Mantle wedge temperatures and their potential relation to volcanic arc location

The mechanisms underpinning the formation of a focused volcanic arc above subduction zones are debated. Suggestions include controls by: (i) where the subducting plate releases water, lowering the solidus in the overlying mantle wedge; (ii) the location where the mantle wedge melts to the highest degree; and (iii) a limit on melt formation and migration imposed by the cool shallow corner of the wedge. Here, we evaluate these three proposed mechanisms using a set of kinematically-driven 2D thermo-mechanical mantle-wedge models in which subduction velocity, slab dip and age, overriding-plate thickness and the depth of decoupling between the two plates are systematically varied. All mechanisms predict, on the basis of model geometry, that the arc-trench distance, D, decreases strongly with increasing dip, consistent with the negative D-dip correlations found in global subduction data. Model trends of sub-arc slab depth, H, with dip are positive if H is wedge-temperature controlled and overriding-plate thickness does not exceed the decoupling depth by more than 50 km, and negative if H is slab-temperature controlled. Observed global H-dip trends are overall positive. With increasing overriding plate thickness, the position of maximum melting shifts to smaller H and D, while the position of the trenchward limit of the melt zone, controlled by the wedge's cold corner, shifts to larger H and D, similar to the trend in the data for oceanic subduction zones. Thus, the limit imposed by the wedge corner on melting and melt migration seems to exert the first-order control on arc position

A multi-proxy investigation of mantle oxygen fugacity along the Reykjanes Ridge

Mantle oxygen fugacity (fO2) governs the physico-chemical evolution of the Earth, however current estimates from commonly used basalt redox proxies are often in disagreement. In this study we compare three different potential basalt fO2 proxies: Fe3+/Fetot, V/Sc and V isotopes, determined on the same submarine lavas from a 700 km section of the Reykjanes Ridge, near Iceland. These samples provide a valuable test of the sensitivities of fO2 proxies to basalt petrogenesis, as they formed at different melting conditions and from a mantle that towards Iceland exhibits increasing long-term enrichment of incompatible elements. New trace element data were determined for 63 basalts with known Fe3+/Fetot. A subset of 19 lavas, covering the geographical spread of the ridge transect, was selected for vanadium isotope analyses. Vanadium is a multi-valence element whose isotopic fractionation is theoretically susceptible to redox conditions. Yet, the VAA composition of basaltic glasses along the Reykjanes Ridge covers only a narrow range (VAA = −1.09 to −0.86‰; 1SD = 0.02–0.09) and does not co-vary with fractionation-corrected Fe3+/Fetot (0.134–0.151; 1SD = 0.005) or V/Sc (6.6–8.5; 1SD = 0.1-1.3) ratios. However, on a global scale, basaltic VAA may be controlled by the extent of melting. The V/Sc compositions of primitive (MgO > 7.5 wt%) basalts show no systematic change along the entire length of the Reykjanes Ridge. Typical peridotite melting models in which source Fe3+/Fetot is constant at 5% and that account for the increased mantle potential temperature nearer the plume center and the fO2 dependent partitioning of V, can reproduce the V/Sc data. However, while these melting models predict that basalt Fe3+/Fetot ratios should decrease with increasing mantle potential temperature towards Iceland, fractionation-corrected Fe3+/Fetot of Reykjanes Ridge lavas remain nearly constant over the ridge length. This discrepancy is explained by source heterogeneity, where an oxidized mantle pyroxenite component contributes to melting with increasing proximity to Iceland. Comparison of observed and modeled Fe3+/Fetot indicate that source variation in fO2 is present under the Reykjanes Ridge, with higher Fe3+/Fetot closer to Iceland. This source variability in fO2 cannot be resolved by V isotopes and redox-sensitive trace element ratios, which instead appear to record magmatic processes

Melt generation beneath Arctic Ridges: Implications from Ule

