Martine Simoes

Late Cenozoic metamorphic evolution and exhumation of Taiwan

The Taiwan mountain belt is composed of a Cenozoic slate belt (Hsuehshan Range units, HR, and Backbone Slates, BS) and of accreted polymetamorphic basement rocks (Tananao Complex, TC). Ongoing crustal shortening has resulted from the collision between the Chinese continental margin and the Luzon volcanic arc, which initiated ~6.5 Ma ago. The grade and age of metamorphism and exhumation are a key record of the development of the orogenic wedge. Because the Taiwan mountain belt is mostly composed by accreted sediments lacking metamorphic index minerals, quantitative constraints on metamorphism are sparse. By contrast, these rocks are rich in carbonaceaous material (CM) and are therefore particularly appropriate for RSCM (Raman Spectroscopy of CM) thermometry. We apply this technique in addition to (U-Th)/He thermochronology on detrital zircons to assess peak metamorphic temperatures (T) and the late exhumational history respectively, along different transects in central and southern Taiwan. In the case of the HR units, we find evidence for high metamorphic T of at least 340°–350°C and locally up to 475°C, and for relative rapid exhumation with zircon (U-Th)/He ages in the range of 1.5–2 Ma. Farther east, the BS were only slightly metamorphosed (T < 330 °C), and zircons are not reset for (U-Th)/He. From the eastern BS to the inner TC schists, T gradually increases from ~350°C up to ~500°C following an inverted metamorphic gradient. Available geochronological constraints and the continuous thermal gradient from the BS to the basement rocks of the TC suggest that the high RSCM T of the TC were most probably acquired during the last orogeny, and were not inherited from a previous thermal event. Zircons yield (U-Th)/He ages of ~0.5–1.2 Ma. Peak metamorphic T and the timing of exhumation do not show along-strike variations over the TC in the studied area. In contrast, exhumation is laterally diachronous and decreases southward in the case of the HR units. In particular, our data imply that the HR units have been exhumed by a minimum of 15 km over the last few Ma. In the case of the BS, they show far less cumulated exhumation and much slower cooling rates. We propose that most of the deformation and exhumation of the Taiwan mountain belt is sustained through two underplating windows located beneath the Hsuehshan Range and the TC. Our data show significant departures from the predictions of the prevailing model in Taiwan, which assumes a homogeneous critical wedge with dominant frontal accretion. Our study sheds new light on how the mountain belt has grown as a possible result of underplating mostly.

Rupture process of the Mw = 7.9 2015 Gorkha earthquake (Nepal): Insights into Himalayan megathrust segmentation

We investigate the rupture process of the 25 April 2015 Gorkha earthquake (Mw = 7.9) using a kinematic joint inversion of teleseismic waves, strong motion data, high-rate GPS, static GPS, and synthetic aperture radar (SAR) data. The rupture is found to be simple in terms of coseismic slip and even more in terms of rupture velocity, as both inversion results and a complementing back projection analysis show that the main slip patch broke unilaterally at a steady velocity of 3.1–3.3 km/s. This feature likely contributes to the moderate peak ground acceleration (0.2 g) observed in Kathmandu. The ~15 km deep rupture occurs along the base of the coupled portion of the Main Himalayan Thrust and does not break the area ranging from Kathmandu to the front. The limitation in length and width of the rupture cannot be identified in the preearthquake interseismic coupling distribution and is therefore discussed in light of the structural architecture of the megathrust.

Mountain building in Taiwan: A thermokinematic model

The Taiwan mountain belt is classically viewed as a case example of a critical wedge growing essentially by frontal accretion and therefore submitted to distributed shortening. However, a number of observations call for a significant contribution of underplating to the growth of the orogenic wedge. We propose here a new thermokinematic model of the Taiwan mountain belt reconciling existing kinematic, thermometric and thermochronological constraints. In this model, shortening across the orogen is absorbed by slip on the most frontal faults of the foothills. Crustal thickening and exhumation are sustained by underplating beneath the easternmost portion of the wedge (Tananao Complex, TC), where the uplift rate is estimated to ~6.3 mm a^(−1), and beneath the westernmost internal region of the orogen (Hsueshan Range units, HR), where the uplift rate is estimated to ~4.2 mm a^(−1). Our model suggests that the TC units experienced a synchronous evolution along strike despite the southward propagation of the collision. It also indicates that they have reached a steady state in terms of cooling ages but not in terms of peak metamorphic temperatures. Exhumation of the HR units increases northward but has not yet reached an exhumational steady state. Presently, frontal accretion accounts for less than ~10% of the incoming flux of material into the orogen, although there is indication that it was contributing substantially more (~80%) before 4 Ma. The incoming flux of material accreted beneath the TC significantly increased 1.5 Ma ago. Our results also suggest that the flux of material accreted to the orogen corresponds to the top ~7 km of the upper crust of the underthrust Chinese margin. This indicates that a significant amount (~76%) of the underthrust material has been subducted into the mantle, probably because of the increase in density associated with metamorphism. We also show that the density distribution resulting from metamorphism within the orogenic wedge explains well the topography and the gravity field. By combining available geological data on the thermal and kinematic evolution of the wedge, our study sheds new light onto mountain building processes in Taiwan and allows for reappraising the initial structural architecture of the passive margin.

The Sumatra subduction zone: A case for a locked fault zone extending into the mantle

[1] A current view is that the portion of the subduction interface that remains locked in the time interval between large interplate earthquakes, hereinafter referred to as the locked fault zone (LFZ), does not extend into the mantle because serpentinization of the mantle wedge would favor stable aseismic sliding. Here, we test this view in the case of the Sumatra subduction zone where the downdip end of the LFZ can be well constrained from the pattern and rate of uplift deduced from coral growth and from GPS measurements of horizontal deformation. These geodetic data are modeled from a creeping dislocation embedded in an elastic half-space and indicate that the LFZ extends 132 ± 10/7 km from the trench, to a depth between 35 and 57 km. By combining this information with the geometry of the plate interface as constrained from two-dimensional gravimetric modeling and seismicity, we show that the LFZ extends below the forearc Moho, which is estimated to lie at a depth of � 30 km, at a horizontal distance of 110 km from the trench. So, in this particular island arc setting, the LFZ most probably extends into the mantle, implying that either the mantle is not serpentinized, or that the presence of serpentine does not necessarily imply stable sliding. From thermal modeling, the temperature at the downdip end of the LFZ is estimated to be 260 ± 100� C. This temperature seems too low for thermally activated ductile flow, so that aseismic slip is most probably due to pressure and/ or temperature induced steady state brittle sliding, possibly favored by fluids released from the subducting slab. INDEX TERMS: 7223 Seismology: Seismic hazard assessment and prediction; 1206 Geodesy and Gravity: Crustal movements—interplate (8155); 1219 Geodesy and Gravity: Local gravity anomalies and crustal structure; 3902 Mineral Physics: Creep and deformation; KEYWORDS: locked fault zone, interseismic deformation, Sumatra

Investigating the kinematics of mountain building in Taiwan from the spatiotemporal evolution of the foreland basin and western foothills

