Dense Seismic Array Research Articles

Three‐Dimensional Teleseismic Elastic Reverse‐Time Migration With Deconvolution Imaging Condition and Its Application to Southwest Japan

AbstractWe have developed a novel deconvolution‐based reverse time migration method to image lithospheric structures using teleseismic data recorded by dense seismic arrays. Unlike traditional approaches that rely on the retrieval of Green's functions (or receiver functions), the new method directly utilizes the recorded three‐component (3‐C) seismic waveforms to reconstruct subsurface wavefields, which has the advantage of avoiding the estimation and removal of source time function. Importantly, it enables the asymptotic estimation of P‐to‐S transmission conversion coefficients at the elastic discontinuities and enhances resolving power for the fine‐scale heterogeneities. Taking these improvements, we have obtained a full three‐dimensional (3‐D) high‐resolution image of the subduction zone beneath southwest Japan. This image provides a comprehensive configuration of the severely deformed subducting plate beneath the Kii Channel and Kinki Peninsula.

The crust-mantle velocity structure beneath the Wudalianchi-Erkeshan-Keluo volcanic belt by joint inversion of receiver functions and ambient noise

The Wudalianchi-Erkeshan-Keluo (WEK) volcanic belt is a significant component of intraplate volcanism in Northeast China and is composed of the Wudalianchi, Erkeshan, and Keluo volcanic clusters. Using joint inversion of receiver functions and ambient noise, we construct a high-resolution 3-D S-wave velocity model of the WEK volcanic belt and its adjacent region, taking advantage of a deployed dense seismic array around this volcanic belt. There is a prominent low-velocity anomaly at 8–15 km depth beneath the Wudalianchi volcanic cluster, suggesting the presence of a crustal magma chamber. Low-velocity anomalies are also observed at 30–35 km depth beneath the Erkeshan volcanic cluster and 30–40 km depth beneath the Keluo volcanic cluster, resulting in discontinuous velocity structures at the Moho discontinuity. We further find a distinct low-velocity anomaly in the uppermost mantle beneath the WEK volcanic belt. Combined with previous geophysical and geochemistry studies, we propose a magma system scenario for the WEK volcanic belt. The upwelling molten material from the asthenosphere accumulated in the uppermost mantle and the magma chamber was formed, which provided the same uppermost mantle magma sources for the WEK volcanic belt.

Three-dimensional VP, VS, and VP/VS imaging based on AI microseismic detection reveals the mechanism of induced earthquakes in the Xiluodu Reservoir Area, China

The Xiluodu Reservoir (XLD), the third hydropower station developed downstream of the Jinsha River, is currently the fourth largest hydropower station worldwide. The reservoir has been impounded for 10 years, and numerous earthquakes have occurred, thus indicating an increased seismic risk. In existing studies, the deep tectonics and fine-scale structure in the reservoir area have not been sufficiently explained. Thus, the mechanism underlying reservoir-induced seismicity in this area remains unclear. In this study, we utilized a microseismic catalogue compiled based on the LOC-FLOW workflow (Li et al., 2023), alongside long-term observation data from a dense seismic array in the Jinsha River downstream, to obtain high-precision 3D VP, VS and VP/VS structures, as well as refined seismic relocations in the XLD Reservoir area. Additionally, we incorporated petrological, geomagnetic, gravity, crustal density, geodetic observation and numerical simulation data, which demonstrate that seismicity in the reservoir head area and near the Manao Fault is influenced by fluids. Moreover, the two M4 earthquakes that occurred east of the Lianfeng Fault are indicative of tectonic activity. Of these earthquakes, the Tianba earthquake cluster in the reservoir head area is located in a high-VP/VS zone, while a low-VS anomaly channel is observed here, which suggests that large amounts of fluid infiltrated along the fault, increasing the pore pressure and reducing the effective stress on the fault. The Wuji earthquake cluster in the head area is also located in a low-VP/VS zone, which indicates a low porosity and possibly undrained conditions. Thus, the earthquakes in this area are the result of coupled elastic stress-pore water pressure interactions.

Fault-scale crustal structure across the Dunhua-Mishan fault (Tanlu northern segment) constrained from teleseismic P-wave receiver functions

