www.nature.com/scientificreports

OPEN

received: 19 April 2016 accepted: 07 July 2016 Published: 04 August 2016

Boron isotope fractionation in magma via crustal carbonate dissolution Frances M. Deegan1,2, Valentin R. Troll1,3, Martin J. Whitehouse2, Ester M. Jolis1 & Carmela Freda3 Carbon dioxide released by arc volcanoes is widely considered to originate from the mantle and from subducted sediments. Fluids released from upper arc carbonates, however, have recently been proposed to help modulate arc CO2 fluxes. Here we use boron as a tracer, which substitutes for carbon in limestone, to further investigate crustal carbonate degassing in volcanic arcs. We performed laboratory experiments replicating limestone assimilation into magma at crustal pressure-temperature conditions and analysed boron isotope ratios in the resulting experimental glasses. Limestone dissolution and assimilation generates CaO-enriched glass near the reaction site and a CO2-dominated vapour phase. The CaO-rich glasses have extremely low δ11B values down to −41.5‰, reflecting preferential partitioning of 10B into the assimilating melt. Loss of 11B from the reaction site occurs via the CO2 vapour phase generated during carbonate dissolution, which transports 11B away from the reaction site as a boron-rich fluid phase. Our results demonstrate the efficacy of boron isotope fractionation during crustal carbonate assimilation and suggest that low δ11B melt values in arc magmas could flag shallowlevel additions to the subduction cycle.

Carbon and boron cycles

Carbon is transferred from the Earth’s interior to the surface mainly by CO2 degassing at arc volcanoes1,2. Arc CO2 emissions represent a mixture from several volatile sources including mantle carbon, subducted carbon from altered ocean crust or oceanic sediments, and possibly also carbon released from the over-riding plate during final magma ascent1–4. Carbonate is the main CO2-carrying mineral phase and its stability at sub-arc depths is debated. Some authors argue that carbonate is unlikely to efficiently break-down under normal subduction conditions5, while others propose that carbonate can be efficiently recycled via e.g., interaction with aqueous fluids6,7. Carbon dioxide outgassing has varied considerably throughout Earth history‒​during the warmer Cretaceous period, for example, it has been calculated that there was as much as 220 to 550% CO2 outgassing relative to present day rates8,9. If this calculation is correct, then two important implications arise; (i) that intensified global continental arc activity may have led to excess degassing at paleo-arcs to explain paleo-climate fluctuations and (ii) that decarbonation reactions in the Earth’s crust may be an important, albeit time-variable, factor in modulating the carbon cycle (cf. refs 3, 4, 8–10). It is the latter aspect in particular that provides the incentive for this study, wherein we investigate how crustal decarbonation can potentially affect volatile systematics, and especially boron cycling, in subduction zones.

Boron as a tracer in subduction zones

Boron is a fluid-mobile trace element and its isotope ratios have been used to evaluate the involvement of dehydrated slab materials in arc magma-genesis for the past ca. 25 years (see ref. 11 for a review of current concepts in boron isotope systematics in subduction zones). Boron isotopes are well suited to this task because 11B is preferentially fractionated into co-existing fluids, giving rise to boron-enriched slab-derived fluids with high δ​ 11 B values that progressively evolve towards low δ​11B values due to continual dehydration reactions and boron loss along the down-going slab11–16. This subduction-driven isotopic evolution is mirrored in the boron isotope 1

Department of Earth Sciences, Centre for Experimental Mineralogy, Petrology, and Geochemistry (CEMPEG), Uppsala University, SE-752 36, Uppsala, Sweden. 2Department of Geosciences, Swedish Museum of Natural History, SE-104 05, Stockholm, Sweden. 3Istituto Nazionale di Geofisica e Vulcanologia (INGV), 00143 Rome, Italy. Correspondence and requests for materials should be addressed to F.M.D. (email: [email protected]) Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

1

www.nature.com/scientificreports/

