PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 10591.
Published online 2017 September 6. doi:  10.1038/s41598-017-09640-0
PMCID: PMC5587744

Growth of the northeastern margin of the Tibetan Plateau by squeezing up of the crust at the boundaries

Abstract

In classic orogenic models, the mountain range is underlain by a deep crustal root. Here we present the crustal and upper mantle structures along two receiver function profiles across Qilian, an orogen experiencing recent growth at the northern margin of the Tibetan plateau. Opposite to an expected crustal root beneath the orogen, the Moho beneath Qilian is arch-like, shallower beneath the center and deepens by up to 10 km beneath its southern and northern boundaries. Additional velocity interfaces sub-parallel to the Moho are observed in the lower crust of the basins south of Qilian, which we interpret as the top of a mechanically strong lower crust thrusting several tens of kilometers underneath Qilian. In the north, the small lateral offset between the surface and mantle traces of the thrust system reveals a steep boundary, indicating that the North China cratonic crust acts as a strong resistance to the northward growth of the plateau, forcing the development of the left-lateral strike-slip Haiyuan fault south of the northern Qilian suture. The young Qilian orogen thus has been rising and growing progressively from the boundaries to the center, squeezed up by more rigid tectonic blocks in the north and south.

Introduction

Classic orogenic models are characterized by a deep crustal root1, which reflects a cumulative history of complex crustal shortening and folding, and subsequent surface erosion and crustal modification2. While the general mechanism of orogeny (crustal shortening) is widely agreed upon, details of the processes, through which orogenic belts develop their crustal roots, remain unclear.

At the northeast margin of the Tibetan plateau, the Qilian terrane consists of complexly deformed earlier Paleozoic arcs3, 4 (Fig. 1). In the Carboniferous (359 to 299 million years ago), the earlier orogenic movement ceased and marine and marine-continental transition sedimentary layers covered the north Qilian4. The Cenozoic reactivation led to the rise of the present-day Qilian mountains, which are undergoing shortening and uplift58, driven presumably mainly by the collision of the Indian plate in the south and confined by the main Asian continental plate in the north. There are debates about whether the building of the Qilian mountains is driven also by southward subduction of the Asian mantle lithosphere912 or not13, 14. While controversy remains as whether deformation in the northeast Tibetan plateau started shortly after the collision of the Indian and Asian plates6, 9, there are mounting evidence showing that most of the orogenic growth in the region occurred since ~15 Ma15, 16. The present pattern of strike-slip dominated crustal deformation in the northeastern Tibetan plateau was established diachronously in late Miocene to Pliocene (11.6–2.6 Ma)5, 17, 18, making the present-day Qilian a relatively young orogen and a good place to understand early evolution of orogens.

Figure 1
Location of the broadband seismic stations on the topographic map of the northeastern Tibetan plateau. Blue and red triangles are our seismic stations deployed in 2013 and in 201012 respectively. Yellow and green triangles are our NETS seismic stations ...

Data

In order to understand the ongoing crustal deformation of the northeast Tibetan plateau and its interaction with the main Asian continental plate, we deployed 47 broadband seismographs along two roughly linear profiles across the northeast Tibetan plateau between 2013 and 2014 (Fig. 1). We used the seismic data from these 47 stations and 33 previous stations in the vicinity of our profiles in this study. The stations were equipped with Guralp 3ESP or 3T broadband seismometers. Both profiles are about 500 km long, running in the SSW–NNE directions. The western profile starts from just north of the Kunlun fault (KF), across the eastern Qaidam basin and the Qilian orogen into the North China platform (NCP, also known as the Gobi Alashan block), with its southern part following our pilot profile of the joint Sino-US North-Eastern Tibetan plateau Seismic (NETS) stations deployed between summers of 2006 and 200719 and the northern part following the Chinese Geological Survey (CGS) stations deployed between the summers of 2010 and 201112. The eastern profile extends from the northern margin of the Songpan-Ganzi terrane, across the KF, the Gonghe basin and the Qilian orogen to the NCP. The eastern profile follows the western flank of the NETS areal stations deployed from 2008 to 2009. Both profiles trend more or less perpendicular to the mountain ranges and the main tectonic boundaries in the region, including the left-lateral Kunlun strike-slip fault, the South Qilian suture (SQS), the North Qilian suture (NQS), and Yumushan fault (YF). These two linear profiles are designed for cross-section imaging of the crustal and upper mantle structure with receiver function profiling techniques (see Method).