The Dunhua-Mishan fault, located in the northern segment of the Tanlu fault zone, experienced multiple tectonic processes associated with the effects of the Pacific Plate subduction and the Indo-Asia collision. The high-resolution fault-scale structure is critical for understanding the fault evolution and potential fault damage. However, the well-defined deep structure of the Dunhua-Mishan fault is still unclear due to the lack of the dense seismic array. In this study, we construct a high-resolution P-wave receiver function imaging based on linear dense seismic array across the fault. Our results reveal the strong Moho depth variation across the Dunhua-Mishan fault zone. The slightly higher Vp/Vs ratio values within the fault zone indicate the presence of a small amount of mafic crust composition. Interestingly, the significant double positive Ps converted phases are observed within the fault zone, which may represent double Moho discontinuities. The double Moho structure may be related to multiple significant tectonic activities in the Tanlu northern segment. These newly observed structures provide new seismic constraints on the formation and evolution of the Tanlu fault zone and probably reflect that the lithospheric structure of the Dunhua-Mishan fault has been modified by a series of tectonic processes.

Multiscale structures of crust-mantle beneath the South China block and their geodynamic implication

The South China block is an ideal location in which to study the impact of paleo-Pacific plate subduction on the upper mantle and crustal structures of the overriding plate, as well as the mechanism of lithospheric thinning. By integrating data from permanent seismic stations of the China Seismic Network and a dense seismic array at Xuefeng Mountain, this study employed P- and S-wave receiver functions (RFs) to image multiscale structures of crust-mantle across the South China block, thereby providing constraints for a comprehensive understanding of its tectonic evolution. The main findings are as follows: (1) The thickness of lithosphere and crust gradually decreases from northwest (NW) to southeast (SE). (2) Beneath the Xuefeng Mountain uplift zone, the P-to-S conversion from the 660 km discontinuity arrives earlier than that predicted from a one-dimensional model, other weak discontinuities can be traced in the upper mantle, including the lithosphere-asthenosphere boundary, the crust becomes thinner, and the Conrad discontinuity becomes deeper. Referring to previous results from tomography and geological reconstructions, we deduce that lithospheric delamination might have occurred beneath the Xuefeng Mountain uplift zone, and the delaminated lithosphere penetrated the upper mantle and reached the mantle transition zone.

Ocean Bottom Distributed Acoustic Sensing for Oceanic Seismicity Detection and Seismic Ocean Thermometry

AbstractA T‐wave is a seismo‐acoustic wave that can travel a long distance in the ocean with little attenuation, making it valuable for monitoring remote tectonic activity and changes in ocean temperature using seismic ocean thermometry (SOT). However, current high‐quality T‐wave stations are sparsely distributed, limiting the detectability of oceanic seismicity and the spatial resolution of global SOT. The use of ocean bottom distributed acoustic sensing (OBDAS), through the conversion of telecommunication cables into dense seismic arrays, is a cost‐effective and scalable means to complement existing seismic stations. Here, we systematically investigate the performance of OBDAS for oceanic seismicity detection and SOT using a 4‐day Ocean Observatories Initiative community experiment offshore Oregon. We first present T‐wave observations from distant and regional earthquakes and develop a curvelet denoising scheme to enhance T‐wave signals on OBDAS. After denoising, we show that OBDAS can detect and locate more and smaller T‐wave events than regional OBS network. During the 4‐day experiment, we detect 92 oceanic earthquakes, most of which are missing from existing catalogs. Leveraging the sensor density and cable directionality, we demonstrate the feasibility of source azimuth estimation for regional Blanco earthquakes. We also evaluate the SOT performance of OBDAS using pseudo‐repeating earthquake T‐waves. Our results show that OBDAS can utilize repeating earthquakes as small as M3.5 for SOT, outperforming ocean bottom seismometers. However, ocean ambient natural and instrumental noise strongly affects the performance of OBDAS for oceanic seismicity detection and SOT, requiring further investigation.

Impact of Ancient Tectonics on Intracontinental Deformation Partitioning: Insights From Crustal Structures of the East Junggar‐Altai Area

AbstractCompressive stress generated at collision fronts can propagate over long distances, inducing deformation within the continent's interior. Nevertheless, the factors governing the partitioning of intracontinental deformation remain enigmatic. The Altai Mountains serve as a type‐example of ongoing intracontinental deformation. Here, we investigate the crustal architecture of the Chinese Altai Mountains, using receiver functions obtained from newly deployed dense seismic nodal arrays. The new seismic results reveal distinct crustal features, including (a) a negative polarity discontinuity beneath Chinese Altai Mountains, suggesting a low‐velocity layer; (b) a north‐dipping mid‐crustal structure beneath the suture zone between East Junggar and Chinese Altai, indicating underthrusting of East Junggar's lower crust beneath the Chinese Altai Mountains; (c) a double Moho structure beneath East Junggar, revealing a high‐velocity lower crustal layer. In conjunction with constraints from previous multi‐disciplinary regional studies, the double Moho structures are interpreted as mafic restite from Late Paleozoic magma underplating. The addition of mafic materials can significantly enhance the rheological strength of East Junggar's crust, causing it to function as an indenter that thrust beneath the Chinese Altai Mountains during the subsequent convergence process. As a consequence, significant deformation occurs in the Chinese Altai region, resulting in the emergence of decollements, as evident by the negative polarity discontinuity. The presence of pre‐existing decollements makes the Altai Mountains region more susceptible to deformation, thereby facilitating the concentration of intracontinental deformation. These findings illuminate the evolution history of the Chinese Altai Mountains and highlight the great impacts of ancient tectonics on intracontinental deformation partitioning.

Imaging Urban Hidden Faults with Ambient Noise Recorded by Dense Seismic Arrays

Abstract The identification of preexisting near-surface faults represents a piece of crucial information needed to correctly assess the seismic hazard of any area. The mapping of these structures is particularly challenging in densely populated and heavily urbanized areas. We use ambient seismic noise recorded by a dense array in Seal Beach, California, to image shallow fault lines via a reflected surface-wave analysis. Our results highlight the presence of previously unknown shallow faults that correlate remarkably well with shallow seismicity and active survey images.

Wind turbines as a metamaterial-like urban layer: an experimental investigation using a dense seismic array and complementary sensing technologies

The deflection and the control of the effects of the complex urban seismic wavefield on the built environment is a major challenge in earthquake engineering. The interactions between the soil and the structures and between the structures strongly modify the lateral variability of ground motion seen in connection to earthquake damage. Here we investigate the idea that flexural and compressional resonances of tall turbines in a wind farm strongly influence the propagation of the seismic wavefield. A large-scale geophysical experiment demonstrates that surface waves are strongly damped in several distinct frequency bands when interacting at the resonances of a set of wind turbines. The ground-anchored arrangement of these turbines produces unusual amplitude and phase patterns in the observed seismic wavefield, in the intensity ratio between stations inside and outside the wind farm and in surface wave polarization while there is no metamaterial-like complete extinction of the wavefield. This demonstration is done by setting up a dense grid of 400 geophones and another set of radial broadband stations outside the wind farm to study the properties of the seismic wavefield propagating through the wind farm. Additional geophysical equipment (e.g., an optical fiber, rotational and barometric sensors) was used to provide essential explanatory and complementary measurements. A numerical model of the turbine also confirms the mechanical resonances that are responsible for the strong coupling between the wind turbines and the seismic wavefield observed in certain frequency ranges of engineering interest.

SItomo – A toolbox for splitting intensity tomography and application in the Eastern Alps

The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave’s non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.

Seismogenic fault of the 2021 Ms 6.0 Luxian induced earthquake in the Sichuan Basin, China constrained by high-resolution seismic reflection and dense seismic array

Accurate documentation of the location and geometry of seismogenic faults are critical for understanding strong seismicity and seismic hazard mitigation. A rising concern worldwide is induced seismicity caused by hydraulic fracturing, yet the geological characteristics and mechanisms of such seismogenic faults remain insufficiently understood. In particular, a first-order issue is whether the strong (M ≥ 5.0) induced earthquakes are caused directly by the slip of shallow faults near shale reservoirs or indirectly by deep basement faults. The September 16, 2021 Ms 6.0 (Mw 5.4) Luxian earthquake in the Sichuan Basin, China is one of the world's greatest induced earthquakes related to shale gas exploitation. Here we constrain the geometry of the seismogenic fault of this earthquake using three-dimensional seismic reflection data and a dense seismic array. Pre-existing faults in the sedimentary cover are clearly visible in seismic coherence slices. However, these shallow faults do not match the seismological parameters of the Luxian earthquake. High-resolution seismic reflection profiles combined with relocated earthquakes indicate that the seismogenic fault is located in the Precambrian basement. The results suggest that the mainshock is most likely caused by the poroelastic effects due to fluid injection. Hydraulic fracturing could have reactivated a large-scale basement fault and triggered the strong earthquakes under a relatively high geo-stress conditions in the study area. The basement fault slipped upward and ruptured the sedimentary cover, resulting in a significant difference between the focal and centroid depths of the mainshock. Our findings emphasize the importance of investigating regional tectonic and geological settings, local stress fields, and pre-existing faults for studying potential induced seismicity in shale gas fields.

Resolvability of Multiple Microseismic P-Wave Source Regions with Two Large Seismic Arrays in China and the United States

