Search

SUMMARYThe Opalinus Clay (OPA) is a clay-rich formation considered as a potential host rock for radioactive waste repositories and as a caprock for carbon storage in Switzerland. Its very low permeability (10−19 to 10−21 m2) makes it a potential sealing horizon, however the presence of faults that may be activated during the lifetime of a repository project can compromise the long-term hydrological confinement, and lead to mechanical instability. Here, we have performed laboratory experiments to test the effect of relative humidity (RH), grain size (g.s.) and normal stress on rate-and-state frictional properties and stability of fault laboratory analogues corresponding to powders of OPA shaly facies. The sifted host rock powders at different grain size fractions ( 35 MPa is dictated by a transition from localized to distributed deformation and (2) water vapour content does not affect the deformation mechanisms and dilation, whereas it decreases steady-state friction (μss), and enhances fault stability.

The earthquake Magnitude-Frequency-Distribution, FMD, is usually modelled with the Gutenberg-Richter relation law, where the b-value controls the relative rate of small and large earthquakes. b-value has been documented to show an inverse dependence on differential stress, it increases with the fault roughness or during fluid-induced earthquakes. For some seismic sequences a near real-time characterization of the b-value has been used to discriminate between foreshocks and aftershocks. Here we examine the influence on b-values of different lithologies hosting earthquakes.In general, seismicity not only localizes along the major structures where mainshocks nucleate, but it can be also distributed within volumes of the seismogenic layer characterized by different lithologies. For the Mw 6.5 2016–2017 Central Italy seismic sequence, the lithology can be properly defined by seismic reflection profiles. Here the fractured carbonate of the Apennines, located at almost 1-2 km and 4-6 km of depth, are characterized by b-values ranging between 1.3 and 1.4 that can be diagnostic of brittle dominated deformation. At 2-4 km and 6-10 km of depth, the Triassic Evaporites showing a bimodal brittle-ductile deformation and compartmentalized fluid overpressure (documented in deep boreholes) are linked to high b-values, in the range of 1.5-1.65 reaching 1.80 for clustered swarms. Between 10-12 km of depth the phyllosilicate rich basement, with its predominant velocity strengthening behaviour, is hosting small magnitude earthquakes with b-values around 1.4. Our results indicate that away from the large earthquake faults, characterized by a stress dependent elasto-frictional rheology, FMD are strongly controlled by rock lithology and style of deformation.

In the last 15 years, activities for geo-energy production are associated to subsurface fluid injection in enhanced geothermal systems, for enhanced oil recovery, for the disposal of wastewater or for carbon dioxide capture and storage. In several regions, M>3 earthquakes occurred following fluid injection, and some of these earthquakes have caused extensive damage, putting geo-energy production projects at risk of being discontinued. Evaluating the conditions under which fluid injection can induce earthquakes is therefore important to safeguard local infrastructures and to ensure continuity of geo-energy projects. To shed light on the effect of fluid injection on a fault located in the proximity of a reservoir, we implemented into the Q-DYN seismic cycle simulator the fluid diffusion equation (one-way coupling). We ran models of seismic cycles on a rate-and-state-dependent fault under a quasi-dynamic approximation, and we developed a systematic study to assess how fault frictional heterogeneities, the stress state of the fault upon injection, the timing of injection relative to the phase of the seismic cycle and factors controlling fluid flow, i.e. permeability, porosity, flow-rate, influence fault slip behavior and earthquake magnitude. Our results show that localized pore-pressure perturbations allow us to gain deeper physical insight into the propagation and arrest of earthquake ruptures and that changes in the fault physical properties can promote a spectrum of fault slip behavior and recorded magnitudes.

Serpentinites are poly-mineralic rocks distributed almost ubiquitously in active tectonic regions worldwide. They are composed of rheologically weak (lizardite and crysotile) and strong (e.g., magnetite and pyroxene) phases. In particular, lizardite typically shows low friction coefficients and is supposed to localise deformation along weak shear zones characterized by aseismic behaviour. Major faults hosting serpentinite lithologies are characterised by seismic activity, tremors, and other slip modes. We advance the hypothesis that low strain domains, which are enriched in rheologically strong phases, can act as potential site of nucleation of unstable slip as the result of the velocity-dependent rheology of magnetite-rich serpentinites. Through an experimental and microstructural approach, we explore the different mechanisms whose interplay controls the complex behaviour of these lithologies.For this study we collected natural samples of lizardite-magnetite rich serpentinites within the low strain domains of the Elba Island ophiolites (Italy). Rocks were characterised, powdered, and deformed in a set of shear experiments at four different normal stresses (25, 50, 75 and 100 MPa) in the biaxial apparatus BRAVA. The experiments consist of an initial phase of sliding at 10 μm/s, a slide-hold-slide test, and two series of velocity stepping (sliding velocity from 0.1 to 300 μm/s). Fundamental parameters to quantify the frictional properties of serpentinites are individuated in the (a-b) value, the critical slip distance Dc, and the critical stiffness kc, which is derived by their combination.The material shows friction values of ~0.4 with velocity weakening behaviour and negative frictional healing. The module of the negative (a-b) parameter increases neatly with decreasing sliding velocity while Dc decreases, causing kc to rise. At low velocities (< 3 μm/s) sliding is unstable and the fault undergoes stick-slip behaviour. This is explained by the increase of the critical stiffness to values higher than the loading system stiffness. This systematic change of mechanical properties and fault slip behaviours with sliding velocity is interpreted to be the result of the time-dependent arrangement of grains in a heterogeneous experimental fault architecture.Back-scattered SEM images of the principal slip zones of recovered samples support this hypothesis. Elasto-frictional behaviour is controlled by the build-up of a partial (granular) load-bearing framework of strong magnetite grains, while visco-frictional rheology is controlled by the (phyllosilicatic) anastomosing and foliated lizardite matrix. At low sliding velocities, the granular phase interacts creating force chains thus promoting frictional instabilities. At higher velocities, dilation promotes the activity of throughgoing weaker phyllosilicate planes thus favouring stable slip.Our experiments shed light on the role of fault rock heterogeneity in nucleating dynamic slip in nature as well as in controlling the slip mode during earthquakes or slow-slip events in serpentinite terrains.

After many large earthquakes aftershocks activity can reach depths greater than the base of the seismogenic zone that is defined by background seismic activity. This observation is generally explained by strain rate induced embrittlement associated with the increase of post-mainshock strain rate, which favors a transition from ductile to brittle behavior. However, the underlying physical processes are not well understood. Here we integrate geological and geophysical data for the 2016–2017 Central Italy seismic sequence with laboratory experiments to provide a geological and physical interpretation for the post-mainshock transient deepening of the base of the seismogenic zone.The base of the seismogenic zone in the central-northern Apennines is set typically at 9-10 kilometers and corresponds to the top of the phyllitic basement. Structural studies on exhumed basement rocks highlight a heterogeneous basement fabric consisting of competent, 10 to 200 m wide, quartz-rich lenses surrounded by an interconnected and frictionally weak phyllosilicate-rich matrix. The matrix controls the overall rheology of the basement due to its interconnectivity, and promotes aseismic deformation because its rate-strengthening behavior.Following each major (Mw > 5.5) event of the 2016–2017 sequence, a dramatic and abrupt increase in seismic rate is observed below 10 km, hence within the basement. Here we document the presence of seismicity clusters made of more than 4 earthquakes and characterized by small magnitude (Mw < 2.5), small dimensions (< 500 meters), small temporal duration (< 14 days) and a swarm-like behavior. Furthermore, these clusters are often represented by multiple or repeating earthquakes with a cross correlation coefficient higher than 0.7 for all the three components. These observations suggest that the increase of shear stressing rate within the basement is responsible for deepening of seismicity. To further explore this idea, we performed laboratory experiments on rocks from exhumed outcrops of basement rocks. We found that shear stressing rate promotes accelerated creep on the phyllosilicate-rich matrix and dynamic instabilities on the quartz-rich gouge belonging to the lenses. Our integrated analysis suggests that the mainshocks of the 2016-2017 seismic sequence promoted an increase of shear stressing rate within the basement allowing the phyllosilicate-rich matrix to creep faster hence favoring the loading and the repeated failures of locked seismogenic patches represented by the quartz-rich lenses.

Pressure, temperature, and infilling fluids influence the petrophysical properties and the associated damaging processes of rocks at all scales. Moreover, each fluid-rock system possesses peculiar mechanical behaviours being these particularly complex in carbonate rocks hosting fluids. In this work, we analyze the laboratory results of deformed clean and hydrocarbon-filled limestones under varying pressure and temperature, providing links between recorded physical properties (seismic velocity), fluid behavior, and damaging. We focus on carbonate-bearing reservoir (Bolognano Formation) rocks, sampled in the Majella massif (Central Italy) that represents a very good analogue for buried carbonate reservoirs. This reservoir is composed by calcarenites with connected porosity of about 20% saturated by hydrocarbon in the solid state at the outcrop conditions. We performed hydrostatic, triaxial and true-triaxial deformation tests up to a temperature of 100º C and a confining pressure up to 100 MPa on both clean and naturally hydrocarbon-filled limestone samples. Results show increasing seismic velocity and Young’s modulus with increasing confining pressures for both clean and saturated samples as expected. However, different results are observed when the temperature is increased. At low temperatures saturated samples show larger seismic velocity and rigidity with respect to clean samples whilst at higher temperatures the opposite occurs. In particular, when temperature is rised up to 100º C the Young’s modulus of the saturated samples dramatically decreases, being this coupled by a clear volume reduction even during hydrostatic tests (no differential stress applied). Accordingly, microstructural observations highlight grain crushing related to a large amount of randomly distributed cracks within saturated samples. On the contrary, clean samples are characterized by few microfractures, pointing out the primary role played by liquid hydrocarbons. These observations are in good agreement with meso and microstructural features observed on outcropping hydrocarbon-filled carbonate-bearing faults. The presence of fluid hydrocarbons (high temperature) severely weakens the rock promoting fracturing whilst at lower temperature the presence of solid hydrocarbons increases the mechanical properties of hydrocarbon-bearing rocks. These observations have a large impact for the petrophysical characterization of reservoirs and for the understanding microscale to mesoscale mechanisms of deformation and fluids movement along deformed rock volumes.

Faults can be reactivated by fluid injection and pore pressure increase in the rock volumes surrounding the fault zone. Induced earthquakes represent only one of the possible responses of active faults to pore pressure perturbations, since other strain transients characterize the spectrum of fault-slip behavior. A series of fluid injection experiments, designed and undertaken in the framework of the ERC Fault Activation and Earthquake Ruptures (FEAR) project, will be conducted in the Bedretto Underground Laboratory for Geosciences and Geoenergies (BULGG, Switzerland) to understand fault reactivation processes on a target well-identified fault zone. Small and accessible faults are to be instrumented to monitor deformation and seismicity during both fluid injection and fault reactivation.The mineralogical, microstructural, and hydraulic properties of the target fault zone are investigated to characterize the fault-slip behavior. Characterization of the frictional response is achieved through a suite of laboratory rock-deformation experiments using both double-direct and rotary experimental apparatuses. Fault stimulation by fluid pressurization was also simulated in laboratory by using an injection protocol reliable for the in-situ hydraulic stimulation and consisting of stepwise pore fluid pressure increase. Experiments undertaken at low velocity with the double-direct apparatus (BRAVA) suggest that the selected fault, composed of mixed phyllosilicate-granular materials, is frictionally stable but yet can be dynamically reactivated by hydraulic stimulation.Experiments were also performed on the fault gouge from the target fault with the rotary experimental apparatus (SHIVA). First, we apply half of the stresses measured at depth in the underground laboratory to accomplish the operating capability of the apparatus: 12 MPa normal stress, 7.5 MPa confining pressure and 1.5 MPa pore fluid pressure. Second, we imposed a slip rate of 10-5 m/s for 0.01 m to have an equally compacted and textured layer. Third, we applied a shear stress so that an equivalent slip tendency of 0.35 is achieved (ca. 2.7 MPa), and kept it constant. We then increased stepwise the pore fluid pressure by 0.1 MPa every 150 s. This allows the spontaneous nucleation of slip events. After fault reactivation, the maximum slip velocity was set to 0.1 m/s. The fluid injection sequence results in a first reactivation (R1). Thanks to the nominally infinite slip available in SHIVA we run a second injection sequence up to a reactivation (R2).Our experiments show two different styles of reactivation between R1 and R2. R1 reactivation is abrupt, with slip rate accelerating up to 0.1 m/s. Instead, R2 has a stage in which slip rate oscillates (0.5-3 mm/s) just before the last step of pore pressure increase leads to acceleration to 0.1 m/s. This would suggest a role for the shear fabric developed during the first reactivation, in which extensive grain size reduction might have led to stiffening of the fault, responsible of the oscillatory slip. This frictional behavior suggests the importance of considering the effect of texture development during multiple cycles of seismic slip. The generalization of our data and observations will contribute to shed light on the mechanics of faults and induced earthquakes by fluid pressure increase.

In the last decades, it has been observed that faults can slip both by slow aseismic creep and seismic events (i.e., earthquakes). Between these two slip modes, a wide variety of fault slip behavior can be observed, including low-frequency earthquakes, slow slip events and tremors. This wide variety of slip modes can radiate seismic energy at different frequencies whose content may be linked to the physical mechanisms at play. In the laboratory, it is possible to reproduce the entire spectrum of fault slip modes by modulating the loading stiffness of the apparatus depending on the critical fault rheologic stiffness (i.e. k/kc). This technique allows us to study, under controlled laboratory conditions, the acoustic signature of different fault slip modes to infer the physical mechanisms at their origin. To shed light on the nucleation mechanisms and seek for reliable precursors to failure of different slip modes, we performed friction experiments on powders that differ for granulometry and grain shape (i.e., glass beads with a grainsize < 150 µm; and quartz powders MinUSil with an average grain size of 10.5 µm), to simulate fault gouge. The experiments were conducted in a double direct shear configuration, instrumented with an array of piezoelectric sensors to record continuously Acoustic Emissions (AEs) at high recording rate (~10MHz). The experiments are performed at a constant displacement rate of 10 µm/s and using a spring to reduce the apparatus stiffness k, to match the critical fault rheological stiffness (kc). Following this procedure we obtain slow slip events (i.e., k = kc) and fast events (i.e. k

Large earthquakes take place on mature faults with hundreds of meters to kilometres of cumulative slip. At shallow depths, the fault zone is generally composed of non-cohesive rock wear products, often referred to as gouge. Seismic and aseismic slip occur in this fault gouge and fracture/brecciation of the wall rock and damage zone can add to the fault gouge as part of the wear process. Gouge thickness generally increases linearly with the cumulative fault shear displacement and laboratory work shows that gouge tends to stabilize fault frictional stability. Previous works show that frictional stability of simulated fault gouge varies as a function of shear displacement. The stability evolution is interpreted as a consequence of the degree of shear localisation within the simulated fault gouge: the more the deformation is localized, the more the fault slip is unstable. This implies that for bare rock surfaces, unstable behaviour is expected as the deformations are forced to be localized at the interface between the two sheared surfaces.On natural faults at large shear displacement (or for faults having a high gouge production rate), a competition must take place between 1) the localization of the deformation at rock-on-rock surfaces, 2) the delocalization of deformation due to gouge production and wall rock brecciation, 3) fault zone lithification and frictional healing and 4) shear localization within the gouge and wear material. The competition and interaction between these phenomena are modulated by cumulative fault slip, temperature and fluid chemistry. In turn, this competition may influence the frictional stability of faults with increasing shear displacement, and thus, their potential seismic activity.To characterise the influence of shear displacement on fault stability, constant velocity and velocity step experiments were performed to large displacement. Two initially intact rocks were chosen as starting material: a high porosity Fontainebleau sandstone and a low porosity quartzite. These samples represent very different resistances to abrasion (i.e., wear production with slip) for the same initial mineral composition (< 95% quartz), which allows us to investigate wear and wear rate on fault stability. Additionally, simulated quartz gouge was tested for comparison. Mechanical data are analysed within the rate-and-state framework, and post-mortem microscopic analyses of the sample were performed. For initially bare surface experiments a threshold shear displacement is required to transition from stable to unstable sliding. Stick-slip events (laboratory earthquakes) evolve systematically as a function of fault zone shear displacement. The inversion of the rate-and-state parameters shows that shear displacement has a dominant influence on both (a-b) and Dc. For all the faults tested, (a-b) decreases with increasing shear displacement. For high wear rates and simulated gouge, Dc decreases with increasing shear displacement. However, for low wear rate faults, Dc is constant within the tested shear displacement. These results demonstrate that, under the tested boundary conditions, fault stability varies systematically with fault maturity and in particular that shear displacement and strain localization are the dominant parameters controlling fault slip stability.

