skip to main content


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
NSF-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 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
  2. 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
  3. SUMMARY

    Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

     
    more » « less
  4. Abstract

    To investigate the controls on continental rifting in the western branch of the East Africa Rift System, we conduct shear velocity imaging of the crust and uppermost mantle beneath the weakly extended Malawi Rift and the Rungwe Volcanic Province (RVP). We use local‐scale measurements of Rayleigh wave phase velocities between 9‐ and 100‐s periods combined with constraints on basin architecture and crustal thickness to invert for shear velocity from the surface to ~135 km. Our resulting 3‐D model reveals a localized low‐velocity anomaly associated with the RVP extending from the crust and through the upper mantle, which can be explained with modestly elevated temperatures. Away from the RVP, velocities within mantle flanking the rift are fast (>4.6 ± 0.1 km/s), suggesting depleted lithospheric mantle to depths of ~100 and >135 km to the west and east of the rift, respectively. The upper mantle beneath the rift axis is characterized by thinned lithosphere with slower velocities than the surrounding plateau, suggestive of thermal and/or chemical modification by the rifting process. Slowest velocities are mildly asymmetric about the rift axis, with the lowest velocities observed beneath the rift and adjacent footwall escarpments. The underlying asthenosphere is only moderately slow (~4.25 ± 0.1 km/s), including beneath the RVP, precluding the presence of significant volumes of partial melt. The positions of localized lithospheric modification and basin‐bounding border faults correlate with the location of Proterozoic mobile belts, suggesting that these sutures provide lithospheric‐scale weakening mechanisms necessary for localizing strain and allowing extension to occur in the Malawi Rift.

     
    more » « less
  5. Abstract

    The Shanxi Rift is one of the world's largest Cenozoic continental rifts. It is one of the most earthquake prone areas in China and has volcanic activities characterized by small volumes of magma. In order to study continental rifting mechanism, intraplate volcanism, and seismicity of the Shanxi Rift, we construct a high‐resolution 3‐D lithospheric shear velocity (Vs) model with a focus of the crust by jointly inverting Rayleigh wave dispersion curves and ellipticity measurements. Compared with published models, our 3‐D Vs model has improved the resolution for shallow structures. Our model reveals several significant velocity features. (1) The sedimentary layer of the Datong Basin is much thicker than those in the southern Shanxi Rift, implying that the sedimentation rate of the Datong Basin is higher, which could be attributed in part to asthenosphere upwelling beneath the Datong Basin area in the past few Ma in addition to passive rifting. (2) A low‐velocity body is observed in the middle crust beneath the Datong Volcanoes, suggesting the presence of a small amount of partial melting in the middle crust. (3) Lower velocities in the middle and lower crust beneath the Datong Basin are coincident with the shallowing of deep earthquakes from the southern to northern Shanxi Rift in the crust, suggesting that higher temperature may elevate the brittle‐to‐ductile transition zone, leading to the shallowing of focal depths of earthquakes in the Shanxi Rift from south to north.

     
    more » « less