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.     

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.     

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     

Distinguishing lower and upper crustal processes in magmas erupted during the buildup to the 7.7 ka climactic eruption of Mount Mazama, Crater Lake, Oregon, using 238U–230Th disequilibria

Uranium-series isotope ratios determined for 35 volcanic rocks and 4 glass separates erupted from ~36 to 4.8 ka at Mt. Mazama, Crater Lake, Oregon, identify both ²³⁰Th-excess and ²³⁸U-excess components. U–Th isotope compositions cover a wide range, exceeding those previously measured for the Cascade arc. Age-corrected (²³⁰Th/²³²Th) and (²³⁸U/²³²Th) activity ratios range from 1.113 to 1.464 and from 0.878 to 1.572 (44.4 % ²³⁰Th-excess to 8.8 % ²³⁸U-excess), respectively. The most distinctive aspect of the data set is the contrast in U–Th isotope ratios between low and high Sr (LSr, HSr) components that have been previously identified in products of the 7.7 ka caldera-forming climactic eruption and preclimactic rhyodacite lavas. The LSr component exclusively contains ²³⁸U-excess, but the HSr component, as well as more primitive lavas, are marked by ²³⁰Th-excess. ²³⁰Th-excesses such as those recorded at Mt. Mazama are commonly observed in the Cascades. Melting models suggest that high ²³⁰Th-excesses observed in the more primitive lavas evolved through mixing of a mantle melt with a partial melt of a mafic lower crustal composition that contained garnet in the residuum that was produced through dehydration melting of amphibolite that was initially garnet free. Dehydration melting in the lower crust offers a solution to the “hot-slab paradox” of the Cascades, where low volatile contents are predicted due to high slab temperatures, yet higher water contents than expected have been documented in erupted lavas. The ²³⁸U-excess observed at Mt. Mazama is rare in Cascade lavas, but occurs in more than half of the samples analyzed in this study. Traditionally, ²³⁸U-excess in arc magmas is interpreted to reflect slab fluid fluxing. Indeed, ²³⁸U-excess in arcs is common and likely masks ²³⁰Th-excess resulting from lower crustal interaction. Isotopic and trace element data, however, suggest a relatively minor role for slab fluid fluxing in the Cascades. We propose that ²³⁸U-excess reflects melting and assimilation of young, hydrothermally altered upper crust. The processes related to generating ²³⁸U-excess are likely important features at Mt. Mazama that accompanied development of a large-scale silicic magma chamber that led to the caldera-forming eruption.     

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.     

U-Th-Pb and Lu-Hf isotopic constraints on the evolution of sub-continental lithospheric mantle, French Massif Central

We have carried out a Pb double-spike and Lu-Hf isotope study of clinopyroxenes from spinel-facies mantle xenoliths entrained in Cenozoic intraplate continental volcanism of the French Massif Central (FMC). U-Th-Pb and Lu-Hf isotope systematics verify the existence of different lithospheric domains beneath the northern and southern FMC. Northern FMC clinopyroxenes have extreme Lu/Hf ratios and ultra-radiogenic Hf (ε_Hf = +39.6 to +2586) that reflect ∼15-25% partial melting in Variscan times (depleted mantle model ages ∼360 Ma). Zr, Hf and Th abundances in these clinopyroxenes are low and unaffected by hydrous/carbonatitic metasomatism that overprinted LILE and light REE abundances and caused decoupling of Lu/Hf-Sm/Nd ratios and Nd-Hf isotopes (ε_Nd = +2.1 to +91.2). Pb isotopes of northern FMC clinopyroxenes are radiogenic (²⁰⁶Pb/²⁰⁴Pb > 19), and typically more so than the host intraplate volcanic rocks. ²³⁸U/²⁰⁴Pb ratios range from 17 to 68, and most samples have distinctively low ²³²Th/²³⁸U (<1) and ²³²Th/²⁰⁴Pb (3-22). Clinopyroxenes from southern FMC lherzolites are generally marked by overall incompatible trace element enrichment including Zr, Hf and Th abundances, and have Pb isotopes that are similar to or less radiogenic than the host volcanic rocks. Hf isotope ratios are less radiogenic (ε_Hf = +5.4 to +41.5) than northern FMC mantle and have been overprinted by silicate-melt-dominated metasomatism that affected this part of FMC mantle. Major element and Lu concentrations of clinopyroxenes from southern FMC harzburgites are broadly similar to northern FMC clinopyroxenes and suggest they experienced similar degrees of melt extraction as northern FMC mantle. ²³⁸U/²⁰⁴Pb (53-111) and ²³²Th/²⁰⁴Pb ratios (157-355) of enriched clinopyroxenes from the southern FMC are extreme and significantly higher than the intraplate volcanic rocks. In summary, mantle peridotites from different parts of the FMC record depletion at ∼360 Ma during Variscan subduction, followed by differing styles of enrichment. Northern FMC mantle was overprinted by a fluid/carbonatitic metasomatic agent that carried elements like U, Pb, Sr and light REE. In contrast, much of the southern FMC mantle was metasomatised by a small-degree partial silicate melt resulting in enrichment of all incompatible trace elements. The extreme mantle ²³⁸U/²⁰⁴Pb (northern and southern FMC), ²³²Th/²³⁸U (northern FMC) and ²³²Th/²⁰⁴Pb ratios (southern FMC), coupled with unremarkable present-day Pb isotope ratios, constrain the timing of enrichment. Mantle metasomatism is a young feature related to melting of the upwelling mantle responsible for Cenozoic FMC volcanism, rather than subduction-related metasomatism intimately associated with mantle depletion during the Variscan orogeny. The varying metasomatic styles relate to pre-existing variations in the thickness of the continental lithospheric lid, which controlled the extent to which upwelling mantle could ascend and melt. In the northern FMC, a thicker and more refractory lithospheric lid (?80 km) only allowed incipient degrees of melting resulting in fluid/carbonatitic metasomatism of the overlying sub-continental lithospheric mantle. The thinner lithospheric lid of the southern FMC (?70 km) allowed larger degrees of melting and resulted in silicate-melt-dominated metasomatism, and also focused the location of the volcanic fields of the FMC above this region.     

