Fracture surface energy of the Punchbowl fault, San Andreas system

Earthquake energy dissipation in a fracture mechanics framework

Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake's potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics.

Multiscale off-fault brecciation records coseismic energy budget of principal fault zone

Breccia and pulverized rock are typical textures in off-fault damage adjacent to a main seismogenic zone. Previously, by estimating the energy required to advance the rupture in this zone using particle size distribution at sub-millimeter/micrometer scales, we could constrain the energy budget during coseismic events. However, whether microscopic estimation is sufficient to capture surface energy fragmentation during an earthquake and the effect of measurement scale variation on calculation of co-seismic energy partitioning remained unclear. Here, we investigated the mechanism of coseismic off-fault damage based on field and microstructural observations of a well-exposed breccia body in Ichinokawa, Japan. We used in situ clast measurements coupled with thin-section analysis of breccia clasts to estimate the energy budget of the damage zone adjacent to the principal slip zone of the Median Tectonic Line (MTL). The total surface energy density and corresponding surface energy per unit fault for a width of ~ 500 m of the dynamical damage zone were estimated. The moment magnitude estimated based on surface energy was 5.8-8.3 Mw. In Ichinokawa, off-fault fragmentation is initiated by coseismic activity and is followed by fluid activity. Under dynamic fragmentation conditions, the scale is important to calculate the surface energy.

Radon concentration in seawater as a geochemical indicator of submarine fault activity in the Yatsushiro Sea, Japan

This study examined the relationship between radon (222Rn) concentrations in seawater and crustal activity in the Yatsushiro Sea by investigating the submarine fault zone situated at the southern end of the Futagawa-Hinagu fault zone, activated by the 2016 Kumamoto earthquake (M7.3). We conducted an analysis of 222Rn concentration in samples of bottom water just above the seafloor and pore water in sediments, utilizing multiple and piston cores from the Hakuho Maru Expedition KH18-3. The findings revealed significantly elevated 222Rn concentrations in the central sites of the Yatsushiro Sea, coinciding with a high-stress field exhibiting dense active faults. Seismicity analysis revealed heightened moment release and a low b-value post the 2016 Kumamoto earthquake, indicative of increased seismic activity and the potential for substantial earthquakes in the Yatsushiro Sea vicinity. Our results indicate that heightened concentrations of 222Rn in seawater can serve as an effective tracer for identifying and estimating submarine fault activities. Moreover, our research highlights the utility of 222Rn concentrations in detecting active submarine faults and assessing their activity. It contributes to a comprehensive understanding of the potential for significant earthquakes in the Yatsushiro Sea in the future.

On-fault earthquake energy density partitioning from shocked garnet in an exhumed seismic midcrustal fault

The energy released during an earthquake is mostly dissipated in the fault zone and subordinately as radiated seismic waves. The on-fault energy budget is partitioned into frictional heat, generation of new grain surface by microfracturing, and crystal-lattice distortion associated with dislocation defects. The relative contribution of these components is debated and difficult to assess, but this energy partitioning strongly influences earthquake mechanics. We use high-resolution scanning-electron-microscopy techniques, especially to analyze shocked garnet in a fault wall-rock, to provide the first estimate of all three energy components for a seismic fault patch exhumed from midcrustal conditions. Fault single-jerk seismicity is recorded by the presence of pristine quenched frictional melt. The estimated value of energy per unit fault surface is ~13 megajoules per square meter for heat, which is predominant with respect to both surface energy (up to 0.29 megajoules per square meter) and energy associated with crystal lattice distortion (0.02 megajoules per square meter).

A split Hopkinson pressure bar for experimental investigation of dynamic pulverization under very high strain rates

