Physical and mechanical rock properties of a heterogeneous volcano; the case of Mount Unzen, Japan
Authors: Kendrick Jackie E. , Schaefer Lauren N. , Schauroth Jenny , Bell Andrew F. , Lamb Oliver D. , Lamur Anthony , Miwa Takahiro , Coats Rebecca , Lavallée Yan , Kennedy Ben M. 
Published: October 2020

Abstract. Volcanoes represent one of the most critical geological settings for hazard modelling due to their propensity to both unpredictably erupt and collapse, even in times of quiescence. Volcanoes are heterogeneous at multiple scales, from porosity which is variably distributed and frequently anisotropic to strata that are laterally discontinuous and commonly pierced by fractures and faults. Due to variable and, at times, intense stress and strain conditions during and post-emplacement, volcanic rocks span an exceptionally wide range of physical and mechanical properties. Understanding the constituent materials' attributes is key to improving the interpretation of hazards posed by the diverse array of volcanic complexes. Here, we examine the spectrum of physical and mechanical properties presented by a single dome-forming eruption at a dacitic volcano, Mount Unzen (Japan) by testing a number of isotropic and anisotropic lavas in tension and compression and using monitored acoustic emission (AE) analysis. The lava dome was erupted as a series of 13 lobes between 1991–1995, and its ongoing instability means much of the volcano and its surroundings remain within an exclusion zone today. During a field campaign in 2015, we selected 4 representative blocks as the focus of this study. The core samples from each block span range in porosity from 9.14 to 42.81 %, and permeability ranges from 1.54 × 10−14 to 2.67 × 10−10 m2 (from 1065 measurements). For a given porosity, sample permeability varies by > 2 orders of magnitude is lower for macroscopically anisotropic samples than isotropic samples of similar porosity. An additional 379 permeability measurements on planar block surfaces ranged from 1.90 × 10−15 to 2.58 × 10−12 m2, with a single block having higher standard deviation and coefficient of variation than a single core. Permeability under confined conditions showed that the lowest permeability samples, whose porosity largely comprises microfractures, are most sensitive to effective pressure. The permeability measurements highlight the importance of both scale and confinement conditions in the description of permeability. The uniaxial compressive strength (UCS) ranges from 13.48 to 47.80 MPa, and tensile strength (UTS) using the Brazilian disc method ranges from 1.30 to 3.70 MPa, with crack-dominated lavas being weaker than vesicle-dominated materials of equivalent porosity. UCS is lower in saturated conditions, whilst the impact of saturation on UTS is variable. UCS is between 6.8 and 17.3 times higher than UTS, with anisotropic samples forming each end member. The Young's modulus of dry samples ranges from 4.49 to 21.59 GPa and is systematically reduced in water-saturated tests. The interrelation of porosity, UCS, UTS and Young's modulus was modelled with good replication of the data. Acceleration of monitored acoustic emission (AE) rates during deformation was assessed by fitting Poisson point process models in a Bayesian framework. An exponential acceleration model closely replicated the tensile strength tests, whilst compressive tests tended to have relatively high early rates of AEs, suggesting failure forecast may be more accurate in tensile regimes, though with shorter warning times. The Gutenberg-Richter b-value has a negative correlation with connected porosity for both UCS and UTS tests which we attribute to different stress intensities caused by differing pore networks. b-value is higher for UTS than UCS, and typically decreases (positive Δb) during tests, with the exception of cataclastic samples in compression. Δb correlates positively with connected porosity in compression, and negatively in tension. Δb using a fixed sampling length may be a more useful metric for monitoring changes in activity at volcanoes than b-value with an arbitrary starting point. Using coda wave interferometry (CWI) we identify velocity reductions during mechanical testing in compression and tension, the magnitude of which is greater in more porous samples in UTS but independent of porosity in UCS, and which scales to both b-value and Δb. Yet, saturation obscures velocity changes caused by evolving material properties, which could mask damage accrual or source migration in water-rich environments such as volcanoes. The results of this study highlight that heterogeneity and anisotropy within a single system not only add uncertainty but also have a defining role in the channelling of fluid flow and localisation of strain that dictate a volcano's hazards and the geophysical indicators we use to interpret them.