Figure 1.  Overview of existing boron data and new data from this study. (a) Literature data for various geochemical reservoirs compared to our experiments. The full range of Mid-Ocean Ridge Basalt (MORB) and Ocean Island Basalt (OIB) literature data is shown and the δ​11B values suggested to represent the primary mantle are highlighted (after refs 41 and 42). (b) Literature δ​11B data for arc volcanoes. (c–e) A four component arc system comprising mantle, slab lithologies, slab fluids, and upper crustal carbonates is shown. Process arrows indicate the effects that shallow slab fluids versus assimilation of degassing crust would have on δ​11B values and B concentration of erupted products. The plotted experimental data correspond to limestonebearing experiments except for panel e) which shows dolostone experimental data. Error bars (1σ​) are smaller than symbol size. Full data sources are provided in Supplementary Table S-5. Abbreviations: TIMS, thermal ionisation mass spectrometry; SIMS, secondary ionisation mass spectrometry.

composition of exhumed rocks and minerals subjected to prograde blueschist facies metamorphism, which generally have negative δ​11B values due to loss of “heavy” boron into a C-O-H fluid phase (Fig. 1a; see also ref. 16). Thus, it is envisaged that the tendency for many arc magmas to display boron concentrations higher than typical Mid-Ocean Ridge Basalt (MORB) mantle is due to input of slab-derived boron and that isotopically diverse fluids produced during subduction can react with the mantle wedge to generate the wide range of δ​11B values recorded in global arc suites (e.g., ref. 11 and Fig. 1a–c). A controlling factor on the boron isotope systematics of arcs is likely related to the thermal structure of a particular subduction zone whereby boron-enriched, high δ​11B fluids will be generated at shallow depths along “hot” slabs and hence low δ​11B components are likely to be subducted to Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

2

www.nature.com/scientificreports/

Figure 2.  Experimental data. Back scattered electron (BSE) images and δ​11B profiles for representative (a,b) Vesuvius and (c,d) Merapi experiments. The solid red line on the BSE images (a,c) represents the SIMS traverse and the red symbols indicate analysis spots. In (b) and (d), average δ​11B values for different glass compositional domains are represented by red horizontal bars and the full range of starting material values measured is indicated by black dashed horizontal lines. The CaO-rich glasses have significantly lower δ​11B values than both the starting material and the CaO-normal glass due to transport of 11B away from the reaction site in the CO2 vapour. Error as in Fig. 1.

greater depth (e.g., ref. 17). In contrast, in “cooler” regimes, high δ​11B components may be preserved and transported to deeper parts of a subduction zone11. Although great progress has been made in understanding the subduction-scale systematics of boron, the processes that may serve to alter the boron isotope composition of magma after segregation from its source, such as crustal assimilation in the over-riding plate, have been addressed less frequently than source-related processes (some notable examples include refs 18–20). Because BO33− (borate) is thought to substitute for the CO32− site in carbonate21, and as carbonate has the potential to strongly degas in magma22–25, we aim to test if boron transport could be tied to carbonate degassing in active subduction zones where the upper plate contains CaCO3-bearing minerals. To do this, we offer high temperature-high pressure experimental simulations of carbonate assimilation at conditions corresponding to the over-riding plate and provide spatially controlled boron isotope analyses of the experimental products by Secondary Ionisation Mass Spectrometry (SIMS).

Results and Discussion

Our experiments were designed to simulate CO2 fluxing in magma at mid- to upper-crustal pressure conditions whereby solid carbonate (mainly CaCO3 but also MgCa(CO3)2 in some cases) was allowed to react with pre-fused, powdered magmatic rock from Mt. Merapi (Indonesia) and Mt. Vesuvius (Italy) at 1200 °C and 0.5 GPa for up to 300 s (refs 24 and 25). These volcanic systems were chosen because both are subduction-related and display evidence for crustal carbonate assimilation in form of erupted calc-silicate xenoliths and chemical signatures in erupted rocks and fumarole gas22,23,26,27. The advantage of our experiments is that they simulate short-term disequilibrium reactions in order to capture the temporal evolution of magma-carbonate interaction as time variable “snapshots”. Major element compositions of starting materials and experimental products are provided in Supplementary Tables S1 and S2. Boron data for the fused starting materials (n =​ 25) and the experimental products (n =​ 147) are provided in Supplementary Tables S3 and S4. Our experimental products comprise a CaO-normal glass similar in composition to the starting materials, a CaO-rich or MgO-rich glass, and a mixing interface between the two domains that shows variable CaO and MgO contents24,25. Incongruent break-down of carbonate produced free CO2 bubbles that permeated all melt domains28 (Fig. 2). The δ​11B values of the starting materials range from −​8.8 to −​3.5‰ for Merapi and from −​14.6 to −​7.6‰ for Vesuvius (Fig. 1a). In the Vesuvius case, the measured δ​11B values overlap with the lower end of the established Vesuvius range (−​7.6 to −​6.3‰; ref. 29), and are similar to literature δ​11B values for other Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