Off-fault damage or pulverized rocks found in large-scale strike-slip faults are of great interest in earthquake research. In order to experimentally investigate rock pulverization, we developed a split Hopkinson pressure bar with compact dimensions and high-speed imaging. The developed experimental setup is capable of generating very high strain rates up to 1320 s-1 with the satisfaction of stress equilibrium, which are essential to reproduce the dynamic pulverization observed in nature and obtain dynamic stress-strain responses accurately. High-speed imaging revealed that cracks initiate and propagate along the grain boundaries at very high speeds, while the dynamic stress-strain response suggested that energy dissipated into the fracture increases with stronger impacts. In addition, we show that the apparatus is capable of producing particle size distributions partly similar to those in naturally pulverized rocks of large-scale strike-slip faults. Thus, our developed system with compact dimensions opens new ways to understand the dynamics of the rock pulverization in off-fault regions of large-scale strike-slip faults.

Fracture-Induced Permeability in Whitby Mudstone

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

Volumetric and shear processes in crystalline rock approaching faulting

Understanding the approach to faulting in continental rocks is critical for identifying processes leading to fracturing in geomaterials and the preparation process of large earthquakes. In situ dynamic X-ray imaging and digital volume correlation analysis of a crystalline rock core, under a constant confining pressure of 25 MPa, are used to elucidate the initiation, growth, and coalescence of microfractures leading to macroscopic failure as the axial compressive stress is increased. Following an initial elastic deformation, microfractures develop in the solid, and with increasing differential stress, the damage pervades the rock volume. The creation of new microfractures is accompanied by propagation, opening, and closing of existing microfractures, leading to the emergence of damage indices that increase as powers of the differential stress when approaching failure. A strong spatial correlation is observed between microscale zones with large positive and negative volumetric strains, microscale zones with shears of opposite senses, and microscale zones with high volumetric and shear strains. These correlations are attributed to microfracture interactions mediated by the heterogeneous stress field. The rock fails macroscopically as the microfractures coalesce and form a geometrically complex 3D volume that spans the rock sample. At the onset of failure, more than 70% of the damage volume is connected in a large fracture cluster that evolves into a fault zone. In the context of crustal faulting dynamics, these results suggest that evolving rock damage around existing locked or future main faults influences the localization process that culminates in large brittle rupture events.

Geochemical Alteration and Mineralogy of Coals under the Influence of Fault Motion: A Case Study of Qi’nan Colliery, China

Geochemical characteristics of rocks in fault zones have been extensively studied, while there are limited studies on coal occurring in fault zones of underground coal mine. In this study, five coal samples were carefully collected from a reverse fault zone in Qi’nan colliery. Systematical detection methods were employed to analyze the different chemical and physical characteristics of fault-related coal samples. Through comparative analysis, the following insights are obtained. Three subdivided fault zones were classified according to the deformation characteristics of coal samples. Frictional heat and strong ductile deformation generated by fault motion led to the dissociation of phenol and carboxyl groups in coal molecules, which sharply decreased the concentrations of elements Co and Mo bound to these functional groups in zone I. The modified pore-cleat system in zone I with higher pore volume and lower permeability allowed solutions containing enriched trace elements to migrate through zone I locally. Concentrations of HREE, MREE and related elements associated with the invasive solutions showed significant positive anomalies in zone I. Precipitation and smearing of clay minerals in zone I led to poorer connectivity. Disruption and delamination of laminar clay minerals by strong compression-shear stress significantly increased the adsorption sites for related elements, especially the HREE and MREE. Nano-scale clay minerals resulting from stress-induced scaly exfoliation also enhanced the retention capability of REE in zone I.

Nanocrystalline Principal Slip Zones and Their Role in Controlling Crustal Fault Rheology

Principal slip zones (PSZs) are narrow (<10 cm) bands of localized shear deformation that occur in the cores of upper-crustal fault zones where they accommodate the bulk of fault displacement. Natural and experimentally-formed PSZs consistently show the presence of nanocrystallites in the <100 nm size range. Despite the presumed importance of such nanocrystalline (NC) fault rock in controlling fault mechanical behavior, their prevalence and potential role in controlling natural earthquake cycles remains insufficiently investigated. In this contribution, we summarize the physical properties of NC materials that may have a profound effect on fault rheology, and we review the structural characteristics of NC PSZs observed in natural faults and in experiments. Numerous literature reports show that such zones form in a wide range of faulted rock types, under a wide range of conditions pertaining to seismic and a-seismic upper-crustal fault slip, and frequently show an internal crystallographic preferred orientation (CPO) and partial amorphization, as well as forming glossy or “mirror-like” slip surfaces. Given the widespread occurrence of NC PSZs in upper-crustal faults, we suggest that they are of general significance. Specifically, the generally high rates of (diffusion) creep in NC fault rock may play a key role in controlling the depth limits to the seismogenic zone.