We present new 238U-230Th-226Ra-210Pb, 235U-231Pa, and Nd, Sr, Hf, and Pb isotope data for the slow- to ultraslow-spreading Mohns, Knipovich, and Gakkel Ridges. Combined with previous work, our data from the Arctic Ridges cover the full range of axial depths from the deep northernmost Gakkel Ridge shallowing upwards to the Knipovich, Mohns, and Kolbeinsey Ridges north of Iceland. Age-constrained samples from the Mohns and Knipovich Ridges have (230Th/238U) activity ratios ranging from 1.165 to 1.30 and 1.101 to 1.225, respectively. The high 230Th excesses of Kolbeinsey, Mohns, and Knipovich mid-ocean ridge basalts (MORB) are erupted from ridges producing relatively thin (Mohns, Knipovich) to thick (Kolbeinsey) oceanic crust with evidence for sources ranging from mostly peridotite (Kolbeinsey) to eclogite-rich mantle (Mohns, Knipovich). Age-constrained lavas from 85ºE on the Gakkel Ridge, on the other hand, overlie little to no crust and range from small (~5%) 230Th excesses to small 238U excesses (~5%). The strong negative correlation between (230Th/238U) values vs. axial ridge depth among Arctic ridge basalts is controlled not only by solidus depth influence on 238U-230Th disequilibria, but also by variations in mantle source lithology and depth to the base of the lithosphere, which is expected to vary at ultra-slow spreading ridges. Small 231Pa excesses (65% excess) in age constrained basalts support the presence of eclogite in the mantle source for this region. Conversely, the ultraslow-spreading Gakkel Ridge basalts are homogeneous, with Sr, Nd, and Hf radiogenic isotopic signatures indicative of a long time-averaged depleted mantle source. The Gakkel samples have minimum (226 Ra/230Th) ratios ranging from 3.07 to 3.65 ± 3%, which lie along and extend the global negative correlation between 226Ra and 230Th excesses observed in MORB. The new 230Th-226Ra data support a model for global MORB production in which deep melts record interaction with shallower materials. This scenario requires either mixing with shallow-derived melts, or melt-rock reaction with shallower rocks in the lithosphere or crust

Equilibrium partitioning and isotopic fractionation of nitrogen between biotite, plagioclase, and K-feldspar during magmatic differentiation

Funding: UK Natural Environment Research Council (NE/R012253/1, NE/V010824/1, NE/P012167/1)A significant portion of the continental crust is composed of plutonic igneous rocks. However, little is known about the geochemical behaviour of N between the different minerals during magmatic differentiation. To provide new constraints for the behaviour of N during crust formation, we have characterised the geochemistry of nitrogen (N) in the compositionally zoned calc-alkaline pluton at Loch Doon, SW Scotland. We present N concentration and N isotope values for whole-rock data alongside biotite, plagioclase and K-feldspar mineral separates and assess the degree to which these data preserve equilibrium partitioning during magmatic differentiation. We show that whole rock likely inherited its N contents and δ15N signatures from the initial source composition and that this signature is homogenous at a pluton scale. Whilst the whole-rock data are best explained as crust-derived N in the source, the degree of homogenisation across a pluton scale is inconsistent with empirical N diffusivities, ruling out syn-emplacement crustal assimilation as the source of N. Instead, our data suggest a crustal signature inherited from depth associated with the Iapetus subduction zone. At a mineral scale, we find that N preferentially partitions into the feldspars over mica in this system in the order K-feldspar > plagioclase ≈ biotite > quartz, with average mineral-mineral distribution coefficients of DNplagioclase-biotite = 1.3 ± 0.6 and DNKspar-biotite = 2.8 ± 0.6. Partitioning is accompanied by a large and near constant equilibrium isotope fractionation factor between biotite and both feldspars (averages are Δ15NPlag-Biotite = +7.8 ± 1.2 ‰ and Δ15NKspar-Biotite = +7.9 ± 1.0 ‰), whereas Δ15NKspar-Plagioclase closely approximates 0 ‰, where both minerals show δ15N overlapping with the bulk rock δ15N values. These results show that mica crystallisation generates in a large negative Δ15N resulting a 15N-depleted reservoir within plutonic rocks. Moreover, our dataset suggests that feldspars might be a more significant host of N in the igneous portion of Earth’s continental and oceanic crust than previous thought.Publisher PDFPeer reviewe

Rhenium elemental and isotopic variations at magmatic temperatures

