We report new experimental data on the composition of magmatic amphiboles synthesised from a variety of granite (sensu lato) bulk compositions at near-solidus temperatures and pressures of ...0.8–10 kbar. The total aluminium content (Al
tot
) of the synthetic calcic amphiboles varies systematically with pressure (
P
), although the relationship is nonlinear at low pressures (<2.5 kbar). At higher pressures, the relationship resembles that of other experimental studies, which suggests of a general relationship between Al
tot
and P that is relatively insensitive to bulk composition. We have developed a new Al-in-hornblende geobarometer that is applicable to granitic rocks with the low-variance mineral assemblage: amphibole + plagioclase (An
15–80
) + biotite + quartz + alkali feldspar + ilmenite/titanite + magnetite + apatite. Amphibole analyses should be taken from the rims of grains, in contact with plagioclase and in apparent textural equilibrium with the rest of the mineral assemblage at temperatures close to the haplogranite solidus (725 ± 75 °C), as determined from amphibole–plagioclase thermometry. Mean amphibole rim compositions that meet these criteria can then be used to calculate
P
(in kbar) from Al
tot
(in atoms per formula unit, apfu) according to the expression:
P
kbar
=
0.5
+
0.331
8
×
Al
tot
+
0.995
4
×
Al
tot
2
This expression recovers equilibration pressures of our calibrant dataset, comprising both new and published experimental and natural data, to within ±16 % relative uncertainty. An uncertainty of 10 % relative for a typical Al
tot
value of 1.5 apfu translates to an uncertainty in pressure estimate of 0.5 kbar, or 15 % relative. Thus the accuracy of the barometer expression is comparable to the precision with which near-solidus amphibole rim composition can be characterised.
Formation and dynamics of magma reservoirs Sparks, R S J; Annen, C; Blundy, J D ...
Philosophical transactions - Royal Society. Mathematical, Physical and engineering sciences/Philosophical transactions - Royal Society. Mathematical, physical and engineering sciences,
02/2019, Volume:
377, Issue:
2139
Journal Article
Peer reviewed
Open access
The emerging concept of a magma reservoir is one in which regions containing melt extend from the source of magma generation to the surface. The reservoir may contain regions of very low fraction ...intergranular melt, partially molten rock (mush) and melt lenses (or magma chambers) containing high melt fraction eruptible magma, as well as pockets of exsolved magmatic fluids. The various parts of the system may be separated by a sub-solidus rock or be connected and continuous. Magma reservoirs and their wall rocks span a vast array of rheological properties, covering as much as 25 orders of magnitude from high viscosity, sub-solidus crustal rocks to magmatic fluids. Time scales of processes within magma reservoirs range from very slow melt and fluid segregation within mush and magma chambers and deformation of surrounding host rocks to very rapid development of magma and fluid instability, transport and eruption. Developing a comprehensive model of these systems is a grand challenge that will require close collaboration between modellers, geophysicists, geochemists, geologists, volcanologists and petrologists. This article is part of the Theo Murphy meeting issue 'Magma reservoir architecture and dynamics'.
Investigations of a variety of continental rifts and margins worldwide have revealed that a considerable volume of melt can intrude into the crust during continental breakup, modifying its ...composition and thermal structure. However, it is unclear whether the cause of voluminous melt production at volcanic rifts is primarily increased mantle temperature or plate thinning. Also disputed is the extent to which plate stretching or thinning is uniform or varies with depth with the entire continental lithospheric mantle potentially being removed before plate rupture. Here we show that the extensive magmatism during rifting along the southern Red Sea rift in Afar, a unique region of sub-aerial transition from continental to oceanic rifting, is driven by deep melting of hotter-than-normal asthenosphere. Petrogenetic modelling shows that melts are predominantly generated at depths greater than 80 kilometres, implying the existence of a thick upper thermo-mechanical boundary layer in a rift system approaching the point of plate rupture. Numerical modelling of rift development shows that when breakup occurs at the slow extension rates observed in Afar, the survival of a thick plate is an inevitable consequence of conductive cooling of the lithosphere, even when the underlying asthenosphere is hot. Sustained magmatic activity during rifting in Afar thus requires persistently high mantle temperatures, which would allow melting at high pressure beneath the thick plate. If extensive plate thinning does occur during breakup it must do so abruptly at a late stage, immediately before the formation of the new ocean basin.
A model for the generation of intermediate and silicic igneous rocks is presented, based on experimental data and numerical modelling. The model is directed at subduction-related magmatism, but has ...general applicability to magmas generated in other plate tectonic settings, including continental rift zones. In the model mantle-derived hydrous basalts emplaced as a succession of sills into the lower crust generate a deep crustal hot zone. Numerical modelling of the hot zone shows that melts are generated from two distinct sources; partial crystallization of basalt sills to produce residual H2O-rich melts; and partial melting of pre-existing crustal rocks. Incubation times between the injection of the first sill and generation of residual melts from basalt crystallization are controlled by the initial geotherm, the magma input rate and the emplacement depth. After this incubation period, the melt fraction and composition of residual melts are controlled by the temperature of the crust into which the basalt is intruded. Heat and H2O transfer from the crystallizing basalt promote partial melting of the surrounding crust, which can include meta-sedimentary and meta-igneous basement rocks and earlier basalt intrusions. Mixing of residual and crustal partial melts leads to diversity in isotope and trace element chemistry. Hot zone melts are H2O-rich. Consequently, they have low viscosity and density, and can readily detach from their source and ascend rapidly. In the case of adiabatic ascent the magma attains a super-liquidus state, because of the relative slopes of the adiabat and the liquidus. This leads to resorption of any entrained crystals or country rock xenoliths. Crystallization begins only when the ascending magma intersects its H2O-saturated liquidus at shallow depths. Decompression and degassing are the driving forces behind crystallization, which takes place at shallow depth on timescales of decades or less. Degassing and crystallization at shallow depth lead to large increases in viscosity and stalling of the magma to form volcano-feeding magma chambers and shallow plutons. It is proposed that chemical diversity in arc magmas is largely acquired in the lower crust, whereas textural diversity is related to shallow-level crystallization.
We present the variation in trace element partition coefficients measured at the interface between rapidly cooled clinopyroxene crystals and co-existing melts. Results indicate that, as the cooling ...rate is increased, clinopyroxene crystals are progressively depleted in Si, Ca and Mg counterbalanced by enrichments in Al (mainly tetrahedral Al
iv
), Na and Ti. Partition coefficients (Ds) for rare earth elements (REE), high field strength elements (HFSE) and transition elements (TE) increase with increasing cooling rate, in response to clinopyroxene compositional variations. The entry of REE into the M2 site is facilitated by a coupled substitution where either Na substitutes for Ca on the M2 site or Al
iv
substitutes for Si in the tetrahedral site. The latter substitution reflects an increased ease of locally balancing the excess charge at M2 as the number of surrounding Al
iv
atoms increases. Due to the lower concentration of Ca in rapidly cooled clinopyroxenes, divalent large ion lithophile elements (LILE) on M2 decrease at the expense of monovalent cations. Conversely, higher concentrations of HFSE and TE on the M1 site are facilitated as the average charge on this site increases with the replacement of divalent-charged cations by Al
vi
. Although crystallization kinetics modify clinopyroxene composition, deviations from equilibrium partitioning are insufficient to change the tendency of a trace element to be compatible or incompatible. Consequently, there are regular relationships between ionic radius, valence of the trace element and D. At both equilibrium and cooling rate conditions, Ds for isovalent cations define parabola-like curves when plotted against ionic radius, consistent with the lattice strain model, demonstrating that the partitioning of trace elements is driven by charge balance mechanisms; cation substitution reactions can be treated in terms of the energetics of the various charge-imbalanced configurations.
Compositions of plagioclase‐melt pairs are commonly used to constrain temperatures (T), dissolved water contents (H2O) and pressures (P) of pre‐eruptive magma storage and transport. However, previous ...plagioclase‐based thermometers, hygrometers, and barometers can have significant errors, leading to imprecise reconstructions of conditions during plagioclase growth. Here, we explore whether we can refine existing plagioclase‐based hygrothermobarometers with either plagioclase‐melt or melt‐only chemistry (±T/H2O), calibrated using random forest machine learning on experimental petrology data (n = 1,152). We find that both the plagioclase‐melt and melt‐only models return similar cross‐validation root‐mean‐square errors (RMSEs), as the melt holds most of the P‐T‐H2O information rather than the plagioclase. T/H2O‐dependent melt models have test set RMSEs of 25°C, 0.70 wt.% and 76 MPa for temperature, H2O content and pressure, respectively, while T/H2O‐independent models have RMSEs of 38°C, 0.97 wt.% and 91 MPa. The melt thermometer and hygrometer are applicable to a wide range of plagioclase‐bearing melts at temperatures between 664 and 1355°C, and with H2O concentrations up to 11.2 wt.%, while the melt barometer is suitable for pressures of ≤500 MPa. An updated plagioclase‐melt equilibrium model has also been calibrated, allowing the equilibrium anorthite content to be predicted with an error of 5.8 mol%. The new P‐T‐H2O‐An models were applied to matrix glasses and melt inclusions from the 1980 Mount St Helens (USA) and 2014–2015 Holuhraun (Iceland) eruptions, corroborating previous independent estimates and observations. Models are available at https://github.com/kyra‐cutler/Plag‐saturated‐melt‐P‐T‐H2O‐An, enabling assessment of plagioclase‐melt equilibrium and characterization of last‐equilibrated P‐T‐H2O conditions of plagioclase‐saturated magmas.
Plain Language Summary
Thermobarometry and hygrometry are common methods for reconstructing magma crystallization conditions (pressure (P), temperature (T) and dissolved water content (H2O)) prior to eruption. Plagioclase is a ubiquitous mineral found in a wide range of volcanic rocks and is often used to estimate P‐T‐H2O conditions. Here, we use machine learning‐based regression to calibrate new models, based on either plagioclase‐melt or melt‐only chemistry, to test whether we can improve the existing range of plagioclase‐based thermobarometers and hygrometers. We also develop an updated model to determine the equilibrium composition of a plagioclase crystal with a given normalized melt composition. We find that the plagioclase‐melt and melt‐only thermobarometers and hygrometers return very similar model errors due to the melt holding nearly all the P‐T‐H2O information. The models can be applied to a wide range of plagioclase‐bearing melts except for the barometer, which is only appropriate for upper crustal pressures of ≤500 MPa (≤13.8 km depth). Models are available at https://github.com/kyra‐cutler/Plag‐saturated‐melt‐P‐T‐H2O‐An, enabling assessment of plagioclase‐melt equilibrium and characterization of last‐equilibrated P‐T‐H2O conditions of plagioclase‐saturated magmas.
Key Points
Investigation into whether existing plagioclase‐based hygrothermobarometers can be refined using random forest machine learning
The algorithm highlights that only melt composition is required to estimate intensive variables (P‐T‐H2O) of plagioclase‐saturated magmas
New melt models show improvement in errors compared to thermodynamic‐based plagioclase‐melt hygrothermobarometers and equilibria models
Volatiles contribute to magma ascent through the sub-volcanic plumbing system. Here, we investigate melt inclusion compositions in terms of major and trace elements, as well as volatiles (H
2
O, CO
2
..., SO
2
, F, Cl, Br, S) for Quaternary Plinian and dome-forming dacite and andesite eruptions in the central and the northern part of Dominica (Lesser Antilles arc). Melt inclusions, hosted in orthopyroxene, clinopyroxene and plagioclase are consistently rhyolitic. Post-entrapment crystallisation effects are limited, and negligible in orthopyroxene-hosted inclusions. Melt inclusions are among the most water-rich yet recorded (≤ 8 wt% H
2
O). CO
2
contents are generally low (< 650 ppm), although in general the highest pressure melt inclusion contain the highest CO
2
. Some low-pressure (< 3 kbars) inclusions have elevated CO
2
(up to 1100–1150 ppm), suggestive of fluxing of shallow magmas with CO
2
-rich fluids. CO
2
-trace element systematics indicate that melts were volatile-saturated at the time of entrapment and can be used for volatile-saturation barometry. The calculated pressure range (0.8–7.5 kbars) indicates that magmas originate from a vertically-extensive (3–27 km depth) storage zone within the crust that may extend to the sub-Dominica Moho (28 km). The vertically-extensive crustal system is consistent with mush models for sub-volcanic arc crust wherein mantle-derived mafic magmas undergo differentiation over a range of crustal depths. The other volatile range of composition for melt inclusions from the central part is F (75–557 ppm), Cl (1525–3137 ppm), Br (6.1–15.4 ppm) and SO
2
(< 140 ppm), and for the northern part it’s F (92–798 ppm), Cl (1506–4428 ppm), Br (not determined) and SO
2
(< 569; one value at 1015 ppm). All MIs, regardless of provenance, describe the same Cl/F correlation (8.3 ± 2.7), indicating that the magma source at depth is similar. The high H
2
O content of Dominica magmas has implications for hazard assessment.
Partial melting in the Earth's mantle plays an important part in generating the geochemical and isotopic diversity observed in volcanic rocks at the surface. Identifying the composition of these ...primary melts in the mantle is crucial for establishing links between mantle geochemical 'reservoirs' and fundamental geodynamic processes. Mineral inclusions in natural diamonds have provided a unique window into such deep mantle processes. Here we provide experimental and geochemical evidence that silicate mineral inclusions in diamonds from Juina, Brazil, crystallized from primary and evolved carbonatite melts in the mantle transition zone and deep upper mantle. The incompatible trace element abundances calculated for a melt coexisting with a calcium-titanium-silicate perovskite inclusion indicate deep melting of carbonated oceanic crust, probably at transition-zone depths. Further to perovskite, calcic-majorite garnet inclusions record crystallization in the deep upper mantle from an evolved melt that closely resembles estimates of primitive carbonatite on the basis of volcanic rocks. Small-degree melts of subducted crust can be viewed as agents of chemical mass-transfer in the upper mantle and transition zone, leaving a chemical imprint of ocean crust that can possibly endure for billions of years.
A model for the generation of intermediate and silicic igneous rocks is presented, based on experimental data and numerical modelling. The model is directed at subduction-related magmatism, but has ...general applicability to magmas generated in other plate tectonic settings, including continental rift zones. In the model mantle-derived hydrous basalts emplaced as a succession of sills into the lower crust generate a deep crustal hot zone. Numerical modelling of the hot zone shows that melts are generated from two distinct sources; partial crystallization of basalt sills to produce residual H2O-rich melts; and partial melting of pre-existing crustal rocks. Incubation times between the injection of the first sill and generation of residual melts from basalt crystallization are controlled by the initial geotherm, the magma input rate and the emplacement depth. After this incubation period, the melt fraction and composition of residual melts are controlled by the temperature of the crust into which the basalt is intruded. Heat and H2O transfer from the crystallizing basalt promote partial melting of the surrounding crust, which can include meta-sedimentary and meta-igneous basement rocks and earlier basalt intrusions. Mixing of residual and crustal partial melts leads to diversity in isotope and trace element chemistry. Hot zone melts are H2O-rich. Consequently, they have low viscosity and density, and can readily detach from their source and ascend rapidly. In the case of adiabatic ascent the magma attains a super-liquidus state, because of the relative slopes of the adiabat and the liquidus. This leads to resorption of any entrained crystals or country rock xenoliths. Crystallization begins only when the ascending magma intersects its H2O-saturated liquidus at shallow depths. Decompression and degassing are the driving forces behind crystallization, which takes place at shallow depth on timescales of decades or less. Degassing and crystallization at shallow depth lead to large increases in viscosity and stalling of the magma to form volcano-feeding magma chambers and shallow plutons. It is proposed that chemical diversity in arc magmas is largely acquired in the lower crust, whereas textural diversity is related to shallow-level crystallization.