Fluids released from subducting hydrated rocks influence volcanism, tectonics, and geochemical cycling, but the mechanisms of fluid escape in subduction zones remain poorly understood. We address this issue by investigating the Erro-Tobbio meta-serpentinites (ET-MS), Italy, exhumed serpentinite rocks that preserve extensive dehydration vein networks formed by the porosity-generating breakdown of antigorite and brucite. We characterized the structure and morphology of these self-organized vein networks and evaluated their hydrodynamic properties using a novel approach. Specifically, we combined X-ray tomography and drone imagery with generative machine learning, electron microscopy, and equilibrium thermodynamics to model and analyze fluid pathways in the ET-MS. In both natural and simulated samples, these dehydration vein networks act as efficient drainage systems, enabling rapid fluid percolation even at porosities below 1%. The maximum network permeability is , several orders of magnitude higher than that of intact serpentinite. Fe-rich olivine and monticellite occur alongside relict brucite and magnetite in these veins. This assemblage indicates that the high permeability arises from porosity localized along brucite- and magnetite-rich veins, where infiltration of reducing fluids enhanced dehydration reactions. These findings demonstrate that serpentinite dehydration in subduction zones can produce flow-optimized vein structures that efficiently channel fluids at low porosity, potentially influencing fluid migration on local to regional scales before widespread dehydration occurs.

Aqueous fluids released by metamorphic dehydration of serpentinites are a key component for seismicity, creep, and geochemical cycling in subduction zones. How these fluids drain and migrate towards the mantle wedge has yet to be fully understood. Here we address the influence of pre-existing structural and mineralogical heterogeneities in serpentinites on dehydration and fluid migration at forearc conditions. We partially dehydrated natural serpentinite containing brucite veins in a piston-cylinder apparatus with a temperature gradient across the conditions of the brucite + antigorite = olivine + fluid reaction (485–520 °C; 1.5 GPa). Micro-tomography, electron microscopy and microstructural analysis of the experimental results, coupled with thermodynamic modelling, show that temperature, mineralogical heterogeneity and variable ingress of external H 2 controlled the dehydration extent. Experimentally formed olivine indicates a topotactic relationship between [100] Ol and [0001] Brc , although the resultant fabric is overall random because brucite was randomly oriented. Olivine forms mono-mineralic aggregates along the walls of brucite veins, displaying very high porosity (up to 32%) and permeability (10 –13 –10 –14 m 2 ). Tracing the pre-existing brucite vein network, these aggregates can form a transient network of interconnected, highly permeable fluid channels that allows drainage and may enhance open-system exchange with neighboring lithologies. Infiltration of reduced external fluids can trigger redox dehydration of magnetite + antigorite to Fe-rich olivine, which renews porosity and propagates focused fluid flow. The distribution of brucite and magnetite, especially as vein networks, therefore has a first-order control on how focused fluid drainage and flow paths develop during subduction of serpentinites.

The dehydration of antigorite is an important reaction in subduction zones with implications on both geochemical and geophysical processes. In this experimental study we focus on the onset of antigorite dehydration and investigate various chemical and physical parameters as possible drivers for the fluid release. We performed hydrostatic and co-axial Griggs experiments on antigorite serpentinites with variable chemical composition and microstructures at high-pressure and high-temperature conditions across the antigorite dehydration (1.5 GPa, 620–670 °C). For these conditions, our thermodynamic models predict the formation of olivine from magnetite decomposition and partial dehydration of antigorite. Detailed analyses of the run products reveal limited magnetite decomposition. Antigorite dehydration is restricted to samples that have been deformed. Nano-sized olivine and orthopyroxene formed locally in oblique dehydration bands and exhibit neither a clear crystallographic preferred orientation nor a topotactic relation with precursor antigorite. We argue that limited local dehydration in our experiments is related to strain and variations in reaction kinetics. Systematic investigation excludes mineralogical and chemical heterogeneities, and temperature gradients as reaction driving potentials. The structural relation of the dehydration bands suggests deformation-related dehydration, which is supported by numerical simulations that couple reaction kinetics with mechanical work rate and self-consistently predict dehydration bands. In this scenario, strain concentration due to applied axial stress locally increases the internal energy of antigorite to reach the activation energy of the dehydration reaction, enabling dehydration. This study highlights the importance of coupled mechanical and chemical processes and provides a mechanistic framework for deformation-induced dehydration of antigorite.

