Nick Varley

Microstructural controls on the physical and mechanical properties of edifice-forming andesites at Volcán de Colima, Mexico

The reliable assessment of volcanic unrest must rest on an understanding of the rocks that form the edifice. It is their microstructure that dictates their physical properties and mechanical behavior and thus the response of the edifice to stress perturbations during unrest. We evaluate the interplay between microstructure and rock properties for a suite of edifice-forming rocks from Volcan de Colima (Mexico). Microstructural analyses expose (1) a pervasive, isotropic microcrack network, (2) a high, subspherical vesicle density, and (3) a wide vesicle size distribution. This complex microstructure severely impacts their physical and mechanical properties. In detail, porosities are high and range from 8 to 29%. As a consequence, elastic wave velocities, Youngs moduli, and uniaxial compressive strengths are low, and permeabilities are high. All of the rock properties demonstrate a wide range. For example, strength decreases by a factor of 8 and permeability increases by 4 orders of magnitude over the porosity range. Below a porosity of 11–14%, the permeability-porosity trend follows a power law with a much higher exponent. Microstructurally, this represents a critical vesicle content that efficiently connects the microcrack population and permits a much more direct path through the sample, rather than restricting flow to long and tortuous microcracks. Values of tortuosity inferred from the Kozeny-Carman permeability model support this hypothesis. However, we find that the complex microstructure precludes a complete description of their mechanical behavior through micromechanical modeling. We urge that the findings of this study be considered in volcanic hazard assessments at andesitic stratovolcanoes.

Long-period earthquakes and co-eruptive dome inflation seen with particle image velocimetry.

Dome growth and explosive degassing are fundamental processes in the cycle of continental arc volcanism. Because both processes generate seismic energy, geophysical field studies of volcanic processes are often grounded in the interpretation of volcanic earthquakes. Although previous seismic studies have provided important constraints on volcano dynamics, such inversion results do not uniquely constrain magma source dimension and material properties. Here we report combined optical geodetic and seismic observations that robustly constrain the sources of long-period volcanic earthquakes coincident with frequent explosive eruptions at the volcano Santiaguito, in Guatemala. The acceleration of dome deformation, extracted from high-resolution optical image processing, is shown to be associated with recorded long-period seismic sources and the frequency content of seismic signals measured across a broadband network. These earthquake sources are observed as abrupt subvertical surface displacements of the dome, in which 20–50-cm uplift originates at the central vent and propagates at ∼50 m s-1 towards the 200-m-diameter periphery. Episodic shifts of the 20–80-m thick dome induce peak forces greater than 109 N and reflect surface manifestations of the volcanic long-period earthquakes, a broad class of volcano seismic activity that is poorly understood and observed at many volcanic centres worldwide. On the basis of these observations, the abrupt mass shift of solidified domes, conduit magma or magma pads may play a part in generating long-period earthquakes at silicic volcanic systems.

Identification of structural controls in an active lava dome with high resolution DEMs: Volcán de Colima, Mexico

Monitoring the topography of active lava domes is critical for detecting changes that may trigger or influence collapse or explosive activity. Internal dome structure and conditions are more difficult to elucidate, but also play vital roles. Here, we describe the exposure (following an explosion) of significant scarps in the active dome at Volcan de Colima, Mexico, that are interpreted as evidence of brittle failure planes and a complex internal dome morphology. In the first use of automated 3D computer vision reconstruction techniques (structure-from-motion and multi-view stereo, SfM-MVS) on an active volcanic dome, we derive high resolution surface models from oblique and archive photographs taken with a consumer camera. The resulting 3D models were geo-referenced using features identified in a web-sourced orthoimage; no ground-based measurements were required. In December 2010, the dome (2.14×106 m3) had a flat upper surface, reflecting an overall ductile emplacement regime. Between then and May 2011, a period of low explosivity was accompanied by a small volume loss (0.4×105 m3) and arcuate steps appeared in the dome surface, suggesting the presence of localized planes of weakness. The complex array of summit scarps was exposed following a significant explosion in June 2011, and is interpreted to be the surface expression of fault planes in the dome. The 1-m resolution DEMs indicated that the region of greatest volume loss was not coincident with the assumed location of the conduit, and that heterogeneity within the dome may have been important during the June explosion.

Particle sizes of andesitic ash fallout from vertical eruptions and co-pyroclastic flow clouds, Volcán de Colima, Mexico

We report particle-size distributions for andesitic ash fall from vertical eruptions and clouds above block-and-ash fl ows at Volcan de Colima, Mexico, during 2004‐2006. We analyzed 17 samples using laser diffraction: 11 from vertical eruptions and 6 from block-and-ash fl ows (copyroclastic fl ow ash). Vertical eruptions produce well-sorted fall deposits, whereas co-pyroclastic fall deposits are poorly sorted, with high proportions of very fi ne grained ash (<30 µm). Statistical analysis shows particle-size distributions of vertical eruptions are more leptokurtic (peaked) than co-pyroclastic samples. Deconvolution of grain-size histograms shows that copyroclastic samples have at least one subpopulation with a mode of 8.3‐8.7 φ (2.4‐3.1 µm). Estimates of the number of particles in different size ranges show that co-pyroclastic samples contain much greater numbers of very fi ne particles than vertical eruption samples. Our results provide no direct evidence that milling or comminution produces hazardous fi ne ash particles in pyroclastic fl ows, but are consistent with that interpretation.

Airborne thermal remote sensing of the Volcán de Colima (Mexico) lava dome from 2007 to 2010

Abstract We investigate high-resolution digital photographs and infrared images of the lava dome eruption at Volcán de Colima, from 2007 to 2010. Qualitative observations provide insight into active volcanic processes (e.g. rockfalls and fracturing) and show that, as the dome advances a substantial cooled talus apron develops, which stabilizes the structure. Progressive collapse of the talus apron as it reaches the crater rim corresponds with the development of a lava lobe, extruding hot lava from deeper within the dome. Quantitative dome surface temperature time-series show that the highest temperature hotspots migrate from the dome sides (250–380 °C) to the top (150–300 °C) and finally to the lava lobe (220–400 °C) as the structurally unstable areas expose fresh material. Net surface heat loss from the dome ranges from 5 to 30 MW, comparable to other dome forming systems. Heat budget calculations confirm that the lava dome retained a hot viscous core throughout the period 2007–2010. We propose that the mechanical stability of the Volcán de Colima dome arises from the shear strength of flanking talus which stabilizes the hot viscous core.

Blowing Off Steam: Tuffisite Formation As a Regulator for Lava Dome Eruptions

Tuffisites are veins of variably sintered, pyroclastic particles that form in conduits and lava domes as a result of localized fragmentation events during gas-and-ash explosions. Those observed in-situ on the active 2012 lava dome of Volcan de Colima range from voids with intra-clasts showing little movement and interpreted to be failure-nuclei, to sub-parallel lenses of sintered granular aggregate interpreted as fragmentation horizons, through to infilled fractures with evidence of viscous remobilization. All tuffisites show evidence of sintering. Further examination of the complex fracture-and-channel patterns reveals viscous backfill by surrounding magma, suggesting that lava fragmentation was followed by stress relaxation and continued viscous deformation as the tuffisites formed. The natural tuffisites are more permeable than the host andesite, and have a wide range of porosity and permeability compared to a narrower window for the host rock, and gauging from their significant distribution across the dome, we posit that the tuffisite veins may act as important outgassing pathways. To investigate tuffisite formation we crushed and sieved andesite from the lava dome and sintered it at magmatic temperatures for different times. We then assessed the healing and sealing ability by measuring porosity and permeability, showing that sintering reduces both over time. During sintering the porosity-permeability reduction occurs due to the formation of viscous necks between adjacent grains, a process described by the neck-formation model of Frenkel (1945). This process leads the granular starting material to a porosity-permeability regime anticipated for effusive lavas, and which describes the natural host lava as well as the most impervious of natural tuffisites. This suggests that tuffisite formation at Volcan de Colima constructed a permeable network that enabled gas to bleed passively from the magma. We postulate that this progressively reduced the lava dome’s ability to seal and build pressure that drives explosions. Indeed, the time interval between explosions during 2007 to 2011 gradually increased before the onset of a period of quiescence starting in June 2011. We suggest that the permeability evolution during tuffisite formation has important consequences for modeling of gas-and-ash explosions, common at dome-forming volcanoes.