The Taiwanese range has resulted from the collision between the Luzon volcanic arc and the Chinese continental margin, which started about 6.5 Myr ago in the north, and has since propagated southward. The building of the range has been recorded in the spatiotemporal evolution of the foreland basin. We analyze this sedimentary record to place some constraints on the kinematics of crustal deformation. The flexure of the foreland under the load of the growing wedge started with a 1.5 Myr long phase of rapid subsidence and sedimentation, which has migrated southward over the last 3.5 Myr at a rate of 31 +10/−5 mm/yr, reflecting the structural evolution of the range and the growth of the topography during the oblique collision. Isopachs from the Toukoshan (~0 to 1.1 Ma) and Cholan (~1.1 to 3.3 Ma) formations, as well as the sedimentation rates retrieved from a well on the Pakuashan anticline, indicate that the foreland basement has been moving toward the center of mass of the orogen by ~45–50 mm/yr during the development of the basin. From there, we estimate the long-term shortening rate across the range to 39.5–44.5 mm/yr. By considering available data on the thrust faults of the foothills of central Taiwan, we show that most (if not all) the shortening across the range is accommodated by the most frontal structures, with little if any internal shortening within the wedge. The range growth appears therefore to have been essentially sustained by underplating rather than by frontal accretion. In addition, only the upper ~7 to 9 km of the underthrusted crust participates to the growth of the orogen. This requires that a significant amount of the Chinese passive margin crust is subducted beneath the Philippine Sea plate.

Kinematics of fault-related folding derived from a sandbox experiment

We analyze the kinematics of fault tip folding at the front of a fold-and-thrust wedge using a sandbox experiment. The analog model consists of sand layers intercalated with low-friction glass bead layers, deposited in a glass-sided experimental device and with a total thickness h = 4.8 cm. A computerized mobile backstop induces progressive horizontal shortening of the sand layers and therefore thrust fault propagation. Active deformation at the tip of the forward propagating basal decollement is monitored along the cross section with a high-resolution CCD camera, and the displacement field between pairs of images is measured from the optical flow technique. In the early stage, when cumulative shortening is less than about h/10, slip along the decollement tapers gradually to zero and the displacement gradient is absorbed by distributed deformation of the overlying medium. In this stage of detachment tip folding, horizontal displacements decrease linearly with distance toward the foreland. Vertical displacements reflect a nearly symmetrical mode of folding, with displacements varying linearly between relatively well defined axial surfaces. When the cumulative slip on the decollement exceeds about h/10, deformation tends to localize on a few discrete shear bands at the front of the system, until shortening exceeds h/8 and deformation gets fully localized on a single emergent frontal ramp. The fault geometry subsequently evolves to a sigmoid shape and the hanging wall deforms by simple shear as it overthrusts the flat ramp system. As long as strain localization is not fully established, the sand layers experience a combination of thickening and horizontal shortening, which induces gradual limb rotation. The observed kinematics can be reduced to simple analytical expressions that can be used to restore fault tip folds, relate finite deformation to incremental folding, and derive shortening rates from deformed geomorphic markers or growth strata.

Kinematic analysis of the Pakuashan fault tip fold, west central Taiwan: Shortening rate and age of folding inception

The Pakuashan anticline is an active fault-tip fold that constitutes the frontalmost zone of deformation along the western piedmont of the Taiwan Range. Assessing seismic hazards associated with this fold and its contribution to crustal shortening across central Taiwan requires some understanding of the fold structure and growth rate. To address this, we surveyed the geometry of several deformed strata and geomorphic surfaces, which recorded different cumulative amounts of shortening. These units were dated to ages 2 ranging from ~ 19 ka to ~ 340 ka using Optical Stimulated Luminescence (OSL). We collected shallow seismic profiles and used previously published seismic profiles to constrain the deep structure of the fold. These data show that the anticline has formed as a result of pure shear with subsequent limb rotation. The cumulative shortening along the direction of tectonic transport (N118E) is estimated to be 1010 +/- 160 m. An analytical fold model derived from a sandbox experiment [Bernard, et al., in press] is used to model growth strata. This yields a shortening rate of 16.3 +/- 4.1 mm/yr and constrains the time of initiation of deformation to 62.2 +/- 9.6 ka. In addition, the kinematic model of Pakuashan is used to assess how tectonics, sedimentation and erosion have sculpted the present-day fold topography and morphology. The fold model, applied here for the first time on a natural example, appears promising in determining the kinematics of fault-tip folds in similar contexts and therefore in assessing seismic hazards on blind thrust faults.

