Search

Inversion and exhumation of crystalline basement units of passive margins is a critical stage of mountain-building processes. With a multi-methodological approach that combines 3D geological modelling and cross-section restoration at different time steps, we unravel the 4D geodynamic evolution of the Aar Massif (central Alps, Switzerland) during the late-stage Alpine orogeny (22 to 0 Ma). Our results demonstrate how a portion of the crystalline basement units of the European continental margin was exhumed in a non-cylindrical framework with an alternation of vertical and horizontal-dominated tectonics. Rapid exhumation (22 to 12 Ma), with deformation along steep reverse/normal faults, alternated with northwest-directed thrusting of the previously exhumed basement units (16 to 0 Ma). We pay special attention to how shear/fault zone patterns change laterally and temporarily to accommodate the non-cylindrical deformation in dome-like basement rises. Published thermochronometric and metamorphic peak temperature data provide additional insights into the temporal sequence of reverse and thrust faulting perpendicular and parallel to the strike of the orogen. Such variations in the fault geometries and associated offsets were constrained by the detailed reconstruction of the shear zone pattern in 3D. In this study, we show how the basement-involved uplift and shortening controlled the late-stage collisional overprint and mechanics of the northern rim of the central Alps.

Fault gouges of granitoid composition represent the principal non-cohesive tectonites within fault zones in the continental crust. The spatial distribution and strength of granitoid fault gouges is therefore crucial for understanding how weak the upper continental crust could be due to the formation of fault zones. Although several laboratory investigations reported the mechanical weakening of granitoid gouges in shear experiments, the deformation mechanism responsible for such behavior remains not well understood.To address this issue, we conducted two series of shear experiments on granitoid gouges by using a ring shear apparatus. The starting gouge powders were derived from crushed granitoid mylonite with a median grain size of 45 μm. In a first set of experiments, gouges were sheared at a sliding velocity of 100 μm/s for a displacement of 15 mm. Temperatures explored ranged from 20°C to 650°C in order to determine the temperature dependence of gouge strength. The second set of experiments is identical to the first ones, except that the applied sliding velocity was set at 1 μm/s to study how fault slip rate influences the strength of gouges.We observe that differences in gouges strengths as a function of sliding velocity and temperature: At a sliding velocity of 100 μm/s, the steady-state shear stress (τ) remains relatively constant at τ=76-82 MPa over the entire temperature range. Contrastingly, at a sliding velocity of 1 μm/s the steady-state shear stress remains temperature-insensitive with τ≈75 MPa up to tempertures of 450°C, but decreases then to τ≈50 MPa at 650°C (Fig.1 a). Furthermore, the amount of decrease of shear stress is strain dependent (Fig.1 b). At even slower sliding velocity of 0.1 µm/s, the shear stresses decrease further to τ≈38 MPa.Microstructurally, all gouges deformed at T≦450°C show typical cataclastic features, where angular clasts with grain size of ~10 μm are surrounded by a fine-grained matrix. Intergranular fracture arrays in Riedel- and Y-shears are well developed over the entire cross section, indicating homogeneous bulk deformation. In contrast, gouges sheared at 650°C with τ≈50 MPa show strain localization in a principal slip zone. It is shear plane parallel with widths up to ~50 µm. Inside the principal slip zone, all grains are dramaticly reduced to nm-size and tightly packed. No intergranular fracture arrays are observed. Outside the principal slip zone, rounded grains with size of ~5 μm are loosely packed, with meniscus cement growing in between. The aforementioned strain localization is enhanced at temperature above 450°C and slip rate below 1μm/s, suggesting that viscous creep mechanisms (e.g. pressure solution) control the deformation process at slow sliding velocities, which is not the case in fast rate experiments. Our results show that the activation of viscous creep mechanisms leads to significant fault zone weakening, while contrasts in grain size keep deformation localized.Figure 1 Shear stress plotted as a function of temperature. Shear stress data collected at (a) 15mm displacement in steady-state, and at (b) 5mm displacement.

Crustal rocks involved in orogenic processes frequently bear evidence for widespread fluid circulation. The hydration history of the granitic continental crust in inverted passive continental margins is of particular interest, as granitoids experience rheological weakening by fluid-rock interaction processes. Regrettably, it is often unclear if hydration occurs during rifting or during tectonic inversion. Hence, it is difficult to appreciate the interplay of pre- and syn-orogenic fluids inside continental crustal segments of rifted margins. The geochemical fingerprint of ancient hydration events is stored in hydrous minerals that crystallized directly from circulating paleo-fluids. Thus, such minerals can shed light on the nature of these ancient fluids, as well as provide temporal constraints if they can be dated. Hence, advances in geochronological methods applied to hydrous minerals may prove pivotal in untangling the history of fluid circulation in the granitic continental crust in orogens.We applied U–Pb geochronology of epidote [i.e., Ca2Al2(Al,Fe3+)Si3O12(OH)] in hydrothermal veins hosted by a late Carboniferous/early Permian calc-alkaline granodiorite in the inverted Adriatic passive continental margin (hereafter “Err nappe”), both located in the eastern Swiss Alps. During Jurassic rifting leading to the break-up of Pangea, the continental crust in the Err nappe was hydrated, as seawater-derived fluids percolated along syn-rift faults. However, geochronological data of epidote reveal that the hydration of the granitic continental crust in the Err nappe occurred also later during inversion. Epidote U–Pb geochronology returned two age clusters: (1) 85.2 ± 9.7 Ma, related to Late Cretaceous compression; and (2) 59.9 ± 2.7 Ma, related to subsequent Paleocene extension. These age clusters unveil two distict events of fluid circulation, which are consistent with the timing of tectonic inversion and deformation proposed in the literature. As confirmed by Pb–Sr–O–H isotope geochemistry of epidote, Late Cretaceous fluid circulation was likely mediated by fluids released by underlying units undergoing metamorphism during Eo-Alpine compression. Notably, the Paleocene fluids circulating during extension were most likely surficial in origin (i.e., meteoric water and/or modified/connate seawater), and they percolated into the granitic continental crust by exploiting extensional faults.In the context of existing data, our results show that the hydration of the granitic continental crust of the Adriatic passive continental margin was mediated by a repeated series of fluid circulation events. Our work advocates that the use of a multi-methodological approach, combining new geochemical and geochronological, tools provides unprecedented insight into complex processes of fluid circulation in the continental crust, and beyond.

On October 05, 2021 an Mw4.0 earthquake struck 6 km south of the village of Arolla, near the tongue of the Arolla Glacier. Almost one year prior to this earthquake, an M3.5 event occurred on November 08, 2020 in the same location. Both earthquakes were followed by a few aftershocks that were detected and located by the Swiss Seismological Service (SED). The unusually shallow depth of 1-2 km of these earthquakes, indications for a mostly thrust-type mechanisms within a region characterized by a predominantly extensional stress regime, and unusual high CLVD (50-70%) components of SED’s routine moment tensor solutions raised questions regarding the triggering mechanism. To understand and explain the possible existence of shallow thrust earthquakes in the area, we perform a thorough seismotectonic analysis that is based on enhancing the existing earthquake catalog of the SED and complementary moment-tensor solutions computed by multiple algorithms. The original SED earthquake catalog contains 83 earthquakes that occur between January 01, 2020 and December 31, 2021 and locate ~5 km around the two mainshocks. Using a deep learning based algorithm (EQTransformer), we detect additionally 253 events, thus the new catalog contains 4 times more earthquakes than the original SED bulletin. Absolute locations for the additional earthquakes are obtained using the probabilistic NonLinLoc method in combination with a recently updated Vp and Vs crustal 3D velocity model. In addition, we compute local magnitudes (MLhc) using SED’s standard procedure, in order to compile a homogeneous catalog consistent with the SED bulletin. The enhanced catalog events are used as templates for a match filtering scheme, which increased the number of detections by at least one magnitude order. Last, we relocate the final catalog using the double difference method towards obtaining a high resolution enhanced earthquake catalog. Spatially, the main cluster shows an intense seismic activity, stretched in N-S direction that matches the strike of the fault planes derived from moment tensor inversion. An additional cluster, that is not present in the SED bulletin locations, is identified next to the area were the aftershock activity of the two main events locates. Furthermore, the enhanced catalog shows a smother temporal evolution with more background events than previously recorded. Overall, we explore the possibility of fluid driven microseismicity that might be related to the nearby glacier. With our study we emphasize the importance of enhanced earthquake catalogs using both machine learning pickers and template matching algorithms. These approaches lead to unravel prior unmapped structures and improve our understanding of the seismotectonic regime in the study area.

Understanding orogen-internal seismic deformation in regions with diffuse spatial earthquake occurrence is challenging. To gain deeper insights into the processes driving seismic fault reactivation, it is crucial to obtain information on the ubiquitous pre-existing fracture patterns. In orogens with long tectonic histories – such as the Alps – such patterns can be complex, and information on their appearance is mainly limited to observations of faults at the surface, while the detailed patterns at depth remain mostly unknown. Moreover, the link between such surface-based fault observations and active seismicity is often ambiguous. However, it has been shown that both earthquake magnitudes (Gutenberg-Richter law) and various fault properties (e.g., length, displacement) follow power-law distributions.In this work, we aim to investigate the potential relationship between the scaling properties of faults and earthquakes, which has been little explored. To this end, we use statistical tools based on field data collected with remote sensing techniques at different scales to quantitatively characterize the length distributions of exposed fault networks at different study sites in the southwestern Swiss Alps. Due to the good outcrop conditions at high elevations, the dense seismic monitoring network, and the enhanced earthquake activity, this region provides an ideal natural laboratory for the study of orogen-internal seismicity. By combining fault trace maps from three different scales, we are able to derive power law parameters and decipher similarities in scaling exponents for the different sites studied. Assuming that the fault networks exist in a similar form at depth and form the pre-existing discontinuities along which recent earthquakes develop, we compare the derived scaling laws with the frequency-magnitude distribution of local seismicity over the past 15 years. Here we find similar scaling properties between the seismicity and fracture networks only at depths below 3 km. However, in shallower regions, the large discrepancy between the scaling laws suggests that partial seismic ruptures of individual fault segments are more common than at greater depths. Such a statistical comparison of fault and earthquake scaling laws provides interesting insights into orogen-internal seismic deformation and fault reactivation processes that also have implications for regional seismic hazard.