This work was funded by Natural Environment Research Council UK Standard Grant to RGH, AJD, and JP (NE/T001119).Recent analytical advances in the measurement of rhenium (Re) isotope ratios allow its potential as a palaeoredox and chemical weathering proxy to be explored. However, a successful isotopic proxy must be grounded by an understanding of its composition and behaviour in the solid Earth. Here, we present Re concentrations and Re isotopic (δ187Re) compositions for a well-characterised sequence of lavas from Hekla volcano, Iceland. The concentration of Re varies from 0.02 to 1.4 ng/g, decreasing from basalt to more evolved lavas. We show that the crystallisation and removal of magnetite is responsible for the Re decrease in this system. By contrast, δ187Re values for the same suite of samples show a relatively narrow range (−0.45 to −0.22 0/1000), suggesting minimal resolvable Re isotope fractionation between magnetite and the silicate melt. Together with other samples, including mid-ocean ridge basalts, these first igneous data can be used to estimate a baseline for terrestrial materials (δ187Re = −0.33 ± 0.15 0/1000, 2 s.d., n = 14), from which low-temperature Re isotope variations in Earth’s surficial environments can be assessed, alongside the global isotope mass balance of Re.Peer reviewe

Nucleosynthetic Heterogeneity Controls Vanadium Isotope Variations in Bulk Chondrites

The vanadium (V) isotope composition of early solar system materials have been hypothesized to be sensitive to high energy irradiation that originated from the young Sun. Vanadium has two isotopes with masses 50 and 51 that have (51)V/(50)V ratio of approximately 410. High energy irradiation produces (50)V from various target isotopes of Ti, Cr and Fe, which would result in light V isotope compositions (expressed as delta (51)V in per mille = 1000 x (((51)V/(50)V(sub sample)/(51)V/(50)V(sub AlfaAesar)) - 1)) relative to a presumably chondritic starting composition. Recently published V isotope data for calcium aluminium inclusions (CAIs) has revealed some very negative values relative to chondrites (by almost -4 per mille) that were indeed interpreted to reflect irradiation processes despite the fact that the studied CAIs all exhibited significant initial abundances of (10)Be, while only a few CAIs displayed light V isotope compositions. It is difficult to relate V isotope variations directly to a singular process because V only possesses two isotopes. Therefore, V isotope variations can principally be produced both mass dependent and independent processes. Mass dependent kinetic stable isotope fractionation is common in CAIs for refractory elements due to partial condensation/evaporation processes. The element strontium (Sr) has an almost identical condensation temperature to V and studies of stable Sr isotope compositions in CAIs reveal both heavy and light values relative to chondrites of several permil. These variations are similar in magnitude to those reported for V isotopes in CAIs, which suggests it is possible that some of the V isotope variation in CAIs could be due to kinetic stable isotope fractionation during condensation/evaporation processes

No V-Fe-Zn isotopic variation in basalts from the 2021 Fagradalsfjall eruption

The Earth’s mantle is chemically heterogeneous in space and time, which is often reflected by variable isotopic compositions of mantle derived basalts. Basalts from the first 40 days of the 2021 Fagradalsfjall eruption, Reykjanes Peninsula, Iceland, display systematic temporal variations in the ratios of incompatible elements alongside resolvable variations in Sr, Nd and Pb radiogenic isotopes. These variations reflect progressive influx of magma derived from melting of a deeper, more enriched and potentially lithologically distinct source. We use this eruptive time series to conduct the first combined V-Fe-Zn isotope study, exploring the sensitivity of the combined isotopic approach, with particular focus on fingerprinting source lithological heterogeneity. We find no analytically resolvable change in V (δ51VAA between −0.95 ± 0.09 ‰ 2 s.d. and −0.86 ± 0.07 ‰ 2 s.d.), Fe (δ56FeIRMM-524 between +0.047 ± 0.042 ‰ 2 s.d. and +0.094 ± 0.049 ‰ 2 s.d.) and Zn (δ66ZnAA-ETH between −0.042 ± 0.003 ‰ 2 s.d. and +0.013 ± 0.027 ‰ 2 s.d.) isotopic compositions. The lack of variability in V-Fe-Zn isotopes, despite the evolving trace element and radiogenic isotope ratios, suggests there is no significant contribution of melts from a lithologically distinct (pyroxenite) mantle component under the Reykjanes Peninsula

