Explore the vast range of titles available.

This introductory article provides a synopsis of our current understanding of the form and dynamics of magma reservoirs in the crust. This knowledge is based on a range of experimental, observational and theoretical approaches, some of which are multidisclipinary and pioneering. We introduce and provide a contextual background for the papers in this issue, which cover a wide range of topics, encompassing magma storage, transport, behaviour and rheology, as well as the timescales on which magma reservoirs operate. We summarize the key findings that emerged from the meeting and the challenges that remain. The study of magma reservoirs has wide implications not only for understanding geothermal and magmatic systems, but also for natural oil and gas reservoirs and for ore deposit formation. This article is part of the Theo Murphy meeting issue ‘Magma reservoir architecture and dynamics’.

There are clear microstructural differences between mafic plutonic rocks that formed in a dynamic liquid-rich environment, in which crystals can be moved and re-arranged by magmatic currents, and those in which crystal nucleation and growth are essentially in situ and static. Crystalline enclaves, derived from deep crustal mushy zones and erupted in many volcanic settings, afford a unique opportunity to use the understanding of microstructural development, established from the study of intrusive plutons, to place constraints on the architecture of sub-volcanic systems. Here, we review the relevant microstructural literature, before applying these techniques to interrogate the crystallization environments of enclaves from the Kameni Islands of Santorini and Rábida Volcano in the Galápagos. Crystals in samples of deep-sourced material from both case studies preserve evidence of at least some time spent in a liquid-rich environment. The Kameni enclaves appear to record an early stage of crystallization during which crystals were free to move, with the bulk of crystallization occurring in a static, mushy environment. By contrast, the Rábida enclaves were sourced from an environment in which hydrodynamic sorting and re-arrangement by magmatic currents were common, consistent with a liquid-rich magma chamber. While presently active volcanoes are thought to be underlain by extensive regions rich in crystal mush, these examples preserve robust evidence for the presence of liquid-rich magma chambers in the geological record. This article is part of the Theo Murphy meeting issue ‘Magma reservoir architecture and dynamics'.

Only a small fraction (approx. 1–20%) of magmas generated in the mantle erupt at the surface. While volcanic eruptions are typically considered as the main exhaust pipes for volatile elements to escape into the atmosphere, the contribution of magma reservoirs crystallizing in the crust is likely to dominate the volatile transfer from depth to the surface. Here, we use multiscale physical modelling to identify and quantify the main mechanisms of gas escape from crystallizing magma bodies. We show that most of the outgassing occurs at intermediate to high crystal fraction, when the system has reached a mature mush state. It is particularly true for shallow volatile-rich systems that tend to exsolve volatiles through second boiling, leading to efficient construction of gas channels as soon as the crystallinity reaches approximately 40–50 vol.%. We, therefore, argue that estimates of volatile budgets based on volcanic activity may be misleading because they tend to significantly underestimate the magmatic volatile flux and can provide biased volatile compositions. Recognition of the compositional signature and volumetric dominance of intrusive outgassing is, therefore, necessary to build robust models of volatile recycling between the mantle and the surface. This article is part of the Theo Murphy meeting issue ‘Magma reservoir architecture and dynamics’.

Crystal mush is rapidly emerging as a new paradigm for the evolution of igneous systems. Mid-ocean ridges provide a unique opportunity to study mush processes: geophysical data indicate that, even at the most magmatically robust fast-spreading ridges, the magma plumbing system typically comprises crystal mush. In this paper, we describe some of the consequences of crystal mush for the evolution of the mid-ocean ridge magmatic system. One of these is that melt migration by porous flow plays an important role, in addition to rapid, channelized flow. Facilitated by both buoyancy and (deformation-enhanced) compaction, porous flow leads to reactions between the mush and migrating melts. Reactions between melt and the surrounding crystal framework are also likely to occur upon emplacement of primitive melts into the mush. Furthermore, replenishment facilitates mixing between the replenishing melt and interstitial melts of the mush. Hence, crystal mushes facilitate reaction and mixing, which leads to significant homogenization, and which may account for the geochemical systematics of mid-ocean ridge basalt (MORB). A second consequence is cryptic fractionation. At mid-ocean ridges, a plagioclase framework may already have formed when clinopyroxene saturates. As a result, clinopyroxene phenocrysts are rare, despite the fact that the vast majority of MORB records clinopyroxene fractionation. Hence, melts extracted from crystal mush may show a cryptic fractionation signature. Another consequence of a mush-dominated plumbing system is that channelized flow of melts through the crystal mush leads to the occurrence of vertical magmatic fabrics in oceanic gabbros, as well as the entrainment of diverse populations of phenocrysts. Overall, we conclude that the occurrence of crystal mush has a number of fundamental implications for the behaviour and evolution of magmatic systems, and that mid-ocean ridges can serve as a useful template for trans-crustal mush columns elsewhere. This article is part of the Theo Murphy meeting issue ‘Magma reservoir architecture and dynamics'.