According to well-established hypotheses based on field observations of natural faults, viscous deformation may localize following pre-existing brittle fractures. The weak behaviour can be explained by brittle grain size reduction and phase mixing, which may activate grain size sensitive processes in the viscous field. To prove this hypothesis, it is necessary to perform experiments to observe the strain and stress evolution in faulted and non-faulted rocks. Pec et al. (2012) performed experiments on granitic rocks by shearing manually crushed granitic powder between coarse solid granitic forcing blocks. However, in their study, there are unavoidable boundary conditions between the forcing blocks and the gouge, and a comparison to an intact rock without fracture is difficult.In our study, we reduce the boundary conditions to a minimum and can directly compare the stresses and microstructural evolution during deformation of intact and fractured granitic ultramylonites at 650°C, confining pressure of 1.2GPa, and a constant displacement rate of 10-8m/s. We perform these experiments on initially solid cylindrical samples in two experimental sets: In set A, we slowly apply the load and confining pressure, to ensure an intact rock sample is deformed. In set B, we create fractures before the experiment starts but already in the closed system of the experimental setup. Once experimental P/T conditions are reached, both experimental sets are deformed to different finite strains to investigate the associated microstructural evolution. The deformation is disseminated in the set A experiments, but localizes strongly along the fracture in experimental set B. The strain is accommodated by viscous granular flow incorporating an impressive grain size reduction of up to 1000x and dissolution/precipitation processes. In addition, the stress records show that in experiments A, initially a 30% higher yield stress has to be overcome before steady state flow, while in set B steady state flow is reached directly without a strain softening increment. In both sets, steady state stresses range around 300MPa, i.e. far below the confining pressure.Applying microstructural observations and mechanical data of our experiments to deformation of granitoid crust in nature reveals that fractures serve to reach mechanical steady state earlier compared to non-fractured crust. As a matter of strain, however, both settings may yield at the same mechanical strengths of resulting shear zones. It is important to note that polymineralic fine-grained ultramylonites are up to four times weaker than monomineralic quartz, presenting an important behaviour of efficient strain localization and rheological properties substantially below those of the end member minerals. Pec, M., Stünitz, H. and Heilbronner, R., 2012. Semi-brittle deformation of granitoid gouges in shear experiments at elevated pressures and temperatures. Journal of Structural Geology, vol. 38, pp. 200-221. https://doi.org/10.1016/j.jsg.2011.09.

Deformation of polymineralic aggregates can be accommodated by viscous granular flow, a process mediated by the interplay among intracrystalline plasticity and dissolution–precipitation, each active in specific minerals under given P–T conditions. Some rock-forming minerals like quartz and feldspars have been intensively studied in terms of deformation processes. Instead, the deformation behavior of epidote and its role during viscous granular flow is not well investigated, although this mineral is ubiquitous in granitic rocks deforming under greenschist-facies conditions. In this contribution, we provide microstructural and geochemical evidence for the occurrence of dissolution–precipitation of epidote during deformation of an epidote–quartz vein. The main part of the vein is deformed, producing a fold, which is visible due to relicts of primary-growth layering inside the vein. The deformation mechanisms active during deformation include dynamic recrystallization of quartz by subgrain rotation recrystallization, producing grain size reduction in the primary vein quartz. Recrystallization occurs contemporaneously with dissolution and (re)precipitation of epidote and quartz grain boundary sliding, leading to a combined process described as viscous granular flow. The combination of grain boundary sliding and dissolution locally and repeatedly produces creep cavities. These represent not only loci for nucleation of new epidote grains at the expense of dissolved ones, but they also allow fluid-mediated transport of elements. The same trace element patterns between old epidote relicts and newly formed grains, with much narrower variability in the latter, indicate a process of chemical homogenization. The nature of the fluid that mediates deformation is investigated using Pb–Sr isotope data of epidote, which suggest that deformation is assisted by internally recycled fluids with the addition of a syn-kinematic external fluid component.

Despite the generally accepted concept that most earthquakes occur along pre-existing faults, the complex 3D geometries of seismically active fault systems at depth often remain unresolved. However, earthquake nucleation and migration processes are heavily influenced by the geometries and properties of such pre-existing structures, which limits our general understanding of earthquake nucleation and fault interactions.Under the assumption that faults are reactivated at spatially and temporally different localities, previous studies have attempted to derive fault geometries from hypocenter locations, but were usually limited by the precision of relocation techniques. Enabled by the recent advances in hypocenter relocation techniques, we present a novel Monte Carlo-based method that uses relatively relocated hypocenters and their uncertainties to image geometries, stress states and kinematics of seismically active fault systems. The application of the developed Python toolbox on a natural earthquake sequence along the Rhone-Simplon fault zone in the northern Valais (Swiss Alps) reveals active strike-slip faults with a contractional stepover. Performed stress analyses indicate varying stress states along the fault system, which has direct implications for fault properties such as the reactivation potential or the fluid transmissivity. Overall, we document the migration of an earthquake swarm across a complex strike-slip fault system at an unprecedented spatiotemporal resolution.Our toolbox can be applied to high-precision hypocenter catalogs of natural earthquake sequences or hydraulic stimulation experiments, which could help to improve our understanding of the role of pre-existing faults on earthquake nucleation and migration processes at various scales.

The susceptibility of catchment rocks to glacial erosion may control the evolution of valley morphology in high-relief mountain ranges such as the Alps. Non-uniform proneness to bedrock erosion may indeed localize knickpoints and overdeepenings characteristic of glacial valleys. Yet, little is known about the explicit influence of bedrock properties (i.e. lithology, hardness, and geological structures) on glacial erosion processes. In this study, we select the Great Aletsch Glacier (Swiss Alps) as a natural laboratory to document and investigate the relationship between bedrock properties and subglacial erosion mechanisms. The Great Aletsch Glacier with a length of more than 20km and an ice thickness of up to 800m is the largest glacier in Central Europe. The underlying bedrock consists of the crystalline basement units of the Aar massif (gneiss, granite, and granodiorite) and is dissected by a large number of steep faults and former ductile shear zones. Geological and remote sensing lineament mapping combined with 3D geological modelling allowed us to make a large-scale characterization of the lithologies and structures’ spatial frequency over the entire length of the glacier. Additionally, we performed field-based rock hardness analyses (Schmidt hammer) along the glacier’s bedrocks (intact rock and faulted/sheared domains) to testify for structure-controlled erosion behaviour. Obtained results demonstrate that: (i) the typology and distribution of faults and shear zones are not uniform over the entire length of the glacier; (ii) high-frequency structure domains correlate with overdeepenings and/or abrupt glacier flow deflection in the direction of the strike of the structures; (iii) low-frequency structure domains correlate to the absence of overdeepenings and a straight glacier trajectory. In terms of erosive resistance domains of intact rock masses show high hardness values for each of the investigated lithologies without substantial variability between the different basement rocks (rebound values ranging from 45 to 60 N/mm2). On the contrary, faulted or sheared domains show a significant drop in hardness value (rebound values ranging from 10 to 40 N/mm2). Based on these results we propose that, for the case of the Great Aletsch Glacier, differences in crystalline basement lithologies do not exert an important role in glacial erosion. We postulate instead that the non-uniform spatial distribution of geological structures imposes a major control on the development of the glacial valley. The substantially reduced bulk hardness within high-frequency structure domains renders indeed the bedrock to be more prone to efficient glacial erosion process at these sites (i.e. glacial quarrying) and therefore to the development of large-scale overdeepenings, local scouring, or changes in the glacier flow direction. By contrast, the more massive undeformed and therefore less erosive low-frequency structures domains coincide with sections with no knickpoints or overdeepenings. In times of global warming and glacial retreat, such structure-controlled bedrock incisions are prone for further enhanced surface weathering and gravitation-controlled erosion processes, such as rockfalls and landslides, providing sites of enhances natural hazard potential.

Fault gouges of granitoid composition represent the principal non-cohesive tectonites within fault zones in the continental crust. Their velocity-dependent friction is crucial for understanding earthquake nucleation and the depth distribution of fault-related seismicity in granitoid shear zones (Wehrens et al. 2016; Blanpied et al. 1998). In the framework of rate-and-state friction laws (RSF), the friction parameter (a-b) is measured in sliding experiments to describe the velocity dependence of friction in fault gouges (Scholz, 1998). A velocity-strengthening system is frictionally stable, (a-b) >0, whereas a velocity-weakening system can be frictionally unstable, (a-b) To address this issue, we conducted velocity stepping sliding experiments under hydrothermal conditions by using a ring shear apparatus. The powdered starting gouge was derived from a granitoid mylonite collected at the NAGRA Grimsel Test Site (Central Swiss Alps). The applied velocity steps were 1-3-10-30-100 μm/s. Pore fluid pressure and the effective normal stress were 100 MPa. Temperatures explored ranged from 20-650 °C. Values of (a-b) were obtained from RSF model inversions of the evolution of friction coefficients at mechanical steady state conditions. Our experiments showed pronounced changes in (a-b) values with across the full range of temperatures up to 650 °C and velocities investigated. At temperatures below ~100 °C and above ~400 °C, we observed mostly velocity strengthening with positive (a-b). In contrast, velocity weakening with negative (a-b) was observed between ~100 °C and ~400 °C. Samples deformed at a sliding velocity of 100 μm/s deviated slightly from this trend, as (a-b) values were negative between ~200 °C and ~400 °C.The presented experimental study demonstrates a significant influence of temperature and sliding velocity on velocity-dependence during deformation of granitoid gouge. We suggest that the observed transitions in velocity dependence reflect an interplay of interactions. In terms of crustal faulting, our data suggest the existence of a seismogenic window that limits the depth distribution of earthquakes on faults in granitoid shear. REFERENCESWehrens, P. C., Berger, A., Peters, M., Spillmann, T., Herwegh, M. 2016: Deformation at the frictional-viscous transition: Evidence for cycles of fluid-assisted embrittlement and ductile deformation in the granitoid crust, Tectonophysics, 693, 66-84.Blanpied M. L., Tullis T. E., Weeks J. D. 1998: Effects of slip, slip rate, and shear heating on the friction of granite.Scholz, C. H. 1998: Earthquakes and friction laws, Nature, 391, 37-42.

