Earthquake rupture properties in presence of thermal -pressurization of pore fluids

<p>Frictional heat generated in the fault core during earthquake rupture can raise the fluid pressure in the slip zone. Such increase of fluid pressure decreases the effective normal stress and thereby lowers the frictional strength of the fault. Therefore, thermal pressurization (TP) of pore fluid affects earthquake rupture processes including nucleation, propagation, and arrest. While the effects of pore pressure and fluid flow rate on dynamic weakening of faults are qualitatively understood, a detailed analysis of how TP affects  earthquake rupture parameters is needed to further deepen our understanding. </p><p>In this study, we investigate the role of two key TP parameters -- hydraulic diffusivity and shear-zone half-width -- earthquake dynamics and kinematic source properties (slip, peak slip-rate, rupture speed and rise time). We conduct  a suite of 3D dynamic rupture simulations applying a rate-and-state dependent friction law (with strong velocity weakening) coupled with thermal-pressurization of pore fluids. Simulations are carried out with the open source software SeisSol (www.seissol.org). The temporal evolution of rupture parameters over ~1’000 randomly  distributed on-fault receivers is statistically analyzed in terms of  mean variations of rupture parameters and correlations among rupture parameters. </p><p>Our simulations reveal that mean slip decreases with increasing hydraulic diffusivity, whereas mean peak slip-rate and rupture speed remain nearly constant. On the other hand, we observe only a slight decrease of mean slip with increasing shear-zone half-width, whereas mean peak slip-rate and rupture speed show clear decrease. The faster diffusion of pore pressure as hydraulic diffusivity increases promotes faster increase of the effective normal stress (and fault strength) behind the main rupture front, reducing the rise time and, therefore also affecting mean slip. An increase in shear-zone half- width represents a heat source distributed over larger fault normal distance causing a second-order effect on mean slip. Additionally, our simulations reveal correlations among rupture parameters: 1) slip has weak negative correlation with peak slip-rate and negligible correlation with rupture speed, but a positive correlation with rise time, 2) peak slip-rate has a strong positive correlation with rupture speed, but a strong negative correlation with rise time, 3) rupture speed has strong negative correlation with rise time. We observe little or negligible effects of variations of hydraulic diffusivity and shear-zone half- width on the correlations between rupture parameters. Overall, our study builds a fundamental understanding on how thermal pressurization of pore fluids affects dynamic and thereby kinematic earthquake rupture properties. Our findings are thus important for the earthquake source modeling community, and particularly, for assessing seismic hazard due to induced events in geo-reservoirs.</p>

Experimental Evidence and Characteristic Recognition of the Nanoweakening of Slip Deformation Zones

A central issue in the study of fault evolution is to identify shear weakening and its mechanism; currently, studies of fault weakening in narrow slip deformation zones, including those of various slipping planes such as schistosity, foliation, cleavage, joints and faults in rocks, are ongoing. To verify the nanoweakening in shear slipping, we carried out experiments: triaxial compression experiments on sandstones and uniaxial compression experiments on granites. Furthermore, on the basis of scanning electron microscopy (SEM) observations and experimental data analyses, we suggested three kinds of nanoweakening in terms of the corresponding strain stages: (1) The slip nanoweakening caused by the strain hardening deformation stage of the shear slip, which creates nanograins with dense coatings that may be due to the nanocoating on the shear planes, can result in rolling friction rather than with sliding friction, and the former is a principal mechanism of sliding nanoweakening. (2) The rheological nanoweakening caused by the strain softening deformation stage; in view of developing weakened deformation due to grain boundary migration (GBM), the flow of synkinematic minerals and melt coating phenomena lead to rheological nanoweakening. (3) The dynamic nanoweakening caused by thermal pressurization and fluid pressurization during the strain softening stage and strain degenerating stage. Thus, when these aspects are considered in defining the relationship between the nanoweakening at the slipping planes and the strain stages, the representative mechanism and its behavior rules can be obtained.

Rupture-dependent breakdown energy in fault models with thermo-hydro-mechanical processes

Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need to reexamine the assumption of predetermined frictional breakdown energy common in dynamic rupture modeling and to better understand the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes.