communications earth & environment
Extremely rapid up-and-down motions of island arc crust during arc-continent collision
Your article has downloaded
Similar articles being viewed by others
Carousel with three slides shown at a time. Use the Previous and Next buttons to navigate three slides at a time, or the slide dot buttons at the end to jump three slides at a time.
Magmatic flare-up causes crustal thickening at the transition from subduction to continental collision
Carlos E. Ganade, Pierre Lanari, … Caroline M. Ribeiro
Long-lived Paleoproterozoic eclogitic lower crust
Sebastian Buntin, Irina M. Artemieva, … Stefan Buske
Post-collisional mantle delamination in the Dinarides implied from staircases of Oligo-Miocene uplifted marine terraces
Philipp Balling, Christoph Grützner, … Kamil Ustaszewski
Episodic construction of the early Andean Cordillera unravelled by zircon petrochronology
José Joaquín Jara, Fernando Barra, … Diego Morata
Arc accretion and crustal reworking from late Archean to Neoproterozoic in Northeast Brazil
Alanielson C. D. Ferreira, Elton L. Dantas, … Ingrid M. Nedel
No mafic layer in 80 km thick Tibetan crust
Gaochun Wang, Hans Thybo & Irina M. Artemieva
Archean crust and metallogenic zones in the Amazonian Craton sensed by satellite gravity data
J. G. Motta, C. R. de SOUZA FILHO, … C. Braitenberg
A 60-Myr record of continental back-arc differentiation through cyclic melting
Lauren C. Wolfram, Roberto F. Weinberg, … Raul Becchio
A likely geological record of deep tremor and slow slip events from a subducted continental broken formation
Francesco Giuntoli & Giulio Viola
Larry Syu-Heng Lai ORCID: orcid.org/0000-0003-4589-97291,
Rebecca J. Dorsey1,
Chorng-Shern Horng2,
Wen-Rong Chi3,4,
Kai-Shuan Shea5 &
Jiun-Yee Yen6
Communications Earth & Environment volume 3, Article number: 100 (2022) Cite this article
Palaeomagnetism
Tectonics
Mountain building and the rock cycle often involve large vertical crustal motions, but their rates and timescales in unmetamorphosed rocks remain poorly understood. We utilize high-resolution magneto-biostratigraphy and backstripping analysis of marine deposits in an active arc-continent suture zone of eastern Taiwan to document short cycles of vertical crustal oscillations. A basal unconformity formed on Miocene volcanic arc crust in an uplifting forebulge starting ~6 Ma, followed by rapid foredeep subsidence at 2.3–3.2 mm yr−1 (~3.4–0.5 Ma) in response to oceanward-migrating flexural wave. Since ~0.8–0.5 Ma, arc crust has undergone extremely rapid (~9.0–14.4 mm yr−1) uplift to form the modern Coastal Range during transpressional strain. The northern sector may have recently entered another phase of subsidence related to a subduction polarity reversal. These transient vertical crustal motions are under-detected by thermochronologic methods, but are likely characteristic of continental growth by arc accretion over geologic timescales.
Vertical crustal motions are fundamental to the creation of topography, development of sedimentary basins, and the rock cycle1,2. Rapid vertical-displacement rates (mm year−1) are often driven by tectonic processes such as crustal thickening in Tibet3 and Taiwan4,5,6, strike-slip deformation along the San Andreas fault7,8, lithospheric thinning in Central Anatolia9 and D'Entrecasteaux island of Papua New Guinea10, and deflections due to changes in surface or subsurface mass loads in Hawaii11 and Antarctica12. Mass redistribution by erosion and sedimentation amplifies vertical crustal motions13 and is important for understanding interactions between tectonic and surface processes involved in mountain building and continental growth14. Among these settings, extreme rates (>10 mm year−1) of long-term (>105–107 years) rock uplift are seldom detected15. Many studies rely on thermochronology16 and petrology-geochemistry17 to assess long-term vertical movements of crustal materials, but these methods are limited by the requirement of suitable minerals, geothermal gradient, temperature, and duration needed to thermally reset thermochronometers. Assumptions of ancient topography and thermal structure are often essential for use of thermochronologic methods to estimate long-term rock uplift18, but it is difficult to constrain these quantities in rapidly growing orogens constructed of unmetamorphosed rocks. Thus, the rates, timescales, and structural controls on vertical motions of shallow crust in active mountain belts remain poorly understood.
Stratigraphic study of inverted syn-orogenic sedimentary basins provides a powerful tool with which to document rapid subsidence and uplift of unmetamorphosed near-surface rocks in zones of mountain building at tectonically active plate margins2. Integrated paleomagnetic and biostratigraphic analysis can yield extremely high resolution [e.g., ± 5–20 thousand years (kyr)] for dating stratigraphic intervals7, and thus, remarkably high-fidelity estimates for rates of vertical crustal motions and tectonic processes. We applied this technique to investigate stratigraphic records of the Coastal Range in eastern Taiwan, which reveal a history of extremely rapid vertical crustal oscillations during accretion of volcanic-arc crust to a continent.
The island of Taiwan has emerged since late Miocene time through active collision between the Luzon island arc on the Philippine Sea plate and the Chinese continental margin of the Eurasian plate19 (Fig. 1). Rapid (~82 mm year−1) oblique convergence20 between the two plates induces rapid exhumation and denudation21, and some studies infer that collision has propagated southward through time22,23,24. Other studies of the metamorphic core (Central Range) and western foreland basin suggest that collision was geologically simultaneous from north to south, with pulses of accelerated exhumation from ~0.1 to 2–4 mm year−1 at 2.0–1.5 million years ago (Ma) and then to 4–8 mm year−1 at ~0.5 Ma25,26. These accelerations correspond to tectonic reorganizations of the overriding Philippine Sea plate27,28 that drive exhumation of high-pressure metamorphosed Miocene arc crust-bearing mélange (Yuli Belt) from >35 km crustal depth5,29,30,31,32 and rapid Quaternary emergence of the unmetamorphosed arc crust in the Coastal Range4.
Fig. 1: The Coastal Range of eastern Taiwan.
a Regional tectonic configuration5,28 (inset), simplified geological map of the Coastal Range34 and analyzed stratigraphic sections (blue boxes), including Hsiukuluan river (HKL), Wulou river (WL), Fungfu (FF), and Fungpin (FP) sections in the north, and Bieh river (BC), Madagida river (MDJ), and Sanshian river (SSS) sections in the south. b Millennial rates of marine terrace uplift and river incision, compiled by Lai et al.99 c Geodetic rates of vertical deformation measured during 2000–200849. Negative values mean subsidence.
The Coastal Range contains a thick succession of Plio-Pleistocene orogen-derived marine flysch, conglomerate (Fanshuliao and Paliwan formations), and olistostromes (Lichi Mélange) that rest unconformably on Miocene volcanic arc basement (Tuluanshan Formation, ca. 15–6 Ma) (Fig. 2)33,34. The basal unconformity (ca. 6–4 Ma) is a broad erosive surface discontinuously capped by thin shallow-marine limestone (Kangkou Limestone, 5.6–3.5 Ma) and limestone-clast-bearing epiclastic deposits (Biehchi Epiclastic Unit, 5–3.5 Ma) that are directly overlain by uncemented deep-water flysch34,35. These relations record slow uplift, erosion, and intermittent sedimentation on arc basement near sea level, followed by rapid subsidence that created accommodation space for thick basin-filling sediment6,33. Basin inversion and uplift must have occurred in the last roughly 1 million years based on observation of the tilted youngest marine flysch and creation of modern Coastal Range topography4. Yet, the timing and rates of these processes remain unclear. Thermochronology and petrological methods are not suitable for solving this problem because of the short duration of burial and very low paleo-geothermal gradient revealed by clay mineralogy (~14 °C km−1)36 and no post-depositional resetting of fission-track detrital thermochronometers37. Early studies used calcareous nannoplankton data38,39 and simplified lithostratigraphic columns to derive a range of subsidence rates (0.8–5 mm year−1) and minimum uplift rate (5.9–7.5 mm year−1) with unknown spatial variability4,6. Many of these studies are now outdated and pre-date recent advances in paleomagnetism and microfossil studies34,40,41,42,43.
Fig. 2: Stratigraphy of the Coastal Range.
Locations of studied sections are shown in Fig. 1. GPTS Geomagnetic polarity timescale. FAD First appearance datums of index microfossils. Tls Tuluanshan Formation; k Kangkou Limestone. See details of stratigraphic correlation in Supplementary Figs. 1, 2.
To document the timing, magnitude, and rates of vertical crustal motions in the Coastal Range of eastern Taiwan, we compiled geologic and magneto-biostratigraphic data to date and accurately reconstruct two composite stratigraphic columns in the northern and southern Coastal Range (Figs. 1, 2). This includes our published data from the southern Coastal Range34 (Supplementary Fig. 1) and new detailed geologic mapping, paleomagnetic measurements, and microfossil identifications of planktonic foraminifera and calcareous nannoplankton in the northern Coastal Range (Supplementary Figs. 2–6; Supplementary Data 1–2). The foraminifera data provide improved paleobathymetry estimates (Supplementary Fig. 7, Supplementary Data 3), and compiled magneto-biostratigraphy constrains depositional age (Fig. 3a). Dense sampling of microfossils and paleomagnetic sites yields high temporal resolution (~1–15 kyr) and high accuracy (2–43 kyr) of age controls that surpass other geologic dating methods (Supplementary Data 4). Because the youngest inverted sediment is consolidated and lithified, we have to account for deposits that accumulated above the top of our measured sections and subsequently were removed by erosion. Porosity-effective stress of sandstone and vitrinite reflectance data indicate that ca. 0.45–1.95 km of strata was eroded off the youngest deposits preserved in our measured sections44. Due to uncertain variability of compaction history and geothermal structure, we used a conservative thickness of 0.5–1.0 km for eroded sediments to reconstruct the deepest subsidence prior to onset of structural inversion (Fig. 3a).
Fig. 3: Sediment-accumulation rates and subsidence-uplift histories of the Coastal Range.
a Stratigraphic age models for the Coastal Range. b, c Geohistory curves of the northern and southern Coastal Range. Symbols and color-fill style follow Fig. 3a. Asterisk in each curve is the reconstructed depth using youngest sediment eroded from top of section. Error bars represent one standard error of the mean. See details of data and calculation results in Supplementary Data 4.
Using high-fidelity constraints on stratigraphy and paleobathymetry, updated eustatic sea level curve (Supplementary Fig. 7), and porosity-depth functions for relevant sediment types (Supplementary Fig. 8), we conducted a modern 1-D backstripping analysis45 to progressively remove decompacted sediment and correct paleo-water depth along two composite sections. This allowed us to reconstruct the history of subsidence and uplift of arc basement in the north and south (Fig. 3b, c). See details in the "Methods" section.
Subsidence-uplift histories of Taiwan's Coastal Range
Our results reveal that >5.48–6.51 km of preserved orogen-derived sediment accumulated between ~3.39 and 0.77 Ma, with a minor increase in sedimentation rate at ~2.0 Ma (Fig. 3a). Volcanic-arc basement subsided to depths of 6.53–7.78 km below modern sea level at rates of 2.26 mm year−1 in the north and 3.24 mm year−1 in the south, with tectonic forces and sediment loads making subequal contributions to the total subsidence (Fig. 3b, c). The age of youngest, now-eroded sediment is estimated to be ~0.61–0.50 Ma, providing a reasonable age estimate for the end of subsidence and onset of structural inversion (Fig. 3a). Since the folded basal unconformity and underlying volcanic basement rocks are exposed well above sea level, the difference between present-day heights of antiformal peaks and the reconstructed depth of arc basement at the end of deposition (0.82–0.77 Ma for the youngest preserved sediment; 0.61–0.50 Ma at the top of eroded sediment) provide a conservative estimate for the amount of vertical displacement during basin inversion. To exhume the volcanic-arc basement to present elevations of the Coastal Range (0–1.33 km in the north, 0–1.68 km in the south) requires extremely rapid rock uplift rates of at least 8.89–14.39 mm year−1 in the north and 10.43–14.24 mm year−1 in the south (Fig. 3b, c; Supplementary Data 4).
The minor increase in sedimentation rate at ~2.0 Ma and onset of basin inversion at ~0.8–0.5 Ma are coeval with proposed episodes of accelerated rock exhumation in the metamorphic Central Range to the west26 (Figs. 3a, 4a). Post-0.5 Ma uplift rates of ca. 9–14 mm year−1 represent the lower bound of long-term average exhumation rates in the Coastal Range. These rates are intermediate between spatially variable millennial uplift rates measured in coastal marine terraces (2.3–11.8 mm year−1)46,47 and fluvial incision rates (15.1–27.3 mm year−1) measured near the western major oblique thrust fault (Longitudinal Valley fault)48 (Fig. 1b). Notably, millennial uplift rates in the northern Coastal Range (2.3–4.7 mm year−1)46,47 are considerably slower than our calculated long-term exhumation rates (~9–14 mm year−1), and geodetic data49 reveal subsidence at rates locally up to ~23.5 mm year−1 in the north (Fig. 1b, c). This implies that the northern part of inverted arc crust may have entered a new subsidence stage very recently.
Fig. 4: Conceptual model for the history of vertical motions (blue arrows) of the Philippine Sea plate arc crust in the past ~6 Ma.
a Forebulge uplift and foredeep subsidence occur on the Luzon arc crust during 6–0.5 Ma, in response to eastward migrating orogenic load33,34. b Post-0.5 Ma rapid uplift due to transpressional deformation and isostatic rebound of the overriding plate4,5,59. c Representative structural cross-section in the modern Coastal Range34,57. Blue lines mark major boundary faults. LVF Longitudinal Valley fault; CRF Central Range fault.
Taken together, the stratigraphic record in the Coastal Range reveals two short cycles of up-and-down crustal motions (Fig. 3b, c). Coeval histories of uplift and subsidence in the north and south support an interpretation of simultaneous crustal dynamics along the collisional suture25,26,28, in contrast to southward-propagating growth of the Taiwan orogen inferred in some models22,23,24. These findings indicate the need for a revised tectonic interpretation to explain the observed rapid vertical oscillations of arc crust in this active arc-continent suture zone.
Drivers of rapid crustal oscillations during oceanic arc accretion
We infer that the first cycle of uplift (~6–3.4 Ma) and subsidence (~3.4–0.5 Ma) of volcanic basement was driven by early flexure and loading of accreting Luzon arc crust. This interpretation is supported by recent results showing that Plio-Pleistocene sediments (Lichi Mélange, Fanshuliao and Paliwan formations) formed on – and were derived from – an east-dipping submarine paleoslope at the steep western margin of the basin34,50. These sediments onlap the westward-inclined basal unconformity on top of Miocene arc basement (Tuluanshan Formation) and display lateral facies changes that record eastward progradation of coarse deposits into the basin34. These observations are best explained by basinward migration of the depocenter in response to eastward migration of a thrust-bounded submarine slope at the retrowedge orogenic front (Fig. 4a). Thus, the first stage of vertical crustal oscillation in the Coastal Range (slow uplift) is interpreted as a signal of forebulge uplift. This stage was followed by subsidence in an evolving retrowedge foredeep basin driven by eastward migration of a flexural wave6,33. These processes took place during development of the prowedge foreland basin in western Taiwan starting ~6.5–5 Ma25,51.
The absence of transgressive deposits at the basal unconformity between shallow-marine limestone (5.6–3.5 Ma) and overlying deep-marine flysch (~3.4–3.1 Ma) requires a sudden increase in water depth at ~3.5 Ma (Fig. 3b, c). This abrupt subsidence likely resulted from local extensional faulting52,53 in the upper dilating crust of a migrating flexural forebulge54, and may be in part related to an increase in orogenic load on the Philippine Sea plate due to accelerated topographic growth in the retrowedge at 3.2 ± 0.6 Ma55. Lithospheric bending of the Philippine Sea plate induced by northward subduction at the Ryukyu trench56 may have also influenced vertical displacements in the accreting arc crust. However, this hypothesis predicts southward migration (relative to the Coastal Range) of vertical-displacement signals, but the observed cycles of uplift and subsidence occurred simultaneously in northern and southern Coastal Range transects (Fig. 3b, c). Therefore, we suggest that northward subduction at the Ryukyu-trench exerted little or no influence on vertical crustal motions during ~6–0.5 Ma basin evolution.
Post-0.5 Ma extremely rapid uplift and basin inversion created the modern topography of the Coastal Range in a small doubly-vergent structural wedge during an abrupt change to wrench-style transpressional deformation34,57 (Fig. 4b, c). While rapid vertical displacements are observed in other oblique-convergent settings7,8, it is unusual to document large-scale exhumation at rates >10 mm year−1 driven solely by upper-crustal deformation. Isostatic adjustments to changes in lithospheric structure provide a mechanism that can explain this behavior58. A tectonic load on the Philippine Sea plate may have been suddenly released when forearc lithosphere was broken and subducted at the collisional suture zone59,60,61. Downward extraction of lithospheric fragments may have caused rapid exhumation of the overriding plate crust near the suture62,63, which may be partially responsible for the post-0.5 Ma uplift of metamorphosed (Yuli Belt) and unmetamorphosed (Coastal Range) arc crust in eastern Taiwan (Fig. 4b, c)5,34,64. Isostatic compensation by erosional unloading may also have contributed to the creation of topography65,66, but its role remains unclear. It is possible that we have underestimated post-0.5 Ma long-term uplift rates in the northern Coastal Range, because present-day topography in the north may currently be influenced by the modern subsidence stage. This youngest and ongoing subsidence likely is driven by northward convergence and reversed subduction polarity of the Philippine Sea plate at the Ryukyu trench67,68 (Fig. 1).
This study reveals a history of extremely rapid vertical oscillations (up-and-down cycles) of near-surface rocks in only ~3 Myr during active accretion of island arc crust to a growing continental margin (Fig. 3). Our results highlight the short-lived episodic nature of arc-continent collision69 and challenge ideas about long timescales that are often invoked for the rise and fall of mountainous topography by coupled tectonic and surface processes (106–107 years)70. We find that growth of topography to form an eroding steep mountain range directly on the footprint of a formerly subsiding marine basin can be accomplished in only ~500–800 kyr. This study shows that the change of direction in vertical motion, and transformation of a subsiding basin to an eroding mountainous source, can occur very quickly to drive extremely rapid rock cycling and mixing at this and other collisional plate boundaries34. Because arc collision and accretion are recognized as a fundamental process in the growth of continental crust71, our results suggest that extremely rapid vertical crustal motions may be characteristic of deformation in arc-continent suture zones through geologic time72.
Existing methods of thermochronology, petrology, and geochemistry provide a useful record of crustal burial and exhumation in these settings, but they cannot detect such short rapid vertical trajectories in shallow unmetamorphosed crust. Integrated magneto-biostratigraphy and basin subsidence analysis thus offers an important tool for documenting high-resolution histories of vertical crustal motions in tectonically active systems. This approach spans a crucial time gap that needs to be bridged to better understand dynamic feedbacks among long-term geologic (tectonic) processes and shorter timescale impacts such as climate change and landscape response1,70,73.
Geological mapping and lithostratigraphy
Detailed geological mapping for this study targeted excellent exposures in road cuts and riverbanks of the Coastal Range (Fig. 1). Marker beds (pebbly mudstone and tuffaceous turbidite) and fault zones were carefully mapped across the study area (Supplementary Figs. 3–4). We also compiled information from previously published geological maps74,75,76,77,78,79,80,81. Lithostratigraphic descriptions in type sections of the southern Coastal Range [Bieh river (BC), Madagida river (MDJ), and Sanshian river (SSS) sections] are from Lai et al.34 (Supplementary Fig. 1). New results in the northern Coastal Range [Hsiukuluan river (HKL), Wulou river (WL), Fungfu (FF), and Fungpin (FP) sections] are compiled in Supplementary Fig. 2. We produced composite columns for the northern and southern Coastal Range (Fig. 2) through correlations based on marker beds and the first appearance datums (FAD) of index fossils (Supplementary Figs. 1, 2) and aided by construction of a balanced geological cross-section (Supplementary Fig. 5). This approach assumes that the thickness (i.e., rock volume) of each unit does not change substantially across the local structures (faults and folds).
Magnetostratigraphy
In the southern Coastal Range, we adapted results and interpretations of paleomagnetic chrons from previous published studies34,42,79,80 (Supplementary Fig. 1). In the northern Coastal Range, we compiled published paleomagnetism data40,41 along with new data from samples collected in coherent strata from continuous sections (Paliwan Formation) for Hsiukuluan river (HKL), Wulou river (WL), and Fungpin (FP) sections, avoiding chaotic mass-transport deposits (slump beds, olistoliths) (Supplementary Figs. 2–4). Paleomagnetic samples were collected using a standard (22 mm diameter) drill core from fresh mudstone exposures, and remanent magnetization was measured with a 2G three-axis cryogenic magnetometer. To remove viscous remanent component of overprinting magnetic signals, we applied stepwise thermal demagnetization (THD, from room temperature up to 800 °C) or alternating-field demagnetization (AFD, from 0 up to 80 mT) to most samples. We applied a combination of THD and AFD (e.g., THD to 360 °C followed by AFD procedure) in cases of some specimens that became thermally unstable at higher temperatures. Through these procedures, we obtained reliable measurements of primary remanent component of the paleomagnetic declination and inclination from Zijderveld-type diagrams at each site (Supplementary Fig. 6). We derived mean paleomagnetic directions by restoring the perturbations of regional folds (bedding dip) (Supplementary Data 1). The age values and their uncertainty ranges of magnetic reversals follow the most recent global geomagnetic polarity timescale82,83.
Calcareous nanoplanktons and planktonic foraminifera biostratigraphy
In the southern Coastal Range, we adapted results and age interpretations from previously published studies34,42,79,80 (Supplementary Fig. 1). In the northern Coastal Range, we digitized and manually georeferenced unpublished calcareous nannoplankton fossil charts (Supplementary Data 2) and sample localities (Supplementary Figs. 3, 4) from original notes and field maps of previous works by W.-R. Chi38,84 (Supplementary Fig. 2). For planktonic foraminifera biostratigraphy, we compiled published data of index fossils85,86,87,88 and collected five new samples in Fungfu (FF) and Fungpin (FP) sections (Supplementary Figs. 2, 4; Supplementary Data 2). Samples were collected from fresh intact exposures of mudstone, and we used 63 μm sieve to extract proper size foraminifera for identification.
Interpretation of depositional ages is based primarily on the first appearance datum (FAD) for index fossils due to the potential of fossil reworking, which is commonly reported in turbidite-dominated deposits of the Coastal Range34,38,89 and confirmed by our work. Age values and uncertainties are based on recent compilations for the Indo-Pacific region90,91 (Fig. 2).
Paleobathymetry
When abundance data of both benthic and planktic foraminifera are available, one can empirically estimate the paleobathymetry (W) of the marine sediment through a regression relating planktic percentage (%Ps) to modern water depth92,93. The planktic percentage (%Ps) is defined as:
$$\begin{array}{c}{ \% P}_{s}=100\cdot \frac{P}{P\,+\,B\,-\,S\,-\,O}\end{array}$$
where P is the number of in situ planktic specimens; B is the number of benthic specimens; S is the number of stress-marker benthic specimens that are more sensitive to other environmental factors (e.g., oxygen level) rather than water depth (in the genera Bolivina, Bulimina, Chilostomella, Fursenkoina, Globobulimina, Uvigerina s.s.); O is the number of organically cemented benthic specimens and reworked fossils. Once (%Ps) is calculated, the paleobathymetry (W) can be estimated through this logistic function:
$$\begin{array}{c}W={{\exp }}\left[{{\alpha }}+\frac{1}{\beta }{{{{{\rm{ln}}}}}}\left(\frac{{ \% P}_{s}}{100\,-\,{ \% P}_{s}}\right)\right]-0.1\end{array}$$
where α and \(\beta\) are empirical constants from modern settings. Despite the lack of modern offshore constraints near Taiwan, we adapted α = 1.45 ± 0.080 and \(\beta\) = 5.23 ± 0.042 (means ± standard errors) regressed from modern analog of 0–4500 m water depth in subtropical Pacific near New Zealand by Hayward et al.92.
We converted qualitative (in ordinal scale) abundance data of from Chang85,86,87 (V: > 50 specimens; A: 21-50 specimens; C: 11–20 specimens; F: 6-10 specimens; R: < 6 specimens; ?: unsurely identified; D: reworked fossils) into quantitative (ratio) scale with ranges of uncertainty (V = 75 ± 25; A = 35 ± 15; C = 15 ± 5; F = 7 ± 3; R = 3 ± 2; ? = 1 ± 1; excluding reworked fossils) to constrain parameters P, B, and S. Parameter O is assumed to be zero because we have excluded reworked fossils. We then used Eq. (2) to calculate the paleobathymetry (W) and its uncertainty (standard error) through Gaussian error propagation at each sample site of Chang78,79,80 along studied sections [Fungfu (FF), Fungpin (FP), Bieh river (BC), Madagida river (MDJ), and Sanshian river (SSS) sections] (Supplementary Data 3).
Fossil abundance data are sometimes absent [e.g., data of Hsiukuluan river (HKL) and Wulou river (WL) sections from Chang and Chen88] or paleobathymetry cannot be calculated using the planktic-percentage method (e.g., %Ps = 0 or 100). In these cases, we used an alternative method from Hohenegger94, which relies solely on the presence/absence of the benthic species and their modern water depth distributions. The basic function of this method can be written as
$$\begin{array}{c}W=\mathop{\sum }\limits_{n=1}^{m}{l}_{n}{{d}_{n}}^{-1}/\mathop{\sum }\limits_{n=1}^{m}{{d}_{n}}^{-1}\end{array}$$
where \({l}_{n}\) and \({d}_{n}\) are the location parameter (i.e., mean water depth) and its dispersion (water depth range) of the nth benthic species, respectively; m is the total amount of benthic species that are considered. This method is based on the idea that species with narrower present-day depth distribution could yield more information about the paleobathymetry of the sediment than species adapted to live in a wide range of water depth. We collected constraints of the minimum and maximum distributed water depths of different benthic species or genera (ωmin and ωmax respectively) from publications to date (see cited data in Supplementary Data 3). We then derived the geometric mean (\({l}_{n}=\sqrt{{\omega }_{{\max }}\cdot {\omega }_{{\min }}}\)) and range (\({d}_{n}={\omega }_{{\max }}-{\omega }_{{\min }}\)) of water depth for each species or genus, and calculated the paleobathymetry (W) through Eq. (3). The uncertainty of paleobathymetry (\({\sigma }_{w}\)) can be estimated as below:
$$\begin{array}{c}{\sigma }_{w}=\sqrt{\mathop{\sum }\limits_{n=1}^{m}{\left({l}_{n}-W\right)}^{2}{{d}_{n}}^{-1}/\mathop{\sum }\limits_{n=1}^{m}{{d}_{n}}^{-1}}\end{array}$$
Lastly, we projected these sample locations along the composite stratigraphic columns to determine changes in paleobathymetry through time (Supplementary Figs. 1, 2, and 7).
Age models and quality of age measurements
Age constraints (paleomagnetic reversals and first appearance datums of index fossils) and their measured present-day stratigraphic heights (before correcting effects of compaction) on the composite sections (Fig. 2) are used to conducted linear regressions between thickness and age (Supplementary Data 4). We visually determine data trends and group data to be fitted in different regression lines (Fig. 3a). These linear models are then used to predict (using extrapolation and interpolation) depositional ages of the section boundaries (base of section, top of preserved section, top of reconstructed section) and paleobathymetry samples (Supplementary Data 3–4).
The resolution of our dating with these integrated magneto-biostratigraphy methods is expected to be as good as the precision [~1–15 thousand years (kyr)] of each astronomically-calibrated age of paleomagnetic reversal83 and first appearance datum (FAD) for index fossils90, once they are accurately placed in the stratigraphic sections. The stratigraphic position of each age-control datum was determined using the two sample sites that define a magnetic polarity reversal or the first appearance of an index fossil. The stratigraphic thickness between two bounding paleomagnetic sample sites mostly lies between 8 to 288 m, except a poor constraint (828 m) at reversal between polarity chrons C2An1.n and C1r.2r (Gauss-Matuyama boundary) near the base of Fungfu (FF) section in the northern Coastal Range (Fig. 2; Supplementary Figs. 1, 2). Using the mean sediment-accumulation rates calculated from height-age linear regressions, the uncertainty (accuracy) for most of our age measurements is estimated to be ± 2–43 kyr (0.1–2%) from targeted age values [±~507 kyr (~10%) for the Gauss–Matuyama boundary in the FF section] (Supplementary Data 4). This demonstrates remarkably high-resolution and high-fidelity controls on depositional ages and associated vertical subsidence and uplift rates.
Subsidence analysis: decompaction and backstripping
We used established numerical methods of basin geohistory analysis45 that correct for loss of pore space during progressive burial and compaction of sediment through time to reconstruct of subsidence history of the sedimentary basin. Stratigraphic thickness, corrected for the effects of compaction and sea level, is therefore used to track the depth of the basement top through time95. This approach assumes that the pore spaces are interconnected (i.e., no overpressure), and the sediment porosity [\({{\varnothing }}({{{{{\rm{z}}}}}})\)] varies as an exponential function of depth (z):
$$\begin{array}{c}{{\varnothing }}\left(z\right)={{\varnothing }}\left(0\right)\cdot {{{{{{\rm{e}}}}}}}^{-z/C}\end{array}$$
where C is a decay constant that varies with lithology. We calculated global average porosity-depth functions from previous publications95,96 to estimate initial porosity [\({{\varnothing }}\left(0\right)\)], decay constant (C), and uncertainties (standard errors) for sandy, muddy, and mean marine sediments (Supplementary Fig. 8; Supplementary Data 4). We assumed that basement rocks below the basal unconformity (i.e., Kangkou Limestone, Tuluanshan Formation) have not compacted. This assumption is supported by the presence of calcite-cemented Miocene volcaniclastic sandstones beneath the unconformity that are directly overlain by uncemented mudstones and sandstones above the contact33,34, showing that rocks beneath the unconformity were cemented prior to deep post-Miocene basin subsidence. Rocks beneath the unconformity may be subject to minor, inconsiderable amounts of compaction that do not affect the results of our subsidence analysis.
For N stratigraphic units, the present-day thickness (\({T}_{0}\)) of each unit was buried at a depth of \({D}_{0}\) when the youngest sediment deposited prior to structural tilting and erosion. There are 14 preserved units (N = 14) in both northern and southern Coastal Range. Stratigraphic positions of all age boundaries are presented in Supplementary Data 4. Uncertainties of \({T}_{0}\) of preserved units (unit 1–14) are limited by the stratal thickness between the two bounding samples that confine the magnetic reversals or the first appearance datums (FAD) of index fossils (Supplementary Figs. 1,2; Supplementary Data 4). The uppermost preserved sediment at the top of unit 1 is consolidated and lithified, which means that there must have been a considerable thickness of strata (called unit 0 in this study) above preserved unit 1 that was subsequently removed by erosion during post-0.5 Ma uplift. Porosity-effective stress of the sandstone44 and compiled vitrinite reflectance97 data yield an estimated range of eroded thickness of ca. 2.2–3.7 km above marker bed tuff Tp12 near the top of the Madagida river (MDJ) section (Supplementary Fig. 1). This implies ca. 0.45–1.95 km sediment that was deposited above unit 1 and has since been removed. For simplicity, we assumed a conservative thickness range of 0.75 ± 0.25 km for the unpreserved unit 0 in both the north and south composite sections, and we used it to reconstruct the top of section prior to structural inversion (Fig. 3; Supplementary Data 4).
To calculate the thickness of each unit at some earlier time (i.e., decompacted thickness, \({T}_{i}\)), when the unit was buried only to a depth of \({D}_{i}\), we can use the mass-balance equation:
$$\begin{array}{c}{\int }_{{D}_{i}}^{{{D}_{i}+{T}_{i}}}\left[1-{{{\varnothing }}}_{i}\left(z\right)\right]{{{{{{\rm{d}}}}}}z}=\,{\int }_{{D}_{0}}^{{{D}_{0}+{T}_{0}}}\left[1-{{{\varnothing }}}_{i}\left(z\right)\right]{{{{{{\rm{d}}}}}}z}\end{array}$$
where the \({{{\varnothing }}}_{i}\left(z\right)\) is the porosity of unit i (i = 0, 1, 2,…, N) at depth z. This approach assumes the volume of sediment grains within the unit does not change. The quantity \(\left[1-{{{\varnothing }}}_{i}\left(z\right)\right]\) represents the volume of sediment grains (per unit volume of the strata) at any level within the unit. After integrating Eq. (5), Eq. (6) is rearranged to iteratively solve \({T}_{i}\):
$$\begin{array}{c}\Longrightarrow {T}_{i}=\mathop{\sum }\limits_{k\,=\,1}^{M}{F}_{i}^{k}\cdot \left[{\varphi }_{1}\cdot {{{{{{\rm{e}}}}}}}^{-{T}_{i}/{C}_{i}^{k}}+{\varphi }_{2}\right]\end{array}$$
with \(\left\{\begin{array}{c}{\varphi }_{1}=-{C}_{i}^{k}\cdot {{{\varnothing }}}_{i}^{k}\left(0\right)\cdot {{{{{{\rm{e}}}}}}}^{-{D}_{i}/{C}_{i}^{k}}\hfill\\ {\varphi }_{2}=-{\varphi }_{1}+{T}_{0}+{C}_{i}^{k}\cdot {{{\varnothing }}}_{i}^{k}\left(0\right)\cdot {{{{{{\rm{e}}}}}}}^{-{D}_{0}/{C}_{i}^{k}}\cdot \left[{{{{{{\rm{e}}}}}}}^{-{D}_{0}/{C}_{i}^{k}}-1\right]\end{array}\right.\)where \({F}_{i}^{k}\) is the fraction of lithology k (k = 1, 2,…, M) in unit i. There are three types of lithology (M = 3) determined in this study. The fraction of sand (k = 1), mud (k = 2), and pebbly mudstone (k = 3) was measured in all sections (Supplementary Figs. 1, 2; Supplementary Data 4). We assumed that unit 0 contains subequal fractions of sand and mud (\({F}_{0}^{1}=0.5\); \({F}_{0}^{2}=0.5\); \({F}_{0}^{3}=0\)). The initial porosity [\({\varnothing }_{i}^{k}\)(0)] and decay constant (\({C}_{i}^{k}\)) for sandy, muddy, and mean marine sediments were applied respectively (Supplementary Fig. 8). After repeating Eqs. (6) and (7) for each unit from top to the base of section, we computed the sum of decompacted unit thickness (\({\epsilon }_{i}\)) at the time when unit i finished deposition:
$$\begin{array}{c}{\epsilon }_{i}=\mathop{\sum }\limits_{i\,=\,0}^{N}{T}_{i}\end{array}$$
These procedures [Eqs. (6–8)] account for the effects of removing incrementally older units from the top and allow the section to expand as underlying units are unloaded.
To evaluate the amounts of total subsidence (i.e., the true depth of the basement, \({{{{{{\rm{\sigma }}}}}}}_{i}\)) at the end of unit i deposition, we need to correct the magnitude of paleo-water depth by assuming the basin was filled with water to sea level:
$$\begin{array}{c}{{{{{{\rm{\sigma }}}}}}}_{i}={\epsilon }_{i}+{W}_{i}-{\Delta }_{i}\end{array}$$
where \({W}_{i}\) and \({\Delta }_{i}\) are the paleobathymetry and sea level relative to modern datum at the time of unit i finishing deposition. We combined means and standard errors of available paleobathymetry estimates (Supplementary Data 3) with reconstructed global sea level from Miller et al.98 using age constraints in each section (within the range of age uncertainty with extra ± 0.2 Myr and ± 0.01 Myr for paleobathymetry and relative sea level respectively) as the values of \({W}_{i}\) and \({\Delta }_{i}\) at each unit boundary (Supplementary Fig. 7; Supplementary Data 4). Paleobathymetry constraint was not available for the reconstructed unit 0 above the preserved strata. In this case, we assumed the paleobathymetry of reconstructed basin top (at the end of unit 0 deposition, \({W}_{0}\)) and its uncertainty range as a half of the \({W}_{1}\) value at the top of preserved strata [\({W}_{0}={W}_{1}/2\pm {W}_{1}/2\)] because basin inversion in the Coastal Range possibly began between the ends of unit 1 and unit 0 (Fig. 3; Supplementary Fig. 7).
Since the folded basal unconformity and underlying volcanic basement rocks are exposed well above sea level (Fig. 4c), we can calculate the minimum and the best-estimated long-term rock exhumation rates using the decompacted depth of basement at depositional ends of unit 1 and unit 0 (\({{{{{{\rm{\sigma }}}}}}}_{1}\) and \({{{{{{\rm{\sigma }}}}}}}_{0}\)) to the present-day elevations of the antiformal peaks in northern and southern Coastal Range (0–1.33 km and 0–1.68 km, respectively) (Fig. 3).
We further estimated the amount of subsidence that was contributed by tectonic forces (so-called "tectonic subsidence," \({\zeta }_{i}\)) at the end of unit i deposition by removing effects of paleo-water depth and local (Airy) isostatic response to applied sediment loads:
$$\begin{array}{c}{\zeta }_{i}={\epsilon }_{i}\frac{\left({\rho }_{a}-{\bar{{\rho }_{s}}}^{i}\right)}{\left({\rho }_{a}-{\rho }_{w}\right)}+{W}_{i}-{\Delta }_{i}\frac{{\rho }_{a}}{\left({\rho }_{a}-{\rho }_{w}\right)}\end{array}$$
where \({\rho }_{w}\) is the density of marine water (\({\rho }_{w}=1025\,{{\mbox{kg}}}\,{{{\mbox{m}}}}^{{{\mbox{-}}}3}\)); \({\rho }_{a}\) is the density of asthenospheric mantle (\({\rho }_{a}=3300\,{{\mbox{kg}}}\,{{{\mbox{m}}}}^{{{\mbox{-}}}3}\)). \({\bar{{\rho }_{s}}}^{i}\) is the average (integral) density of the unit i with various lithologies k at each time frame of sedimentation:
$$\begin{array}{c}{\bar{{\rho }_{s}}}^{i}=\left(\mathop{\sum }\limits_{i\,=\,0}^{N}\mathop{\sum }\limits_{k\,=\,1}^{M}{F}_{i}^{k}\cdot {\Omega }_{i}^{k}\right)/{\epsilon }_{i}\end{array}$$
\({\Omega }_{i}^{k}\) represents the average (integral) 1-D mass of unit i for given lithology k at some earlier time, assuming the pore spaces were filled by marine water,
$${\Omega }_{i}^{k} ={{\int}_{{D}_{i}}^{{{{D}_{i}+{T}_{i}}}}}\left\{{{{\varnothing }}}_{i}^{k}\left(z\right)\cdot {\rho }_{w}+\left[1-{{{\varnothing }}}_{i}^{k}\left(z\right)\right]\cdot {({\rho }_{g})}_{k}\right\}{{{{{{\rm{d}}}}}}z}\\ ={({\rho }_{g})}_{k}\cdot {T}_{i}+\left[{\rho }_{w}-{({\rho }_{g})}_{k}\right]\cdot {C}_{i}^{k}\cdot {{{\varnothing }}}_{i}^{k}\left(0\right)\cdot {{{{{{\rm{e}}}}}}}^{-{D}_{i}/{C}_{i}^{k}}\cdot \left[1-{{{{{{\rm{e}}}}}}}^{-{T}_{i}/{C}_{i}^{k}}\right]$$
where \({({\rho }_{g})}_{k}\) represents the averaged density of sediment grain for each lithology k in unit i. We applied the means and standard errors of the grain densities for sandy, muddy, and mean marine sediments (for pebbly mudstone) respectively from previously published data95 (Supplementary Fig. 8; Supplementary Data 4).
Uncertainties of all input parameters were considered in the analysis. Thus, we were able to estimate the uncertainties of calculated tectonic subsidence (\({\zeta }_{i}\)) and decompacted total subsidence (\({{{{{{\rm{\sigma }}}}}}}_{i}\)) through Gaussian error propagation or delta method (Supplementary Data 4). Throughout the analysis, we found that the primary source of error is the present-day (non-decompacted) thickness of each unit (\({T}_{0}\)), followed by uncertainties in paleobathymetry estimates (\({W}_{i}\)).
Lastly, we plotted the subsidence history (i.e., geohistory) diagrams for northern and southern composite sections, which show changes in depth of basement [i.e., decompacted total subsidence (\({{{{{{\rm{\sigma }}}}}}}_{i}\))] and calculated corresponding amount of tectonic subsidence (\({\zeta }_{i}\)) through time (Fig. 3b, c). The amount of total subsidence prior to decompaction correction was also plotted as a standard convention. The position of basement top at the beginning of subsidence for each composite section was reconstructed by projecting along the oldest segment of subsidence curve to the sea level (0 km). This indicates the end of formation of the basal unconformity6,33, which is characterized by a broad erosive surface that formed near sea level and is capped by discontinuous thin (ca. 50–200 m) deposits of shallow-marine limestone (Kangkou Limestone, from 5.57–4.37 Ma to 4.31–3.47 Ma) and limestone-clast bearing epiclastic rocks (Biehchi Epiclastic Unit, from 5.53–4.31 Ma to 3.82–3.35 Ma)34,35 (Supplementary Figs. 1, 2).
Data presented in this paper are provided in Supplementary Data 1–4 (Microsoft Excel spreadsheets) and permanently stored at https://doi.org/10.6084/m9.figshare.19350530.
Code availability
Analyses and production of figures were conducted using R version 4.0.3. All scripts and data in required formats are freely available at https://doi.org/10.5281/zenodo.5823613.
Champagnac, J.-D., Molnar, P., Sue, C. & Herman, F. Tectonics, climate, and mountain topography. J. Geophys. Res.: Solid Earth 117, https://doi.org/10.1029/2011JB008348 (2012).
Turner, J. P. & Williams, G. A. Sedimentary basin inversion and intra-plate shortening. Earth-Sci. Rev. 65, 277–304 (2004).
Taylor, M., Forte, A. M., Laskowski, A. & Ding, L. Active uplift of southern Tibet revealed. GSA Today 31, 4–10 (2021).
Lundberg, N. & Dorsey, R. J. Rapid Quaternary emergence, uplift, and denudation of the Coastal Range, eastern Taiwan. Geology 18, 638–641 (1990).
Sandmann, S., Nagel, T. J., Froitzheim, N., Ustaszewski, K. & Münker, C. Late Miocene to Early Pliocene blueschist from Taiwan and its exhumation via forearc extraction. Terr. Nova 27, 285–291 (2015).
Lundberg, N. & Dorsey, R. J. in New Perspectives in Basin Analysis Frontiers in Sedimentary Geology Ch. 13, (eds. Kleinspehn, K. & Paola, C.) 265–280 (Springer, 1988).
Dorsey, R. J., Housen, B. A., Janecke, S. U., Fanning, C. M. & Spears, A. L. F. Stratigraphic record of basin development within the San Andreas fault system: Late Cenozoic Fish Creek–Vallecito basin, southern California. Geol. Soc. Am. Bull. 123, 771–793 (2011).
Spotila, J. A., Farley, K. A., Yule, J. D. & Reiners, P. W. Near-field transpressive deformation along the San Andreas fault zone in southern California, based on exhumation constrained by (U-Th)/He dating. J. Geophys. Res.: Solid Earth 106, 30909–30922 (2001).
Öğretmen, N. et al. Evidence for 1.5 km of uplift of the central Anatolian plateau's southern margin in the last 450 kyr and Implications for Its multiphased uplift history. Tectonics 37, 359–390 (2018).
Baldwin, S. L. et al. Pliocene eclogite exhumation at plate tectonic rates in eastern Papua New Guinea. Nature 431, 263–267 (2004).
Watts, A. B. & Cochran, J. R. Gravity anomalies and flexure of the lithosphere along the Hawaiian-Emperor seamount chain. Geophys. J. Int. 38, 119–141 (1974).
Barletta, V. R. et al. Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability. Science 360, 1335–1339 (2018).
Willett, S. D. Orogeny and orography: The effects of erosion on the structure of mountain belts. J. Geophys. Res.: Solid Earth 104, 28957–28981 (1999).
Lister, G. & Forster, M. Tectonic mode switches and the nature of orogenesis. Lithos 113, 274–291 (2009).
Gingerich, P. D. Rates of geological processes. Earth-Sci. Rev. 220, 103723 (2021).
Reiners, P. W. & Brandon, M. T. Using thermochronology to understand orogenic erosion. Annu. Rev. Earth Planet. Sci. 34, 419–466 (2006).
Brown, M. The contribution of metamorphic petrology to understanding lithosphere evolution and geodynamics. Geosci. Frontiers 5, 553–569 (2014).
Braun, J., van der Beek, P. & Batt, G. Quantitative Thermochronology: Numerical Methods for the Interpretation of Thermochronological Data (Cambridge University Press, 2006).
Byrne, T. B. et al. in Arc-Continent Collision Frontiers in Earth Sciences Ch. 8 (eds. Brown, D. & Ryan, P. D.) 213-245 (Springer, 2011), https://doi.org/10.1007/978-3-540-88558-0_8.
Yu, S.-B., Chen, H.-Y. & Kuo, L.-C. Velocity field of GPS stations in the Taiwan area. Tectonophysics 274, 41–59 (1997).
Dadson, S. J. et al. Links between erosion, runoff variability and seismicity in the Taiwan orogen. Nature 426, 648–651 (2003).
Suppe, J. Mechanics of mountain building and metamorphism in Taiwan. Mem. Geol. Soc. China 4, 67–89 (1981).
Teng, L. S. Geotectonic evolution of late Cenozoic arc-continent collision in Taiwan. Tectonophysics 183, 57–76 (1990).
Huang, C.-Y., Yuan, P. B. & Tsao, S.-J. Temporal and spatial records of active arc-continent collision in Taiwan: A synthesis. Geol. Soc. Am. Bull. 118, 274–288 (2006).
Lee, Y.-H. et al. Simultaneous mountain building in the Taiwan orogenic belt. Geology 43, 451–454 (2015).
Hsu, W.-H. et al. Pleistocene onset of rapid, punctuated exhumation in the eastern Central Range of the Taiwan orogenic belt. Geology 44, 719–722 (2016).
Wu, J., Suppe, J., Lu, R. & Kanda, R. Philippine Sea and East Asian plate tectonics since 52 Ma constrained by new subducted slab reconstruction methods. J. Geophys. Res.: Solid Earth 121, 4670–4741 (2016).
Sibuet, J.-C., Zhao, M., Wu, J. & Lee, C.-S. Geodynamic and plate kinematic context of South China Sea subduction during Okinawa trough opening and Taiwan orogeny. Tectonophysics 817, 229050 (2021).
Chen, W.-S. et al. A reinterpretation of the metamorphic Yuli belt: Evidence for a middle-late Miocene accretionary prism in eastern Taiwan. Tectonics 36, 188–206 (2017).
Mesalles, L. et al. A Late-Miocene Yuli belt? New constraints on the eastern Central Range depositional ages. Terrestrial Atmos. Ocean. Sci. 31, 403–414 (2020).
Lo, Y.-C., Chen, C.-T., Lo, C.-H. & Chung, S.-L. Ages of ophiolitic rocks along plate suture in Taiwan orogen: Fate of the South China Sea from subduction to collision. Terrestrial Atmos. Ocean. Sci. 31, 383–402 (2020).
Yui, T.-F. et al. Dating thin zircon rims by NanoSIMS: the Fengtien nephrite (Taiwan) is the youngest jade on Earth. Int. Geol. Rev. 56, 1932–1944 (2014).
Dorsey, R. J. Collapse of the Luzon volcanic arc during onset of arc-continent collision: evidence from a Miocene-Pliocene unconformity, eastern Taiwan. Tectonics 11, 177–191 (1992).
Lai, L. S.-H. et al. Polygenetic mélange in the retrowedge foredeep of an active arc-continent collision, Coastal Range of eastern Taiwan. Sedimentary Geol. 418, 105901 (2021).
Huang, C.-Y. & Yuan, P. B. Stratigraphy of the Kangkou Limestone in the Coastal Range, eastern Taiwan. J. Geol. Soc. China 37, 585–605 (1994).
Dorsey, R. J., Buchovecky, E. J. & Lundberg, N. Clay mineralogy of Pliocene-Pleistocene mudstones, eastern Taiwan: Combined effects of burial diagenesis and provenance unroofing. Geology 16, 944–947 (1988).
Kirstein, L. A., Carter, A. & Chen, Y.-G. Impacts of arc collision on small orogens: new insights from the Coastal Range detrital record, Taiwan. J. Geol. Soc. London 171, 5–8 (2014).
Chi, W.-R., Namson, J. & Suppe, J. Stratigraphic record of plate interactions in the Coastal Range of eastern Taiwan. Mem. Geol. Soc. China 4, 155–194 (1981).
Chang, S. S. L. & Chi, W.-R. Neogene nannoplankton biostratigraphy in Taiwan and the tectonic implications. Petroleum Geol. Taiwan 19, 93–147 (1983).
Lee, T.-Q. Study of the polarity transition record of the upper Olduvai event from Wulochi sedimentary sequence of the Coastal Range, eastern Taiwan. Terrestrial Atmos. Ocean. Sci. 3, 503–518 (1992).
Lee, T.-Q. Evolution tectonique et geodynamique neogene et quaternaire de la chaine cotiere de Taiwan: apport du paleomagnetisme Doctorate thesis, Universite Pierre et Marie Curie (1989).
Horng, C.-S. & Shea, K.-S. Dating of the Plio-Pleistocene rapidly deposited sequence based on integrated magneto-biostratigraphy: a case study of the madagida-chi section, Coastal Range, eastern Taiwan. J. Geol. Soc. China 39, 31–58 (1996).
Chen, W.-H. et al. Depleted deep South China Sea δ13C paleoceanographic events in response to tectonic evolution in Taiwan–Luzon Strait since Middle Miocene. Deep Sea Res. II: Topical Stud. Oceanogr. 122, 195–225 (2015).
Hong, S.-M. Using Porosity-effective Stress Relationship Curve to Evaluate the Erosion Amount of the Fore-arc Basin in the Southern Part of Coastal Range Master thesis, National Central University (2020).
Allen, P. A. & Allen, J. R. Basin Analysis: Principles and Application to Petroleum Play Assessment 3rd edn (John Wiley & Sons, Ltd., 2013).
Chen, W.-S., Yang, C.-Y., Chen, S.-T. & Huang, Y.-C. New insights into Holocene marine terrace development caused by seismic and aseismic faulting in the Coastal Range, eastern Taiwan. Quaternary Sci. Rev. 240, 106369 (2020).
Hsieh, M.-L. & Rau, R.-J. Late Holocene coseismic uplift on the Hua-tung coast, eastern Taiwan: evidence from mass mortality of intertidal organisms. Tectonophysics 474, 595–609 (2009).
Shyu, J. B. H., Sieh, K., Avouac, J.-P., Chen, W.-S. & Chen, Y.-G. Millennial slip rate of the Longitudinal Valley fault from river terraces: Implications for convergence across the active suture of eastern Taiwan. J. Geophys. Res. 111, B08403 (2006).
Ching, K.-E. et al. Modern vertical deformation rates and mountain building in Taiwan from precise leveling and continuous GPS observations, 2000–2008. J. Geophys. Res.: Solid Earth 116, B08406 (2011).
Page, B. M. & Suppe, J. The Pliocene Lichi Mélange of Taiwan: its plate-tectonic and olistostromal origin. Am. J. Sci. 281, 193–227 (1981).
Lin, A. T. & Watts, A. B. Origin of the West Taiwan basin by orogenic loading and flexure of a rifted continental margin. J. Geophys. Res.: Solid Earth 107, ETG 2-1–ETG 2-19 (2002).
Barrier, E. & Angelier, J. Active collision In eastern Taiwan: the Coastal Range. Mem. Geol. Soc. China 7, 135–159 (1986).
Huang, C.-Y. et al. Tectonics of short-lived intra-arc basins in the arc-continent collision terrane of the Coastal Range, eastern Taiwan. Tectonics 14, 19–38 (1995).
Di Martire, D., Ascione, A., Calcaterra, D., Pappalardo, G. & Mazzoli, S. Quaternary deformation in SE Sicily: Insights into the life and cycles of forebulge fault systems. Lithosphere 7, 519–534 (2015).
Mesalles, L. et al. From submarine continental accretion to arc-continent orogenic evolution: The thermal record in southern Taiwan. Geology 42, 907–910 (2014).
Wang, W.-H. & Lee, Y.-H. 3-D plate interactions in central Taiwan: Insight from flexure and sandbox modeling. Earth. Planet. Science Lett. 308, 1–10 (2011).
Hsieh, Y.-H., Liu, C.-S., Suppe, J., Byrne, T. B. & Lallemand, S. The Chimei submarine canyon and fan: A record of Taiwan arc-continent collision on the rapidly deforming over-riding plate. Tectonics 39, e2020TC006148 (2020).
Watts, A. B. Isostasy and Flexure of the Lithosphere (Cambridge University Press, 2001).
Shyu, J. B. H., Wu, Y.-M., Chang, C.-H. & Huang, H.-H. Tectonic erosion and the removal of forearc lithosphere during arc-continent collision: Evidence from recent earthquake sequences and tomography results in eastern Taiwan. J. Asian Earth Sci. 42, 415–422 (2011).
Chemenda, A. I., Yang, R. K., Stephan, J. F., Konstantinovskaya, E. A. & Ivanov, G. M. New results from physical modelling of arc–continent collision in Taiwan: evolutionary model. Tectonophysics 333, 159–178 (2001).
Malavieille, J. et al. in Geology and Geophysics of an Arc-continent Collision, Taiwan Vol. 358 187-211, https://doi.org/10.1130/0-8137-2358-2.187 (Geological Society of America, 2002).
Froitzheim, N., Pleuger, J. & Nagel, T. J. Extraction faults. J. Struct. Geol. 28, 1388–1395 (2006).
Majka, J. et al. Microdiamond discovered in the Seve Nappe (Scandinavian Caledonides) and its exhumation by the "vacuum-cleaner" mechanism. Geology 42, 1107–1110 (2014).
Zhang, Y., Tsai, C.-H., Froitzheim, N. & Ustaszewski, K. The Yuli Belt in Taiwan: Part of the suture zone separating Eurasian and Philippine Sea plates. Terrestrial Atmos. Ocean. Sci. 31, 415–435 (2020).
Abbott, L. D. et al. Measurement of tectonic surface uplift rate in a young collisional mountain belt. Nature 385, 501–507 (1997).
Molnar, P. Isostasy can't be ignored. Nat. Geoscience 5, 83–83 (2012).
Teng, L. S., Lee, C. T., Tsai, Y. B. & Hsiao, L. Y. Slab breakoff as a mechanism for flipping of subduction polarity in Taiwan. Geology 28, 155–158 (2000).
Yen, J. Y. et al. Insights into seismogenic deformation during the 2018 Hualien, Taiwan, earthquake sequence from InSAR, GPS, and modeling. Seismol. Res. Lett. 90, 78–87 (2018).
Draut, A. E. & Clift, P. D. Differential preservation in the geologic record of intraoceanic arc sedimentary and tectonic processes. Earth-Sci. Rev. 116, 57–84 (2013).
Allen, P. A. Time scales of tectonic landscapes and their sediment routing systems. Geol. Soc. London, Special Publ. 296, 7–28 (2008).
Albarède, F. The growth of continental crust. Tectonophysics 296, 1–14 (1998).
Clift, P. D., Schouten, H. & Draut, A. E. A general model of arc-continent collision and subduction polarity reversal from Taiwan and the Irish Caledonides. Geol Soc London Special Publications 219, 81–98 (2003).
Romans, B. W., Castelltort, S., Covault, J. A., Fildani, A. & Walsh, J. P. Environmental signal propagation in sedimentary systems across timescales. Earth-Sci Rev 153, 7–29 (2016).
Hsu, T. L. Geology of the Coastal Range, eastern Taiwan. Bull Central Geol Surv. Taiwan 8, 39–63 (1956).
Wang, Y. & Chen, W.-S. Geological map of eastern Coastal Range (1:100,000). Central Geological Survey, Ministry of Economic Affairs of Taiwan, New Taipei, Taiwan, 1993.
Chen, W.-S. & Wang, Y. Geological map of Taiwan - Fengpin sheet (1:50,000). Central Geological Survey, Ministry of Economic Affairs of Taiwan, New Taipei, Taiwan, 1997.
Yi, D.-C., Chen, C.-Y. & Lin, C.-W. Geological map of Taiwan - Guangfu sheet (1:50,000). Central Geological Survey, Ministry of Economic Affairs of Taiwan, New Taipei, Taiwan, 2012.
Teng, L. S., Tsai, Y.-l & Kuo, S.-T. On the Chimei Fault zone of the Coastal Range, eastern Taiwan. Bull. Central Geol. Surv. Ministry. Econom. Affairs. Taiwan 29, 1–44 (2016).
Lai, L. S.-H. & Teng, L. S. Stratigraphy and structure of the Tai-Yuan basin, southern Coastal Range, eastern Taiwan. Bull. Central Geol. Surv. Ministry. Econom. Affairs.Taiwan 29, 45–76 (2016).
Lai, L. S.-H., Ng, T.-W. & Teng, L. S. Stratigraphic correlation of tuffaceous and psephitic strata in the Paliwan formation, southern Coastal Range of eastern Taiwan. Bull. Central Geol. Surv. Ministry. Econom. Affairs Taiwan 31, 1–32 (2018).
Hsu, C.-W., Liu, Y.-C. & Yen, Y.-C. A study on the fault trace and the recent drilling data of the Chimei Fault, eastern Taiwan. Special Publication. Central Geol. Surv. Ministry. Econom. Affairs Taiwan 31, 91–116 (2017).
Ogg, J. G. in Geologic Time Scale 2020 Ch. 5 (eds Gradstein, F. M. Ogg, J. G., Schmitz, M. D. & Ogg G. M.) 159–192 (Elsevier, 2020) https://doi.org/10.1016/B978-0-12-824360-2.00005-X.
Channell, J. E. T., Singer, B. S. & Jicha, B. R. Timing of Quaternary geomagnetic reversals and excursions in volcanic and sedimentary archives. Quaternary Sci. Rev. 228, 106114 (2020).
Chi, W.-R., Namson, J. & Mei, W. W. Calcareous nannoplankton biostratigraphy of the late Neogene Sediments exposed along the Hsiukuluanchi In the Coastal Range, eastern Taiwan. Petroleum Geol. Taiwan 17, 75–87 (1980).
Chang, L.-S. A biostratigraphic study of the Tertiary in the Coastal Range, eastern Taiwan, based on smaller foraminifera (I: southern part). Proc. Geol. Soc. China 10, 64–76 (1967).
Chang, L.-S. A biostratigraphic study of the Tertiary in the Coastal Range, eastern Taiwan, based on smaller foraminifera (III: middle part). Proc. Geolog. Soc. China 12, 89–101 (1969).
Chang, L.-S. A biostratigraphic study of the Tertiary in the Coastal Range, eastern Taiwan, based on smaller foraminifera (II: northern part). Proc. Geol. Soc. China 11, 19–33 (1968).
Chang, L.-S. & Chen, T. H. A biostratigraphic study of the Tertiary along the Hsiukuluanchi in the Coastal Range, eastern Taiwan, based on smaller foraminifera. Proc. Geol. Soc. China 13, 115–128 (1970).
Chen, W.-S. Tectonostratigraphic framework and age of the volcanic-arc and collision basins in the Coastal Range, eastern Taiwan. Western Pac. Earth Sci. 9, 67–98 (2009).
Chuang, C.-K. et al. Integrated stratigraphy of ODP Site 1115 (Solomon Sea, southwestern equatorial Pacific) over the past 3.2 Ma. Marine Micropaleontol. 144, 25–37 (2018).
Raffi, I. et al. in Geologic Time Scale 2020 (eds Felix M. Gradstein, James G. Ogg, Mark D. Schmitz, & Gabi M. Ogg) Ch. 29, 1141–1215 (Elsevier, 2020), https://doi.org/10.1016/B978-0-12-824360-2.00029-2.
Hayward, B. W. & Triggs, C. M. Using multi-foraminiferal-proxies to resolve the paleogeographic history of a lower Miocene, subduction-related, sedimentary basin (Waitemata basin, New Zealand). J. Foraminiferal Res. 46, 285–313 (2016).
van Hinsbergen, D. J. J., Kouwenhoven, T. J. & van der Zwaan, G. J. Paleobathymetry in the backstripping procedure: Correction for oxygenation effects on depth estimates. Palaeogeogr. Palaeoclimatol. Palaeoecol. 221, 245–265 (2005).
Hohenegger, J. Estimation of environmental paleogradient values based on presence/absence data: a case study using benthic foraminifera for paleodepth estimation. Palaeogeogr. Palaeoclimatol. Palaeoecol. 217, 115–130 (2005).
Lee, E. Y., Novotny, J. & Wagreich, M. Subsidence Analysis and Visualization: For Sedimentary Basin Analysis and Modelling (Springer International Publishing, 2019).
Kominz, M. A., Patterson, K. & Odette, D. Lithology dependence of porosity In slope and deep marine sediments. J. Sediment. Res. 81, 730–742 (2011).
Chim, L. K., Yen, J.-Y., Huang, S.-Y., Liou, Y.-S. & Tsai, L. L.-Y. Using Raman spectroscopy of carbonaceous materials to track exhumation of an active orogenic belt: an example from eastern Taiwan. J. Asian Earth Sci. 164, 248–259 (2018).
Miller, K. G. et al. Cenozoic sea-level and cryospheric evolution from deep-sea geochemical and continental margin records. Sci. Adv. 6, eaaz1346 (2020).
Lai, L. S.-H., Roering, J. J., Finnegan, N. J., Dorsey, R. J. & Yen, J.-Y. Coarse sediment supply sets the slope of bedrock channels in rapidly uplifting terrain: field and topographic evidence from eastern Taiwan. Earth Surf. Processes Landforms 46, 2671–2689 (2021).
This manuscript was greatly improved by two valuable reviews. The authors also appreciate feedback from two anonymous reviewers on early versions of this manuscript. We thank Li Lo, Edward Davis, Jhih-Jyun Zeng, and Kevin Gardner for discussions that improved the quality of this manuscript and data analysis. We also thank Chien-Hao Wang, Kuo-Hang Chen, Chun-Hung Lin for field and laboratory assistance. This research was funded by the Geological Society of America (2018 Graduate Student Research Grant) and Ministry of Education of Taiwan (2019 Government Scholarship for Study Abroad) to L.S.-H.L. and the Ministry of Science and Technology of Taiwan to C.-S.H. (NSC102-2116-M-001-005), J.-Y.Y., and R.J.D (MOST107-2811-M-259-002; MOST108-2116-M-259-005-MY2).
Department of Earth Sciences, University of Oregon, Eugene, OR, 97403, USA
Larry Syu-Heng Lai & Rebecca J. Dorsey
Institute of Earth Sciences, Academia Sinica, Taipei, 11529, Taiwan
Chorng-Shern Horng
Department of Earth Sciences, National Cheng Kung University, Tainan, 70101, Taiwan
Wen-Rong Chi
Department of Earth Sciences, National Central University, Taoyuan, 32001, Taiwan
Central Geological Survey, Ministry of Economic Affairs, New Taipei, 235055, Taiwan
Kai-Shuan Shea
Department of Natural Resources and Environmental Studies, National Dong Hwa University, Hualien, 97401, Taiwan
Jiun-Yee Yen
Larry Syu-Heng Lai
Rebecca J. Dorsey
L.S.-H.L., R.J.D., J.-Y.Y. conceptualized the study and conducted geological surveys. C.-S.H. and L.S.-H.L. were responsible for paleomagnetic analysis. W.-R.C. and K.-S.S. performed microfossil identifications. All authors contributed to the manuscript and figure preparation.
Correspondence to Larry Syu-Heng Lai.
Communications Earth & Environment thanks Kamil Ustaszewski and Yuan-Hsi Lee for their contribution to the peer review of this work. Primary Handling Editors: Derya Gürer, Joe Aslin, Clare Davis and Heike Langenberg. Peer reviewer reports are available.
Description of Additional Supplementary Files
Supplementary Data 1
Peer Review File
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Lai, L.SH., Dorsey, R.J., Horng, CS. et al. Extremely rapid up-and-down motions of island arc crust during arc-continent collision. Commun Earth Environ 3, 100 (2022). https://doi.org/10.1038/s43247-022-00429-2
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Communications Earth & Environment (Commun Earth Environ) ISSN 2662-4435 (online)