The Aar Massif is a mid-crustal basement section of the European plate and it was intensely deformed during the Alpine orogeny. Alpine deformation of Aar Massif granitoids is expressed by a pervasive network of ductile shear zones consisting of fine-grained polymineralic (ultra)mylonites dominated by viscous granular flow processes (Wehrens et al., 2016). Fluid circulation and hydration reactions are recorded by both granitic protoliths and shear zones. In particular, they are made evident by the alteration of feldspar into hydrous minerals like epidote and mica. The timing of this hydration event and, consequently, whether Alpine deformation was initiated in already altered granitoids or in fresh ones were unclear. This lack of time constraints led to a pivotal question: did deformation initiate in rocks that were altered by pre-kinematic hydration, or was hydration syn-kinematic and driven by the formation of Alpine shear zones? Laser ablation ICP-MS U–Pb geochronology applied to epidote in hydrothermal veins provides new evidence for pre-Alpine hydration of granitoids in the Aar Massif. Two veining events are recognized: 1) one at ca. 276 Ma occurring during Permian transtension and rifting, and 2) another at ca. 14 Ma related to late Alpine exhumation phases. Initial 207Pb/206Pb ratios of all Permian epidote samples overlap within uncertainty, suggesting only one fluid source and equilibration path. Also, these ratios are more radiogenic than those of the host rocks at the time of vein formation. The hydrogen isotopic composition of the Permian fluids was calculated from measurements in bulk epidote separates by high-temperature conversion elemental analyzer. With a temperature range of epidote crystallization estimated between 200–300 °C, the calculated δDfluid value is -57 to -44 ‰. An external source for the Permian fluids is suggested by the disequilibrium of Pb isotopes between hydrothermal epidote and host granitoids. Percolation of meteoric water along transtensional faults and interaction with syn-rift sediments before reaching the granitoids underneath is suggested by the Permian transtensional geodynamics and supported the hydrogen isotopic composition of the Permian fluids. The occurrence of Permian fluid circulation in the Aar Massif granitoids indicates that these rocks were altered before the onset of Alpine deformation. In fact, it can be inferred that fluid circulation caused not only veining, but also pervasive flow and thus the alteration of magmatic feldspar into fine-grained hydrous minerals throughout the granitoids of the present-day Aar Massif. This enabled pre-Alpine storage of water and the creation of numerous additional grain boundaries, both favoring viscous granular flow during Alpine deformation. In this context, the localization of strain in polymineralic aggregates containing hydrous minerals can recycle stored pre-kinematic (i.e., Permian) water. Thus, the initiation of Alpine deformation did not necessarily require the addition of syn-kinematic (i.e., Alpine) fluids, although their presence is confirmed by this and previous studies (e.g., Ricchi et al., 2019; Peverelli et al., 2021).

The low permeability and excellent sealing properties of mudstones places such sedimentary rocks into focus of geo-engineering applications using clay formations as natural barriers for contaminant transport. Here we investigate microscale deformation structures in the Opalinus Clay in northern Switzerland, which is currently under investigation as a host rock for radioactive waste confinement. We aim to characterize paleo-faulting/fracturing as well as subsequent mineralization events. For this purpose, drill core samples were investigated macroscopically, as well as by low and high-resolution optical light microscopy and scanning electron microscopy (SEM) in combination with energy-dispersive X-ray spectroscopy (EDX). These data were combined with high-resolution trace element maps obtained by Synchrotron X-ray Fluorescence Microscopy (SXFM). By combining the observed microstructures with the micro-chemistry of the associated mineralization events we yield a grouping between different processes and, in combination with cross-cutting relationships, a relative timing of the different deformation events.Commonly the Opalinus Clay is weakly deformed, with only few localized deformation structures. The latter include: calcite veins (mm thickness) as well as mineralized (calcite and celestite) thrusts, normal faults and strike-slip faults (all cm thickness). Micro-textural analysis shows that low-angle thrust faults with calcite slickensides on their dip-slip surfaces localize on pre-existing horizontal fibrous calcite veins. The horizontal veins imply an early deformation stage with temporarily high pore fluid pressures under sub-horizontal max. principle stresses within the highly anisotropic mudstone. The first order analysis of the major element chemistry between the calcite forming slickensides and the fibrous veins shows significant differences in Mg, Fe and Mn contents. From a fluid-mechanical perspective, this finding implies that generation of fibrous veins during a first fluid event provides the mechanical discontinuity, which is reused during later fluid assisted thrusting.The overprinting relationships between fibrous veins and slickensides indicate that deformation/precipitation events occurred in a cyclic fashion. Such information is a key towards the understanding of fluid assisted deformation and mineralization processes in compacted and anisotropic clay formations. On a regional scale, variations in paleo-deformation-mineralization events in the Opalinus Clay imply regional differences likely related to a gradually varying intensity of compressional (thrusting) and extensional (normal faulting) tectonics throughout northern Switzerland.

Rheological models of Earth’s granitoid mid- to upper crust are commonly based on the physico-chemical properties of the most abundant rock forming minerals quartz and feldspar. However, there is increasing field evidence that deformation in these rocks localizes in ultrafine-grained polymineralic shear zones, which are weaker than any of the end member minerals. Especially at the brittle to viscous transition, the localization and deformation mechanisms, i.e. the role of incipient brittle deformation vs. continuous viscous strain localization, is not yet fully understood.To fill this gap in knowledge, ultramylonite samples with granitic composition from the Central Aar Granite (Aar Massif, Central Switzerland) were deformed using a Griggs type apparatus. The foliation of the ultramylonitic starting material was oriented 45° to the compression direction, to investigate the influence of grain size and composition on strain localization in the different mylonite bands. Two types of coaxial experiments were conducted at 650°C, and 1.2 GPa confining pressure: A) Discrete fractures were created before the shear deformation starts; B) No fractures were induced during an early stage of the experiment.All experiments have in common that strain is accommodated in 20-100 µm wide viscous shear zones with elongated grains and minor grain size reduction. In these shear zones, most strain is further localized in 10-20 µm wide zones, showing dramatic grain size reduction down to few tens of nanometres. In the experimentally generated shear zones, both, brittle and viscous processes are active. In terms of overall rock strength, all newly formed ultrafine-grained shear zones are up to three times weaker than comparable experiments on pure quartz or coarser grained granites – which agrees well with field observations. Furthermore, pre-fractured type A) is up to two times weaker than the non-fractured type B), and the orientation and number of shear zones is also fundamentally different between the two experiment types.This study confirms two weakening factors promoting different types of strain localization at the brittle to viscous transition: 1) The existence of fractures and their interconnectivity – facilitating highly-localized grain size reduction; 2) Initial sample heterogeneity by polymineralic composition and ultrafine grain size – generating grain size reduction along strain gradients by activating viscous processes. Further quantitative microstructural analyses will reveal the role of chemistry and the deformation mechanisms on the localization behaviour.