Aqueous fluids released by metamorphic dehydration reactions are key components for magmatism, seismicity, creep, and geochemical cycling in subduction zones. How these fluids drain and migrate towards the mantle wedge is not fully understood, partly because the recognition and interpretation of deep fluid pathways in the exhumed rock record is challenging. Serpentinites are among the most important H2O carriers in subducted slabs, with dehydration occurring for example during the lizardite to antigorite transition (ca. 300 – 350 °C), brucite breakdown (ca. 500 °C) and antigorite dehydration (620 – 670 °C). Fluid production and the related formation of interconnected porosity allowing fluid migration are influenced by pre-existing chemical and mineralogical heterogeneities, as well as microstructure and porosity. In oceanic and low-grade metamorphic serpentinites such heterogeneities are very common. To investigate their effect on metamorphic dehydration and fluid migration during subduction to forearc conditions, we experimentally dehydrated natural serpentinite that contains abundant brucite formed during the prograde lizardite–antigorite transformation [1]. In the starting material, brucite occurs as veins and as intergrowths with serpentine in the matrix. We performed piston-cylinder experiments at conditions of brucite dehydration in subduction zones (520 – 570 °C; 1.5 GPa), coupled with micro-tomography (µ-CT), Raman, electron microscopy and microstructural analysis. The experimental results show the formation of olivine as (i) veinlets along the rims of brucite veins, (ii) surrounding and replacing Fe-oxides, and as (iii) tabular grains growing in the serpentinite matrix at the hot spot of the sample cylinder. All olivine types are related to newly formed porosity visible at the resolution of the µ-CT (1.2 µm voxel size). Broad-ion beam polished FE-SEM analysis of the starting material indicates that veins of brucite (± serpentine, Fe-oxides) have significantly more nano-porosity than the serpentine matrix. This observation and the formation of olivine veinlets along the previous brucite vein walls in the experiment suggest that the presence of brucite veins –formed early during shallow forearc metamorphism of serpentinite– will influence fluid production and migration pathways during brucite and antigorite dehydration at deep forearc conditions. Incidentally, our results further demonstrate that preferential dehydration occurs when an external reducing agent (in the case of the experiment H2 most likely derived from the graphite heater) triggers the replacement of serpentine + magnetite by olivine, in line with previous experimental and natural observations [2,3]. Incipient fluid release from serpentinite is thus heterogeneous at the microscale, and will cause local fluid pressure variations that may lead to flow and ultimately drainage. This process, possibly in combination with deformation, deviatoric stress and/or external fluid flux, may favour the development of commonly mono-mineralic olivine veins, which are inferred to form by an interplay of brucite dehydration and reactive fluid flow [4,5]. [1] Menzel et al., 2018, Lithos[2] Eberhard et al., 2023, Journal of Petrology[3] Padrón-Navarta et al., 2023, Nature Geoscience[4] Plümper et al., 2017, Nature Geoscience[5] Huber et al., 2022, G-cubed Funding: M.D.M: Junta de Andalucía (Postdoc_21_00791) and project “RUSTED”, MCIU Spain (PID2022-136471N-B-C21 & 22). LE: NWO (VI.Vidi.193.030), EXCITE (TNA-C3-2023-13).