(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.     

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.     

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.     

High U/Th partitioning by clinopyroxene from alkali silicate and carbonatite metasomatism: an origin for Th/U disequilibrium in mantle melts?

Clinopyroxene/melt pairs in strongly potassic silicate and carbonatite melts exhibit unusually high U/Th partitioning ratios of ˜ 3 and ˜ 2, respectively. These values are much higher than those found for aluminous clinopyroxenes in peridotite, and have the potential to cause significant (²³⁰Th)/(²³⁸U) isotope enrichment in volcanics. The potassic silicate (lamproite) and carbonatite melts correspond closely to the main agents of mantle metasomatism, indicating that clinopyroxene in metasomatized regions of the mantle may greatly affect U/Th disequilibria. Recycling of alkali pyroxenite veins in the oceanic lithosphere formed by solidification of melt in the extremities of the MORB melting region presents an alternative to eclogite recycling in MORB and OIB genesis.     

Dynamics of melting beneath a small-scale basaltic system: a U-Th–Ra study from Rangitoto volcano,Auckland volcanic field,New Zealand

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.     

Determining melt productivity of mantle sources from U-Th and U-Pa disequilibria; an example from Pico Island, Azores

We use coupled ²³⁸U-²³⁰Th and ²³⁵U-²³¹Pa disequilibria measurements from Pico Island, Azores to examine the melting behavior of the underlying mantle. U-series disequilibria in young, mafic lavas are dependent on the melting rate of their source, which in most cases is primarily controlled by its melt productivity. Mafic lithologies such as eclogite and pyroxenite have much higher melt productivities than peridotite and so U-series measurements may provide constraints on the mineralogy of the melting mantle. Recent Pico Mountain lavas show limited geochemical variations and a restricted range of U-series disequilibria with (²³⁰Th/²³⁸U) = 1.22-1.25 and (²³¹Pa/²³⁵U) = 1.46-1.50. Using a simple, dynamic melting model of a homogeneous source, these results can be reproduced with melting rates of <1 × 10⁻⁴ kg/m³/a and melt porosities of <0.7% near the onset of melting. For a plausible range of upwelling rates, this implies that the melt productivity is <6%/GPa. This value is consistent with a garnet peridotite source but not with more highly productive mafic lithologies. Given independent evidence for the involvement of mafic lithologies such as recycled oceanic crust in Pico magmagenesis, we suggest a scenario in which initial eclogitic melts are dispersed through melt-rock reaction into a larger volume of surrounding peridotite. Subsequent re-melting of the resultant incompatible element enriched peridotite carries a geochemical signature of the mafic lithologies but not necessarily a record of their high melt productivity.     

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.     

Geochemistry and Sr–Nd–Pb isotopic composition of the Harrat Al-Madinah Volcanic Field,Saudi Arabia

Whole-rock geochemical and Sr, Nd and Pb isotope data are presented for the Harrat Al-Madinah volcanic field, in the north western part of the Arabian plate, aiming to understand their origin and the composition of their mantle source. This area is an active volcanic field characterized by the occurrence of two historic eruptions approximately in 641 and 1256 A.D. Field investigation of the main volcanic landforms indicates dominantly monogenetic strombolian eruptions, in addition to local phreatomagmatic eruption style. The lavas consist mainly of alkali olivine basalt, olivine transitional basalt, and hawaiite with ocean island basalt (OIB)-like characteristics. Evolved rocks, represented by mugearites, benmoreites, and trachytes, occur mainly as domes, tuff cones and occasionally as lava flows. Chemical variations in the evolved rocks indicated their evolution by low pressure crystal fractionation of olivine, plagioclase, clinopyroxene, and Fe–Ti oxides from the relatively primitive basalts. The isotopic compositions of ¹⁴³Nd/¹⁴⁴Nd (0.512954–0.512995), ⁸⁷Sr/⁸⁶Sr (0.702899 to–0.702977) and Pb (²⁰⁶Pb/²⁰⁴Pb = 18.5515–18.7446, ²⁰⁷Pb/²⁰⁴Pb = 15.5120–15.5222, ²⁰⁸Pb/²⁰⁴Pb = 38.1347–38.4468), show restricted variations suggesting only minor crustal contamination. They defined an array consistent with mixing of two geochemically distinct components of depleted MORB-mantle (DMM) and high ²³⁸U/²⁰⁴Pb ratio (HIMU). The variations in Tb/Yb, La/Yb and Sm/Yb ratios in the relatively primitive basalts (MgO > 6 wt.%) indicated garnet peridotite source. However, the positive Nb, Sr, Ba and Ti anomalies in the primitive mantle-normalized incompatible element patterns and the significant variation between Zr/Nb vs. Ce/Y and La/Yb vs. Yb suggest contribution of an amphibole-bearing spinel lherzolite source. Moreover, the negative correlations between SiO₂ vs. ⁸⁷Sr/⁸⁶Sr and Th vs. ¹⁴³Nd/¹⁴⁴Nd are interpreted as an indication of mixing melts derived from two end-members; one is garnet bearing asthenospheric source with OIB characteristic and the other is amphibole-bearing spinel lherzolite. The Harrat Al-Madinah volcanic field occurs near the Red Sea Rift System and its origin reflects a strong lithospheric control on the loci of partial melting. The dominantly NNW alignment patterns of the volcanoes, which is similar to the regional Red Sea trend, may suggest that the magmas were produced by decompression partial melting triggered by lithospheric extension related to the Red Rift.     

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.     

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.     

Consequences of diffuse and channelled porous melt migration on uranium series disequilibria

Magmas erupted at mid-ocean ridges (MORB) result from decompression melting of upwelling mantle. However, the mechanism of melt transport from the source region to the surface is poorly understood. It is debated whether melt is transported through melt-filled conduits or cracks on short time scales (<∼ 10³ yrs), or whether there is a significant component of slow, equilibrium porous flow on much longer time scales (>∼ 10³-10⁴ yrs). Radiogenic excess ²²⁶Ra in MORB indicates that melt is transported from the melting region on time scales less than the half life of ²²⁶Ra (∼1600 yrs), and has been used to argue for fast melt transport from the base of the melting column. However, excess ²²⁶Ra can be generated at the bottom of the melt column, during the onset of melting, and at the top of the melt column by reactive porous flow. Determining the depth at which ²²⁶Ra is generated is critical to interpreting the rate and mechanism of magma migration. A recent compilation of high quality U-series isotope data show that in many young basalts, ²²⁶Ra excess in MORB is negatively correlated with ²³⁰Th excess. The data suggest that ²²⁶Ra excess is generated independently of ²³⁰Th excess, and cannot be explained by “dynamic” or fractional melting, where observed radiogenic excesses are all generated at the base of the melt column. One explanation is that the negative correlation of activity ratios is a result of mixing of slow moving melt that has travelled through reactive, low-porosity pathways and relatively fast moving melt that has been transported in unreactive high-porosity channels. We investigate this possibility by calculating U-series disequilibria in a melting column in which high-porosity, unreactive channels form within a low-porosity matrix that is undergoing melting. The results show that the negative correlation of ²²⁶Ra and ²³⁰Th excesses observed in MORB can be produced if ∼60% of the total melt flux travels through the low-porosity matrix. This melt maintains ²²⁶Ra excesses via chromatographic fractionation of Ra and Th during equilibrium transport. Melt that travels through the unreactive, high-porosity channels is not able to maintain significant ²²⁶Ra excesses because Ra and Th are not fractionated from each other during transport and the transport time for melt in the channels to reach the top of the melt column is longer than the time scale for ²²⁶Ra excesses to decay. Mixing of melt from the high porosity channels with melt from the low-porosity matrix at the top of the melting column can produce a negative correlation of ²²⁶Ra and ²³⁰Th excesses with the slope and magnitude observed in MORB. This transport process can also account for other aspects of the geochemistry of MORB, such as correlations between La/Yb, α_Sm/Nd, and Th/U and ²²⁶Ra and ²³⁰Th excess.     

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.     

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.     

Late Precambrian OIB magmatism in the Kuznetsk Alatau,Siberia: Geochemical features of the Kulbyurstyug Formation volcanics

Basaltoids of the Vendian–Cambrian Kulbyurstyug volcanic complex in the eastern part of Kuznetsk Alatau have a high content of titanium (TiO₂ 3–4 wt %). They are relatively enriched with LILE and HFSE (Ba 360–900, Zr 160–726, Nb 66–101, ΣREE up to 225–329 ppm), and demonstrate a fractionated spectrum of REEs (La/Yb ～ 13–17), high Nb/U (44–66), and low Th/Ta and Th/Nb. These features are comparable with derivatives of OIB (Ocean Island Basalts) magmatism. The rock geochemistry suggests the possible formation of an initial mafic melt in the garnet-bearing peridotite mantle with 2–4% degree of melting and the presence of a small amount of spinel. By the age and composition, the studied volcanics are correlated with the OIB magmatism products, occurring in the adjacent Gorny Altai and also considered to be derivatives of intraplate magmatism of Riphean–Early Cambrian age. This magmatism was caused by plume activity during formation of the Paleo-Asian Ocean.