Search

Meaningful analysis of uranium‐series isotopic disequilibria in basaltic lavas relies on the use of complex forward numerical models like dynamic melting (McKenzie, 1985, https://doi.org/10.1016/0012-821x(85)90001-9) and equilibrium porous flow (Spiegelman & Elliott, 1993, https://doi.org/10.1016/0012-821x(93)90155-3). Historically, such models have either been solved analytically for simplified scenarios, such as constant melting rate or constant solid/melt trace element partitioning throughout the melting process, or have relied on incremental or numerical calculators with limited power to solve problems and/or restricted availability. The most public numerical solution to reactive porous flow, UserCalc (Spiegelman, 2000, https://doi.org/10.1029/1999gc000030) was maintained on a private institutional server for nearly two decades, but that approach has been unsustainable in light of modern security concerns. Here, we present a more long‐lasting solution to the problems of availability, model sophistication and flexibility, and long‐term access in the form of a cloud‐hosted, publicly available Jupyter notebook. Similar to UserCalc, the new notebook calculates U‐series disequilibria during time‐dependent, equilibrium partial melting in a one‐dimensional porous flow regime where mass is conserved. In addition, we also provide a new disequilibrium transport model which has the same melt transport model as UserCalc, but approximates rate‐limited diffusive exchange of nuclides between solid and melt using linear kinetics. The degree of disequilibrium during transport is controlled by a Damköhler number, allowing the full spectrum of equilibration models from complete fractional melting (Da=0) to equilibrium transport (Da=∞). Access the full executable Jupyter notebook version of this manuscript here. Key Points: Cloud‐based Jupyter notebook presents an open source, reproducible tool for modeling U‐series in basaltsEquilibrium and pure disequilibrium porous flow U‐series models with 1D conservation of massScaled porous flow model introduces incomplete equilibrium scenario with reaction rate limitations [ABSTRACT FROM AUTHOR]

We address the question of whether melting of the mafic oceanic crust occurs beneath ordinary volcanic arcs using constraints from U-Series ( 238 U/ 232 Th, 230 Th/ 232 Th and 226 Ra/ 230 Th) measurements. Alteration of the top few hundred meters of the mafic crust leads to strong U enrichment. Via decay of 238 U to 230 Th, this results in elevated ( 230 Th/ 232 Th) (where brackets indicate activity ratios) over time-scales of ∼350 ka. This process leads to the high ( 230 Th/ 232 Th), between 2.6 and 11.0 in the mafic altered oceanic crust (AOC) sampled at ODP Sites 801 and 1149 near the Izu-Bonin-Mariana arc. Th activity ratios in the Izu arc lavas range from ( 230 Th/ 232 Th) = 1.2–2.0. These values are substantially higher than those in bulk sediment subducting at the Izu trench and also extend to higher values than in mid-ocean ridge basalts and the Mariana arc. We show that the range in Th isotope ratios in the Izu arc lavas is consistent with the presence of a slab melt from a mixed source consisting of AOC and subducted sediments with an AOC mass fraction of up to approximately 80 wt.% in the component added to the arc lava source. The oceanic plate subducting at the Izu arc is comparatively cold which therefore indicates that temperatures high enough for fluid-saturated melting of the AOC are commonly achieved beneath volcanic arcs. The high ratio of AOC/sediments of the slab melt component suggested for the Izu arc lavas requires preferential melting of the AOC. This can be achieved when fluid-saturated melting of the slab is triggered by fluids derived from underlying subducted serpentinites. Dehydration of serpentinites and migration of the fluid into the overlying crust causes melting to start within the AOC. The absence of a significant sediment melt component suggests there was insufficient water to flux both AOC and overlying sediments. [ABSTRACT FROM AUTHOR]

The island of São Miguel, Azores piques geochemical interest due to a strikingly large range in the long-lived radiogenic isotope ratios of its lavas. The 'enriched' signatures (e.g., radiogenic 87Sr/86Sr and unradiogenic 143Nd/144Nd) observed in lavas found in the east of the island have been proposed to originate from recycled crustal components in their mantle sources. Such fertile lithologies should have higher melt productivities (amount of melt generated per decrement in pressure) than peridotitic mantle and this heritage ought to be evident in the U-series signatures of melts derived from such materials. Specifically, combined 238U-230Th and 235U-231Pa disequilibria in erupted lavas should identify significant differences in the melting rate and/or the mineralogical makeup of their sources (e.g., the amount of garnet). To this end, new 238U-230Th-226Ra and 235U-231Pa disequilibria are presented for sixteen mafic lavas (< 50 wt.% SiO2) that encompass the full range of São Miguel's long-lived radiogenic isotope variability. However, primary U-series signatures of mafic lavas can be influenced by assimilation of more differentiated products. Notably, syenite nodules and cumulate alkali feldspar xenocrysts have been identified in São Miguel basalts. We argue that secondary contamination of basalts is most evident in the hugely variable 230Th-226Ra disequilibria (~ 18% deficits to > 300% excesses of 226Ra) that strongly correlate with Ba/Th ratios. The longer-lived U-series nuclides are only slightly perturbed in a number of the more enriched lavas, resulting in increased noise in the dataset. In apparently uncontaminated samples, isotopically depleted lavas from the west of São Miguel display slightly higher 230Th and 231Pa excesses than the enriched samples from the east. The U-series signature of the enriched lavas can be modeled with a combination of higher absolute values of D¯U U = and D¯Th,, coupled with higher melting rates compared to the depleted western source. These parameters are consistent with a peridotite source enriched with contributions from more fusible lithologies such as eclogite and/or pyroxenite in the eastern source, but not with direct derivation from such fusible lithologies. The only subtle variations in U-series isotopes coupled with large, systematic variations in radiogenic isotope compositions suggest that enriched components and/or their melts are well-mixed with mantle peridotite, thus masking their distinctive melting behavior. Finally, basalts from Sete Cidades volcano in the west of São Miguel have very similar long-lived isotope compositions and U-series disequilibria to basalts from the island of Pico, Azores. The geochemical similarity of lavas from islands at the center and periphery of the archipelago argue against systematic differences in melting behavior related to sampling different portions of an underlying plume structure. [ABSTRACT FROM AUTHOR]

Abstract: No young, off-axis, mid-ocean ridge lavas have yet been directly linked to underlying off-axis melt bodies. In this study, we present new measurements of 238U–230Th–226Ra–210Pb isotope compositions for a suite of lavas from the overlapping spreading center (OSC) at 9°03′N on the East Pacific Rise (EPR). These lavas span a large range of compositions, from basalt to dacite, and include both axial and off-axis samples recovered from a prominent, axis-parallel pillow ridge and a flat-topped seamount that overlie the westernmost extent of a 4-km-wide melt lens (Kent et al., 2000). We report 210Pb excesses in axial basalts and basaltic andesites, which we suggest results from gas-magma fractionation of 222Rn from 226Ra beneath dacite magmas. In addition, our U-series ages agree with visual observations, indicating that while most recent volcanic activity occurs at the spreading axis, active volcanism also occurs away from the axis. Specifically, the off-axis pillow ridge and seamount samples overlying the off-axis subsurface melt body have eruption ages of less than 8ka, and likely as young as 1ka, despite being located on crust that has a spreading age of ~75ka. The young ages of these lavas, combined with existing geological, geochemical and geophysical constraints, provide evidence for a genetic link between the pillow ridge and seamount lavas and the seismically imaged, underlying off-axis melt lens. This link demonstrates that off-axis volcanism does not necessarily come from a sub-axial magma body and can be sourced directly from off-axis magma bodies. [Copyright &y& Elsevier]

The Auckland volcanic field is a Quaternary monogenetic basaltic field of 50 volcanoes. Rangitoto is the most recent of these at ~500 year BP and may mark a change in the behaviour of the field as it is the largest by an order of magnitude and is unusual in that it erupted magmas of alkalic then subalkalic basaltic composition in discrete events separated by ≤50 years. Major and trace element geochemistry together with Sr-Nd and U-Th-Ra isotopes provides the basis for modelling the melting conditions that brought about the eruption of two chemically different lavas with very little spatial or temporal change. Sr-Nd isotopes suggest that the source for both eruptions is similar with a slight degree of heterogeneity. The basalts show high Th-excess compared with comparable continental volcanic fields. We show that the alkalic basalts give evidence for lower degrees of partial melting, higher amounts of residual garnet, a longer melting column and lower melting and upwelling rates compared with the subalkalic basalts. The low upwelling rates (0.1-1.5 cm/year) modelled for both magmas do not suggest a plume or major upwelling in the mantle region beneath Auckland; therefore, we suggest localised convection due to relict movement from the active subduction system situated 400 km to the southeast. A higher porosity for the initial alkalic basalt is based on Ra-excesses, suggesting movement of melt by two different porosities: the initial melt travelling in fast high porosity channels from greater depths preserving a high Th-excess and the subsequent subalkalic magma travelling from a shallower depth through lower porosity diffuse channels preserving a high Ra-excess; this creates a negative array in (Ra/Th) versus (Th/U) space previously only seen in mid ocean ridge Basalt data. This mechanism suggests the Auckland volcanic field may operate by the presence of discrete melt batches that are able to move at different depths and speeds giving the field its erratic spatial and temporal pattern of eruptions, a type of behaviour that may have implications for the evolution of other continental volcanic fields worldwide. [ABSTRACT FROM AUTHOR]

The discovery of chemically and isotopically enriched mid-ocean ridge basalts (E-MORB) has offered substantial insight into the origin, time scales, and length scales of mantle heterogeneity. However, the exact processes involved in producing this E-MORB enrichment are vigorously debated. Additionally, because the ages of E-MORB are not well constrained, the petrogenetic, temporal, and geological relationships between E-MORB and normal (N)-MORB are not known. To investigate these relationships and to explore how melting and melt transport processes contribute to or modify enriched mantle source compositions and generate E-MORB melts beneath mid-ocean ridges, we measured major and trace elements, and Sr, Nd, Hf, Pb, and U–Th–Ra isotopes for a suite of lavas that were collected off-axis, including several E-MORB, at 9–10°N along the East Pacific Rise (EPR). These data show coherent mixing trends among long-lived radiogenic isotopes, U-series nuclides, and incompatible trace elements, implying that mixing of melts from different sources occurs at different depths. Our results are consistent with previous studies that show that melting occurs in a two-porosity melting regime, with high-porosity channels forming deeply in the presence of garnet and transporting enriched melts with large 230Th excesses to the crust, whereas low-porosity channels transport melts more slowly, allowing them to equilibrate at shallow depths and develop large 226Ra excesses at the expense of diminished 230Th excesses. Forward modeling of the trace element data also is consistent with mixing of melts in a two-porosity melting regime. U-series age constraints suggest that E-MORB neither erupt at systematically different times from N-MORB, nor necessarily through different pathways. Previous studies of E-MORB at 9–10°N have suggested that E-MORB compositions could be explained by off-axis eruption. However, when considered in light of previously published magnetic paleointensity and U-series age constraints, recent geological studies, and the most widely accepted contemporary understanding of volcanic construction at 9–10°N EPR, the asymmetric, off-axis distribution of E-MORB at 9–10°N EPR is consistent with, and more simply explained by, a model in which E-MORB erupted within the axial summit trough (AST) and flowed down the ridge flanks (∼0–3 km). These E-MORB subsequently spread away from the AST, and, finally, were preserved on the seafloor through asymmetric construction of the extrusive layer. Taken together, the range of ages of E-MORB at 9–10°N EPR and the geochemical and isotopic mixing trends suggest that enriched melts are continuously supplied to the ridge axis, but because of their small proportions relative to the volumetrically and volcanically dominant N-MORB, E-MORB preservation and exposure is comparatively scarce. [ABSTRACT FROM AUTHOR]

Abstract: In historical basalts from a wide range of tectonic settings, 210Pb is often found to have an activity deficit relative to its predecessor 226Ra. Several processes have been hypothesized as causes of 210Pb deficits in basalts. In subduction zone and ocean island environments, 210Pb deficits have often been attributed to shallow magmatic degassing. At mid-ocean ridges, 210Pb deficits have been inferred to result from mantle melting, limiting the time between melt production and eruption to 100 years or less. Here we present an alternative mechanism for producing 210Pb deficits in magmas, by diffusive exchange between a melt and cumulate minerals (plagioclase and/or clinopyroxene) in the crust. The deficit in 210Pb develops in response to its decay toward secular equilibrium with 226Ra within the mineral grains; decay provides an internal sink for 210Pb that drives continuous diffusive exchange with the melt. Deficits in 210Pb develop under a broad range of conditions, in enriched and depleted melts and during interaction with young or old cumulate minerals. The magnitude of the deficit depends mainly on the equilibrium mineral/melt partition coefficients for Pb and Ra and on the melt/rock ratio during diffusive interaction, and is only weakly dependent on the relative diffusivities of Ra and Pb in the minerals or the trace element disequilibrium between the melt and cumulate minerals. Plagioclase in the crust has greater leverage on the 210Pb–226Ra system than any silicate mineral present during mantle melting, and is capable of inducing significant 210Pb deficits in the melt even at melt fractions above 50%. Its influence on the melt is also rapid, with a substantial 210Pb deficit developing in less than a year and approaching a steady state value after several decades or less. The strong control crustal cumulates are capable of exerting on 210Pb–226Ra fractionation in melts indicates that they may have a significant role in a wide range of tectonic environments, and suggests caution in interpreting 210Pb deficits as a signature of mantle melting, or as a product of 222Rn degassing. [Copyright &y& Elsevier]

Rhyolite pumices and co-erupted granophyric (granite) xenoliths yield evidence for rapid magma generation and crystallization prior to their eruption at 15·2 ± 2·9 ka at the Alid volcanic center in the Danikil Depression, Eritrea. Whole-rock U and Th isotopic analyses show 230Th excesses up to 50% in basalts <10 000 years old from the surrounding Oss lava fields. The 15 ka rhyolites also have 30–40% 230Th excesses. Similarity in U–Th disequilibrium, and in Sr, Nd, and Pb isotopic values, implies that the rhyolites are mostly differentiated from the local basaltic magma. Given the (230Th/232Th) ratio of the young basalts, and presumably the underlying mantle, the (230Th/232Th) ratio of the rhyolites upon eruption could be generated by in situ decay in about 50 000 years. Limited (∼5%) assimilation of old crust would hasten the lowering of (230Th/232Th) and allow the process to take place in as little as 30 000 years. Final crystallization of the Alid granophyre occurred rapidly and at shallow depths at ∼20–25 ka, as confirmed by analyses of mineral separates and ion microprobe data on individual zircons. Evidently, 30 000–50 000 years were required for extraction of basalt from its mantle source region, subsequent crystallization and melt extraction to form silicic magmas, and final crystallization of the shallow intrusion. The granophyre was then ejected during eruption of the comagmatic rhyolites. [ABSTRACT FROM PUBLISHER]

Abstract: Trace element and U-series isotope analyses are presented for a suite of recent (<10 ka) basalts from the axial portion of the Kenya rift. Samples from throughout the rift have LREE-enriched patterns with HREE >10×chondrite and the LREE between 60 and 200×chondrite. REE fractionation is consistent with melting a garnet lherzolite source region with between 2% and 6% modal garnet. Other trace element ratios are distinct from OIB, notably Zr/Hf which ranges from 43 to 48, whilst at a given Zr content the Zr/Hf ratio is significantly greater than that found in OIB. (238U/232Th) range from 0.362 to 1.036, (230Th/232Th) from 0.503 to 1.109, with (230Th/238U) ranging from 0.783 to 2.966. All but two samples are in 230Th excess or in secular equilibrium. Samples with elevated (238U/232Th), also have Rb/Cs>120, but unexceptional 208Pb⁎/206Pb⁎ and hence κ Pb values. These samples have experienced U and Cs loss and are excluded from further consideration. Of the unaltered samples, all have (238U/232Th) generally lower than OIB, with maximum values of <0.8, and some <0.6. The maximum (230Th/238U) is 1.39, similar to OIB. Although none of the basalts has a primary composition, (230Th/238U) does not vary systematically with indices of fractionation, and comparison with evolved rocks from Kenya indicates that 238U–230Th disequilibrium in the basalts is not the product of fractionation and crustal residence, but a product of melt generation. The maximum (230Th/238U) that can be generated by batch melting, assuming a source mineralogy consistent with the REE variation is 1.05 and so the variation in (230Th/238U) is attributed to more complex models of melt generation and/or transport. Both dynamic melting and equilibrium porous flow suggest mantle upwelling rates of ≤2 cm year−1. It is suggested that the Kenya basalts represent melts derived from lithospheric mantle that has been thermally reactivated by and incorporated into the underlying (East African) mantle plume. [Copyright &y& Elsevier]

Abstract: We report U–Th disequilibria data for a suite of 13 young basaltic samples from the Samoan Islands, which represent the end-member mantle component EM2, and 4 historic lavas from Mt. Erebus, typifying young HIMU. The Samoan samples have low 230Th/232Th and 238U/232Th, consistent with the enriched nature (EM2) of the Samoan mantle source, whereas the Mt. Erebus samples have high 206Pb/204Pb and intermediate (230Th/232Th) and (238U/232Th). When considered in the context of the global oceanic basalt database, the Samoan samples'' low 230Th/232Th and 238U/232Th greatly extend the global correlations between 230Th/232Th and 238U/232Th with 87Sr/86Sr and 143Nd/144Nd and change the functional form of these correlations from linear to hyperbolic. Using a maximum likelihood non-linear inversion method, we show that these correlations of 230Th/232Th and 238U/232Th with 87Sr/86Sr and 143Nd/144Nd can be approximated by two-component mixing. However, the global oceanic basalt data also show considerable scatter about the best-fit mixing curves. This scatter is attributed to additional source components and melting processes influencing the lavas (230Th/232Th) and (238U/232Th). Additional source components are supported by the Pb isotope data, which clearly require more than two endmember mantle components. We also show that the extent and variability in (230Th/238U) decreases as a function of source enrichment, with “depleted” MORB showing the largest extents and greatest variability in (230Th/238U) and “enriched” OIB, like Samoa, showing smaller extents and less variability of (230Th/238U). We interpret this observation in terms of differences in the melting regimes beneath mid-ocean ridges and ocean islands. [Copyright &y& Elsevier]

Abstract: U-series studies in oceanic basalts are critical for understanding melting and melt extraction in the Earth''s mantle. The combined results of a detailed geochemical study of melting and melt extraction at Theistareykir, northern Iceland, provide a strong case for melt extraction via channeled melt flow at an active spreading ridge. It has often been argued, however, that widely used melting and melt extraction models, which simulate channeled melt extraction (i.e. fractional and/or dynamic melting), can only partially explain the global U-series systematics in oceanic basalts. As a consequence, more complicated models have been invoked, which suggest different styles of melt extraction at different depths/pressures in the mantle, so-called “two-porosity models”. Alternatively, diffusion-controlled mechanisms have been proposed. Here we show that U–Th–Pa–Ra systematics in oceanic basalts can indeed be explained by models where melt transport occurs without chemical equilibrium between melt and solid when variations in all three critical melting parameters (residual porosity, upwelling rate of the solid mantle and melt velocity) are taken into account. Melting at ridges requires systematic variation of at least two critical melting parameters, most likely upwelling and melt extraction rate. Melts generated with increasing lateral distance to the ridge axis are generated with slower upwelling rates and are extracted with lower velocities than melts created closer to the ridge axis. Melting at ocean islands, on the other hand, can successfully be explained by variations in upwelling rate only. Global U-series systematics in OIB originate from superimposed global variations in upwelling velocity due to different buoyancy fluxes and from local variation in upwelling velocity as a function of radial distance to the plume center. The model proposed here is consistent with other geochemical data for oceanic basalts and strongly supports melt extraction via high-porosity channels as a general means of melt extraction from the Earth''s upper mantle. [Copyright &y& Elsevier]

Abstract: We report U-series systematics (238U-230Th-231Pa-226Ra) for basalts from the Azores islands and the nearby mid-Atlantic ridge with the aim of constraining melting processes in this region of plume–ridge interaction. Both 230Th and 231Pa excess show rough negative trends with indices of source enrichment (Sr, Nd and Pb isotopes) which cannot be explained by melting of a source variably enriched by mafic components (e.g. pyroxenite). However, there is evidence for a distinct source component, possibly sediment, beneath Sao Miguel which may give rise to these trends. 231Pa excess and Sm-Nd fractionation are positively and negatively correlated (respectively) with distance from the centre of the plume, while 230Th excesses show no clear trend. This suggests that mantle upwelling velocity exerts a strong control on U-Pa fractionation and that melting is initiated deeper in the centre of the plume. These observations are used to infer melting conditions beneath the mid-Atlantic ridge and Azores islands. The mantle upwelling velocity is constrained to be 3–4 cm/a which is smaller than that inferred for Hawaii or Iceland. Predicted 226Ra-excesses are similar to, or lower than those measured, suggesting high velocity melt ascent in channels and rapid differentiation prior to eruption. The effect of the presence of volatiles in the source is also investigated and we show that the presence of water could explain the 230Th and 231Pa trends as a function of distance across the plume. A 3-D model for the Azores plume can reproduce the U-series isotope observations. [Copyright &y& Elsevier]

Basaltic lavas provide a window into their mantle sources and the associated geodynamic processes associated with them. Obtaining absolute chronology on young basalt flows is often challenging. Some gypsum crusts in pyro-duct (lava tube) caves crystallize during the last "cooling" stages of tube formation probably at relatively low temperatures (i.e., <100 °C). Here we propose uranium-series dating of these 'speleogenetic' gypsum crusts as a way to establish the ages of young basalt flows. We demonstrate the technique using cave crusts from four young basalt flows in El Malpais National Monument, New Mexico, three of which are moderately well-dated by alternative methods (14C, 36Cl, and 3He ages). They all exhibit well-formed lava tube systems that contain gypsum crust speleothems interpreted to be remnant of the volcanic sulfur or sulfate processes when the tube systems were active during their terminal stages. δ18O SO4 , δ34S, 87Sr/86Sr, and δ234U values for the gypsum crusts have near-basaltic signatures, and the crusts contain uranium concentrations suitable for uranium-series dating in three of these four basalt flows. Analyses of gypsum from caves of two well-dated basalt flows yielded uranium-series ages that agree with the ages determined by other techniques showing that absolute ages of young basalt flows can be achieved by this method. The fourth flow, the Hoya de Cibola basalt, previously undated by other absolute means has a U–Th gypsum crust isochron age of 16.9 ± 0.9 ka. The uranium-series age and isotopic analyses of these gypsum crusts should be useful in geologic applications involving young basalt flows. [ABSTRACT FROM AUTHOR]