Results

Our two receiver function profiles (Figs 2 and and3)3) reveal a crustal structure beneath Qilian that is opposite to the expected deep crustal root in classic orogenic models1, 20, 21. At about 50~58 km below the sea level (BSL), the Moho beneath Qilian is arch-like, more obvious on the western profile, deepening from ~50 km BSL at the center to 55–60 km BSL at the northern and southern boundaries of the orogen, respectively (Fig. 2, Figs S3 and S5). To highlight this arch-like feature, we split the southern, middle and northern sections of the Qilian orogen and stacked the receiver functions within each section to enhance the signal of the Moho. Figure 4 shows that the Moho beneath the center of the orogen is unequivocally shallower than its northern and southern sections along the western profile. On the eastern profile, the Moho arch is less pronounced, deepening by about 3 km from the center of the Qilian orogen to the boundaries (Figs 3 and and44 and Fig. S4). Along a similar receiver function profile ~200 km east of our eastern profile (AA’ in Fig. 1), the Moho under the less-developed east Qilian orogen is nearly flat at ~48 km BSL22.

Figure 2
Seismic image and interpretations of crustal structure along the western profile. (A) Topography and station positions along the profile. (B) Seismic cross-section constructed with 4005 P-wave receiver functions along the profile. Positive and negative ...
Figure 3
Seismic image and interpretations of crustal structure along the eastern profile. (A) Topography and station positions along the profile. (B) Seismic cross-section constructed with 3693 P-wave receiver functions along the profile. Horizontal distances ...
Figure 4
Stacked receiver functions along sections of the western and eastern profiles. (A) From left to right, stacked receive functions for the southern (x = 40–60 km in Fig. 2), middle (x = 70–110) ...

There are major differences in the way the Qilian crust interacts with the tectonic blocks north and south of the orogen. In the north, the Moho beneath the NCP is characterized by consistently high amplitudes of the converted seismic waves on the two profiles (Figs 2 and and3).3). The small perturbations to the Moho near ~39°N on the western profile and near 38.5°N on the eastern profile can be attributed to multiple waves converted and reflected at the base of sedimentary rocks (marked ‘D’ beneath the NCP in Figs 2 and and3).3). The slight northward deepening of the Moho at the northern ends of the two profiles is most likely due to the under-corrected effects of the increasingly thick sedimentary layer toward north, because the shallow velocity near the northern ends of the profiles is not well constrained due to a lack of seismic stations there23. The Moho beneath the NCP remains relatively flat, at about ~48 km BSL on the western profile and ~45 km BSL on the eastern profile. It deepens south of the YF and reaches ~58 km BSL beneath the NQS. The left-lateral strike-slip Haiyuan fault, which accommodates the eastward motion of the crust in the northeastern Tibetan plateau, is tens of kilometers south of the NQS near the western profile and merges with the NQS near the eastern profile (Fig. 1). Beneath both profiles, the most depressed points of the Moho near the southern end of the NCP section are directly beneath the surface trace of the NQS (Figs 2, ,33 and Fig. S5). The lateral distance between the surface trace of the thrust fault and the most depressed points of the Moho are about 25 km on the two profiles, suggesting a steep thrust system between the Qilian orogen and the NCP, in contrast to the shallow underthrust at the southern margin of the Tibetan plateau24. This observation supports the notion of southward under-thrusting of the NCP crust11, at least at the crustal level, as under-thrusting of the crust and mantle lithosphere places the southern edge of the NCP Moho south of the surface trace of the thrust fault, as observed.

South of the Qilian orogen, the Qaidam Moho is seen clearly deepening to north from ~48 km BSL north of the Kunlun mountains to ~63 km BSL beneath the SQS on the western profile (Fig. 2 and Fig. S3). From this point to farther north, the Moho becomes less certain to discern on this profile, but for reasons given below, we infer that the Qaidam Moho extends north of the surface trace of the SQS. In contrast, the corresponding Moho on the eastern profile is nearly horizontal south of 36°N, at ~50 km BSL (Fig. 3). Farther north, it deepens slightly, to ~57 km BSL, and can be traced continuously to ~30 km north of the SQS (Fig. 3). The Moho is somewhat blurred beneath the Gonghe basin (~35.7–36.3°N) on the eastern profile, probably due to multiple waves converted and reflected at the bottom of the Gonghe basin since the directly converted-wave (marked ‘S’ beneath the Gonghe basin in Fig. 3) is very strong and clearly seen on the profile.