Slip along pre-existing faults in the Earth’s crust occurs whenever the shear stress resolved on the fault plane overcomes fault frictional strength, potentially generating catastrophic earthquakes. The coupling between shear stress and normal stress during fault loading depends on 1) the orientation of the fault within the stress field and 2) the tectonic setting. In compressional settings, a load-strengthening path occurs because along thrust faults the increase in shear stress is coupled with an increase in effective normal stress. On the contrary, in extensional settings, the increase in shear stress is coupled with a decrease in effective normal stress, resulting in load-weakening paths for normal faults.Analytical approaches to evaluate the potential for fault reactivation are generally based on the assumption that faults are ideal planes, characterized by zero thickness and constant friction, embedded in homogeneous isotropic elastic media. However, natural faults typically host thick fault cores and highly fractured damage zones, which can compact or dilate under different loading paths (i.e., different coupling between normal and shear stress). In addition, in most laboratory friction experiments, the fault is loaded under constant or increasing normal stress and at optimal orientation for reactivation. Here, we present laboratory experiments simulating reactivation of thick gouge-bearing faults that experienced different loading paths.Our results show that the differential stress required for reactivation strongly differs from theoretical predictions, and unfavourably oriented faults appear systematically weaker, especially when a thick gouge layer is present. Before reactivation fault zone compacts in load-strengthening paths whereas dilation is observed in load-weakening path. Upon fault reactivation at comparable normal stress, load-strengthening promotes stable creep whereas load-weakening results in accelerated slip. Our study highlights the importance of fault thickness and loading path in fault hydromechanical coupling and stability with significant implications for fluid circulation within fault zones and earthquake mechanics.

