Lupei Zhu

Moho depth variation in southern California from teleseismic receiver functions

The number of broadband three-component seismic stations in southern California has more than tripled recently. In this study we use the teleseismic receiver function technique to determine the crustal thicknesses and V_p/V_s ratios for these stations and map out the lateral variation of Moho depth under southern California. It is shown that a receiver function can provide a very good “point” measurement of crustal thickness under a broadband station and is not sensitive to crustal P velocity. However, the crustal thickness estimated only from the delay time of the Moho P-to-S converted phase trades off strongly with the crustal V_p/V_s ratio. The ambiguity can be reduced significantly by incorporating the later multiple converted phases, namely, the PpPs and PpSs+PsPs. We propose a stacking algorithm which sums the amplitudes of receiver function at the predicted arrival times of these phases by different crustal thicknesses H and Vp/Vs ratios. This transforms the time domain receiver functions directly into the H-V_p/V_s domain without need to identify these phases and to pick their arrival times. The best estimations of crustal thickness and V_p/V_s ratio are found when the three phases are stacked coherently. By stacking receiver functions from different distances and directions, effects of lateral structural variation are suppressed, and an average crustal model is obtained. Applying this technique to 84 digital broadband stations in southern California reveals that the Moho depth is 29 km on average and varies from 21 to 37 km. Deeper Mohos are found under the eastern Transverse Range, the Peninsular Range, and the Sierra Nevada Range. The central Transverse Range, however, does not have a crustal root. Thin crusts exist in the Inner California Borderland (21–22 km) and the Salton Trough (22 km). The Moho is relatively flat at the average depth in the western and central Mojave Desert and becomes shallower to the east under the Eastern California Shear Zone (ECSZ). Southern California crust has an average V_p/V_s ratio of 1.78, with higher ratios of 1.8 to 1.85 in the mountain ranges with Mesozoic basement and lower ratios in the Mojave Block except for the ECSZ, where the ratio increases.

Crustal structure across the San Andreas Fault, southern California from teleseismic converted waves

Crustal structure along a 140 km long profile across the San Andreas Fault (SAF) in southern California was imaged by stacking teleseismic P–S converted phases recorded by a dense, short-period seismic array. The crust/mantle discontinuity (Moho) is visible as a continuous feature at a depth around 30 km but is offset 6 to 8 km beneath the SAF. A small Moho disruption can also be seen under the Eastern California Shear Zone (ECSZ). These results suggest that the SAF and ECSZ extend through the entire crust. The Moho upwarp under the San Gabriel Mountains indicates that the mountain ranges were lifted en masse as a result of crustal buckling under horizontal compression.

Intermediate depth earthquakes beneath the India-Tibet Collision Zone

We report on three intermediate depth earthquakes in the India‐Tibet collision zone, two under the Himalayan Thrust Belt (HTB) and one beneath the Indus Zangbo suture. The mb magnitudes of these three events are from 4.3 to 4.9, and are too small to be well located by conventional means. However, from modeling their broadband waveforms recorded at near‐regional distances on a temporary PASSCAL array, we can confidently confine the sources to be below the crust, between 70 and 80 km deep. The existence of these intermediate depth earthquakes in this area suggests relatively low temperatures in the mantle lithosphere. The two events under the HTB display strike‐slip mechanisms with some normal faulting component; this is quite different from the shallow thrust events typical of the same area. The principal P and T axes of all 3 fault plane solutions show roughly NS compression and EW extension, consistent with a regional stress field produced by the indenting of the India continent.

Mushy magma beneath Yellowstone

A recent prospective on the Yellowstone Caldera discounts its explosive potential based on inferences from tomographic studies which suggests a high degree of crystallization of the underlying magma body. In this study, we show that many of the first teleseismic P-wave arrivals observed at seismic stations on the edge of the caldera did not travel through the magma body but have taken longer but faster paths around the edge. After applying a number of waveform modeling tools, we obtain much lower seismic velocities than previous studies, 2.3 km/sec (V_p) and 1.1 km/sec (V_s). We estimate the physical state of the magma body by assuming a fluid-saturated porous material consisting of granite and a mixture of rhyolite melt and water and CO_2 at a temperature of 800°C and pressure at 5 km (0.1 GPa). We found that this relatively shallow magma body has a volume of over 4,300 km^3 and is about 32% melt saturated with about 8% water plus CO_2 by volume.

Fault-Plane Determination of the 18 April 2008 Mount Carmel, Illinois, Earthquake by Detecting and Relocating Aftershocks

We developed a sliding-window cross-correlation (SCC) detection technique and applied the technique to continuous waveforms recorded by the Cooperative New Madrid Seismic Network stations following the 18 April 2008 Illinois earthquake. The technique detected more than 120 aftershocks down to M_L 1.0 in the 2 week time window following the mainshock, which is three times more than the number of aftershocks reported by the seismic network. Most aftershocks happened within 24 hrs of the mainshock. We then relocated all events by the double-difference relocation algorithm. Accurate P- and S-wave differential arrival times between events were obtained by waveform cross correlation. After relocation, we used the L1 norm to fit all located events by a plane to determine the mainshock fault plane. The best-fit plane has a strike of 292°±11° and dips 81°±7° to the northeast. This plane agrees well with the focal mechanism solutions of the mainshock and four largest aftershocks. By combining the aftershock locations and focal mechanism solutions, we conclude that the 18 April earthquake occurred on a nearly vertical left-lateral strike-slip fault orienting in the west-northwest–east-southeast direction. The fault coincides with the proposed left-stepping Divide accommodation zone in the La Salle deformation belt and indicates reactivation of old deformation zone by contemporary stresses in the Midcontinent.

