The invasion of active magma chambers by primitive magma of deeper provenance is a frequent occurrence in volcanic systems, and it is commonly associated with pressurization. Chamber replenishment is driven by pressure and buoyancy forces that cause magma ascent towards shallow depths. We examine the end‐member case of pure buoyancy‐driven (natural) convection in crustal reservoirs deriving from the presence of degassed, dense magma at shallow level, that can originate a gravitational instability. Space‐time‐dependent numerical simulations of magma dynamics in composite underground systems reveal highly nonlinear pressure evolution dominated by decompression at shallow depths. This counterintuitive result originates from the compressible nature of multiphase magmas and their complex convection and mixing dynamics. Shallow magma chamber decompression on replenishment is favored by large volatile contents of the uprising magma, resulting in large density contrasts among the resident and the incoming components. These results show that the intuitive concept of magma chamber pressurization upon replenishment may not always hold in real situations dominated by buoyancy, and provide new perspectives for the interpretation of geophysical records at active volcanoes.

Arrival of magma from depth into shallow reservoirs has been documented as one of the possible processes leading to eruption. Magma intruding and rising to the surface interacts with the already emplaced, degassed magmas residing at shallower depths, leaving chemical signatures in the erupted products.We performed two-dimensional numerical simulations of the arrival of gasrich magmas into shallow reservoirs. We solve the fluid dynamics for the two interacting magmas, evaluating the space-time evolution of the physical properties of the mixture.Convection and mingling develop quickly into the chamber and feeding conduit/dyke, leading on longer timescales to a density stratification with the lighter, gas-richer magma, mixed with different proportions of the resident magma, rising to the top of the chamber due to buoyancy. Over timescales of hours, the magmas in the reservoir appear to have mingled throughout, and convective patterns become harder to identify.Our simulations have been performed changing the geometry of the shallow reservoir and the gas content of the initial end-member magmas. Horizontally elongated magma chambers, as well as higher density contrasts between the two magmas, cause faster ascent velocities and also increase the mixing efficiency.

Many volcanic eruptions are shortly preceded by injection of new magma into a pre-existing, shallow (<10 km) magma chamber, causing convection and mixing between the incoming and resident magmas. These processes may trigger dyke propagation and further magma rise, inducing long-term (days to months) volcano deformation, seismic swarms, gravity anomalies, and changes in the composition of volcanic plumes and fumaroles, eventually culminating in an eruption. Although new magma injection into shallow magma chambers can lead to hazardous event, such injection is still not systematically detected and recognized. Here, we present the results of numerical simulations of magma convection and mixing in geometrically complex magmatic systems, and describe the multiparametric dynamics associated with buoyant magma injection. Our results reveal unexpected pressure trends and pressure oscillations in the Ultra-Long-Period (ULP) range of minutes, related to the generation of discrete plumes of rising magma. Very long pressure oscillation wavelengths translate into comparably ULP ground displacements with amplitudes of order 10 −4 -10 −2 m. Thus, new magma injection into magma chambers beneath volcanoes can be revealed by ULP ground displacement measured at the surface.

The most violent and catastrophic volcanic eruptions on Earth have been triggered by the refilling of a felsic volcanic magma chamber by a hotter more mafic magma. Examples include Vesuvius 79 AD, Krakatau 1883, Pinatubo 1991, and Eyjafjallajökull 2010. Since the first hypothesis, plenty of evidence of magma mixing processes, in all tectonic environments, has accumulated in the literature allowing this natural process to be defined as fundamental petrological processes playing a role in triggering volcanic eruptions, and in the generation of the compositional variability of igneous rocks. Combined with petrographic, mineral chemistry and geochemical investigations, isotopic analyses on volcanic rocks have revealed compositional variations at different length scales pointing to a complex interplay of fractional crystallization, mixing/mingling and crustal contamination during the evolution of several magmatic feeding systems. But to fully understand the dynamics of mixing and mingling

[1] Volcanic eruptions are often accompanied by spatiotemporal migration of ground deformation, a consequence of pressure changes within magma reservoirs and pathways. We modeled the propagation of pressure variations through the east rift zone (ERZ) of K" ılauea Volcano, Hawai'i, caused by magma withdrawal during the early eruptive episodes (1983)(1984)(1985) of the ongoing Pu'u ' " O'" o-Kupaianaha eruption. Eruptive activity at the Pu'u ' " O'" o vent was typically accompanied by abrupt deflation that lasted for several hours and was followed by a sudden onset of gradual inflation once the eruptive episode had ended. Similar patterns of deflation and inflation were recorded at K" ılauea's summit, approximately 15 km to the northwest, albeit with time delays of hours. These delay times can be reproduced by modeling the spatiotemporal changes in magma pressure and flow rate within an elastic-walled dike that traverses K" ılauea's ERZ. Key parameters that affect the behavior of the magma-dike system are the dike dimensions, the elasticity of the wall rock, the magma viscosity, and to a lesser degree the magnitude and duration of the pressure variations themselves. Combinations of these parameters define a transport efficiency and a pressure diffusivity, which vary somewhat from episode to episode, resulting in variations in delay times. The observed variations in transport efficiency are most easily explained by small, localized changes to the geometry of the magma pathway.