Abstract P-wave microseisms are useful for understanding ocean waves. Resolving and locating multiple P-wave source regions using seismic data can provide valuable information about ocean waves. The resolvability of multiple microseismic P-wave source regions depends on the location accuracy and resolution, which can be improved using multiple large seismic arrays. In this article, we investigate the source locations of P-wave microseisms at the period of 5–10 s by combining the backprojection results from two large dense seismic arrays located in China (ChinArray) and the United States (USArray). We independently process data recorded by ChinArray and USArray during a two-year period (2014–2015) that border both the North Pacific and North Atlantic. Then the results are normalized and summed or intersected in the source region to improve the accuracy of the P-wave microseism source locations by reducing the deviation from the velocity structure model and the array response function. The results show that we can resolve two to three sources with a scale of ∼500–1000 km within one large P-wave source region. We also investigate how array parameters such as aperture, interstation spacing, and geographic position affect the detectability and accuracy of the P-wave microseism sources. The discrepancy in P-wave microseism source locations between backprojection observation and ocean model predictions in source number, source scale, and source region scope imply that the ocean model needs to be improved.

Coherent Wavefield Reconstruction Improves Event Location With Dense Seismic Arrays

Coherent Wavefield Reconstruction Improves Event Location With Dense Seismic Arrays

Earthquake monitoring and high-resolution velocity tomography for the central Longmenshan fault zone by a temporary dense seismic array

Earthquake monitoring and high-resolution velocity tomography for the central Longmenshan fault zone by a temporary dense seismic array

Detailed crustal deformation in the Ailaoshan Orogenic Belt revealed by receiver functions from a dense array

The Ailaoshan Orogenic Belt (AOB), located at the southeastern boundary of the Tibetan Plateau, is an ideal place for investigating the mechanisms of lateral growth of Tibet. Using the data recorded by a dense seismic array across the Ailaoshan belt, we investigate the detailed lateral variations of crustal anisotropy on the basis of Pms phase of receiver functions. Remarkable crustal anisotropy is observed throughout this study region with a mean delay time of 0.33 ± 0.19 s, indicating the anisotropy primarily originates in the middle-lower crust. The fast directions beneath the AOB including the Ailaoshan-Red River shear zone (ARRSZ) and its western low-grade metamorphic unit generally align with the NW-SE strike of ARRSZ. The weak anisotropy in the South China Block (SCB) argues that the block is relatively stable, with limited internal deformation. Meanwhile, the anisotropy beneath the western boundary of the SCB is strong, and the N-S oriented fast direction is influenced by both the crustal stress and Xiaojiang Fault. Combining the high Vp/Vs and significant lateral variations of crustal anisotropy parameters, we suggest that the strike-slip motion along the ARRSZ induces the partial melting and pronounced anisotropy in the middle-lower crust of AOB, without the presence of crustal flow. The differences between crustal and mantle anisotropy indicate crust-mantle decoupling deformation of the AOB, supporting the block extrusion model occurring only in the crustal scale as the primary deformation pattern.

Determining the location of the Bayan Obo rare earth elements mineralization body by the transfer learning method

Bayan Obo is the largest rare earth element (REE) deposit in the world. The occurrence of REE is closely related to the dolomite in this area. Dolomite serves as the mother rock of REE mineralization and the ore body. How to accurately locate and characterize dolomite is the key to determining the distribution of REE and estimating its reserves. A large amount of geophysical work has been conducted in this area, including a dense seismic array, various electromagnetic methods, gravity and aeromagnetic surveys, as well as numerous petrophysical property measurements. To fully leverage the results obtained by these geophysical methods and develop an understanding of the physical property structure, a multisource geophysical data fusion technology is developed. First, various physical property profiles obtained from inversion on the same profile are converted into images with identical resolution and dimension. Then, an image adaptive feature extraction technique based on transfer learning is used to extract features of different scales from multisource images. Subsequently, the fusion image is reconstructed based on the local nearest neighbor weighted average feature fusion rule to obtain the final fusion result. This aids in identifying the spatial appearance pattern of the target for detection. Given the physical characteristics of the mineralized dolomite, which has high density, high resistivity, and high magnetic susceptibility, its location and shape can be defined in the fusion image. The results indicate that the occurrence depth of dolomite can extend up to 1500 m and the dolomite has a southward tilt as one of its primary structural characteristics. The predicted range of dolomite distribution is consistent with the formation range revealed by drilling, making it a reliable basis for predicting the distribution of rare earth ore bodies.