To assess the seismogenic potential of fault zones it is crucial to understand fluid-rock interactions in these zones, because alteration affects the fault strength and stability, as well as the deformation mechanisms.The San Andreas fault (SAF) system is known for infrequent large magnitude (M≥7) earthquakes, whereas some segments lack such strong seismic events [1]. Here, strain is largely accommodated by creep motion. Aseismic creep can be enhanced by the presence of fluids, which may additionally drive mineral reactions. For example, fluid composition and magnesite deposits in the SAF segment between San Juan Bautista and Parkfield suggests carbonation due to infiltration of CO2-bearing fluids into the fault [2]. Carbonation of ultramafic rocks leads to the formation of talc, which is known to be frictionally weak and promotes creep when wet [3]. However, our thermodynamic fluid-infiltration calculations show that carbonation will not produce pure talc but lizardite-talc-magnesite (LTM) and talc-magnesite rocks (soapstone) and, with increasing extent of reactive fluid flow, talc-magnesite-quartz (TMQ) and magnesite-quartz rocks (listvenite). The strength and seismogenic potential of serpentinite fault zones undergoing carbonation thus may change dynamically as the mineral proportions and assemblages change, but the respective frictional behaviour of these assemblages is unknown.We performed rotary-shear experiments on gouge layers with compositions ranging from lizardite-serpentinite to LTM, soapstone, TMQ and listvenite at pressure, temperature and pore fluid pressures corresponding to a depth of about 10 km (300 °C, 250 MPa normal stress and 100 MPa pore pressure). We measured the frictional strength within the velocity range of 0.002 µm/s to 10 µm/s.Our data show that lizardite gouges are relatively strong and slightly velocity-weakening. The friction coefficient dropped from 0.45 at 0.002 µm/s to 0.42 at 10 µm/s. A similar velocity-dependence is observed for soapstone gouges, although at lower absolute friction coefficients of 0.3 to 0.28. Interestingly, listvenite gouges show the opposite behavior, with friction coefficients increasing from 0.25 at 0.002 µm/s to 0.48 at 10 µm/s. Stick-slips were only observed in serpentinite and soapstone gouges at low velocities. Increasing velocities and progressing carbonation causes stable slip behavior. Microtextural observations indicate strong grain-size reduction and basal cleavage in serpentinite gouges. On the contrary, soapstone and listvenite gouges show a fine-grained magnesite matrix surrounding the silicates.Our results suggest that serpentinized fault zones have the potential to nucleate unstable slip. The results further confirm the strong weakening effect of carbonation. CO2-fluid-rock interaction in ultramafic fault gouges may effectively suppress the nucleation of earthquakes. Since also listvenite gouges deformed aseismic and are found to be frictionally weak at low velocities, we suggest that besides talc also magnesite plays an important role in the deformation behavior of carbonated ultramafic fault zones. [1] Jolivet et al. 2015. Geophys. Res. Lett. doi:10.1002/2014GL062222.[2] Klein et al. 2022. Geophys. Res. Lett. doi:10.1029/2022GL099185.[3] Moore et al. 2008. Tectonophysics. doi:10.1016/j.tecto.2007.11.039 FundingLE: NWO (VI.Vidi.193.030)M.D.M: Junta de Andalucía (Postdoc_21_00791) and MCIU, Spain (PID2022-136471N-B-C22)

Dehydration reactions play a pivotal role in the dynamics and seismicity at subduction zones and in the deep water cycle. These reactions often occur during rock deformation. The dehydration of antigorite serpentinite is particularly important at subduction zones. This dehydration has been investigated with laboratory experiments of serpentinite deformation. Yet, the reproduction of such laboratory deformation experiments of serpentinite dehydration with mathematical models is still a major challenge. Here, we test a two-dimensional (2D) hydro-mechanical-chemical (HMC) numerical model for serpentinite dehydration by comparing the numerical results with the results of laboratory experiments.The laboratory experiments are performed with a Griggs apparatus. Natural antigorite serpentinites with and without preferred orientation are deformed by vertical compression for a confining pressure of 1.5 GPa and maximum differential stresses between 350 and 700 MPa. For comparison, also experiments with hydrostatic stress are performed. The experimental temperature is between 620 and 650 °C. The applied confining pressure and temperature are in the olivine stability field according to the measured chemical composition of the serpentinite and thermodynamic Perple_X calculations. However, olivine only forms locally in the serpentinite if the serpentinite is deformed under differential stress. Olivine does not form in serpentinite under hydrostatic stress. Hence, we hypothesize that olivine formation is controlled by reaction kinetics and that the kinetics are locally faster in serpentinite that deforms under differential stress.We elaborate a 2D HMC numerical algorithm that can simulate dehydration and olivine generation in a deforming serpentinite [Schmalholz et al., 2023]. The algorithm is based on a staggered finite difference discretization and employs a matrix-free, pseudo-transient iterative solver. Furthermore, the algorithm is programmed in the Julia language, employs the ParallelStencil package, and runs on GPUs. We discuss three major numerical challenges: First, the treatment of large changes in solid density during the generation of olivine by serpentinite dehydration. Second, the treatment of large temporal and spatial gradients in the unknowns, such as porosity and fluid pressure. Third, the treatment of strongly nonlinear relations between unknowns and parameters, such as the relations between density and fluid pressure, porosity and permeability, and porosity and rock viscosity. We implement several mathematical formulations for the reaction kinetics and discuss which formulation can explain the laboratory results best. We further discuss potential numerical benchmarks of such HMC algorithms for modelling the coupling of chemical reactions, fluid flow and rock deformation. ReferencesSchmalholz, S. M., E. Moulas, L. Räss, and O. Müntener (2023), Serpentinite Dehydration and Olivine Vein Formation During Ductile Shearing: Insights From 2D Numerical Modeling on Porosity Generation, Density Variations, and Transient Weakening, Journal of Geophysical Research: Solid Earth, 128(11), e2023JB026985, doi:https://doi.org/10.1029/2023JB026985.

Detailing the relationship between stress and reactions in metamorphic rocks has been controversial, and much of the debate has centered on theory. Here, we add to this discussion and make a major advance by showing in time-resolved synchrotron microtomography experiments that a reacting and deforming sample experiencing an elastic differential stress produces a fabric orthogonal to the largest principal stress. This fabric forms very early in the reaction and can be shown to be unrelated to strain. The consequences of this are significant because a non-hydrostatic stress state is a very common geological occurrence. Our data provide the basis for new interpretations of the classical, and enigmatic, serpentine fabrics of Val Malenco, Italy, and Cerro del Almirez, Spain, where we relate the reported fabrics to transient, and cyclical, differential stresses from magma intrusion and the earthquake cycle.

To investigate the effect of carbon-bearing phases on the release of fluids in subducted serpentinites, we performed high-pressure multi-anvil experiments on representative ophicarbonate assemblages over a pressure range from 2.5 GPa to 5 GPa and from 450 °C to 900 °C, across the antigorite-out reaction. Parallel experiments were performed on carbonate-free serpentinites. In all experiments, we monitored and/or controlled the oxygen fugacity. The addition of 20 wt. % CaCO 3 to a serpentinite assemblage at 2.5 GPa is found to decrease the onset of the serpentine dehydration by over 100 °C, in comparison to carbonate-free assemblages. Similarly, the final disappearance of serpentine is also affected by the presence of CaCO 3 . For a bulk CaCO 3 content of 20 wt. %, this causes a decrease in maximum stability of antigorite by 50 °C. For a bulk CaCO 3 content exceeding 25 wt. %, this difference can be as high as 100 °C in warm and 150 °C in cold subduction zones, causing antigorite to be completely dehydrated at 500 °C. This results from the reaction of CaCO 3 with serpentine to form clinopyroxene and Mg-rich carbonates. This reaction, however, causes no discernible decrease in the proportion of carbonate, indicating that the amount of released carbon is insignificant. Whilst CaCO 3 , therefore, influences serpentine stability, there is no significant effect of hydrous phases on the carbonate stability. On the other hand, a MgCO 3 -bearing system shows no significant effects on the serpentinite stability field. Further experiments and oxygen fugacity calculations indicate that graphite is not stable in typical magnetite-bearing serpentinites. The reduction of carbonates to graphite would require oxygen fugacities that are 1–2 log units below those of magnetite-bearing serpentinites. This confirms earlier studies and indicates that reduction of carbonates can only occur through the infiltration of external H 2 -rich fluids.

The trigger mechanism of intermediate depth earthquakes (30 - 300 km) is a long-standing debate. Many studies showed that these seismic events nucleate along a double-seismic zone within the subducting slab. The seismic events of the lower plane coincide with the depth of major dehydration reactions in the lithospheric mantle. Consequently, it is thought that these events are related to the release of fluids. Several scenarios are currently discussed that might lead to brittle deformation. Among these are dehydration embrittlement and dehydration-driven stress transfer.Antigorite is one of the most important candidates for fluid release due to its high H2O content and stability limits within the lower Wadati-Benioff zone. The release of water through antigorite dehydration can be calculated by equilibrium thermodynamics and is mainly a function of temperature. This does, however, not account for deformation (e.g., stored internal strain energy) leading to local variations in the free energy of minerals. In this study we aim to explore the effect of shear stress on the stability of antigorite.We performed high-pressure and high-temperature experiments in a Griggs rig. We used intact drill cores of two different starting materials for our experiments: a foliated and an isotropic antigorite-serpentinite. Both starting materials did not contain relict olivine and/or orthopyroxene. We run our experiments at 620 to 650 °C with a confining pressure of 1.5 GPa and a strain rate of 10-6 s-2. Subsequent analyses of the experimental runs revealed no dehydration products within the bulk sample. However, we observed the formation of ultra-fine grained (< 100 nm) olivine and orthopyroxene along narrow zones, which are orientated 30 to 40 ° with respect to the compression axis. These zones are similar in all runs and independent of the starting material microstructure. We thus propose that shear stress localization within our cylindrical sample triggered the dehydration. Within subduction zones local variations in stress field due to mineralogical or textural heterogeneities could promote dehydration, eventually leading to seismic events through stress transfer.

On Earth, subduction zones facilitate the cycling of volatiles between the Earth’s surface and interior. Volatile cycling has significant effects on the long-term state of the Earth’s climate and tectono-magmatic events, including volcanism and earthquakes. A key stage in the volatile cycle is the devolatilization of the subducting oceanic lithosphere, in which volatiles can escape the previously hydrated rocks. However, it is not well known how efficiently volatiles are transported at this stage. To better understand how volatiles escape at these conditions, we have analyzed the dehydration-related vein networks of the Erro-Tobbio meta-serpentinites (ET-MS), Italy. The ET-MS display well preserved networks of metamorphic olivine veins. These veins are the result of the dehydration reaction of antigorite and brucite to produce H2O and olivine. However, due to the low permeability of serpentinite at depth, the dehydration reaction requires the formation of self-organizing vein networks to allow the produced fluid to escape [1]. Thus, the metamorphic olivine veins in ET-MS may be used as a proxy for fluid flow pathways. We took a multiscale approach to analyzing the network architectures. For microscale (~16 µm voxel size) and mesoscale (~200 µm voxel size) resolutions, X-ray tomography methods are sufficient to visualize the three-dimensional structure of the networks. However, for large scale observations these methods are inapplicable. To solve this, we apply a novel workflow to analyze outcrop scale (~10 m) network systems in three dimensions using only two-dimensional data. By training a generative adversarial network (GAN) with two-dimensional data conditioned by spatial orientation, we can generate statistically representative three-dimensional networks that mimic those of the ET-MS. These representations also display similar characteristics in their respective pore-network-models. With this method, it is possible to produce reasonable three-dimensional approximations of the ET-MS vein networks using only photogrammetry data of the outcrops. In turn, this allows us to extract metrics, such as permeability, that describe the volatile transport efficiency of the ET-MS, and further, how these characteristics change at a broad range of scales. [1] Plümper et al. (2017) Nature Geoscience 10(2), 150-156.

