Shallow Velocity Research Articles

Probing a planet from the subsurface to the atmosphere with infrasound data

The analysis of infrasound waves has a significant potential for retrieving a range of geophysical parameters across various scales, such as atmospheric structures, characteristics of surface and buried sources, and seismic velocity structures. This potential was recently showcased in infrasound studies that illustrated the ability to retrieve Earth's shallow velocities and Mars' near-surface atmospheric winds from remote acoustic observations. Consequently, infrasound is becoming an essential component of planetary missions, providing insights into the interiors and atmospheres of other terrestrial planets in the solar system. However, utilizing infrasound data requires efficient forward wave simulation techniques and an accurate description of waveform sensitivity to source and medium parameters. Even under ideal conditions, many inverse problems associated with infrasound-based probing are inherently ill-posed, necessitating regularization or other approaches to ensure reliable results. In this contribution, we highlight recent seismo-acoustic research findings, addressing atmospheric probing at both smaller and global scales, as well as innovative methods to accelerate infrasonic wave propagation modeling.

Imaging the Garlock Fault Zone With a Fiber: A Limited Damage Zone and Hidden Bimaterial Contrast

AbstractThe structure of fault zones and the ruptures they host are inextricably linked. Fault zones are narrow, which has made imaging their structure at seismogenic depths a persistent problem. Fiber‐optic seismology allows for low‐maintenance, long‐term deployments of dense seismic arrays, which present new opportunities to address this problem. We use a fiber array that crosses the Garlock Fault to explore its structure. With a multifaceted imaging approach, we peel back the shallow structure around the fault to see how the fault changes with depth in the crust. We first generate a shallow velocity model across the fault with a joint inversion of active source and ambient noise data. Subsequently, we investigate the fault at deeper depths using travel‐time observations from local earthquakes. By comparing the shallow velocity model and the earthquake travel‐time observations, we find that the fault's low‐velocity zone below the top few hundred meters is at most unexpectedly narrow, potentially indicating fault zone healing. Using differential travel‐time measurements from earthquake pairs, we resolve a sharp bimaterial contrast at depth that suggests preferred westward rupture directivity.

Isolation management to protect threatened native galaxiid fish species: Lessons from Aotearoa New Zealand

Abstract The use of barriers in freshwater systems to mitigate invasions, known as isolation management, has been increasingly implemented as a protection strategy for native fish. Barriers are crucial for securing vulnerable populations, but much remains to be learnt about the drivers of success and the implications of deliberate fragmentation. In Aotearoa‐New Zealand, henceforth Aotearoa, native non‐migratory galaxiid (NMG) fishes are characterized by a disproportionate number of threatened species with introduced salmonids (trout) identified as a key driver of NMG declines. We conducted a qualitative review of barriers protecting NMG in Aotearoa, illustrating varying outcomes and current knowledge gaps with case studies. Barrier types ranged from natural bedrock waterfalls, shallow velocity chutes and drying reaches to deliberately engineered structures (known commonly as exclusion barriers). Barriers were typically located in headwater streams or spring systems. NMG species secured by fragmentation often have high conservation value and barriers are essential for their protection. Despite the critical role of these barriers, the isolated NMG populations above them are at risk of stochastic extirpation, invasion, or re‐invasion. Exclusion barriers may restrict access to habitat for native migratory fish species and can sometimes alter or remove ideal mesohabitat from protected populations. Overall, when using exclusion barriers, managers must assess extinction risk for individual populations by isolation or invasion, balancing benefits and consequences. However, these assessments are frequently based on incomplete data and hindered by lack of general understanding of the underpinning processes. Despite well‐documented benefits to NMG population survival, key knowledge gaps such as long‐term viability of isolated populations, combinations of physical barrier parameters conferring exclusion under all flow regimes and extent of habitat loss to other native species still exist.

High-Resolution 3D Shallow S-Wave Velocity Structure Revealed by Ambient-Noise Double Beamforming with a Dense Array in Guangzhou Urban Area, China

Abstract Shallow velocity structure surveys are very important for urban seismic hazard monitoring and risk assessment. Ambient-noise tomography provides an ideal way to obtain urban fine structure. In this study, we obtained a high-resolution 3D VS model of the metropolitan areas of the Pearl River Delta (PRD) using the ambient-noise double-beamforming method with a dense nodal array. The new model reveals shallow structures that correlate well with surface geological features, with low-velocity anomalies in fault depressions and high-velocity anomalies in fault uplifts. Our findings reveal detailed fault geometries and basin characteristics of the PRD. The Guangzhou–Conghua fault emerges as a prominent velocity boundary, playing a significant role in controlling the development and subsidence of the Longgui basin. The Xinhui–Shiqiao fault and Shougouling fault are identified as major faults that control the formation and evolution of depressions in the PRD. The basin structures in the PRD are classified as semigraben basins controlled by synsedimentary faults. The long axes of the sub-basins align with the strike of the major faults, and the deposit centers are located in close proximity to these faults. Furthermore, our investigation reveals low-velocity anomalies along the faults, suggesting the existence of pre-existing faults facilitating heat transfer and fluid/melt migration from the deep crust. Our results provide new constraints on the geometric structure of the sedimentary basins and fault systems in the PRD area, thereby contributing to urban seismic hazard assessment and offering valuable insights into potential geothermal resources.

Full Dispersion‐Spectrum Inversion of Surface Waves

AbstractNowadays, the most successful applications of full‐waveform inversion (FWI) involve marine seismic data under acoustic approximations. Elastic FWI of land seismic data is still challenging in theory and practice. Here, we propose a full dispersion spectrum inversion method and apply it to seismic data acquired in West Antarctica. Inspired by the conventional surface wave dispersion curve inversion method, we propose to invert the surface wave dispersion spectrum instead of the complicated waveforms. We compare the frequency‐velocity, frequency‐slowness, and frequency‐wavenumber spectra in terms of their ability to resolve dispersion modes and the feasibility of their adjoint updates and conclude that the frequency‐slowness spectrum is the best for our inversion objectives. We test four objective functions, subtraction, zero‐lag crosscorrelation, optimal transport, and the local‐crosscorrelation to quantify the spectrum mismatch and provide the corresponding adjoint source. We then theoretically analyze the convexity of the proposed objective functions and examine their convergence behavior using numerical examples. We also compare the proposed method with the classic FWI method and the traditional surface wave dispersion curve inversion method and discuss the strengths and weaknesses of each method. This technique is employed to evaluate the shallow velocity structures beneath a seismic array stationed in West Antarctica. Our proposed inversion scheme is also useful for more general applications such as imaging the shallow subsurface of the critical zones, like geothermal reservoirs, and CO2 storage sites.

Examination of shallow and deep S-wave velocity structures from microtremor array measurements and receiver function analysis at strong-motion stations in Kathmandu basin, Nepal

The Himalayan collision zone, where the Indian Plate subducts beneath the Eurasian Plate at a low angle, has caused many devastating earthquakes. The Kathmandu basin, situated in this region, is surrounded by mountains on all sides and is filled with distinct soft lake sediments with a highly undulating bedrock topography. The basin has been experiencing rapid urbanization, and the growing population in its major cities has increased the vulnerability to seismic risk during future earthquakes. Several strong-motion stations have recently been deployed in the Kathmandu basin. It is expected that the data captured by this strong-motion station array will further enhance our understanding of site amplification in sedimentary basins. Clear P-to-S converted waves have been observed in the strong-motion records. In this study, we investigate the medium boundary that generated these converted waves. First, we estimate the shallow velocity structures, which correspond to the topographic slopes or surface geology, beneath the strong-motion stations. We then apply a receiver function analysis to the strong-motion records. The receiver function indicates that the interface between the soft sediment and seismic bedrock serves as a boundary that generates converted waves. The obtained results can be used for tuning three-dimensional velocity structures.Graphical

Imaging the shallow velocity structure of the slow‐spreading ridge of the South China Sea with downward continued multichannel seismic data

AbstractHigh‐resolution shallow oceanic crust velocity models provide crucial information on the tectonothermal history of the oceanic crust. The ocean bottom seismometers record wide‐angle seismic reflection and refraction data to image deeper structures compared with streamer data set. However, most ocean bottom seismometers experiments produce low‐resolution velocity models with limited shallow crustal structure due to sparse ocean bottom seismometers spacing. Multichannel seismic data recorded by towed streamers provide complementary seismic images of the oceanic crust but yield little information on subseafloor velocity because most subseafloor refractions are masked by seafloor reflections. Therefore, it is difficult to obtain fine‐scale velocity structure of shallow upper oceanic crust with both ocean bottom seismometers and multichannel seismic data. Downward continuation technique redatumed the shots and receivers to the seafloor to collapse the seafloor reflections to the zero offset and extract refractions as first arrivals from nearly zero offset, enabling dense ray coverage at the shallow crust. We applied the downward continuation and traveltime tomography methods to two synthetic models, Marmousi and SEAM Phase I Salt models, to demonstrate the performance of the strategy in the situations of flat seafloor and rough seafloor topography. We conducted the first‐arrival traveltime tomography on downward continued towed‐streamer multichannel seismic data across a slow‐spreading ridge of the South China Sea, providing unprecedented details of shallow velocity structure in the sediments. The low velocity sediments revealed by traveltime tomography match well with the prestack depth migration profile.

Shear-Wave Velocity Model for the Dead Sea Transform from Multimode Inversion of Surface Waves Excited by the February 2023 Southeast Türkiye Earthquake Sequence

Abstract The dispersive nature of surface waves can be used for shear-wave velocity inversion at different scales. We show that four large earthquakes from the 2023 Türkiye earthquake sequence generate visible surface waves recorded by a dense strong-motion network deployed along the Dead Sea Transform (DST) in Israel. Thanks to favorable geometrical conditions and source radiation patterns, we observe both Rayleigh and Love waves that travel predominantly parallel to the network. We can reliably compute the dispersion of three Love-wave modes and two Rayleigh-wave modes. Using these dispersion curves, we invert for a 1D S-wave velocity model of the entire DST, outperforming an existing model. Statistical and kernel sensitivity analysis show high certainty down to a depth of 30 km thanks to the multimode joint inversion in the frequency band of 0.03–0.3 Hz. Using a multiwindow approach, we invert for an along-strike laterally varying velocity model of the DST. Although it is limited to the fundamental Love-wave mode, using the 1D model as a constraint allows us to recover a shallow (10 km) velocity structure in agreement with previous studies of the area. Despite the simplicity of our used approach, it can be used as a basis for more advanced studies.

Constraints and insights into the tectonics of the Bayan Obo deposit from the shallow velocity structure

The Bayan Obo deposit in Inner Mongolia is a crucial source of rare earth elements-niobium-iron, with significant economic importance for China. The ore-bearing dolomite is believed to have originated from mantle-derived magma, indicating high prospectivity in the deeper sections. However, the deep distribution of the ore-bearing dolomite remains poorly understood. Here, we develop a 3D S-wave velocity ([Formula: see text]) structure down to 2.5 km depth beneath the Bayan Obo using ambient noise tomography techniques, which used 25 days of continuous waveform data recorded by 312 short-period seismometers. Combined with previous geological, geochemical, and geophysical studies, our results reveal the possible 3D spatial distribution of the ore-bearing dolomite beneath the Bayan Obo. The results indicate that the [Formula: see text] structures beneath the main, east, and west pits change from relatively high-velocity (high-V) zones (with [Formula: see text] of approximately 2.8–3.2 km/s) at shallow depths (depth to approximately 1.5–1.8 km) to significant low-velocity (low-V) zones (with [Formula: see text] < 2.6 km/s) at deeper depths (depth below approximately 1.8 km). We suggest that the content of ore-bearing dolomite may decrease significantly at deeper depths (depth below approximately 1.8 km) beneath these pits, which possibly indicates a diminished prospecting potential there. Our results also indicate a significant high-V zone (with [Formula: see text] > 3.4 km/s) at depths of approximately 1–2.5 km between the main and west pits. Its lithology remains unclear and is proposed to be a key scientific issue in the future to improve our understanding of the ore-controlling structure associated with the tectonic events in the Bayan Obo deposit.

Array analysis of seismic noise at the Sos Enattos mine, the Italian candidate site for the Einstein Telescope