At ~15–20 km above the Moho south and near the SQS, there is a prominent phase corresponding to a positive velocity contrast interface (increasing velocity with depth) at ~40–60 km BSL on the western profile and at ~30–40 km BSL on the eastern profile (marked ‘1’ in Figs 2 and and3).3). The southward extension of this velocity interface is obscured by the multiple waves from the basin sedimentary layer. Feng et al. (2014) also observed such an interface12 nearly along our western profile and stressed its similarity with the doublet phase observed beneath southern Tibet11, 25, 26. The doublet phase beneath southern Tibet appears more or less horizontal within a narrow depth range of ~55–58 km BSL, consistent with the interpretation of an eclogitization phase change boundary11, 26, 27. In contrast, the lower crustal velocity interface beneath the Qaidam and Gonghe basins runs parallel to the Moho and is clearly not horizontal (thus at very different confining pressures). If this interface is the eclogite phase boundary, the deepening interface requires a substantially colder lower crust at a greater depth towards the SQS, an unlikely circumstance. We thus interpret this lower crustal velocity interface as the top of the lower crust, which is likely mafic to have the high velocity needed to create a large enough velocity contrast with the crust above. In either case, the Moho, the bottom of this lower crust or an unlikely eclogitized lower crust, must follow the top interface to 20–40 km north of the SQS beneath the western profile (Fig. 2), similar to the Moho on the eastern profile (Fig. 3). A possible explanation for the lack of a clear arrival from the Moho north of the SQS on the western profile is that a strong lower crustal layer at the bottom of the Qaidam crust at greater depths may significantly reduce the impedance contrast of this part of the Moho, therefore reducing the amplitude of the converted seismic waves and making this part of the Moho more difficult to identify.

Discussion

The western Qilian terrane appears to be a mechanically weak and recently reactivated orogen bounded by two mechanically strong tectonic blocks: The NCP and Qaidam basin. The west Qilian is characterized by a low velocity zone in its mid-to-lower crust and its Vp/Vs ratio is typical for more felsic crustal composition23, 28. In the north, the Qilian crust abuts the western block of the North China craton consisting of a basement of Precambrian metamorphic rocks29, which is generally considered to be older, colder and mechanically stronger than the Qilian’s basement that experienced oceanic subduction and arc magmatism between 520–400 Ma3. GPS data show that the strain rate within the Qilian orogen is an order of magnitude greater than the rate within the western block of the north China craton7, 8. In the south, the basement of the Qaidam basin is composed of Precambrian-Silurian ( >419 Ma) metamorphic rocks30. The Qaidam crust and the mantle lithosphere appears to be exceptionally strong mechanically, with an effective elastic thickness of ~70 km, significantly greater than that of 10–30 km for the rest of the Tibetan plateau31. Seismic studies show that the lower crust beneath the Qaidam and Gonghe basins has a much higher velocity than in the surrounding tectonic units23, also indicating a mechanically strong lower crust.

Taken together, the receiver function profiles suggest that the recent thickening of the Qilian orogenic crust is most developed at its northern and southern boundaries, squeezed up by under-thrusting of a mechanically strong Qaidam crust in the south and blocked by the North China cratonic crust in the north. Analysis of surface geology and reconstruction of seismic reflection profiles reveal that the shortening along the northern Qilian Shan frontal thrust system is much higher than in the interior of Qilian32, consistent with our observation of a thicker crust near the northern boundary. The arch-like Moho across Qilian may characterize crustal deformation of young orogens sandwiched between strong tectonic blocks. The surface topography and Moho display mirror images, particularly on the eastern profile (Fig. 3), reflecting a component of the “pure-shear” mode of crustal deformation33 on the cross-sections, in a region otherwise dominated by strike-slip crustal deformation34.

