Search for: All records

Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical‐component GSC velocity as an approximation for the incident source wavelet. The MohoP‐to‐sconversions are clearly visible on DAS receiver functions. We then derive meter‐scale arrival time measurements along the entire 20‐km‐long array. We are also able to measure path‐averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future.

 

The 2010M_W7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events (M≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus onPnand its depth phasespPnandsPn, which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stackPnwaveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of thePn, directP, and directSphases recorded for event pairs. Focal depths of 93 out of 122M≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.

 

The structure of the mantle transition zone (MTZ) in subduction zones is essential for understanding subduction dynamics in the deep mantle and its surface responses. We constructed the P (V_p) and SH velocity (V_s) structure images of the MTZ beneath Northeast Asia based on two‐dimensional (2‐D) triplicated waveform modeling. In the upper MTZ, a normalV_pbut 2.5% lowV_slayer compared with IASP91 are required by the triplication data. In the lower MTZ, our results show a relatively higher‐velocity layer (+2%V_pand −0.5%V_scompared to IASP91) with a thickness of ~140 km and length of ~1,200 km atop the 660‐km discontinuity. Taking this anomaly as the stagnant slab and considering the plate convergence rate of 7–10 cm/year in the western Pacific region during the late Cenozoic, we deduced that the stagnant slab has a subduction age of less than 30 Ma. This suggests that the observed stagnancy of the slab in the MTZ beneath Northeast Asia may have occurred no earlier than the Early Oligocene. From the constraints derived individually onV_pandV_sstructures, highV_p/V_sratios are obtained for the entire MTZ beneath Northeast Asia, which may imply a water‐rich and/or carbonated environment. Within the overall higher‐velocity stagnant slab, a low‐velocity anomaly was further detected, with a width of ~150 km,V_pandV_sreductions of 1% and 3% relative to IASP91. Such a gap may have provided a passage for hot deep mantle materials to penetrate through the thick slab and feed the Changbaishan volcano.

 

SUMMARY Receiver functions (RFs) estimated on dense arrays have been widely used for the study of Earth structures across multiple scales. However, due to the ill-posedness of deconvolution, RF estimation faces challenges such as non-uniqueness and data overfitting. In this paper, we present an array-based RF deconvolution method in the context of emerging dense arrays. We propose to exploit the wavefield coherency along a dense array by joint inversions of waveforms from multiple events and stations for RFs with a minimum number of phases required by data. The new method can effectively reduce the instability of deconvolution and help retrieve RFs with higher fidelity. We test the algorithm on synthetic waveforms and show that it produces RFs with higher interpretability than those by the conventional RF estimation practice. Then we apply the method to real data from the 2016 Incorporated Research Institutions for Seismology (IRIS) community wavefield experiment in Oklahoma and are able to generate high-resolution RF profiles with only three teleseismic earthquakes recorded by the temporary deployment. This new method should help enhance RF images derived from short-term high-density seismic profiles.