Phase-transformation-induced lubrication of earthquake sliding

Frictional failure is not possible at depth in Earth, hence earthquakes deeper than 30-50 km cannot initiate by overcoming dry friction. Moreover, the frequency distribution of earthquakes with depth is bimodal, suggesting another change of mechanism at about 350 km. Here I suggest that the change at 30-50 km is from overcoming dry friction to reduction of effective stress by dehydration embrittlement and that the change at 350 km is due to desiccation of slabs and initiation by phase-transformation-induced faulting. High-speed friction experiments at low pressure indicate that exceeding dry friction provokes shear heating that leads to endothermic reactions and pronounced weakening. Higher-pressure studies show nanocrystalline gouge accompanying dehydration and the highest pressure experiments initiate by exothermic polymorphic phase transformation. Here I discuss the characteristic nanostructures of experiments on high-speed friction and high-pressure faulting and show that all simulated earthquake systems yield very weak transformation-induced lubrication, most commonly nanometric gouge or melt. I also show that phase-transformation-induced faulting of olivine to spinel can propagate into material previously transformed to spinel, apparently by triggering melting analogous to high-speed friction studies at low pressure. These experiments taken as a whole suggest that earthquakes at all depths slide at low frictional resistance by a self-healing pulse mechanism with rapid strength recovery.This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

Experimental constraints on dynamic fragmentation as a dissipative process during seismic slip

Various fault damage fabrics, from gouge in the principal slip zone to fragmented and pulverized rocks in the fault damage zone, have been attributed to brittle deformation at high strain rates during earthquake rupture. Past experimental work has shown that there exists a critical threshold in stress-strain rate space through which rock failure transitions from failure along a few discrete fracture planes to intense fragmentation. We present new experimental results on Arkansas Novaculite (AN) and Westerly Granite (WG) in which we quantify fracture surface area produced by dynamic fragmentation under uniaxial compressive loading and examine the controls of pre-existing mineral anisotropy on dissipative processes at the microscale. Tests on AN produced substantially greater new fracture surface area (approx. 6.0 m² g^-1) than those on WG (0.07 m² g^-1). Estimates of the portion of energy dissipated into brittle fracture were significant for WG (approx. 5%), but appeared substantial in AN (10% to as much as 40%). The results have important implications for the partitioning of dissipated energy under extreme loading conditions expected during earthquakes and the scaling of high-speed laboratory rock mechanics experiments to natural fault zones.This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

Heating, weakening and shear localization in earthquake rupture

Field and borehole observations of active earthquake fault zones show that shear is often localized to principal deforming zones of order 0.1-10 mm width. This paper addresses how frictional heating in rapid slip weakens faults dramatically, relative to their static frictional strength, and promotes such intense localization. Pronounced weakening occurs even on dry rock-on-rock surfaces, due to flash heating effects, at slip rates above approximately 0.1 m s^-1 (earthquake slip rates are typically of the order of 1 m s^-1). But weakening in rapid shear is also predicted theoretically in thick fault gouge in the presence of fluids (whether native ground fluids or volatiles such as H₂O or CO₂ released by thermal decomposition reactions), and the predicted localizations are compatible with such narrow shear zones as have been observed. The underlying concepts show how fault zone materials with high static friction coefficients, approximately 0.6-0.8, can undergo strongly localized shear at effective dynamic friction coefficients of the order of 0.1, thus fitting observational constraints, e.g. of earthquakes producing negligible surface heat outflow and, for shallow events, only rarely creating extensive melt. The results to be summarized include those of collaborative research published with Nicolas Brantut (University College London), Eric Dunham (Stanford University), Nadia Lapusta (Caltech), Hiroyuki Noda (JAMSTEC, Japan), John D. Platt (Carnegie Institution for Science, now at *gramLabs), Alan Rempel (Oregon State University) and John W. Rudnicki (Northwestern University).This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion'.

