Battaglia, Maurizio

Understanding whether fractures and faults impact the CO2 migration through the overburden is critical in the evaluation and monitoring of CO2 geological storage sites. We present a numerical model and workflow to describe the hydraulic behaviour of a fault located in the shallow part of the overburden. This helps to evaluate the sealing potential of the system in case of unwanted CO2 migration toward the surface and to design an efficient monitoring plan. The model configuration is representative of several experiments performed at real sites under quite shallow conditions (50–500 m). The model results, applied to a selected fault outcropping in the Apennines (Italy), show that most of the gas migrates through the high permeable footwall damage zone. A significant amount of gas then dissolves into the water, emphasizing the importance of accurate modelling to assess the hazard of CO2 leakage into near-surface aquifers or to the surface. The occurrence of pressure buildup close to the fault core points out that detailed modelling of the migration conditions is required to predict gas path through a fault zone.

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.

The ability of Synthetic Aperture Radar (SAR) to image the Earth's surface, even through dense cloud cover and in night-and-day conditions, can facilitate the evaluation and monitoring of natural hazards and the management of natural disasters. The family of SAR satellite sensors orbits the Earth at an altitude ranging from 500 to 800 km, following sun-synchronous, near-polar orbits, slightly inclined with respect to Earth meridians. The most commonly used bands in SAR applications are the C-band (5–6 GHz, ~5, 6 cm wavelength), the X-band (8–12 GHz, ~3, 1 cm wavelength), and the L-band (1–2 GHz ~23 cm wavelength) with a temporal resolution depending on the satellite revisiting time. The availability of SAR has made a new spectrum of measurements possible on a global and spatial scale not attainable by ground-based studies, revealing critical insights into remote or poorly understood areas (e.g., Biggs et al., 2014). This Research Topics presents a review of articles on the state-of-art in the application of SAR sensors to study surface deformation in different geologic environments and triggered by a variety of processes. The topics discussed range from the analysis of co-seismic deformation (Merryman Boncori) to studies of volcanic unrest (Dzurisin et al.; Garthwaite et al.), monitoring of landslides (Bianchini et al.) and ground subsidence in urban areas (Solari et al.).

Calderas are collapse structures related to the emptying of magmatic reservoirs, often associated with large eruptions from long-lived magmatic systems. Understanding how magma is transferred from a magma reservoir to the surface before eruptions is a major challenge. Here we exploit the historical, archaeological and geological record of Campi Flegrei caldera to estimate the surface deformation preceding the Monte Nuovo eruption and investigate the shallow magma transfer. Our data suggest a progressive magma accumulation from ~1251 to 1536 in a 4.6 ± 0.9 km deep source below the caldera centre, and its transfer, between 1536 and 1538, to a 3.8 ± 0.6 km deep magmatic source ~4 km NW of the caldera centre, below Monte Nuovo; this peripheral source fed the eruption through a shallower source, 0.4 ± 0.3 km deep. This is the first reconstruction of pre-eruptive magma transfer at Campi Flegrei and corroborates the existence of a stationary oblate source, below the caldera centre, that has been feeding lateral eruptions for the last ~5 ka. Our results suggest: 1) repeated emplacement of magma through intrusions below the caldera centre; 2) occasional lateral transfer of magma feeding non-central eruptions within the caldera. Comparison with historical unrest at calderas worldwide suggests that this behavior is common.

Surface tension plays an important role in the nucleation of H2O gas bubbles in magmatic melts and in the time-dependent rheology of bubble-bearing magmas. Despite several experimental studies, a physics based model of the surface tension of magmatic melts in contact with H2O is lacking. This paper employs gradient theory to develop a thermodynamical model of equilibrium surface tension of silicate melts in contact with H2O gas at low to moderate pressures. In the last decades, this approach has been successfully applied in studies of industrial mixtures but never to magmatic systems. We calibrate and verify the model against literature experimental data, obtained by the pendant drop method, and by inverting bubble nucleation experiments using the Classical Nucleation Theory (CNT). Our model reproduces the systematic decrease in surface tension with increased H2O pressure observed in the experiments. On the other hand, the effect of temperature is confirmed by the experiments only at high pressure. At atmospheric pressure, the model shows a decrease of surface tension with temperature. This is in contrast with a number of experimental observations and could be related to microstructural effects that cannot be reproduced by our model. Finally, our analysis indicates that the surface tension measured inverting the CNT may be lower than the value measured by the pendant drop method, most likely because of changes in surface tension controlled by the supersaturation.

Analysis of microgravity and surface displacement data collected at the summit of Kīlauea Volcano, Hawaii (USA), between December 2009 and November 2012 suggests a net mass accumulation at ~1.5 km depth beneath the northeast margin of Halema‘uma‘u Crater, within Kīlauea Caldera. Although residual gravity increases and decreases are accompanied by periods of uplift and subsidence of the surface, respectively, the volume change inferred from the modeling of interferometric synthetic aperture radar deformation data can account for only a small portion (as low as 8%) of the mass addition responsible for the gravity increase. We propose that since the opening of a new eruptive vent at the summit of Kīlauea in 2008, magma rising to the surface of the lava lake outgasses, becomes denser, and sinks to deeper levels, replacing less dense gas-rich magma stored in the Halema‘uma‘u magma reservoir. In fact, a relatively small density increase (<200 kgm_3) of a portion of the reservoir can produce the positive residual gravity change measured during the period with the largest mass increase, between March 2011 and November 2012. Other mechanisms may also play a role in the gravity increase without producing significant uplift of the surface, including compressibility of magma, formation of olivine cumulates, and filling of void space by magma. The rate of gravity increase, higher than during previous decades, varies through time and seems to be directly correlated with the volcanic activity occurring at both the summit and the east rift zone of the volcano.