3

www.nature.com/scientificreports/ Italian magmatic systems (e.g., −​13.7‰ at Stromboli; ref. 30). The boron concentration of the Vesuvius starting glass ranges from 12 to 14 μ​g/g, and is hence close to the reported range of 14 to 36 μ​g/g for Vesuvius erupted products29. To the best of our knowledge, there are currently no published δ​11B data for Merapi. Boron concentration of our Merapi starting glass ranges from 15 to 18 μ​g/g, consistent with recorded Merapi whole-rock values of 12 to 20 μ​g/g (ref. 31), but slightly lower than reported concentrations for Merapi clinopyroxene-hosted melt inclusions (35 to 109 μ​g/g; ref. 32). Boron isotope profiles were analysed across the interface between CaO-rich and CaO-normal glass (Fig. 2). The δ​11B values and B concentration of the CaO-normal glasses range from −​5.3 to  +​ 1.2‰ and from 10 to 15 μ​g/g for Merapi and from −​14.7 to −​4.9‰ and 9 to 241 μ​g/g for Vesuvius. In contrast, the δ​11B values of the CaO-rich and (MgO)CaO-rich glasses range from −​21.9 to −​8.6‰ for Merapi and from −​41.5 to −​13.6‰ for Vesuvius. These values fall considerably below the δ​11B values of many subduction systems globally and they are significantly lower than the experimental starting materials (Fig. 1d,e; see also Supplementary Table S5 for sources of all literature data presented). The δ​11B values of the CaO-rich glasses are also considerably lower than those of carbonate in the literature. In general, biogenic carbonate has an average δ​11B value of +​19.1‰ and variable B concentration (Fig. 1). In contrast, lithified carbonate (i.e. lime/dolostone) has lower but still positive δ​11B values, which range from +​1.5 to +​8.4‰ and B concentration of 2 to 18 μ​g/g (refs 33 and 34; Fig. 1, data sources in Supplementary Table S5). We therefore conclude that simple binary mixing between the reactants cannot explain the low δ​11B melts in the experiments and implies additional processes at work. Boron concentration in the CaO-rich and CaO(MgO)-rich glasses ranges from 4 to 7 μ​g/g and 2 to 85 μg​ /g for Merapi and Vesuvius, respectively. A negative correlation between δ​11B value and B concentration is observed in all experiments (Fig. 1d,e) and hence the low δ​11B values recorded must be related to B degassing from the melt(s) under the experimental conditions (see Methods). We also note that the experimental data (n =​ 172) for the most part overlap the range of natural magmatic values (Fig. 1), supporting the view that the experiments provide a useful analogue for natural magmatic processes. We hence argue that boron isotope fractionation during decarbonation in the upper plate represents another means to generate low δ​11B values in arc magmas in addition to slab dehydration during subduction (cf. refs 3, 8, 9, 10 and 22). As our experiments were heated to 1200 °C over a period of six minutes (see Methods), the fractionation factor between tetrahedrally and trigonally coordinated boron, α​, varied with experiment duration, since α​is inversely proportional to temperature13. Realising that α​was not constant in our experiments precludes calculation of a unique Rayleigh model for our experimental data. To explain our lowest δ​11B values, α​probably varied from ca. 1.012 to 1.002 (the latter at 1200 °C), but these values have large uncertainties due to the fact that CO2 was continually fluxing in the experiments and that fractionation commenced already below 1200 °C. Nonetheless, a similar range of α​values may apply to situations where a large thermal gradient exists, such as across metamorphic aureoles in plutonic complexes, or along down-going slabs in subduction zones. Boron isotope fractionation in the experimental melts would also have been mirrored by evolving δ​11B values in the co-existing fluid, with δ​11B values becoming lower over time in a similar fashion to the slab dehydration models performed by ref. 15. Assuming that boron isotope fractionation is independent of pressure35, our data are in-line with these and similar models that predict values lower than −​30‰ in dehydrated subducted materials, and others as low as −​35‰ for phengite-free dehydrated assemblages at subduction depths15. The fundamental implication here is that extremely low δ​11B values can be generated in subducted material at ~3 GPa as well as in magma in the upper arc crust at ~0.5 GPa due to the presence of a coexisting fluid phase that serves to scavenge boron from the rock or silicate melt. Our experiments also reveal several dynamic aspects of boron transport in magma. In particular, the relatively elevated δ​11B values of some of the CaO-normal glass domains compared to the δ​11B values of the starting material (Fig. 2) leads us to propose a conceptual model (Fig. 3). In our model, carbonate dissolution and degassing at the onset of magma-carbonate interaction is the catalyst for boron isotope fractionation. The newly formed CO2 phase scavenges 11B from the carbonate and silicate melt, causing 11B(OH)3 to enter the newly generated volatile phase by substitution for CO2 and assimilation of 10B-rich material to occur at the decarbonation reaction site. The highly mobile volatile phase then rapidly migrates away from the reaction site, resulting in the generation of a fluid with a relatively high δ​11B value that progressively evolves towards lower δ​11B values, similar to some arc fluids and models thereof, as discussed above (Fig. 1c). Conversely, the relatively unaffected, CaO-normal melt further away from the carbonate dissolution site would be undersaturated with respect to CO2, which could facilitate coupled 11B and CO2 reabsorption in the melt and hence relatively high δ​11B melt values (Fig. 3). The extent to which this process is expressed as “low” or “high” δ​11B melts in individual volcanic systems depends on several factors, including the amount of carbonate assimilated, the viscosity of the magma and, particularly, the solubility of CO2 in the melt, since low melt solubility of CO2 will promote bubble formation and thus boron extraction from the co-existing melt. This process would be most effective under low pressures, since CO2 solubility in silicate melts decreases with decreasing pressure (see discussion in refs 24, 25, 28 and references therein), which would make boron extraction into a CO2-bearing phase most efficient in the upper parts of the crust. In conclusion, our data demonstrate that short time-scale (syn-eruptive) carbonate assimilation can result in heterogeneous and locally very low δ​11B melt values, similar to predictions for subducted materials. Distinguishing between these processes may be aided by, e.g., the presence or absence of crustal xenolith suites and by thermobarometry to constrain crystallisation depth of the main mineral phases. Boron isotope fractionation in magma via crustal carbonate dissolution is likely to be most pertinent for volcanoes sited on volatile-bearing sedimentary crust, as for instance found in continental subduction settings, and may help to identify upper crustal additions to the carbon cycle.

Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

