Modeling episodic fluid-release events in the ductile carbonates of the Glarus thrust

The exposed Glarus thrust displays midcrustal deformation with tens of kilometers of displacement on an ultrathin layer, the principal slip zone (PSZ). Geological observations indicate that this structure resulted from repeated stick-slip events in the presence of highly overpressured fluids. Here we show that the major characteristics of the Glarus thrust movement (localization, periodicity, and evidence of pressurized fluids) can be reconciled by the coupling of two processes, namely, shear heating and fluid release by carbonate decomposition. During this coupling, slow ductile creep deformation raises the temperature through shear heating and ultimately activates the chemical decomposition of carbonates. The subsequent release of highly overpressurized fluids forms and lubricates the PSZ, allowing a ductile fault to move tens of kilometers on millimeter-thick bands in episodic stick-slip events. This model identifies carbonate decomposition as a key process for motion on the Glarus thrust and explains the source of overpressured fluids accessing the PSZ.

Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: reconstruction of differential stress magnitudes and pre-deformation fluid properties

A well developed theoretical framework is available in which paleofluid properties, such as chemical composition and density, can be reconstructed from fluid inclusions in minerals that have undergone no ductile deformation. The present study extends this framework to encompass fluid inclusions hosted by quartz that has undergone weak ductile deformation following fluid entrapment. Recent experiments have shown that such deformation causes inclusions to become dismembered into clusters of irregularly shaped relict inclusions surrounded by planar arrays of tiny, new-formed (neonate) inclusions. Comparison of the experimental samples with a naturally sheared quartz vein from Grimsel Pass, Aar Massif, Central Alps, Switzerland, reveals striking similarities. This strong concordance justifies applying the experimentally derived rules of fluid inclusion behaviour to nature. Thus, planar arrays of dismembered inclusions defining cleavage planes in quartz may be taken as diagnostic of small amounts of intracrystalline strain. Deformed inclusions preserve their pre-deformation concentration ratios of gases to electrolytes, but their H2O contents typically have changed. Morphologically intact inclusions, in contrast, preserve the pre-deformation composition and density of their originally trapped fluid. The orientation of the maximum principal compressive stress (σ1σ1) at the time of shear deformation can be derived from the pole to the cleavage plane within which the dismembered inclusions are aligned. Finally, the density of neonate inclusions is commensurate with the pressure value of σ1σ1 at the temperature and time of deformation. This last rule offers a means to estimate magnitudes of shear stresses from fluid inclusion studies. Application of this new paleopiezometer approach to the Grimsel vein yields a differential stress (σ1–σ3σ1–σ3) of ∼300 MPa∼300 MPa at View the MathML source390±30°C during late Miocene NNW–SSE orogenic shortening and regional uplift of the Aar Massif. This differential stress resulted in strain-hardening of the quartz at very low total strain (<5%<5%) while nearby shear zones were accommodating significant displacements. Further implementation of these experimentally derived rules should provide new insight into processes of fluid–rock interaction in the ductile regime within the Earth's crust.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

The tectonometamorphic evolution of the Sesia–Dent Blanche nappes (internal Western Alps): review and synthesis

This study reviews and synthesizes the present knowledge on the Sesia–Dent Blanche nappes, the highest tectonic elements in the Western Alps (Switzerland and Italy), which comprise pieces of pre-Alpine basement and Mesozoic cover. All of the available data are integrated in a crustal-scale kinematic model with the aim to reconstruct the Alpine tectono-metamorphic evolution of the Sesia–Dent Blanche nappes. Although major uncertainties remain in the pre-Alpine geometry, the basement and cover sequences of the Sesia–Dent Blanche nappes are seen as part of a thinned continental crust derived from the Adriatic margin. The earliest stages of the Alpine evolution are interpreted as recording late Cretaceous subduction of the Adria-derived Sesia–Dent Blanche nappes below the South-Alpine domain. During this subduction, several sheets of crustal material were stacked and separated by shear zones that rework remnants of their Mesozoic cover. The recently described Roisan-Cignana Shear Zone of the Dent Blanche Tectonic System represents such a shear zone, indicating that the Sesia–Dent Blanche nappes represent a stack of several individual nappes. During the subsequent subduction of the Piemonte–Liguria Ocean large-scale folding of the nappe stack (including the Roisan-Cignana Shear Zone) took place under greenschist facies conditions, which indicates partial exhumation of the Dent Blanche Tectonic System. The entrance of the Briançonnais micro-continent within the subduction zone led to a drastic change in the deformation pattern of the Alpine belt, with rapid exhumation of the eclogite-facies ophiolite bearing units and thrust propagation towards the foreland. Slab breakoff probably was responsible for allowing partial melting in the mantle and Oligocene intrusions into the most internal parts of the Sesia–Dent Blanche nappes. Finally, indentation of the Adriatic plate into the orogenic wedge resulted in the formation of the Vanzone back-fold, which marks the end of the pervasive ductile deformation within the Sesia–Dent Blanche nappes during the earliest Miocene.

Valle d’Aosta section of the Sesia Zone: multi-stage HP metamorphism and assembly of a rifted continental margin

The spectrum characteristic of the EMC ranges from eclogites (containing omphacite and/or jadeite, garnet, phengite, glaucophane, zoisite, chloritoid, rutile) to phengite schists, calcschists, and marbles, as well as a variety of orthogneisses. Despite the intense polyphase deformation and HP-metamorphic recrystallization, it is possible in some locations to recognize pre-Alpine characteristics in some of the protoliths. For instance, two types of felsic orthogneiss can be distinguished in the Aosta Valley, one derived from Permian granitoids (with local preservation of intrusive contacts, magmatic inclusions, leucocratic veins and other magmatic structures; Stop 3), the other derived from pre-Variscan leuco-monzogranite, such as the building stone mined at the “Argentera” quarry near Settimo Vittone / Montestrutto (Stop 2; so-called “Verde Argento” contains jadeite, phengite, K-feldspar, quartz). Polycyclic and more rarely monocyclic metasediments contain evidence of a complex Alpine PTDt-evolution, locally including relics of their prograde history from blueschist, one or more stages at eclogite facies. Recent petrochronological studies have dated this HP-evolution of the Sesia Zone in some detail. In the area visited, clear evidence of HP-cycling has been identified in one km-size tectonic slice (Stop 1), but not in adjacent parts of the EMC, indicating “yo-yo tectonics”. Partial retrogression and attendant ductile to brittle deformation of the HP-rocks is evident in one of the outcrops (Stop 4). Apart from the four localities in the Sesia Zone, a final outcrop introduces HP-rocks of the adjacent Piemonte oceanic unit, specifically calc-schists and ophiolite members of the “Zermatt-Saas” zone. The hilltop outcrop (Stop 5) displays foliated antigorite schist with peridotite relics (clinopyroxene, spinel) containing lenses derived from doleritic dykes. These fine-grained metarodingites and the folded veins containing Mg-chlorite and titanoclinohumite within serpentinite once again indicate equilibration under low-temperature eclogite facies conditions. However, these units reached that HP stage more than 20 Ma after the youngest eclogite facies imprint recognized in the Sesia Zone. Despite nearly half a century of intense study in the Sesia Zone, the complex assembly of its HP-terranes and their relation to more external parts of the Western Alps remains incompletely understood. This field guide merely introduces a few of the classic outcrops and discusses some of the critical evidence they contain, but it could not incorporate details on each stage of the evolution recognized so far.

Inherited structural controls on fault geometry, architecture and hydrothermal activity: an example from Grimsel Pass, Switzerland

Exhumed faults hosting hydrothermal systems provide direct insight into relationships between faulting and fluid flow, which in turn are valuable for making hydrogeological predictions in blind settings. The Grimsel Breccia Fault (Aar massif, Central Swiss Alps) is a late Neogene, exhumed dextral strike-slip fault with a maximum displacement of 25–45 m, and is associated with both fossil and active hydrothermal circulation. We mapped the fault system and modelled it in three dimensions, using the distinctive hydrothermal mineralisation as well as active thermal fluid discharge (the highest elevation documented in the Alps) to reveal the structural controls on fluid pathway extent and morphology. With progressive uplift and cooling, brittle deformation inherited the mylonitic shear zone network at Grimsel Pass; preconditioning fault geometry into segmented brittle reactivations of ductile shear zones and brittle inter-shear zone linkages. We describe ‘pipe’-like, vertically oriented fluid pathways: (1) within brittle fault linkage zones and (2) through alongstrike- restricted segments of formerly ductile shear zones reactivated by brittle deformation. In both cases, low-permeability mylonitic shear zones that escaped brittle reactivation provide important hydraulic seals. These observations show that fluid flow along brittle fault planes is not planar, but rather highly channelised into sub-vertical flow domains, with important implications for the exploration and exploitation of geothermal energy.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Tectono-metamorphic evolution of a nappe stack: A case study of the Swiss Alps

Fold-and-thrust belts are prominent structures that occur at the front of compressional orogens. To unravel the tectonic and metamorphic evolution of such complexes, kinematic investigations, quantitative microstructural analysis and geothermometry (calcite–graphite, calcite–dolomite) were performed on carbonate mylonites from thrust faults of the Helvetic nappe stack in Central Switzerland. Paleo-isotherms of peak temperature conditions and cooling stages (fission track) of the nappe pile were reconstructed in a vertical section and linked with the microstructural and kinematic evolution. Mylonitic microstructures suggest that under metamorphic conditions close to peak temperature, strain was highly localized within thrust faults where deformation temperatures spatially continuously increased in both directions, from N to S within each nappe and from top–down in the nappe stack, covering a temperature range of 180–380 °C. Due to the higher metamorphic conditions, thrusting of the lowermost nappe, the Doldenhorn nappe, was accompanied by a much more pronounced nappe internal ductile deformation of carbonaceous rock types than was the case for the overlying Wildhorn- and Gellihorn nappes. Ongoing thrusting brought the Doldenhorn nappe closer to the surface. The associated cooling resulted in a freezing in of the paleo-isotherms of peak metamorphic conditions. Contemporaneous shearing localized in the basal thrust, initially still in the ductile deformation regime and finally as brittle faulting and cataclasis inducing ultimately an inverse metamorphic zonation. With ongoing exhumation and the formation of the Helvetic antiformal nappe stack, a bending of large-scale tectonic structures (thrusts, folds), peak temperature isotherms and cooling isotherms occurred. While this local bending can directly be attributed to active deformation underneath the section investigated up to times of 2–3 ma, a more homogeneous uplift of the entire region is suggested for the very late and still active exhumation stage.