Does hydrologic circulation mask frictional heat on faults after large earthquakes?

Knowledge of frictional resistance along faults is important for understanding the mechanics of earthquakes and faulting. The clearest in situ measure of fault friction potentially comes from temperature measurements in boreholes crossing fault zones within a few years of rupture. However, large temperature signals from frictional heating on faults have not been observed. Unambiguously interpreting the coseismic frictional resistance from small thermal perturbations observed in borehole temperature profiles requires assessing the impact of other potentially confounding thermal processes. We address several issues associated with quantifying the temperature signal of frictional heating including transient fluid flow associated with the earthquake, thermal disturbance caused by borehole drilling, and heterogeneous thermal physical rock properties. Transient fluid flow is investigated using a two‐dimensional coupled fluid flow and heat transport model to evaluate the temperature field following an earthquake. Simulations for a range of realistic permeability, frictional heating, and pore pressure scenarios show that high permeabilities (>10−14 m2) are necessary for significant advection within the several years after an earthquake and suggest that transient fluid flow is unlikely to mask frictional heat anomalies. We illustrate how disturbances from circulating fluids during drilling diffuse quickly leaving a robust signature of frictional heating. Finally, we discuss the utility of repeated borehole temperature profiles for discriminating between different interpretations of thermal perturbations. Our results suggest that temperature anomalies from even low friction should be detectable at depths >1 km 1 to 2 years after a large earthquake and that interpretations of low friction from existing data are likely robust.Keywords: Frictional heat on faults, Hydrologic circulationKeywords: Frictional heat on faults, Hydrologic circulatio

Asymmetric brittle deformation at the Pāpaku Fault, Hikurangi Subduction Margin, NZ, IODP Expedition 375

Quantifying fault damage zones provides a window into stress distribution and rheology around faults. International Ocean Discovery Program (IODP) Expeditions 372/375 drilled an active thrust splay fault within the Hikurangi subduction margin. The fault, which is hosted in Pleistocene clastic sediments, is surrounded by brittle fractures and faults as well as ductile deformation features. We find that fracture density in the damage zone enveloping the fault is asymmetric, with the hanging wall showing greater overall fracture density and at greater distances from the fault than the footwall. Furthermore, the peak in fracture density occurs within an area of mesoscale folding and localized slip in the hanging wall rather than adjacent to the main fault zone. We attribute the asymmetry in damage to disparate deformation histories between the hanging wall and footwall, greater ductile deformation within the footwall, and/or dynamic stress asymmetry around a propagating rupture. Damage asymmetry is common at shallow depths in subduction zones and influences the mechanical and hydrological properties of the fault, such as channelized fluid flow and fault stability. Finally, we demonstrate that subduction zone faults show similar damage-displacement scaling as continental faults

International ocean discovery program expedition 375 preliminary report: Hikurangi subduction margin coring and observatories unlocking the secrets of slow slip through drilling to sample and monitor the forearc and subducting plate, 8 March - 5 May 2018

Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expedition 375 was undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough by (1) coring at four sites, including an active fault near the deformation front, the upper plate above the high-slip SSE sourc e region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll Seamount, and (2) installing borehole observatories in an active thrust near the deformation front and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372 (26 November 2017-4 January 2018; see th e Expedition 372 Preliminary Report for further details on the LWD acquisition program). Northern Hikurangi subduction margin SSEs recur every 1-2 years and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. Sampling of material from the sedimentary section and oceanic basement of the subducting plate reveals the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate all the way to the trench, indicating that the shallow thrust fault zone targeted during Expedition 375 may also lie in the SSE rupture area. Hence, sampling at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expedition 375 (together with the Hikurangi subduction LWD component of Expedition 372) was designed to address three fundamental scientific objectives: (1) characterize the state and composition of the incoming plate and shallow plate boundary fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; (2) characterize material properties, thermal regime, and stress conditions in the upper plate above the core of the SSE source region; and (3) install observatories at an active thrust near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface

Slow slip source characterized by lithological and geometric heterogeneity

Slow slip events (SSEs) accommodate a significant proportion of tectonic plate motion at subduction zones, yet little is known about the faults that actually host them. The shallow depth (<2 km) of well-documented SSEs at the Hikurangi subduction zone offshore New Zealand offers a unique opportunity to link geophysical imaging of the subduction zone with direct access to incoming material that represents the megathrust fault rocks hosting slow slip. Two recent International Ocean Discovery Program Expeditions sampled this incoming material before it is entrained immediately down-dip along the shallow plate interface. Drilling results, tied to regional seismic reflection images, reveal heterogeneous lithologies with highly variable physical properties entering the SSE source region. These observations suggest that SSEs and associated slow earthquake phenomena are promoted by lithological, mechanical, and frictional heterogeneity within the fault zone, enhanced by geometric complexity associated with subduction of rough crust

Topographically driven groundwater flow and the San Andreas heat flow paradox revisited

Evidence for a weak San Andreas Fault includes borehole heat flow measurements that show no evidence for a frictionally generated heat flow anomaly and the inferred orientation of σ1 nearly perpendicular to the fault trace. Interpretations of the stress orientation data remain controversial, at least in close proximity to the fault, leading some researchers to hypothesize that the San Andreas Fault is, in fact, strong and that its thermal signature may be removed or redistributed by topographically driven groundwater flow in areas of rugged topography, such as typify the San Andreas Fault system. To evaluate this scenario, we use a steady state, two-dimensional model of coupled heat and fluid flow within cross sections oriented perpendicular to the fault and to the primary regional topography. Our results show that existing heat flow data near Parkfield, California, do not readily discriminate between the expected thermal signature of a strong fault and that of a weak fault. In contrast, for a wide range of groundwater flow scenarios in the Mojave Desert, models that include frictional heat generation along a strong fault are inconsistent with existing heat flow data, suggesting that the San Andreas Fault at this location is indeed weak. In both areas, comparison of modeling results and heat flow data suggest that advective redistribution of heat is minimal. The robust results for the Mojave region demonstrate that topographically driven groundwater flow, at least in two dimensions, is inadequate to obscure the frictionally generated heat flow anomaly from a strong fault. However, our results do not preclude the possibility of transient advective heat transport associated with earthquakes

(Table T6) Consolidation properties and permeabilities of sediments from ODP Leg 205 sites

Vertical permeability and sediment consolidation measurements were taken on seven whole-round drill cores from Sites 1253 (three samples), 1254 (one sample), and 1255 (three samples) drilled during Ocean Drilling Program Leg 205 in the Middle America Trench off of Costa Rica's Pacific Coast. Consolidation behavior including slopes of elastic rebound and virgin compression curves (Cc) was measured by constant rate of strain tests. Permeabilities were determined from flow-through experiments during stepped-load tests and by using coefficient of consolidation (Cv) values continuously while loading. Consolidation curves and the Casagrande method were used to determine maximum preconsolidation stress. Elastic slopes of consolidation curves ranged from 0.097 to 0.158 in pelagic sediments and 0.0075 to 0.018 in hemipelagic sediments. Cc values ranged from 1.225 to 1.427 for pelagic carbonates and 0.504 to 0.826 for hemipelagic clay-rich sediments. In samples consolidated to an axial stress of ~20 MPa, permeabilities determined by flow-through experiments ranged from a low value of 7.66 x 10**-20 m**2 in hemipelagic sediments to a maximum value of 1.03 x 10**-16 m**2 in pelagic sediments. Permeabilities calculated from Cv values in the hemipelagic sediments ranged from 4.81 x 10**-16 to 7.66 x 10**-20 m**2 for porosities 49.9%-26.1%