4

www.nature.com/scientificreports/

Figure 3.  Conceptual model. Tetrahedrally coordinated boron is present in carbonate and silicate melt and decarbonation at the onset of assimilation triggers boron isotope fractionation as follows: CaCO3(BOH)4 (carbonate) +​SiO2(BOH)4 (silicate melt) →​CaO-rich silicate melt +​10B(OH)4 (in melt) +​  CO2 (fluid) +​ 11 B(OH)3 (in CO2-rich fluid). In other words, assimilation of carbonate into the melt gives rise to Ca-rich melt and a co-existing CO2 phase that mingles with CaO-normal melt. Transport of trigonally coordinated 11B in CO2 bubbles away from the reaction site and subsequent partial reabsorption in CO2-undersaturated melt at the distal parts of the capsule gives rise to relatively high δ​11B values in portions of the CaO-normal glass.

Methods

Experimental methods.  The experiments presented in this paper were designed to simulate assimilation of

carbonate crust by magma using Merapi and Vesuvius volcanoes as type examples. All experiments were carried out using the end-loaded piston cylinder apparatus at the HP-HT Laboratory of Experimental Volcanology and Geophysics, at the Istituto Nazionale di Geofisica e Vulcanologia (INGV), Rome, Italy (see http://www.roma1. ingv.it/laboratori/laboratorio-hp-ht/). The end-loaded piston cylinder is calibrated for use in the pressure range 0.5–2 GPa (±​50 MPa) and the lowest end of this range was selected for the experiments (0.5 GPa), which corresponds to approximately mid-crustal depths. The experiments were carried out under super-liquidus conditions for the magmatic starting material by employing an experimental temperature of 1200 °C (±​5  °C). The following starting materials were used in the experiments (see Supplementary Table S1):

(i). Merapi series. A sample of 1994 Merapi basaltic-andesite whole rock powder was fused and hydrated with ultrapure Milli-Q water to ca. 2.5 wt. % H2O. This glass sample was powdered again for insertion into the experimental capsules. The crustal carbonate reactant used in the Merapi series was a sample of local Java platform carbonate, sourced from a quarry to the south of Merapi volcano. The carbonate was inserted into the capsules as small, solid fragments of rock, weighing ~9 to 10 mg. Petrographic descriptions and major element compositions of the starting materials used in the Merapi series are provided in ref. 24. (ii). Vesuvius series. The starting materials used in the Vesuvius experiments were a shoshonitic lava flow from Vesuvius29 fused and hydrated with ultrapure Milli-Q water to ca. 2.0 wt. % H2O and limestone and dolostone from the local Procida carbonate formation36,37. The limestone and dolostone added to the experimental charges were small solid fragments of rock, weighing ~6 to 8 mg. Two sets of experiments were performed containing a limestone and shoshonite hydrated glass and a dolostone and shoshonitic hydrated glass. Detailed descriptions and compositions of the experimental products are provided in ref. 25. The experiment starting materials were loaded into platinum capsules with 3 mm outer diameter before insertion into a 19.1mm NaCl–crushable alumina–pyrophyllite–pyrex assembly38. The capsules were pressurised to 0.5 GPa and then heated from ambient temperature to the target temperature of 1200 °C using a ca. 6 min duration heat-up ramp. Thereafter, the experiments were either i) quenched immediately on reaching the target temperature; these experiments constituted “zero-time” runs; or ii) held at the target temperature for durations of 60 s, 90 s, 150 s (Merapi only), and 300 s before quenching. A detailed description of the experimental method is provided in refs 24 and 25 and references therein. At the end of the runs, the capsules were retrieved from the piston cylinder and mounted in low-volatility Struers EpoFix epoxy resin under vacuum to impregnate the void spaces and prevent sample loss that may otherwise occur due to the high vesicularity of the samples. The long axis of the capsules were placed parallel to the surface of the epoxy block and the capsules were subsequently polished along their long axes until the enclosed experimental products were revealed on the surface at a depth corresponding to approximately the middle of the capsule thickness. The samples were then carbon coated for electron microprobe (EMP) imaging and analysis at Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

