0000000000208759
AUTHOR
Stefan M. Schmalholz
Control of 3D tectonic inheritance on fold-and-thrust belts: insights from 3D numerical models and application to the Helvetic nappe system
Abstract. Fold-and-thrust belts and associated tectonic nappes are common in orogenic regions. They exhibit a wide variety of geometries and often a considerable along-strike variation. However, the mechanics of fold-and-thrust belt formation and the control of the initial geological configuration on this formation are still incompletely understood. Here, we apply three-dimensional (3D) thermo-mechanical numerical simulations of the shortening of the upper crustal region of a passive margin to investigate the control of 3D laterally variable inherited structures on the fold-and-thrust belt evolution and associated nappe formation. We consider tectonic inheritance by applying an initial mode…
Quantifying the impact of mechanical layering and underthrusting on the dynamics of the modern India-Asia collisional system with 3-D numerical models
The impact of mechanical layering and the strength of the Indian lower crust on the dynamics of the modern India-Asia collisional system are studied using 3-D thermomechanical modeling. The model includes an Indian oceanic domain, Indian continental domain, and an Asian continental domain. Each domain consists of four layers: upper/lower crust, and upper/lower lithospheric mantle. The Tarim and Sichuan Basins are modeled as effectively rigid blocks and the Quetta-Chaman and Sagaing strike-slip faults as vertical weak zones. The geometry, densities, and viscosities are constrained by geophysical data sets (CRUST2.0, gravity, and seismology). Both static (no horizontal movement of model bound…
On backflow associated with oceanic and continental subduction
SUMMARY A popular idea is that accretion of sediment at a subduction zone commonly leads to the formation of a subduction channel, which is envisioned as a narrow zone located above a subducting plate and filled with vigorously circulating accreted sediment and exotic blocks. The circulation can be viewed as a forced convection, with downward flow in the lower part of the channel due to entrainment by the subducting plate, and a ‘backflow’ in the upper part of the channel. The backflow is often cited as an explanation for the exhumation of high-pressure/low-temperature metamorphic rocks from depths of 30 to 50 km. Previous analyses of this problem have mainly focused on the restricted case …
2D Hydro-Mechanical-Chemical Modeling of (De)hydration Reactions in Deforming Heterogeneous Rock: The Periclase-Brucite Model Reaction
Deformation at tectonic plate boundaries involves coupling between rock deformation, fluid flow, and metamorphic reactions, but quantifying this coupling is still elusive. We present a new two-dimensional hydro-mechanical-chemical numerical model and investigate the coupling between heterogeneous rock deformation and metamorphic (de)hydration reactions. We consider linear viscous compressible and power-law viscous shear deformation. Fluid flow follows Darcy's law with a Kozeny-Carman type permeability. We consider a closed isothermal system and the reversible (de)hydration reaction: periclase and water yields brucite. Fluid pressure within a circular or elliptical inclusion is initially bel…
Nonlithostatic pressure during subduction and collision and the formation of (ultra)high-pressure rocks
The mechanisms that result in the formation of high-pressure (HP) and ultrahigh-pressure (UHP) rocks are controversial. The usual interpretation assumes that pressure is close to lithostatic, petrological pressure estimates can be transferred to depth, and (U)HP rocks have been exhumed from great depth. An alternative explanation is that pressure can be larger than lithostatic, particularly in continental collision zones, and (U)HP rocks could thus have formed at shallower depths. To better understand the mechanical feasibility of these hypotheses, we performed thermomechanical numerical simulations of a typical subduction and collision scenario. If the subducting crust is laterally homogen…
Quantification and visualization of finite strain in 3D viscous numerical models of folding and overthrusting
Abstract Finite strain analysis and three-dimensional (3D) numerical modeling are important methods to understand the deformation history of rocks. Here, we analyze finite strain in 3D numerical simulations of power-law viscous folding and overthrusting. Simulations with different and laterally varying detachment geometries cause a lateral transition from folding (for thicker detachments) to overthrusting. We compute the 3D finite strain tensor, the principal strain values and their orientations. We compute the Nadai strain, e S , and the Lode’s ratio, ν , representing the strain symmetry (constriction or flattening). We apply Hsu diagrams to visualize strain distribution in e S - ν space, …
Transition from thin- to thick-skinned tectonics and consequences for nappe formation: Numerical simulations and applications to the Helvetic nappe system, Switzerland
Abstract We use two-dimensional numerical simulations to study the deformation style and the formation of nappes during the shortening of basement–cover systems with half-grabens. We perform simulations with only linear and power-law viscous rheology, and also thermo-mechanical simulations with viscoelastoplastic rheology and temperature-dependent viscosities. The accumulated, area-averaged strain is calculated for the sediments, es, and for the basement, eb, to quantify the deformation style. The results show that the transition from thin-skinned dominated (es/eb > 1) to thick-skinned dominated (es/eb ≤ 1) deformation is controlled by two effective viscosity ratios: the ratio of viscosity …
Widening of Hydrous Shear Zones During Incipient Eclogitization of Metastable Dry and Rigid Lower Crust—Holsnøy, Western Norway
The partially eclogitized crustal rocks on Holsnøy in the Bergen Arcs, Norway, indicate that eclogitization is caused by the interplay of brittle and ductile deformation promoted by fluid infiltration and fluid‐rock interaction. Eclogitization generated an interconnected network of millimeter‐to‐kilometer‐wide hydrous eclogite‐facies shear zones, which presumably caused transient weakening of the mechanically strong lower crust. To decipher the development of those networks, we combine detailed lithological and structural mapping of two key outcrops with numerical modeling. Both outcrops are largely composed of preserved granulite with minor eclogite‐facies shear zones, thus representing th…