Magma reservoirs play a key role in controlling numerous processes in planetary evolution, including igneous differentiation and degassing, crustal construction, and volcanism. For decades, scientists have tried to understand what happens in these reservoirs, using an array of techniques such as field mapping/petrology/geochemistry/geochronology on plutonic and volcanic lithologies, geophysical imaging of active magmatic provinces, and numerical/experimental modeling. This review paper tries to follow this multi-disciplinary framework while discussing past and present ideas. We specifically focus on recent claims that magma columns within the Earth's crust are mostly kept at high crystallinity (\"mush zones\"), and that the dynamics within those mush columns, albeit modulated by external factors (e.g., regional stress field, rheology of the crust, pre-existing tectonic structure), play an important role in controlling how magmas evolve, degas, and ultimately erupt. More specifically, we consider how the chemical and dynamical evolution of magma in dominantly mushy reservoirs provides a framework to understand: (1) the origin of petrological gradients within deposits from large volcanic eruptions (\"ignimbrites\"); (2) the link between volcanic and plutonic lithologies; (3) chemical fractionation of magmas within the upper layers of our planet, including compositional gaps noticed a century ago in volcanic series (4) volatile migration and storage within mush columns; and (5) the occurrence of petrological cycles associated with caldera-forming events in long-lived magmatic provinces. The recent advances in understanding the inner workings of silicic magmatism are paving the way to exciting future discoveries, which, we argue, will come from interdisciplinary studies involving more quantitative approaches to study the crust-reservoir thermo-mechanical coupling as well as the kinetics that govern these open systems.

One of the outstanding problems in understanding the behavior of intermediate-to-silicic magmatic systems is the mechanism(s) by which large volumes of crystal-poor rhyolite can be extracted from crystal-rich mushy storage zones in the mid-deep crust. The mechanisms commonly invoked are hindered settling, micro-settling, and compaction. The concept of micro-settling involves extraction of grains from a crystal framework during Ostwald ripening and has been shown to be non-viable in the metallic systems for which it was originally proposed. Micro-settling is also likely to be insignificant in silicic mushes, because ripening rates are slow for quartz and plagioclase, contact areas between grains in a crystal mush are likely to be large, and abundant low-angle grain boundaries promote grain coalescence rather than ripening. Published calculations of melt segregation rates by hindered settling (Stokes settling in a crystal-rich system) neglect all but fluid dynamical interactions between particles. Because tabular silicate minerals are likely to form open, mechanically coherent, frameworks at porosities as high as ~ 75%, settling of single crystals is only likely in very melt-rich systems. Gravitationally-driven viscous compaction requires deformation of crystals by either dissolution–reprecipitation or dislocation creep. There is, as yet, no reported microstructural evidence of extensive, syn-magmatic, internally-generated, viscous deformation in fully solidified silicic plutonic rocks. If subsequent directed searches do not reveal clear evidence for internally-generated buoyancy-driven melt segregation processes, it is likely that other factors, such as rejuvenation by magma replenishment, gas filter-pressing, or externally-imposed stress during regional deformation, are required to segregate large volumes of crystal-poor rhyolitic liquids from crustal mushy zones.

Textures and compositions of minerals can be used to infer the physiochemical conditions present within magmatic systems. Given that plagioclase is an abundant phase in many magmatic systems, understanding the link between texture and process is vital. Here, we present a database of textural and compositional data for > 1800 plagioclase crystals in mid-ocean ridge basalt from the Gakkel Ridge (Arctic Ocean) to investigate the physiochemical conditions and processes that govern the formation of plagioclase textures and compositions. The Gakkel basalts have high modal crystal contents (up to 50%). The crystal cargo is complex, with both individual plagioclase and glomerocrysts showing large variations in crystal habit, zoning and resorption. The most common types of zoning are reverse and patchy; we attribute patchy zoning to infilling following either skeletal growth or resorption. Resorption is abundant, with multiple resorption events commonly present in a single crystal, and results from both magmatic recharge and decompression. Periods of strong undercooling, distinct to quench crystallisation, are indicated by matured skeletal crystals and thin normally zoned melt inclusion-rich bands following resorption. Individual samples often contain diverse textural and compositional plagioclase groups. Furthermore, most plagioclase is not in equilibrium with its host melt. Finally, the porous open structures of some glomerocrysts suggest that they represent pieces of entrained disaggregated mush. We interpret this to indicate that the crystal cargo is not generally phenocrystic in origin. Instead, plagioclase crystals that formed in different parts of a mush-dominated plumbing system were entrained into ascending melts. The textures of individual crystals are a function of their respective histories of (under)cooling, magma mixing and decompression. The morphologies of melt inclusion trapped in the plagioclase crystals are associated with specific host crystal textures, suggesting a link between plagioclase crystallisation processes and melt inclusion entrapment. The database of plagioclase presented herein may serve as a template for the interpretation of plagioclase textures in magmatic systems elsewhere.