5

www.nature.com/scientificreports/ INGV Rome. Representative analyses of the CaO-normal and (MgO)CaO-rich experimental glasses are given in Supplementary Table S2. Prior to SIMS analysis, the carbon coating was gently removed from the experiments by polishing with 1 μ​m diamond-paste, after which the samples were cleaned with pure ethanol and coated with gold for Secondary Ionisation Mass Spectrometry (SIMS) analysis.

SIMS protocol.  Boron isotope analyses were performed using the CAMECA IMS 1280 ion microprobe at the Nordsim facility in Stockholm, Sweden, employing an analytical protocol based on ref. 39. An O2− primary beam of 6 nA with an accelerating voltage of 13.0 kV and imaging a 200 μ​m aperture was used to produce a 20 μ​ m analysis spot. The primary beam was employed to sputter a 25 μ​m square raster for 150 seconds prior to data acquisition to eliminate surface contamination. Secondary ions were transferred to the mass spectrometer using a nominal potential of 10 kV and measured at a mass resolution (M/Δ​M) of 2860, which is sufficient to eliminate isobaric interference of 10B1H+ on 11B. Signals were measured by magnet peak switching into an ion counting electron multiplier, and the sequence comprised a measurement at mass 9.33 (2 s integration time), 10B+ (8 s), 11 + B (4 s), and 30Si++ (2 s). Instrumental mass fractionation was determined and corrected for by employing a reference volcanic glass from Lipari Island (B6) with δ​11BNBS and B concentration values of −​3.3‰ and 197 μ​g/g, respectively34. We note that in an inter-laboratory comparison study, the δ​11B values of B6 showed a large degree of variance of ±​3.6‰ (2 SD)34. Despite this, B6 behaves extremely well on the micron-scale, as it is homogeneous and boron-rich. The value chosen here of −​3.3‰ for B6 is based on analyses by G.D. Layne (pers. comm.) calibrated to reference material GB-4 that has an accepted value of −​12.8‰ (e.g., ref. 40). This value of −​3.3‰ also corresponds to one of the values obtained by positive thermal ionisation mass spectrometry (PTIMS) in the inter-laboratory comparison study34. As our experimental and starting material glasses were both analysed using the same method and employing the same value for standard B6 (−​3.3‰), our data are internally consistent and the relative enrichment and/or depletion patterns in our boron isotope data are valid. However, there may be a systematic bias in our data when compared to data generated by different SIMS laboratories using different values for the B6 glass standard. The 11B yields ranged from 500 to 700 cps/μg​ g−1/nA for the standard. Internal precision (1σ​mean) on 11B/10B ranged from ±​0.5‰ to ±​2.9‰ (semean) for both B6 and the unknowns. The data were acquired over several analytical sessions, with the external error (reproducibility) on B6 ranging from ±​0.42 to 1.29‰ (RSD; n =​  69). External errors were propagated onto the overall analytical uncertainty for each analysis. The uncertainty on the boron concentration measurements (determined from 10B/30Si++ ratios) is estimated at 4.3% (RSD) which, when propagated with the internal errors, gives an absolute uncertainty of ca. 10% (RSD) at the 2σ​level. Boron data for the fused starting materials (n =​ 25) and the experimental products (n =​ 147) are provided in Supplementary Tables S3 and S4.