The area surrounding the dismissed mine of Sos Enattos (Sardinia, Italy) is the Italian candidate site for hosting Einstein Telescope (ET), the third-generation gravitational wave (GW) observatory. One of the goals of ET is to extend the sensitivity down to frequencies well below those currently achieved by GW detectors, i.e. down to 2 Hz. In the bandwidth [1,10] Hz, the seismic noise of anthropogenic origin is expected to represent the major perturbation to the operation of the infrastructure, and the site that will host the future detector must fulfill stringent requirements on seismic disturbances. In this paper we describe the operation of a temporary, 15-element, seismic array deployed in close proximity to the mine. Signals of anthropogenic origin have a transient nature, and their spectra are characterized by a wide spectral lobe spanning the [3,20] Hz frequency interval. Superimposed to this wide lobe are narrow spectral peaks within the [3,8] Hz frequency range. Results from slowness analyses suggest that the origin of these peaks is related to vehicle traffic along the main road running east of the mine. Exploiting the correlation properties of seismic noise, we derive a dispersion curve for Rayleigh waves, which is then inverted for a shallow velocity structure down to depths of approx 150 m. This data, which is consistent with that derived from analysis of a quarry blast, provide a first assessment of the elastic properties of the rock materials at the site candidate to hosting ET.

Interplay Between Seismic and Aseismic Deformation on the Central Range Fault During the 2013 Mw 6.3 Ruisui Earthquake (Taiwan)

AbstractThe 2013 Ruisui earthquake represents the first unequivocal evidence of the activity of the Central Range fault (CRF) in central Longitudinal Valley, Taiwan. Using a joint Bayesian finite‐fault source inversion of Global Navigation Satellite System and strain time series, we infer that coseismic rupture occurred between 4 and 19 km depth with maximum slip of 0.5 m located near the hypocenter. We then apply a variational Bayesian independent component analysis approach to displacement signals to infer a 3‐months long afterslip located in the near‐source region. This observation represents the first evidence of aseismic slip on the CRF. Combining geodetic and seismological analysis with simulations based on rate‐and‐state friction mechanics, we analyze the interplay between seismic and aseismic deformation during the earthquake sequence. We observe that afterslip is the dominant postseismic deformation mechanism, with >95% of the moment being released aseismically in the postseismic phase and also likely represents the driving force controlling aftershock productivity. Finally, we infer the presence of a shallow velocity strengthening zone (∼0–4 km depth) associated with spatially heterogeneous slip during the postseismic phase with maximum slip of 0.18 m located above the zone of maximum coseismic deformation.

Along-Strike Variation in the Shallow Velocity Structure beneath the Chenghai Fault Zone, Yunnan, China, Constrained from Methane Sources and Dense Arrays

Abstract The fine structure of the fault zone and the surrounding area is the basis for understanding the process of earthquake nucleation and rupture propagation. To obtain the high-resolution structure of the Chenghai fault (CHF) and the nearby basins, we deployed two dense arrays and excited eight methane sources across the fault from October to November 2020. Based on the 611 P-wave travel times, we obtained the shallow velocity structure beneath the arrays using the simul2000 travel-time inversion program, and the results are as follows: (1) the shallow velocity structure beneath the CHF is very complex, with obvious velocity contrasts on both the sides of the regional fault; (2) low-velocity zones (LVZs) beneath the CHF show clear along-strike variations. The LVZs extend to ∼500 m in depth with widths of ∼2 km and ∼5 km beneath the Qina and Pianjiao arrays, respectively, which are consistent with the Quaternary sediments, and the velocity contrasts along the interface of the LVZ can reach 20%–50%; and (3) the distribution of shallow surface tectonic geomorphology is mainly controlled by regional fault activities that are formed under the combined action of regional near-east–west stretching and clockwise rotation of microblocks. Our results can help improve cognition and seismic hazard assessment for potential earthquakes on the CHF, as well as lay the foundation for understanding the seismic wave velocity variation mechanism in the fault zone.

Surface wave dispersion curve inversion using mixture density networks