Slip rates on the Chelungpu and Chushiang thrust faults inferred from a deformed strath terrace along the Dungpuna river, west central Taiwan

The Chelungpu fault produced the September 1999 M_w = 7.6 Chi-Chi earthquake, central Taiwan. The shortening rate accommodated by this structure, integrated over several seismic cycles, and its contribution to crustal shortening across the Taiwanese range have remained unresolved. To address the issues, we focus our study on the Chelungpu and Chushiang thrust faults within the southernmost portion of the Chi-Chi rupture area. Structural measurements and available seismic profiles are used to infer the subsurface geometry of structures. The Chushiang and Chelungpu faults appear as two splay faults branching onto a common ramp that further north connects only to the Chelungpu surface trace. We survey a deformed strath terrace along the Dungpuna river, buried under a 11,540 ± 309 years old fill deposit. Given this age, the dip angles of the faults, and the vertical throw determined from the offset of the strath terrace across the surface fault traces, we estimate slip rates of 12.9 ± 4.8 and 2.9 ± 1.6 mm/yr on the Chelungpu and Chushiang faults, respectively. These yield a total shortening rate of 15.8 ± 5.1 mm/yr to be absorbed on their common decollement at depth. This total value is an upper bound for the slip rate on the Chelungpu fault further north, where the Chushiang fault disappears and transfers shortening to adjacent faults. Combining these results with the recently constrained shortening rate on the Changhua blind thrust reveals that all these frontal faults presently absorb most of the long-term horizontal shortening across the Taiwanese range. They thus stand as the major sources of seismic hazards in this heavily populated area. The return period of earthquakes similar to the Chi-Chi event over a ∼80 km long stretch of the Western Foothills is estimated to be ~64 years. This value is an underestimate because it assumes that all the faults locked during the interseismic period slip only during such large events. Comparison with historical seismicity suggests that episodic aseismic deformation might also play a major role in accommodating shortening.

Erosion influences the seismicity of active thrust faults

Assessing seismic hazards remains one of the most challenging scientific issues in Earth sciences. Deep tectonic processes are classically considered as the only persistent mechanism driving the stress loading of active faults over a seismic cycle. Here we show via a mechanical model that erosion also significantly influences the stress loading of thrust faults at the timescale of a seismic cycle. Indeed, erosion rates of about ~0.1-20 mm yr(-1), as documented in Taiwan and in other active compressional orogens, can raise the Coulomb stress by ~0.1-10 bar on the nearby thrust faults over the inter-seismic phase. Mass transfers induced by surface processes in general, during continuous or short-lived and intense events, represent a prominent mechanism for inter-seismic stress loading of faults near the surface. Such stresses are probably sufficient to trigger shallow seismicity or promote the rupture of deep continental earthquakes up to the surface.

Kinematics of the active West Andean fold‐and‐thrust belt (central Chile): Structure and long‐term shortening rate

West-verging thrusts, synthetic with the Nazca - South America subduction interface, have been recently discovered at the western front of the Andes. At ~33°30’S, the active San Ramon fault stands as the most frontal of these west-verging structures, and represents a major earthquake threat for Santiago, capital city of Chile. Here we elaborate a detailed 3D structural map and a precise cross-section of the West Andean fold-and-thrust belt based on field observations, satellite imagery and previous structural data, together with digital topography. We then reconstruct the evolution of this frontal belt using a trishear kinematic approach. Our reconstruction implies westward propagation of deformation with a total shortening of 9-15 km accumulated over the last 25 Myr. An overall long-term shortening rate of 0.1-0.5 mm/yr is deduced. The maximum value of this shortening rate compares well with the rate that may be inferred from recent trench data across the San Ramon fault and the slip associated with the two past Mw > 7 earthquakes. This suggests that the San Ramon fault is most probably the only presently active fault of the West Andean fold-and-thrust-belt and that most - if not all - the deformation is to be released seismically.