Porosity of the melting zone and variations in the solid mantle upwelling rate beneath Hawaii: inferences from 238U-230Th-226Ra and 235U-231Pa disequilibria

Measurements of ²³⁸U-²³⁰Th-²²⁶Ra and ²³⁵U-²³¹Pa disequilibria in a suite of tholeiitic-to-basanitic lavas provide estimates of porosity, solid mantle upwelling rate and melt transport times beneath Hawaii. The observation that (²³⁰Th/²³⁸U) > 1 indicates that garnet is required as a residual phase in the magma sources for all of the lavas. Both chromatographic porous flow and dynamic melting of a garnet peridotite source can adequately explain the combined U-Th-Ra and U-Pa data for these Hawaiian basalts. For chromatographic porous flow, the calculated maximum porosity in the melting zone ranges from 0.3–3% for tholeiites and 0.1–1% for alkali basalts and basanites, and solid mantle upwelling rates range from 40 to 100 cm yr⁻¹ for tholeiites and from 1 to 3 cm yr⁻¹ for basanites. For dynamic melting, the escape or threshold porosity is 0.5–2% for tholeiites and 0.1–0.8% for alkali basalts and basanites, and solid mantle upwelling rates range from 10 to 30 cm yr⁻¹ for tholeiites and from 0.1 to 1 cm yr⁻¹ for basanites. Assuming a constant melt productivity, calculated total melt fractions range from 15% for the tholeiitic basalts to 3% for alkali basalts and basanites.     

Subduction and melting processes inferred from U-Series, Sr-Nd-Pb isotope, and trace element data, Bicol and Bataan arcs, Philippines

We present U-series, Sr-Nd-Pb isotope, and trace element data from the two principal volcanic chains on Luzon Island, developed over oppositely dipping subduction zones, to explore melting and mass transfer processes beneath arcs. The Bataan (western) and Bicol (eastern) arcs are currently subducting terrigenous and pelagic sediments, respectively, which have different trace element and isotopic compositions. The range of (²³⁰Th/²³⁸U) disequilibria for both arcs is 0.85-1.15; only lavas from Mt. Mayon (Bicol arc) have ²³⁰Th activity excesses. Bataan lavas have higher ⁸⁷Sr/⁸⁶Sr and lower ¹⁴³Nd/¹⁴⁴Nd than Bicol lavas (⁸⁷Sr/⁸⁶Sr = 0.7042-0.7046, ¹⁴³Nd/¹⁴⁴Nd = 0.51281-0.51290 vs. ⁸⁷Sr/⁸⁶Sr = 0.70371-0.70391, ¹⁴³Nd/¹⁴⁴Nd = 0.51295-0.51301) and both arcs show steep linear arrays towards sediment values on ²⁰⁷Pb/²⁰⁴Pb vs. ²⁰⁶Pb/²⁰⁴Pb diagrams. Analysis of incompatible element and isotopic data allows identification of a sediment component that, at least in part, was transferred as a partial melt to the mantle wedge peridotite. Between 1% and 5% sediment melt addition can explain the isotopic and trace element variability in the rocks from both arcs despite the differences in sediment supply. We therefore propose that sediment transfer to the mantle wedge is likely mechanically or thermally limited. It follows that most sediments are either accreted, reside in the sub-arc lithosphere, or are recycled into the convecting mantle. However, whole-sale sediment recycling into the upper mantle is unlikely in light of the global mid-ocean ridge basalt data. Fluid involvement is more difficult to characterize, but overall the Bicol arc appears to have more fluid influence than the Bataan arc. Rock suites from each arc can be related by a dynamic melting process that allows for ²³⁰Th ingrowth, either by dynamic or continuous flux melting, provided the initial (²³⁰Th/²³²Th) of the source is ∼0.6-0.7. The implication of either model is that inclined arrays on the U-Th equiline diagram may not have chronologic significance. Modeling also suggests that U-series disequilibria are influenced by the tectonic convergence rate, which dictates mantle matrix flow. Thus with slower matrix flow there is a greater degree of ²³⁰Th ingrowth. While other factors such as prior mantle depletion and addition of a subducted component may explain some aspects of U-series data, an overall global correlation between tectonic convergence rate and the extent of U-Th disequilibria may originate from melting processes.     

Plume-lithosphere interaction beneath Mt. Cameroon volcano, West Africa: Constraints from U-Th-Ra and Sr-Nd-Pb isotope systematics

Precise measurements of ²³⁸U-²³⁰Th-²²⁶Ra disequilibria in lavas erupted within the last 100 yr on Mt. Cameroon are presented, together with major and trace elements, and Sr-Nd-Pb isotope ratios, to unravel the source and processes of basaltic magmatism at intraplate tectonic settings. All samples possess ²³⁸U-²³⁰Th-²²⁶Ra disequilibria with ²³⁰Th (18-24%) and ²²⁶Ra (9-21%) excesses, and there exists a positive correlation in a (²²⁶Ra/²³⁰Th)-(²³⁰Th/²³⁸U) diagram. The extent of ²³⁸U-²³⁰Th-²²⁶Ra disequilibria is markedly different in lavas of individual eruption ages, although the (²³⁰Th/²³²Th) ratio is constant irrespective of eruption age. When U-series results are combined with Pb isotope ratios, negative correlations are observed in the (²³⁰Th/²³⁸U)-(²⁰⁶Pb/²⁰⁴Pb) and (²²⁶Ra/²³⁰Th)-(²⁰⁶Pb/²⁰⁴Pb) diagrams. Shallow magma chamber processes like magma mixing, fractional crystallization and wall rock assimilation do not account for the correlations. Crustal contamination is not the cause of the observed isotopic variations because continental crust is considered to have extremely different Pb isotope compositions and U/Th ratios. Melting of a chemically heterogeneous mantle might explain the Mt. Cameroon data, but dynamic melting under conditions of high D_U and D_U/D_Th, long magma ascent time, or disequilibrium mineral/melt partitioning, is required. The most plausible scenario to produce the geochemical characteristics of Mt. Cameroon samples is the interaction of melt derived from the asthenospheric mantle with overlying sub-continental lithospheric mantle which has elevated U/Pb (>0.75) and Pb isotope ratios (²⁰⁶Pb/²⁰⁴Pb > 20.47) due to late Mesozoic metasomatism.     

