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 gas-rich 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.Supplementary material:Videos showing the evolution of magma composition with time in the shallow chamber for simulations 1–4 are available at www.geolsoc.org.uk/SUP18819

Large-volume, high-crystallinity, chemically homogeneous ignimbrites, dubbed ‘monotonous intermediates’, provide a unique opportunity to investigate the evolution of crustal magmatic reservoirs. We present the results of hydrothermal experiments on a dacite from Fish Canyon Tuff (FCT) in Colorado (USA), a classic example of a monotonous intermediate deposit, in order to characterize the variations in chemical and physical properties of hydrous dacite magmas as a function of temperature. The experiments (200 MPa, 720–1100°C) span the inferred pre-eruptive conditions of FCT magmas, and are shown to provide the best match to the chemical and physical properties of the erupted magmas at 790±10°C under conditions at or close to water-saturation. The results show the important effect of water content in controlling the chemical and physical evolution of magma, and the contrasted behaviour of water-saturated v. water-undersaturated magmas. In both cases, however, there is a broad interval of temperature (200°C) over which crystal fraction changes little. By recasting this behaviour in terms of enthalpy, rather than temperature, as the independent variable we show that this interval corresponds to a minimum in the change in crystallinity per unit of energy added or subtracted from the system, such that small perturbations to the heat content of the system (e.g. by cooling or new magma injections) results in very little change in magma properties. The crystal content in this interval is 55–65 wt%, which is close to the phenocryst content (40–55 wt%) of monotonous intermediates. We propose that crystal-rich magmas tend to settle in this ‘petrological trap’, changing little in physical and chemical properties over time as the system grows. Petrological trapping enables very large volumes of intermediate magma to accumulate in the shallow crust until such time as the net buoyancy force of these crystal-rich magma is sufficient to overcome the strength of the roof rocks, leading to a potentially very large eruption.

Determining the timescales of magma degassing is essential for understanding the mechanisms controlling the eruption style and the dynamics of magmatic systems. Towards this end, we measured 210Pb–226Ra disequilibria in andesite lavas erupted from Volcán de Colima between 1998 and 2010. (210Pb/226Ra)0 activity ratios range from 0.86 to 1.09, and are best explained in terms of 222Rn degassing and accumulation. The range in 210Pb deficits indicates that the timescales of 222Rn degassing did not exceed 11 years. 210Pb excesses are rare and small (<10%), which signifies that 222Rn degassing is more effective than 210Pb accumulation in this intermediate system despite the relatively low gas output at the surface. The absence of significant 210Pb excesses strongly suggests that the volcanic activity results from episodic ascent of small magma batches through the vapour-saturated section of the magmatic system. Overall, the degassing models based on 210Pb–226Ra disequilibrium suggest an open and complex subvolcanic magmatic system comprising several conduits in which multiple magma batches reside for up to 10 years. Shifts from effusive to explosive Vulcanian eruptive phases are not related to changes in degassing mode on timescales resolvable using 210Pb–226Ra disequilibria.

Understanding the magmatic processes that drive unrest at silicic calderas remains a major goal in Volcanology. Rabaul in Papua New Guinea is an exceptional location because after two decades of unrest and a peak in seismicity and deformation in 1983–85, eruptive activity began in 1994 and is still ongoing. A particularly large sub-Plinian eruption occurred from Tavurvur in October 2006. Whole-rock compositions are andesitic and reflect mixing/mingling between basaltic and dacitic magmas from the same system. The magmas that fed the 2006 eruption were stored at about 930°C, with 1–3 wt% H2O, 25–520 ppm CO2, and 50–2500 ppm SO2 in the melt. Melt inclusions hosted in pyroxene, and plagioclase phenocrysts record fractional crystallization at ≤200 MPa under relatively dry and poorly oxidizing conditions. Magma mixing/mingling is expressed as heterogeneous glass compositions, strongly zoned phenocrysts, and mafic crystal aggregates. A textural maturation from fine, acicular to large, blocky crystal clots implies different relative ages of formation. Modelling the chemical zoning of plagioclase shows that mafic–silicic interactions started a couple of decades prior to the 2006 eruption and continued until days to weeks prior to eruption. Basaltic replenishments have been driving unrest and eruption at the Rabaul caldera since the 1970s.Supplementary material:Tables and figures reporting the composition of the Tavurvur 2006, Kombiu and 1.4 ka BP caldera samples and showing thermodynamic modelling with MELTS are available at http://www.geolsoc.org.uk/SUP18816