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: High‐Resolution 3‐D Shear Wave Velocity Model of Northern Taiwan via Bayesian Joint Inversion of Rayleigh Wave Ellipticity and Phase Velocity With Formosa Array
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
Award ID(s):
1753362
PAR ID:
10361499
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
5
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We construct a new shear velocity model for the San Gabriel, Chino and San Bernardino basins located in the northern Los Angeles area using ambient noise correlation between dense linear nodal arrays, broadband stations, and accelerometers. We observe Rayleigh and Love waves in the correlation of vertical (Z) and transverse (T) components, respectively. By combining Hilbert and Wavelet transforms, we obtain the separated fundamental and first higher mode of the Rayleigh wave dispersion curves based on their distinct particle motion polarization. Basin depths constrained by receiver functions, gravity, and borehole data are incorporated into the prior model. Our 3D shear wave velocity model covers the upper 3–5 km of the crust in the San Gabriel, Chino and San Bernardino basin area. The Vs model is in agreement with the geological and geophysical cross‐sections from other studies, but discrepancies exist between our model and a Southern California Earthquake Center community velocity model. Our shear wave velocity model shows good consistency with the CVMS 4.26 in the San Gabriel basin, but predicts a deeper and slower sedimentary basin in the San Bernardino and Chino basins than the community model. 
    more » « less
  2. Abstract We construct a 3D shear velocity model of the Salt Lake Valley using Rayleigh waves excited by the 31 March 2020 Mw 6.5 central Idaho earthquake recorded on a 168-station temporary nodal geophone network and the 49-station permanent regional network. The temporary array—deployed in response to the March 18 Mw 5.7 Magna earthquake—serendipitously recorded clear surface waves between 10 and 20 s period from the Idaho event at ∼500 km epicentral distance, from which we measure both Rayleigh wave phase velocity and ellipticity (H/V ratio). In addition, we employ multicomponent earthquake coda cross correlation to extend the measurements down to 5 s period. Because Rayleigh wave ellipticity features outstanding shallow sensitivity, we invert for a 3D upper crust VS model of the Salt Lake Valley. Our model shows basin structure in general agreement with and complements the current Community Velocity Model, which is mostly constrained by borehole and gravity measurements. Our model thus provides critical information for future earthquake hazard assessment studies, which require detailed shallow velocity structure. 
    more » « less
  3. SUMMARY Taiwan, one of the most active orogenic belts in the world, undergoes orogenic processes that can be elucidated by the doubly vergent wedge model, explaining the extensive island-wide geological deformation. To provide a clearer depiction of its cross-island orogenic architecture, we apply ambient noise tomography across an east–west linear seismic array in central Taiwan, constructing the first high-resolution 2-D shear velocity model of the upper crust in the region. We observe robust fundamental- and higher-mode Rayleigh waves, with the latter being mainly present in the western Coastal Plain. We develop a multimode double-beamforming method to determine local phase velocities across the array between 2- and 5-s periods. For each location, we jointly invert all available fundamental- and higher-mode phase velocities using a Bayesian-based inversion method to obtain a 1-D model. All 1-D models are then combined to form a final 2-D model from the surface to ∼10 km depth. Our newly developed 2-D model clearly delineates major structural boundaries and fault geometries across central Taiwan, thereby corroborating the previously proposed pro-wedge and retro-wedge models while offering insight into regional seismic hazards. 
    more » « less
  4. Abstract We apply ambient noise tomography to a seismic array from the Trans‐Haiti project to obtain a 2‐D shear wave velocity (Vs) across Haiti. We perform multi‐component noise cross‐correlation, measure Rayleigh wave phase velocity and its horizontal‐to‐vertical amplitude ratio (H/V) between periods of 3–18 s, and jointly invert both measurements into Vs for the crustal structures of Haiti. Both H/V and phase velocity measurements exhibit consistent patterns related to the geologic units. Sedimentary basins—CSE and Plateau Central basins—show higher H/V values, while mountain areas—Massif de la Selle, Chaine des Matheux, Montagnes Noires and Massif de Nord—exhibit lower H/V. Regarding phase velocity, higher velocities are observed in northern and southern Haiti, likely reflecting the thinner crust compared to the thicker crust showing lower velocities in the central part. While our Vs model is consistent with previous model that suggested thinner crustal thickness in the northern and southern Haiti, with thickening in the center, the Moho interface in the central domain might be shallower than previously thought. 
    more » « less
  5. SUMMARY Long-period (T > 10 s) shear wave reflections between the surface and reflecting boundaries below seismic stations are useful for studying phase transitions in the mantle transition zone (MTZ) but shear-velocity heterogeneity and finite-frequency effects complicate the interpretation of waveform stacks. We follow up on a recent study by Shearer & Buehler (hereafter SB19) of the top-side shear wave reflection Ssds as a probe for mapping the depths of the 410-km and 660-km discontinuities beneath the USArray. Like SB19, we observe that the recorded Ss410s-S and Ss660s-S traveltime differences are longer at stations in the western United States than in the central-eastern United States. The 410-km and 660-km discontinuities are about 40–50 km deeper beneath the western United States than the central-eastern United States if Ss410s-S and Ss660s-S traveltime differences are transformed to depth using a common-reflection point (CRP) mapping approach based on a 1-D seismic model (PREM in our case). However, the east-to-west deepening of the MTZ disappears in the CRP image if we account for 3-D shear wave velocity variations in the mantle according to global tomography. In addition, from spectral-element method synthetics, we find that ray theory overpredicts the traveltime delays of the reverberations. Undulations of the 410-km and 660-km discontinuities are underestimated when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ. Therefore, modelling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes. 
    more » « less