Timescales of magma differentiation from basalt to andesite beneath Hekla Volcano, Iceland: Constraints from U-series disequilibria in lavas from the last quarter-millennium flows

Measurements of ²³⁸U-²³⁰Th-²²⁶Ra disequilibria, Sr-Nd-Pb-Hf isotopes and major-trace elements have been conducted for lavas erupted in the last quarter-millennium at Hekla volcano, Iceland. The volcanic rocks range from basalt to dacite. Most of the lavas (excluding dacitic samples) display limited compositional variations in radiogenic Sr-Nd-Pb-Hf isotopes (⁸⁷Sr/⁸⁶Sr = 0.70319-0.70322; ¹⁴³Nd/¹⁴⁴Nd = 0.51302-0.51305; ²⁰⁶Pb/²⁰⁴Pb = 19.04-19.06; ²⁰⁷Pb/²⁰⁴Pb = 15.53-15.54; ²⁰⁸Pb/²⁰⁴Pb = 38.61-38.65; ¹⁷⁶Hf/¹⁷⁷Hf = 0.28311-0.28312). All the samples possess (²³⁰Th/²³⁸U) disequilibrium with ²³⁰Th excesses, and they show systematic variations in (²³⁰Th/²³²Th) and (²³⁸U/²³²Th) ratios. The highest ²²⁶Ra excesses occur in the basalt and most differentiated andesite lavas, while some basaltic-andesite lavas have (²²⁶Ra/²³⁰Th) ratio that are close to equilibrium. The ²³⁸U-²³⁰Th-²²⁶Ra disequilibria variations cannot be produced by simple closed-system fractional crystallization with radioactive decay of ²³⁰Th and ²²⁶Ra in a magma chamber. A closed-system fractional crystallization model and assimilation and fractional crystallization (AFC) model indicate that the least differentiated basaltic andesites were derived from basalt by fractional crystallization with a differentiation age of ∼24 ± 11 kyr, whereas the andesites were formed by assimilation of crustal material and fractionation of the basaltic-andesites within 2 kyr. Apatite is inferred to play a key role in fractionating the parent-daughter nuclides in ²³⁰Th-²³⁸U and ²²⁶Ra-²³⁰Th to make the observed variations. Our proposed model is that several batches of basaltic-andesite magmas that formed by fractional crystallization of a basaltic melt from a deeper reservoir, were periodically injected into the shallow crust to form individual magma pockets, and subsequently modifying the original magma compositions via simultaneous assimilation and fractional crystallization. The assimilant is the dacitic melt, which formed by partial melting of the crust.     

Understanding melt generation beneath the slow-spreading Kolbeinsey Ridge using U, Th, and Pa excesses

To examine the petrogenesis and sources of basalts from the Kolbeinsey Ridge, one of the shallowest locations along the global ridge system, we present new measurements of Nd, Sr, Hf, and Pb isotopes and U-series disequilibria on 32 axial basalts. Young Kolbeinsey basalts (full-spreading rate = 1.8 cm/yr; 67°05′-70°26′N) display (²³⁰Th/²³⁸U) < 1 and (²³⁰Th/²³⁸U) > 1 with (²³⁰Th/²³⁸U) from 0.95 to 1.30 and have low U (11.3-65.6 ppb) and Th (33.0 ppb-2.40 ppm) concentrations. Except for characteristic isotopic enrichment near the Jan Mayen region, the otherwise depleted Kolbeinsey basalts (e.g. ⁸⁷Sr/⁸⁶Sr = 0.70272-0.70301, ε_Nd = 8.4-10.5, ε_Hf = 15.4-19.6 (La/Yb)_N = 0.28-0.84) encompass a narrow range of (²³⁰Th/²³²Th) (1.20-1.32) over a large range in (²³⁸U/²³²Th) (0.94-1.32), producing a horizontal array on a (²³⁰Th/²³²Th) vs. (²³⁸U/²³²Th) diagram and a large variation in (²³⁰Th/²³⁸U). However, the (²³⁰Th/²³⁸U) of the Kolbeinsey Ridge basalts (0.96-1.30) are inversely correlated with (²³⁴U/²³⁸U) (1.001-1.031). Samples with low (²³⁰Th/²³⁸U) and elevated (²³⁴U/²³⁸U) reflect alteration by seawater or seawater-derived materials. The unaltered Kolbeinsey lavas with equilibrium ²³⁴U/²³⁸U have high (²³⁰Th/²³⁸U) values (?1.2), which are consistent with melting in the presence of garnet. This is in keeping with the thick crust and anomalously shallow axial depth for the Kolbeinsey Ridge, which is thought to be the product of large degrees of melting in a long melt column. A time-dependent, dynamic melting scenario involving a long, slowly upwelling melting column that initiates well within the garnet peridotite stability zone can, in general, reproduce the (²³⁰Th/²³⁸U) and (²³¹Pa/²³⁵U) ratios in uncontaminated Kolbeinsey lavas, but low (²³¹Pa/²³⁵U) ratios in Eggvin Bank samples suggest eclogite involvement in the source for that ridge segment.     

U-TH-PA-RA study of the Kamchatka arc: new constraints on the genesis of arc lavas

The ²³⁸U-²³⁰Th-²²⁶Ra and ²³⁵U-²³¹Pa disequilibria have been measured by mass spectrometry in historic lavas from the Kamchatka arc. The samples come from three closely located volcanoes in the Central Kamchatka Depression (CKD), the most active region of subducted-related volcanism in the world. The large excesses of ²²⁶Ra over ²³⁰Th found in the CKD lavas are believed to be linked to slab dehydration. Moreover, the samples show the uncommon feature of (²³⁰Th/²³⁸U) activity ratios both lower and higher than 1. The U-series disequilibria are characterized by binary trends between activity ratios, with (²³¹Pa/²³⁵U) ratios all >1. It is shown that these correlations cannot be explained by a simple process involving a combination of slab dehydration and melting. We suggest that they are more likely to reflect mixing between two end-members: a high-magnesia basalt (HMB) end-member with a clear slab fluid signature and a high-alumina andesite (HAA) end-member reflecting the contribution of a slab-derived melt. The U-Th-Ra characteristics of the HMB end-member can be explained either by a two-step fluid addition with a time lag of 150 ka between each event or by continuous dehydration. The inferred composition for the dehydrating slab is a phengite-bearing eclogite. Equilibrium transport or dynamic melting can both account for ²³¹Pa excess over ²³⁵U in HMB end-member. Nevertheless, dynamic melting is preferred as equilibrium transport melting requires unrealistically high upwelling velocities to preserve fluid-derived ²²⁶Ra/²³⁰Th. A continuous flux melting model is also tested. In this model, ²³¹Pa-²³⁵U is quickly dominated by fluid addition and, for realistic extents of melting, this process cannot account for (²³¹Pa/²³⁵U) ratios as high as 1.6, as observed in the HMB end-member.The involvement of a melt derived from the subducted oceanic crust is more likely for explaining the HAA end-member compositions than crustal assimilation. Melting of the oceanic crust is believed to occur in presence of residual phengite and rutile, resulting in no ²²⁶Ra-²³⁰Th disequilibrium and low ²³¹Pa excess over ²³⁵U in the high-alumina andesites. Consequently, it appears that high-alumina andesites and high-magnesia basalts have distinct origins: the former being derived from melting of the subducted oceanic crust and the latter from hydrated mantle. It seems that there is no genetic link between these two magma types, in contrast with what was previously believed.     

U-series disequilibria in suspended river sediments and implication for sediment transfer time in alluvial plains: The case of the Himalayan rivers

²³⁸U-²³⁴U-²³⁰Th radioactive disequilibria were analyzed in suspended sediments (collected at different depths) from the Ganges River and one of its main tributaries: the Narayani-Gandak River. Results associated with bedload sediment data suggest that uranium-series (U-series) disequilibria in river sediments of the Ganges basin vary with grain size and sampling location. The range of observed U-series disequilibria is explained by a mixing model between a coarse-grained sediment end-member, represented by bedload and bank sediments, and a fine-grained end-member that both originate from Himalaya but undergo different transfer histories within the plain. The coarse-grained sediment end-member transits slowly (i.e. >several 100’s ky) in the plain whereas the fine-grained sediment end-member is transferred much faster (<20-25 ky), as indicated by the absence of significant variations in Th isotope composition of the fine-grained sediment end-members. These results show that U-series isotopes can be used to quantify the various transfer times of river sediments of different sizes and infer that there can be an order of magnitude of difference, or more, between the transfer time of suspended and bedload sediments. This underlines that a good knowledge of the proportion of suspended vs. bedload sediments transported in the river is required to accurately assess how fast erosion products are transferred in a catchment and how fast a catchment is likely to respond to external forcing factors.     

U-series disequilibria in Kick’em Jenny submarine volcano lavas: A new view of time-scales of magmatism in convergent margins

We present data for U and its decay series nuclides ²³⁰Th, ²²⁶Ra, ²³¹Pa, and ²¹⁰Po for 14 lavas from Kick’em Jenny (KEJ) submarine volcano to constrain the time-scales and processes of magmatism in the Southern Lesser Antilles, the arc having the globally lowest plate convergence rate. Although these samples are thought to have been erupted in the last century, most have (²²⁶Ra)/(²¹⁰Po) within ±15% of unity. Ten out of 14 samples have significant ²²⁶Ra excesses over ²³⁰Th, with (²²⁶Ra)/(²³⁰Th) up to 2.97, while four samples are in ²²⁶Ra-²³⁰Th equilibrium within error. All KEJ samples have high (²³¹Pa)/(²³⁵U), ranging from 1.56 to 2.64 and high ²³⁸U excesses (up to 43%), providing a global end-member of high ²³⁸U and high ²³¹Pa excesses. Negative correlations between Sr, sensitive to plagioclase fractionation, and Ho/Sm, sensitive to amphibole fractionation, or K/Rb, sensitive to open system behavior, indicate that differentiation at KEJ lavas was dominated by amphibole fractionation and open-system assimilation. While (²³¹Pa)/(²³⁵U) does not correlate with differentiation indices such as Ho/Sm, (²³⁰Th)/(²³⁸U) shows a slight negative correlation, likely due to assimilation of materials with slightly higher (²³⁰Th)/(²³⁸U). Samples with ²²⁶Ra excess have higher Sr/Th and Ba/Th than those in ²²⁶Ra-²³⁰Th equilibrium, forming rough positive correlations of (²²⁶Ra)/(²³⁰Th) with Sr/Th and Ba/Th similar to those observed in many arc settings. We interpret these correlations to reflect a time-dependent magma differentiation process at shallow crustal levels and not the process of recent fluid addition at the slab-wedge interface.The high ²³¹Pa excesses require an in-growth melting process operating at low melting rates and small residual porosity; such a model will also produce significant ²³⁸U-²³⁰Th and ²²⁶Ra-²³⁰Th disequilibrium in erupted lavas, meaning that signatures of recent fluid addition from the slab are unlikely to be preserved in KEJ lavas. We instead propose that most of the ²³⁸U-²³⁰Th, ²²⁶Ra-²³⁰Th, and ²³⁵U-²³¹Pa disequilibria in erupted KEJ lavas reflect the in-growth melting process in the mantle wedge (reflecting variations in U/Th, daughter-parent ratios, fO₂, and thermal structure), followed by modification by magma differentiation at crustal depths. Such a conclusion reconciles the different temporal implications from different U-series parent-daughter pairs and relaxes the time constraint on mass transfer from slab to eruption occurring in less than a few thousand years imposed by models whereby ²²⁶Ra excess is derived from the slab.     

Partial melting processes above subducting plates: Constraints from Pa-U disequilibria

The processes involved in the formation and transport of partial melts above subducting plates remain poorly constrained relative to those at mid-ocean ridges. In particular, ²³⁸U-²³⁰Th-²²⁶Ra disequilibria, that might normally be used to constrain melting dynamics, tend to be swamped by the effects of fluid addition from the down-going plate. The ²³¹Pa-²³⁵U system provides an exciting exception to this because the highly incompatible nature of Pa means that fractionation and in-growth during partial melting overwrite the effects of fluid U addition. We present ²³¹Pa-²³⁵U data on 50 well-characterised lavas from seven subduction zones in order to examine partial melting processes. Measured (²³¹Pa/²³⁵U) ratios are all >1 and 15% are >2. Overall (²³¹Pa/²³⁵U) shows broad positive correlations with (²³⁰Th/²³⁸U) and La/Yb and negative trends against Ba/Th and (²²⁶Ra/²³⁰Th). These systematics can differ from arc to arc but suggest that (²³¹Pa/²³⁵U) tends to be higher in sediment-rich arc lavas where the effects of fluid addition are muted and there is less of a ²³¹Pa deficit for melting to overprint. We have explored the effects of decompression melting, frictional drag dynamic melting with and without ageing subsequent to fluid U addition to the wedge as well as flux melting models. Globally, average (²³¹Pa/²³⁵U) appears to correlate negatively with convergence rate and so in the numerical models we use the local subduction rate for the rate of matrix flow through the melting zone. Using this assumption and reasonable values for other parameters, the melting models can simulate the overall range of (²³¹Pa/²³⁵U) and some of the data trends. However, it is clear that local variations in some parameters, especially source composition and extent of melting, exert a major influence on ²³¹Pa-²³⁵U disequilibria. Some data, which lie at a high angle to the modelled trends, may be explained by mixing between small degree hydrous melts formed near the slab and larger degree, decompression melts produced at shallow depth.     

Chemical and isotopic constraints on the generation and transport of magma beneath the East Pacific Rise

Interpretation of U-series disequilibria in midocean ridge basalts is highly dependent on the bulk partition coefficients for U and Th and therefore the mineralogy of the mantle source. Distinguishing between the effect of melting processes and variable source compositions on measured disequilibria (²³⁸U-²³⁰Th-²²⁶Ra and ²³⁵U-²³¹Pa) requires measurement of the radiogenic isotopes Hf, Nd, Sr, and Pb. Here, we report measurements of ²³⁸U-²³⁰Th-²²⁶Ra and ²³⁵U-²³¹Pa disequilibria; Hf, Nd, Sr, and Pb isotopic; and major and trace element compositions for a suite of 20 young midocean ridge basalts from the East Pacific Rise axis between 9°28′ and 9°52′N. All of the samples were collected within the axial summit trough using the submersible Alvin. The geological setting and observational data collected during sampling operations indicate that all the rocks are likely to have been erupted from 1991 to 1992 or within a few decades of that time. In these samples, ²³⁰Th excesses and ²²⁶Ra excesses are variable and inversely correlated. Because the eruption ages of the samples are much less than the half-life of ²²⁶Ra, this inverse correlation between ²³⁰Th and ²²⁶Ra excesses can be considered a primary feature of these lavas. For the lava suite analyzed in this study, ²²⁶Ra and ²³⁰Th excesses also vary with lava composition: ²²⁶Ra excesses are negatively correlated with Na₈ and La/Yb and positively correlated with Mg#. Conversely, ²³⁰Th excesses are positively correlated with Na₈ and La/Yb and negatively correlated with Mg#. Th/U, ²³⁰Th/²³²Th, and ²³⁰Th excesses are also variable and correlated to one another. ²³¹Pa excesses are large but relatively constant and independent of Mg#, La/Yb, Th/U, and Na₈. The isotope ratios ¹⁴³Nd/¹⁴⁴Nd, ¹⁷⁶Hf/¹⁷⁷Hf, ⁸⁷Sr/⁸⁶Sr, and ²⁰⁸Pb/²⁰⁶Pb are constant within analytical uncertainty, indicating that they were derived from a common source. The source is homogeneous with respect to parent/daughter ratios Lu/Hf, Sm/Nd, Rb/Sr, and Th/U; therefore, the measured variations of Th/U, ²³⁰Th, and ²²⁶Ra excesses and major and trace element compositions in these samples are best explained by polybaric melting of a homogeneous source, not by mixing of compositionally distinct sources.     

Time scale and conditions of weathering under tropical climate: Study of the Amazon basin with U-series