SUMMARYIn many seismological, environmental and engineering applications a detailed S-wave velocity model of the shallow subsurface is required. This is generally achieved by the inversion of surface wave dispersion curves using various inversion methods. The classical inversion approaches suffer from several shortcomings, such as inaccurate solutions due to local minima or large computation times in case of a wide parameter space. A number of machine learning (ML) approaches have been suggested to tackle these problems, which however do not provide probabilistic solutions and/or constrain layer number and layer thickness to a fixed value. In this study, we develop a novel neural network (NN) approach in order to characterize the shallow velocity structure from Love and Rayleigh wave dispersion curves. The novelty of our method lies in the simultaneous estimation of layer numbers, layer depth and a complete probability distribution of the S-wave velocity structure. This is achieved by a two-step ML approach, where (1) a regular NN classifies the number of layers within the upper 100 m of the subsurface and (2) a mixture density network outputs the depth estimates together with a fully probabilistic solution of the S-wave velocity structure. We show the advantages of our ML approach compared to a conventional neighbourhood inversion and a Markov chain Monte Carlo algorithm. Our ML approach is then applied to dispersion curves extracted from recorded noise data in Munich, Germany. The resulting velocity profile is in accordance with lithologic information at the site, which highlights the potential of our approach.

The 3-D shallow velocity structure and sedimentary structure of 2017 Ms6.6 Jinghe earthquake source area derived from dense array observations of ambient noise

To understand the shallow structure and complex sedimentary environment of the 2017 Jinghe Ms6.6 earthquake focal area, we used one mouth of continuous sesimic data from a dense seismic array of 208 short period stations around the earthquake focal area, and applied the ambient noise tomography (ANT) method to image the three-dimensional Shear wave velocity structure at the depth less than 4 km. The shear wave velocity shown clear lateral variations and vertical variations from the surface to the deeper regions and has a tight correlation with surface geological and tectonic features in the study area. Obvious low-velocity anomalies has been presented throughout most of the Jinghe depression, whereas the Borokonu Mountains presented high-velocity anomalies. The thickness of the Cenozoic sedimentary basement in the study area is approximately 1–4 km, and the distribution of thickness is highly uneven. The crystalline basement in the study area has strong bending deformation, and the non-uniform Cenozoic sediments are related to the strong bending deformation of the crystalline basement. The Kusongmuchik piedmont fault is a high-angle thrust fault cutting through the base. There are also many medium low-angle faults, which do not penetrate the surface, which has indicated that they are in a concealed state at present. The results have provided a shallow high-resolution 3D velocity model that can be used in the simulation of strong ground motion and for evaluating potential seismic hazards.

The Magnitude 4.5 Earthquake in the Karst Area of Guizhou on 21 August 2021: An Extremely Shallow Earthquake with a Closing-Crack Source

AbstractGuizhou province is located in a karst development area in China, and there are many moderate earthquakes, among which extremely shallow earthquakes are likely to cause damage to buildings, transportation lines, and other engineering facilities. Accurate focal depth and moment tensor can provide vital information for earthquake disaster assessment and karst collapse monitoring in karst areas. The network in the study area is relatively sparse, and the velocity structure is complex. For shallow earthquakes, the shallow velocity structure may cause certain errors in the inversion of the focal depth and mechanisms. In this study, based on the surface-wave information of regional stations, the shallow velocity structure is inverted. Based on the updated velocity model, we obtain the focal depth and mechanism of the Mw 4.5 earthquake on 21 August 2021 in Guizhou province. The results show that the optimal focal depth is about 2 km, and the moment of horizontal closing crack (M​crack) accounts for about 31% of the full moment tensor. The first motions of the P wave at the near-seismic stations of the earthquake are negative, similar to the typical collapses. However, compared to mining collapse earthquakes, the proportion of closing-crake components is relatively low. The dominant double-couple source of the earthquake is a thrust fault, which is consistent with the distribution of faults in the study area. It is speculated that a collapse in the karst area, resulting in negative P-wave polarities and about 31% closing-crake source, induced the shallow earthquake.

Resolving minute temporal seismic velocity changes induced by earthquake damage: the more stations, the merrier?