References

1. Burton, M. R., Sawyer, G. M. & Granieri, D. Deep carbon emissions from volcanoes. Rev. Mineral. Geochem. 75, 323–354 (2013). 2. Dasgupta, R. Ingassing, storage, and outgassing of terrestrial carbon through geologic time. Rev. Mineral. Geochem. 75, 183–229 (2013). 3. Lee, C.-T. A., Thurner, S., Paterson, S. & Cao, W. The rise and fall of continental arcs: interplays between magmatism, uplift, weathering, and climate. Earth Planet. Sci. Lett. 425, 105–119 (2015). 4. Zhang, M. et al. Magma-derived CO2 emissions in the Tengchong volcanic field, SE Tibet: implications for deep carbon cycle at intra-continent subduction zone. J. Asian Earth Sci. doi: org/10.1016/j.jseaes.2016.06.009 (2016). 5. Kerrick, D. & Connolly, J. Metamorphic devolatilization of subducted marine sediments and the transport of volatiles into the Earth’s mantle. Nature 36, 293–296 (2001). 6. Ague, J. J. & Nicolescu, S. Carbon dioxide released from subduction zones by fluid-mediated reactions. Nat. Geosci. 7, 2–7 (2014). 7. Keleman, P. B. & Manning, C. E. Reevaluating carbon fluxes in subduction zones, what goes down, mostly comes up. PNAS 112, E3997–E4006 (2015). 8. Johnston, F. K. B., Turchyn, A. V. & Edmonds, M. Decarbonation efficiency in subduction zones: implications for warm Cretaceous climates. Earth Planet. Sci. Lett. 303, 143–152 (2011). 9. Lee, C.-T. A. et al. Continental arc–island arc fluctuations, growth of crustal carbonates, and long-term climate change. Geosphere 9, 21–36 (2013). 10. Lee, C.-T. A. & Lackey, J. S. Global continental arc flare-ups and their relation to long-term greenhouse conditions. Elements 11, 125–130 (2015). 11. Ryan, J. G. & Chauvel, C. The subduction-zone filter and the impact of recycled materials on the evolution of the mantle in Treatise On Geochemistry 2nd edn (ed. Carlson, R.), Ch. 3.13, 479–508 (Elsevier, 2014). 12. Peacock, S. M. & Hervig, R. L. Boron isotopic composition of subduction-zone metamorphic rocks. Chem. Geol. 160, 281–290 (1999). 13. Hervig, R. et al. Isotopic and elemental partitioning of boron between hydrous fluid and silicate melt. Am. Mineral. 87, 769–774 (2002). 14. Wunder, B., Meixner, A., Romer, R. L., Wirth, R. & Heinrich, W. The geochemical cycle of boron: constraints from boron isotope partitioning experiments between mica and fluid. Lithos 84, 206–216 (2005). 15. Marschall, H. R., Altherr, R. & Rüpke, L. Squeezing out the slab—modelling the release of Li, Be and B during progressive highpressure metamorphism. Chem. Geol. 239, 323–335 (2007). 16. Pabst, S. et al. The fate of subducted oceanic slabs in the shallow mantle: insights from boron isotopes and light element composition of metasomatized blueschists from the Mariana forearc. Lithos 132–133, 162–179 (2012). 17. Leeman, W. P., Tonarini, S., Chan, L. H. & Borg, L. E. Boron and lithium isotopic variations in a hot subduction zone-the southern Washington Cascades. Chem. Geol. 212, 101–124 (2004). 18. Rosner, M., Erzinger, J., Franz, G. & Trumbull, R. B. Slab-derived boron isotope signatures in arc volcanic rocks from the central Andes and evidence for boron isotope fractionation during progressive slab dehydration. Geochem. Geophys. Geosyst. 4, 9005 (2003). 19. Tonarini, S., Armienti, P., D’Orazio, M. & Innocenti, F. Subduction-like fluids in the genesis of My. Etna magmas: evidence from boron isotopes and fluid mobile elements. Earth Planet. Sci. Lett. 192, 471–483 (2001). 20. Tonarini, S. et al. B/Nb amd δ​11B systematics in the Phlegrean volcanic district, Italy. J. Volcanol. Geoth. Res. 133, 123–139 (2004). 21. Hemming, N. G. & Hanson, G. N. Boron isotopic composition and concentration in modern marine carbonates. Geochim. Cosmochim. Acta 56, 537–543 (1992).

Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