The development of the Moho arch appears to be correlated with the development of the orogen from the planar Moho under the less-developed eastern Qilian orogen 200 km east of our eastern profile22, to the small Moho arch under our eastern profile (Figs 3 and and4)4) and finally the most pronounced Moho arch beneath the more developed western Qilian (Figs 2 and and4).4). We speculate that as the tectonic blocks continue to converge, crustal thickening will propagate from the boundaries to the center and the arch-like Moho will eventually disappear, forming a classic orogenic crust root beneath the orogen.

Methods

In this study, we obtained the images of the crustal and upper-mantle structure of the NE Tibetan plateau using common-conversion-point (CCP) stacks of P-wave receiver functions (RFs), a well-established method that has been applied in many similar studies of the Tibetan plateau11, 26. This method enhances converted S waves produced by P waves of distant earthquakes impinging on seismic interfaces. The interfaces and their properties can be determined from the converted waves based on the time delays between the direct and converted waves, which are mainly proportional to the depths of the interfaces, and on the amplitudes of the converted waves, which depend on the magnitudes and signs of the velocity contrasts.

We used teleseismic events with magnitudes greater than 4.5 and epicentral distances greater than 30°. Most of them occurred to the east of the profiles (Fig. S1). To ensure that only good-quality waveforms are used, we also applied a visual selection process to all the data. We performed time-domain iterative deconvolution35 in the calculation of RFs, in which the Gaussian filter factor was set to 2.5 to retain sufficient high frequency (up to 1 Hz) signals to image the crustal and upper-mantle structure. We isolated the converted waves from the incident waves using the wave-vector method36 which decomposes P- and S-waves and removes the effect of the free surface by multiplying the inverse matrix of free-surface response with the observed vertical, radial and transverse components of data. A free-surface response matrix is calculated for surface velocities of Vp = 5.34 km/s and Vs = 3.00 km/s, which removed most of the effect of the free surface. We used variable reference models for the individual stations (Fig. S2) to migrate the RFs of different seismic stations from time to space in a 3-D volume by tracing the rays from the location of each seismic station. The reference models were obtained by interpolations at each station of the shear-wave velocity model derived from joint inversion of receiver functions and Rayleigh wave dispersion23, and the Vp/Vs ratios derived from h-k stacking of RFs obtained previously23 and in this study (Table S1). Then all the RFs were projected and stacked onto cross-sections aligned to due north with a fixed horizontal stacking bin width of 1 km.

Data Availability

The original seismic data from the NETS and ASCENT experiments are available at the IRIS DMC. Receiver functions from all the stations used in this study, including non-IRIS data, are available upon request from D.S.

Electronic supplementary material

Acknowledgements

This research was funded by the Chinese Geological Survey (1212011220903), the National Natural Science Foundation (41374109) of China, the Ministry of Science and Technology of China (2006DFA21350) and U.S. National Science foundation (0738779). Dr. Eric Sandvol provided assistants in the field. We thank Drs. Xiaoming Xu, Xinfu Li, Xiaofeng Liang, Meng Gong, Jinsheng Jia and Hongyi Li for the participation of fieldwork. Comments by Dr. An Yin and two anonymous reviewers greatly helped improving the manuscript.

Author Contributions

Author Contributions

