Geomechanical rock properties of a basaltic volcano

Abstract

In volcanic regions, reliable estimates of mechanical properties for specific volcanic events such as cyclic inflation-deflation cycles by magmatic intrusions, thermal stressing, and high temperatures are crucial for building accurate models of volcanic phenomena. This study focuses on the challenge of characterizing volcanic materials for the numerical analyses of such events. To do this, we evaluated the physical (porosity, permeability) and mechanical (strength) properties of basaltic rocks at Pacaya Volcano (Guatemala) through a variety of laboratory experiments, including: room temperature, high temperature (935°C), and cyclically-loaded uniaxial compressive strength tests on as-collected and thermally-treated rock samples. Knowledge of the material response to such varied stressing conditions is necessary to analyze potential hazards at Pacaya, whose persistent activity has led to 13 evacuations of towns near the volcano since 1987. The rocks show a non-linear relationship between permeability and porosity, which relates to the importance of the crack network connecting the vesicles in these rocks. Here we show that strength not only decreases with porosity and permeability, but also with prolonged stressing (i.e., at lower strain rates) and upon cooling. Complimentary tests in which cyclic episodes of thermal or load stressing showed no systematic weakening of the material on the scale of our experiments. Most importantly, we show the extremely heterogeneous nature of volcanic edifices that arise from differences in porosity and permeability of the local lithologies, the limited lateral extent of lava flows, and the scars of previous collapse events. Input of these process-specific rock behaviors into slope stability and deformation models can change the resultant hazard analysis. We anticipate that an increased parameterization of rock properties will improve mitigation power.

Introduction

Several studies of non-volcanic slopes using numerical models and geomechanical material properties have been adapted for the investigation of volcanic collapses, slope stability, and ground deformation. These include applying Limit Equilibrium Methods (LEM) (Voight et al., 1983; Iverson, 1995; Elsworth and Voight, 1996; Hürlimann et al., 2000; Reid et al., 2000; Donnadieu et al., 2001; Okubo, 2004; Apuani and Corazzato, 2010) and numerical modeling by Distinct Element Methods (DEM) (Zimbelman et al., 2004), Finite Difference Methods (FDM) (Hürlimann et al., 2000; Apuani et al., 2005a,b, 2007; Apuani and Corazzato, 2010) and Finite Element Methods (FEM) (Sousa and Voight, 1995; Apuani and Corazzato, 2010; del Potro et al., 2013; Schaefer et al., 2013). However, numerical data on the relevant mechanical properties for volcanic rocks and rock masses used in these types of studies remain a major source of uncertainty. Quantifying these properties is challenging due to the diversity of materials that compose volcanoes, the unknown interior structures and geometry, and the often inaccessible locations of volcanic rock masses (Thomas et al., 2004; del Potro and Hürlimann, 2008).

Experimental rock mechanics has traditionally been executed in the fields of civil engineering and mining as well as structural geology, with a focus on sedimentary and intrusive igneous rock (see Paterson and Wong, 2005 and references therein). Early efforts were made to understand rock behavior and strength under high pressure (Scholz, 1968; Tapponnier and Brace, 1976; Batzle et al., 1980), high temperature (Griggs et al., 1960; Simmons and Cooper, 1978; Handin and Carter, 1979; Fredrich and Wong, 1986), and in water saturated conditions (Bauer et al., 1981). Additional work was done to connect mineralogy and microstructures with rock strength (Merriam et al., 1970; Přikryl, 2001) and to model the mechanical properties of porous material; that is, the pore-emanated crack model (Sammis and Ashby, 1986) and the wing-crack model (Ashby and Sammis, 1990). Only recently have laboratory studies begun to systematically investigate the microstructural, physical, and mechanical properties of volcanic rocks, a material infamously known for its formation under disequilibrium conditions, thereby rich in heterogeneities at all scales. These include experiments on a range of material under relevant volcanic conditions such as thermal stressing (Vinciguerra et al., 2005; Kendrick et al., 2013a; Heap et al., 2014a), cyclic inflation-deflation cycles by intrusions (Heap et al., 2009, 2010; Kendrick et al., 2013a), fragmentation (Spieler et al., 2004; Kueppers et al., 2006; Scheu et al., 2008), and flow or fracture at high temperatures and/or pressures (Balme et al., 2004; Rocchi et al., 2004; Smith et al., 2005; Lavallée et al., 2007, 2008; Benson et al., 2008; Cordonnier et al., 2009; Loaiza et al., 2012; Kendrick et al., 2013b).

These investigations have shown that porosity is a fundamental control on rock properties. At volcanoes, the porous network comprises of vesicles and fractures that build the permeable network and which controls degassing (Ashwell et al., 2015). Porosity also controls the mechanical response of materials whereby their strength decreases and elasticity increases with porosity (e.g., Al-Harthi et al., 1999; Spieler et al., 2004; Paterson and Wong, 2005; Heap et al., 2014a,b). Thermal stressing of rock has also been argued to affect its strength, although damage imparted by thermal cracking has been shown to lead to unsystematic changes in mechanical response (Jones et al., 1997; Rocchi et al., 2003; Balme et al., 2004; Kendrick et al., 2013a; Heap et al., 2014a); the ambiguity of this claim has been linked to the extent of pre-existing microcrack damage in a rock (Vinciguerra et al., 2005; Heap et al., 2014a). Other rocks can debilitate upon thermal stressing, as is the case during decarbonation (Heap et al., 2013b) and dehydroxylization (Heap et al., 2012) of certain rock-forming minerals. Cyclic loading can lead to mechanical weakening (Kendrick et al., 2013a) and an increase in heterogeneity of elastic modulus values, likely due to an increase in anisotropic crack damage within the samples (Heap and Faulkner, 2008; Heap et al., 2010).

This paper focuses on the challenge of characterizing volcanic materials for numerical analyses and addressing the scarcity of published geomechanical data for volcanic materials. To do this, we evaluated the mechanical properties of rocks at Pacaya Volcano (Guatemala) using a variety of laboratory experiments, including room temperature, high temperature, and cyclically-loaded uniaxial compressive strength tests (UCS) on as-collected and thermally treated rock samples. Pacaya is an ideal case study for material characterization for use in deformation analysis due to its persistent activity, significant morphological changes, and history of slope instability (Vallance et al., 1995; Matías Gómez et al., 2012; Schaefer et al., 2013). By considering the non-unique geomechanical behaviors when the materials are subjected to this wide range of stressing conditions, we may allow the determination of more encompassing and realistic models that will help better predict volcanic processes and imminent catastrophes.

