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.
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.
more » « less- NSF-PAR ID:
- 10492933
- 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. 506-524
- Size(s):
- ["p. 506-524"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
Abstract Non-invasive surface wave methods are increasingly being used as the primary technique for estimating a site’s small-strain shear wave velocity (Vs). Yet, in comparison to invasive methods, non-invasive surface wave methods suffer from highly variable standards of practice, with each company/group/analyst estimating surface wave dispersion data, quantifying its uncertainty (or ignoring it in many cases), and performing inversions to obtain Vs profiles in their own unique manner. In response, this work presents a well-documented, production-tested, and easy-to-adopt workflow for developing estimates of experimental surface wave dispersion data with robust measures of uncertainty. This is a key step required for propagating dispersion uncertainty forward into the estimates of Vs derived from inversion. The paper focuses on the two most common applications of surface wave testing: the first, where only active-source testing has been performed, and the second, where both active-source and passive-wavefield testing has been performed. In both cases, clear guidance is provided on the steps to transform experimentally acquired waveforms into estimates of the site’s surface wave dispersion data and quantify its uncertainty. In particular, changes to surface wave data acquisition and processing are shown to affect the resulting experimental dispersion data, thereby highlighting their importance when quantifying uncertainty. In addition, this work is accompanied by an open-source Python package, swprocess , and associated Jupyter workflows to enable the reader to easily adopt the recommendations presented herein. It is hoped that these recommendations will lead to further discussions about developing standards of practice for surface wave data acquisition, processing, and inversion.more » « less
-
more » « less
Abstract The Formosa array, with 137 broadband seismometers and ∼5 km station spacing, was deployed recently in Northern Taiwan. Here by using eight months of continuous ambient noise records, we construct the first high‐resolution three‐dimensional (3‐D) shear wave velocity model of the crust in the area. We first calculate multi‐component cross‐correlations to extract robust Rayleigh wave signals. We then determine phase velocity maps between 3 and 10 s periods using Eikonal tomography and measure Rayleigh wave ellipticity at each station location between 2 and 13 s periods. For each location, we jointly invert the two types of Rayleigh wave measurements with a Bayesian‐based inversion method for a one‐dimensional shear wave velocity model. All piecewise continuous one‐dimensional models are then used to construct the final 3‐D model. Our 3‐D model reveals upper crustal structures that correlate well with surface geological features. Near the surface, the model delineates the low‐velocity Taipei and Ilan Basins from the adjacent fast‐velocity mountainous areas, with basin geometries consistent with the results of previous geophysical exploration and geological studies. At a greater depth, low velocity anomalies are observed associated with the Linkou Tableland, Tatun Volcano Group, and a possible dyke intrusion beneath the Southern Ilan Basin. The model also provides new geometrical constraints on the major active fault systems in the area, which are important to understand the basin formation, orogeny dynamics, and regional seismic hazard. The new 3‐D shear wave velocity model allows a comprehensive investigation of shallow geologic structures in the Northern Taiwan.
-
more » « less
Abstract I present a 3‐D isotropic shear wave velocity model of the crust and uppermost mantle beneath the Alaska‐Aleutian subduction zone offshore of the Alaska Peninsula, based on seismic data recorded by the Alaska Amphibious Community Seismic Experiment (AACSE) array and some other networks. The model derives from Rayleigh wave phase speed measurements extracted from ambient seismic noise. A new three‐station interferometry (Zhang et al., 2020) approach is applied to improve the data coverage of ambient noise surface waves. Based on the ambient noise Rayleigh wave dispersion data, a Bayesian Monte Carlo inversion is performed to produce the shear wave velocity model. There are several prominent structures captured by the model, including: (1) The major sedimentary basins across the study region are identified by model. (2) Crustal thickness estimates are related with the geological structures. (3) The imaged slab edge is consistent with both the Slab 2.0 model (Hayes et al., 2018) and earthquake locations. And lots of geological and tectonic features related to subduction zone are captured, including the serpentinized forearc and partial melting zone beneath the Aleutian arc volcanoes. (4) Near the Shumagin gap, reduction in Vs is observed at the uppermost part of the incoming Pacific plate, consistent with the active source study of Shillington et al. (2015). The Vs reduction reflects hydration of the oceanic plate which could be related to local seismicity variation.
-
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
