Readings in Geophysics Ge 169277 Vito Rubino Heating

Readings in Geophysics – Ge 169/277 Vito Rubino • Heating and weakening of faults during earthquake slip J. R. Rice, JGR, 2006 • Thermal decomposition of carbonates J. Sulem and V. Famin, JGR, 2009 6 th March 2012

Heating and weakening of faults during earthquake slip J. R. Rice, JGR, 2006

Heating and weakening of faults during earthquake slip • Earthquakes occur because fault strength weakens with increasing slip or slip rate • What physical processes determine how weakening occurs? • Relevant weakening processes in large crustal events are suggested to be thermal and to involve: 1) thermal pressurization of pore fluid (at highly stressed frictional microcontacts during rapid slip) 2) flash heating (at highly stressed frictional microcontacts during rapid slip) reduces effective normal stress and shear strength for a given friction coefficient reduces friction coefficient • Macroscopic melting, or possibly gel formation in silica-rich lithologies, may become important at large slip. Predictions How to verify predictions • Strength drop should often be nearly complete at large slip • Onset of melting should be precluted over much (and, for small enough slip, all) of the seismogenic zone implied by by these thermal mechanisms Shear fracture energies Compare with seismic data of crustal events

Extreme localization of slip on mature faults (Punchbowl Fault) [Chester and Chester, 1998; Chester et al. , 2003, 2004; Chester and Goldsby, 2003] 1 mm 5 mm

Field studies that support above morphology • Field observations of mature crustal faults suggest that slip in individual events occurs primarily within a thin shear zone < 1 -5 mm, within a finely granulated, ultracataclastic fault core. • Evidence for such morphology emerged from studied of exhumed and now inactive faults Fault Location North Branch San Gabriel Punchbowl San Andreas system Southern California Papers Chester et al. , 1993 Chester and Chester, 1998; Chester 2004 Median Tectonic Line Japan Wibberly and Shimamoto, 2003 Hanaore Japan Noda and Shimamoto, 2005 Nojima Japan Lockner et al. , 2000

Thermal pressurization of pore fluid • This mechanisms assumes that fluids (water, typically) are present within the fault gouge which shears. • Shear strength t during seismic slip can still be represented by the classical effective law: Normal stress Pressure Frictional heating would cause the fluid (if it was unconstrained rather then caged by the densely packed solid particles) to expand in volume much more than the solid cage A pressure increase must be induced in the pore fluid during slip Unless: • Shear-induced dilatancy of the gouge cage overwhelms thermal expansion effects • Gouge is highly permeable Since sn can typically be assumed to remain constant during slip Strength t is reduced, ultimately, towards zero as shear heating continues to raise temperature so that p approaches sn.

[Slide modified from J. Rice, Caltech GALCIT-Seismo. Lab seminar, 29 March 2007] Specific heat per unit volume of fault gouge Work of shear stress First law of thermodynamics Fourier’s law Fluid mass flux Darcy’s law Energy flux k - permeabilty hf - fluid viscosity (Segall and Rice, 1995, 2004) Assume constant friction coefficient Thermal diffusivity Current mass of pore fluid It describes pore pressure evolution Pore pressure change per unit temperature change Hydraulic diffusivity Increments dp in pore pressure can be related to increments dm in fluid mass and increments d. T in temperature

Predictions of shear strength vs. slip V = 1 m/s Scaled friction strength f = 0. 25 based on flash heating Parameter L* includes the dependence of scaled friction strength on: • porothermoelastic properties of the solid and fluid • slip rate and friction Parameter L* has the unit of length diffusivity divided by slip rate Slip • Continued weakening over a very broad range of of slip d (scaled by L*) [Slide modified from J. Rice, Caltech GALCIT-Seismo. Lab seminar, 29 March 2007]

[Slide from J. Rice, Caltech GALCIT-Seismo. Lab seminar, 29 March 2007]

[Slide modified from J. Rice, Caltech GALCIT-Seismo. Lab seminar, 29 March 2007] • According to the predictions, G increases substantially with slip during an event • The average G in an event (as inferred from seismic slip inversions) also shows a clear tendency to increase with the average slip in the event (Abercombie and Rice, 2005; Rice et al. , 2005; Tinti et al. , 2005).

[Slide from J. Rice, Caltech GALCIT-Seismo. Lab seminar, 29 March 2007]

Other weakening mechanisms • Macroscopic melting (i. e. , when a coherent melt layer has formed along the whole sliding surface) may be relevant for sufficiently large combination of slip and initial effective normal stress. Melts, if sufficiently hot, have a low viscosity and may lubricate faults reducing dynamic friction (Sibson, 1975; Spray, 1993; Brodsky and Kanamori, 2001) • Gel formation. This mechanism has been identified with silica-rich lithologies, when: Sliding produces fine wear debris in presence of water (Goldsby and Tullis, 2002; Di Toro et al. , 2004) • Thermal pressurization caused by rapid emission of CO 2 during the decarbonation. (J. Sulem and V. Famin, JGR, 2009)

Thermal decomposition of carbonates in fault zones: Slip-weakening and temperature-limiting effects J. Sulem and V. Famin, JGR, 2009