Many metamorphic rocks have a fabric. What is often not clear is how much deformational or metamorphic processes contributed to the formation of these fabrics. Are foliations always the result of strain? When does intrinsic crystallographic anisotropy alone lead to the formation of structural elements? Understanding the relative contributions of deformation and metamorphism in rock fabrics is fundamentally important because it is foundational to understanding the role of stress in reacting and deforming rocks.To this end, we make a major advance in our understanding of fabric development in reacting rocks by showing in time-resolved (4D) synchrotron microtomography (µCT) experiments that when a gypsum dehydration reaction occurs in a differentially stressed sample the reaction products develop orthogonally to the largest principal stress. This is an important finding because we can show with our µCT data that this preferred orientation forms early in the reaction and at very small strains (<1%). Using a simple kinematic model we can demonstrate that it cannot have formed because of reorientation during mechanical compaction. It remains to be established if it is nucleation or growth of bassanite that is being affected by the stress or both. Our experiments suggest that metamorphic transformations may be inherently anisotropic when reacting under the influence of a non-hydrostatic stress state. The consequences of this are many. For example, there will be cases in natural rocks where the interpretation of a lineation, foliation or crystallographic preferred orientation as formed by strain may be incorrect. Moreover, the physical properties (e.g. hydraulic and mechanics) of metamorphic rocks could also be significantly anisotropic from early in a transformation. Mass transport pathways might initialise as channelled or partitioned conduits which would have an impact during subduction and in thin-skinned tectonics. Our data reveal a critical new finding related to the very common geological occurrence of reacting rocks experiencing a differential stress.

We have performed in situ time‐of‐flight neutron diffraction experiments to examine the uptake of deuterium in iron monosulfide at pressures up to 11.4 GPa and temperatures to 1300 K. A D 2 fluid was formed in the experiments through the decomposition of ND 3 BD 3 , resulting in an oxygen fugacity of approximately 1.2 log units below the iron‐wüstite buffer. Deuterium positions and site occupancies were determined in FeS V, using Rietveld refinements of the powder neutron diffraction patterns. Our structural model indicates that two normally unoccupied sites in the P 6 3 / mmc FeS V structure, at Wyckoff positions 6h and 4f , are partially occupied by D atoms, with the latter being more dominant. The deuterium content D x in FeSD X increases with both pressure and temperature over the experimental conditions explored, from 0.126 (14) at 2.3 GPa and 787 K to 1.20 (16) at 9.7 GPa and 1300 K. The unit‐cell volume expansion per deuterium atom is 1.53 ± 0.16 Å 3 at 6.9 GPa and 960 K, which is smaller than that determined for metallic iron phases at similar conditions. The variation in unit‐cell volume indicates that most deuterium is lost from FeS V upon temperature quenching at high‐pressures. By fitting the obtained FeS V deuterium site occupancies to a thermodynamic model, estimates for the hydrogen contents of iron monosulfide at conditions and oxygen fugacities consistent with the base of the cratonic lithosphere can be made. This results in values in the range of 1,700–2,700 ppm, which contribute to approximately 2–3 ppm hydrogen in the bulk mantle.