Thallium concentration and thallium isotope composition of lateritic terrains

Continental weathering plays a key role in modifying the geochemical budget of terrestrial reservoirs. Laterites are the products of extreme sub-aerial continental weathering. This study presents the first investigation of thallium (Tl) abundances and stable isotopic compositions of lateritic terrains. Two laterite profiles from India of differing protolith and age are studied. Thallium concentrations range between 7 and 244 ng/g for a basalt-based lateritic profile and 37–652 ng/g within a greywacke lateritic profile. The average Tl stable isotope composition of the two profiles is similar to many typical igneous materials, however, the intense tropical weathering causes a small but resolvable fractionation of Tl stable isotopes towards heavy values in the residual soils. The profiles are dominated by significant positive isotope excursions (reported as ε205Tl relative to standard NBS997) of +3.5 ± 0.5 2sd and +6.2 ± 0.5 2sd at the inferred palaeowater tables within both laterite profiles. These signatures likely reflect combined changes in redox state and mineralogy. Extensive mineral dissolution under through-flowing fluids alongside the formation of new phases such as phyllosilicates and Mn- and Fe-oxides and hydroxides likely account for some of the Tl mobilisation, sorption and coprecipitation. In the case of laterites, the formation of the new phases and role of surface sorption likely contribute to stable Tl isotope fractionation. The identification of strong isotope excursions at inferred palaeowater tables encourages future research to determine specific mineral phases that may drive significant fractionation of Tl stable isotopes. This study showcases the magnitude of natural variation possible in terrestrial soils. Such information is key to the nascent application of Tl isotope compositions as tracers of anthropogenic pollution

A new method for the determination of low-level actinium-227 in geological samples

Author Posting. © The Author(s), 2012. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Journal of Radioanalytical and Nuclear Chemistry 296 (2013): 279-283, doi:10.1007/s10967-012-2140-0.We developed a new method for the determination of 227Ac in geological samples. The method uses extraction chromatographic techniques and alpha-spectrometry and is applicable for a range of natural matrices. Here we report on the procedure and results of the analysis of water (fresh and seawater) and rock samples. Water samples were acidified and rock samples underwent total dissolution via acid leaching. A DGA (N,N,N’,N’-tetra-n-octyldiglycolamide) extraction chromatographic column was used for the separation of actinium. The actinium fraction was prepared for alpha spectrometric measurement via cerium fluoride micro-precipitation. Recoveries of actinium in water samples were 80±8 % (number of analyses n=14) and in rock samples 70±12 % (n=30). The minimum detectable activities (MDA) were 0.017-0.5 Bq kg-1 for both matrices. Rock sample 227Ac activities ranged from 0.17 to 8.3 Bq kg-1 and water sample activities ranged from below MDA values to 14 Bq kg-1of 227Ac. From the analysis of several standard rock and water samples with the method we found very good agreement between our results and certified values

Cumulate causes for the low contents of sulfide-loving elements in the continental crust

Despite the economic importance of chalcophile (sulfide-loving) and siderophile (metal-loving) elements (CSEs), it is unclear how they become enriched or depleted in the continental crust, compared with the oceanic crust. This is due in part to our limited understanding of the partitioning behaviour of the CSEs. Here I compile compositional data for mid-ocean ridge basalts and subduction-related volcanic rocks. I show that the mantle-derived melts that contribute to oceanic and continental crust formation rarely avoid sulfide saturation during cooling in the crust and, on average, subduction-zone magmas fractionate sulfide at the base of the continental crust prior to ascent. Differentiation of mantle-derived melts enriches lower crustal sulfide- and silicate-bearing cumulates in some CSEs compared with the upper crust. This storage predisposes the cumulate-hosted compatible CSEs (such as Cu and Au) to be recycled back into the mantle during subduction and delamination, resulting in their low contents in the bulk continental crust and potentially contributing to the scarcity of ore deposits in the upper continental crust. By contrast, differentiation causes the upper oceanic and continental crust to become enriched in incompatible CSEs (such as W) compared with the lower oceanic and continental crust. Consequently, incompatible CSEs are predisposed to become enriched in subduction-zone magmas that contribute to continental crust formation and are less susceptible to removal from the continental crust via delamination compared with the compatible CSEs