Extreme hydrothermal conditions at an active plate-bounding fault.

Temperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre. At temperatures above 300-450 degrees Celsius, usually found at depths greater than 10-15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional-mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults.

Seismic Slip-Pulse Experiments Simulate Induced Earthquake Rupture in the Groningen Gas Field.

Rock materials show dramatic dynamic weakening in large-displacement (m), high-velocity (∼1 m/s) friction experiments, providing a mechanism for the generation of large, natural earthquakes. However, whether such weakening occurs during induced M3-4 earthquakes (dm displacements) is unknown. We performed rotary-shear experiments on simulated fault gouges prepared from the source-, reservoir- and caprock formations present in the seismogenic Groningen gas field (Netherlands). Water-saturated gouges were subjected to a slip pulse reaching a peak circumferential velocity of 1.2-1.7 m/s and total displacements of 13-20 cm, at 2.5-20 MPa normal stress. The results show 22%-81% dynamic weakening within 5-12 cm of slip, depending on normal stress and gouge composition. At 20 MPa normal stress, dynamic weakening from peak friction coefficients of 0.4-0.9 to 0.19-0.27 was observed, probably through thermal pressurization. We infer that similar effects play a key role during induced seismic slip on faults in the Groningen and other reservoir systems.

Adverse events of exercise therapy in randomised controlled trials: a systematic review and meta-analysis.