Regional waveform calibration in the Pamir-Hindu Kush region

Twelve moderate-magnitude earthquakes (m_b 4–5.5) in the Pamir-Hindu Kush region are investigated to determine their focal mechanisms and to relocate them using their regional waveform records at two broadband arrays, the Kyrgyzstan Regional Network (KNET), and the 1992 Pakistan Himalayas seismic experiment array (PAKH) in northern Pakistan. We use the “cut-and-paste” source estimation technique to invert the whole broadband waveforms for mechanisms and depths, assuming a one-dimensional velocity model developed for the adjacent Tibetan plateau. For several large events the source mechanisms obtained agree with those available from the Harvard centroid moment tensor (CMT) solutions. An advantage of using regional broadband waveforms is that focal depths can be better constrained either from amplitude ratios of Pnl to surface waves for crustal events or from time separation between the direct P and the shear-coupled P wave (sPn + sPmP) for mantle events. All the crustal events are relocated at shallower depths compared with their International Seismological Centre bulletin or Harvard CMT depths. After the focal depths are established, the events are then relocated horizontally using their first-arrival times. Only minor offsets in epicentral location are found for all mantle events and the bigger crustal events, while rather large offsets (up to 30 km) occur for the smaller crustal events. We also tested the performance of waveform inversion using only two broadband stations, one from the KNET array in the north of the region and one from the PAKH array in the south. We found that this geometry is adequate for determining focal depths and mechanisms of moderate size earthquakes in the Pamir-Hindu Kush region.

Crustal Velocity Structure in Italy from Analysis of Regional Seismic Waveforms

Abstract In this article, we use regional seismic waveforms recorded by the recently installed Istituto Nazionale di Geofisica e Vulcanologia (INGV) national network and the Mediterranean Very Broadband Seismographic Network (MedNet) stations to develop 1D crustal velocity models for the Italian peninsula. About 55,000 P -wave and 35,000 S -wave arrival times from 4727 events are used to derive average seismic parameters in the crust and uppermost mantle. We define four regions, according to geological constraints and recent travel-time tomography results. Based on the average seismic parameters, we combine broadband seismic waveforms and travel times of regional phases to model crustal structures for the four regions by applying the genetic algorithm. Our results indicate smooth velocity gradients with a depth beneath the Apennines and a deep Moho beneath the central Alps. Green’s functions from the regionalized 1D velocity models are used to determine source depths and focal mechanisms for 37 events with a magnitude larger than 3.5 by a grid search technique. Our results show that normal and strike-slip faulting source mechanisms dominate the Apenninic belt and that most thrust faulting events occur in the Adriatic Sea and the outer margin of the northern Apennines.

Determination of earthquake focal depths and source time functions in central Asia using teleseismic P waveforms

We developed a new method to determine earthquake source time functions and focal depths. It uses theoretical Greens function and a time-domain deconvolution with positivity constraint to estimate the source time function from the teleseismic P waveforms. The earthquake focal depth is also determined in the process by using the time separations of the direct P and depth phases. We applied this method to 606 earthquakes between 1990 and 2005 in Central Asia. The results show that the Centroid Moment Tensor solutions, which are routinely computed for earthquake larger than M5.0 globally using very long period body and surface waves, systematically over-estimated the source depths and durations, especially for shallow events. Away from the subduction zone, most of the 606 earthquakes occurred within the top 20 km of crust. This shallow distribution of earthquakes suggests a high geotherm and a weak ductile lower crust in the region.

Deformation in the lower crust and downward extent of the San Andreas Fault as revealed by teleseismic waveforms

High resolution images of crustal structure across the San Andreas Fault (SAF) were obtained by using the common conversion point stacking of teleseismic P-to-S converted waves recorded during the Los Angeles Region Seismic Experiments (LARSE-I and II). In the upper crust, several sedimentary basins were delineated in the images, including the San Fernando and the Santa Clarita Basins. The San Fernando Basin reaches a depth of 8 km under the northern edge of the San Fernando Valley. On the LARSE-I profile, the downward projection of the SAF truncates several lower crustal interfaces including the Moho on both sides. The Moho is vertically offset by as much as 8 km. Along the LARSE-II profile, the impedance contrast and slope of the Moho are seen to change across the fault. These results indicate that the fault penetrates into the lower crust and probably uppermost mantle as a narrow (<10 km) feature. The Moho beneath the San Gabriel Mountains is shallower (∼26 km) than under the San Gabriel Valley to the south and the Mojave Desert to the north, suggesting that the mountain ranges were lifted en masse by horizontal compression. On the northeast side of the SAF, the Mojave Desert has a sharp and essentially flat Moho at a depth of ∼32 km. The lower crustal structure beneath the San Fernando and Santa Clarita Valleys along the LARSE-II profile south of the SAF is complicated as indicated by the large undulation and low impedance contrast of the Moho. These observations suggest that the deformation in the lower crust is localized and often concentrates near boundaries of crustal blocks or beneath those places which have experienced intensive faulting and deformation in the upper crust.

Preliminary Study of Crust-Upper Mantle Structure of the Tibetan Plateau by Using Broadband Teleseismic Body Waveforms

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure.Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km.