6

www.nature.com/scientificreports/ 22. Troll, V. R. et al. Crustal CO2 liberation during the 2006 eruption and earthquake events at Merapi volcano, Indonesia. Geophys. Res. Lett. 39, L11302 (2012). 23. Jolis, E. M. et al. Skarn xenolith record crustal CO2 liberation during Pompeii and Pollena eruptions, Vesuvius volcanic system, central Italy. Chem. Geol. 415, 17–36 (2015). 24. Deegan, F. M. et al. Magma-carbonate interaction processes and associated CO2 release at Merapi volcano, Indonesia: insights from experimental petrology. J. Petrol. 51, 1027–1051 (2010). 25. Jolis, E. M. et al. Experimental simulation of magma–carbonate interaction beneath Mt. Vesuvius, Italy. Contrib. Mineral. Petr. 166, 1335–1353 (2013). 26. Dallai, L., Cioni, R., Boschi, C. & D’Oriano, C. Carbonate-derived CO2 purging magma at depth: influence on the eruptive activity of Somma-Vesuvius, Italy. Earth Planet. Sci. Lett. 310, 84–95 (2011). 27. Troll, V. R. et al. Magmatic differentiation processes at Merapi volcano: inclusion petrology and oxygen isotopes. J. Volcanol. Geoth. Res. 261, 38–49 (2013). 28. Blythe, L. S. et al. CO2 bubble generation and migration during magma-carbonate interaction. Contrib. Mineral. Petr. 169, 1–16 (2015). 29. Di Renzo, V. et al. Magmatic history of Somma-Vesuvius on the basis of new geochemical and isotopic data from a deep borehole (Camaldoli della Torre). J. Petrol. 48, 753–784 (2007). 30. Schiavi, F., Kobayashi, K., Nakamura, E., Tiepolo, M. & Vannucci, R. Trace element and Pb–B–Li isotope systematics of olivinehosted melt inclusions: insights into source metasomatism beneath Stromboli (southern Italy). Contrib. Mineral. Petr. 163, 1011–1031 (2011). 31. Gertisser, R. & Keller, J. Trace element and Sr, Nd, Pb and O isotope variations in medium-K and high-K volcanic rocks from Merapi volcano, Central Java, Indonesia: evidence for the involvement of subducted sediments in Sunda arc magma genesis. J. Petrol. 44, 457–489 (2003). 32. Preece, K. et al. Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia. Contrib. Mineral. Petr. 168, 1061 (2014). 33. Vengosh, A., Kolodny, Y., Starinsky, A., Chivas, A. R. & McCulloch, M. T. Coprecipitation and isotopic fractionation of boron in modern biogenic carbonates. Geochim. Cosmochim. Acta 55, 2901–2910 (1991). 34. Gonfiantini, R. et al. Intercomparison of boron isotope and concentration measurements. Part II: evaluation of results. Geostand. Geoanal. Res. 27, 41–57 (2003). 35. Sanchez-Valle, C. et al. Boron isotopic fractionation between minerals and fluids: new insights from in situ high pressure-high temperature vibrational spectroscopic data. Geochim. Cosmochim. Acta 69, 4301–4313 (2005). 36. Barberi, F., Bizouard, H., Clocchiatti, R., Metrich, N., Santacroce, R. & Sbrana, A. The Somma-Vesuvius magma chamber: A petrological and vulcanological approach. Bull. Volcanol. 44, 295–315 (1981). 37. Civetta, L., Galati, R. & Santacroce, R. Magma mixing and convective compositional layering within the Vesuvius magma chamber. Bull. Volcanol. 53, 287–300 (1991). 38. Freda, C., Baker, D. & Ottolini, L. Reduction of water loss from gold-palladium capsules during piston cylinder experiments by use of pyrophyllite powder. Am. Mineral. 86, 234–237 (2001). 39. Straub, S. M. & Layne, G. D. The systematics of boron isotopes in Izu arc front volcanic rocks. Earth Planet. Sci. Lett. 198, 25–39 (2002). 40. Chaussidon, M., Robert, F., Mangin, D., Hanon, P. & Rose, E. F. Analytical procedures for the measurement of boron isotope compositions by ion microprobe in meteorites and mantle rocks. Geostandard. Newslett. 21, 7–17 (1997). 41. Chaussidon, M. & Marty, B. Primitive boron isotope composition of the mantle. Science 269, 383–386 (1995). 42. Marschall, H. R. & Monteleone, B. D. Boron isotope analysis of silicate glass with very low boron concentrations by secondary ion mass spectrometry. Geostand. Geoanal. Res. 39, 31–46 (2014).