OBJECTIVE: To evaluate the relative risk (RR) of serious and non-serious adverse events in patients treated with exercise therapy compared with those in a non-exercising control group. DESIGN: Systematic review and meta-analysis. DATA SOURCES: Primary studies were identified based on The Cochrane Database of Systematic Reviews investigating the effect of exercise therapy. ELIGIBILITY CRITERIA: At least two of the authors independently evaluated all identified reviews and primary studies. Randomised controlled trials were included if they compared any exercise therapy intervention with a non-exercising control. Two authors independently extracted data. The RR of serious and non-serious adverse events was estimated separately. RESULTS: 180 Cochrane reviews were included and from these, 773 primary studies were identified. Of these, 378 studies (n=38 368 participants) reported serious adverse events and 375 studies (n=38 517 participants) reported non-serious adverse events. We found no increase in risk of serious adverse events (RR=0.96 (95%CI 0.90 to 1.02, I2: 0.0%) due to exercise therapy. There was, however, an increase in non-serious adverse events (RR=1.19 (95%CI 1.09 to 1.30, I2: 0.0%). The number needed to treat for an additional harmful outcome for non-serious adverse events was 6 [95%CI 4 to 11). CONCLUSION: Participating in an exercise intervention increased the relative risk of non-serious adverse events, but not of serious adverse events. Exercise therapy may therefore be recommended as a relatively safe intervention.PROSPERO registration numberCRD42014014819.

Fault zone fabric and fault weakness.

Geological and geophysical evidence suggests that some crustal faults are weak compared to laboratory measurements of frictional strength. Explanations for fault weakness include the presence of weak minerals, high fluid pressures within the fault core and dynamic processes such as normal stress reduction, acoustic fluidization or extreme weakening at high slip velocity. Dynamic weakening mechanisms can explain some observations; however, creep and aseismic slip are thought to occur on weak faults, and quasi-static weakening mechanisms are required to initiate frictional slip on mis-oriented faults, at high angles to the tectonic stress field. Moreover, the maintenance of high fluid pressures requires specialized conditions and weak mineral phases are not present in sufficient abundance to satisfy weak fault models, so weak faults remain largely unexplained. Here we provide laboratory evidence for a brittle, frictional weakening mechanism based on common fault zone fabrics. We report on the frictional strength of intact fault rocks sheared in their in situ geometry. Samples with well-developed foliation are extremely weak compared to their powdered equivalents. Micro- and nano-structural studies show that frictional sliding occurs along very fine-grained foliations composed of phyllosilicates (talc and smectite). When the same rocks are powdered, frictional strength is high, consistent with cataclastic processes. Our data show that fault weakness can occur in cases where weak mineral phases constitute only a small percentage of the total fault rock and that low friction results from slip on a network of weak phyllosilicate-rich surfaces that define the rock fabric. The widespread documentation of foliated fault rocks along mature faults in different tectonic settings and from many different protoliths suggests that this mechanism could be a viable explanation for fault weakening in the brittle crust.

Cathodoluminescence as a tracing technique for quartz precipitation in low velocity shear experiments.

Two simulated gouges (a pure quartz and a quartz-muscovite mixture) were experimentally deformed in a ring shear apparatus at a constant low velocity under hydrothermal conditions favourable for dissolution-precipitation processes. Microstructural analysis using scanning electron microscope cathodoluminescence imaging and cathodoluminescence spectroscopy combined with chemical analysis showed that quartz dissolution and precipitation occurred in both experiments. The starting materials and deformation conditions were chosen so that dissolution-precipitation microstructures could be unambiguously identified from their cathodoluminescence signal. Precipitated quartz was observed as blue luminescent fracture fills and overgrowths with increased Al content relative to the original quartz. In the pure quartz gouge, most of the shear deformation was localized on a boundary-parallel slip surface. Sealing of fractures in a pulverized zone directly adjacent to the slip surface may have helped keeping the deformation localized. In the quartz-muscovite mixture, some evidence was observed of shear-accommodating precipitation of quartz in strain shadows, but predominantly in fractures, elongating the original grains. Precipitation of quartz in fractures implies that the length scale of diffusive mass transfer in frictional-viscous flow is shorter than the length of the quartz domains. Additionally, fracturing might play a more important role than generally assumed. Our results show that cathodoluminescence, especially combined with chemical analysis, is a powerful tool in microstructural analyses of experimentally deformed quartz-bearing material and visualizing quartz precipitation.

Surface microstructures developed on polished quartz crystals embedded in wet quartz sand compacted under hydrothermal conditions.

Intergranular pressure solution plays a key role as a deformation mechanism during diagenesis and in fault sealing and healing. Here, we present microstructural observations following experiments conducted on quartz aggregates under conditions known to favor pressure solution. We conducted two long term experiments in which a quartz crystal with polished faces of known crystallographic orientation was embedded in a matrix of randomly oriented quartz sand grains. For about two months an effective axial stress of 15 MPa was applied in one experiment, and an effective confining pressure of 28 MPa in the second. Loading occurred at 350 °C in the presence of a silica-saturated aqueous solution. In the first experiment, quartz sand grains in contact with polished quartz prism ([Formula: see text]) faces became ubiquitously truncated against these faces, without indenting or pitting them. By contrast, numerous sand-grain-shaped pits formed in polished pyramidal ([Formula: see text]) and ([Formula: see text]) crystal faces in the second experiment. In addition, four-leaved and (in some cases) three-leafed clover-shaped zones of precipitation formed on these prism faces, in a consistent orientation and pattern around individual pits. The microstructures observed in both experiments were interpreted as evidence for the operation of intergranular pressure solution. The dependence of the observed indentation/truncation microstructures on crystal face orientation can be explained by crystallographic control of stress-induced quartz dissolution kinetics, in line with previously published experimental and petrographic data, or possibly by an effect of contact orientation on the stress-induced driving force for pressure solution. This should be investigated in future experiments, providing data and microstructures which enable further mechanism-based analysis of deformation by pressure solution and the effect of crystallographic control on its kinetics in quartz-rich sands and sandstones.

Healing Behavior of Simulated Fault Gouges From the Groningen Gas Field and Implications for Induced Fault Reactivation.

We investigated the frictional strength recovery (healing) and subsequent reactivation and slip-weakening behavior of simulated fault gouges derived from key stratigraphic units in the seismogenic Groningen gas field (N. E. Netherlands). Direct-shear, slide-hold-slide (SHS) experiments were performed at in situ conditions of 100 °C, 40 MPa effective normal stress and 10-15 MPa pore fluid pressure (synthetic formation brine). Sheared gouges were allowed to heal for periods up to 100 days before subsequent reshearing. The initial coefficient of (steady) sliding friction µ was highest in the Basal Zechstein caprock (µ = 0.65 ± 0.02) and Slochteren sandstone reservoir (µ = 0.61 ± 0.02) gouges, and the lowest in the Ten Boer claystone at the reservoir top (µ = 0.38 ± 0.01) and in the Carboniferous shale substrate (µ ≈ 0.45). Healing and subsequent reactivation led to a marked increase (∆µ) in (static) friction coefficient of up to ~0.16 in Basal Zechstein and ~0.07 in Slochteren sandstone gouges for the longest hold periods investigated, followed by a sharp strength drop (up to ~25%) and slip-weakening trajectory. By contrast, the Ten Boer and Carboniferous gouges showed virtually no healing or strength drop. Healing rates in the Basal Zechstein and Slochteren sandstone gouges were significantly affected by the stiffness of different machines used, in line with the Ruina slip law, and with a microphysical model for gouge healing. Our results point to marked stratigraphic variation in healed frictional strength and healing rate of faults in the Groningen system, and high seismogenic potential of healed faults cutting the reservoir and Basal Zechstein caprock units, upon reactivation.

Microphysical Model Predictions of Fault Restrengthening Under Room-Humidity and Hydrothermal Conditions: From Logarithmic to Power-Law Healing.

The maximum fault strength and rate of interseismic fault strengthening ("healing") are of great interest to earthquake hazard assessment studies, as they directly relate to event magnitude and recurrence time. Previous laboratory studies have revealed two distinct frictional healing behaviors, referred to as Dieterich-type and non-Dieterich-type healing. These are characterized by, respectively, log-linear and power-law increase in the strength change with time. To date, there is no physical explanation for the frictional behavior of fault gouges that unifies these seemingly inconsistent observations. Using a microphysical friction model previously developed for granular fault gouges, we investigate fault strengthening analytically and numerically under boundary conditions corresponding to laboratory slide-hold-slide tests. We find that both types of healing can be explained by considering the difference in grain contact creep rheology at short and long time scales. Under hydrothermal conditions favorable for pressure solution creep, healing exhibits a power-law evolution with hold time, with an exponent of ~1/3, and an "apparent" cutoff time (α) of hundreds of seconds. Under room-humidity conditions, where grain contact deformation exhibits only a weak strain-rate dependence, the predicted healing also exhibits a power-law dependence on hold time, but it can be approximated by a log-linear relation with α of a few seconds. We derive analytical expressions for frictional healing parameters (i.e., healing rate, cutoff time, and maximum healing), of which the predictions are consistent with numerical implementation of the model. Finally, we apply the microphysical model to small fault patches on a natural carbonate fault and interpret the restrengthening during seismic cycles.

On interevent time distributions of avalanche dynamics.

Physical systems characterized by stick-slip dynamics often display avalanches. Regardless of the diversity of their microscopic structure, these systems are governed by a power-law distribution of avalanche size and duration. Here we focus on the interevent times between avalanches and show that, unlike their distributions of size and duration, the interevent time distributions are able to distinguish different mechanical states of the system. We use experiments on granular systems and numerical simulations of emulsions to show that systems having the same probability distribution for avalanche size and duration can have different interevent time distributions. Remarkably, these interevent time distributions look similar to those for earthquakes and, if different from an exponential, are indirect evidence of non trivial space-time correlations among avalanches. Our results therefore indicate that interevent time statistics are essential to characterise the dynamics of avalanches.

Microscale cavitation as a mechanism for nucleating earthquakes at the base of the seismogenic zone.

Major earthquakes frequently nucleate near the base of the seismogenic zone, close to the brittle-ductile transition. Fault zone rupture at greater depths is inhibited by ductile flow of rock. However, the microphysical mechanisms responsible for the transition from ductile flow to seismogenic brittle/frictional behaviour at shallower depths remain unclear. Here we show that the flow-to-friction transition in experimentally simulated calcite faults is characterized by a transition from dislocation and diffusion creep to dilatant deformation, involving incompletely accommodated grain boundary sliding. With increasing shear rate or decreasing temperature, dislocation and diffusion creep become too slow to accommodate the imposed shear strain rate, leading to intergranular cavitation, weakening, strain localization, and a switch from stable flow to runaway fault rupture. The observed shear instability, triggered by the onset of microscale cavitation, provides a key mechanism for bringing about the brittle-ductile transition and for nucleating earthquakes at the base of the seismogenic zone.