Many volcanic edifices are subject to flank failure, usually produced by a combination of events, rather than any single process. From a dynamic point of view, the cause of collapse can be divided into factors that contribute to an increase in shear stress, and factors that contribute to the reduction in the friction coefficient μ of a potential basal failure plane. We study the potential for flank failure at Mount Etna considering a schematic section of the eastern flank, approximated by a wedge-like block. For such geometry, we perform a (steady state) limit equilibrium analysis: the resolution of the forces parallel to the possible basal failure plane allows us to determine the total force acting on the potentially unstable wedge. An estimate of the relative strength of these forces suggests that, in first approximation, the stability is controlled primarily by the balance between block weight, lithostatic load and magmatic forces. Any other force (sea load, hydrostatic uplift, and the uplift due to mechanical and thermal pore-fluid pressure) may be considered of second order. To study the model sensitivity, we let the inferred slope α of the basal surface failure vary between −10° and 10°, and consider three possible scenarios: no magma loading, magmastatic load, and magmastatic load with magma overpressure. We use error propagation to include in our analysis the uncertainties in the estimates of the mechanics and geometrical parameters controlling the block equilibrium. When there is no magma loading, the ratio between destabilizing and stabilizing forces is usually smaller than the coefficient of friction of the basal failure plane. In the absence of an initiating mechanism, and with the nominal values of the coefficient of friction μ = 0.7 ± 0.1 proposed, the representative wedge will remain stable or continue to move at constant speed. In presence of magmastatic forces, the influence of the lateral restraint decreases. If we consider the magmastatic load only, the block will remain stable (or continue to move at constant speed), unless the transient mechanical and thermal pressurization significantly decrease the friction coefficient, increasing the instability of the flank wedge for α > 5° (seaward dipping decollement). When the magma overpressure contribution is included in the equilibrium analysis, the ratio between destabilizing and stabilizing forces is of the same order or larger than the coefficient of friction of the basal failure plane, and the block will become unstable (or accelerate), especially in the case of the reduction in friction coefficient. Finally, our work suggests that the major challenge in studying flank instability at Mount Etna is not the lack of an appropriate physical model, but the limited knowledge of the mechanical and geometrical parameters describing the block equilibrium.

The Long Valley caldera (California) formed ~760,000 yr ago following the massive eruption of the Bishop Tuff. Postcaldera volcanism in the Long Valley volcanic fi eld includes lava domes as young as 650 yr. The recent geological unrest is characterized by uplift of the resurgent dome in the central section of the caldera (75 cm in the past 33 yr) and earthquake activity followed by periods of relative quiescence. Since the spring of 1998, the caldera has been in a state of low activity. The cause of unrest is still debated, and hypotheses range from hybrid sources (e.g., magma with a high percentage of volatiles) to hydrothermal fl uid intrusion. Here, we present observations of surface deformation in the Long Valley region based on differential synthetic aperture radar interferometry (InSAR), leveling, global positioning system (GPS), two-color electronic distance meter (EDM), and microgravity data. Thanks to the joint application of InSAR and microgravity data, we are able to unambiguously determine that magma is the cause of unrest.

Time-dependent gravimetric measurements can detect subsurface processes long before magma flow leads to earthquakes or other eruption precursors. The ability of gravity measurements to detect subsurface mass flow is greatly enhanced if gravity measurements are analyzed and modeled with ground-deformation data. Obtaining the maximum information from microgravity studies requires careful evaluation of the layout of network benchmarks, the gravity environmental signal, and the coupling between gravity changes and crustal deformation. When changes in the system under study are fast (hours to weeks), as in hydrothermal systems and restless volcanoes, continuous gravity observations at selected sites can help to capture many details of the dynamics of the intrusive sources. Despite the instrumental effects, mainly caused by atmospheric temperature, results from monitoring at Mt. Etna volcano show that continuous measurements are a powerful tool for monitoring and studying volcanoes. Several analytical and numerical mathematical models can be used to fit gravity and deformation data. Analytical models offer a closed-form description of the volcanic source. In principle, this allows one to readily infer the relative importance of the source parameters. In active volcanic sites such as Long Valley caldera (California, U.S.A.) and Campi Flegrei (Italy), careful use of analytical models and high-quality data sets has produced good results. However, the simplifications that make analytical models tractable might result in misleading volcanological interpretations, particularly when the real crust surrounding the source is far from the homogeneous/isotropic assumption. Using numerical models allows consideration of more realistic descriptions of the sources and of the crust where they are located (e.g., vertical and lateral mechanical discontinuities, complex source geometries, and topography). Applications at Teide volcano (Tenerife) and Campi Flegrei demonstrate the importance of this more realistic description in gravity calculations.

In the last four decades, Campi Flegrei caldera has been the world’s most active caldera characterized by intense unrest episodes involving huge ground deformation and seismicity, but, at the time of writing, has not culminated in an eruption. We present a careful review, with new analyses and interpretation, of all the data and recent research results. We deal with three main problems: the tentative reconstruction of the substructure; the modelling of unrest episodes to shed light on possible pre-eruptive scenarios; and the probabilistic estimation of the hazards from explosive pyroclastic products. The results show, for the first time at a volcano, that a very peculiar mechanism is generating episodes of unrest, involving mainly activation of the geothermal system from deeper magma reservoirs. The character and evolution of unrest episodes is strongly controlled by structural features, like the ring-fault system at the borders of the caldera collapse. The use of detailed volcanological, mathematical and statistical procedures also make it possible to obtain a detailed picture of eruptive hazards in the whole Neapolitan area. The complex behaviour of this caldera, involving interaction between magmatic and geothermal phenomena, sheds light on the dynamics of the most dangerous types of volcanoes in the world.