Search

Fluid-rock interactions play a critical role in Earth’s lithosphere and environmental subsurface systems. In the absence of chemical mass transport, mineral-hydration reactions would be accompanied by a solid-volume increase that may induce differential stresses and associated reaction-induced deformation processes, such as dilatant fracturing to increase fluid permeability. However, the magnitudes of stresses that manifest in natural systems remain poorly constrained. We used optical and electron microscopy to show that one of the simplest hydration reactions in nature [MgO + H2O = Mg(OH)2] can induce stresses of several hundred megapascals, with local stresses of as much as ∼1.5 GPa. We demonstrate that these stresses not only cause fracturing but also induce plastic deformation with dislocation densities (1015 m−2) exceeding those typical of tectonically deformed rocks. If these reaction-induced stresses can be transmitted across larger length scales, they may influence the bulk stress state of reacting regions. Moreover, the structural damage induced may be the first step toward catastrophic rock failure, triggering crustal seismicity.

We present results from a natural deformed shear zone in the Turon de Técouère massif of the French Pyrenees that directly addresses the processes involved in strain localization, a topic that has been investigated for the last 40 years by structural geologists. Paleopiezometry indicates that differential stresses are variable both spatially across the zone, and temporally during exhumation. We have, however, also calculated strain rate, which remains constant despite changes in stress. This result appears to be at odds with recent experimental deformation on monophase (olivine) rocks, which indicate that strain localization occurs dominantly as a result of constant stress. We hypothesize that in the Turon de Técouère massif—and many natural shear zones—strain localization occurs as a result of reactions, which decrease the grain size and promote the activation of grain size sensitive deformation mechanisms. From a tectonics perspective, this study indicates that the deformation rate in a particular plate boundary is relatively uniform. Stress, however, varies to accommodate this deformation. This viewpoint is consistent with deformation at a plate boundary, but it is not the typical way in which we interpret strain localization.

Fluid flow in crystalline rocks in the absence of fractures or ductile shear zones dominantly occurs by grain boundary diffusion, as it is faster than volume diffusion. It is, however, unclear how reactive fluid flow is guided through such pathways. We present a microstructural, mineral chemical, and thermodynamic analysis of a static fluid-driven reaction from dry granulite to ‘wet’ eclogite. Fluid infiltration resulted in re-equilibration at eclogite-facies conditions, indicating that the granulitic protolith was out of equilibrium, but unable to adjust to changing P–T conditions. The transformation occurred in three steps: (1) initial hydration along plagioclase grain boundaries, (2) complete breakdown of plagioclase and hydration along phase boundaries between plagioclase and garnet/clinopyroxene, and (3) re-equilibration of the rock to an eclogite-facies mineral assemblage. Thermodynamic modelling of local compositions reveals that this reaction sequence is proportional to the local decrease of the Gibbs free energy calculated for ‘dry’ and ‘wet’ cases. These energy differences result in increased net reaction rates and the reactions that result in the largest decrease of the Gibbs free energy occur first. In addition, these reactions result in a local volume decrease leading to porosity formation; i.e., pathways for new fluid to enter the reaction site thus controlling net fluid flow. Element transport to and from the reaction sites only occurs if it is energetically beneficial, and enough transport agent is available. Reactive fluid flow during static re-equilibration of nominally impermeable rocks is thus guided by differences in the energy budget of the local equilibrium domains.

Complex microstructures of zircon recrystallisation have been discovered in igneous and metamorphic zircons from the Jack Hills Metasedimentary Belt in Western Australia, and the granitic injection complex of Harris in north-west Scotland. Cathodoluminescence (CL), electron backscatter diffraction (EBSD) and wavelength dispersive spectroscopy (WDS) analysis reveal that recrystallisation in the studied zircons is often characterised by (1) translation, bending and fading of the oscillatory zoning; (2) structural recovery of the zircon; (3) an enrichment in Hf and depletion in Y in recrystallised zircon; and (4) the formation of a recrystallisation interface with a minor (

Melt transport mechanisms have an important impact on the chemical composition of the percolated host rock and the migrating melts. Melt migration is usually assumed to occur at grain boundaries. However, microstructural studies revealed the occurrence of polyphase inclusions along dislocations, subgrain boundaries and microcracks in single mineral grains. The inclusions are interpreted as crystallized melt pockets suggesting that melts can migrate within deformed crystals. Intracrystalline melt migration and diffusive re-equilibration can lead to significant mineral trace element enrichments when associated with dissolution–precipitation reactions. In this contribution, we study a body of replacive troctolites associated with the Erro-Tobbio ophiolitic mantle peridotites (Ligurian Alps, Italy). The replacive formation of the olivine-rich troctolite involved extensive impregnation of a dunitic matrix, i.e. partial dissolution of olivine and concomitant crystallization of interstitial phases. The olivine matrix is characterized by two distinct olivine textures: (i) coarse deformed olivine, representing relicts of the pre-existing mantle dunite matrix (olivine1), and (ii) fine-grained undeformed olivine, a product of the melt–rock interaction process (olivine2). Previous studies documented a decoupling between olivine texture and trace element composition, namely enriched trace element compositions in olivine1 rather than in olivine2, as would be expected from the dissolution–precipitation process. Notably, the trace element enrichments in deformed olivines are correlated with the occurrence of elongated 10 µm size polyphase inclusions (clinopyroxene, Ti-pargasite, chromite) preferentially oriented along olivine crystallographic axes. These inclusions show irregular contacts and have no crystallographic preferred orientation with the host olivine, and the phases composing the inclusions show similar chemical compositions to the vermicular phases formed at the grain boundaries during late-stage reactive crystallization of the troctolite. This suggests that the investigated inclusions did not form as exsolutions of the host olivine but rather by input of metasomatic fluids percolating through the deformed olivine grains during closure of the magmatic system. We infer that strongly fractionated volatile-rich melts were incorporated in oriented microfractures within olivine1 and led to the crystallization of the polyphase inclusions. The presence of intracrystalline melt greatly enhanced diffusive re-equilibration between the evolved melt and the percolated olivine1, in turn acquiring the enriched character expected in neoformed olivine crystals. Intracrystalline melt percolation can have strong geochemical implications and can lead to efficient re-equilibration of percolated minerals and rocks.

Seismological observations indicate strong variations in wave velocities around faults both co-seismically during earthquakes, and post-seismically. Recovery is commonly associated with a reduction in crack damage. Here, we explore the recovery associated with time-dependent mechanical closure of cracks. We report results from laboratory experiments conducted on dry cores of Carrara marble at room temperature. We deformed cylindrical samples in the semi-brittle regime to induce crack damage before subjecting them to hydrostatic and triaxial stress conditions for extended periods of time while recording dilatancy and wave speeds repeatedly. We report wave speed increases of up to 40% of the damage-induced wave speed drop in samples subject to hydrostatic loading. Moreover, we report the occurrence of significant wave speed increases contemporaneously with time-dependent creep in triaxially loaded samples. Wave speed recovery during creep is only observed below a threshold creep strain rate, a result we interpret as a transition from brittle to plastic creep with decreasing strain rate. We interpret the wave speed increase in terms of reduced crack density and increased contact area within the crack array, and show that around 40% of the total crack surface has to be closed to justify the observed wave speed recoveries. We propose that mechanical crack closure is driven by the viscous relaxation of the bulk rock under the influence of locked-in stresses at low confining pressure, and of the external stresses at higher confining pressure. Our study shows that mechanical crack closure is a significant source of time-dependent wave speed recovery.

The high-pressure metamorphic Nevado-Filábride Complex (NFC) in the Betics mountain range of southeastern Spain exhibits continental and ocean-derived tectonic units, which are key for understanding the geodynamic evolution of the Western Mediterranean. We address the current debate in the definition of tectonic units, the emplacement of (ultra)mafic rocks, and the timing of burial metamorphism by conducting a structural study combined with single grain fusion 40Ar/39Ar dating of white micas in structurally critical outcrops of the eastern Sierra de Los Filábres. One older 40Ar/39Ar age population (38–27 Ma) is found at distance from the main shear zones in the relics of an early foliation, while a younger 40Ar/39Ar population (22–12 Ma) is dominant in the vicinity of these shear zones, where the early foliation is obliterated. Both age groups are interpreted as the record of deformation or fluid-induced recrystallization during distinct fabric-forming events, while alternative scenarios are discussed. A key observation is the presence of an ophiolitic mélange, which—together with new and published geochronological data—allows for a new tectonic hypothesis. This considers Paleogene subduction beneath a Jurassic oceanic lithosphere, followed by the continued subduction of NFC and overlying ophiolites below the Alpujárride Complex. Exhumation during westward slab roll-back led to the formation of an extensional detachment system that obliquely cut nappe contacts. Although the timing constraints for high pressure-low temperature (HP-LT) metamorphism in the NFC remain inconclusive, the new tectonic hypothesis provides a solution that can account for both Paleogene and Miocene ages of HP-LT metamorphism.

The Multi-scale Laboratories (MSL) are a network of European laboratories bringing together the scientific fields of analogue modeling, paleomagnetism, experimental rock and melt physics, geochemistry and microscopy. MSL is one of nine (see below) Thematic Core Services (TCS) of the European Plate Observing System (EPOS). The overarching goal of EPOS is to establish a comprehensive multidisciplinary research platform for the Earth sciences in Europe. It aims at facilitating the integrated use of data, models, and facilities, from both existing and new distributed pan European Research Infrastructures, allowing open access and transparent use of data. The TCS MSL network allows researchers to collaborate with other labs and scientists. By becoming part of the rapidly growing TCS MSL network, new laboratories are offered a platform to showcase their research data output, laboratory equipment and information, and the opportunity to open laboratories to guest researchers through the Transnational Access (TNA) program. The EPOS Multi-scale laboratories offer researchers a fully operational data publication chain tailored to the specific needs of laboratory research, from a bespoke metadata editor, through dedicated, (domain-specific) data repositories, to the MSL Portal showcasing these citable data publications. During this process the data publications are assigned with digital object identidiers (DOI), published with open licenses (e.g. CC BY 4.0) and described with standardized and machine-readable rich metadata (following the FAIR Principles to make research data Findable, Accessible, Interoperable and Reusable. The TCS MSL is currently working on linking these data publications to the EPOS Central Portal1, the main discovery and access point for European multi-disciplinary data, and on increasing the number of connected data repositories. info:eu-repo/semantics/publishedVersion

SUMMARY Particle image velocimetry (PIV), a method based on image cross-correlation, is widely used for obtaining velocity fields from time-series of images of deforming objects. Rather than instantaneous velocities, we are interested in reconstructing cumulative deformation, and use PIV-derived incremental displacements for this purpose. Our focus is on analogue models of tectonic processes, which can accumulate large deformation. Importantly, PIV provides incremental displacements during analogue model evolution in a spatial reference (Eulerian) frame, without the need for explicit markers in a model. We integrate the displacements in a material reference (Lagrangian) frame, such that displacements can be integrated to track the spatial accumulative deformation field as a function of time. To describe cumulative, finite deformation, various strain tensors have been developed, and we discuss what strain measure best describes large shape changes, as standard infinitesimal strain tensors no longer apply for large deformation. PIV or comparable techniques have become a common method to determine strain in analogue models. However, the qualitative interpretation of observed strain has remained problematic for complex settings. Hence, PIV-derived displacements have not been fully exploited before, as methods to qualitatively characterize cumulative, large strain have been lacking. Notably, in tectonic settings, different types of deformation—extension, shortening, strike-slip—can be superimposed. We demonstrate that when shape changes are described in terms of Hencky strains, a logarithmic strain measure, finite deformation can be qualitatively described based on the relative magnitude of the two principal Hencky strains. Thereby, our method introduces a physically meaningful classification of large 2-D strains. We show that our strain type classification method allows for accurate mapping of tectonic structures in analogue models of lithospheric deformation, and complements visual inspection of fault geometries. Our method can easily discern complex strike-slip shear zones, thrust faults and extensional structures and its evolution in time. Our newly developed software to compute deformation is freely available and can be used to post-process incremental displacements from PIV or similar autocorrelation methods.

Seismic faulting causes wall rock damage driven by both mechanical stress and thermal energy. In the lower crust, coseismic damage has important implications for wall rock permeability, the progress of subsequent fluid-driven metamorphic reactions, and rock rheology. Wall rock microstructures reveal high-stress conditions near the slip surface during lower crustal earthquakes, however, there is limited documentation on the thermal effect. Here, we present a transmission electron microscopy study of coseismic microfractures in plagioclase feldspar from lower crustal granulites from the Bergen Arcs, Western Norway. Focused ion beam foils are collected 1.25 mm and 1.8 cm from a 2 mm thick eclogite facies pseudotachylyte vein. Dislocation-free plagioclase aggregates fill the microfractures and record a history of recovery from a short-lived high stress-temperature (σ-T) state caused by seismic slip and frictional melting along the nearby fault surface. The plagioclase aggregates retain the crystallographic orientation of the host rock and shape preferred orientation relative to the fault slip surface. We propose that plagioclase partially amorphized along the microfractures at peak stress conditions followed by repolymerization to form dislocation-free grain aggregates within the timeframe of pseudotachylyte formation. The heat from the slip surface dissipated into the wall rock causing a short-lived temperature peak. Subsequent cooling led to exsolution of intermediate plagioclase compositions by spinodal decomposition within a few millimeters distance to the fault surface. Our findings provide microstructural evidence for the high σ-T conditions that are expected in the proximity of seismic faults, highlighting the importance of micro- and nanostructures for the understanding of earthquakes ruptures.

The strength of the lithosphere is typically modelled based on constitutive equations for steady-state flow. However, strain hardening may cause significant evolution of strength in the colder load-bearing portion of the lithosphere. Recent rheological data from low-temperature deformation experiments on olivine suggest that strain hardening occurs due to the presence of temperature-independent back stresses generated by long-range elastic interactions among dislocations. These interpretations provided the basis for a flow law that incorporates hardening by the development of back stress. Here, we test this dislocation-interaction hypothesis by examining the microstructures of olivine samples deformed plastically at room temperature either in a deformation-DIA apparatus at differential stresses of ≤ 4.3 GPa or in a nanoindenter at applied contact stresses of ≥ 10.2 GPa . High-angular resolution electron backscatter diffraction maps reveal the presence of geometrically necessary dislocations with densities commonly above 10 14 m − 2 and intragranular heterogeneities in residual stress on the order of 1 GPa in both sets of samples. Scanning transmission electron micrographs reveal straight dislocations aligned in slip bands and interacting with dislocations of other types that act as obstacles. The resulting accumulations of dislocations in their slip planes, and associated stress heterogeneities, are consistent with strain hardening resulting from long-range back-stresses acting among dislocations and thereby support the form of the flow law for low-temperature plasticity. Based on these observations, we predict that back stresses among dislocations will impart significant mechanical anisotropy to deformed lithosphere by enhancing or reducing the effective stress. Therefore, strain history, with associated microstructural and micromechanical evolution, is an important consideration for models of lithospheric strength. The microstructural observations also provide new criteria for identifying the operation of back-stress induced strain hardening in natural samples and therefore provide a means to test the applicability of the flow law for low-temperature plasticity.

The Cap de Creus peninsula in NE Spain exposes Variscan amphibolite-facies rocks, transected by greenschist-facies mylonitic shear zones accommodating displacements of up to 1 km. One of these shear zones contains pseudotachylytes with 1–4 mm thick fault veins and cm-scale high-angle injection veins and breccias. The pseudotachylyte matrix encloses mm-to μm-scale wallrock fragments and consists of ultrafine biotite and feldspar, often with a microporphyritic structure suggesting crystallization from a melt. Recent dating of the host shear zone yielded a Jurassic (170-160 Ma) age. Here we present 40Ar/39Ar dating of the pseudotachylytes yielding Early Eocene ages between 52.76 ± 1.64 Ma and 44.55 ± 0.77 Ma, indicating Alpine brittle faulting at times that Jurassic low-angle ductile shears had been tilted to their present vertical position following Alpine thrusting in the Pyrenees. The kinematics of the pseudotachylytes suggest strike-slip reactivation along the host shear zone, consistent with plate-kinematic analyses indicating coeval Iberia - Europe motion dominated by strike-slip. The dimensions of the pseudotachylyte-bearing zone and estimated brittle displacements suggest a moment magnitude of ~4.7–4.9. The average thicknesses of the fault veins indicate deformation at ≥4 km depth, an energy density of ~9.106 J/m2 and a frictional resistance of ~23 MPa.

Shale host rock and containment potential are largely determined by the connected pore network in the rock, and the connection between the pore network and the naturally present or mechanically induced fracture network together determines the total bulk permeability. Pore connectivity in shales is poorly understood because most of the porosity is present in sub-micrometer-sized pores that are connected through nanometer-sized pore throats. We have used a number of different techniques to investigate the microstructure and permeability of Early Jurassic shales from the UK (Whitby Mudstone), under intact and fractured conditions. Whitby Mudstone is a clay matrix-rich rock (50-70%), with different mineralogical layers on the sub-millimeter scale and very low natural permeability (10-19 to 10-22 m2), representative of many gas shales and caprocks present in Europe. Artificial fracturing of this shale increases its permeability by 2-5 orders of magnitude at low confining pressure (5 MPa). At high confining pressures (30 MPa), permeability changes were more sensitive to the measuring direction with respect to the bedding orientation. Given the distinct lack of well-defined damage zones, most of the permeability increase is controlled by fracture permeability, which is sensitive to the coupled hydro-chemo-mechanical response of the fractures to fluids.

The ophiolitic peridotite and gabbro of Moncuni (Southern Lanzo Massif, Western Alps) retain pre-subduction mantle-to-oceanic, high-temperature ( > 700 C ∘ ) ductile fabrics. These fabrics are overprinted by seismic fracturing and faulting associated with pseudotachylytes. Within the gabbro, the pseudotachylytes preserve dry glass and pristine microlites. The occurrence of rare, minute garnet and the static development of eclogite-facies assemblages in local hydrated domains indicate that pseudotachylytes experienced subduction conditions of 600 C ∘ and 2.1 GPa. The exceptional survival of glass and the absence of post-oceanic ductile deformation demonstrate prevailing dry conditions during the entire Alpine subduction and exhumation path. Dry conditions inhibited reaction kinetics and viscous flow. In contrast, the majority of the Alpine ophiolites, derived from the upper hydrated portions of the oceanic lithosphere, show pervasive fluid-assisted metamorphism and ductile deformation. The Moncuni body can, therefore, be regarded as representative for the rheological behaviour during subduction of seismic, dry, deeper oceanic lithosphere that is rarely exhumed to the Earth's surface. In Moncuni, the brittle-ductile transition of dry oceanic rocks is constrained to be between 600 and 750 C ∘ . This temperature range corresponds to the observed cut-off of intermediate-depth seismicity within subducting slabs. We infer that the base of the seismic layer corresponds to the brittle-ductile transition of a dry slab rather than the locus of antigorite breakdown triggering earthquakes by dehydration embrittlement.

Borosilicate glasses are currently used for the immobilization of highly radioactive waste and are materials of choice for many biomedical and research industries. They are metastable materials that corrode in aqueous solutions, reflected by the formation of silica-rich surface alteration layers (SAL). Until now, there is no consensus in the scientific community about the reaction and transport mechanism(s) and the rate-limiting steps involved in the formation of SALs. Here we report the results of multi-isotope tracer (2H,18O,10B, 30Si, 44Ca) corrosion experiments that were performed with precorroded and pristine glass monoliths prepared from the six-component international simple glass and a quaternary aluminum borosilicate glass. Results of transmission electron microscopy and nanoscale analyses by secondary ion mass spectrometry reveal a nanometer-sharp interface between the SAL and the glass, where decoupling of isotope tracer occurs, while proton diffusion and ion exchange can be observed within the glass. We propose a unifying mechanistic model that accounts for all critical observations so far made on naturally and experimentally corroded glasses. It is based on an interface-coupled glass dissolution-silica precipitation reaction as the main SAL forming process. However, a diffusion-controlled ion exchange front may evolve in the glass ahead of the dissolution front if SAL formation at the reaction interface significantly slows down due to transport limitations. A unifying mechanistic model has been developed for silicate glass corrosion that can explain all critical observations made to-date. Borosilicate glasses are often used in biomedical devices and to store and dispose of radioactive waste. They decay in aqueous solution via the generation of a porous ‘surface alteration layer’ (SAL), the structure of which is different to the bulk. How the SAL forms is still not clear, however, an international team lead by Thorsten Geislern at the University of Bonn, Germany, has now used multi-isotope tracer experiments, to provide detailed insight into the distinct chemical and transport steps occurring during SAL formation. Their results suggest that an ‘interface-coupled dissolution-precipitation’ reaction is the main mechanism at play during SAL formation, but that, if slowed by transport limitations, it may be replaced by an ‘interdiffusion’ process.

Serpentinization may provide a unique environment for the abiotic formation of condensed carbonaceous matter. This could support the deep biosphere and contribute to the deep carbon cycle, and may have provided the first building blocks for life. However, thus far, condensed carbonaceous matter has been found only in association with the minor mineral constituents of serpentinites. In contrast, here we show the direct association between carbonaceous matter and the dominant Fe oxide in serpentinites, magnetite. Our samples were recovered from the Yap Trench, western Pacific Ocean, with a human-occupied vehicle at a depth of 6413 m below sea level. The carbonaceous matter coincides with some micron-sized magnetite grains, but particularly with nanosized Fe oxides within serpentinite nanopores. Vibrational spectroscopy reveals that the condensed carbonaceous matter contains both aliphatic and aromatic compounds, but there is no evidence for functional groups typical for biological organics. Based on these observations, we suggest that physicochemical phenomena in serpentinite nanopores and nanosized catalytically active minerals may play a key role in the abiotic synthesis of complex carbonaceous matter.

Amino acids and polycyclic aromatic hydrocarbons (PAHs) belong to the range of organic compounds detected in meteorites. In this study, we tested empirically and theoretically if PAHs are precursors for amino acids in carbonaceous chondrites, as previously suggested. We conducted experiments to synthesize amino acids from fluoranthene (PAH), with ammonium bicarbonate as a source for ammonia and carbon dioxide under mimicked asteroidal conditions. In our thermodynamic calculations, we extended our analysis to additional PAH-amino acid combinations. We explored 36 reactions involving the PAHs naphthalene, anthracene, fluoranthene, pyrene, triphenylene, and coronene and the amino acids glycine, alanine, valine, leucine, phenylalanine, and tyrosine. Our experiments do not show the formation of amino acids, whereas our theoretical results hint that PAHs could be precursors of amino acids in carbonaceous chondrites at low temperatures.

Nanoplatelets offer many possibilities to construct advanced materials due to new properties associated with their (semi)two-dimensional shapes. However, precise control of both positional and orientational order of the nanoplatelets in three dimensions, which is required to achieve emerging and collective properties, is challenging to realize. Here, we combine experiments, advanced electron tomography and computer simulations to explore the structure of supraparticles self-assembled from nanoplatelets in slowly drying emulsion droplets. We demonstrate that the rich phase behaviour of nanoplatelets, and its sensitivity to subtle changes in shape and interaction potential can be used to guide the self-assembly into a wide range of different structures, offering precise control over both orientation and position order of the nanoplatelets. Our research is expected to shed light on the design of hierarchically structured metamaterials with distinct shape- and orientation- dependent properties. Nanoplatelets can be used as anisotropic building blocks for constructing novel optoelectronic materials. Here, Wang et al. show a route of assembling nanoplatelets with controllable positional and orientational order in three dimensions facilitated by the surface tension of drying emulsion droplets.

In this study, we examine the effect of transient mantle creep on the prediction of glacial isostatic adjustment (GIA) signals. Specifically, we compare predictions of relative sea level (RSL) change from GIA from a set of Earth models in which transient creep parameters are varied in a simple Burgers model to a reference case with a Maxwell viscoelastic rheology. The model predictions are evaluated in two ways: first, relative to each other to quantify the effect of parameter variation, and second, for their ability to reproduce well-constrained sea level records from selected locations. Both the resolution and geographic location of the RSL observations determine whether the data can distinguish between model cases. Model predictions are most sensitive to the inclusion of transient mantle deformation in regions that are near-field and peripheral relative to former ice sheets. This sensitivity appears particularly true along the North American west coast in the region of the former Cordilleran Ice Sheet, which experienced rapid sea-level fall following deglaciation between 14 and 12 kyr BP. Relative to the Maxwell case, Burgers models better reproduce this rapid phase of regional postglacial sea-level fall. As well, computed goodness-of-fit values in this region show a clear preference for models where transient deformation is present in the whole or lower mantle, and for models where the rigidity of the Kelvin element is weakened relative to the rigidity of the Maxwell element. In contrast, model predictions of relative sea-level change in the far-field show weak sensitivity to the inclusion of transient deformation.