J.S. conducted the data processing, D.S. led the data collection and analysis. J.S. and D.S. was involved in earlier drafts of the manuscript. Y.S. was involved in the data analysis and led the writing of the submitted manuscript. W.Z., G.X., H.S. and Y. Song contributed to the data collection, understanding of the results and editing of the manuscript.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-09640-0

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Dewey JF, Bird JM. Mountain belts and the new global tectonics. J. Geophys. Res. 1970;75:2625–2647. doi: 10.1029/JB075i014p02625. [Cross Ref]
2. Fischer KM. Waning buoyancy in the crustal roots of old mountains. Nature. 2002;417:933–936. doi: 10.1038/nature00855. [PubMed] [Cross Ref]
3. Yin A, et al. Early Paleozoic Tectonic and Thermomechanical Evolution of Ultrahigh-Pressure (UHP) Metamorphic Rocks in the Northern Tibetan Plateau, Northwest China. International Geology Review. 2007;49:681–716. doi: 10.2747/0020-6814.49.8.681. [Cross Ref]
4. Song S, Niu Y, Su L, Xia X. Tectonics of the north Qilian orogen, NW China. Gondwana Res. 2013;23:1378–1401. doi: 10.1016/j.gr.2012.02.004. [Cross Ref]
5. Chen Z, et al. Global positioning system measurements from eastern Tibet and their implications for India/Eurasia intercontinental deformation. J. Geophys. Res. 2000;105(16):215–16,227.
6. Yin A, Harrison TM. Geologic evolution of the Himalayan–Tibetan orogen. Annu. Rev. Earth Planet. Sci. 2000;28:211–280. doi: 10.1146/annurev.earth.28.1.211. [Cross Ref]
7. Wang Q, et al. Present-day crustal deformation in China constrained by global positioning system measurements. Science. 2001;294:574–577. doi: 10.1126/science.1063647. [PubMed] [Cross Ref]
8. Zhang PZ, et al. Continuous deformation of the Tibetan plateau from global positioning system data. Geology. 2004;32:809–812. doi: 10.1130/G20554.1. [Cross Ref]
9. Meyer B, et al. Crustal thickening in Gausu-Qinghai, lithospheric mantle subduction, and oblique, strike-slip controlled growth of the Tibet plateau. Geophys. J. Int. 1998;135:1–47. doi: 10.1046/j.1365-246X.1998.00567.x. [Cross Ref]
10. Tapponnier P, et al. Oblique stepwise rise and growth of the Tibet Plateau. Science. 2001;294:1671–1677. doi: 10.1126/science.105978. [PubMed] [Cross Ref]
11. Kind R, et al. Seismic images of crust and upper mantle beneath Tibet: Evidence for Eurasian plate subduction. Science. 2002;298:1219–1221. doi: 10.1126/science.1078115. [PubMed] [Cross Ref]
12. Feng M, et al. Structure of the crust and mantle down to 700km depth beneath the East Qaidam basin and Qilian Shan from P and S receiver functions. Geophys. J. Int. 2014;199:1416–1429. doi: 10.1093/gji/ggu335. [Cross Ref]
13. Yue H, et al. Lithospheric and upper mantle structure of the northeastern Tibetan Plateau. J. Geophys. Res. 2012;117
14. Liang X, et al. A complex Tibetan upper mantle: a fragmented Indian slab and no south-verging subduction of Eurasian lithosphere. Earth Planet. Sci. Lett. 2012;333–334:101–111. doi: 10.1016/j.epsl.2012.03.036. [Cross Ref]
15. Fang X, Li J, Song C, Yan M. Late Miocene-Quaternary synchronous-but-magnitude-different episodic rapid uplifts of the NE Tibetan plateau: A synthesis from flexural basins. ACTA Geol. Sinica (English ed.) 2013;87(supp.):207–209.
16. Yuan D-Y, et al. The growth of northeastern Tibet and its relevance to large-scale continental geodynamics: A review of recent studies. Tectonics. 2013;32:1358–1370. doi: 10.1002/tect.20081. [Cross Ref]
17. Tapponnier P, et al. Active thrusting and folding in the Qilian Shan, and decoupling between upper crust and mantle in northeastern Tibet. Earth Planet. Sci. Lett. 1990;97:382–403. doi: 10.1016/0012-821X(90)90053-Z. [Cross Ref]
18. Zheng D, Clark MK, Zhang P, Zheng W, Farley KA. Erosion, fault initiation and topographic growth of the North Qilian Shan (northern Tibetan Plateau) Geosphere. 2010;6:937–941. doi: 10.1130/GES00523.1. [Cross Ref]
19. Shi D, Shen Y, Zhao W, Li A. Seismic evidence for a Moho offset and south-directed thrust at the easternmost Qaidam–Kunlun boundary in the Northeast Tibetan plateau. Earth planet. Sci. Lett. 2009;288:329–334. doi: 10.1016/j.epsl.2009.09.036. [Cross Ref]
20. Oxburgh ER. The deep structure of orogenic belts - the root problem. Geological Society, London, Special publications. 1969;3:251–273. doi: 10.1144/GSL.SP.1969.003.01.13. [Cross Ref]
21. Vanderhaeghe O, Medvedev S, Fullsack P, Beaumont C, Jamieson RA. Evolution of orogenic wedges and continental plateau: insights from crustal thermal-mechanical models overlying subducting mantle lithosphere. Geophys. J. Int. 2003;153:27–51. doi: 10.1046/j.1365-246X.2003.01861.x. [Cross Ref]
22. Ye Z, et al. Seismic evidence for the North China plate underthrusting beneath northeastern Tibet and its implications for plateau growth. Earth Planet. Sci. Lett. 2015;426:109–117. doi: 10.1016/j.epsl.2015.06.024. [Cross Ref]
23. Zheng D, et al. Crustal and upper mantle structure beneath the northeastern Tibetan Plateau by joint analysis of receiver functions and Rayleigh wave dispersion. Geophys. J. Int. 2016;204:583–590. doi: 10.1093/gji/ggv469. [Cross Ref]
24. Schulte-Pelkum V, et al. Imaging the Indian subcontinent beneath the Himalaya. Nature. 2005;435:1222–1225. doi: 10.1038/nature03678. [PubMed] [Cross Ref]
25. Yuan X, Ni J, Kind R, Mechie J, Sandvol E. Lithospheric and upper mantle structure of southern Tibet from a seismological passive source experiment. J. Geophys. Res. 1997;102(27):491–27,500.
26. Nábělek J, et al. Underplating in the Himalaya-Tibet collision zone revealed by the Hi-CLIMB experiment. Science. 2009;325:1371–1374. doi: 10.1126/science.1167719. [PubMed] [Cross Ref]
27. Shi D, et al. Receiver Function Imaging of Crustal Suture, Steep Subduction, and Mantle Wedge in the Eastern India-Tibet Continental Collision Zone. Earth Planet. Sci. Lett. 2015;414:6–15. doi: 10.1016/j.epsl.2014.12.055. [Cross Ref]
28. Li H, et al. The distribution of the mid-to-lower crustal low-velocity zone beneath the northeastern Tibetan Plateau revealed from ambient noise tomography. J. Geophys. Res. Solid Earth. 2014;119:1954–1970. doi: 10.1002/2013JB010374. [Cross Ref]
29. Vincent SJ, Allen MB. Evolution of the Minle and Chaoshui Basins, China: Implications for Mesozoic strike-slip basin formation in central Asia. Geol. Soc. Am. Bull. 1999;111:725–742. doi: 10.1130/0016-7606(1999)111<0725:EOTMAC>2.3.CO;2. [Cross Ref]
30. Wu, C., Yin, A., Zhang, J.Y., Liu, W.C. & Ding, L. Geologic history of the Paleozoic-Mesozoic Eastern Kunlun Arc, central Tibetan Plateau: Implications for the tectonic evolution of the Tethyan orogenic system. Lithosphere, in press.
31. Braitenberg C, Wang Y, Fang J, Hsu HT. Spatial variations of flexure parameters over the Tibet-Qinghai plateau. Earth Planet. Sci. Lett. 2003;205:211–224. doi: 10.1016/S0012-821X(02)01042-7. [Cross Ref]
32. Zuza, A.V., Cheng, X. & Yin, A. Testing models of Tibetan plateau formation with Cenozoic shortening estimates across the Qilian Shan-Nan Shan thrust belt. Geosphere 12 (2), doi:10.1130/GES01254.1 (2016).
33. England PC, Houseman GA. Finite strain calculations of continental deformation, 2, Comparison with the India-Asia collision zone. J. Geophys. Res. 1986;91:3664–3676. doi: 10.1029/JB091iB03p03664. [Cross Ref]
34. Tapponnier, P., Peltzer, G., Le Dain, A.Y., Armijo, R. & Cobbold, P. Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine. Geology10, 611–616, doi:10.1130/0091-7613(1982) 10<611:PETIAN>2.0.CO;2 (1982).
35. Ligorria, J. & Ammon, C.J. Iterative deconvolution and receiver function estimation. Bull. Seismol. Soc. Am. 89, 1395–1400, doi:10.1.1.473.9379 (1999).
36. Reading A, Kennett B, Sambridge M. Improved inversion for seismic structure using transformed, S-wave vector receiver functions: Removing the effect of the free surface. Geophys. Res. Lett. 2003;30 doi: 10.1029/2003GL018090. [Cross Ref]
37. Wessel P, Smith WHF. New, improved version of generic mapping tools released. EOS Trans. 1998;79 doi: 10.1029/98EO00426. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group