Geologic setting

Pacaya is a basaltic stratovolcano complex in the Central American Volcanic Arc (Figure 1A) that has been persistently active since 1961 with over 250 lava flows, intermittent strombolian activity, and ash and fumarolic plumes (Matías Gómez et al., 2012; Rose et al., 2013). Between 0.6 and 1.5 ka B.P., the SW sector of the initial cone failed in a major edifice collapse, forming a debris avalanche that traveled 25 km SW, containing approximately 0.65 km³ of debris (Kitamura and Matías Gomez, 1995; Vallance et al., 1995). The edifice re-built within the collapse amphitheater until the nineteenth century, when it went into a period of repose. Activity renewed once again in 1961, loading the SW flank of the edifice with 100–150 m thick of primarily effusive lava flows. Although the volcano has grown substantially, the collapse scarp is still exposed today on the north and east sides of the cone on account of the asymmetric loading of material on the SW flank.

Figure 1

Explosive and effusive eruptions have caused constant morphometric changes to the edifice: in 1962, an oval-shaped area, approximately 300 m by 200 m, subsided near the summit. This subsidence has since been in-filled by erupted material. In 2010, a large collapse trough 50–80 m deep and 100 m wide, extending 600 m from the summit, developed concurrently with explosive activity on May 27, 2010 (Figure 1B). Although the exact causes of these collapses are unknown, shallow magma systems and high growth rates leading to over-steepening of the cone have been postulated as factors affecting edifice stability (Eggers, 1983; Vallance et al., 1995; Matías Gómez, 2009; Schaefer et al., 2013). The rock characterization determined in this study can help to better understand events such as these and assess potential future intrusions, collapses, and deformation events.

Materials and experimental methodology

Experimental materials

Sample characterization

Pacaya volcano is made of diverse lithologies (lavas and volcaniclastic deposits) which vary in extent and character significantly, though non-systematically (Wunderman and Rose, 1984; Matías Gómez et al., 2012). This heterogeneity at the large scales calls for a characterization of equally diverse rock types. For this experimental study, eight sample blocks were collected from different eruptive deposits (Figure 1C, see Supplementary Image for detailed volcanological map from Matías Gómez et al., 2012), which cover a range of macroscopically representative physical properties of the edifice-forming basalt at Pacaya (Table 1). These samples, which vary in porosity, crystal size, and vesicle size distribution include: two samples north of the summit from lava flows in January and August–October of 2013; two samples of blocks from rock fall collected from on top of the 2010 deposits; two samples from lava flows associated with the May-June eruption in 2010; one sample from the brecciated top of a lava flow that occurred in February of 1974; and one sample from the (>1.5 ka B.P.) southern portion of the ancestral scarp (Figures 2A–H). Typically, basalt at Pacaya is porphyritic with olivine and plagioclase phenocrysts (45 wt%), with lesser amounts of clinopyroxene and iron-titanium oxides and plagioclase dominated microlites (Bardintzeff and Deniel, 1992). No significant petrographic differences have been found between prehistoric and modern lava flows. Each of the sample blocks was drilled to prepare a series of cylindrical cores with a diameter of 26 mm and a length of 52 mm to maintain a length-to-diameter ratio of 2:1. Although other samples were collected in the field (see Supplementary Image), the eight sample blocks were chosen from which cores could be drilled to provide the widest range of representative samples while ensuring reliable and repeatable measurements of strength. Based on the abundance of effusive products, and the similarity of eruptive product composition over time, we feel as though our samples represent the majority of edifice-building material variability.

Figure 2

Porosity and permeability

Before each test, porosity was determined for each core using a 100cc helium pycnometer from Micromeritics, providing measurements with an accuracy of 0.1%. First, the density of a cylindrical core (ρ_rock) was estimated and used to calculate the connected porosity (ϕ_c). Total porosity (ϕ_T) was also determined for each sample type by measuring the density of a powdered sample (ρ_powder), and computing:

From this, isolated porosity was calculated by subtracting connected porosity from total porosity.

Permeability was determined for each sample block using a TinyPerm II Permeameter from New England Research Inc., which was also used to measure the permeability of 97 in-situ lavas and breccias (unconsolidated air-fall tephras to brecciated surfaces of lava flows) at Pacaya during fieldwork in October 2013.

Thermal stressing

Thermal stressing occurs in volcanic edifices due to heat radiating from the magma in the conduit that may be upwards of 1100°C. Rocks in direct contact with conduit magma may experience partial melting, and further afield temperatures exceeding 900°C are likely common. Thermal analysis was performed on cores of samples in a box furnace by heating them to 1050°C and then re-cooling to room temperature at a rate of 10°C/min to assess physical changes, including oxidation and the onset of melting. This data was used to constrain the maximum experimental temperature for uniaxial compressive strength tests, as it revealed partial melting between 1000 and 1050°C for some samples; therefore a temperature of 935°C was chosen for the high-temperature tests as well as for the thermal stressing tests.

For thermal stressing, three cores from each block were placed in a box furnace and heated at a rate of 4°C/min to the target temperature, left to dwell for 1 h to allow thermal equilibration and completion of any reactions which might have taken place, and then cooled to room temperature using the same rate. A rate of 4°C/min was chosen to minimize damage that could result from temperature gradients within the sample with the presumption that heating in volcanic scenarios occurs relatively steadily, except for rocks in direct contact with magma. Porosity of samples was measured before and after thermal treatment. Thin sections were made from each block with and without thermal stressing (Figures 2A–P).

Experimental methods

Uniaxial compressive strength (UCS) experiments were conducted at constant strain rates in a 100 kN Instron uniaxial loading frame (model 1362) in the Experimental Volcanology Laboratory at the University of Liverpool. A summary of the experiment type, temperature, and strain rate is shown in Table 2. Applied load (in kN) and displacement (mm) were measured throughout at a rate of 10–100 Hz. Following calibration, to take into account the compliance of the experimental set-up at the given conditions, axial stress (σ_a; MPa) and strain (ε_a) were calculated. The Static Young's modulus (E_s) was calculated from the linear portion of the stress/strain curve (at 30% of the peak stress) via:

Room temperature experiments

UCS experiments were performed on dry specimens at room temperature by testing three cores (52 mm length and 26 mm diameter) from each of the eight original sample blocks as well as three thermally stressed cores from seven sample blocks. The cylindrical cores were axially loaded at strain rates of 10⁻¹, 10⁻³, and 10⁻⁵ s⁻¹ (5, 0.05, and 0.0005 mm/s, respectively) until failure occurred, as defined by a sharp stress drop.

Cyclic experiments

In cyclic experiments, the uniaxial compressive stress was applied at a constant strain rate of 10⁻⁵ s⁻¹ until 15% of the sample's value of UCS (as determined from previous room temperature UCS tests at a strain rate of 10⁻⁵ s⁻¹), then unloaded to 5% of this stress. Subsequently, a further loading/unloading cycle was imposed, but this time reaching an additional 15% load increment (i.e., 30% of UCS), then again unloaded to 5%. These cycles were repeated by increasing the peak stress in 15% stress increments of the UCS until failure. All cyclic tests were conducted at room temperature on dry, as-collected specimens at a strain rate of 10⁻⁵ s⁻¹.

High temperature experiments

High temperature tests were conducted by placing a core in the same uniaxial press surrounded by a three-zone split cylinder furnace, and heated at a rate of 4°C/min until a temperature of 935 ± 3°C was reached. Before undertaking the experimental program, the compliance of the setup at high temperature was measured and considered to remove any of the mechanical artifacts of the machine on the dataset. Samples were held to thermally equilibrate at this temperature for an hour before the experiment began. Temperature was consistently monitored during heating, cooling, and while the compressive load was applied by a thermocouple in contact with the outer surface of the specimen. Samples at high temperature were tested at strain rates of 10⁻¹ and 10⁻⁵ s⁻¹ to determine end-member strain-rate behavior.

Results

Porosity and permeability

Samples range from 2 to 48% total porosity, and isolated pores makes up a varying but sometimes substantial amount of total porosity (Figure 3) as a result of the varying abundance of microcracks and vesicles (Figure 2). Although porosity varied greatly between samples, it remained relatively constant between cores from the same sample, with standard deviations of <3.3% for five or more measurements. Pre- and post-thermal stressing pycnometer measurements showed only minor changes in porosity of −0.4 to +0.2% despite visible differences in sample structure that resulted from thermal cracking and/or oxidation (Figure 2). Permeability measurements showed a bi-modal distribution (Figure 4A) that approximately split into coherent lavas which were sampled for this study (peak at 10⁻¹³), and brecciated lava flow surfaces and unconsolidated air-fall deposits (peak at 10⁻¹⁰). The sample's permeability was observed to scale to porosity via the relationship for effusive lavas constrained by Mueller et al. (2005) (Figure 4B).

Figure 3

Figure 4

Uniaxial compressive strength results

Here we present the results of 64 UCS experiments and the effects of (1) strain rate, (2) porosity, (3) cyclic loading, (4) thermal stressing, and (5) temperature on the strength of the tested rocks. For all UCS stress-strain curves, see Supplementary Material.

A comparison of UCS stress-strain curves for all as-collected, room temperature cores at 10⁻¹ and 10⁻⁵ s⁻¹ strain rates is shown in Figure 5. These curves display the shape of typical brittle failure under compression, outlined in Heap et al. (2014a) as follows: (1) convex shape as microcracks aligned sub-perpendicular to the direction of loading begin to close (leading to a large strain despite little stress accumulation); (2) linear shape as the rock deforms elastically; (3) concave shape as the rock deforms inelastically and experiences strain hardening; and (4) rock failure marked by a substantial stress drop. Higher strain rate tests fail in a more brittle manner (marked by sharper curve peaks at the stress drop) than slower strain rates (Figure 5) and reach both a higher axial strain and a higher peak stress before failure occurs (Figure 6A). Indeed, room temperature, thermally stressed, and high temperature tests all show a clear pattern of apparent strengthening with increasing strain rates, which is more significant at lower porosity (Figure 6B, Table 3). On average, room temperature cores tested at a strain rate of 10⁻¹ s⁻¹ failed at axial stresses 5% higher than strain rates of 10⁻³ s⁻¹, and strain rates of 10⁻³ s⁻¹ failed 29% higher than strain rates of 10⁻⁵ s⁻¹. In a similar pattern, thermally stressed cores tested at a strain rate of 10⁻¹ s⁻¹ failed at axial stresses 4% higher than strain rates of 10⁻³ s⁻¹, and strain rates of 10⁻³ s⁻¹ failed at a stress 22% higher than the stress at strain rates of 10⁻⁵ s⁻¹. Finally, high-temperature cores failed at an average 52% higher axial stress at 10⁻¹ s⁻¹ as compared to strain rates of 10⁻⁵ s⁻¹ (Table 3).

Figure 5

Figure 6

The UCS results show that failure stress is highly dependent on porosity, with maximum stress decreasing non-linearly as porosity increases (Figure 6B). Samples can be divided into three porosity groups, with samples having <10% porosity averaging a UCS of 138 MPa, samples with 10–25% porosity with a strength decreasing by 55% to an average of 62 MPa, and samples with >25% porosity with a strength decreasing another 66% to an average of 21 MPa. Overall, axial compressive strength for Pacaya basalt ranges widely from 12 to 174 MPa (Table 3), with the lower end falling into strength categories as low as chalk (Hoek, 1994). Given that the porosity of the rocks scales with permeability, we chose to assess the relationship of UCS with permeability, as this measures the strong control that cracks have on the physical properties of rocks. The data shows a remarkable linearity in UCS decrease with logarithmic increase in rock permeability (Figure 6C). Cyclic loading tests in which the peak load of each successive cycle was increased had a variable effect on strength (Table 3, Figure 6). Each cycle peak was 15% of the failure strength of the room temperature tests higher (as-collected cores at 10⁻⁵ s⁻¹). The highest stress was usually reached on the final cycle, exceeding all previous stress cycles (Figures 7A,B). However, some cyclically loaded samples failed at a lower stress than the peak stress experienced in the penultimate cycle (Figures 7C,D). Compared to the mechanical behavior of the rocks during standard UCS tests at strain rates of 10⁻⁵ s⁻¹, four out of seven rocks failed at a lower stress during the cyclic tests, and five out of seven rocks failed at lower strain (Figure 8A). Continuous stress-strain curves (Figures 7B,D) show a degree of hysteresis, showing that at a given stress, strain is higher during unloading than during loading of a particular cycle. Each successive cycle led to further accumulated (non-recoverable) strain noted when the samples were unloaded to 5% of their value of their UCS (Figure 8B).