ABSTRACT Serpentinites play an important role in the delivery of water into subduction zones. In addition, serpentinites also contain ferric Fe and can transport significant redox potential. We present high-pressure and high-temperature experiments and Mössbauer spectroscopy measurements on natural lizardite and antigorite samples equilibrated at various oxygen fugacities in order to quantify the relationship between the oxygen fugacity f(O2) and the Fe3+/Fetot ratio in these two phases. In antigorite, Fe3+ partitions into the octahedral site and is charge balanced by tetrahedral Al. In lizardite, tetrahedral Fe3+ is observed only at low temperature as well as under high f(O2), whereas Fe3+ prefers the octahedral site at temperatures exceeding 500 °C and at 3 to 5 GPa. Although metastable, lizardite remains in redox equilibrium in our experiments at conditions above the lizardite to antigorite phase transformation at 300 °C and demonstrates a similar stability to antigorite. The Al concentration of lizardite is found to be temperature dependent, and it was possible to reequilibrate the Fe3+/Fetot ratio of lizardite from 0.1 to 0.9 by using redox buffers such as Fe metal, graphite, graphite–calcite, Re–ReO2 and Ru–RuO2. Our experiments on antigorite demonstrate that antigorite does not adjust its Al concentration on experimental time scales. Since Fe3+ is charge balanced by Al, it was also not possible to manipulate the Fe3+/Fetot ratio of antigorite. The coexisting phases, however, show chemical equilibration with this antigorite composition. We have retrieved the standard Gibbs energy for Fe3+- and Al-endmembers of antigorite and lizardite and calculated the metamorphic evolution of subducting serpentinites. The lizardite to antigorite transformation does not cause a decrease in the bulk Fe3+/Fetot ratio under f(O2) buffered conditions, in contrast to observations from some natural settings, but does result in the formation of additional magnetite due to antigorite having a lower Fe3+/Fetot ratio than lizardite at equilibrium. If the f(O2) of antigorite serpentinite is buffered during subduction, such as due to the presence of graphite and carbonate, the bulk Fe3+/Fetot ratio decreases progressively. On the other hand, in a closed system where the bulk serpentinite Fe3+/Fetot ratio remains constant, the f(O2) increases during subduction. In this scenario, the f(O2) of an antigorite serpentinite with a typical Fe3+/Fetot ratio of 0.4 increases from the fayalite–magnetite–quartz to the hematite–magnetite f(O2) buffer during dehydration. These f(O2) results confirm earlier inferences that fluids produced by antigorite dehydration may not contain sufficient oxidised sulphur species to oxidise the mantle wedge. Sufficiently high levels of f(O2) to mobilise oxidised sulphur species may be reached upon antigorite dehydration, however, if closed system behaviour maintains a high bulk redox potential across the lizardite to antigorite phase transformation. Alternatively, oxidation of the mantle wedge might be achieved by oxidising agents from sources in subducted oceanic crust and sediments.

For decades, millions of people have been waiting for the "big one" in either Tokyo, Istanbul or Los Angeles. No one is able to predict where and when such a disaster will happen first. People expect it as a self-evident fact, will experience it as a stroke of fate, and will speak of it as a “tragedy”.Our aim cannot be (yet?) to predict where and when major earthquakes will occur, because that would amount to claiming to be able to announce in advance where and when the lightning strikes. Nevertheless, we believe that our understanding of the seismic process and associated risk should greatly benefit from the following question: what parameters control whether a dynamic rupture nucleates, grows or stops?It is crucial to understand the processes and conditions causing the initial stages of catastrophic rock tearing under pressure, and the interplay between mineral- and tectonic-scale factors. Both fluid percolation events and transformation-driven stress transfers can trigger mechanical instabilities ultimately causing rupture nucleation. And once a rupture has nucleated, similar processes should also occur within the damage zone and modulate the ability of small ruptures to “self-propagate” towards large seismic events.Our lack of understanding is considerable about the exact conditions for rupture nucleation and dynamic propagation. While observational methods help image mechanical instabilities, laboratory experiments provide insights on the physics of the lubrication processes enabling seismic faults to grow under pressure. Unfortunately, there remains a significant gap in scientific communication between researchers using different analytical methods or conceptual views.Here we outline transdisciplinary connections between the contributions to the session. From seismology to electrical conductivity measurements in the laboratory, from field geology to numerical modelling, from machine learning to mineralogy, from geodesy to mineral physics, here we walk on the frontier of knowledge in order to reshape the central questions that we need to ask to further investigate the rupture phenomenon.