Joint inversion of Rayleigh wave ellipticity and phase velocity for crustal structure in Taiwan
Introduction
Due to the convergence between the Eurasian Plate (EUP) and the Philippine Sea Plate (PSP), the Taiwan Orogen is experiencing violent deformation as manifested by the significant seismicity in and around its surroundings. East to the Taiwan Orogen, the PSP is subducting northwestward beneath the EUP along the Ryukyu trench, whereas to the south, the PSP is overriding the EUP along the Manila trench. Overall, the PSP is advancing towards the EUP at ~90 mm/yr (Sella et al., 2002). Acting as one of the youngest orogens of the world, the Taiwan Orogen has attracted much attention from geoscientists, especially since the deployments of the broadband seismological array and GPS network in this region (e.g., Roecker et al., 1987 Wu et al., 2007 Kuo-Chen et al., 2012 Huang et al., 2014a, Huang et al., 2014b Huang, 2017).
With regard to the slab interactions near Taiwan, a series of geophysical observations advocate southwestward propagating slab tearing along the Chinese continental margin, i.e. following the geometry of the continental-ocean boundary (e.g., Lallemand et al., 2001 Teng et al., 2000 Ustaszewski et al., 2012). Specifically, Teng et al. (2000) inferred that the Eurasian slab has broken off all the way from under the southern Ryukyu to north Taiwan. Correspondingly, the north Taiwan has undergone a post-collisional lithospheric extension resulted from the Eurasian slab breakoff (Teng, 1996 & Teng et al., 2000; Ustaszewski et al., 2012 Wang et al., 1999), as manifested by the horizontal GPS measurements (Bos et al., 2003). On the other hand, the geodetic leveling indicates that the Taiwan mountain belt uplifted at 0.2–18 mm/yr in the most recent 20 years, and the maximum uplift rates, 22.9 mm/yr and 25.8 mm/yr, are observed at the Central Range and Coastal Range respectively (Ching et al., 2011b). It is therefore expected that the central Taiwan is still in the process of mountain building, in contrast to the northeastern Taiwan suffering the post-collisional lithospheric extension. In this tectonic scenario, the mountain building process in Taiwan is a typical example that begins with the formation of the collisional orogen and ends with subsequent collisional collapse (Teng, 1996).
To understand the mountain building in Taiwan, several conceptual models have been proposed based on the present-day deformation pattern and the observed seismicity distribution. For example, the thin-skinned model suggests that the EUP subducts beneath the PSP with an accretionary wedge overlying a planar decollement fault (Suppe, 1981 Dahlen, 1990). Similarly, the lithospheric collision hypothesis advocates strong crust thickening and shortening occur on both sides of the Longitudinal Valley which is the suture zone that separates the PSP and EUP (Wu et al., 1997). Essentially, these models give only two-dimensional views of the Taiwan orogeny. Recently, new results on slab interactions beneath Taiwan are reported (e.g., Huang et al., 2014a, Huang et al., 2014b Huang et al., 2015 Ai et al., 2019b). Especially, we have provided a crustal and uppermost mantle model beneath the Taiwan Island, Taiwan Strait and southwestern Ryukyu Arc from ambient noise tomography (Ai et al., 2019b). This model, however, has a poor resolution in the crustal part limited to the frequency band of the surface wave dispersions. To decipher how the mountain building in Taiwan is linked to the subsurface slab interactions, a high-resolution crust model is desired. In the present study, we present three 2-D shear wave velocity transects that are located in the northern, central and southern Taiwan respectively (Fig. 1) by joint inversion of teleseismic Rayleigh wave H/V ratios and phase velocities. By introducing Rayleigh wave ellipticity as an additional type of observable into shear wave velocity inversion, we can overcome the low depth-resolution of the shallow structure that obtained from inversion with phase velocity alone (Ai et al., 2019b), as the H/V ratio is more sensitive to the shallow structure than the phase velocity at the same period (Boore and Nafi Toksöz, 1969; Lin et al., 2012 Lin and Schmandt, 2014 & Lin et al., 2014).
Similar seismological studies aimed at resolving the arc-continent collision structure were also conducted along the three transects in which we apply the joint inversion (Van Avendonk et al., 2014 & Van Avendonk et al., 2015; Kuo et al., 2016). Compared to the P wave tomography studies, the present joint inversions constrain the crust in western Taiwan well, which is dominated by relatively lower seismicity, though the lateral resolution in our results may be limited by the inter-station distance. Overall, the shear wave velocity images obtained here are complementary to the P wave velocity profiles and the 3-D lithosphere-scale models and provide more details on the evolution of the Taiwan Orogen associated with the slab interactions.
Section snippets
The TAIGER data
After the occurrence of the 1999 Chi-Chi earthquake, the TAiwan Integrated GEodynamic Research (TAIGER) project (2004–2009) was carried out to better understand the seismogenic environment in Taiwan (Wu et al., 2014). One of the main objectives of the TAIGER is to refine the tomographic images of the lithosphere by extending the aperture of the existing broadband networks (Kuo-Chen et al., 2012). In this study, to measure the teleseismic Rayleigh wave H/V ratios, we use the land-based broadband
Rayleigh wave H/V ratios
H/V maps at selected periods from 12 to 80 s are presented in Fig. 7. In general, the H/V measurements in western Taiwan are obviously greater than that observed in the east, especially at short periods. At ~12 s, in central segment of the Coastal Plain (Fig. 7a), the H/V values are as high as 3.0, which is much higher than that observed at the rest parts of Taiwan. At intermediate periods (Fig. 7c & d), the highest H/V ratios appear in the Longitudinal Valley. In the longer period band (Fig. 7
Deflection of Rayleigh wave packet from the great-circle path at the continental margin
The radial components of Rayleigh wave packets used for H/V measurement in this study are obtained by rotating the N′ and E' components by α degrees. As demonstrated in Fig. 2, both the non-great-circle effect (deflection of Rayleigh wave packet from the great-circle path) ϕ and the sensor misorientation θ contribute to α (α = β + ϕ – θ). We thus cannot assess the non-great-circle effects in Taiwan quantitatively without knowledge on the misorientations of the stations. Here, we define ϕ – θ as
Conclusion
In this study, broadband Rayleigh wave H/V measurements are obtained from teleseismic waveforms collected from the TAIGER project. By performing joint inversion of the H/V measurements and phase velocities obtained previously, three shear wave velocity transects across the Taiwan Orogen are obtained. We summarize the main findings as follows:
- (1)
A global optimization strategy for joint inversion of Rayleigh wave H/V ratios and phase velocities based on Particle Swarm Optimization algorithm is
Credit of authors
Sanxi Ai: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Software; Supervision; Validation; Visualization; Writing - original draft; Writing - review & editing
Yong Zheng: Funding acquisition; Investigation; Project administration; Resources; Software; Supervision; Validation; Visualization; Writing - original draft; Writing - review & editing
Sixue Wang: Formal analysis; Writing - review & editing.
Lipeng He: Formal analysis; Writing - review & editing.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgment
Teleseismic waveform data for H/V measurements in this study is provided by the TAIGER project. The waveform data are archived at Data Management Center of IRIS (www.iris.edu). The H/V ratio measurements and shear wave velocities obtained in this study are publicly deposited at https://github.com/aisanxi/Vs-1D-Taiwan. This work is jointly supported by the CEA project (2019CSES0109), the NSFC grants (41874053, 42030108), and the MOST Special Fund from the State Key Laboratory of Geological
References (68)
- et al.
Joint vp and vs tomography of Taiwan: implications for subduction-collision orogeny
Earth & Planetary Science Letters
(2014) - et al.
Penetration of mid-crustal low velocity zone across the Kunlun Fault in the NE Tibetan Plateau revealed by ambient noise tomography
Earth Planet. Sci. Lett.
(2014) - et al.
Crustal structures from the Wuyi-Yunkai orogen to the Taiwan orogen: the onshore-offshore wide-angle seismic experiments of the TAIGER and ATSEE projects
Tectonophysics
(2016) - et al.
New insights on 3-D plates interaction near Taiwan from tomography and tectonic implications
Tectonophysics
(2001) - et al.
Inversion of the Paleogene Chinese continental margin and thick-skinned deformation in the Western Foreland of Taiwan
J. Struct. Geol.
(2006) - et al.
Seismic imaging of subducting continental lower crust beneath the Pamir
Earth Planet. Sci. Lett.
(2013) - et al.
Application of particle swarm optimization to interpret Rayleigh wave dispersion curves
J. Appl. Geophys.
(2012) The particle swarm optimization algorithm: convergence analysis and parameter selection
Inf. Process. Lett.
(2003)- et al.
Crust–mantle boundaries in the Taiwan–Luzon arc-continent collision system determined from local earthquake tomography and 1D models: implications for the mode of subduction polarity reversal
Tectonophysics
(2012) - et al.
Post-collisional magmatism around northern Taiwan and its relation with opening of the Okinawa Trough
Tectonophysics
(1999)
Taiwan orogeny: thin-skinned or lithospheric collision?
Tectonophysics
Subsurface imaging, TAIGER experiments and tectonic models of Taiwan
J. Asian Earth Sci.
Seismic evidence on different rifting mechanisms in southern and northern segments of the Fenhe-Weihe Rift zone
Journal of Geophysical Research: Solid Earth
Ambient noise tomography across the Taiwan Strait, Taiwan Island, and southwestern Ryukyu Arc: Implications for subsurface slab interactions
Tectonics
Backazimuth estimation reliability using surface wave polarization
Geophys. Res. Lett.
Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements
Geophysical Journal International
Closing the gap between regional and global travel time tomography
Journal of Geophysical Research: Solid Earth
Rayleigh wave particle motion and crustal structure
Bull. Seismol. Soc. Am.
Surface deformation and tectonic setting of Taiwan inferred from a GPS velocity field
Journal of Geophysical Research: Solid Earth
Analysis of Rayleigh-wave multipath propagation at LASA
Bull. Seismol. Soc. Am.
An automated Rayleigh-wave detection algorithm
Bull. Seismol. Soc. Am.
Modern vertical deformation rates and mountain building in Taiwan from precise leveling and continuous GPS observations, 2000–2008
Journal of Geophysical Research: Solid Earth
Present-day kinematics of active mountain building in Taiwan from GPS observations during 1995–2005
Journal of Geophysical Research: Solid Earth
Critical taper model of fold-and-thrust belts and accretionary wedges
Annu. Rev. Earth Planet. Sci.
Direction of approach of Rayleigh waves and related problems (Part II)
Bull. Seismol. Soc. Am.
P-wave tomography and azimuthal anisotropy of the Manila-Taiwan-southern Ryukyu region
Tectonics
Array-analysis of two-dimensional variations in surface wave phase velocity and azimuthal anisotropy in the presence of multi-pathing interference
Phase velocities of Rayleigh waves in the MELT experiment on the East Pacific rise
Science
Geological ages of Taiwan stratigraphy and tectonic events (in Chinese)
Sci. Sin. Terrae
Investigating the lithospheric velocity structures beneath the Taiwan region by nonlinear joint inversion of local and teleseismic P wave data: Slab continuity and deflection
Geophys. Res. Lett.
Layered deformation in the Taiwan orogen
Science
Rayleigh-wave multipathing along the west coast of North America
Bull. Seismol. Soc. Am.
Particle swarm optimization
Three-dimensional P velocity structures of the lithosphere beneath Taiwan from the analysis of TAIGER and related seismic data sets
Journal of Geophysical Research: Solid Earth
Cited by (3)
Inversion of Rayleigh Wave Dispersion Curves Via BP Neural Network and PSO
2023, Research SquareLayered Evolution of the Oceanic Lithosphere Beneath the Japan Basin, the Sea of Japan
2023, Journal of Geophysical Research: Solid Earth