Earthquake nucleation has been understood as controlled by the frictional properties of fault zones. Mature fault zones host abrasive wear products, such as gouges, which result from the frictional sliding occurring in successive slip events. Shear localization in fault gouges is strongly dependent on, among others, fault mineralogical composition and grain size distribution, originating a wide variety of microstructural textures that may be related to different types of fault motion from aseismic creep, slow earthquakes to fast slip events. Additionally, within a fault, one can encounter different stages of maturity, ranging from an incipient and poorly-developed fault zone (i.e. discontinuous and thin gouge layer) to a mature fault zone that has experienced a lot of wear from previous sliding events (i.e. well-developed gouge layer). The localization of deformation within a mature gouge layer has been identified as possibly responsible for mechanical weakening and as an indicator of a change in stability within the fault.To gain insights on the role of grain size distribution, and thus fault maturity, in slip behavior and fault rheology, we performed friction experiments on quartz fault gouge in a double direct shear configuration using a biaxial apparatus (BRAVA at INGV in Rome, Italy). The experiments were performed at a constant normal stress of 40 MPa and under 100% humidity. We investigated different sliding velocities, from 10 µm/s to 1 mm/s, to assess time-dependent physical processes. Different bi-disperse mixtures of quartz were sheared to reproduce different initial grain size distributions within the fault (F110, average grain size and Min-u-sil, average grain size ). Samples were carefully collected at the end of the experiments to prepare thin sections for microstructural analyses.A first set of experiments was performed increasing the proportion between smaller and larger particles within a homogeneous blend. The friction evolves from a strain-hardening behavior for a sample with only F110 to a slip-weakening one for the one with only Min-u-sil. The difference in rheology is observable in the analyzed microstructures. Particularly, the two end members clearly show comminution and localization along boundary shear planes, whereas mixtures of the two sizes of particles only present a more diffused deformation. In the second set of experiments, we sheared gouges with a horizontal layering of the two grain sizes and observed different behaviors in terms of friction and rheology. These layered gouges present strain hardening behavior, with a strengthening part corresponding to the material of the layer in contact with the sliding block and a steady-state part with slightly higher friction than for the homogeneous mixtures.These results give important information on the connection between grain size distribution, shear localization, and the resulting fault slip behavior. In this context, the proportion between small/large particles and their distribution and percentages within the fault plays an important role in controlling fault rheology. We also complete our knowledge by using Discrete Element Method, simulating gouge sliding with different grain scale properties (size, distribution, cementation…), and observing a detailed evolution of shear localizations.

