Editorial: Cutting-Edge Analogue Modeling Techniques Applied to Study Earth Systems

Our understanding of Earth systems is built on field observations, geological and geophysical investigations and modeling. For over two hundred years, geologists are building analog models to test theories and understand the physics leading to field observations. Analog models do not aim to reproduce nature but rather to simplify the system so that parameters like geometry, kinematics, or dynamics can be isolated and investigated. Analog models allow to investigate complex three-dimensional problems at high-resolution. In addition to deciphering outcrop observations, analog models offer the opportunity to predict structures not accessible for direct observation. Analog models provide a full 4-D view of geological processes, allowing for investigating the time evolution of structures

Combining low-temperature thermochronology with 3-D probabilistic kinematic modeling including uncertainties in the Eastern Alps

To understand the exhumation history of the Alps and its foreland, it is important to accurately reconstruct its time-temperature evolution. This is often done employing thermokinematic models. However, one problem of many current approaches is that they rely on prescribed geometric structures at depth without considering their uncertainty. Therefore, the aim of this work is to compare low-temperature thermochronological data with a 3-D probabilistic kinematic model. To this end, we combine 3-D kinematic forward modeling with a systematic random sampling approach to automatically generate an ensemble of kinematic models in the range of assigned uncertainties. These can later be used to obtain a 3-D probabilistic exhumation map, from which exhumation values for the sample positions of thermochronological data can be interpolated, and compared to estimates made solely from thermochronology. In a next step, the uncertainties assigned to the kinematic model can be updated with the thermochronological data, to obtain an even more robust model. We apply this approach to the Bavarian Subalpine Molasse, which is particularly suited as a test case, as it connects the Alpine orogen with its foreland, and should shed light on the strain distributions during the latest stages of Alpine mountain building. Preliminary results using previously published data show that the estimated exhumation from the modeling can serve as a constraint to thermochronological interpretations, leading to an uncertainty reduction. In a next step, we will use our own (U-Th)/He measurements to obtain an integrated picture of foreland evolution and associated uncertainties over space and time

Dilatant normal faulting in jointed cohesive rocks: a physical model study

Dilatant faults often form in rocks containing pre-existing joints, but the effects of joints on fault segment linkage and fracture connectivity are not well understood. We present an analogue modeling study using cohesive powder with pre-formed joint sets in the upper layer, varying the angle between joints and a rigid basement fault. We analyze interpreted map-view photographs at maximum displacement for damage zone width, number of connected joints, number of secondary fractures, degree of segmentation and area fraction of massively dilatant fractures. Particle imaging velocimetry provides insight into the deformation history of the experiments and illustrates the localization pattern of fault segments. Results show that with increasing angle between joint-set and basement-fault strike the number of secondary fractures and the number of connected joints increase, while the area fraction of massively dilatant fractures shows only a minor increase. Models without pre-existing joints show far lower area fractions of massively dilatant fractures while forming distinctly more secondary fractures

Dating the youngest deformation in the Alps with ESR thermochronometry

Low-temperature thermochronology is a useful tool to reconstruct tectonic deformation and landscape evolution within the first 2 km of the crust. It is a suitable tool to investigate deformation associated with cooling and exhumation of the lower crust in orogenic settings. Low temperature thermochronology is applied here to understand the Neogenic post-collisional extensional event that occurred in the Alps, because a gap in previous age dating exists between a thousand and a million years. Quartz is the most common mineral in the crust; occurring in magmatic as well as sedimentary and metamorphic rocks. The potential of quartz electron-spin resonance (ESR) as a radiation dosimeter has been well documented, and many studies applied the method to date sediments and heated rocks (e.g. tephra). In this study, we apply quartz ESR dating as an ultralow-temperature thermochronometer, characterized by a closure temperature of 30°-90°, and dating range of 103-107 years. We show the results of ESR thermochronometry on quartz applied to rocks from crustal-scale faults in the Central (Simplon Fault) and Eastern Alps (Brenner and Salzachtal Faults). Here, the lower crust has been tectonically exhumed, associated with exhumation of the Lepontine Dome and Tauern Window, respectively. Thermochronological data are available from this area, such as fission tracks or U-Th/He data on zircon and apatite. Results of the ESR measurements of 15 samples crossing the Brenner and Salzachtal faults (northern and western border of the Tauern Window) show that the ESR ages of quartz get younger (<1Ma) inside the western part of the Tauern Window, in accordance with fission track and (U-Th)/He ages. In general, younger ages (between 200 and 500 ka) are also obtain closer to the fault zone, localized near (e.g. Simplon Fault) or at the bottom of the valley (e.g. Brenner Fault), compared with the protolithic rocks (600-900 ka). We interpret the trend of the ESR ages as an exhumation of the isotherms due to both recent uplift of the footwall of the fault and for erosion of the valley, where the later overprints the former. These results promise to establish ESR as an ultra-low thermochronometer using quartz for the Quaternary landscape reconstruction of the Alpine chain

The Memory of a Fault Gouge: An Example from the Simplon Fault Zone (Central Alps)

Faut gouge forms at the core of the fault as the result of a slip in the upper brittle crust. Therefore, the deformation mechanisms and conditions under which the fault gouge was formed can document the stages of fault movement in the crust. We carried out a microstructural analysis on a fault gouge from a hanging-wall branch fault of the Simplon Fault Zone, a major low-angle normal fault in the European Alps. We use thin-section analysis, together with backscattered electron imaging and X-ray diffractometry (XRD), to show that a multistage history from ductile to brittle deformation within the fault gouge. We argue that this multistage deformation history is the result of continuous exhumation history from high to low temperature, along the Simplon Fault Zone. Because of the predominance of pressure solution and veining, we associated a large part of the deformation in the fault gouge with viscous-frictional behaviour that occurred at the brittle-ductile transition. Phyllosilicates and graphite likely caused fault lubrication that we suggested played a role in localizing slip along this major low-angle normal fault

The Memory of a Fault Gouge: An Example from the Simplon Fault Zone (Central Alps)

A fault gouge forms at the core of the fault as the result of a slip in the upper brittle crust. Therefore, the deformation mechanisms and conditions under which the fault gouge was formed can document the stages of fault movement in the crust. We carried out a microstructural analysis on a fault gouge from a hanging-wall branch fault of the Simplon Fault Zone, a major low-angle normal fault in the Alps. We use thin-section analysis, together with backscattered electron imaging and X-ray diffractometry (XRD), to show that a multistage history from ductile to brittle deformation, together with a continuous exhumation history from high to low temperature, took place within the fault gouge. Because of the predominance of pressure solution and veining, we associated a large part of the deformation in the fault gouge with viscous-frictional behaviour that occurred at the brittle-ductile transition. Phyllosilicates and graphite likely caused fault lubrication that we suggested played a role in the formation of this major low-angle normal fault. © 2022 by the authors. Licensee MDPI, Basel, Switzerland

Critical taper analysis reveals lithological control of variations in detachment strength: An analysis of the Alpine basal detachment (Swiss Alps)

Although evidence for weak detachments underlying foreland thrust belts exists, very little is known about the lateral variations in effective strength, as well as the geological nature of such variations. Using critical taper analysis, we show that a detailed and systematic measurement of surface slope of the Central European Alps reveals variations in strength parameter F along the detachment, based on the argument that the Alps are close to the critical state. We show that the basal detachment is very weak near the deformation front but strengthens toward the hinterland. Very low F (effective coefficient of friction plus normalized cohesion) values of  0.54) but may also require additional mechanisms of dynamic weakening

Foreland dynamics as a measure of mountain building processes

Forelands record the uplift and exhumation history of mountain belts. The alpine foreland basin is particularly exciting, as is shows late-orogenic exhumation, possibly as a reaction to mantle-driven, plate convergence, or climatic forcings. However, inferring the contribution of the individual drivers to exhumation from stratigraphic or thermochronological data is challenging. The reason for this are along strike variability basin of stratigraphy, different degree of exhumation, as well as structural style of the Subalpine Molasse (i.e., the fold-thrust belt at the southern fringe of the basin). Furthermore, the influence of fluid flow on the thermochronological ages is unknown. Exhumation estimates in the central part of the basin are mostly based on stratigraphic arguments. Thermochronological data is scarce and limited to local studies. As the Molasse has also been uplifted in the central part of the basin since the Miocene, it is probable that it also responds to deep-seated processes, but to a lesser extent than the western part of the basin. This may be a result of different slab dynamics along strike the orogen. To test this, we used detrital and in situ low-temperature thermochronological age dating to shed light on the surface expression of the underlying geodynamic process (Figure 1). Data shows that most ages in the central part of the basin are unreset, while resetting occurs in the southernmost tectonic slices of the Subalpine Molasse. Generally, Miocene shortening in the Subalpine Molasse progressively decreases from west to east. The pattern coincides with slab geometries at depth (Mock et al., 2020). A general trend of lesser erosion from west to east is also visible in the flat lying Molasse based on vitrinite reflectance data. This suggests that a geodynamic driver is required for explaining basin exhumation on basin scale. Locally, the pattern is more complex. Particularly in the Subalpine Molasse, exhumation may be associated with plate convergence. To test the influence of faulting on exhumation, we constrained the geometries of the fold-thrust belt. Using a new compilation of stratigraphy and structures along the entire Alpine deformation front (Ortner et al., 2023), we identified two key regions: the Bregenzerach south of the eastward termination of the Jura Mountains, and the Hausham Syncline southeast of Munich. The Bregenzerach region lies at the surface boundary between Eastern and Western Alps. Furthermore, previously published thermochronological data indicate thrust activity in the mid-Miocene. Structures at depth are reasonably well-constrained due to good outcrop conditions and seismic data. The Hausham Syncline represents the region where structures at depth are less well constrained, and additionally the frontal triangle zone of the Subalpine Molasse tapers out. Structural modeling shows that it is possible to quantify the uncertainty of structures at depth, paving towards thermo-kinematic modeling including structural uncertainty (Brisson et al., 2023; Frings et al., 2023). The extensive thermochronological dataset offers the opportunity to identify local particularities not in line with the general trends observed in the data. Using thermal springs as proxy for heat flow (Luijendijk et al., 2020), we show that fluid flow may at least locally influence the cooling pattern. This is important for translating cooling into exhumation, particularly in regions where less data is available and thus outliers may be overlooked

Time scale bias in erosion rates of glaciated landscapes

Deciphering erosion rates over geologic time is fundamental for understanding the interplay between climate, tectonic, and erosional processes. Existing techniques integrate erosion over different time scales, and direct comparison of such rates is routinely done in earth science. On the basis of a global compilation, we show that erosion rate estimates in glaciated landscapes may be affected by a systematic averaging bias that produces higher estimated erosion rates toward the present, which do not reflect straightforward changes in erosion rates through time. This trend can result from a heavy-tailed distribution of erosional hiatuses (that is, time periods where no or relatively slow erosion occurs). We argue that such a distribution can result from the intermittency of erosional processes in glaciated landscapes that are tightly coupled to climate variability from decadal to millennial time scales. In contrast, we find no evidence for a time scale bias in spatially averaged erosion rates of landscapes dominated by river incision. We discuss the implications of our findings in the context of the proposed coupling between climate and tectonics, and interpreting erosion rate estimates with different averaging time scales through geologic time