Magma intrusion beneath long valley caldera confirmed by temporal changes in gravity

A New Analysis of Caldera Unrest through the Integration of Geophysical Data and FEM Modeling: The Long Valley Caldera Case Study

The Long Valley Caldera, located at the eastern edge of the Sierra Nevada range in California, has been in a state of unrest since the late 1970s. Seismic, gravity and geodetic data strongly suggest that the source of unrest is an intrusion beneath the caldera resurgent dome. However, it is not clear yet if the main contribution to the deformation comes from pulses of ascending high-pressure hydrothermal fluids or low viscosity magmatic melts. To characterize the nature of the intrusion, we developed a 3D finite element model which includes topography and crust heterogeneities. We first performed joint numerical inversions of uplift and Electronic Distance Measurement baseline length change data, collected during the period 1985–1999, to infer the deformation-source size, position, and overpressure. Successively, we used this information to refine the source overpressure estimation, compute the gravity potential and infer the intrusion density from the inversion of deformation and gravity data collected in 1982–1998. The deformation source is located beneath the resurgent dome, at a depth of 7.5 ± 0.5 km and a volume change of 0.21 ± 0.04 km3. We assumed a rhyolite compressibility of 0.026 ± 0.0011 GPa−1 (volume fraction of water between 0% and 30%) and estimated a reservoir compressibility of 0.147 ± 0.037 GPa−1. We obtained a density of 1856 ± 72 kg/m3. This density is consistent with a rhyolite melt, with 20% to 30% of dissolved hydrothermal fluids.

Basal nucleation and the prevalence of ascending swarms in Long Valley caldera

Earthquake swarms are ubiquitous in volcanic systems, being manifestations of underlying nontectonic processes such as magma intrusions or volatile fluid transport. The Long Valley caldera, California, is one such setting where episodic earthquake swarms and persistent uplift suggest the presence of active magmatism. We quantify the long-term spatial and temporal characteristics of seismicity in the region using cluster analysis on a 25-year high-resolution earthquake catalog derived using leading-edge deep-learning algorithms. Our results show that earthquake swarms beneath the caldera exhibit enlarged families with statistically significant tendency for upward migration patterns. The ascending swarms tend to nucleate at the base of the seismogenic zone with a spatial footprint that is laterally constrained by the southern rim of the caldera. We suggest that these swarms are driven by the transport of volatile-rich fluids released from deep volcanic processes. The observations highlight the potential for extreme spatial segmentation of earthquake triggering processes in magmatic systems.

Thoughts on the criteria to determine the origin of volcanic unrest as magmatic or non-magmatic

As our ability to detect volcanic unrest improves, we are increasingly confronted with the question of whether the unrest has a magmatic origin (magma on the move) or a non-magmatic origin from a change in the hydrothermal system (fluids that are not magma on the move) or tectonic processes. The cause of unrest has critical implications for the potential eruptive hazard (e.g. used in constructing Bayesian Event Trees), but is frequently the subject of debate, even at well-studied systems. Here, we propose a set of multi-disciplinary observations and numerical models that could be used to evaluate conceptual models about the cause of unrest. These include measurements of gas fluxes and compositions and the isotopic signature of some components (e.g. H₂, He, C, SO₂, H₂ O, CH₄ and CO₂), the spatial and temporal characteristics of ground deformation, thermal output, seismicity, changes in gravity, and whether there is topographic uplift or subsidence spanning hundreds to thousands of years. In several volcanic systems, both magmatic and non-magmatic unrest is occurring at the same time. While none of these observations or models is diagnostic on its own, we illustrate several examples where they have been used together to make a plausible conceptual model of one or more episodes of unrest and whether eruptions did or did not follow the unrest. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.

Tree-ring width reveals the preparation of the 1974 Mt. Etna eruption

Reduced near-infrared reflectance observed in September 1973 in Skylab images of the western flank of Mt. Etna has been interpreted as an eruption precursor of the January 1974 eruption. Until now, it has been unclear when this signal started, whether it was sustained and which process(es) could have caused it. By analyzing tree-ring width time-series, we show that the reduced near-infrared precursory signal cannot be linked to a reduction in annual tree growth in the area. However, comparing the tree-ring width time-series with both remote sensing observations and volcano-seismic activity enables us to discuss the starting date of the pre-eruptive period of the 1974 eruption.

Insensitivity of Tree-Ring Growth to Temperature and Precipitation Sharpens the Puzzle of Enhanced Pre-Eruption NDVI on Mt. Etna (Italy)

On Mt. Etna (Italy), an enhanced Normalized Difference in Vegetation Index (NDVI) signature was detected in the summers of 2001 and 2002 along a distinct line where, in November 2002, a flank eruption subsequently occurred. These observations suggest that pre-eruptive volcanic activity may have enhanced photosynthesis along the future eruptive fissure. If a direct relation between NDVI and future volcanic eruptions could be established, it would provide a straightforward and low-cost method for early detection of upcoming eruptions. However, it is unclear if, or to what extent, the observed enhancement of NDVI can be attributed to volcanic activity prior to the subsequent eruption. We consequently aimed at determining whether an increase in ambient temperature or additional water availability owing to the rise of magma and degassing of water vapour prior to the eruption could have increased photosynthesis of Mt. Etna's trees. Using dendro-climatic analyses we quantified the sensitivity of tree ring widths to temperature and precipitation at high elevation stands on Mt. Etna. Our findings suggest that tree growth at high elevation on Mt. Etna is weakly influenced by climate, and that neither an increase in water availability nor an increase in temperature induced by pre-eruptive activity is a plausible mechanism for enhanced photosynthesis before the 2002/2003 flank eruption. Our findings thus imply that other, yet unknown, factors must be sought as causes of the pre-eruption enhancement of NDVI on Mt. Etna.

ROCK PHYSICS. Rock physics of fibrous rocks akin to Roman concrete explains uplifts at Campi Flegrei Caldera

Uplifts in the Campi Flegrei caldera reach values unsurpassed anywhere in the world (~2 meters). Despite the marked deformation, the release of strain appears delayed. The rock physics analysis of well cores highlights the presence of two horizons, above and below the seismogenic area, underlying a coupled process. The basement is a calc-silicate rock housing hydrothermal decarbonation reactions, which provide lime-rich fluids. The caprock above the seismogenic area has a pozzolanic composition and a fibril-rich matrix that results from lime-pozzolanic reactions. These findings provide evidence for a natural process reflecting that characterizing the cementitious pastes in modern and Roman concrete. The formation of fibrous minerals by intertwining filaments confers shear and tensile strength to the caprock, contributing to its ductility and increased resistance to fracture.

Accelerated uplift and magmatic intrusion of the Yellowstone caldera, 2004 to 2006

The Yellowstone caldera began a rapid episode of ground uplift in mid-2004, revealed by Global Positioning System and interferometric synthetic aperture radar measurements, at rates up to 7 centimeters per year, which is over three times faster than previously observed inflation rates. Source modeling of the deformation data suggests an expanding volcanic sill of approximately 1200 square kilometers at a 10-kilometer depth beneath the caldera, coincident with the top of a seismically imaged crustal magma chamber. The modeled rate of source volume increase is 0.1 cubic kilometer per year, similar to the amount of magma intrusion required to supply the observed high heat flow of the caldera. This evidence suggests magma recharge as the main mechanism for the accelerated uplift, although pressurization of magmatic fluids cannot be ruled out.