A robust, in situ estimate of shear-wave velocity VS and the small-strain damping ratio DS (or equivalently, the quality factor QS) is crucial for the design of buildings and geotechnical systems subjected to vibrations or earthquake ground shaking. A promising technique for simultaneously obtaining both VS and DS relies on the Multichannel Analysis of Surface Waves (MASW) method. MASW can be used to extract the Rayleigh wave phase velocity and phase attenuation data from active-source seismic traces recorded along linear arrays. Then, these data can be inverted to obtain VS and DS profiles. This paper introduces two novel methodologies for extracting the phase velocity and attenuation data. These new approaches are based on an extension of the beamforming technique which can be combined with a modal filter to isolate different Rayleigh propagation modes. Thus, the techniques return reliable phase velocity and attenuation estimates even in the presence of a multimode wavefield, which is typical of complex stratigraphic conditions. The reliability and effectiveness of the proposed approaches are assessed on a suite of synthetic wavefields and on experimental data collected at the Garner Valley Downhole Array and Mirandola sites. The results reveal that, under proper modelling of wavefield conditions, accurate estimates of Rayleigh wave phase velocity and attenuation can be extracted from active-source MASW wavefields over a broad frequency range. Eventually, the estimation of soil mechanical parameters also requires a robust inversion procedure to map the experimental Rayleigh wave parameters into soil models describing VS and DS with depth. The simultaneous inversion of phase velocity and attenuation data is discussed in detail in the companion paper.
This paper deals with in situ characterization of the small-strain shear-wave velocity VS and damping ratio DS from an advanced interpretation of Multi-channel Analysis of Surface Waves (MASW) surveys. A new approach based on extracting Rayleigh wave data using the CFDBFa method has been discussed in the companion paper. This paper focuses on mapping the experimental Rayleigh wave phase velocity and attenuation into profiles of VS and DS versus depth, which is achieved through a joint inversion procedure. The joint inversion of phase velocity and attenuation data utilizes a newly developed Monte Carlo global search algorithm, which implements a smart sampling procedure. This scheme exploits the scaling properties of the solution of the Rayleigh eigenvalue problem to modify the trial earth models and improve the matching with the experimental data. Thus, a reliable result can be achieved with a limited number of trial ground models. The proposed algorithm is applied to the inversion of synthetic data and of experimental data collected at the Garner Valley Downhole Array site, as described in the companion paper. In general, inverted soil models exhibit well-defined VS profiles, whereas DS profiles are affected by larger uncertainties. Greater uncertainty in the inverted DS profiles is a direct result of higher variability in the experimental attenuation data, the limited wavelength range at which reliable values of attenuation parameters can be retrieved, and the sensitivity of attenuation data to both DS and VS. Nonetheless, the resulting inverted earth models agree well with alternative in situ estimates and geological data. The results stress the feasibility of retrieving both stiffness and attenuation parameters from active-source MASW testing and the effectiveness of extracting in situ damping ratio estimates from surface wave data.
more » « less- NSF-PAR ID:
- 10492932
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 237
- Issue:
- 1
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 525-539
- Size(s):
- ["p. 525-539"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
null (Ed.)SUMMARY We determine crustal shear wave velocity structure and crustal thickness at recently deployed seismic stations across West Antarctica, using a joint inversion of receiver functions and fundamental mode Rayleigh wave phase velocity dispersion. The stations are from both the UK Antarctic Network (UKANET) and Polar Earth Observing Network/Antarctic Network (POLENET/ANET). The former include, for the first time, four stations along the spine of the Antarctic Peninsula, three in the Ellsworth Land and five stations in the vicinity of the Pine Island Rift. Within the West Antarctic Rift System (WARS) we model a crustal thickness range of 18–28 km, and show that the thinnest crust (∼18 km) is in the vicinity of the Byrd Subglacial Basin and Bentley Subglacial Trench. In these regions we also find the highest ratio of fast (Vs = 4.0–4.3 km s–1, likely mafic) lower crust to felsic/intermediate upper crust. The thickest mafic lower crust we model is in Ellsworth Land, a critical area for constraining the eastern limits of the WARS. Although we find thinner crust in this region (∼30 km) than in the neighbouring Antarctic Peninsula and Haag-Ellsworth Whitmore block (HEW), the Ellsworth Land crust has not undergone as much extension as the central WARS. This suggests that the WARS does not link with the Weddell Sea Rift System through Ellsworth Land, and instead has progressed during its formation towards the Bellingshausen and Amundsen Sea Embayments. We also find that the thin WARS crust extends towards the Pine Island Rift, suggesting that the boundary between the WARS and the Thurston Island block lies in this region, ∼200 km north of its previously accepted position. The thickest crust (38–40 km) we model in this study is in the Ellsworth Mountain section of the HEW block. We find thinner crust (30–33 km) in the Whitmore Mountains and Haag Nunatak sectors of the HEW, consistent with the composite nature of the block. In the Antarctic Peninsula we find a crustal thickness range of 30–38 km and a likely dominantly felsic/intermediate crustal composition. By forward modelling high frequency receiver functions we also assess if any thick, low velocity subglacial sediment accumulations are present, and find a 0.1–0.8-km-thick layer at 10 stations within the WARS, Thurston Island and Ellsworth Land. We suggest that these units of subglacial sediment could provide a source region for the soft basal till layers found beneath numerous outlet glaciers, and may act to accelerate ice flow.more » « less
-
more » « less
Abstract Near‐surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near‐surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short‐period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near‐surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation.
-
Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100 km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models.more » « less
-
SUMMARY Measurements of various physical properties of oceanic sediment and crustal structures provide insight into a number of geological and geophysical processes. In particular, knowledge of the shear wave velocity (VS) structure of marine sediments and oceanic crust has wide ranging implications from geotechnical engineering projects to seismic mantle tomography studies. In this study, we propose a novel approach to nonlinearly invert compliance signals recorded by colocated ocean-bottom seismometers and high-sample-rate pressure gauges for shallow oceanic shear wave velocity structure. The inversion method is based on a type of machine learning neural network known as a mixture density neural network (MDN). We demonstrate the effectiveness of the MDN method on synthetic models with a fixed deployment depth of 2015 m and show that among 30 000 test models, the inverted shear wave velocity profiles achieve an average error of 0.025 km s−1. We then apply the method to observed data recorded by a broad-band ocean-bottom station in the Lau basin, for which a VS profile was estimated using Monte Carlo sampling methods. Using the mixture density network approach, we validate the method by showing that our VS profile is in excellent agreement with the previous result. Finally, we argue that the mixture density network approach to compliance inversion is advantageous over other compliance inversion methods because it is faster and allows for standardized measurements.more » « less