Thermal decomposition of carbonates in fault zones: Slip-weakening and temperature-limiting effects 1) Endothermic reaction of calcite decomposition Limits Numerical simulations of seismic slip at depths of 5 -8 km Decarbonation has two critical consequences on seismic slip 2) Rapid emission of CO 2 by decarbonation Increases Decarbonation Coseismic temperature increase to less ≈ 800 °C Chemical decomposition of calcite Ca. CO 3 Ca. O + CO 2 Slip-weakening effect of thermal pressurization

Why studying thermal decomposition of carbonates? Field observations (Pizzino et al. , 2004; Lewicki and Brantley, 2000; Lewicki et al. , 2003; Famin et al. , 2008; Italiano et al. , 2008; Hirono et al. , 2006, 2007) Production of CO 2 attested important Experimental observations (Han et al. , 2007) Additional mechanism to be investigated among possible fault weakening processes.

Evidence of CO 2 release in several active crustal faults • Corinth rift, Greece. Chemical analyses of water springs near the seismogenic Heliki and Aegion faults revealed an anomalously high content of dissolved CO 2 compared with the regional values (Pizzino et al. , 2004) • San Andreas Fault. The surface trace of the SAF displays a positive anomaly of CO 2 fluxes (Lewicki and Brantley, 2000) and this CO 2 comes from a shallow source, not from the mantle (Lewicki et al. , 2003). • Nojima Fault, Japan. Kobe earthquake, 1995. A microinfrared analysis of exhumed pseudotachylites (i. e. friction induced melts produced by seismic slip) from the Nojima fault revealed a carbon supersaturation in the melts, and the quantity of CO 2 released by friction melting during the 1995 Kobe earthquake was evaluated to 1. 8 to 3. 4 10 3 tons (Famin et al. , 2008). • Central Apennines. Italiano et al. , 2008 reported enhanced fluxes of crustal CO 2 (i. e. not mantellic) during the 1997 -1998 seismic crisis of major faults, and proposed that coseismic decarbonation was responsible for the CO 2 emission. • Chelungpu Fault, Taiwan. Chi earthquake. The fault core was depleted in carbon relative to the damage zone, and the deplation was attributed to a decarbonation induced by frictional heat (Hirono et al. , 2006, 2007) • High velocity friction experiments on Carrara marble have shown that thermal decomposition of calcite due to frictional heating induces a pronounced fault weakening (Han et al. , 2007).

Evidence of CO 2 release in several active crustal faults Field observations Experimental observations Fault Location Reference Heliki and Aegion Corinth rift, Greece Pizzino et al. , 2004 San Andreas California Lewicki and Brantley, 2000; Lewicki et al. , 2003) Nojima Japan Famin et al. , 2008 Central Appenienes Italy Italiano et al. , 2008 Chelungpu Fault Taiwan Hirono et al. , 2006, 2007 High velocity friction experiments Carrara marble Han et al. , 2007

Deforming shear band with heat and fluid fluxes Governing equations Energy Balance Mass Balance

Evolution of temperature, pore pressure and shear stress Decomposition of carbonate rock begins at ~ 700 °C Temperature increase drastically slowed due to the energy consumed in the endothermic chemical reaction 1 2 fault weakening Also evolution of T and p with slip, if considering constant slip rate V = 1 m/s (slip is: d = V t). fault strengthening Two competing effects act on the evolution of the pore pressure: 1) CO 2 production induces an additional fluid mass and thus a pressurization of the pore fluid 2) Increase of porosity due to the solid decomposition induces an increase of the permeability of the medium which limits the pressurization Shear stress decreases rapidly during initial pressurization and then increases again Mineral decomposition of the rock is a mechanism of fault weakening in a first stage and of fault strengthening in a second stage.

Correction to J. Sulem & V. Famin, JGR, 2009 Account for mass of Ca. O produced in the chemical reaction of calcite decomposition (decarbonation) in the mass balance of the solid phase

Influence of some parameters • Shear band thickness • Friction coefficient • Slip velocity • Initial permeability • Initial state of stress and initial temperature

How much can these models explain ? Rice, JGR, 2006 • In absence of a strong weakening mechanism, temperature rise would lead to widespread melting, yet evidence of glass (pseudotachylyte) that would be left from rapid recooling is not pervasive on most exhumed faults. • Relevant weakening processes in large crustal events are therefore likely to be thermal. Two possible weakening processes are: 1) thermal pressurization of pore fluid 2) flash heating at highly stressed frictional microcontacts • For sufficiently large slip, a macroscopic melt layer may form too at high enough normal stress • In silica-rich lithologies, weakening by gel formation may occur instead. Sulem & Famin, JGR, 2006 • Decarbonation is a source of CO 2 that significantly increases the slip-weakening effect of thermal pressurization • An important result of this model is that the endothermic reaction of calcite decomposition limits the coseismic temperature increase to less than ~ 800°C within a carbonate shear band under rapid slip • Large earthquakes are more prone to temperature limitation • This model provides another explanation to the anomalous absence of positive heat flow on active crustal faults, such as San Andreas: a large part of the heat produced by friction would be consumed by endothermic reactions • Friction melting hampered by endothermic calcite decomposition in carbonate fault rocks (and probably in other faults containing a sufficiently high proportion of volatile-rich secondary minerals) • This is consistent with relative scarcity of pseudotachylytes in mature faults such as the Punchbowl fault, and their occurrence in less evolved faults such as Nojima.