Fault zone properties can change significantly during the seismic cycle in response to stress changes, microcracking and wall rock damage. Lab experiments show consistent changes in elastic properties prior to and after lab earthquakes (EQ) and previous works show that machine learning/deep learning (ML/DL) techniques are successful for capturing such changes. Here, we apply DL techniques to assess whether similar changes occur during the seismic cycle of tectonic EQ. The main motivation is to generalize lab-based findings to tectonic faulting, to predict failure and identify precursors. The novelty is that we use EQ traces as probing signals to estimate the fault state.We train DL model to distinguish foreshocks, aftershocks and time to failure of the Mw 6.5 2016 Norcia EQ in central Italy, October 30th 2016. We analyze a 25-second window of 3-component data around the P- and S-wave arrivals for events near the Norcia fault with M>0.5 and ±2 months before/after the Norcia mainshock. Normalized waveforms are used to train a Convolutional Neural Network (CNN). As a first task we divide events into two classes (foreshocks/aftershocks), and then refine the classification as a function of time-to-failure (TTF) for the mainshock. Our DL model perform very well for TTF classification into 2, 4, 8, or 9-classes for the 2 months before/after the mainshock. We explore a range of seismic ray paths near, through, and away from the Norcia mainshock fault zone. Model performance exceeds 90% for most stations. Waveform investigations show that wave amplitude is not the key factor; other waveform properties dictate model performance. Models derived from seismic spectra, rather than time-domain data, are equally good. We challenged the model in several ways to confirm the results. We found reduced performance in training the model with the wrong mainshock time and by omitting data immediately before/after the mainshock. Foreshock/aftershock identification is significantly degraded also by removing high frequencies (filtering seismic data above 25 Hz). We tested data from different years to understand seasonality at individual stations for the time period September to December and removed these effects. Comparing these seasonality effects defined from noise with our EQ results shows that foreshocks/aftershocks for the 2016 Norcia mainshock are well resolved. Training with data containing EQ offers a huge increase in classification performance over noise only, proving that EQ signals are the sole that enable assessing timing as a function of the fault status. To confirm our results and understand which stations are able to detect changes of fault properties we perform a further test cleaning the signals from the seasonality by confounding the DL with a shuffled noise (adversarial training).We conclude that DL is able to recognize variations in the stress state and fracture during the seismic cycle. The model uses EQ-induced changes in seismic attenuation to distinguish foreshocks from aftershocks and time to failure. This is an important step in ongoing efforts to improve EQ prediction and precursor identification through the use of ML and DL.

Summary of triaxial tests on anhydrites cylindrical samples.

Tectonic fault zones are subject to normal stress variations with a wide range of spatiotemporal scales, resulting in stress field alteration. These perturbations can spread over a wide range of frequencies and amplitudes from the high frequency passage of seismic waves generated by earthquakes, to the low frequency of solid earth tides and underground fluid injection cycles. As a result of these normal stress perturbations, critically stressed faults can be reactivated. The resulting slip mode is then controlled by fault friction and elastic properties of the surrounding rock. Existing works show that complex behaviors may arise from the interplay between friction changes with slip and slip rate and stress perturbations.To shed light on the mechanics of fault dynamic triggering we performed experiments in a Biaxial Apparatus in a Double Direct Shear configuration under critically stable stiffness conditions (K/Kc~1). We used powdered quartz gouge (Min-U-Sil 40) as starting material, and conducted experiments at reference normal stress of σn = 10-13.5 MPa. After shearing the material and reaching a steady state sliding, normal stress oscillations were applied with various amplitudes, varying from A = 0.5-2 MPa, and periods, T = 0.5-50 s. In addition, we used the laboratory derived friction parameters as input for forward modeling using Rate-and-State friction laws in order to assess if these laws can explain our data. Our results show that creeping faults, under critical stiffness conditions, are sensitive to normal stress perturbations showing a variety of slip behaviors depending on amplitude and frequency of the oscillations:Oscillation frequency has a major effect on fault stability. Low and high frequencies cause a Coulomb-like response of the shear stress, that is accompanied by a complex frictional response with slow events and period doubling. At the critical frequency predicted by the Rate-and-State friction, we observe dynamic weakening resulting in regular stick-slip events. Oscillation amplitude also plays a role with the main effect depending on the magnitude of the perturbation. Using a modified Rate-and-State equation (Linker and Dieterich, 1992), we are able to accurately model the laboratory data. Our results show that normal stress perturbation on a laboratory creeping fault, at critical stiffness condition, can reproduce the entire spectrum of fault slip behavior depending on the oscillation properties.

Analysis of seismicity can illuminate active fault zone structures but also pervasive deformation occurring within large volumes of the seismogenic zone. Here we show that, for the Mw 6.5 2016–2017 Central Italy seismic sequence, seismicity not only localizes along SW-dipping normal faults, but also occurs within larger volumes of Triassic Evaporites, TE, composed of alternated anhydrites and dolostones. This off-fault and diffuse microseismicity shows a different frequency-magnitude distribution than on-fault seismicity along the structures hosting the largest events. We interpret that, during the sequence, shear strain-rate increase promoted widespread ductile deformation within TE that light-up with diffuse microseismicity. This interpretation is supported by field and laboratory observations showing that TE background ductile deformation is complex and dominated by distributed failure and folding of the anhydrites associated with boudinage fracturing and faulting of dolostones. Our results indicate that ductile crustal deformation can cause off-fault and diffuse microseismicity, which obeys to different scaling laws than on-fault seismicity on structures characterized by elasto-frictional stick-sip behaviour.

Giant rockslides are widespread and sensitive to hydrological forcing, especially in climate change scenarios. They creep slowly for centuries and then can fail catastrophically posing major threats to society. However, the mechanisms regulating the slow-to-fast transition toward their catastrophic collapse remain elusive. We couple laboratory experiments on natural rockslide shear zone material and in situ observations to provide a scale-independent demonstration that short-term pore fluid pressure variations originate a full spectrum of creep styles, modulated by slip-induced undrained conditions. Shear zones respond to pore pressure increments by impulsive acceleration and dilatancy, causing spontaneous deceleration followed by sustained steady-rate creep. Increasing pore pressure results in high creep rates and eventual collapse. Laboratory experiments quantitatively capture the in situ behavior of giant rockslides and lay physically-based foundations to understand the collapse of giant rockslides., Giant rockslides creep slowly for centuries and then can fail catastrophically, posing major threats to society. Here, the authors use observational and experimental evidence to quantitatively capture the full spectrum of giant rockslide behaviour until collapse, that is modulated by hydro-mechanical response to short-term fluid pressure perturbations.

