Attenuation tomography. Modelling regional love waves: Imperial Valley to Pasadena

Abstract to Part I

The inverse problem of seismic wave attenuation is solved by an iterative back-projection method. The seismic wave quality factor, Q, can be estimated approximately by inverting the S-to-P amplitude ratios. Effects of various uncertain ties in the method are tested and the attenuation tomography is shown to be useful in solving for the spatial variations in attenuation structure and in estimating the effective seismic quality factor of attenuating anomalies.

Back-projection attenuation tomography is applied to two cases in southern California: Imperial Valley and the Coso-Indian Wells region. In the Coso-Indian Wells region, a highly attenuating body (S-wave quality factor (Q_β ≈ 30) coincides with a slow P-wave anomaly mapped by Walck and Clayton (1987). This coincidence suggests the presence of a magmatic or hydrothermal body 3 to 5 km deep in the Indian Wells region. In the Imperial Valley, slow P-wave travel-time anomalies and highly attenuating S-wave anomalies were found in the Brawley seismic zone at a depth of 8 to 12 km. The effective S-wave quality factor is very low (Q_β ≈ 20) and the P-wave velocity is 10% slower than the surrounding areas. These results suggest either magmatic or hydrothermal intrusions, or fractures at depth, possibly related to active shear in the Brawley seismic zone.

No-block inversion is a generalized tomographic method utilizing the continuous form of an inverse problem. The inverse problem of attenuation can be posed in a continuous form , and the no-block inversion technique is applied to the same data set used in the back-projection tomography. A relatively small data set with little redundancy enables us to apply both techniques to a similar degree of resolution. The results obtained by the two methods are very similar. By applying the two methods to the same data set, formal errors and resolution can be directly computed for the final model, and the objectivity of the final result can be enhanced.

Both methods of attenuation tomography are applied to a data set of local earthquakes in Kilauea, Hawaii, to solve for the attenuation structure under Kilauea and the East Rift Zone. The shallow Kilauea magma chamber, East Rift Zone and the Mauna Loa magma chamber are delineated as attenuating anomalies. Detailed inversion reveals shallow secondary magma reservoirs at Mauna Ulu and Puu Oo, the present sites of volcanic eruptions. The Hilina Fault zone is highly attenuating, dominating the attenuating anomalies at shallow depths. The magma conduit system along the summit and the East Rift Zone of Kilauea shows up as a continuous supply channel extending down to a depth of approximately 6 km. The Southwest Rift Zone, on the other hand, is not delineated by attenuating anomalies, except at a depth of 8-12 km, where an attenuating anomaly is imaged west of Puu Kou. The Ylauna Loa chamber is seated at a deeper level (about 6-10 km) than the Kilauea magma chamber. Resolution in the Mauna Loa area is not as good as in the Kilauea area, and there is a trade-off between the depth extent of the magma chamber imaged under Mauna Loa and the error that is due to poor ray coverage. Kilauea magma chamber, on the other hand, is well resolved, according to a resolution test done at the location of the magma chamber.

Abstract to Part II

Long period seismograms recorded at Pasadena of earthquakes occurring along a profile to Imperial Valley are studied in terms of source phenomena (e.g., source mechanisms and depths) versus path effects. Some of the events have known source parameters, determined by teleseismic or near-field studies, and are used as master events in a forward modeling exercise to derive the Green's functions (SH displacements at Pasadena that are due to a pure strike-slip or dip-slip mechanism) that describe the propagation effects along the profile. Both timing and waveforms of records are matched by synthetics calculated from 2-dimensional velocity models. The best 2-dimensional section begins at Imperial Valley with a thin crust containing the basin structure and thickens towards Pasadena. The detailed nature of the transition zone at the base of the crust controls the early arriving shorter periods (strong motions), while the edge of the basin controls the scattered longer period surface waves. From the waveform characteristics alone, shallow events in the basin are easily distinguished from deep events, and the amount of strike-slip versus dip-slip motion is also easily determined. Those events rupturing the sediments, such as the 1979 Imperial Valley earthquake, can be recognized easily by a late-arriving scattered Love wave that has been delayed by the very slow path across the shallow valley structure.

Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan

The fine-scale seismic structure of the central Mexico, southern Peru, and southwest Japan subduction zones is studied using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. The morphology of the transition from flat to normal subduction is explored in central Mexico and southern Peru, while in southwest Japan the spatial coincidence of a thin ultra-slow velocity layer (USL) atop the flat slab with locations of slow slip events (SSEs) is explored. This USL is also observed in central Mexico and southern Peru, where its lateral extent is used as one constraint on the nature of the flat-to-normal transitions.

In western central Mexico, I find an edge to this USL which is coincident with the western boundary of the projected Orozco Fracture Zone (OFZ) region. Forward modeling of the 2D structure of the subducted Cocos plate using a finite-difference algorithm provides constraints on the velocity and geometry of the slab’s seismic structure in this region and confirms the location of the USL edge. I propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, by a process analogous to that which occurred when the Rivera plate separated from the proto-Cocos plate 10 Ma.

In eastern central Mexico, observations of a sharp transition in slab dip near the abrupt end of the Trans Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of the USL is found to be coincident with these features and with the western boundary of a zone of decreased seismicity, indicating a change in structure which I interpret as evidence of a possible tear. Analysis of intraslab seismicity patterns and focal mechanism orientations and faulting types provides further support for a possible tear in the South Cocos slab. This potential tear, together with the tear along the projection of the OFZ to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments.

In southern Peru, observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest a smooth contortion of the slab. Concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry. The lateral extent of the USL atop the horizontal Nazca slab is found to be coincident with the margin of the projected linear continuation of the subducting Nazca Ridge, implying a causal relationship, but not a slab tear. Waveform modeling of the 2D structure in southern Peru provides constraints on the velocity and geometry of the slab’s seismic structure and confirms the absence of any tears in the slab.

In southwest Japan, I estimate the location of a possible USL along the Philippine Sea slab surface and find this region of low velocity to be coincident with locations of SSEs that have occurred in this region. I interpret the source of the possible USL in this region as fluids dehydrated from the subducting plate, forming a high pore-fluid pressure layer, which would be expected to decrease the coupling on the plate interface and promote SSEs.

Electromagnetic decay of the 1.7- and 2.43-MeV levels in ^9Be

The 1.7- and 2.43-MeV levels in ⁹Be were populated with the reaction ¹¹B(d, α)⁹Be* by bombarding thin boron on carbon foils with 1.7-MeV deuterons. The alpha particles were analyzed in energy with a surface-barrier counter set at the unique kinematically determined angle and the recoiling ⁹Be nuclei at 90^o were analyzed in rigidity with a magnetic spectrometer, in energy by a surface-barrier counter at the spectrometer focus, and in velocity by the time delay between an alpha and a ⁹Be count. When a pulse from the spectrometer counter was in the appropriate delayed coincidence with a pulse from the alpha counter, the two pulses were recorded in a two-dimensional pulse height analyzer. Most of the ⁹Be* decay by particle breakup. Only those that gamma decay are detected by the spectrometer counter. Thus the experiment provides a direct measurement of Γ_rad/Γ. Analysis of 384 observed events gives Γ_rad/Γ = (1.16 ± 0.14) X 10^-4 for the 2.43-MeV level. Combining this ratio with the value of Γ_rad = 0.122 ± 0.015 eV found from inelastic electron scattering gives Γ = (1.05 ± 0.18) keV. For the 1.7-MeV level, an upper limit, Γ_rad/Γ ≤ 2.4 = 10^-5, was determined.

An investigation of the three-body reaction He^3(He^3, 2p)He^4

Reactions produced by the He³ bombardment of the He³ have been investigated for bombarding energies from 1 to 20 MeV using a tandem Van de Graaff accelerator. Proton spectra from the three-body reaction He³(He³, 2p)He⁴ have been measured with a counter telescope at 13 angles for 9 bombarding energies between 3 and 18 MeV. The results are compared with a model for the reaction which includes a strong p-He⁴ final-state interaction. Alpha-particle spectra have been obtained at 12 and 18 MeV for forward angles with a magnetic spectrometer. These spectra indicate a strongly forward-peaked mechanism involving the ¹S₀ p-p interaction in addition to the p-He⁴ interaction. Measurements of p-He⁴ and p-p coincidence spectra at 10 MeV confirm these features of the reaction mechanism. Deuteron spectra from the reaction of He³(He³, d)pHe³ have been measured at 18 MeV. A triton spectrum from the reaction He³(He³, t)3p at 20 MeV and 4⁰ is interpreted in terms of a sequential decay through an excited state of the alpha particle at 20.0 MeV. No effects are observed which would indicate an interaction in the residual (3p) system. Below 3 MeV the He³(He³, 2p)He⁴ reaction mechanism is observed to be changing and further measurements are suggested in view of the importance of this reaction in stellar interiors.

A study of the ^9 Be nucleus

The reaction ⁷Li(³He, p)⁹Be has been used to measure excitations and intrinsic widths of levels in ⁹Be below the ⁷Li + d threshold. Previously unreported levels have been found at excitations of (13.78 ± .03) MeV and (16.671 ± .008) MeV with widths of (590 ± 60) keV and (41 ± 4) keV respectively. Two overlapping levels have been found at (11.81 ± .02) MeV and (11.29 ± .03) MeV with widths of (400 ± 30) keV and (620 ± 70) keV respectively. Branching ratios from ⁹Be levels populated in this reaction to the ground and first excited states of ⁸Be have been measured by observing the associated protons in coincidence with the decay neutrons. Branching ratios were found to be:

Excitation in ⁹Be .... Branching Ratio.......... Final Nucleus.........

(MeV) .......................... (percent) .....................................

.. 2.43 ........................... 7.5 ± 1.5 .............. ⁸Be(g.s.)

.. 3.03 ........................... 87 ± 13......................................

.. 4.65 ........................... 13 ± 4.......................................

.. 6.76 .............................. ≤ 2 ......................................

.. 11.29 ...............................≤ 2 ......................................

.. 11.81 ...............................≤ 3 ......................................

.. 6.76 ........................... .41 ≤ B.R. ≤ .69 ....... ⁸Be(2⁺)

.. 11.29 ........................... 14 ± 4 .......................................

.. 11.81 ........................... 12 ± 4 .......................................

Corresponding reduced widths for neutron emission are calculated and a comparison of the results with the expectations of current nuclear models is made. In particular the measured branching ratio to ⁸Be(g.s.) from ⁹Be(2.43 MeV) corresponds to an f-wave reduced with θ²_f = 2.1 x 10^-2, in units of ħ²/mR², with R = 4.35 fm. A comparison of this value with that predicted by a Nilsson model calculation, in which ⁹Be is taken to be a deformed nucleus, is discussed. The measured value for θ²_f is found to be consistent with that expected on the basis of measured E2-transition rates between rotational levels in ⁹Be.

Firing Patterns of Cerebellar Purkinje Cells During Locomotion and Sleep

The cerebellum is a major supraspinal center involved in the coordination of movement. The principal neurons of the cerebellar cortex, Purkinje cells, receive excitatory synaptic input from two sources: the parallel and climbing fibers. These pathways have markedly different effects: the parallel fibers control the rate of simple sodium spikes, while the climbing fibers induce characteristic complex spike bursts, which are accompanied by dendritic calcium transients and play a key role in regulating synaptic plasticity. While many studies using a variety of species, behaviors, and cerebellar regions have documented modulation in Purkinje cell activity during movement, few have attempted to record from these neurons in unrestrained rodents. In this dissertation, we use chronic, multi-tetrode recording in freely-behaving rats to study simple and complex spike firing patterns during locomotion and sleep. Purkinje cells discharge rhythmically during stepping, but this activity is highly variable across steps. We show that behavioral variables systematically influence the step-locked firing rate in a step-phase-dependent way, revealing a functional clustering of Purkinje cells. Furthermore, we find a pronounced disassociation between patterns of variability driven by the parallel and climbing fibers, as well as functional differences between cerebellar lobules. These results suggest that Purkinje cell activity not only represents step phase within each cycle, but is also shaped by behavior across steps, facilitating control of movement under dynamic conditions. During sleep, we observe an attenuation of both simple and complex spiking, relative to awake behavior. Although firing rates during slow wave sleep (SWS) and rapid eye movement sleep (REM) are similar, simple spike activity is highly regular in SWS, while REM is characterized by phasic increases and pauses in simple spiking. This phasic activity in REM is associated with pontine waves, which propagate into the cerebellar cortex and modulate both simple and complex spiking. Such a temporal coincidence between parallel and climbing fiber activity is known to drive plasticity at parallel fiber synapses; consequently, pontocerebellar waves may provide a mechanism for tuning synaptic weights in the cerebellum during active sleep.

The Physiology and Computation of Pyramidal Neurons

A variety of neural signals have been measured as correlates to consciousness. In particular, late current sinks in layer 1, distributed activity across the cortex, and feedback processing have all been implicated. What are the physiological underpinnings of these signals? What computational role do they play in the brain? Why do they correlate to consciousness? This thesis begins to answer these questions by focusing on the pyramidal neuron. As the primary communicator of long-range feedforward and feedback signals in the cortex, the pyramidal neuron is set up to play an important role in establishing distributed representations. Additionally, the dendritic extent, reaching layer 1, is well situated to receive feedback inputs and contribute to current sinks in the upper layers. An investigation of pyramidal neuron physiology is therefore necessary to understand how the brain creates, and potentially uses, the neural correlates of consciousness. An important part of this thesis will be in establishing the computational role that dendritic physiology plays. In order to do this, a combined experimental and modeling approach is used.

This thesis beings with single-cell experiments in layer 5 and layer 2/3 pyramidal neurons. In both cases, dendritic nonlinearities are characterized and found to be integral regulators of neural output. Particular attention is paid to calcium spikes and NMDA spikes, which both exist in the apical dendrites, considerable distances from the spike initiation zone. These experiments are then used to create detailed multicompartmental models. These models are used to test hypothesis regarding spatial distribution of membrane channels, to quantify the effects of certain experimental manipulations, and to establish the computational properties of the single cell. We find that the pyramidal neuron physiology can carry out a coincidence detection mechanism. Further abstraction of these models reveals potential mechanisms for spike time control, frequency modulation, and tuning. Finally, a set of experiments are carried out to establish the effect of long-range feedback inputs onto the pyramidal neuron. A final discussion then explores a potential way in which the physiology of pyramidal neurons can establish distributed representations, and contribute to consciousness.