Acknowledgements

We are grateful for technical support from V. Misiti and A. Cavallo at INGV Rome and from L. Ilyinsky and K. Lindén at NRM Stockholm. We also thank L. Blythe for help with the Vesuvius experiments, C. Harris for reviewing an earlier version of the manuscript, and M. Burton, D.R. Hilton, G.D. Layne, and H. Marschall for insightful discussions. This work was supported by the Swedish Research Council (F.M.D., V.R.T.), the Royal Swedish Academy of Science (F.M.D.), Project FIRBMIUR “Development of innovative technologies for environmental protection from natural events” (C.F.), and the Centre for Natural Disaster Science (CNDS) at UU (V.R.T.). The NordSIM facility is supported by the research funding agencies of Denmark, Iceland, Norway and Sweden, the Geological Survey of Finland, and the Swedish Museum of Natural History. This is NordSIM publication #463.

Author Contributions

This study was initiated by F.M.D., V.R.T. and M.J.W. The experimental protocol was designed by V.R.T. and C.F. and the experimental runs were carried out by F.M.D. and C.F. Boron isotope analyses by SIMS were performed by F.M.D., M.J.W. and E.M.J. Data compilation and presentation was led by F.M.D. The manuscript was jointly written by F.M.D. and V.R.T. with input from all authors.

Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep Competing financial interests: The authors declare no competing financial interests. How to cite this article: Deegan, F. M. et al. Boron isotope fractionation in magma via crustal carbonate dissolution. Sci. Rep. 6, 30774; doi: 10.1038/srep30774 (2016). This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ © The Author(s) 2016 Scientific Reports | 6:30774 | DOI: 10.1038/srep30774

7

Boron isotope fractionation in magma via crustal carbonate dissolution.

Carbon dioxide released by arc volcanoes is widely considered to originate from the mantle and from subducted sediments. Fluids released from upper ar...
1MB Sizes 0 Downloads 7 Views