Groundwater Interacting at Depth With Hot Plastic Magma Triggers Phreatic Eruptions at Yugama Crater Lake of Kusatsu-Shirane Volcano (Japan)

Interpreting the triggering mechanisms for phreatic eruptions is a key to improving the hazard assessment of crater lakes. Yugama Crater Lake at Kusatsu-Shirane volcano, Japan, is the site of frequent phreatic eruptions with the recent eruptions in 1982–83, 1989, and 1996, as well as volcanic unrest, including earthquake swarms in 2014 and 2018. To understand the magma–hydrothermal interaction beneath Yugama Crater Lake, we analyzed lake waters from November 2005 to May 2021. From 2005 to 2012, Cl and SO4 concentrations decreased slowly, suggesting the development of a self-sealing zone surrounding the crystallizing magma. We focused on Ca, Al, and Si concentrations as representatives of the breach and dissolution of minerals comprising the self-sealing zone and the Mg/Cl ratio as an indicator for enhanced interaction between groundwater and hot plastic rock within the self-sealing zone. In 2006–2007, the Ca, Al, Si concentrations and the Mg/Cl ratio increased. No Cl and SO4 increase during this period suggests the self-sealing zone was leached by deep circulating groundwater rather than by magmatic fluids injection. After the 2014 earthquakes, Ca, Al, and Si increased again but were associated with a significant Cl increase and a pH decrease. We believe that the HCl-rich magmatic fluids breached the self-sealing zone, leading to fluids injection from the crystallizing magma to the Yugama crater. During this period, the Mg/Cl ratio did not increase, meaning that magmatic fluids ascending from the breached area of the self-sealing zone inhibited deep intrusion of groundwater into the hot plastic rock region. In 2018, magmatic fluids ascended through the self-sealing zone again with less intensity than in 2014. All eruptions since 1982 have been accompanied by a Mg/Cl ratio increase and a Cl decrease, whereas, when a significant HCl input occurs, as in 2014, no eruptions and no Mg/Cl ratio increase occurred. This demonstrates that the groundwater–hot plastic rock interaction, rather than the magmatic fluids input, played an essential role in triggering phreatic eruptions; i.e., phreatic eruptions can potentially occur without clear signs of fresh magma intrusions.

Effect of Joint Density on Rockburst Proneness of the Elastic-Brittle-Plastic Rock Mass

The prediction of rockburst proneness is the basis of preventing and controlling rockburst disasters in rock engineering. Based on energy theory and damage mechanics, the quantitative functional relationship between joint density and energy density was derived. Then, the theoretical results were verified by numerical simulation and uniaxial compression test, and the effect of joint density on rockburst proneness of the elastic-brittle-plastic rock mass was discussed. The results show that the relationship between the joint density and the dissipated energy index of the jointed rock mass is a logarithmic function. With the same total input energy, the higher the joint density, the more the damage dissipation energy. Even in the case of high joint density, the rock mass still has limited resistance to external failure. Under the same joint density, the strength of parallel jointed rock mass is better than that of the cross-jointed rock mass, and the parallel jointed rock mass can accumulate more elastic strain energy and has higher rockburst proneness. The joint density is closely related to the ability of the rock mass to store high strain energy. The higher the joint density is, the weaker the ability to accumulate the elastic strain energy of rock mass is and the lower the rockburst proneness is. It is helpful to predict rockburst proneness by investigating and studying the properties of geological discontinuities. The research results have some theoretical and engineering guiding significance for the prediction of rockburst proneness of the jointed rock mass.

The spatial–temporal total friction coefficient of the fault viewed from the perspective of seismo-electromagnetic theory

Recently, it has been shown theoretically how the lithospheric stress changes could be linked with magnetic anomalies, frequencies, spatial distribution and the magnetic-moment magnitude relation using the electrification of microfractures in the semibrittle–plastic rock regime (Venegas-Aravena et al., 2019). However, this seismo-electromagnetic theory has not been connected with the fault's properties in order to be linked with the onset of the seismic rupture process itself. In this work we provide a simple theoretical approach to two of the key parameters for seismic ruptures which are the friction coefficient and the stress drop. We use sigmoidal functions to model the stress changes in the nonelastic regime within the lithosphere. We determine the temporal changes in frictional properties of faults. We also use a long-term friction coefficient approximation that depends on the fault dip angle and four additional parameters that weigh the first and second stress derivative, the spatial distribution of the nonconstant stress changes, and the stress drop. We found that the friction coefficient is not constant in time and evolves prior to and after the earthquake occurrence regardless of the (nonzero) weight used. When we use a dip angle close to 30∘ and the contribution of the second derivative is more significant than that of the first derivative, the friction coefficient increases prior to the earthquake. During the earthquake event the friction drops. Finally, the friction coefficient increases and decreases again after the earthquake occurrence. It is important to mention that, when there is no contribution of stress changes in the semibrittle–plastic regime, no changes are expected in the friction coefficient.

The spatial-temporal total friction coefficient of the fault viewed from the seismo-electromagnetic theory

Recently, it has been shown theoretically how the lithospheric stress changes could be linked with magnetic anomalies, frequencies, spatial distribution and the magnetic-moment magnitude relation using the electrification of microfractures in the semi brittle-plastic rock regimen [Venegas-Aravena et al. Nat. Hazards Earth Syst. Sci. 19, 1639–1651 (2019)]. However, this Seismo-electromagnetic Theory still has not shown any relation, approach or changes in the fault's properties in order to be linked with the beginning of seismic rupture process itself. In this work we show the first and simple theoretical approach to one of the key parameters for seismic ruptures as is the friction coefficient and the stress drop. We use sigmoidal stress changes in the non-elastic regimen within lithosphere described before to figure out the temporal changes in frictional properties of faults. We also use a long term friction coefficient approximation that can depend on the fault dip angle, four parameters that weight the first and second stress derivative, the spatial distribution of the non-constant stress changes and the stress drop. It is found that the friction coefficient is not constant in time and evolve previous and after the earthquake occurs regardless of the (non-zero) weight used. When we use a dip angle close to 30 degrees and the contribution of the second derivative is more significant than the first derivative, the friction coefficient increase previous the earthquake. Then, the earthquake occurs and the friction drop. Finally, the friction coefficient increases and decreases after the earthquake. When there is no contribution of stress changes in the semi brittle-plastic regimen, no changes are expected in the friction coefficient.