Owing to still ongoing convergence within the Europe-Adria collision zone, Switzerland is affected by heterogeneously distributed moderate seismic activity. The project SeismoTeCH aims to improve the understanding of the links between the seismic activity, existing fault structures and geodynamics in Switzerland and its close vicinity. We started with compiling existing databases on faults (fault densities, lengths and orientations), seismic activity (spatial hypocenter and magnitude distributions, detection of seismic lineaments, focal mechanisms), orientations of mean principal stress axes and recent crustal movements (GNSS, high precision levelling) in order to establish potential correspondences as well as regional variations.Due to the long-lasting Alpine deformation, fault-orientation patterns as well as fault-densities vary between specific tectonic domains (Jura/North-Alpine foreland, Alpine frontal sediment nappe systems, External Crystalline Massifs, inner-Alpine domains and Southern Alps). Despite this variability, the fault patterns show first order correlations with the spatial arrangement of newly mapped seismic lineaments, earthquake focal planes and associated focal mechanisms. This correlation indicates a regional geodynamics-controlled reactivation of the specific fault networks during current crustal movements. In terms of recent surface movements, variations in (i) horizontal GNSS movements with respect to stable Europe and (ii) vertical uplift (from levelling and GNSS data) have to be discriminated. (i) From E to W in southern Switzerland (S-Grisons–Ticino–Valais, S of Rhone-Simplon line), horizontal movements change from NW to SW directions (velocities >0.5-0.8mm/yr). The southern Adria crustal block shows minimal to no lateral motions in the W-part and a clear NE-directed motion that is progressively increasing towards the E. This motion can be correlated with the so-called counter-clockwise rotation of the Adriatic plate. North of aforementioned domain, N- to NW-directed movements dominate but velocities decrease progressively from the central Alpine domains (Overall recent surface movements and seismicity in and along Central Alpine crustal blocks are affected by buoyancy-driven vertical combined with transpressional/-tensional horizontal movements indicating a lithosphere-scale geodynamic forcing.

Studies on moment tensors (MT) and focal mechanisms are of great importance for assessing regional and local seismotectonic processes, especially when a high-quality, dense network is in operation. However, common MT inversion methods are largely restricted to magnitudes > 3.5. In order to lower the completeness of MT catalogs, improved Green’s functions and/or hybrid inversion techniques are needed. In this study, we revisit small-to-moderate earthquakes, which occurred in Switzerland and surrounding regions by means of various MT inversion methods and assess the potential to improve completeness of MT catalogs in Central Alps region. To accomplish this, we implement state-of the art methods for MT inversion using either full waveform data or combinations of first-motion polarities with amplitudes and amplitude ratios. Methods based on full waveform inversion considered in this study are ISOLA (Sokos & Zahradnik 2013) and Grond (Heimann et al. 2018), as well as techniques based on amplitudes and/or polarities (HybridMT (Kwiatek et al. 2016), MTfit (Pugh & White 2018)), which can solve MTs for smaller magnitude earthquakes. Hence, the combination of multiple techniques allows to compute full or deviatoric MTs for a broader range of magnitudes and enrich the existing catalogs. We first apply these methods to recent earthquake sequences occurred in the Central Alps between 2019 and 2021. During that period, several earthquake sequences, like the one associated with the 2021 M4.1 Arolla earthquake, occurred and show complexity on the waveforms, due to their shallow focal depths. In addition, several of the standard MT solutions calculated by the Swiss Seismological Service (SED) for these earthquakes indicate complex moment tensors with unusually high percentage of the CLVD component. To check whether such CLVD component is real and not an artifact caused, for instance, by unmodeled heterogeneities, we invert for full and deviatoric MTs using multiple 1D velocity models and algorithms. Additionally, we perform MT inversions for several earthquakes either within selected earthquake sequences or regional background seismicity. The resulting MT solutions are compared to existing high-quality focal mechanisms computed using first motion polarities as well as to high-precision double difference locations. Uncertainties of MT solutions are estimated using bootstrap-based methods. This work contributes towards an enriched high-quality focal mechanisms database for Switzerland, which could be used to revisit the regional to local stress field at unprecedented resolution and provides new insights into the complexities of active fault systems in the Central Alps region., EGUsphere

Accuracy of hypocenter location, in particular focal depth, is a precondition for high-resolution seismotectonic analysis of natural and induced seismicity. For instance, linking seismicity with mapped fault segments requires hypocenter accuracy at the sub-kilometer scale. In this study, we demonstrate that inaccurate velocity models and improper phase selection can bias absolute hypocenter locations and location uncertainties, resulting in errors larger than the targeted accuracy. To avoid such bias in densely instrumented seismic networks, we propose a coupled hypocenter-velocity inversion procedure restricted to direct, first-arriving, mainly upper-crustal Pg and Sg phases. On the basis of synthetic tests and selected ground-truth events we demonstrate that a sub-kilometer hypocenter accuracy can be achieved by regional-scale, three-dimensional Pg and Sg velocity models combined with dynamic phase selection and a non-linear location algorithm. The tomographic inversion uses about 60’000 Pg and 30’000 Sg quality-checked phases of local earthquakes in the Central Alps region. The derived models image the VP and VS structure of the Central Alps upper crust at unprecedented resolution, including small-scale anomalies such as those caused by Subalpine Molasse units below the Alpine front. The relocation procedure is applied to more than 18’000 earthquakes and the relocated hypocenters reveal previously unrecognized seismogenic structures, for instance in the Swiss Molasse basin south of Bern. The ML 4.6 Urnerboden earthquake of 2017 is used as an example to demonstrate how the derived 3-D velocity structure and relocated hypocenters can be jointly interpreted to constrain the lithology hosting upper-crustal seismicity in the Central Alps., Journal of Geophysical Research: Solid Earth, 126 (12), ISSN:2169-9313, ISSN:0148-0227, ISSN:2169-9356

Ice-driven mechanical weathering in mountainous environment is considered as an efficient process for slow but cyclical mechanical preconditioning of rockfall events. In this study, we simulate subcritical microfracture propagation under frost wedging conditions along pre-existing mechanical weaknesses of intact rock bridges with an innovative experimental approach. Two series of freeze–thaw experiments conducted in an environmental chamber were carried out to investigate and monitor the propagation of artificially induced fractures (AIF) in two twin gneiss samples. A displacement sensor recorded the sample’s in situ deformation in an environmental chamber during the experiments. 3D X-ray CT scans, performed before and after the experiments, as well as thin sections showing the post-experiment state of the deformed samples allowed tracking and quantification of fracture propagation. Our results demonstrate that frost wedging propagated the AIFs 1.25 cm2 and 3.5 cm2 after 42 and 87 freeze–thaw cycles, respectively. The experiments show that volumetric expansion of water upon freezing, cooperating with volumetric thermal expansion and contraction of the anisotropic rock, plays a key role in fracture widening and propagation. Based on these results, this study proposes that: (1) frost wedging exploits intrinsic pre-existing mechanical anisotropies of the rock; (2) the fracturing process is not continuous but alternates between stages of fast propagation and more quiet stages of stress accumulation; and (3) downward migration of “wedging grains,” stuck between the walls of the fracture, increases the tensile stress at the tip, widening and propagating the fractures with each freeze–thaw cycle. The experimental design developed in this study offers the chance to visualize and quantify the long-term efficiency of frost wedging in near-natural scenarios.

Epidote – here defined as minerals belonging to the epidote–clinozoisite solid solution – is a low-μ (μ=238U/204Pb) mineral occurring in a variety of geological environments and participating in many metamorphic reactions that is stable throughout a wide range of pressure–temperature conditions. Despite containing fair amounts of U, its use as a U−Pb geochronometer has been hindered by the commonly high contents of initial Pb, with isotopic compositions that cannot be assumed a priori. We present a U−Pb geochronology of hydrothermal-vein epidote spanning a wide range of Pb (3.9–190 µg g−1), Th (0.01–38 µg g−1), and U (2.6–530 µg g−1) contents and with μ values between 7 and 510 from the Albula area (eastern Swiss Alps), from the Grimsel area (central Swiss Alps), and from the Heyuan fault (Guangdong Province, China). The investigated epidote samples show appreciable fractions of initial Pb contents (f206=0.7–1.0) – i.e., relative to radiogenic Pb – that vary to different extents. A protocol has been developed for in situ U−Pb dating of epidote by spot-analysis laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) with a magmatic allanite as the primary reference material. The suitability of the protocol and the reliability of the measured isotopic ratios have been ascertained by independent measurements of 238U/206Pb and 207Pb/206Pb ratios, respectively, with quadrupole and multicollector ICP-MS applied to epidote micro-separates digested and diluted in acids. For age calculation, we used the Tera–Wasserburg (207Pb/206Pb versus 238U/206Pb) diagram, which does not require corrections for initial Pb and provides the initial 207Pb/206Pb ratio. Petrographic and microstructural data indicate that the calculated ages date the crystallization of vein epidote from a hydrothermal fluid and that the U−Pb system was not reset to younger ages by later events. Vein epidote from the Albula area formed in the Paleocene (62.7±3.0 Ma) and is related to Alpine greenschist-facies metamorphism. The Miocene (19.2±4.3 and 16.9±3.7 Ma) epidote veins from the Grimsel area formed during the Handegg deformation phase (22–17 Ma) of the Alpine evolution of the Aar Massif. Identical initial 207Pb/206Pb ratios reveal homogeneity in Pb isotopic compositions of the fluid across ca. 100 m. Vein epidote from the Heyuan fault is Cretaceous in age ( 107.2±8.9 Ma) and formed during the early movements of the fault. In situ U−Pb analyses of epidote returned reliable ages of otherwise undatable epidote–quartz veins. The Tera–Wasserburg approach has proven pivotal for in situ U−Pb dating of epidote, and the decisive aspect for low age uncertainties is the variability in intra-sample initial Pb fractions.

In the framework of the SeismoTeCH project, which aims at advancing our understanding of seismotectonic processes in Switzerland, we present the first 3-D attenuation model of the upper crust for the Central Alps and their northern foreland. The 3-D inversions derive the quality factor Q (1/attenuation) using path attenuation t∗ observations for 4,192 distributed earthquakes recorded on permanent and temporary stations, including both velocity and acceleration records for the period 2002-2019. We followed a procedure of gradational inversions, in which a series of inversions with increasingly grid complexity are performed, with the goal of obtaining a useful Q model everywhere despite the varied data distribution. The Qs and Qp results show large-scale features in the upper crust, which are consistent with a recently improved high-resolution velocity models of the same region and serve to refine the interpretations of crustal structures from Vp and Vp/Vs. For example, the foreland region of southern Germany and northern Switzerland show a low Q crustal block bounded by high Q regions in the uppermost layer between -2.5 and 2.0 km depth. This markedly correlates with the overlying surface geology, where low Q areas coincide with the Molasse Basin, and the transition between low and high Q regions outline the geological boundary between the Molasse and the Mesozoic sediments towards north and the Alpine front to the south. At depths ranging between 2.0 - 6.5 km, low Q is imaged along the Rhone valley in the Valais in southwest Switzerland. This region presents the transition between the Centrals and Western Alps and hosts the presently seismically most active fault zones. As the attenuation of fractured areas is enhanced by fluids, low Q values may relate here to distributed microfractures that produce greater fracture connectivity and permeability in a relatively higher strain-rate zone. These geophysical constraints seem to support crustal scale fluid flow along fracture networks as manifest by the prominent occurrence of hot springs in this area. On the other hand, the moderate-to-high Qs and Qp (400-800) along with low Vp/Vs ratio and high Vs observed in the external Aar Massif could be indicative of metamorphic processes leading to different Vp/Vs ratios compared to the basement in the northern foreland (Black Forest Massif), and possibly image the continuation of the massif 20-30 km further to the northeast. In combination with recently developed Vp and Vs velocity models, the developed 3-D attenuation models provide additional constraints in terms of composition and physical properties of the uppermost crust of the central Alps as well as crucial input for next generation seismic hazard models of Switzerland, allowing for a more realistic prediction of earthquake related ground motions.

Recently, the application of LA–ICP–MS has enabled U–Pb dating of epidote minerals within the epidote–clinozoisite solid solution series (Peverelli et al., 2020). Epidote crystallization ages can provide an absolute time frame of deformation sequences when combined with detailed microstructural and metamorphic P–T analysis. Epidote deformation occurs in a brittle manner over a wide range of conditions below its closure temperature for Pb diffusion (685–750 °C; Dahl, 1997); hence, such deformation will not affect its formation U–Pb age. Nevertheless, the possibility of isotopically resetting epidote via fluid–mineral interaction has to be taken into account even at low deformation temperatures.We investigated the geochemical and Sr–Pb isotopic characteristics of epidote in one hydrothermal vein in the Aar Massif (central Swiss Alps). The vein is associated with an Alpine shear zone and it is composed of aggregates of 0.1–1 mm anhedral to subhedral epidote grains (epidote-A) + green biotite within a quartz matrix. This quartz dynamically recrystallized by subgrain rotation at temperatures above 400 °C (Stipp et al., 2002) along with crystallization of a second epidote generation (epidote-B) made of tiny (< 0.1 mm) anhedral epidote grains in part mantling epidote-A and defining a fold. We address whether interaction with the fluid that precipitated epidote-B chemically affected epidote-A, i.e. whether the U–Pb age measured by LA–ICP–MS in epidote-A still dates its crystallization upon vein formation or displays age disturbance.LA–ICP–MS Sr and Pb concentration data overlap between epidote-A and epidote-B, as do their REE patterns, with (La/Yb)N ratios of 0.03–0.92. Lead and Sr isotopic signatures were measured respectively by solution MC–ICP–MS and by TIMS in epidote-A and in separates mixing different proportions of epidote-A and -B (no pure mechanical separates of epidote-B possible), and they are different. This requires open-system conditions during deformation, i.e., introduction of an external fluid with higher 87Sr/86Sr and 208Pb/206Pb ratios during crystallization of epidote-B. Despite the presence of an external fluid and the incorporation of external Sr and Pb in epidote-B, LA–ICP–MS U–Pb isotopic data for epidote-A define a regression in a Tera–Wasserburg plot indicating an age of 19.2 ± 4.3 Ma, consistent with epidote-A crystallization during original vein opening. The preservation of the crystallization age in epidote-A indicates that interaction with the fluid that formed epidote-B did not geochemically and isotopically affect epidote-A. The consistency in trace element contents between epidote-A and -B hints that the epidote-forming cations were inherited by the fluid from epidote-A, and thus suggests dissolution-precipitation as the formation process for epidote-B. Dahl, Earth Planet. Sci. Lett. 150, 277–290, 1997.Peverelli et al., Geochronology Discuss. [preprint], https://doi.org/10.5194/gchron-2020-27, in review, 2020.Stipp et al., Geological Society, London, Special Publications, 200(1), 171-190, 2002

Inversion of passive margins and their transportation into fold-and-thrust belts is a critical stage of mountain building processes and their structural interpretation is fundamental for understanding collisional orogens. Due to the multitude of parameters that influence their formation (e.g. the interaction between sedimentary cover and basement, the mechanical stratigraphy or the rheology of different rock types) as well as along-strike internal variations, a single cross-sectional view is insufficient in exploring the 3D evolution of a fold-and-thrust belt. Hence, a 3D geological characterization is required to better comprehend such complex systems. Based on a detailed digital map, a 3D structural model of the current tectonic situation and sequential retrodeformation, we elaborate the 3D evolution of a part of the former European passive continental margin. In this setting, we focus on the Doldenhorn Nappe (DN) and the underlying western Aar massif (external Central Alps, Switzerland). The DN is part of the Helvetic nappe system and consists of a large-scale recumbent fold with a thin inverted limb of intensively deformed sediments (Herwegh and Pfiffner 2005). The sedimentary rocks of the DN were deposited in Mesozoic-Cenozoic times in a small-sized basin, which has been inverted during the compression of the Alpine orogeny (Burkhard 1988). Along NNW-SSE striking geological cross-sections, restoration techniques reveal the original asymmetric triangular shape of the DN basin and how the basin has been exhumed from ~ -12 km (Berger et al. 2020) to its present position at 4km elevation above sea level throughout several Alpine deformation stages. Moreover, the model allows to visualize the current structural position of the DN and the massif as well as the geometric and overprinting relationships of the articulated deformation sequence that shaped the investigated area throughout the Alpine evolution. Here we document that: (i) the DN is a strongly non-cylindrical recumbent fold that progressively pinches out toward the NE; (ii) significant along-strike (W-E) stratigraphy thickness variations are reflected in structural variations from a single basal thrust deformation (W) to an in-sequence thrust deformation (E); and (iii) the progressive exhumation of the basement units towards the E and thrusting towards the N. In this context, special emphasis is given to illustrate how three-dimensional geometry of inherited pre-orogenic structures (e.g., Variscan-Permian and rifting related basement cover structures) play a key role in the structural style of fold-and-thrust belts. In summary, today’s structural position of the DN is the result of the inversion of a small basin in an early stage of thrusting, which was followed by sub-vertical buoyancy driven exhumation of the Aar massif and subsequent thrust related shortening. All three stages are deeply coupled with an original non-cylindrical shape of the former European passive continental margin.

The Rawil depression north of the Rhone Simplon fault zone (southwestern Swiss Alps) was host of the Mw = 5.8 Sion earthquake in 1946 (Fäh et al., 2011). It is nowadays one of the seismically most active regions in Switzerland and seismicity forms a cluster, which is elongated approximately in WSW-ENE direction over 40-50 km. In November 2019, a remarkable earthquake sequence occurred within the center of this cluster north of the village of Anzère, with more than 300 earthquakes up to ML = 3.3 recorded by the Swiss Seismological Service within 20 days.Detecting associated full-scale 3D fault patterns solely based on earthquake hypocenters is challenging because of commonly too limited spatial resolution and insufficient number of seismic events. Within the framework of SeismoTeCH, we aim to improve these limitations by a combination of high-precision hypocenter relocation techniques, reconstruction of subsurface fault patterns and correlative links between surface and subsurface data. Assuming that a fault is seismically active multiple times and that the seismic stress-release is initiated at different locations along the fault, we can calculate 3D fault plane orientations from the hypocenter locations. Together with the 17 focal mechanisms derived for the Anzère sequence, we are able to gain geometrical and kinematic information of the seismic faults in 3D. Our analysis reveals a seismically active transpressional step-over structure within a dextral strike-slip fault zone. With remote sensing and field observations, we detect exhumed faults with similar orientations and kinematics that presumably represent step-over structures, interconnecting previously known strike-slip fault zones.Although seismic activity occurs at depths between 3-5 km, we conclude that the observed surface fault systems in the Rawil depression can be correlated in terms of fault patterns with those assumed at depth. The linkage of the recent seismicity with structural observations of exhumed, potentially paleo-seismic faults in combination with recent hypocenter relocation techniques therefore have great potential to provide further insights into fault linkage and earthquake rupturing processes. ReferencesFäh, D., Giardini, D., Kästli, P., Deichmann, N., Gisler, M., Schwarz-Zanetti, G., Alvarez-Rubio, S., Sellami, S., Edwards, B., Allmann, B., Bethmann, F., Wössner, J., Gassner-Stamm, G., Fritsche, S., Eberhard, D., 2011. ECOS-09 Earthquake Catalogue of Switzerland Release 2011. Report and Database. Public catalogue, 17.4.2011. Swiss Seismological Service ETH Zürich, Report SED/RISK/R/001/20110417.

Accuracy of hypocenter location, in particular focal depth, is a precondition for high-resolution seismotectonic analysis of natural and induced seismicity. For instance, linking seismicity with mapped fault segments requires hypocenter accuracy at the sub-kilometer scale. In this study, we demonstrate that inaccurate velocity models and improper phase selection can bias absolute hypocenter locations and location uncertainties, resulting in errors larger than the targeted accuracy. To avoid such bias in densely instrumented seismic networks, we propose a coupled hypocenter-velocity inversion restricted to direct, upper-crustal Pg and Sg phases. The derived three-dimensional velocity models, combined with dynamic phase selection and non-linear location algorithms result in a highly accurate earthquake catalog, including consistent hypocenter uncertainties. We apply this procedure to about 60’000 Pg and 30’000 Sg quality-checked phases of local earthquakes in the Central Alps region. The derived tomographic models image the Vp and Vs velocity structure of the Central Alps’ upper crust at unprecedented resolution, including small-scale anomalies such as those caused by a Permo-Carboniferous trough in the northern foreland, Subalpine Molasse below the Alpine front or crystalline basement units within the Penninic nappes. The external Aar Massif is characterized by low Vp/Vs ratios of about 1.625-1.675 in the depth range of 2-6.5 km, which we relate to a felsic composition of the uplifted crustal block, possibly with increased quartz content. Finally, we discuss along-strike variations imaged by relocated seismicity in the Central Alps and demonstrate how joint interpretation of velocity structure and hypocenters provides additional constraints on lithologies of upper-crustal seismicity.

During the accretion of foreland basin sediments into an accretionary or orogenic wedge, the sediments dehydrate and deform. Both dehydration and deformation intensity increase from the outer to the inner wedge and are a function of metamorphic processes and strain. Here, we study the microstructural evolution of slates from the exhumed Flysch units making up a paleo accretionary wedge in the European Alps. With classic SEM imaging and synchrotron X-ray fluorescence microscopy, we document the evolution of slate fabrics and calcite veins and aim at understanding the role of the evolving slate fabrics for strain localisation and fluid flow at the micro-scale.The investigated slate samples are from a NW-SE transect covering the outer and inner wedge from 200 to 330 °C. The metamorphic gradient directly correlates with an increasing background strain gradient. With the use of the autocorrelation function, we quantify the evolution of the microfabrics along the metamorphic gradient and document deformation stages from a weakly deformed slate without foliation in the outer wedge to a strongly deformed slate with a dense spaced foliation in the inner wedge. The foliation mainly forms by dissolution-precipitation processes, which increase with increasing metamorphic gradient.The slate matrix reveals two main sets of veins. The first vein set includes micron-scaled calcite veinlets with very high spatial densities. The second vein set includes layer parallel calcite veins that form vein-arrays (couple of metres thick) in the inner wedge. Both vein sets could have moved large amounts of fluids through the wedge. The spatial distribution of the micron-veinlets reveals that fluids were moved pervasively. In the case of the layer parallel veins forming vein-arrays, fluid flow was localized, supported by the dense spaced foliation formed in the slate matrix in the inner wedge. This way we now establish a direct link between the microstructural evolution in the slate matrix and associated dehydration, where fluids become increasingly channelled towards the inner wedge. Knowing that the vein-arrays have length dimensions in the order or hundreds of metres to kilometres, these structures are important for larger-scale fluid flow, the feeding of fluids into megathrusts and for related seismic activity in the wedge.

Journal of structural geology 150, 104404 (2021). doi:10.1016/j.jsg.2021.104404, Published by Elsevier Science, Amsterdam [u.a.]

White mica has been widely used to date microstructures and tectonic events in faults, shear-zones and folds because of its suitability for radiogenic dating. However, complex (i) microstructural evolution, (ii) individual chemical evolution of the K-bearing phases, (iii) mixing of ‘detrital grains’ with newly formed and/or recrystallized or chemically reset grains as well as (iv) volume diffusion may result in apparent K-Ar ages. Here, specimens from a prograde sediment sequence of the exhumed fossil European Alpine accretionary wedge were used to investigate resetting processes of white mica by the type and intensity of deformation as well as peak metamorphic conditions. We combine the K-Ar system with mass and mineral quantities from grain size fractions to calculate the amount of recrystallized white mica in each grain size fraction along the metamorphic gradient. Increasing recrystallization with increasing metamorphic grade is related to thermally activated pressure solution and dissolution-precipitation creep, as seen by the formation of a spaced cleavage of recrystallized phyllosilicates documented through Synchrotron X-ray Fluorescence Microscopy and Scanning Electron Microscope imaging techniques. Increasing recrystallization by dissolution-precipitation processes induces chemical resetting of the isotopic system, resulting in a prograde decrease of apparent K-Ar ages. We demonstrate that Ar volume diffusion does not play a significant role for the low-temperature samples, promoting recrystallization as the important physico-chemical process for age resetting. However, white mica chemistry reveals that no simple relation between isotopic resetting and grain size exists along the prograde path. Reliable age information can therefore only be obtained in the case of (nearly) complete resetting, which accounts only for the smallest grain size fraction at the highest metamorphic temperature. These findings could shed new light on accurate dating of mica-rich fault rocks, where the time constraints depend not only on the temperature, but also on the amount and type of deformation.

Monoclinic epidote is a low-µ (µ = 283U / 204Pb) mineral occurring in a variety of geological environments, participating in many metamorphic reactions and stable throughout a wide range of pressure–temperature conditions. Despite containing fair amounts of U, its use as a U–Pb geochronometer has been hindered by the commonly high contents of initial Pb with isotopic compositions that cannot be assumed a priori. We present U–Pb geochronology of hydrothermal-vein epidote spanning a wide range of Pb (3.9–190 µg g−1), Th (0.009–38 µg g−1) and U (2.6–530 µg g−1) contents and with µ values between 7–510 from the Albula area (eastern Swiss Alps), from the Grimsel area (central Swiss Alps) and from the Heyuan fault (Guangdong province, China). The investigated epidote samples show appreciable fractions of initial Pb that vary to different extents. A protocol has been developed for in situ U–Pb dating of epidote by spot-analysis laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) with a magmatic allanite as primary reference material. The suitability of the protocol and the reliability of the measured isotopic ratios have been ascertained by independent measurements of 238U / 206Pb and 207Pb / 206Pb ratios respectively by quadrupole and multicollector ICP–MS applied to epidote micro-separates digested and diluted in acids. For age calculation, we used the Tera–Wasserburg (207Pb / 206Pb–238U / 206Pb) diagram, which does not require corrections for initial Pb and provides the initial 207Pb / 206Pb ratio if all intra-sample analyses are co-genetic. Petrographic and microstructural data indicate that the calculated ages date the crystallization of vein epidote from a hydrothermal fluid and that the U–Pb system was not reset to younger ages by later events. Vein epidote from the Albula area formed in the Paleocene (62.7 ± 3.0 Ma) and is related to Alpine greenschist-facies metamorphism. The Miocene (19.1 ± 4.0 Ma and 16.9 ± 3.7 Ma) epidote veins from the Grimsel area formed during the Handegg phase (22–17 Ma) of the Alpine evolution of the Aar Massif. Identical initial 207Pb / 206Pb ratios reveal homogeneity in Pb isotopic compositions of the fluid across ca. 200 m. Vein epidote from the Heyuan fault is Cretaceous in age (108.1 ± 8.4 Ma) and formed during the early movements of the fault. In situ U–Pb analyses of epidote returned reliable ages of otherwise undatable epidote-quartz veins. The Tera–Wasserburg approach has proven pivotal for in situ U–Pb dating of epidote and the decisive aspect for low age uncertainties is the variability in intra-sample initial Pb fractions.

In experiments designed to understand deep shear zones, we show that periodic porous sheets emerge spontaneously during viscous creep, forming a hydro-mechanical anisotropy that influences mass transfer. These findings challenge the current paradigm of viscosity in solid rocks. In particular, they showcase how shear zones may actively focus mass transport and highlight the possibility that viscous rocks could locally transition from flow to fracture. Our work demonstrates that viscosity in solids is not directly comparable to viscosity in fluids and this is consequential for a range of important solid Earth topics, like slow earthquakes, the flow of glacial ice and the tectonics of exoplanets.

The present-day deformation in the Central Alps is dominated by vertical uplift, at rates up to 1.5 mm/yr as indicated by high-precision levelling and GPS data. Understanding the driving mechanisms of this neotectonic uplift and its link to seismicity in the Central Alps requires accurate locations of current deformation processes within the upper crust. Especially the question if and how deformation in the crystalline basement is coupled with deformation in the overlaying nappe systems is key to understand the neotectonic processes. Seismicity provides important information on deformation in the uppermost crust, however, an accuracy of focal depths in the order of few kilometers and less is required to distinguish sources in the basement from sources in the sedimentary cover.In this study, we demonstrate how insufficient crustal velocity models and inconsistent seismic phase selection can lead to biased hypocenter solutions, which hamper such high-resolution seismotectonic interpretations. We propose a relocation procedure combining a new high-resolution Pg and Sg 3D crustal model of the Central Alps with a dynamic seismic phase selection to overcome this bias and to improve accuracy of hypocenter solutions. The new tomographic model is based on more than 60,000 Pg and 30,000 quality-checked Sg phases of earthquakes, which occurred in the greater Central Alpine region between 1996 and 2019. In combination with a nonlinear, probabilistic earthquake location algorithm, the model was used to relocate more than 18’000 earthquakes, which occurred in this region over the past 36 years. The derived catalog includes a consistent error and quality assessment, calibrated against ground-truth events like quarry blasts.The relocated seismicity in the Central Alps is interpreted together with additional information from the tomographic model, focal mechanisms, geophysical, geological and geodetic data. We focus our interpretation on the eastern Aar massif as well as on the Rawil depression, located in-between the outcropping Aar and Aiguilles-Rouge massifs. Both regions were recently affected by remarkable seismic events. The ML4.6 Urnerboden earthquake of 2017 occurred near the eastern termination of the Aar massif, while a sequence of about 350 events occurred in the Rawil earthquake lineament near the Sanetschpass in November 2019. Both sequences provide unique insights into active faults in the Central Alps and we image systems of sub-vertically oriented strike-slip faults of variable strike, which root in the crystalline basement in both regions. Our results document the existence of active strike-slip fault systems in the External Crystalline Massifs of the Central Alps in regions of maximum change in uplift rates. We therefore discuss possible models relating the observed strike-slip kinematics to buoyancy-driven vertical tectonic processes.

The seismic-interseismic cycle strongly relates to the interplay between dilation owing to fracturing and frictional granular flow on one hand side and hydrothermal cementation processes on the other side. This study investigates different fault rocks of a crustal-scale fault zone in the Central Alps (Switzerland). We combine microstructural with geochemical approaches to decipher the interaction of grain size reduction via frictional processes with precipitation and resulting particle size increases. The three major fault rocks, i.e. (1) cockade-bearing breccias, (2) cataclasites, and (3) fault gouges, differ in their microstructure. The chemical data clearly demonstrate a decreasing gain of volume along this group of tectonites. Their different precipitation volumes most likely relate to dynamic changed of the local permeability of these rocks. The fluid pathways control the precipitation at different localities and times, which affect the healing of these fault rocks inducing a gain in rock strength. During the next deformation event, the extent of healing therefore directly controls the mechanical behavior of the rock. The estimated volume gain (~+110%) in cockade-bearing breccias is consistent with the seismic dilatant behavior of these frictional rocks as already proposed from other arguments (Berger and Herwegh 2019). This is in contrast to the fault-gouges with only minor gains in volume and mass resulting in a predominantly non-cohesive deformation style. This example indicates that permeability evolution (and related hydrothermal processes) strongly influences the mechanical behavior of such faults. This shows the highly dynamic behavior with time in long-lived fault systems. These dynamic changes in precipitation and resulting different strengths occur at different timescales from minutes (seismic events) to thousands of years.Ref.: Berger, A., Herwegh, M., 2019. Cockade structures as a paleo-earthquake proxy in upper crustal hydrothermal systems. Nature Scientific Reports, 9, 9209.