Rock mechanics. Superplastic nanofibrous slip zones control seismogenic fault friction

Understanding the internal mechanisms controlling fault friction is crucial for understanding seismogenic slip on active faults. Displacement in such fault zones is frequently localized on highly reflective (mirrorlike) slip surfaces, coated with thin films of nanogranular fault rock. We show that mirror-slip surfaces developed in experimentally simulated calcite faults consist of aligned nanogranular chains or fibers that are ductile at room conditions. These microstructures and associated frictional data suggest a fault-slip mechanism resembling classical Ashby-Verrall superplasticity, capable of producing unstable fault slip. Diffusive mass transfer in nanocrystalline calcite gouge is shown to be fast enough for this mechanism to control seismogenesis in limestone terrains. With nanogranular fault surfaces becoming increasingly recognized in crustal faults, the proposed mechanism may be generally relevant to crustal seismogenesis.

Rapid acceleration leads to rapid weakening in earthquake-like laboratory experiments

After nucleation, a large earthquake propagates as an expanding rupture front along a fault. This front activates countless fault patches that slip by consuming energy stored in Earth's crust. We simulated the slip of a fault patch by rapidly loading an experimental fault with energy stored in a spinning flywheel. The spontaneous evolution of strength, acceleration, and velocity indicates that our experiments are proxies of fault-patch behavior during earthquakes of moment magnitude (M(w)) = 4 to 8. We show that seismically determined earthquake parameters (e.g., displacement, velocity, magnitude, or fracture energy) can be used to estimate the intensity of the energy release during an earthquake. Our experiments further indicate that high acceleration imposed by the earthquake's rupture front quickens dynamic weakening by intense wear of the fault zone.

Origin, separation and identification of environmental nanoparticles: a review

The biogeochemical and ecological impacts of environmental nanoparticles (ENPs) are some of the fastest growing areas of research today. However, efficient separation and collection of ENPs in natural systems remains difficult. This review article is focused on experimental investigation of separation and identification of ENPs, including nanoparticles with size fractions in the range of <2000, 450 to 2000, 100 to 450 and 1 to 100 nm. An automated ultrafiltration device (AUD) was used successfully to overcome the problem of efficiently collecting ENPs in large quantities in red soils. A significant amount of hematite nanoparticles was present on the surface coating of kaolinite nanoparticles and aggregated hematite nanoparticles overlapping the edge of a kaolinite flake in a size range of 5 to 8 nm. Synchrotron XRD technique is more straightforward and powerful than conventional XRD with oriented specimens and random powder methods for identifying nanoparticles, crystallinity, and particle size in red soils, particularly for the illite, kaolinite, goethite and hematite nanoparticles. The AUD apparatus can be employed to efficiently collect large quantities of soil and related ENPs for investigation of their structural characteristics and surface properties, which have significant impact on weathering reaction pathways, catalysis, the fate of vital elements and environmental pollutants, and ecosystem restoration.

Slip zone and energetics of a large earthquake from the Taiwan Chelungpu-fault Drilling Project

Determining the seismic fracture energy during an earthquake and understanding the associated creation and development of a fault zone requires a combination of both seismological and geological field data. The actual thickness of the zone that slips during the rupture of a large earthquake is not known and is a key seismological parameter in understanding energy dissipation, rupture processes and seismic efficiency. The 1999 magnitude-7.7 earthquake in Chi-Chi, Taiwan, produced large slip (8 to 10 metres) at or near the surface, which is accessible to borehole drilling and provides a rare opportunity to sample a fault that had large slip in a recent earthquake. Here we present the retrieved cores from the Taiwan Chelungpu-fault Drilling Project and identify the main slip zone associated with the Chi-Chi earthquake. The surface fracture energy estimated from grain sizes in the gouge zone of the fault sample was directly compared to the seismic fracture energy determined from near-field seismic data. From the comparison, the contribution of gouge surface energy to the earthquake breakdown work is quantified to be 6 per cent.