Figure 7

Figure 8

The effect of thermal stressing of cores varied across our investigated suite of samples: 12 out of 20 cores failed at a lower stress than the original samples (Table 3). The difference in porosity before and after thermal treatment was, however, minimal (−0.4 to +0.2%) and overall, the changes in strength by thermal stressing are within the natural variability of each sample (Figures 6A,B). Additionally, there were no significant variations in the magnitude of strain to failure or style of failure upon thermal stressing (see Supplementary Information).

The deformation of rock at high-temperature exhibited an increase in strength. On average, strength at 935°C was 19 and 5% higher than the rock strength at room temperature for deformation at strain rates of 10⁻¹ and 10⁻⁵ s⁻¹, respectively. These tests revealed that rock failure at high temperature requires a higher magnitude of strain than at room temperature; this is especially the case for the lower porosity (strongest) samples (Figure 6A) The peak stress of these tests also scale inversely with porosity as do the tests at room temperature (Figure 6B).

Almost all the samples had a Young's Modulus that scaled approximately linearly to peak strength (Figure 9A), the only outliers were the high-temperature experiments, which had a lower Young's Modulus for a given compressive strength (Figure 9A). This difference becomes particularly prominent at higher strength, and corresponds to a slightly higher strain to failure, especially at low porosity (Figure 9B).

Figure 9

Discussion

Examination of the mechanical properties of eruptive products at Pacaya volcano provide us with a landscape of data from which to improve our understanding of volcano edifice properties. Of the attributes measured, porosity has the strongest control on rock strength, followed by strain rate, then temperature, which show non-negligible, second order effects. Cyclic loading and thermal stressing seem to have lesser, unsystematic influences on the mechanical behavior of the rocks tested (Figure 6B). The mechanical data show that the strength of the tested Pacaya basalts vary by slightly over one order of magnitude from 12 to 174 MPa across a range of porosities from 2 to 48%, corresponding to a permeability range of 10⁻¹²–10⁻¹⁵ m² for the experimental samples. Given that the samples tested here only represent coherent rocks forming the volcano (Figure 4A), it is likely that incoherent lithologies such as unconsolidated breccias (not studied experimentally, but for which we measured high permeabilities in the field; Figure 4A) have negligible strength, leading to presumably very weak layers which can be destabilized and jeopardize the stability of the volcano.

The wide range of strength in volcanic rocks is commonly attributed to varying compositions and microstructures inherent in volcanic rock (e.g., Byerlee, 1978; Heap et al., 2014b). Strength of materials has long been described as a function of their porosity (see Paterson and Wong, 2005). In volcanic rock, the porosity is made up of vesicles and microcracks and as such, the relationship presented between strength and permeability, though not commonly investigated, deserves consideration. Early theoretical efforts to model the mechanical properties of porous material have proposed models describing the process of failure for pore-bearing rocks, including the pore-emanated crack model (Sammis and Ashby, 1986), as well as the process of failure of crack-bearing rocks, including the wing-crack model (Ashby and Sammis, 1990). A recent micromechanical study of volcanic rock behavior suggested that the abundance of both vesicles and cracks in volcanic rocks thwart our ability to predict failure using these models (e.g., Heap et al., 2014b). Pores and cracks are fundamental to volcanic systems; the importance of permeability and its influence on gas release vs. build-up in a volcanic systems play an important role on the explosivity of an eruption (e.g., Mueller et al., 2005) but also structural stability (Thomas et al., 2004). Here, we suggest that the strength of rock and the character of the porous and permeable network are intrinsically related (Figure 4B) and influence the mechanical properties of volcanic edifices (Figures 6B,C).

Strain-rate strengthening is observed and may be a result of increased elasticity as well as crack coalescence lagging behind axial deformation at faster rates (e.g., Jaeger and Cook, 1976). Strain-rate strengthening is observed across all conditions similar to previous studies on volcanic rocks (e.g., Rocchi et al., 2004), but is more effective at high temperature, as was suggested by Lindholm et al. (1974). The strain-rate effect implies that slow strain events, such as downslope movement, creep, or shallow magma pooling, could cause more pervasive damage and be more detrimental to the stability of an edifice than fast strain events such as dike intrusions or ground motion produced during earthquakes.

