skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Estimation of Vs30 at the EarthScope Transportable Array Stations by Inversion of Low‐Frequency Seismic Noise
Abstract One of the main sources of seismic noise below 0.05 Hz is the atmospheric pressure variation, especially when surface pressure variations are large. When surface broadband seismic stations are equipped with pressure sensors, there is high coherence between pressure and seismic signals at low frequencies. The amount of ground deformation under surface pressure variations reflects the characteristics of near‐surface elastic structure and allows us to estimate near‐surface shear‐modulus structure using an inversion method. In the inversion method, we have the surface observable η(f) = Sz/Sp, where f is a frequency between 0.01 and 0.05 Hz, and Sz and Sp are the power spectral densities of vertical seismic data and of surface pressure data. We derive depth sensitivity kernels for η(f) with which we invert for elastic moduli of the shallow structure. Between 0.01 and 0.05 Hz, sensitivity kernels typically have peaks at depths within the uppermost 100 m. Based on vertically heterogeneous 1‐D structures, we estimate Vs30 at 744 USArray Transportable Array stations. Vs30 is the time‐averaged shear‐wave velocity from the surface to the 30‐m depth. We compare our results with various surficial geology maps. Although Vs30 has high horizontal variability over a short distance on the scale of hundreds of meters, we find correlations between Vs30 and large‐scale geological structures, such as mapped units and surficial materials. We find good agreement between estimated Vs30 and mapped Quaternary sediment depths, where stations with thicker underlying sediment tend to have slower Vs30.  more » « less
Award ID(s):
2225216
PAR ID:
10610304
Author(s) / Creator(s):
;
Publisher / Repository:
Journal of Geophysical Research: Solid Earth
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
127
Issue:
4
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Vertical records of ocean-bottom seismographs (OBSs) are usually noisy at low frequencies, and one important noise source is the varying ocean-bottom pressure that results from ocean-surface water waves. The relation between the ocean-bottom pressure and the vertical seafloor motion, called the compliance pressure transfer function (PTF), can be derived using background seismic data. During an earthquake, earthquake signals also generate ocean-bottom pressure fluctuations, and the relation between the ocean-bottom pressure and the vertical seafloor motion is named the seismic PTF in this article. Conventionally, we use the whole pressure records and the compliance PTF to remove the compliance noise; the earthquake-induced pressure and the seismic PTF are ignored, which may distort the original signals. In this article, we analyze the data from 24 OBSs with water depth ranging from 107 to 4462 m. We find that for most stations, the investigated frequency range (0.01–0.2 Hz) can be divided into four bands depending on the water depth. In band (I) of lowest frequencies (<0.11, <0.05, and <0.02  Hz for water depth of 107, 1109, and 2650 m, respectively), the vertical seafloor acceleration is composed mostly of pressure compliance noise, which can be removed using the compliance PTF. The compliance PTF is much smaller than the seismic PTF, so distortion of earthquake signals is negligible. In band (II) of higher frequencies (0.11–0.20, 0.05–0.11, and 0.02–0.05 Hz for water depth of 107, 1109, and 2650 m, respectively), the vertical acceleration and ocean-bottom pressure are largely uncorrelated. In bands (III) and (IV) of even higher frequencies (>0.11 and >0.08  Hz for water depth of 1109 and 2650 m, respectively), the compliance noise is negligible, and the ocean-bottom pressure is mostly caused by the seafloor motion. Thus, the compliance can be safely ignored in frequency band (I). 
    more » « less
  2. SUMMARY The near-surface seismic structure (to a depth of about 1000 m), particularly the shear wave velocity (VS), can strongly affect the propagation of seismic waves and, therefore, must be accurately calibrated for ground motion simulations and seismic hazard assessment. The VS of the top (<300 m) crust is often well characterized from borehole studies, geotechnical measurements, and water and oil wells, while the velocities of the material deeper than about 1000 m are typically determined by tomography studies. However, in depth ranges lacking information on shallow lithological stratification, typically rock sites outside the sedimentary basins, the material parameters between these two regions are typically poorly characterized due to resolution limits of seismic tomography. When the alluded geological constraints are not available, models, such as the Southern California Earthquake Center (SCEC) Community Velocity Models (CVMs), default to regional tomographic estimates that do not resolve the uppermost VS values, and therefore deliver unrealistically high shallow VS estimates. The SCEC Unified Community Velocity Model (UCVM) software includes a method to incorporate the near-surface earth structure by applying a generic overlay based on measurements of time-averaged VS in top 30 m (VS30) to taper the upper part of the model to merge with tomography at a depth of 350 m, which can be applied to any of the velocity models accessible through UCVM. However, our 3-D simulations of the 2014 Mw 5.1 La Habra earthquake in the Los Angeles area using the CVM-S4.26.M01 model significantly underpredict low-frequency (<1 Hz) ground motions at sites where the material properties in the top 350 m are significantly modified by the generic overlay (‘taper’). On the other hand, extending the VS30-based taper of the shallow velocities down to a depth of about 1000 m improves the fit between our synthetics and seismic data at those sites, without compromising the fit at well-constrained sites. We explore various tapering depths, demonstrating increasing amplification as the tapering depth increases, and the model with 1000 m tapering depth yields overall favourable results. Effects of varying anelastic attenuation are small compared to effects of velocity tapering and do not significantly bias the estimated tapering depth. Although a uniform tapering depth is adopted in the models, we observe some spatial variabilities that may further improve our method. 
    more » « less
  3. Abstract Records of pressure variations on seismographs were historically considered unwanted noise; however, increased deployments of collocated seismic and acoustic instrumentation have driven recent efforts to use this effect induced by both wind and anthropogenic explosions to invert for near-surface Earth structure. These studies have been limited to shallow structure because the pressure signals have relatively short wavelengths (<∼300 m). However, the 2022 eruption of Hunga Tonga–Hunga Ha’apai (also called “Hunga”) volcano in Tonga generated rare, globally observed, high-amplitude infrasound signals with acoustic wavelengths of tens of kilometers. In this study, we examine the acoustic-to-seismic coupling generated by the Hunga eruption across 82 Global Seismographic Network (GSN) stations and show that ground motion amplitudes are related to upper (0 to ∼5 km) crust material properties. We find high (>0.8) correlations between pressure and vertical component ground motion at 83% of the stations, but only 30% of stations show this on the radial component, likely due to complex tilt effects. We use average elastic properties in the upper 5.2 km from the CRUST1.0 model to estimate vertical seismic/acoustic coupling coefficients (SV/A) across the GSN network and compare these to recorded observations. We exclude many island stations from these comparisons because the 1° resolution of the CRUST1.0 model places a water layer below these stations. Our simple modeling can predict observed SV/A within a factor of 2 for 94% of the 51 non-island GSN stations with high correlations between pressure and ground motion. These results indicate that analysis of acoustic-to-seismic coupling from the eruption could be used to place additional constraints on crustal structure models at stations with collocated seismic and pressure sensors. Ultimately, this could improve tomographic imaging models, which rely on methods that are sensitive to local structure. 
    more » « less
  4. Abstract. Large subduction earthquakes can rupture the shallow part of the megathrust with unusually large displacements and tsunamis. The long duration of the seismic source and high upper-plate compliance contribute to large and protracted long-period motions of the outer upper plate. The resulting shear stress at the sediment–water interface in, for example, the Mw 9.0 2011 Tohoku–Oki earthquake could account for surficial sediment remobilization on the outer margin. We test this hypothesis by simulating in physical tank experiments the combined effects of high- and low-frequency seismic motions on sediment of different properties (chemistry, grain size, water content, and salinity). Our results show that low-frequency motion during a 2011-like earthquake can entrain several centimeters of surficial sediment and that entrainment can be enhanced by high-frequency vertical oscillations. These experiments validate a new mechanism of co-seismic sediment entrainment in deep-water environments. 
    more » « less
  5. Abstract We present observations and modeling of spatial eigen‐functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5‐Hz geophones that cross the SJFZ with ~10–30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after theSwave arrival are consistent for more than 50 near‐fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio show three main frequency peaks at ~1.3, ~2.0, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen‐functions of resonance waves in a low‐velocity layer (fault zone) between two quarter‐spaces. Using a grid search‐based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65%Swave velocity reduction compared to the host rock. The SW edge of the low‐velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together. 
    more » « less