Empirical scaling relationships between fault or slip dimensions and earthquake magnitudes are often used to assess the maximum possible earthquake magnitude of a territory. Upon the assumption of the reactivability of any fault, the earthquake magnitudes derived from the surface fault length (FLEM) are compared at the national scale in Italy against catalogued magnitudes. FLEMs are obtained by considering a comprehensive fault dataset regardless of fault age, stress field orientation, strain rate, etc. In particular, (1) a comprehensive catalogue of all known faults is compiled by merging the most complete databases available; (2) FLEM is then derived from fault length; and (3) the resulting FLEMs are compared (i.e. the mathematical difference) with catalogued earthquake magnitudes. Results show that the largest FLEMs as well as the largest differences between FLEMs and catalogued magnitudes are observed for poorly constrained faults, mainly inferred from subsurface data. It is suggested that these areas have to be further characterized to better estimate fault dimension and segmentation and hence properly assess the FLEM. Where, in contrast, the knowledge of faults is geologically well constrained, the calculated FLEM is often consistent with the catalogued seismicity, with the 2σ value of the distribution of differences being 1.47 and reducing to 0.53 when considering only the Mw≥6.5 earthquakes. Our work highlights areas, in Italy, where further detailed studies on faults are required.

We performed a suite of experiments aimed at examining the frictional properties of unaltered basalts at conditions considered to be representative of slip at shallow depths in volcano-tectonic environments and in-situ geo-energy basaltic sites. Scientific drilling and field studies on exhumed subsurface faults and fractures analogues suggest that, frictional sliding in basalts can occur in shear zones within a volume of wear debris or along localized joint surfaces. To illuminate how microstructural heterogeneities affect the nucleation of slip instabilities in basalts, we sheared simulated fault gouge and bare rock surfaces at low normal stresses (4–30 MPa) at ambient temperature, under room-dry and wet conditions. We performed velocity steps (0.1–300 μm/s) and slide-hold-slides (30–3000 s holds) to determine the frictional stability and healing properties of basalts. In all the tests, we observed high friction coefficient associated with important frictional restrengthening. Overall, our results show that microstructural heterogeneities strongly affect the friction velocity dependence of basalts: while for normal stresses ≥10 MPa, shear localization accompanied by cataclasis and grain size reduction favors the transition to velocity weakening behavior of powdered samples, on bare surfaces gouge production during shearing promotes a transition to a velocity strengthening behavior. Our results imply that at the tested conditions, friction instabilities may promptly nucleate in shear zones where deformation within (unaltered) basaltic gouge layers is localized, such as those located along volcanic flanks, while joint surfaces characterized by rough rock-on-rock contacts are less prone to unstable slip, which is suppressed at velocities ≥10 μm/s.

Many rock deformation experiments used to characterize the frictional properties of tectonic faults are performed on powdered fault rocks or on bare rock surfaces. These experiments have been fundamental to document the frictional properties of granular mineral phases and provide evidence for crustal faults characterized by high friction. However, they cannot entirely capture the frictional properties of faults rich in phyllosilicates. Numerous studies of natural faults have documented fluid-assisted reaction softening promoting the replacement of strong minerals with phyllosilicates that are distributed into continuous foliations. To study how these foliated fabrics influence the frictional properties of faults we have: 1) collected foliated phyllosilicate-rich rocks from natural faults; 2) cut the fault rock samples to obtain solid wafers 0.8-1.2 cm thick and 5 cm x 5 cm in area with the foliation parallel to the 5x5cm face of the wafer; 3) performed friction tests on both solid wafers sheared in their in situ geometry and powders, obtained by crushing and sieving and therefore disrupting the foliation of the same samples; 4) recovered the samples for microstructural studies from the post experiment rock samples; and 5) performed microstructural analyses via optical microscopy, scanning and transmission electron microscopy. Mechanical data show that the solid samples with well-developed foliation show significantly lower friction in comparison to their powdered equivalents. Micro- and nano-structural studies demonstrate that low friction results from sliding along the foliation surfaces composed of phyllosilicates. When the same rocks are powdered, frictional strength is high, because sliding is accommodated by fracturing, grain rotation, translation and associated dilation. Friction tests indicate that foliated fault rocks may have low friction even when phyllosilicates constitute only a small percentage of the total rock volume, implying that a significant number of crustal faults are weak.

Tectonic faults typically break in a single rupture mode within the range of styles from slow slip to dynamic earthquake failure. However, in increasingly well‐documented instances, the same fault segment fails in both slow and fast modes within a short period, as in the sequences that culminated in the 2011 Mw = 9.0 Tohoku‐Oki, Japan, and 2014 Mw = 8.2 Iquique, Chile, earthquakes. Why slow slip alternates with dynamic rupture in certain regions but not in others is not well understood. Here, we integrate laboratory experiments and numerical simulations to investigate the physical conditions leading to cycles where the two rupture styles alternate. We show that a bifurcation takes place near the stability transition with sequences encompassing various rupture modes under constant loading rate. The range of frictional instabilities and slip cycles identified in this study represents important end‐members to understand the interaction of slow and fast slip on the same fault segment. Ministry of Education (MOE) Published version This work is funded by the Earth Observatory of Singapore and by the Singapore Ministry of Education and also the 2017 Stephen Riady Funding from the Earth Observatory of Singapore (M4430260.B50.500000).

Understanding the factors controlling fracture frequency distribution can greatly improve the assessment of fluid circulation in fault damage zones, with evident implications for fault mechanics, hydrogeology and hydrocarbon exploration. This is particularly important for relay zones that are usually characterized by strong damage and structural complexity. We investigated the fracture frequency within an outcrop adjacent to the front fault segment of a relay ramp, hosted within peritidal carbonates that forms part of the Tre Monti fault (Central Italy). We analysed the distribution of fracture frequency in the outcrop through (1) scanlines measured in the field, (2) oriented rock samples, and (3) scan-areas performed on a virtual outcrop model. Fracture frequency increases with distance from the front segment of the relay ramp. Moreover, supratidal and intertidal carbonate facies exhibit higher fracture frequency than subtidal limestones. This trend of increased fracture frequency has two main explanations. (1) The number of subsidiary faults and their associated damage zones increases moving away from the front segment. (2) The supratidal and intertidal carbonate facies content increases toward the centre of the relay ramp. Our results indicate that the fracture frequency pattern is very complex in relay ramps hosted in shallow-water limestones and that its prediction necessitates a good control on structures and sedimentary facies distribution.

Dynamic changes in the stress field during the seismic cycle of tectonic faults can control frictional stability and the mode of fault slip. Small perturbation in the stress field, like those produced by tidal stresses can influence the evolution of frictional strength and fault stability with the potential of triggering a variety of slip behaviors. However, an open question that remains still poorly understood is how amplitude and frequency of stress changes influence the triggering of an instability and the associated slip behavior, i.e. slow or fast slip.Here we reproduce in the laboratory the spectrum of fault slip behaviors, from slow-slip to dynamic stick-slip, by matching the critical fault rheologic stiffness (kc) with the surrounding stiffness (k). We investigate the influence of normal stress variations on the slip style of a quartz rich fault gouge at the stability boundary, i.e. k/kc slightly less than one, by adopting two techniques: 1) instantaneous step-like changes and 2) sinusoidal variations in normal stress. For the latter case, modulations of normal stress were chosen to have amplitudes greater, less or equal to the typical stress drop observed during unperturbed experiments. Also, the period was varied to be greater, less or equal to the typical recurrence time of laboratory slow-slip events. During the experiments, we continuously record ultrasonic wave velocity to monitor the microphysical state of the fault. We find that frictional stability is profoundly affected by variation in normal stress giving rise to a variety of slip behaviors. Furthermore, during strain accumulation and fabric development, changes in normal stress permanently influence the microphysical state of the fault gouge increasing kc and producing a switch from slow to fast stick-slip. Our results indicate that perturbations in the stress state can trigger a variety of slip behaviors along the same fault patch. These results have important implications for the formulation of constitutive laws in the framework of rate- and state- friction, highlighting the necessity to account for the microphysical state of the fault in order to improve our understanding of frictional stability.