Detection of Low-Velocity Layer Using Passive Surface-Wave Data and Deep Learning: An Application in The Urban Area of Hangzhou, China

Passive surface-wave methods using dense seismic arrays have gained growing attention in near-surface high-resolution imaging in urban environments. Deep learning (DL) can release a tremendous workload brought by dense seismic arrays. We presented a case study of shear-wave velocity (Vs) structure imaging in the Hangzhou urban area (eastern China) using DL inversion. Noise data were recorded by dense linear arrays with approximately 5 m spacing deployed along two crossing roads for investigating the top 80 m of the subsurface. Phase-velocity dispersion curves are extracted from virtual shot gathers using multichannel analysis of surface waves. We divided the area where the low-velocity layer (LVL) may exist into three layers with a thickness of 5 m. We gave the four layers weak constraints to generate training dataset and adopted a convolutional neural network to directly invert fundamental-mode Rayleigh-wave phase velocity for 1D Vs models. To improve the accuracy, we further applied the sensitivities to weight the loss function in DL inversion. The obtained pseudo-2D Vs profiles correspond to the velocities estimated from logging data and previous survey. The well-trained neural network successfully identified that the LVL is located at 50-60 m deep. And this network was also achieved accurately the inversion of a dense seismic network nearby. The results of this survey demonstrate the accuracy and efficiency of delineating near-surface structures from traffic-induced noise using the DL technique, which has great potential for monitoring subsurface changes in urban areas.

Physics‐Informed Neural Networks for Elliptical‐Anisotropy Eikonal Tomography: Application to Data From the Northeastern Tibetan Plateau

AbstractWe develop a novel approach for multi‐frequency, elliptical‐anisotropic eikonal tomography based on physics‐informed neural networks (pinnEAET). This approach simultaneously estimates the medium properties controlling anisotropic Rayleigh waves and reconstructs the traveltimes. The physics constraints built into pinnEAET's neural network enable high‐resolution results with limited inputs by inferring physically plausible models between data points. Even with a single source, pinnEAET can achieve stable convergence on key features where traditional methods lack resolution. We apply pinnEAET to ambient noise data from a dense seismic array (ChinArray‐Himalaya II) in the northeastern Tibetan Plateau with only 20 quasi‐randomly distributed stations as sources. Anisotropic phase velocity maps for Rayleigh waves in the period range from 10–40 s are obtained by training on observed traveltimes. Despite using only about 3% of the total stations as sources, our results show low uncertainties, good resolution and are consistent with results from conventional tomography.

Application of Horizontal-to-vertical Spectral Ratio with a Dense Linear Array on Prospecting Urban Underground Active Faults: a Case in The Mingguang City, Anhui Province

Reliable knowledge of sedimentary structures is beneficial for the study of active faults covered by sediments and can also provide important information for seismic hazard assessment, earthquake emergency response and urban planning. Based on the dense seismic array with a spacing of ∼ 60 m between stations and a length of ∼ 8 km, the horizontal-to-vertical spectral ratio (HVSR) curves were obtained at 133 stations with the ambient noise, and the two-dimensional sedimentary structures of the Quaternary and Neogene strata were imaged over two main faults of the Tanlu fault zone in the eastern Jiashan Basin by frequency-depth conversion. The variations (∼ 130 m) in the depth of the Quaternary and Neogene sedimentary basement directly beneath the fault zone agree well with the spatial characteristics of the fault revealed by the shallow seismic reflection profile with active sources at the same locations.

First‐Order Transition in Appalachian Orogenic Processes Revealed by Along‐Strike Variation of the Moho Geometry

AbstractAlong‐strike variation of the Laurentian rifted margin and the Appalachian orogen has long been recognized in the geologic record. We investigated the manifestation of this along‐strike variation at depth by generating scattered wavefield migration profiles from four dense seismic arrays deployed across the Appalachian orogen at different latitudes. All profiles exhibit a similar crustal thickness decrease of 15–20 km from the Mesoproterozoic Grenville Province to the Paleozoic Appalachian accreted terranes, but the Moho architecture differs dramatically along strike. The profiles beneath the central and southern Appalachians show a smoothly varying Moho geometry; in contrast, there is an abrupt Moho depth offset beneath the New England Appalachians. This contrast in Moho geometry may result from variations in the Laurentian rifted margin architecture, changes in Taconic orogeny subduction polarity, and greater crustal shortening during the Acadian‐Neoacadian orogeny in southern New England and the Alleghanian orogeny in the central and southern Appalachians. A first‐order along‐strike transition in the behavior of Appalachian orogenic processes is located between the central and New England Appalachians.