SUMMARYGround shaking induced by earthquakes often introduces transient changes in seismic velocity monitored with ambient noise. These changes are usually attributed to relaxation behaviour following the coseismic damage in the subsurface and are of relevance for post-seismic hazard mitigation. However, the velocity evolution associated with this phenomenon can occur at very small timescales and amplitudes that are not resolved with seismic interferometry and are therefore challenging to link to laboratory experiments. A way to improve the temporal resolution of the velocity time-series is to test whether the estimation of the relative seismic velocity changes dv/v obeys the ergodic hypothesis in which the joint use of colocated stations would lead to better resolved measurements. In this study, we present results from a dense seismic array that was deployed for 2 weeks at the remarkable Patache site in Chile. Thanks to high temporal averaging capabilities, we are able to resolve seismic velocity changes in the 3–6 Hz frequency band at a 10-min resolution around the occurrence of a moderate earthquake (PGV ∼1 cm s–1). We report a velocity drop of ∼0.4 per cent in the first 10 min after ground shaking. Half of this initial drop was recovered within the 2 following days. The shape of the recovery follows a log-linear shape over the whole observed recovery phase, analogous to slow dynamics experiments. When normalized by the total amount of processed data, we show that the ergodic hypothesis almost perfectly holds in our network: the dv/v signal-to-noise ratio (SNR) obtained when averaging a few observation with large stacking durations for the correlation functions is almost equal to the SNR when using a large number of observations with small stacking durations. To understand if the ergodicity is linked to a particular site property, we use the array capabilities to identify the surf at the shoreline as the source of the noise and to derive a 1-D shear velocity profile with the focal spot imaging technique and a transdimensional Bayesian inversion framework. The inversion shows that hard rocks lie close to the surface indicating that this material hosts the observed shallow velocity changes. We discuss our high-resolution measurements and attribute them to a stable noise source excited by the shore, the ergodicity property and an ideal subsurface structure. Finally, we discuss the effect of moderate earthquakes on subsurface damage and the potential relaxation processes in hard rocks.

Characterizing the shallow structure with the multimodal dispersion curves and the body wave refraction traveltimes from deep seismic sounding data

Deep seismic sounding (DSS) profiles are one of the most powerful tools for detecting crustal structures, and they have been deployed worldwide. Generally, the analysis of DSS data mainly focuses on body waves, while the surface waves are considered noise. We suggest that the surface waves in DSS data can be used to constrain subsurface structures. In this study, we use a DSS profile in the Piedmont and Atlantic Coastal Plain as an example to present the usage of the DSS surface wave. Multimodal dispersion curves were extracted from the DSS data with the Frequency-Bessel transform method, and were used in Monte Carlo joint inversions with body wave refraction traveltimes to constrain the shallow structures. Through the inversion, a horizontal stratum on the surface was identified in the Piedmont, and a two-layer sedimentary structure was identified in the Atlantic Coastal Plain. Comparisons with existing studies verified the accuracy of the shallow structures obtained in this study, demonstrating that the shallow velocity structure could be well constrained with the additional constraints provided by the multimodal dispersion curves. Thus, we believe that further research on the surface waves recorded in DSS surveys is warranted.

Estimation of Shallow Shear Velocity Structure in a Site with Weak Interlayer Based on Microtremor Array

A site condition survey is extremely important for the seismic fortification of major projects. The distribution of underlying weak interlayer in sites is extremely harmful to buildings. However, it is a technical problem to find out the distribution of weak interlayer in the overburden. The shallow velocity structure can directly reflect the change characteristics of a stratigraphic structure. In this paper, acquisition of background noise is conducted using a microtremor linear array method, and the distribution characteristics of two typical stratigraphic structures in Wuhan, Hubei Province, are obtained through an inversion of the apparent S-wave velocity; meanwhile, the equivalent shear-wave velocity and the overburden thickness are estimated, which provides a basis for site classification. The research results are as follows: (1) The two-dimensional profile of the apparent S-wave velocity obtained by the microtremor linear array method can be used for fine imaging of the stratum with weak interlayer, and its distribution form and velocity structure characteristics are highly consistent with those of the drilling data. (2) Compared to the borehole data obtained through in situ test, the error of the overburden thickness and the equivalent shear-wave velocity estimated by the inversion of the apparent S-wave velocity is only about 10%, and the estimated parameters can be directly used for site classification. These results can provide important parameters for seismic fortification of major projects, and also provide reference for the exploration of unfavorable geological bodies, such as weak interlayer in complex urban areas, in the future, which can have good scientific significance and popularization value.

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

AbstractSubsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity (VS) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface VS variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred VS variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The VS variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

AbstractDistributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.