The Rio Solimões/Amazonas (Amazon River) and its major tributaries have been analyzed for U-series nuclides. ²³⁸U-²³⁴U-²³⁰Th-²²⁶Ra disequilibria have been measured in the dissolved (<0.2 μm) and suspended loads (>0.2 μm) as well as bed sands. U-series disequilibria are closely related to major and trace element compositions and therefore reflect elemental fractionation during chemical weathering. Moreover, while the dissolved load records present-day weathering, suspended particles integrate the erosion history over much longer time scales (>100 ka). Lowland rivers are characterized by long time scales of chemical erosion (?100 ka) resulting in a high weathering intensity. Moreover, exchange between suspended particles and the dissolved load may explain the U-series signature for these rivers. By combining U-series and Pb isotopes in suspended particles, we show that erosion in the Rio Madeira basin occurred as a multi-step process, whereby the pristine continental crust was eroded several hundreds of Ma ago to produce sediments that have then been integrated in the Cordillera by crustal shortening and are currently eroded. In contrast, recent erosion of a pristine crust is more likely for the Rio Solimões/Amazonas (<10 ka). The suspended particles of the rivers draining the Andes (Solimões/Amazonas, Madeira) suggest time scales of weathering ranging between 4 and 20 ka. This indicates that suspended particles transported by those rivers are not stored for long periods in the Andean foreland basin and the tropical plain. The sediments delivered to the ocean have resided only a few ka in the Amazon basin (6.3 ± 1 ka for the Rio Amazonas at Óbidos). Nevertheless, a large fraction of the sediments coming out from the Andes are trapped in the foreland basin and may never reach the ocean. Erosion in the Andes is not operating in steady state. U-series systematics shows unambiguously that rivers are exporting a lot more sediments than predicted by steady-state erosion and that is a consequence of soil destruction greater than production. By relating this observation to the short time scales of weathering inferred for the Andes (a few ka), it appears that the erosion regime has been recently perturbed, resulting in high denudation rates. A possible explanation would be the increase in precipitation less than 5 ka proposed by recent paleoclimatic studies. Our results indicate that erosion responds rapidly to high-frequency climatic fluctuations.     

U-series disequilibria in volcanic rocks from the Canary Islands: Plume versus lithospheric melting

U-series disequilibria are presented for Holocene samples from the Canary Islands and interpreted with special emphasis on the separate roles of plume vs. lithospheric melting processes. We report Th and U concentrations and (²³⁸U)/(²³²Th), (²³⁰Th)/(²³²Th), (²³⁰Th)/(²³⁸U) and (²³⁴U)/(²³⁸U) for 43 samples, most of which are minimally differentiated, along with (²²⁶Ra)/(²³⁰Th) and (²³¹Pa)/(²³⁵U) for a subset of these samples, measured by thermal ionization mass spectrometry (TIMS). Th and U concentrations range between 2 and 20 ppm and 0.5 and 6 ppm, respectively. Initial (²³⁰Th)/(²³⁸U) ranges from 1.1 to 1.6. (²²⁶Ra)/(²³⁰Th)_o ranges between 0.9 and 1.8 while (²³¹Pa)/(²³⁵U)_o ranges between 1.0 and 2.0.Our interpretation of results is based on a three-fold division of samples as a function of incompatible element ratio, such as Nb/U. The majority of samples have Nb/U = 47 ± 10, similar to most MORB and OIB. Higher ratios are found exclusively in alkali basalts and tholeiites from the eastern Canary Islands whereas lower ratios are exclusively found in differentiated rocks from the western Canary Islands. Those with ordinary Nb/U ratios are attributed to melting within the slowly ascending HIMU-dominated Canary plume.Higher Nb/U, generally found in more silica rich basalts from the eastern islands, is attributed to lithospheric contamination. Based on their trace element characteristics, two possible contaminants are amphibole veins (± other minerals) crystallized in the mantle from previous plume-derived basanite or re-melted plume-derived intrusive rocks. The high Nb/U signature of these materials is imparted on a melt of the lithosphere created either by the diffusive infiltration of alkalis or by direct reaction between basanites and peridotite. Mixing between plume-derived basanite and lithospheric melt accounts for the U-series systematics of most eastern island magmas including the well-known Timanfaya eruption. Lower Nb/U ratios in differentiated rocks from the western islands are attributed to fractional crystallization of amphibole ± phlogopite ± sphene from basanite during its ascent through the lithosphere. Based on changes in disequilibria, phonolites and tephrites are interpreted to result from rapid differentiation of primitive parents within millennia.     

The extent of U-series disequilibria produced during partial melting of the lower crust with implications for the formation of the Mount St. Helens dacites

The extent to which U-series disequilibria can be produced during partial melting of mafic lower crust is quantified using a simple batch melting model and both experimental and theoretical partition coefficients for U, Th and Ra. We show that partial melting of mafic lower crust can only produce small disequilibria between ²³⁸U, ²³⁰Th and ²²⁶Ra. Crystallisation of basalt and mixing between young basalt and crustal derived melts will have a similar or smaller effect. Consequently, U-series disequilibrium in arc andesites and dacites can generally only be an inherited feature derived from a mantle parent, unless the timescales of silicic magma production within the crust are short compared to the half-life of ²²⁶Ra. Our results have profound implications for several recent models of silicic magma production by thermal incubation and partial melting of the lower crust. We show that the ²²⁶Ra excess observed in most arc andesites and dacites requires extremely rapid differentiation and/or the involvement of mantle derived basalts less than a few thousand years old. Application to Mount St. Helens suggests that crystallisation of young mantle-derived magma is likely to be the dominant process in the formation of these dacites.     

(Pa/U)-(Th/U) of young mafic volcanic rocks from Nicaragua and Costa Rica and the influence of flux melting on U-series systematics of arc lavas