Despite the numerous advantages of storing CO2 into basalts by dissolving carbon dioxide into water prior to its injection, the large amount of H2O required for this operation poses an increased risk of fluid overpressure into the fault/fracture networks, and renders the seismicity analysis pivotal to upscale this storage method to voluminous basaltic occurrences diffused worldwide.To deepen our knowledge on the frictional strength, stability, and the healing properties of basalt-built faults, we carried out friction tests on basalts from Mt. Etna using the biaxial deformation machine, BRAVA, and the rotary-shear apparatus, SHIVA (HP-HT laboratory of INGV-Rome, Italy). Specimens were selected for their relative abundance of olivine and pyroxene crystals, i.e. the main sources of divalent cations in silicate rocks necessary to trap CO2 into basalts.Experiments were performed both on synthetic powdered samples and bare surfaces, at room-dry and water drained-saturated conditions, at room temperature and pressure. Bare surfaces consisted in basalt slabs and hollowed cylinders, which were mounted on BRAVA and SHIVA apparatus, respectively. Samples were subjected to 5 to 30 MPa normal stress (σn) for powdered samples and in the range 5 to 10 MPa for bare surfaces.At the investigated normal stresses, frictional sliding data obey Byerlee’s law for friction, with the friction coefficient µ = 0.59 – 0.78. Differences in μ mainly reflect sample variability, different experimental configurations, sample geometry, and, to a lesser extent, the boundary conditions (dry/wet). However, in detail, basalt slabs are generally characterized by the highest friction coefficient and hollow cylinders exhibit a slight increase in friction coefficient with increasing shear displacement, due to the progressive slip hardening resulting from gouge production during frictional sliding.Velocity step increases were conducted on BRAVA after steady values of friction were attained (~ 6.5-7.5 mm for gouge and ~ 3 mm slip for bare surfaces) and consisted in velocity sequences from 0.1 to 300 µm s-1, with ~ 500 μm displacement for each step. Rate-and-state friction experiments show opposite mechanical behavior between bare surfaces and synthetic fault gouge: while bare surfaces experience a transition from rate-weakening at low sliding velocity (V) to rate-strengthening behavior at higher V without any clear dependence on the applied σn, gouge revealed a negative trend of (a-b) with shear velocity at σn > 5 MPa and a velocity-weakening behavior at V ≥ 30 µm s-1, regardless of experimental conditions. We ascribe this different behavior to shear delocalization owing to frictional wear production in bare surfaces, and to shear localization accompanied by grain size reduction along the Riedel R1 and boundary B shear zones in fault gouges, as also confirmed by microstructural analysis.The velocity weakening behavior of fault gouge, coupled with the fast healing rates retrieved from slide-hold-slide experiments (500 µm displacement cumulated at V = 10 µm/s followed by hold times from 30 to 3000 s), define high strength zones that are potentially seismogenic. Conversely, velocity strengthening behavior of bare surfaces promotes aseismic creep at V ≥ 100 µm s-1, regardless of the faster restrengthening compared to fault gouge.

Giant rockslides creep for centuries and then can fail catastrophically posing major threats to society. There is growing evidence that creeping landslides are widespread worldwide and extremely sensitive to hydrological forcing, especially in climate change scenarios. Rockslide creep is the results of progressive rock failure processes, leading to rock damage accumulation, permeability enhancement and strain localization within basal shear zones similar to tectonic faults. As shear zone accumulate strain, they become thicker and less permeable, favoring the development of perched aquifers. Since then, the creep behavior of mature rockslides becomes dominated by hydro-mechanical interaction with external triggers, e.g. rainfall and snowmelt. However, the mechanisms regulating the slow-to-fast transition toward their catastrophic collapse remain elusive, and statistical and simplified mathematical models used for collapse prediction are usually unable to account for the full spectrum of observed slip behaviors.Here we couple laboratory experiments on natural rockslide shear zone material, sampled from high quality drillcores, and in situ observations (groundwater level and surface displacement) to investigate the mechanism of rockslide response to short-term pore pressure variations within basal shear zones at the Spriana rockslide (Italy). Using a biaxial apparatus within a pressure vessel, we characterized the strength and permeability of the phyllosilicate-rich shear zone material at in situ stress, as well as the rate and state frictional properties for shear rates typical of the slow-to-fast transition of real rockslides. Then we carried out non-conventional pore pressure-step creep experiments, in which shear stress is maintained at subcritical levels and pore pressure is increased stepwise while monitoring shear zone slip and dilatancy until runaway failure.Our results, that are quantitatively consistent with in situ monitoring observations, provide a scale-independent demonstration that short-term pore pressure variations originate a full spectrum of creep styles, modulated by slip-induced undrained conditions. Shear zones respond to fluid pressure increments by impulsive acceleration and dilatancy, causing spontaneous deceleration followed by sustained steady-rate creep. Increasing fluid pressure results in high creep rates and eventual collapse. Laboratory experiments quantitatively capture the in situ behavior of giant rockslides, providing physically-based basis to improve forecasting models for giant mature rockslides in crystalline rocks.

Pore fluid overpressures in active fault systems can drive fluid flow and cause fault weakening and seismicity. In return, deformation accommodated by different modes of failure (e.g. brittle vs. ductile) also affects fault zone permeability and, hence, fluid flow and pore fluid pressure distribution. Current numerical simulation techniques model how fluid flow controls fault reactivation and associated seismicity. However, the control exerted by pore fluid pressure on the transition from slow aseismic fault sliding to fast seismic sliding, during the earthquake nucleation phase, is still poorly understood. Here, we model overpressured, supercritical CO2 fluid flow in natural faults, where non-linear, complex feedback between fluid flow, fluid pressure and fault deformation controls the length of the nucleation phase of an earthquake and the duration of the interseismic period. The model setup is an analogue for recent seismic source events in the Northern Apennines of Italy (e.g. M w 6.0 1997-98 Colfiorito and M w 6.5 2016 Norcia earthquakes). Our modelling results of Darcy fluid flow show that the duration of the nucleation phase can be reduced by orders of magnitude, when realistic models of fault zone architecture and pore pressure- and deformation-dependent permeability are considered. In particular, earthquake nucleation phase duration can drop from more than 10 years to a few days/minutes, while the seismic moment can decrease by a factor of 6. Notably, the moment of aseismic slip ( M 0 = 10 9 N m ) obtained during the nucleation phase modelled in our study is of the same order as the detection limit of local strain measurements using strain meters. These findings have significant implications for earthquake early warning systems, as the duration and moment of the nucleation phase will affect the likelihood of timely precursory signal detection. Interestingly, aseismic slip has been measured up to a few months before some recent large earthquakes, although in a different tectonic context than the model developed here, rekindling interest in the nucleation phase of earthquakes. In addition, our results have important implications for short and long term earthquake forecasting, as crustal fluid migration during the interseismic period may control fault strength and earthquake recurrence intervals.