Hafnium, U, Th, and REE content of zircons from the Spirit Mountain batholith in southern Nevada correlate with calculated temperatures from the Ti-in-zircon thermometer to support field and petrologic evidence of rejuvenation of crystal mush and melt extraction events during the 2-million year accumulation of the granitoid batholith. Marked variation in zircon composition from sample to sample, from grain to grain within individual samples, and from zone to zone within individual grains documents in detail a history of fluctuating conditions with repeated episodes of replenishment, reheating, crystal mush rejuvenation, fractional crystallization, and melt segregation. The zircons exhibit compositional and thermal variability indicative of variations in host melt composition due to (1) melt rejuvenation, mixing, and fractionation (2) coeval growth of other REE-rich accessory minerals, and possibly (3) fluctuation in f O 2 .

Percolation through magma mush is a key transport mechanism for melts in the crust and is influenced by the permeability of the crystal framework. Existing models for mush permeability do not account for the range of microstructures that can evolve as mushes crystallize or compact to low melt fractions. Here, we use numerically generated domains of cuboids at the random maximum packing as a starting geometry for a loose magma mush. We then expand the cuboid edges into the pore spaces sequentially, representing a geometrical simulation of crystal overgrowth and crystallization. At each iterative step, we measure the melt fraction, specific surface area, and melt permeability via 3D fluid flow simulations. We find that (a) the permeability drops proportional to the drop in surface area as the melt fraction reduces, (b) the permeability falls to zero at a percolation threshold ϕc=0.0187±0.0010 ${\\phi }_{c}=0.0187\\mathit{\\pm }0.0010$ that is independent of scale and insensitive to the starting cuboid geometry, and (c) once the percolation threshold is determined, our data match a universal percolation model without requiring any free fitting parameters. We show how this percolation model accounts for any 3D shape of the crystals that comprise the evolving mush. Importantly, this approach demonstrates that mush permeability can remain non‐zero in texturally unequilibrated mushes, down to very low melt fractions. Our model outperforms previous models, which overestimate mush permeability by up to three orders of magnitude, and our model can be used to accurately predict how mush permeability changes as mushes mature and crystallize, with implications for quantifying melt extraction, percolation rates, and melt reservoir assembly.

Global Synthetic Aperture Radar (SAR) measurements made over the past decades provide insights into the lateral extent of magmatic domains, and capture volcanic process on scales useful for volcano monitoring. Satellite-based SAR imagery has great potential for monitoring topographic change, the distribution of eruptive products and surface displacements (InSAR) at subaerial volcanoes. However, there are challenges in applying it routinely, as would be required for the reliable operational assessment of hazard. The deformation detectable depends upon satellite repeat time and swath widths, relative to the spatial and temporal scales of volcanological processes. We describe the characteristics of InSAR-measured volcano deformation over the past two decades, highlighting both the technique’s capabilities and its limitations as a monitoring tool. To achieve this, we draw on two global datasets of volcano deformation: the Smithsonian Institution Volcanoes of the World database and the Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics volcano deformation catalogue, as well as compiling some measurement characteristics and interpretations from the primary literature. We find that a higher proportion of InSAR observations capture non-eruptive and non-magmatic processes than those from ground-based instrument networks, and that both transient (< month) and long-duration (> 5 years) deformation episodes are under-represented. However, satellite radar is already used to assess the development of extended periods of unrest and long-lasting eruptions, and improved spatial resolution and coverage have resulted in the detection of previously unrecognised deformation at both ends of the spatial scale (~ 10 to > 1000 km 2 ). ‘Baseline’ records of past InSAR measurements, including ‘null’ results, are fundamental for any future interpretation of interferograms in terms of hazard‚ both by providing information about past deformation at an individual volcano, and for assessing the characteristics of deformation that are likely to be detectable (and undetectable) using InSAR. More than half of all InSAR deformation signals attributed to magmatic processes have sources in the shallow crust (< 5 km depth). While the depth distribution of InSAR-derived deformation sources is affected by measurement limitations, their lateral distribution provides information about the extent of active magmatic domains. Deformation is common (24% of all potentially magmatic events) at loci ≥5 km away from the nearest active volcanic vent. This demonstrates that laterally extensive active magmatic domains are not exceptional, but can comprise the shallowest part of trans-crustal magmatic systems in a range of volcanic settings.