We present U, Th, and Pa isotope data for young lavas from Costa Rica and Nicaragua in the Central American arc. Thorium isotopic ratios for Costa Rica and Nicaragua differ dramatically: Costa Rican lavas are characterized by low (²³⁰Th/²³²Th) (1 to 1.2) and, for four out of five lavas, (²³⁰Th/²³⁸U) greater than unity. Nicaraguan lavas have high (²³⁰Th/²³²Th) (2.2 to 2.7) and, for five of six samples, (²³⁰Th/²³⁸U) less than unity. All lavas have (²³¹Pa/²³⁵U) greater than unity, with initial values ranging from 1.27 to 1.77, but those from Costa Rica have larger ²³¹Pa excesses. There is a broad positive correlation between (²³¹Pa/²³⁵U) and (²³⁰Th/²³⁸U) similar to the worldwide trend for arcs outlined by Pickett and Murrell (1997), although many of the Nicaraguan lavas skirt the high end of that trend. In greater detail, the Central American data appear to divide into separate high-(²³¹Pa/²³⁵U) and low-(²³¹Pa/²³⁵U) tiers. These tiers may be different because of either different residence times in the crust or different proportions of sedimentary components from the slab.Substantial (²³¹Pa/²³⁵U) excesses (>1.5) in both Costa Rica and Nicaragua require a melting process that allows for enhanced daughter (²³¹Pa) ingrowth. With increasing U addition, (²³¹Pa/²³⁰Th) increases in a manner that cannot be explained adequately by aging of fluid components before partial melting and eruption. Thus, either some ²³¹Pa is added from the slab, or melting-enhanced ²³¹Pa ingrowth is greater in sources that have experienced a larger amount of slab-derived flux and a higher extent of melting. These observations can be explained if regions that have undergone greater extents of fluxing and melting have experienced these processes over a longer time interval than those that have had little flux added and little melt extracted. We propose a flux-ingrowth melting model in which corner flow in the mantle wedge supplies fresh hot mantle into a zone of slab fluid addition. Partial melting occurs in response to this fluxing. We assume critical melting at low porosity (∼10⁻³), rapid fluid flux to the melting region, and rapid melt transport. Solid mantle traverses the melting region over 10⁵ to 10⁶ yr, thereby allowing ²³¹Pa and ²³⁰Th ingrowth from U retained in the residues of melt extraction. Magmas are aggregated from all parts of the melting regime, mixing melts from incipiently fluxed regions with those from sources that have experienced more extensive fluid addition, partial melting, and daughter nuclide ingrowth. With suitable assumptions about component addition from the slab, this flux-ingrowth model matches a wide range of U-series and trace element data from Costa Rican and Nicaraguan lavas, with required average extents of melting of ∼1 to 3% and 7 to 15%, respectively. Upwelling and/or extensive melt-rock reaction are not required to explain large (²³¹Pa/²³⁵U) excesses in Central America or other arcs. On Th isotope equiline plots, the model produces linear arrays that resemble isochrons but that have no age significance. Instead, these arrays are generated by mixing of melts from sources that have experienced fluid addition and partial melting over a range of time intervals, as seems likely in arc source regions. Finally, the flux-ingrowth model predicts considerable ²²⁶Ra excesses for integrated magmas. If we assume that ²²⁶Ra is added continuously with the slab-derived fluid, the model predicts large and increasing (²²⁶Ra/²³⁰Th) with increasing melting and slab-component addition, without requiring the addition of a distinct late fluid.     

Geochemical evolution of a shallow magma plumbing system during the last 500 years, Miyakejima volcano, Japan: Constraints from U-Th-Ra systematics

In order to unravel magma processes and the geochemical evolution of shallow plumbing systems beneath active volcanoes, we investigated U-series disequilibria of rocks erupted over the past 500 years (1469-2000 AD) from Miyakejima volcano, Izu arc, Japan. Miyakejima volcanic rocks show ²³⁸U-²³⁰Th-²²⁶Ra disequilibria with excess ²³⁸U and ²²⁶Ra, due to the addition of slab-derived fluids to the mantle wedge. Basaltic bombs of the 2000 AD eruption have the lowest (²³⁰Th/²³²Th) ratio compared to older Miyakejima eruptives, yielding the youngest ²³⁸U-²³⁰Th model age of 2 kyr. This reinforces our previous model that fluid release from the slab and subsequent magma generation in the mantle wedge beneath Miyakejima occur episodically on a several-kyr timescale. In the last 500 years, Miyakejima eruptives show: (1) a vertical trend in a (²³⁰Th/²³²Th)-(²³⁸U/²³²Th) diagram and (2) a positive linear correlation in a (²²⁶Ra/²³⁰Th)₀ − 1/²³⁰Th diagram, which is also observed in lavas from some of the single eruptions (e.g., 1940, 1962, and 1983 AD). The variations cannot be produced by simple fractional crystallization in a magma chamber with radioactive decay of ²³⁰Th and ²²⁶Ra, but it is possibly produced by synchronous generation of melts in the mantle wedge with different upwelling rate or addition of multiple slab-derived fluids. A much more favorable scenario is that some basaltic magmas were intermittently supplied from deep in the mantle and injected into the crust, subsequently modifying the original magma composition and producing variations in (²³⁰Th/²³²Th) and (²²⁶Ra/²³⁰Th)₀ ratios via assimilation and fractional crystallization (AFC). The assimilant of the AFC process would be a volcanic edifice of previous Miyakejima magmatism. Due to the relatively short timescales involved, the interaction between the assimilant and recent Miyakejima magmatism has not been recorded by the Sr-Nd-Pb isotopic systems. In such cases, Th isotopes and (²²⁶Ra/²³⁰Th) ratio are excellent geochemical tracers of magmatic evolution.     

Diffusive fractionation of U-series radionuclides during mantle melting and shallow-level melt-cumulate interaction

U-series radioactive disequilibria in basaltic lavas have been used to infer many important aspects of melt generation and extraction processes in Earth’s mantle and crust, including the porosity of the melting zone, the solid mantle upwelling rate, and the melt transport rate. Most of these inferences have been based on simplified theoretical treatments of the fractionation process, which assume equilibrium partitioning of U-series nuclides among minerals and melt. We have developed a numerical model in which solid-state diffusion controls the exchange of U-series nuclides among multiple minerals and melt. First the initial steady-state distribution of nuclides among the phases, which represents a balance between diffusive fluxes and radioactive production and decay, is calculated. Next, partial melting begins, or a foreign melt is introduced into the system, and nuclides are again redistributed among the phases via diffusion. U-series nuclides can be separated during this stage due to differences in their diffusivity; radium in particular, and possibly protactinium as well, can be strongly fractionated from slower-diffusing thorium and uranium. We show that two distinct processes are not required for the generation of ²²⁶Ra and ²³⁰Th excesses in mid-ocean ridge basalts, as has been argued previously; instead the observed negative correlations of the (²²⁶Ra/²³⁰Th) activity ratio with (²³⁰Th/²³⁸U) and with the extent of trace element enrichment may result from diffusive fractionation of Ra from Th during partial melting of the mantle. Alternatively, the (²²⁶Ra/²³⁰Th) disequilibrium in mid-ocean ridge basalts may result from diffusive fractionation during shallow-level interaction of mantle melts with gabbroic cumulates, and we show that the results of the interaction have a weak dependence on the age of the cumulate if both plagioclase and clinopyroxene are present.     

Geochemical evolution of historical lavas from Askja Volcano, Iceland: Implications for mechanisms and timescales of magmatic differentiation

The mechanisms and the timescales of magmatic evolution were investigated for historical lavas from the Askja central volcano in the Dyngjufjöll volcanic massif, Iceland, using major and trace element and Sr, Nd, and Pb isotopic data, as well as ²³⁸U-²³⁰Th-²²⁶Ra systematics. Lavas from the volcano show marked compositional variation from magnesian basalt through ferrobasalt to rhyolite. In the magnesian basalt-ferrobasalt suite (5-10 wt% MgO), consisting of lavas older than 1875 A.D., ⁸⁷Sr/⁸⁶Sr increases systematically with increasing SiO₂ content; this suite is suggested to have evolved in a magma chamber located at ∼600 MPa through assimilation and fractional crystallization. On the other hand, in the ferrobasalt-rhyolite suite (1-5 wt% MgO), including 1875 A.D. basalt and rhyolite and 20th century lavas, ⁸⁷Sr/⁸⁶Sr tends to decrease slightly with increasing SiO₂ content. It is suggested that a relatively large magma chamber occupied by ferrobasalt magma was present at ∼100 MPa beneath the Öskjuvatn caldera, and that icelandite and rhyolite magmas were produced by extraction of the less and more evolved interstitial melt, respectively, from the mushy boundary layer along the margin of the ferrobasalt magma chamber, followed by accumulation of the melt to form separate magma bodies. Ferrobasalt and icelandite lavas in the ferrobasalt-rhyolite suite have a significant radioactive disequilibrium in terms of (²²⁶Ra/²³⁰Th), and its systematic decrease with magmatic evolution is considered to reflect aging, along with assimilation and fractional crystallization processes. Using a mass-balance model in which simultaneous fractional crystallization, crustal assimilation, and radioactive decay are taken into account, the timescale for the generation of icelandite magma from ferrobasalt was constrained to be <∼3 kyr which is largely dependent on Ra crystal-melt partition coefficients we used.     

Mantle dynamics and mantle melting beneath Niuafo’ou Island and the northern Lau back-arc basin

A suite of young volcanic basaltic lavas erupted on the intra-plate island of Niuafo’ou and at active rifts and spreading centres (the King’s Triple Junction and the Northeastern Lau Spreading Centre) in the northern Lau Basin is used to examine the pattern of mantle flow and the dynamics of melting beneath this complex back-arc system. All lavas contain variable amounts of a subduction related component inherited from the Tonga subduction zone to the east. All lavas have higher ⁸⁷Sr/⁸⁶Sr, lower ¹⁴³Nd/¹⁴⁴Nd and more radiogenic Pb isotope compositions than basalts erupted at the Central Lau Spreading Centre in the central Lau Basin, and are interpreted as variable mixtures of subduction-modified, depleted upper mantle, and mantle residues derived from melting beneath the Samoan Islands which has leaked through a tear in the subducting Pacific Plate beneath the Vitiaz Lineament at the northern edge of the Lau Basin. Our data can be used to map out the present-day distribution of Samoan mantle in this region, and show that it influences the compositions of lavas erupted as far as 400 km from the Samoan Islands. The distribution of Samoan-influenced lavas implies south- and southwest-wards mantle flow rates of >4 cm/year. U-series disequilibria in historic Niuafo’ou lavas have average (²³⁰Th/²³⁸U) = 1.13, (²³¹Pa/²³⁵U) = 2.17, (²²⁶Ra/²³⁰Th) = 2.11, and together with major and trace element data require ∼5% partial melting of mantle at between 2 and 3 GPa, with a residual porosity of 0.002 and an upwelling rate of 1 cm year⁻¹. We suggest that intraplate magmatism in the northern Lau Basin results from decompression melting during southward flow of mantle from beneath old (110–120 Ma), relatively thick Pacific oceanic lithosphere to beneath young (<5 Ma), thinner oceanic lithosphere beneath the northern Lau Basin.     

A possible role for garnet pyroxenite in the origin of the “garnet signature” in MORB

 Geochemical data have been interpreted as requiring that a significant fraction of the melting in MORB source regions takes place in the garnet peridotite field, an inference that places the onset of melting at ≥80 km. However, if melting begins at such great depths, most models for melting of the suboceanic mantle predict substantially more melting than that required to produce the 7±1 km thickness of crust at normal ridges. One possible resolution of this conflict is that MORBs are produced by melting of mixed garnet pyroxenite/spinel peridotite sources and that some or all of the “garnet signature” in MORB is contributed by partial melting of garnet pyroxenite layers or veins, rather than from partial melting of garnet peridotite. Pyroxenite layers or veins in peridotite will contribute disproportionately to melt production relative to their abundance, because partial melts of pyroxenite will be extracted from a larger part of the source region than peridotite partial melts (because the solidus of pyroxenite is at lower temperature than that of peridotite and is encountered along an adiabat 15–25 km deeper than the solidus of peridotite), and because melt productivity from pyroxenite during upwelling is expected to be greater than that from peridotite (pyroxenite melt productivity will be particularly high in the region before peridotite begins melting, owing to heating from the enclosing peridotite). For reasonable estimates of pyroxenite and peridotite melt productivities, 15–20% of the melt derived from a source region composed of 5% pyroxenite and 95% peridotite will come from the pyroxenite. Most significantly, garnet persists on the solidus of pyroxenite to much lower pressures than those at which it is present on the solidus of peridotite, so if pyroxenite is present in MORB source regions, it will probably contribute a garnet signature to MORB even if melting only occurs at pressures at which the peridotite is in the spinel stability field. Partial melting of a mixed spinel peridotite/garnet pyroxenite mantle containing a few to several percent pyroxenite can explain quantitatively many of the geochemical features of MORB that have been attributed to the onset of melting in the stability field of garnet lherzolite, provided that the pyroxenite compositions are similar to the average composition of mantle-derived pyroxene-rich rocks worldwide or to reasonable estimates of the composition of subducted oceanic crust. Sm/Yb ratios of average MORB from regions of typical crustal thickness are difficult to reconcile with derivation by melting of spinel peridotite only, but can be explained if MORB sources contain ∼5% garnet pyroxenite. Relative to melting of spinel peridotite alone, participation of model pyroxenite in melting lowers aggregate melt Lu/Hf without changing Sm/Nd ratios appreciably. Lu/Hf-Sm/Nd systematics of most MORB can be accounted for by melting of a spinel peridotite/garnet pyroxenite mantle provided that the source region contains 3–6% pyroxenite with ≥20% modal garnet. However, Lu/Hf-Sm/Nd systematics of some MORB appear to require more complex melting regimes and/or significant isotopic heterogeneity in the source. Another feature of the MORB garnet signature, (²³⁰Th)/(²³⁸U)>1, can also be produced under these conditions, although the magnitude of (²³⁰Th)/(²³⁸U) enrichment will depend on the rate of melt production when the pyroxenite first encounters the solidus, which is not well-constrained. Preservation of high (²³⁰Th)/(²³⁸U) in aggregated melts of mixed spinel peridotite/garnet pyroxenite MORB sources is most likely if the pyroxenites have U concentrations similar to that expected in subducted oceanic crust or to pyroxenite from alpine massifs and xenoliths. The abundances of pyroxenite in a mixed source that are required to explain MORB Sm/Yb, Lu/Hf, and (²³⁰Th)/(²³⁸U) are all similar. If pyroxenite is an important source of garnet signatures in MORB, then geochemical indicators of pyroxenite in MORB source regions, such as increased trace element and isotopic variability or more radiogenic Pb or Os, should correlate with the strength of the garnet signature. Garnet signatures originating from melts of the garnet pyroxenite components of mixed spinel peridotite/garnet pyroxenite sources would also be expected to be stronger in regions of thin crust. Received: 15 February 1995/Accepted: 7 February 1996     

Combined U-Th/He and Ar/Ar geochronology of post-shield lavas from the Mauna Kea and Kohala volcanoes, Hawaii

Late Quaternary, post-shield lavas from the Mauna Kea and Kohala volcanoes on the Big Island of Hawaii have been dated using the ⁴⁰Ar/³⁹Ar and U-Th/He methods. The objective of the study is to compare the recently demonstrated U-Th/He age method, which uses basaltic olivine phenocrysts, with ⁴⁰Ar/³⁹Ar ages measured on groundmass from the same samples. As a corollary, the age data also increase the precision of the chronology of volcanism on the Big Island. For the U-Th/He ages, U, Th and He concentrations and isotopes were measured to account for U-series disequilibrium and initial He. Single analyses U-Th/He ages for Hamakua lavas from Mauna Kea are 87 ± 40 to 119 ± 23 ka (2σ uncertainties), which are in general equal to or younger than ⁴⁰Ar/³⁹Ar ages. Basalt from the Polulu sequence on Kohala gives a U-Th/He age of 354 ± 54 ka and a ⁴⁰Ar/³⁹Ar age of 450 ± 40 ka. All of the U-Th/He ages, and all but one spurious ⁴⁰Ar/³⁹Ar ages conform to the previously proposed stratigraphy and published ¹⁴C and K-Ar ages. The ages also compare favorably to U-Th whole rock-olivine ages calculated from ²³⁸U-²³⁰Th disequilibria. The U-Th/He and ⁴⁰Ar/³⁹Ar results agree best where there is a relatively large amount of radiogenic ⁴⁰Ar (>10%), and where the ⁴⁰Ar/³⁶Ar intercept calculated from the Ar isochron diagram is close to the atmospheric value. In two cases, it is not clear why U-Th/He and ⁴⁰Ar/³⁹Ar ages do not agree within uncertainty. U-Th/He and ⁴⁰Ar/³⁹Ar results diverge the most on a low-K transitional tholeiitic basalt with abundant olivine. For the most alkalic basalts with negligible olivine phenocrysts, U-Th/He ages were unattainable while ⁴⁰Ar/³⁹Ar results provide good precision even on ages as low as 19 ± 4 ka. Hence, the strengths and weaknesses of the U-Th/He and ⁴⁰Ar/³⁹Ar methods are complimentary for basalts with ages of order 100-500 ka.