Fractures constitute the main pathway for fluids in fault damage zones hosted in low-porosity rocks. Understanding the factors controlling fracture distribution is hence fundamental to better assess fluids circulation in fault damage zones, with evident implications for fault mechanics, hydrogeology and hydrocarbon exploration. Being usually characterized by a strong damage and structural complexity, this is of particularly importance for relay zones.We integrated classical and modern structural geology techniques to investigate the factors controlling fracture distribution within a portion of a relay ramp damage zone pertaining to the Tre Monti fault (Central Italy). The damage zone is hosted within peritidal carbonates and located at the footwall of the relay ramp front segment. We analysed the distribution of the fracture density in the outcrop through (1) scanlines measured in the field, (2) oriented rock samples, and (3) scan-areas performed on a virtual outcrop model obtained by aerial structure-from-motion.Our results highlight structural and lithological control on fracture distribution. Scanlines and virtual scan-areas show that fracture density increases with the distance from the front segment of the relay ramp. Moreover, all the methods highlight that supratidal and intertidal carbonate facies exhibit higher fracture density than subtidal limestones.This apparently anomalous trend of fracture density, that increases moving away from a main fault segment, has two main explanations. (1) The damage is associated with the relay ramp development: approaching the centre of the relay ramp (i.e., moving away from the front segment) an increase in the number of subsidiary faults with their associated damage zones promotes high fracture densities. (2) The increase in fracture density can be attributed to the increasing content in supratidal and intertidal carbonate facies that are more abundant in the centre of the relay ramp.Our results provide important suggestions for factors controlling fracture distribution and fluid flow within relay ramps hosted by shallow water limestones. We show that the trend of fracture distribution with respect to a main fault is not easily predictable in presence of a relay ramp, because it can be modulated by the subsidiary faults formation and slip during the relay ramp development. Moreover, carbonate facies play a non-negligible role in fracture distribution within fault zones hosted in shallow-water carbonates.

A great number of earthquakes occur within thick carbonate sequences in the shallow crust. At the same time, carbonate fault rocks exhumed from a depth < 6 km (i.e., from seismogenic depths) exhibit the coexistence of structures related to brittle (i.e., cataclasis) and ductile deformation processes (i.e., pressure-solution and granular plasticity). We performed friction experiments on water-saturated simulated carbonate-bearing faults for a wide range of normal stresses (from 5 to 120 MPa) and slip velocities (from 0.3 to 100 μm/s). At high normal stresses (σn > 20 MPa) fault gouges undergo strain-weakening, that is more pronounced at slow slip velocities, and causes a significant reduction of frictional strength, from μ = 0.7 to μ = 0.47. Microstructural analysis show that fault gouge weakening is driven by deformation accommodated by cataclasis and pressure-insensitive deformation processes (pressure solution and granular plasticity) that become more efficient at slow slip velocity. The reduction in frictional strength caused by strain weakening behaviour promoted by the activation of pressure-insensitive deformation might play a significant role in carbonate-bearing faults mechanics.

One of the major challenges regarding the disposal of radioactive waste in geological formations is to ensure isolation of radioactive contamination from the environment and the population. Shales are suitable candidates as geological barriers. However, the presence of tectonic faults within clay formations put the long-term safety of geological repositories into question. In this study, we carry out frictional experiments on intact samples of Opalinus Clay, i.e. the host rock for nuclear waste storage in Switzerland. We report experimental evidence suggesting that scaly clays form at low normal stress (≤20 MPa), at sub-seismic velocities (≤300 μm/s) and is related to pre-existing bedding planes with an ongoing process where frictional sliding is the controlling deformation mechanism. We have found that scaly clays show a velocity-weakening and -strengthening behaviour, low frictional strength, and poor re-strengthening over time, conditions required to allow the potential nucleation and propagation of earthquakes within the scaly clays portion of the formation. The strong similarities between the microstructures of natural and experimental scaly clays suggest important implications for the slip behaviour of shallow faults in shales. If natural and anthropogenic perturbations modify the stress conditions of the fault zone, earthquakes might have the potential to nucleate within zones of scaly clays controlling the seismicity of the clay-rich tectonic system, thus, potentially compromising the long-term safeness of geological repositories situated in shales.

Graphitization, or the progressive maturation of carbonaceous material, is considered an irreversible process. Thus, the degree of graphite crystallinity, or its structural order, has been calibrated as an indicator of the peak metamorphic temperatures experienced by the host rocks. However, discrepancies between temperatures indicated by graphite crystallinity versus other thermometers have been documented in deformed rocks. To examine the possibility of mechanical modifications of graphite structure and the potential impacts on graphite thermometry, we performed laboratory deformation experiments. We sheared highly crystalline graphite powder at normal stresses of 5 and 25 megapascal (MPa) and aseismic velocities of 1, 10 and 100 µm s−1. The degree of structural order both in the starting and resulting materials was analyzed by Raman microspectroscopy. Our results demonstrate structural disorder of graphite, manifested as changes in the Raman spectra. Microstructural observations show that brittle processes caused the documented mechanical modifications of the aggregate graphite crystallinity. We conclude that the calibrated graphite thermometer is ambiguous in active tectonic settings.

The physical characterization of carbonate–bearing normal faults is fundamental for resource development and seismic hazard. Here we report laboratory measurements of density, porosity, Vp, Vs, elastic moduli and permeability for a range of effective confining pressures (0.1-100 MPa), conducted on samples representing different structural domains of a carbonate-bearing fault. We find a reduction in porosity from the fault breccia (11.7% total and 6.2% connected) to the main fault plane (9% total and 3.5% connected), with both domains showing higher porosity compared to the protolith (6.8% total and 1.1% connected). With increasing confining pressure, P-wave velocity evolves from 4.5 km/s to 5.9 km/s in the fault breccia, is constant at 5.9 km/s approaching the fault plane and is low (4.9 km/s) in clay-rich fault domains. We find that while the fault breccia shows pressure sensitive behaviour (a reduction in permeability from 2*10−16 m2 to 2*10−17 m2), the cemented cataclasite close to the fault plane is characterized by pressure independent behaviour (permeability 4*10−17 m2). Our results indicate that the deformation processes occurring within the different fault structural domains influence the physical and transport properties of the fault zone. In-situ Vp profiles match well the laboratory measurements demonstrating that laboratory data are valuable for implications at larger scale. Combining the experimental values of elastic moduli and frictional properties it results that at shallow crustal levels M ≤ 1 earthquakes are less favoured, in agreement with earthquake-depth distribution during the L'Aquila 2009 seismic sequence that occurred on carbonates.

We studied the frictional behaviour and deformation products of large (20 cm × 20 cm bare surfaces) experimental limestone faults. We sheared samples in a direct shear configuration, with an imposed normal force of 40–200 kN and shear velocity of 10 μm/s. The steady-state shearing of these surfaces yielded a coefficient of friction 0.7

