The use of surface wave testing for near‐surface engineering applications has increased in recent years. Typical surface wave analysis is based on the dispersion of surface waves in one‐dimensional layered models. One‐dimensional models are inappropriate for measurements at sites with appreciable lateral variability, a likely scenario in many engineering applications. Use of such models can subsequently undermine the reliability and accuracy of the surface wave results. Full waveform inversion (FWI) is a high‐resolution imaging technique that is proven to outperform the conventional dispersion‐based analysis of surface waves. Much of the near‐surface literature has focused on full waveform inversion of Rayleigh waves developed by the interaction of primary‐ and vertically polarized shear waves (P‐SV), and the capabilities of surface waves generated by horizontally polarized shear waves (Love waves) in mapping near‐surface spatial variabilities have not been fully explored. In this numerical study, full waveform inversion of Rayleigh and Love waves was performed on two different spatially correlated Gaussian random fields (mean
Corneal stiffness plays a critical role in shaping the cornea with respect to intraocular pressure and physical interventions. However, it remains difficult to measure the mechanical properties noninvasively. Here, we report the first measurement of shear modulus in human corneas in vivo using optical coherence elastography (OCE) based on surface elastic waves. In a pilot study of 12 healthy subjects aged between 25 and 67, the Rayleigh-wave speed was 7.86 ± 0.75 m/s, corresponding to a shear modulus of 72 ± 14 kPa. Our data reveal two unexpected trends: no correlation was found between the wave speed and IOP between 13–18 mmHg, and shear modulus decreases with age (− 0.32 ± 0.17 m/s per decade). We propose that shear stiffness is governed by the interfibrillar matrix, whereas tensile strength is dominated by collagen fibrils. Rayleigh-wave OCE may prove useful for clinical diagnosis, refractive surgeries, and treatment monitoring.
more » « less- NSF-PAR ID:
- 10217991
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 10
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT VS of 200 and 500 m/s) to mimic the natural spatial variability of geologic materials. Each soil structure was produced at a low and high level of stiffness variability. Two sources with different frequency contents, 25 and 50 Hz, were used to evaluate the effects of source characteristics on the resolution of Rayleigh and Love waveform inversions. The inverted results from the high‐velocity domain demonstrated that Love waveform inversion using high‐frequency seismic sources outperforms Rayleigh full waveform inversion in detecting the shape and the velocities of horizontally deposited geologic materials. Results from the low‐velocity domain also confirmed that Love full waveform inversion was comparable or superior to Rayleigh full waveform inversion, though the performance difference was less significant. However, the 25‐Hz frequency inversions yielded superior results than the 50‐Hz frequency inversions for the low‐velocity domain because the dominant wavelength of the high‐frequency signals becomes so small that it offers an impractically small investigation depth. -
more » « less
Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.
-
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 This study presents an azimuthally anisotropic shear wave velocity model of the crust and uppermost mantle beneath Alaska, based on Rayleigh wave phase speed observations from 10 to 80 s period recorded at more than 500 broadband stations. We test the hypothesis that a model composed of two homogeneous layers of anisotropy can explain these measurements. This “Two‐Layer Model” confines azimuthal anisotropy to the brittle upper crust along with the uppermost mantle from the Moho to 200 km depth. This model passes the hypothesis test for most of the region of study, from which we draw two conclusions. (a) The data are consistent with crustal azimuthal anisotropy being dominantly controlled by deformationally aligned cracks and fractures in the upper crust undergoing brittle deformation. (b) The data are also consistent with the uppermost mantle beneath Alaska and surroundings experiencing vertically coherent deformation. The model resolves several prominent features. (1) In the upper crust, fast directions are principally aligned with the orientation of major faults. (2) In the upper mantle, fast directions are aligned with the compressional direction in compressional tectonic domains and with the tensional direction in tensional domains. (3) The mantle fast directions located near the Alaska‐Aleutian subduction zone and the surrounding back‐arc area form a toroidal pattern that is consistent with mantle flow directions predicted by recent geodynamical models. Finally, the mantle anisotropy is remarkably consistent with SKS fast directions, but to fit SKS split times, anisotropy must extend below 200 km depth across most of the study region.
-
more » « less
SUMMARY 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.