Geophysical characterization of hydrothermal systems and intrusive bodies, El Chichon volcano (Mexico)

[1] The 1982 explosive eruptions of El Chichon volcano (Chiapas, Mexico) destroyed the inner dome and created a 1-km-wide and 180-m-deep crater within the somma crater. A shallow hydrothermal system was exposed to the surface of the new crater floor and is characterized by an acid crater lake, a geyser-like Cl-rich spring (soap pool), and numerous fumarole fields. Multiple geophysical surveys were performed to define the internal structure of the volcanic edifice and its hydrothermal system. We carried out a high-resolution ground-based geomagnetic survey in the 1982 crater and its surroundings and 38 very low frequency (VLF) transects around the crater lake. A 3-D inversion of the ground-based magnetic data set highlighted three high-susceptibility isosurfaces, interpreted as highly magnetized bodies beneath the 1982 crater floor. Inversion of a digitized regional aeromagnetic map highlighted four major deeply rooted cryptodomes, corresponding to major topographic highs and massive lava dome outcrops outside and on the somma rim. The intracrater magnetic bodies correspond closely to the active hydrothermal vents and their modeled maximum basal depth matches the elevation of the springs on the flanks of the volcano. Position, dip, and vertical extent of active and extinct hydrothermal vents identified by VLF-EM surveys match the magnetic data set. We interpret the shallow lake spring hydrothermal system to be mostly associated with buried remnants of the 550 BP dome, but the Cl-rich soap pool may be connected to a small intrusion emplaced at shallow depth during the 1982 eruption.

Pulsed Vulcanian explosions : A characterization of eruption dynamics using Doppler radar

Understanding the dynamics of ongoing volcanic eruptions is essential for predicting the input and transport of volcanic ash in the atmosphere. To constrain near-vent dynamic processes of explosive Vulcanian events, we used Doppler radar measurements, providing tephra velocities and a proxy of the mass flux, in two field experiments at Volcan de Colima (Mexico) and Santiaguito (Guatemala). We find that explosive eruptions at both volcanoes consist of individual degassing pulses. The analysis of the timing of such pulses shows that both volcanoes have preferred interpulse times of 3 s (Santiaguito) and 2–5 s (Volcan de Colima). The interpulse time during one event may change, but it often returns to the preferred interpulse time. This behavior is similar at both volcanoes and the interpulse time roughly follows a log-logistic distribution indicating the interplay of two competing processes. These could be short-duration degassing of the uppermost conduit versus decreasing permeability due to progressive gas loss and compaction during an eruption.

Seismic and experimental insights into eruption precursors at Volcán de Colima

Abstract We combine geophysical and experimental observations to interpret preeruptive unrest at Volcán de Colima in 1998. 17,893 volcanic earthquakes were detected between 1 October and 31 December 1998, including 504 clusters. Using seismic ambient noise interferometry, we observe a drop in velocity prior to the eruption linked to damage accumulation during magma ascent. This is supported by experimental observations where static stress causes a velocity decrease prior to failure. Furthermore, we observe acoustic emission clusters during the experiments, with lower porosity samples producing higher numbers of repeaters. This behavior introduces tensile failure as an additional viable mechanism for clusters during magma ascent. The findings suggest that preeruptive magma ascent may be monitored to variable degrees of accuracy via descriptions of damage accumulation and associated seismic velocity changes.

Tremor and its duration‐amplitude distribution at Popocatépetl volcano, Mexico

The duration-amplitude distribution was calculated for the tremor observed at Popocatepetl volcano during episodes of activity in 2000, 2012, 2013 and 2014. An exponential function was used to obtain a good fit for the duration-amplitude distribution, and the source of volcanic tremor is probably generated by the transportation of magmatic fluids and its coupling with the host rock within the volcanic conduit. In particular, harmonic tremor has shown large amplitudes, durations, and mean values of amplitude, more than spasmodic or pulsalting tremor. This is due to different generation mechanisms: in the case of harmonic tremor, it is produced during magma ascent and lava dome growth, whilst spasmodic and pulsalting tremor are associated with fragmentation of the lava dome and gas emissions. This paper presents the duration-amplitude distribution as a method to estimate the intensity of the tremor at Popocatepetl, a volcano with the major risk in all Mexico.