High temperature strengthening occurred across the sample set. While in some cases it is possible that chemical alteration, partial melting or viscous relaxation of interstitial glass may heal pre-existing fractures, we suggest that under the conditions investigated here, it is more plausible that increased strength relates to expansion and densification (via closure of cracks) and higher elasticity in molecular bonds (see Young's Modulus values, Table 3). For example, the experiment temperature (at 935°C) is above the order-disorder transition of plagioclase, which is known to increase seismic velocities (Kono et al., 2008) and may correlate to increased elasticity. Indeed, here we show that rock deformation in a fully brittle regime at high temperatures experience significantly more strain to failure than at room temperature. This characteristic is not unique to the basalts at high temperature, though many other rock types often become weaker with temperature, especially if they undergo alteration (e.g., decarbonation, dehydroxylization; Heap et al., 2012, 2013a,b) or significant thermal expansion (e.g., due to anisotropic expansion in crystals; Meredith and Atkinson, 1982; Zhang et al., 2009).

Experiments conducted at high temperatures cause rocks to be stronger and more elastic, suggesting that areas near magma intrusions, or areas heated over prolonged periods of time, could experience more deformation before failure than previously inferred from knowledge of rock behavior at low temperatures. According to Ebmeier et al. (2013), a lack of surface deformation before eruptions at Pacaya suggests that magma rises directly from the base of the crust without a period of crustal storage. The high temperature effects may thus explain the trivial amount of pre-eruptive ground deformation associated with dike propagation and magma transport in cold rocks prior to eruption at Pacaya. When also considering the strain-rate strengthening effect, these relatively fast dike injection events could also lead to an increase in edifice elasticity, resulting in less pervasive damage and thus less detriment to the edifice forming rock. However, we must note that the results from cyclic experiments show that if dikes are continually intruded, the country rock could be mechanically weakened over time through strain accumulation, particularly in areas with high elastic modulus values.

Thermally stressing rocks has a complex effect on strength: some cores were slightly weakened, while others were slightly strengthened post-heat treatment. The variable effect of thermal stressing at similar temperatures has also been noted by changes in crack content in other rocks such as andesite (Heap et al., 2014a) and granite (Wang et al., 2013). Although heating and cooling can both induce cracking, there seems to be no systematic response as a function of starting porosity. Volcanic rocks form by cooling and as such, their microstructures and importantly here, microcracks, preserve a record of the cooling history. In our study, the cores were heated and cooled at a rate likely lower than that experienced during cooling upon eruption; it is possible that such mild thermal conditions would minimize the extent of crack damage overprint. According to Heap et al. (2009), so long as the mineralogical assemblage remains intact, microcracks (and subsequent weakening) induced by thermal stressing are likely a function of the original thermal stressing imposed by natural cooling; that is, thermal stressing must be applied at a rate that exceeds natural cooling rate to induce significant damage. Recent studies assessing the effect of fractures induced by thermal stressing concluded that its occurrence would likely not (significantly) affect volcanic edifice stability (Vinciguerra et al., 2005; Heap et al., 2009; Kendrick et al., 2013b). However, we point out that additional mechanical stresses resulting from magmatic forcing during transport would create a zone of damage, increasing the overall heterogeneity of the edifice.

Cyclic loading tests had a variable effect on failure stress and strain, however, they do suggest that repeated dike intrusion or inflation and deflation has the potential to cause edifice weakening by imparting a permanent strain on rocks. The stress-strain curves show hysteresis, and during unloading to 5% of the peak failure stress, the strain remains above zero and gradually increases from cycle to cycle (Figure 8). Such behavior has been observed in other volcanic rocks, including andesite (Kendrick et al., 2013a) and basalt (Heap et al., 2009, 2010), along with sedimentary rocks (e.g., sandstone; Song et al., 2013). Our work further suggests that material with porosities above 5 vol. % can typically absorb stress with fewer permanent damage accumulation (Figure 5B). This elastic response may result from compressibility of the rocks into the porous space, attempting to shut the cracks and vesicles. Elastic response may also be due to changes in porosity isolation as pores are connected through microcracking, with the resultant higher connected porosity allowing for more fluid stress compensation between pores, simultaneously preventing stiffening. This effect would be more abundant in rocks with higher porosity. According to Gurevich (2003), approximately 10% porosity or greater provides enough space for fluid to flow into pores and escape fracture. The implications of cyclic loading/ unloading suggest that volcanic rocks may not behave entirely elastically and accumulate damage. Once imparted, damage in turn affects the resultant material response to stress (changes in strength and Young's modulus), which may favor stress rotation (e.g., Faulkner et al., 2006), with potentially more widespread consequences on the distribution of stress in the edifice.

Conclusions

An experimental investigation of the mechanical properties of basalts at Pacaya volcano has been undertaken to increase our knowledge of contributing variables on slope stability in volcanic edifices. The rocks forming Pacaya volcano are heterogeneous at a range of scales: at the edifice scale, we find the main lithologies are locally heterogeneous due to the limited spread of lava flows as well as to flow emplacement inside earlier collapse scars; at the rock scale, we find that the rocks are macroscopically and microscopically heterogeneous, showing little variation in mineralogy, but strong variation in textures, porosities (formed of pores and cracks) and permeability. The permeability is constrained to increase non-linearly with porosity. Laboratory testing of the mechanical behavior shows that the strength decreases with porosity, increases with strain rate, increases with ambient (deformation) temperature, and is unsystematically varied by both thermal stressing, and due to prolonged stressing histories where cycles of loading/ unloading can induce mechanical hysteresis. We anticipate that integration of such process-specific rock properties in modeling of volcanic events will improve hazard analysis interpretation at active volcanoes such as Pacaya.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

• 1

  Al-HarthiA. A.Al-AmriR. M.ShehataW. M. (1999). The porosity and engineering properties of vesicular basalt in Saudi Arabia. Eng. Geo. 54, 313–320. 10.1016/S0013-7952(99)00050-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  ApuaniT.CorazzatoC. (2010). Etna flank dynamics: a sensitivity analysis by numerical modeling in ISRM International Workshop on Rock Mechanics and Geoengineering in Volcanic Environments (Puerto de la Cruz: International Society for Rock Mechanics).

  • Google Scholar
  • View reference in article

• 3

  ApuaniT.CorazzatoC.CancelliA.TibaldiA. (2005a). Physical and mechanical properties of rock masses at Stromboli: a dataset for volcano instability evaluation. Bull. Eng. Geol. Environ. 64, 419–431. 10.1007/s10064-005-0007-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  ApuaniT.CorazzatoC.CancelliA.TibaldiA. (2005b). Stability of a collapsing volcano (Stromboli-Italy): limit equilibrium analysis and numerical modelling. J. Volcanol. Geotherm. Res. 144, 191–210. 10.1016/j.jvolgeores.2004.11.028

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  ApuaniT.MerriA.MasettiM. (2007). Effects of volcanic seismic events on the Stromboli stability by finite difference numerical modeling in Proceedings of the International Workshop on Volcanic Rocks. 11th ISRM Congress (Ponta Delgada), 101–109.

  • Google Scholar
  • View reference in article

• 6

  AshbyM. F.SammisC. G. (1990). The damage mechanics of brittle solids in compression. Pure Appl. Geophys. 133, 489–521. 10.1007/BF00878002

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  AshwellP. A.KendrickJ. E.LavalléeY.KennedyB. M.HessK. U.von AulockF. W.et al. (2015). Permeability of compacting porous lavas. J. Geophys. Res. Solid. Earth120, 1605–1622. 10.1002/2014jb011519

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BalmeM. R.RocchiV.JonesC.SammondsP. R.MeredithP. G.BoonS. (2004). Fracture toughness measurements on igneous rocks using a high pressure, high temperature rock fracture mechanics cell. J. Volcanol. Geotherm. Res. 132, 159–172. 10.1016/S0377-0273(03)00343-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BardintzeffJ.DenielC. (1992). Magmatic evolution of Pacaya and Cerro Chiquito volcanological complex, Guatemala. Bull. Volcanol. 54, 267–283. 10.1007/BF00301482

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BatzleM. L.SimmonsG.SiegfriedR. W. (1980). Microcrack closure in rocks under stress: direct observation. J. Geophys. Res. Solid. Earth85, 7072–7090. 10.1029/JB085iB12p07072

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  BauerS. J.FriedmanM.HandinJ. (1981). Effects of water-saturated on strength and ductility of three igneous rocks at effective pressures to 50 MPa and temperatures to partial melting, in Rock Mechanics from Research to Application, eds EinsteinH. H.ScandriatoD. P. (Rotterdam: Balkema), 79–84.

  • Google Scholar
  • View reference in article

• 12

  BensonP. M.VinciguerraS.MeredithP. G.YoungR. P. (2008). Laboratory simulation of volcano seismicity. Science322, 249–252. 10.1126/science.1161927

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  ByerleeJ. (1978). Friction of rocks. Pure Appl. Geophys. 116, 615–626. 10.1007/BF00876528

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  CordonnierB.HessK. U.LavalléeY.DingwellD. B. (2009). Rheological properties of dome lavas: case study of Unzen volcano. Earth Planet Sci. Lett. 279, 263–272. 10.1016/j.epsl.2009.01.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  del PotroR.HürlimannM. (2008). Geotechnical classification and characterization of materials for stability analyses of large volcanic slopes. Eng. Geol. 98, 1–17. 10.1016/j.enggeo.2007.11.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  del PotroR.HürlimannM.PinkertonH. (2013). Modelling flank instabilities on stratovolcanoes: parameter sensitivity and stability analyses of Teide, Tenerife. J. Volcanol. Geotherm. Res. 256, 50–60. 10.1016/j.jvolgeores.2013.02.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  DonnadieuF.MerleO.BessonJ. C. (2001). Volcanic edifice stability during cryptodome intrusion. Bull. Volcanol. 63, 61–72. 10.1007/s004450000122

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  EbmeierS. K.BiggsJ.MatherT. A.AmelungF. (2013). Applicability of InSAR to tropical volcanoes: insights from Central America. Geol. Soc. Lond. Spec. Publ. 380, 15–37. 10.1144/SP380.2

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  EggersA. (1983). Temporal Gravity and Elevation Changes at Pacaya Volcano, Guatemala. Ph.D. dissertation, Dartmouth College.

  • Google Scholar
  • View reference in article

• 20

  ElsworthD.VoightB. (1996). Evaluation of volcano flank instability triggered by dyke intrusion. Geol. Soc. Lond. Spec. Publ. 110, 45–53. 10.1144/GSL.SP.1996.110.01.03

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  FaulknerD. R.MitchellT. M.HealyD.HeapM. J. (2006). Slip on ‘weak’ faults by the rotation of regional stress in the fracture damage zone. Nature444, 922–925. 10.1038/nature05353

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  FredrichJ. T.WongT. F. (1986). Micromechanics of thermally induced cracking in three crustal rocks. J. Geophys. Res. 91, 12743–12764. 10.1029/JB091iB12p12743

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  GriggsD. T.TurnerF. J.HeardH. C. (1960). Chapter 4: deformation of Rocks at 500° to 800°C. Geol. Soc. Am. Mem. 79, 39–104. 10.1130/mem79-p39

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  GurevichB. (2003). Elastic properties of saturated porous rocks with aligned fractures. J. Appl. Geophys. 54, 203–218. 10.1016/j.jappgeo.2002.11.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  HandinJ.CarterN. (1979). Rheological properties of rock at high temperatures. Proc. Congr. Int. Soc. Rock. Mech. 4, 97–106.

  • Pubmed Abstract
  • Google Scholar
  • View reference in article

• 26

  HeapM. J.FaulknerD. R. (2008). Quantifying the evolution of static elastic properties as crystalline rock approaches failure. Int. J. Rock Mech. Min. 45, 564–573. 10.1016/j.ijrmms.2007.07.018

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HeapM. J.FaulknerD. R.MeredithP. G.VinciguerraS. (2010). Elastic moduli evolution and accompanying stress changes with increasing crack damage: implications for stress changes around fault zones and volcanoes during deformation. Geophys. J. Int. 183, 225–236. 10.1111/j.1365-246X.2010.04726.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  HeapM. J.LavalléeY.LaumannA.HessK.-U.MeredithP. G.DingwellD. B. (2012). How tough is tuff in the event of fire. Geology40, 311–314. 10.1130/G32940.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  HeapM. J.LavalléeY.MeredithP. G.DingwellD. B.HuismannS.WeiseF. (2013b). The influence of thermal-stressing on the physical, mechanical and chemical properties of siliceous-aggregate high-strength concrete. Constr. Build. Mater. 42, 248–265. 10.1016/j.conbuildmat.2013.01.020

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HeapM. J.LavalléeY.PetrakovaL.BaudP.ReuschléT.VarleyN. R.et al. (2014a). Microstructural controls on the physical and mechanical properties of edifice-forming andesites at Volcán de Colima, Mexico. J. Geophys. Res. Solid. Earth119, 2925–2963. 10.1002/2013JB010521

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  HeapM. J.MolloS.VinciguerraS.LavalléeY.HessK.-U.DingwellD. B.et al. (2013a). Thermal weakening of the carbonate basement under Mt. Etna volcano (Italy): implications for volcano instability. J. Volcanol. Geotherm. Res. 250, 42–60. 10.1016/j.jvolgeores.2012.10.004

  • CrossRef
  • Google Scholar

• 32

  HeapM. J.VinciguerraS.MeredithP. G. (2009). The evolution of elastic moduli with increasing crack damage during cyclic stressing of a basalt from Mt. Etna volcano. Tectonophysics471, 153–160. 10.1016/j.tecto.2008.10.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  HeapM.XuT.ChenC.-F. (2014b). The influence of porosity and vesicle size on the brittle strength of volcanic rocks and magma. Bull. Volcanol. 76, 1–15. 10.1007/s00445-014-0856-0

  • CrossRef
  • Google Scholar

• 34

  HoekE. (1994). Strength of rock and rock masses. ISRM News J. 2, 4–16.

  • Google Scholar

• 35

  HürlimannM.MartJ.LedesmaA. (2000). Mechanical relationship between catastrophic volcanic landslides and caldera collapse. Geophys. Res. Lett. 27, 2393–2396. 10.1029/2000GL011564

  • CrossRef
  • Google Scholar

• 36

  IversonR. (1995). Can magma-injection and groundwater forces cause massive landslides on Hawaiian volcanoes?J. Volcanol. Geotherm. Res. 66, 295–308.

  • Google Scholar

• 37

  JaegerJ. C.CookN. G. W. (1976). Fundamentals of Rock Mechanics. London: Chapman and Hall.

  • Google Scholar

• 38

  JonesC.KeaneyG.MeredithP. G.MurrellS. A. F. (1997). Acoustic emission and fluid permeability measurements on thermally cracked rocks. Phys. Chem. Earth22, 13–17.

  • Google Scholar

• 39

  KendrickJ. E.LavalléeY.HessK. U.HeapM. J.GauntH. E.MeredithP. G.et al. (2013b). Tracking the permeable porous network during strain-dependent magmatic flow. J. Volcanol. Geotherm. Res. 260, 117–126. 10.1016/j.jvolgeores.2013.05.012

  • CrossRef
  • Google Scholar

• 40

  KendrickJ. E.SmithR.SammondsP.MeredithP. G.DaintyM.PallisterJ. D. (2013a). The influence of thermal and cyclic stressing on the strength of rock from Mount St. Helens, Washington. Bull. Volcanol. 75, 1–12. 10.1007/s00445-013-0728-z

  • CrossRef
  • Google Scholar

• 41

  KitamuraS.Matías GomezO. (1995). Tephra Stratigraphic Approach to the Eruptive History of Pacaya Volcano, Guatemala. Science Reports–Tohoku University, Seventh Series: Geography, 45, 1–41.

  • Google Scholar

• 42

  KonoY.MiyakeA.IshikawaM.ArimaM. (2008). Temperature derivatives of elastic wave velocities in plagioclase (An51±1) above and below the order-disorder transition temperature. Am Mineral. 93, 558–564. 10.2138/am.2008.2591

  • CrossRef
  • Google Scholar

• 43

  KueppersU.ScheuB.SpielerO.DingwellD. B. (2006). Fragmentation efficiency of explosive volcanic eruptions: a study of experimentally generated pyroclasts. J. Volcanol. Geotherm. Res. 153, 125–135. 10.1016/j.jvolgeores.2005.08.006

  • CrossRef
  • Google Scholar

• 44

  LavalléeY.HessK. U.CordonnierB.DingwellD. B. (2007). Non-Newtonian rheological law for highly crystalline dome lavas. Geology35, 843–846. 10.1130/G23594A.1

  • CrossRef
  • Google Scholar

• 45

  LavalléeY.MeredithP. G.DingwellD. B.HessK. U.WassermannJ.CordonnierB.et al. (2008). Seismogenic lavas and explosive eruption forecasting. Nature453, 507–510. 10.1038/nature06980

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  LindholmU. S.YeakleyL. M.NagyA. (1974). The dynamic strength and fracture properties of dresser basalt. Int. J. Rock Mech. Min. 11, 181–191. 10.1016/0148-9062(74)90885-7

  • CrossRef
  • Google Scholar

• 47

  LoaizaS.FortinJ.SchubnelA.GueguenY.VinciguerraS.MoreiraM. (2012). Mechanical behavior and localized failure modes in a porous basalt from the Azores. Geophys. Res. Lett. 39, 6. 10.1029/2012gl053218

  • CrossRef
  • Google Scholar

• 48

  Matías GómezR. O. (2009). Volcanological Map of the 1961–2009 Eruption of Volcán de Pacaya, Guatemala. M.S. thesis, Michigan Technological University, Guatemala.

  • Google Scholar

• 49

  Matías GómezR. O.RoseW. I.PalmaJ. L.Escobar-WolfR. (2012). Notes on a map of the 1961–2010 eruptions of Volcán de Pacaya, Guatemala. Geol. Soc. Am. Digit. Map Chart Ser. 10, 10. 10.1130/2012.DMCH010

  • CrossRef
  • Google Scholar

• 50

  MeredithP. G.AtkinsonB. K. (1982). High-temperature tensile crack propagation in quartz: experimental results and application to time-dependent earthquake rupture, Earthquake. Predic. Res. 1, 377–391.

  • Google Scholar

• 51

  MerriamR.RiekeH. H.KimY. C. (1970). Tensile strength related to mineralogy and texture of some granitic rocks. Eng. Geo. 4, 155–160. 10.1016/0013-7952(70)90010-4

  • CrossRef
  • Google Scholar

• 52

  MuellerS.MelnikO.SpielerO.ScheuB.DingwellD. B. (2005). Permeability and degassing of dome lavas undergoing rapid decompression: an experimental determination. Bull. Volcanol. 67, 526–538. 10.1007/s00445-004-0392-4

  • CrossRef
  • Google Scholar

• 53

  OkuboC. H. (2004). Rock mass strength and slope stability of the Hilina slump, Kilauea volcano, Hawaii. J. Volcanol. Geotherm. Res. 138, 43–76. 10.1016/j.jvolgeores.2004.06.006

  • CrossRef
  • Google Scholar

• 54

  PatersonM. S.WongT. F. (2005). Experimental Rock Deformation-the Brittle Field. Berlin; Heidelberg: Springer Science and Business Media.

  • Google Scholar

• 55

  PřikrylR. (2001). Some microstructural aspects of strength variation in rocks. Int. J. Rock Mech. Min. 38, 671–682. 10.1016/S1365-1609(01)00031-4

  • CrossRef
  • Google Scholar

• 56

  ReidM. E.ChristianS. B.BrienD. L. (2000). Gravitational stability of three-dimensional stratovolcano edifices. J. Geophys. Res. 105, 6043–6056. 10.1029/1999JB900310

  • CrossRef
  • Google Scholar

• 57

  RocchiV.SammondsP. R.KilburnC. R. (2003). Flow and fracture maps for basaltic rock deformation at high temperatures. J. Volcanol. Geotherm. Res. 120, 25–42. 10.1016/S0377-0273(02)00343-8

  • CrossRef
  • Google Scholar

• 58

  RocchiV.SammondsP. R.KilburnC. R. J. (2004). Fracturing of Etnean and Vesuvian rocks at high temperatures and low pressures. J. Volcanol. Geotherm. Res. 132, 137–157. 10.1016/S0377-0273(03)00342-1

  • CrossRef
  • Google Scholar

• 59

  RoseW. I.PalmaJ. L.Escobar WolfR.Matías GomezR. O. (2013). A 50 yr eruption of a basaltic composite cone: Pacaya, Guatemala. Geol. S. Am. S. 498, 1–21. 10.1130/2013.2498(01)

  • CrossRef
  • Google Scholar

• 60

  SammisC. G.AshbyM. F. (1986). The failure of brittle porous solids under compressive stress states. Acta. Metall. 34, 511–526. 10.1016/0001-6160(86)90087-8

  • CrossRef
  • Google Scholar

• 61

  SchaeferL. N.OommenT.CorazzatoC.TibaldiA.Escobar-WolfR.RoseW. I. (2013). An integrated field-numerical approach to assess slope stability hazards at volcanoes: the examples of Pacaya Volcano, Guatemala. Bull. Volcanol. 75, 1–18. 10.1007/s00445-013-0720-7

  • CrossRef
  • Google Scholar

• 62

  ScheuB.KueppersU.MuellerS.SpielerO.DingwellD. B. (2008). Experimental volcanology on eruptive products of Unzen volcano. J. Volcanol. Geotherm. Res. 175, 110–119. 10.1016/j.jvolgeores.2008.03.023

  • CrossRef
  • Google Scholar

• 63

  ScholzC. H. (1968). Microfracturing and the inelastic deformation of rock in compre ssion. J. Geophys. Res. 73, 1417–1432. 10.1029/JB073i004p01417

  • CrossRef
  • Google Scholar

• 64

  SimmonsG.CooperH. W. (1978). Thermal cycling cracks in three igneous rocks. Int. J. Rock Mech. Min. 15, 145–148. 10.1016/0148-9062(78)91220-2

  • CrossRef
  • Google Scholar

• 65

  SmithR.SammondsP. R.KilburnC. R. (2005). Experimental studies of lava dome fracture, in AGU Fall Meeting Abstracts (San Francisco, CA), 1.

  • Google Scholar

• 66

  SongH.ZhangH.KangY.HuangG.FuD.QuC. (2013). Damage evolution study of sandstone by cyclic uniaxial test and digital image correlation. Tectonophysics608, 1343–1348. 10.1016/j.tecto.2013.06.007

  • CrossRef
  • Google Scholar

• 67

  SousaJ.VoightB. (1995). Multiple-pulsed debris avalanche emplacement at Mount St. Helens in 1980: evidence from numerical continuum flow simulations. J. Volcanol. Geotherm. Res. 66, 227–250. 10.1016/0377-0273(94)00067-Q

  • CrossRef
  • Google Scholar

• 68

  SpielerO.KennedyB.KueppersU.DingwellD. B.ScheuB.TaddeucciJ. (2004). The fragmentation threshold of pyroclastic rocks. Earth Planet Sci. Lett. 226, 139–148. 10.1016/j.epsl.2004.07.016

  • CrossRef
  • Google Scholar

• 69

  TapponnierP.BraceW. F. (1976). Development of stress-induced microcracks in Westerly granite. Int. J. Rock Mech. Min. 13, 103–112. 10.1016/0148-9062(76)91937-9

  • CrossRef
  • Google Scholar

• 70

  ThomasM. E.PetfordN.BromheadE. N. (2004). Volcanic rock-mass properties from Snowdonia and Tenerife: implications for volcano edifice strength. J. Geol. Soc. Lond. 161, 939–946. 10.1144/0016-764903-166

  • CrossRef
  • Google Scholar

• 71

  VallanceJ. W.SiebertL.RoseW. I.Raul GirónJ.BanksN. G. (1995). Edifice collapse and related hazards in Guatemala. J. Volcanol. Geotherm. Res. 66, 337–355. 10.1016/0377-0273(94)00076-S

  • CrossRef
  • Google Scholar

• 72

  VinciguerraS.TrovatoC.MeredithP. G.BensonP. M. (2005). Relating seismic velocities, thermal cracking and permeability in Mt. Etna and Iceland basalts. Int. J. Rock Mech. Min. Sci. 42, 900–910. 10.1016/j.ijrmms.2005.05.022

  • CrossRef
  • Google Scholar

• 73

  VoightB.JandaR. J.GlickenH.DouglassP. M. (1983). Nature and mechanics of the Mount St Helens rockslide-avalanche of 18 May 1980. Géotechnique33, 243–273. 10.1680/geot.1983.33.3.243

  • CrossRef
  • Google Scholar

• 74

  WangX. Q.SchubnelA.FortinJ.GuéguenY.GeH. K. (2013). Physical properties and brittle strength of thermally cracked granite under confinement. J. Geophys. Res. Solid Earth118, 6099–6112. 10.1002/2013JB010340

  • CrossRef
  • Google Scholar

• 75

  WundermanR. L.RoseW. I. (1984). Amatitlán, an actively resurging cauldron 10 km south of Guatemala City. J. Geophys. Res. Solid Earth89, 8525–8539. 10.1029/JB089iB10p08525

  • CrossRef
  • Google Scholar

• 76

  ZhangL.MaoX.LuA. (2009). Experimental study on the mechanical properties of rocks at high temperature. Sci. China Ser. E52, 641–646. 10.1007/s11431-009-0063-y

  • CrossRef
  • Google Scholar

• 77

  ZimbelmanD. R.WattersR. J.FirthI. R.BreitG. N.Carrasco-NunezG. (2004). Stratovolcano stability assessment methods and results from Citlaltépetl, Mexico. Bull. Volcanol. 66, 66–79. 10.1007/s00445-003-0296-8

  • CrossRef
  • Google Scholar