A paleohydrothermal giant quartz reef (at least 75 m wide, 40 km long) and abundant hot springs at the Heyuan fault, South China, provide an excellent opportunity to investigate hydrothermal flows from the Mesozoic through to present-day.The giant quartz reef has formed in the extensional regime initiated in the Mesozoic, while a change to compressional stress on the Heyuan in the Cenozoic led to the development of cross-cutting strike-slip faults and associated vertical fracture network. Here, we present multiscale observations and analyses from the earlier long-term extensional phase.Detailed microstructural analyses identified a 'quartz-reef window' of formation occurring between ~200-350˚C, linking in both quasi-static criteria (accommodation space; massive fluid sources; and a cap rock/seal) and dynamic mechanisms (episodic-dynamic permeability; brittle-ductile cycles; and fluid injection though brittle-ductile equivalent of Sibson's 'fault-valve' behaviour.This oscillatory brittle-ductile fault-valve is recorded in the field through its apparent contradiction between idiomorphic 5 cm long quartz crystal growth in mode-I fractures, embedded at large-scale inside far from equilibrium fault zones with mylonitic and cataclastic microstructures. Another characteristic feature is the increasing quartz vein frequency towards the core shown by enrichment of SiO2, with depletion of K2O and Na2O in tectonites during alteration from the host granite; a reaction partly sourcing the SiO2 for the quartz reef.We present a first theoretical model compatible with the observation of oscillatory macroscale far from equilibrium conditions, followed by long periods of micro-scale local equilibrium. The model can in particular describe mechanisms of abundant SiO2 dominated fluid release reaching episodically above hydrostatic pressures followed by long periods of SiO2 precipitation, allowing growth of up to 5 cm long idiomorphic quartz crystals in subparallel open channels, which presumably were held open by high fluid pressures. In this interpretation, the observations instabilities are seen to stem from the multiscale and multiphysics of the mineral reactions at the brittle-ductile transition, promoted by a slow extensional geodynamic driver at the Heyuan fault.The new approach allows interpretation of rock physics properties in terms of recently discovered Thermo-Hydro-Mechanical-Chemical (THMC) multiscale wave-like instabilities. In the model short wavelength chemical dissolution-precipitation reaction waves are bouncing between the phyllonitic cap rock and the mylonitic shear zone below. A resonance phenomenon of constructive interference in a finite width around the future quartz-reef triggers the long-time scale steady-state attractor allowing quartz reef growth over geodynamic time scales. We show that this solitary wave limit forms a standing wave matching the characteristic periodic pattern of mode-I quartz veining around the reef and also explaining the fluid overpressures leading to local hydro-fracturing.

Ice-driven mechanical weathering in mountainous environment is considered an efficient process for slow preconditioning of rockfalls. In this study (Musso Piantelli et al., 2020), we simulate with an innovative experimental approach subcritical fracture-propagation under frost-wedging conditions through pre-existing weaknesses of intact rock bridges. Two series of freeze-thaw experiments in an environmental chamber have been designed to investigate and monitor the propagation of artificially-induced fractures (AIF) in two twin gneiss samples. By employing 3D X-Ray Computed Tomography and a displacement sensor, an accurate characterization and new insights into the fracture-propagation mechanism are provided. Our results demonstrate that frost wedging propagated the AIFs of 1.25 cm2 and 3.5 cm2 after 42 and 87 freeze-thaw cycles, respectively. The experiments show that volumetric expansion of water upon freezing, cooperating with volumetric thermal expansion and contraction of the rock, plays a key role in fracture widening and propagation. Based on these results, this study proposes that: (i) frost wedging exploits intrinsic pre-existing weaknesses of the rock; (ii) the fracturing process is not continuous but alternates propagation stages to phases of tensile stress accumulation; and (iii) downward migration of “wedging grains”, stuck between the walls of the fracture, increases the tensile stress at the tip, widening and propagating the fractures with each freeze-thaw cycle. The experimental design developed in this study offers the chance to visualize fracture-propagation in natural joints quantifying the long-term efficiency of this process in near-natural scenarios.REFERENCESMusso Piantelli, F., Herwegh, M., Anselmetti, F.S., Waldvogel, M., Gruner, U., (2020). Microfracture propagation in gneiss through frost wedging: insights from an experimental study. Natural Hazards, 1-18. https://doi.org/10.1007/s11069-019-03846-3

Late stage Alpine collision is the result of collision of the European and Adriatic continental plates. For this stage, particularly the External Crystalline Massifs (ECM) and their forelands provide information on deformation style/kinematics, exhumation history and geodynamic driving forces. Using the ECMs as marker for the non-thinned European passive continental margin, the margin’s paleogeography with its curved geometries controls in parts later compressional tectonics. During closure of the Valaisan in Eocene times, a major NW-SE trending sinistral transfer zone must have acted as lateral ramp between the NNW migrating Penninic front and the westerly situated European margin units (Argentera and Maritime Alps, stage 1 of Schmid et al. 2017). Hence, first transpressive movements were documented by transpressional strike-slip faults in the case of Argentera and southern Belledonne s.l. Massifs. The SW-NE trending ECMs became affected during a first stage of horizontal tectonics in Oligocene, when first the margin sediments were scrapped off their substratum (Helvetics, Chaînes Subalpines, Dauphinois) followed by thick-skinned thrusting. Transport directions gradually changed from W to N from the Western towards the Eastern ECMs. This spread in transport direction is the consequence of an Alp-internal vertical uplift (Internal Crystalline Massifs (ICM), Lepontine dome) induced by indentation of Adriatic mantle into European lithosphere (stage 2 of Schmid et al. 2017). With further down-bending of the European crust, a major change to a vertical tectonic deformation style occurred in Mid to Late Miocene. Steep reverse faults in the ECMs, in parts with oblique slip components (Mont Blanc), witness this stage with its enhanced vertical rock uplift component. The latter is most pronounced in the Aar Massif and gradually decreases towards the West (Belledonne s.l.). With continuous ICM exhumation in Late Miocene, deformation style switches again to horizontal tectonics, leading to ‘en bloc’ exhumation above basement thrusts of all massifs (Belledonne to Aar Massifs). Progressive shortening induced thrust propagation into the foreland sediments as well as the Jura mountains. In the case of the Argentera Massif, oroclinal bending probably led to a substantial anticlockwise rotation, which goes in hand with the rotation of the entire SW Alpine arc (stage 3 of Schmid et al. 2017). In a geodynamic context, the ‘Adriatic push model’ could explain aforementioned stages of classical horizontal tectonics. Not so, however, the observed severe components of vertical tectonics in the case of ECMs. Here, lower crustal delamination, with a loss of lithospheric negatively buoyant forcing and consequent strong increase in positive buoyancy explains generation and activity of steep reverse faulting in the ECMs. In this context, the ‘orogeny slab rollback’ model provides a physically more consistent framework to explain the observed deformation sequences of late-stage continent-continent collision.Schmid, S.M., Kissling, E., Diehl, T., van Hinsbergen, D.J., Molli, G., 2017. Ivrea mantle wedge, arc of the Western Alps, and kinematic evolution of the Alps–Apennines orogenic system. Swiss Journal of Geosciences 110, 581-612.

In plate-convergent settings, fluid-saturated sediments dehydrate during subduction. The sediments are subsequently accreted to the upper plate. Along their dehydration-deformation path, the initial unconsolidated soft marine sediments become thick, foliated, impermeable meta-sedimentary sequences. Fluid flow through such ‘non’-porous low-permeability rocks is concentrated in fracture networks, ranging from the mm- to the km-scale. We study the interplay between ductile and brittle deformation processes and fluid flow by investigating calcite veins in slates from the exhumed European Alpine accretionary wedge across scales (µm to km). These slates experienced peak metamorphic temperatures between 200°C and 330°C and represent the transition between the upper aseismic and seismic zone. With the use of Synchrotron X-ray Fluorescence Microscopy (SXFM), we investigate the slates by visualizing trace-element distributions. This technique shows that alternating cycles of slow pressure-dissolution processes and brittle fracturing persist over long time scales. At the micron-scale, pressure solution of the initial carbonate-rich slates is indicated by an enrichment of newly recrystallized phyllosilicates on cleavage planes and in pressure shadows. These ductile deformation features are mutually overprinted by calcite veins (aperture 10 µm), which are nicely visualized with Sr-SXFM maps. Increasing compaction and recrystallization in the slate-rich matrix leads to progressed dehydration resulting in an increased pore fluid pressure and subsequent hydrofracturing. The micron-sized fractures are immediately filled in with minerals, which are oversaturated at that time in the fluid, resulting in the formation of (i) micron-veinlets. Micron-veinlets collect (ii) into mm-cm sized veins, which themselves form (iii) vein arrays and (iv) mega-arrays, respectively at the 50-100 m and 300-400 m scale. This upscaling of fluid pathways indicates a localised fluid transport through the accretionary wedge, which has important implications for the understanding of the mechanical stability of the accretionary wedge and related seismic activity.

3D modelling of complex and irregular geological bodies is an expanding discipline that combines two-dimensional cartographic and structural data managed with GIS technology. This study presents a complete workflow developed to process geological information to build a 3D model of major stratigraphic, structural and tectonic boundaries. The investigated area is located in the western part of the Aar Massif (external Central Alps, Switzerland) characterized by pronounced topographic (600–(1) Generation of 2D polylines in a map view: a two-dimensional dataset of sequences of polylines has been generated in ArcGIS (10.3.1) defining the starting dataset for the major stratigraphic and tectonic boundaries of the bedrock units. This dataset has been compiled and integrated by using: (i) GeoCover vector datasets 1:25 000 of the Swiss Geological Survey; (ii) The Geological Special Map 1:100 000 of the Aar Massif and the Tavetsch and Gotthard Nappes of the Swiss Geological Survey; (iii) data from literature; and (iv) additional field work conducted for this study in key-locations.(2) Projection of 2D information onto 3D digital elevation model: with the 3D structural modelling software Move (Petex/Midland Valley; 2019.1) the boundaries have then been projected on a digital elevation model (swissALTI3D) with 2 m resolution.(3) Construction of tectonic cross sections: the use of geometric arguments as well as structural measurements allows for projection of these boundaries into a dense regularly spaced network of 2D cross-sections.(4) Interpolation of 3D surfaces: the surface and cross-sections boundaries can be interpolated by applying 3D projection and meshing techniques resulting in a final 3D structural model.Generally, steps (2–4) require iterative adaptations particularly in the case of surface areas being covered by glaciers or unconsolidated Quaternary sediments. In the model, special emphasis is given to visualize the current structural disposition of the western Aar Massif as well as the relative geometric and overprinting relationships of the deformation sequence that shaped the investigated area throughout the Alpine deformation. Finally, since in the investigated area underground data are scarce, an assessment of the relative uncertainties related to input data and is intended to be performed following the approach proposed by Baumberger (2014) and Ferńandez (2005). The workflow presented here offers the chance to gain validation approaches for domains only weakly constrained or with no surface data available, by generating a 3D model that integrates all accessible geological information and background knowledge. REFERENCESBaumberger, R. (2014): Quantification of Lineaments: Link between internal 3D structure and surface evolution 328 of the Hasli valley (Aar massif, central alps, Switzerland), University of Bern, PhD Thesis, unpublished.Ferńandez, O. (2005): Obtaining a best fitting plane through 3D georeferenced data, Journal of Structural Geology 27, pp. 855–858

We revisit large shear strain deformation experiments on Carrara marble and observe that anisotropic porous domains develop spontaneously during shearing. Specifically, as samples are deformed periodic porous sheets are documented to emerge and are found to transfer mass. These results imply that viscous shear zones may naturally partition fluids into highly anisotropic bands. As this hydro-mechanical anisotropy is produced by creep, each porous sheet is interpreted to represent a transient dynamic pathway for fluid transport. It is unclear how long each porous domain is uniquely sustained but it is clear that sheets are persistently present with increasing strain. Our results forward the idea that viscous shear zones have dynamic transport properties that are not related to fracturing or chemical reaction. We believe these new results provide experimental foundation for changing the paradigm of viscosity in rocks to include dynamic permeability. In our view making this change in perspective could alter many classical interpretations in natural banded mylonite zones, for example shear zone parallel syn-kinematic veining may be the result of pore sheet instability and ductile fracturing.

Seismotectonic models that combine all the relevant seismotectonic data (e.g., hypocenter locations and velocity models, focal mechanisms and moment-tensors, faults, geodetic data, and in-situ/regional stress data) constitute a pre-requisite to better understand the interplay between stress, faulting and seismicity of a region. This study is a contribution to the multiannual project SeismoTeCH funded by the Swiss Geophysical Commission (SGPK) and coordinated by the University of Bern to produce an integrative seismotectonic model for the entire territory of Switzerland. In this context, our aim is to provide an up-to-date, high-quality, and consistent catalog of first-motion focal mechanisms computed by the Swiss Seismological Service (SED) since 1976. For this purpose, we developed a quality classification scheme for existing mechanisms based on a priori independent information (mainly applied to the oldest mechanisms in the catalog) combined with statistical methods such as HASH (Hardebeck and Shearer, 2002) and probabilistic source mechanisms scanner algorithms (Massin and Malcom, 2018) to probe the solution space and translate probability density functions to a discrete quality rating.Tests on selected problematic mechanisms are also carried out in order to assess the sensitivity of the focal mechanisms to the velocity models used to calculate location and take-off angles. In particular, we compare existing solutions using the standard 3D P-wave model of the SED with solutions based on recently derived high-resolution 3D Pg+Sg models. These tests are functional to understand the benefits of developing an updated full crustal velocity model for first-motion focal mechanisms calculations, in particular in relations to the focal depths and the accuracy of take-off angles.Finally, to improve the completeness of the existing catalog, we explore new methodologies that would incorporate automated (possibly real-time) and semi-automated techniques for expanding the calculation of first-motion focal mechanisms (and moment tensors) to events of smaller magnitude. The Anzere/Sanetschpass sequence of November 2019 is used to assess and develop these new methods. As a preliminary result of these combined efforts, we present here a high-resolution map of strain-based deformation across Switzerland. This work represents a benchmark for future regional-scale stress inversion and sets the basis for the development of a consistent, fully accessible, and dynamic focal mechanisms database for Switzerland.

The widespread presence of epidote-bearing veins and hydrous minerals such as micas in meta-granitoid rocks attests to the large extent of hydration of the exhuming continental crust. The ability of epidote (Ca2Al3Si3O12(OH) – Ca2Al2Fe3+Si3O12(OH)) to incorporate a wide variety of trace elements renders this mineral a promising geochemical tracer of circulating fluid(s).We report trace element and microstructural data on epidote-bearing veins from the Aar Massif (Central Alps) and from the Albula Pass (Eastern Alps). We characterized and classified the epidote-bearing veins based on their extent of deformation, shape and size of the epidote grains, coexisting minerals, and degree of dynamic recrystallization of associated quartz. Laser ablation ICP-MS data of individual epidote crystals reveal prominent zoning, confirmed by electron probe maps for Sr and Mn. Overall, low to very low Th/U ratios (down to 0.0005 in the Aar Massif veins and 0.001 in the Albula ones) with Th often below limits of detection (< 0.1 µg/g at 16 µm beam size) go along with variably LREE-depleted patterns (and CI Chondrite-normalized LaN/YbN ~0.35 in the Aar Massif veins and ~0.60 in the Albula Pass veins). Strontium contents are variable (hundreds to thousands of µg/g) and mostly high (up to 10100 µg/g in the Aar Massif samples and 12800 µg/g in the Albula Pass samples). The in-situ geochemical data are linked to the microstructures in order to assess whether microstructures can be related to variations in trace elements, also considering the role of coexisting phases. Moreover, trace element data of samples from the Aar Massif are compared to metamorphic host-rock epidotes and cleft epidotes from the same massif.We find that REE patterns of Aar Massif vein epidotes are clearly different than those of metamorphic host-rock epidotes and of cleft epidotes. In addition, REE patterns vary based on the microstructural characteristics of veins. Overall REE patterns of the Albula Pass vein epidotes resemble those from the Aar Massif. Different veins and microstructures define clusters in Sr vs. Y, Eu anomaly vs. Th/U ratios, and Eu anomaly vs. U values. Geochemical heterogeneities are observed among sampling localities within the Aar Massif.The fact that the geochemical characteristics of retrograde hydrothermal vein epidotes are clearly different than those of high-grade metamorphic and metamorphic host-rock epidotes, and the relationship between geochemical characteristics and microstructures support the hypothesis that the deformation did not alter the original geochemical record through neomineralization. Our data argue for the potential of epidote as a powerful fluid tracer in the granitoid continental crust.

Creep cavities are increasingly recognized as an important syn-kinematic feature of shear zones, but much about this porosity needs investigation. Largely, observations of creep cavities are restricted to very fine grained mature ultramylonites, and it is unclear when they developed during deformation. Specifically, a question that needs testing is should grain size reduction during deformation produce creep cavities? To this end, we have reanalyzed the microstructure of a large shear strain laboratory experiment that captures grain size change by dynamic recrystallization during mylonitization. We find that the experiment does contain creep cavities. Using a combination of scanning electron microscopy and spatial point statistics, we show that creep cavities emerge with, and because of, subgrain rotation recrystallization during ultramylonite formation. As dynamic recrystallization is ubiquitous in natural shear zones, this observation has important implications for the interpretation of concepts such as the Goetze criterion, paleopiezometery, and phase mixing.

Groundwater flow in granitic bedrock is of major interest for underground projects such as radioactive waste disposal. It is generally accepted that granitic rocks of the upper crust are characterized as faulted low-porosity rocks showing fault-related permeability. In this study, the influence of existing faults on the present-day water flow in the Grimsel Test Site (Switzerland), an underground rock laboratory situated in granitoid rocks, was investigated by mapping water discharges. As a result, the link between water flow and faults considering slip-tendency analysis and fault intersections is evaluated. Water-conducting features were combined in a structural-permeability favorability map. Faults and dykes occur as three orientation groups, NE–SW, E–W, and NW–SE trending, all steeply dipping southwards with fault intersections also steeply plunging southwards. In total, 100 water discharges were mapped in summer 2014 and 85 in winter 2016, which are located along faults or fault intersections. A comparison of water discharges with structures showed that high-slip-tendency metabasic dykes and fault or dyke intersections represent the dominant flow paths. Further, it could be demonstrated that higher slip-tendency tends to lead to enhanced average hydraulic conductivity and therefore more constant water flow. Based on water fluxes, fault intersections are inferred to represent first-order locations of water percolation followed by high-slip-tendency metabasic dykes. The combination of all water-conducting features into a structural-permeability map results in covering all water discharges. Therefore, the structural-permeability favorability map can serve as suitable representation for constraining water inflow in fractured granitoid host rocks.

International audience; Low-temperature chronometers offer potential to gain insights into the temporal evolution of hydrothermal systems. The long-lived fault-bound Grimsel pass hydrothermal system (including a fossil and an active part) in the Central European Alps serves here as a key site to test such an application. Zircon and apatite grains were separated from samples collected along a fault transect. The resulting zircon (U-Th)/He ages are homogenous along the profile at 8-9 Ma and thus record the regional cooling evolution, remaining unaffected by the younger hydrothermal activity. In contrast, the apatite (U-Th)/He ages show three age groups: One Group (1) of ca. 5 Ma inside and outside the hydrothermal zone matches the low-temperature part of the regional cooling trend, while group (2) with ages as young as 1-2 Ma occurs in a central narrow zone associated with hydrothermal activity. One sample (group 3) displays older apparent ages compared to the regional cooling trend. Group (2) apatite samples reveal a different cathodoluminescence texture and trace-element chemistry, which we interpret together with the young age as apatite growth or re-crystallization within the hydrothermal system. Forward 1D modelling of He diffusion indicates that apatite (U-Th)/He ages should always be reset when exposed to hot thermal waters (up to ~140 • C) present over ka timescales or to intermediate temperature waters (~90 • C) over Ma timescales. Combining our measured apatite (U-Th)/He ages with forward modelling results highlight that, besides regional cooling trends, local heat anomalies within hydrothermal zones are very variable in space and time. Combined trace-element geochemistry and (U-Th)/He dating shows local occurrence of newly-formed apatites crystals, which are best described as geochronometers rather than thermochronometers. Such information is important to explore the longevity of hydrothermal systems and associated spatial distributions of heat anomalies.