Some faults are considered strong because their strength is consistent with the Coulomb criterion under Byerlee's friction, 0.6 μ 0.85 . In marked contrast, numerous studies have documented significant fault weakening induced by fluid-assisted reaction softening that generally takes place during the long-term evolution of the fault. Reaction softening promotes the replacement of strong minerals with phyllosilicates. Phyllosilicate development within foliated and interconnected fault networks has been documented at different crustal depths, in different tectonic regimes and from a great variety of rock types, nominating fluid-assisted reaction softening as a general weakening mechanism within the seismogenic crust. This weakening originates at the grain-scale and is transmitted to the entire fault zone via the interconnectivity of the phyllosilicate-rich zones resulting in a friction as low as 0.1 μ 0.3 . Collectively, geological data and results from laboratory experiments provide strong supporting evidence for structural and frictional heterogeneities within crustal faults. In these structures, creep along weak and rate-strengthening fault patches can promote earthquake nucleation within adjacent strong and locked, rate-weakening portions. Some new frontiers on this research topic regard: 1) when and how a seismic rupture nucleating within a strong patch might propagate within a weak velocity strengthening fault portion, and 2) if creep and slow slip can be accurately detected within the earthquake preparatory phase and therefore represent a reliable earthquake precursor.

Minor fault geometry and kinematics within relay ramps is strongly related to the stress field perturbations that can be produced when two major fault segments overlap and interact. Here we integrate classical fieldwork and interpretation of a virtual outcrop to investigate the geometry and kinematics of subsidiary faults within a relay ramp along the Tre Monti normal fault in the Central Apennines. Although the Tre Monti fault strikes parallel to the regional extension (NE-SW) it shows predominant dip-slip kinematics, suggesting a NW-SE oriented extension acting at sub-regional scale (1–10 km). Conversely, the slickenlines collected on the front segment of the relay ramp highlight right-lateral kinematics. The subsidiary faults in the relay ramp show a complex geometry (variable attitudes) and slickenlines describe multiple kinematics (left-lateral, dip-slip, right-lateral), independently of their orientation. Our fault slip analysis indicates that a local stress field retrieved from the kinematic inversion of the slickenlines collected on the front segment, and likely promoted by the interaction between the overlapping fault segments that bound the relay zone, can explain most of the geometry and kinematics of the subsidiary faults. Further complexity is added by the temporal interaction with both the regional and sub-regional stress fields.

Faults in the brittle crust constitute preexisting weakness zones that can be reactivated depending on their friction, orientation within the local stress field, and stress field magnitude. Analytical approaches to evaluate the potential for fault reactivation are generally based on the assumption that faults are ideal planes characterized by zero thickness and constant friction. However, natural faults are complex structures that typically host thick fault rocks. Here we experimentally investigate the reactivation of gouge-bearing faults and compare the resulting data with theoretical predictions based on analytical models. We simulate preexisting faults by conducting triaxial experiments on sandstone cylinders containing saw-cuts filled with a clay-rich gouge and oriented at different angles, from 30 degrees to 80 degrees, to the maximum principal stress. Our results show the reactivation of preexisting faults when oriented at 30 degrees, 40 degrees, and 50 degrees to the maximum principal stress and the formation of a new fracture for fault orientations higher than 50 degrees. Although these observations are consistent with the fault lock-up predicted by analytical models, the differential stress required for reactivation strongly differs from theoretical predictions. In particular, unfavorable oriented faults appear systematically weaker, especially when a thick gouge layer is present. We infer that the observed weakness relates to the rotation of the stress field within the gouge layer during the documented distributed deformation that precedes unstable fault reactivation. Thus, the assumption of zero-thickness planar fault provides only an upper bound to the stress required for reactivation of misoriented faults, which might result in misleading predictions of fault reactivation.

The presence of weak phyllosilicates in mature carbonate fault zones has been invoked to explain weak faults. However, the relation between frictional strength, fault stability, mineralogical composition, and fabric of fault gouge, composed of strong and weak minerals, is poorly constrained. We used a biaxial apparatus to systematically shear different mixtures of shale (68% clay, 23% quartz and 4% plagioclase) and calcite, as powdered gouge, at room temperature, under constant normal stresses of 30, 50, 100 MPa and under room-dry and pore fluid-saturated conditions, i.e. CaCO3-equilibrated water. We performed 30 friction experiments during which velocity-stepping and slide-hold-slide tests were employed to assess frictional stability and to measure frictional healing, respectively. Our frictional data indicate that the mineralogical composition of fault gouges significantly affects frictional strength, stability, and healing as well as the presence of CaCO3-equilibrated water. Under room-dry condition, the increasing shale content determines a reduction in frictional strength, from μ = 0.71 to μ = 0.43, a lowering of the healing rates and a transition from velocity-weakening to velocity-strengthening behavior. Under wet condition, with increasing shale content we observe a more significant reduction in frictional strength (μ = 0.65–0.37), a near-zero healing and a velocity strengthening behavior. Microstructural investigations evidence a transition from localized deformation promoted by grain size reduction, in calcite-rich samples, to a more distributed deformation with frictional sliding along clay-enriched shear planes in samples with shale content greater than 50%. For faults cutting across sedimentary sequences composed of carbonates and clay-rich sediments, our results suggest that clay concentration and its ability to form foliated and interconnected networks promotes important heterogeneities in fault strength and slip behavior.

The presence of calcite in and near faults, as the dominant material, cement, or vein fill, indicates that the mechanical behaviour of carbonate-dominated material likely plays an important role in shallow- and mid-crustal faulting. To better understand the behaviour of calcite, under loading conditions relevant to earthquake nucleation, we sheared powdered gouge of Carrara Marble, >98 per cent CaCO3, at constant normal stresses between 1 and 100 MPa under water-saturated conditions at room temperature. We performed slide-hold-slide tests, 1–3000 s, to measure the amount of static frictional strengthening and creep relaxation, and velocity-stepping tests, 0.1–1000 μm s–1, to evaluate frictional stability. We observe that the rates of frictional strengthening and creep relaxation decrease with increasing normal stress and diverge as shear velocity is increased from 1 to 3000 μm s–1 during slide-hold-slide experiments. We also observe complex frictional stability behaviour that depends on both normal stress and shearing velocity. At normal stresses less than 20 MPa, we observe predominantly velocity-neutral friction behaviour. Above 20 MPa, we observe strong velocity-strengthening frictional behaviour at low velocities, which then evolves towards velocity-weakening friction behaviour at high velocities. Microstructural analyses of recovered samples highlight a variety of deformation mechanisms including grain size reduction and localization, folding of calcite grains and fluid-assisted diffusion mass transfer processes promoting the development of calcite nanograins in the highly deformed portions of the experimental fault. Our combined analyses indicate that calcite fault gouge transitions from brittle to semi-brittle behaviour at high normal stress and slow sliding velocities. This transition has important implications for earthquake nucleation and propagation on faults in carbonate-dominated lithologies.