An experimental study of amphibole stability in low-pressure granitic magmas and a revised Al-in-hornblende geobarometer

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$^\text{tot}$) of the synthetic calcic amphiboles varies systematically with pressure ($\small \textit{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$^\text{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 $\small \textit{P}$ (in kbar) from Al$^\text{tot}$ (in atoms per formula unit, apfu) according to the expression: $P$ (kbar) = 0.5 + 0.331(8) × Al$^\text{tot}$ + 0.995(4) × (Al$^\text{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$^\text{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.BHP Billiton, Royal Society (Wolfson Research Merit Award), California Institute of Technology (Moore Scholarship

Role of f02 in fluid saturation of oceanic basalt.

Assessing the conditions under which magmas become fluid-saturated has important bearings on the geochemical modelling of magmas because volatile exsolution may profoundly alter the behaviour of certain trace elements that are strongly partitioned in the coexisting fluid1. Saal et al.2 report primitive melt inclusions from dredged oceanic basalts of the Siqueiros transform fault, from which they derive volatile abundances of the depleted mantle, based on the demonstration that magmas are not fluid-saturated at their eruption depth and so preserve the mantle signature in terms of their volatile contents. However, in their analysis, Saal et al.2 consider only fluid−melt equilibria, and do not take into account the homogeneous equilibria between fluid species, which, as we show here, may lead to a significant underestimation of the pressure depth of fluid saturation

Quantifying garnet-melt trace element partitioning using lattice-strain theory: New crystal-chemical and thermodynamic constraints

Many geochemical models of major igneous differentiation events on the Earth, the Moon, and Mars invoke the presence of garnet or its high-pressure majoritic equivalent as a residual phase, based on its ability to fractionate critical trace element pairs (Lu/Hf, U/Th, heavy REE/light REE). As a result, quantitative descriptions of mid-ocean ridge and hot spot magmatism, and lunar, martian, and terrestrial magma oceans require knowledge of garnet-melt partition coefficients over a wide range of conditions. In this contribution, we present new crystal-chemical and thermodynamic constraints on the partitioning of rare earth elements (REE), Y and Sc between garnet and anhydrous silicate melt as a function of pressure (P), temperature (T), and composition (X). Our approach is based on the interpretation of experimentally determined values of partition coefficients D using lattice-strain theory. In this and a companion paper (Draper and van Westrenen this issue) we derive new predictive equations for the ideal ionic radius of the dodecahedral garnet X-site,

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melts: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO-MgO-A

Amphibole and apatite insights into the evolution and mass balance of Cl and S in magmas associated with porphyry copper deposits

Chlorine and sulfur are of paramount importance for supporting the transport and deposition of ore metals at magmatic–hydrothermal systems such as the Coroccohuayco Fe–Cu–Au porphyry–skarn deposit, Peru. Here, we used recent partitioning models to determine the Cl and S concentration of the melts from the Coroccohuayco magmatic suite using apatite and amphibole chemical analyses. The pre-mineralization gabbrodiorite complex hosts S-poor apatite, while the syn- and post-ore dacitic porphyries host S-rich apatite. Our apatite data on the Coroccohuayco magmatic suite are consistent with an increasing oxygen fugacity (from the gabbrodiorite complex to the porphyries) causing the dominant sulfur species to shift from S2− to S6+ at upper crustal pressure where the magmas were emplaced. We suggest that this change in sulfur speciation could have favored S degassing, rather than its sequestration in magmatic sulfides. Using available partitioning models for apatite from the porphyries, pre-degassing S melt concentration was 20–200 ppm. Estimates of absolute magmatic Cl concentrations using amphibole and apatite gave highly contrasting results. Cl melt concentrations obtained from apatite (0.60 wt% for the gabbrodiorite complex; 0.2–0.3 wt% for the porphyries) seems much more reasonable than those obtained from amphibole which are very low (0.37 wt% for the gabbrodiorite complex; 0.10 wt% for the porphyries). In turn, relative variations of the Cl melt concentrations obtained from amphibole during magma cooling are compatible with previous petrological constraints on the Coroccohuayco magmatic suite. This confirms that the gabbrodioritic magma was initially fluid undersaturated upon emplacement, and that magmatic fluid exsolution of the gabbrodiorite and the pluton rooting the porphyry stocks and dikes were emplaced and degassed at 100–200 MPa. Finally, mass balance constraints on S, Cu and Cl were used to estimate the minimum volume of magma required to form the Coroccohuayco deposit. These three estimates are remarkably consistent among each other (ca. 100 km3) and suggest that the Cl melt concentration is at least as critical as that of Cu and S to form an economic mineralization

The evolution and storage of primitive melts in the Eastern Volcanic Zone of Iceland: the 10 ka Grímsvötn tephra series (i.e. the Saksunarvatn ash)

Major, trace and volatile elements were measured in a suite of primitive macrocrysts and melt inclusions from the thickest layer of the 10 ka Grímsvötn tephra series (i.e. Saksunarvatn ash) at Lake Hvítárvatn in central Iceland. In the absence of primitive tholeiitic eruptions (MgO > 7 wt.%) within the Eastern Volcanic Zone (EVZ) of Iceland, these crystal and inclusion compositions provide an important insight into magmatic processes in this volcanically productive region. Matrix glass compositions show strong similarities with glass compositions from the AD 1783–84 Laki eruption, confirming the affinity of the tephra series with the Grímsvötn volcanic system. Macrocrysts can be divided into a primitive assemblage of zoned macrocryst cores (An_78–An_92, Mg#_cpx = 82–87, Fo_79.5–Fo_87) and an evolved assemblage consisting of unzoned macrocrysts and the rims of zoned macrocrysts (An_60–An_68, Mg#_cpx = 71–78, Fo_70–Fo_76). Although the evolved assemblage is close to being in equilibrium with the matrix glass, trace element disequilibrium between primitive and evolved assemblages indicates that they were derived from different distributions of mantle melt compositions. Juxtaposition of disequilibrium assemblages probably occurred during disaggregation of incompatible trace element-depleted mushes (mean La/Yb_melt = 2.1) into aphyric and incompatible trace element-enriched liquids (La/Yb_melt = 3.6) shortly before the growth of the evolved macrocryst assemblage. Post-entrapment modification of plagioclase-hosted melt inclusions has been minimal and high-Mg# inclusions record differentiation and mixing of compositionally variable mantle melts that are amongst the most primitive liquids known from the EVZ. Coupled high field strength element (HFSE) depletion and incompatible trace element enrichment in a subset of primitive plagioclase-hosted melt inclusions can be accounted for by inclusion formation following plagioclase dissolution driven by interaction with plagioclase-undersaturated melts. Thermobarometric calculations indicate that final crystal-melt equilibration within the evolved assemblage occurred at ~1140°C and 0.0–1.5 kbar. Considering the large volume of the erupted tephra and textural evidence for rapid crystallisation of the evolved assemblage, 0.0–1.5 kbar is considered unlikely to represent a pressure of long-term magma accumulation and storage. Multiple thermometers indicate that the primitive assemblage crystallised at high temperatures of 1240–1300°C. Different barometers, however, return markedly different crystallisation depth estimates. Raw clinopyroxene-melt pressures of 5.5–7.5 kbar conflict with apparent melt inclusion entrapment pressures of 1.4 kbar. After applying a correction derived from published experimental data, clinopyroxene-melt equilibria return mid-crustal pressures of 4±1.5 kbar, which are consistent with pressures estimated from the major element content of primitive melt inclusions. Long-term storage of primitive magmas in the mid-crust implies that low CO_2 concentrations measured in primitive plagioclase-hosted inclusions (262–800 ppm) result from post-entrapment CO_2 loss during transport through the shallow crust. In order to reconstruct basaltic plumbing system geometries from petrological data with greater confidence, mineral-melt equilibrium models require refinement at pressures of magma storage in Iceland. Further basalt phase equilibria experiments are thus needed within the crucial 1–7 kbar range.D.A.N. was supported by a Natural Environment Research Council studentship (NE/1528277/1) at the start of this project. SIMS analyses were supported by Natural Environment Research Council Ion Microprobe Facility award (IMF508/1013).This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s00410-015-1170-

The History, Relevance, and Applications of the Periodic System in Geochemistry

Geochemistry is a discipline in the earth sciences concerned with understanding the chemistry of the Earth and what that chemistry tells us about the processes that control the formation and evolution of Earth materials and the planet itself. The periodic table and the periodic system, as developed by Mendeleev and others in the nineteenth century, are as important in geochemistry as in other areas of chemistry. In fact, systemisation of the myriad of observations that geochemists make is perhaps even more important in this branch of chemistry, given the huge variability in the nature of Earth materials – from the Fe-rich core, through the silicate-dominated mantle and crust, to the volatile-rich ocean and atmosphere. This systemisation started in the eighteenth century, when geochemistry did not yet exist as a separate pursuit in itself. Mineralogy, one of the disciplines that eventually became geochemistry, was central to the discovery of the elements, and nineteenth-century mineralogists played a key role in this endeavour. Early “geochemists” continued this systemisation effort into the twentieth century, particularly highlighted in the career of V.M. Goldschmidt. The focus of the modern discipline of geochemistry has moved well beyond classification, in order to invert the information held in the properties of elements across the periodic table and their distribution across Earth and planetary materials, to learn about the physicochemical processes that shaped the Earth and other planets, on all scales. We illustrate this approach with key examples, those rooted in the patterns inherent in the periodic law as well as those that exploit concepts that only became familiar after Mendeleev, such as stable and radiogenic isotopes

Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

The 2010 eruption of Merapi (VEI 4) was the volcano’s largest since 1872. In contrast to the prolonged and effusive dome-forming eruptions typical of Merapi’s recent activity, the 2010 eruption began explosively, before a new dome was rapidly emplaced. This new dome was subsequently destroyed by explosions, generating pyroclastic density currents (PDCs), predominantly consisting of dark coloured, dense blocks of basaltic andesite dome lava. A shift towards open-vent conditions in the later stages of the eruption culminated in multiple explosions and the generation of PDCs with conspicuous grey scoria and white pumice clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxene-hosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H2O, CO2) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have re-equilibrated during shallower storage and/or ascent, at depths of ~0.6–9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform “buffered” Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km3 in order to produce the mass of SO2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing