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.     

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.     

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.     

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.     

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

Uranium series disequilibrium in a young surficial uranium deposit,northeastern Washington,U.S.A.

A recently discovered ore-grade accumulation of U in organic-rich sediments of late Quaternary age provides an opportunity for studying the early association of U, U-daughters, and organic matter in a natural setting. The U occurs in valley-fill sediments of peat, peaty clay, silt, and sand along the north fork of Flodelle Creek, Stevens County, Washington. Radiometric techniques (delayed neutron, high-resolution gamma-ray spectrometry, thin-source alpha spectrometry) were employed to determine the abundance and distribution of U-series nuclides, the extent of secular equilibrium within the U decay series, and the apparent U-series ages of U incorporation.Sixteen lithologically distinct intervals were sampled from a 292 cm core. Uranium contents range from 140 to 2790 ppm and are positively correlated with organic contents. Measured alpha activity ratios of ²³⁴U/²³⁸U (1.31–1.38) are very similar to those reported in coexisting waters, suggesting a rather constant isotopic composition of introduced U. Much lower Th contents of <10–40 ppm are controlled by the type and abundance of silicate detritus. The youth of the host sediments (<15 000 a) and the paucity of associated radioactivity suggested large excesses of U relative to radioactive daughters and such excesses were observed, particularly in the shallowest intervals. Apparent ages of U emplacement determined by the (alpha) activity ratio of ²³⁰Th daughter to ²³⁴U parent show a general increase with depth and fair agreement with estimated depositional ages. This observation suggests dominantly syndepositional or early post depositional emplacement of U followed by decay-generated buildup of ²³⁰Th daughter with time. However, interval by interval comparisons of the relative abundances of other daughters, particularly ²²⁶Ra and ²¹⁰Pb, indicate variability caused by processes other than closed-system growth and decay, probably because chemically diverse daughters that are decay-generated in situ have differing mobilities and because upwelling ground water continuously adds more U and minor amounts of daughters. If ²³⁰Th is considered the least susceptible to these modifications, the data suggest some addition of ²³⁴U in the deepest intervals and some loss of ²²⁶Ra and/or gain of ²²²Rn throughout the studied core.     

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.     

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.     

Extreme U-Th Disequilibrium in Rift-Related Basalts, Rhyolites and Granophyric Granite and the Timescale of Rhyolite Generation, Intrusion and Crystallization at Alid Volcanic Center, Eritrea

Rhyolite pumices and co-erupted granophyric (granite) xenolithsyield evidence for rapid magma generation and crystallizationprior to their eruption at 15·2 ± 2·9 kaat the Alid volcanic center in the Danikil Depression, Eritrea.Whole-rock U and Th isotopic analyses show ²³⁰Th excesses upto 50% in basalts <10 000 years old from the surroundingOss lava fields. The 15 ka rhyolites also have 30–40%²³⁰Th excesses. Similarity in U–Th disequilibrium, andin Sr, Nd, and Pb isotopic values, implies that the rhyolitesare mostly differentiated from the local basaltic magma. Giventhe (²³⁰Th/²³²Th) ratio of the young basalts, and presumablythe underlying mantle, the (²³⁰Th/²³²Th) ratio of the rhyolitesupon eruption could be generated by in situ decay in about 50000 years. Limited (5%) assimilation of old crust would hastenthe lowering of (²³⁰Th/²³²Th) and allow the process to takeplace in as little as 30 000 years. Final crystallization ofthe Alid granophyre occurred rapidly and at shallow depths at20–25 ka, as confirmed by analyses of mineral separatesand ion microprobe data on individual zircons. Evidently, 30000–50 000 years were required for extraction of basaltfrom its mantle source region, subsequent crystallization andmelt extraction to form silicic magmas, and final crystallizationof the shallow intrusion. The granophyre was then ejected duringeruption of the comagmatic rhyolites. KEY WORDS: U-series; zircon; ion microprobe; volcano; geochronology     

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.     

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.     

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.     

Mineralogy,geochemistry and petrogenesis of Kurile island-arc basalts

Whole-rock (major- and trace-element) and mineral chemical data are presented for basaltic rocks from the main evolutionary stages of the Kurile island arc, NW Pacific. An outer, inactive arc contains a Cretaceous-Lower Tertiary sequence of tholeiitic, calcalkaline and shoshonitic basalts. The main arc (Miocene-Quaternary) is dominated by weakly tholeiitic, with lesser, alkalic basalts. The mineralogy of Kuriles basalts is characterised by An-rich plagioclases, a continuous transition from chromites to titanomagnetites, pyroxenes with low Fe³⁺ contents and without strong Fe-enrichment, abundance of groundmass pigeonites and the absence of amphiboles. There is an increase in K₂O contents both along-arc (northwards) and towards the reararc side. The basalts show an exceptionally wide but continuous range of K₂O contents (0.1–4.7%) which correlate with other LIL element contents. Tholeiitic basalts with low LIL element contents, La/Yb and Th/U, but high K/ Rb, P₂O₅/La and Zr/Nb were derived from depleted, lherzolitic mantle which had suffered fluid metasomatism by K, Rb, Cs, Sr, Ba, Pb and H₂O only. Alkali basalts are also thought to be derived from depleted mantle but melt metasomatism involved addition of all LIL elements to a garnet lherzolite mantle. The Kuriles basalts and their mantle sources range continuously between these two end-member compositions. The metasomatic fluids/melts were probably released by early dehydration and later melting within subducted oceanic lithosphere though the process is not adequately constrained.     

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.     

Review of melting experiments on carbonated eclogite and peridotite: insights into mantle metasomatism

Experimental studies on the partial melting of eclogite and peridotite provide important clues on mantle metasomatism. Here, we review results from some of the recent experiments and show that melting of carbonated eclogite and peridotite can produce carbonatitic to carbonated silicate melt, in which carbonates melt preferentially before Ti oxides and silicates. Low-degree melting results in carbonatitic melt coexisting with Ti oxides and silicates. This process also leads to the fractionation between some high-field strength elements (Nb, Ta, Zr, Hf, and HREE) and highly incompatible elements (U and Th) in the melt. When Ti oxides are nearly exhausted in eclogite, extremely high TiO₂ contents (e.g. 19 wt.%) are present in the melt with marked concentration of Nb and Ta. These results help to explain the features of carbonatitic metasomatism and the Nb–Ta spike in oceanic island basalts as identified in experimental studies. These studies also explain the reducing conditions that stabilize diamond in the deep mantle (>150 km) as well as the occurrence of diamond at different depths reported in various studies. Melting in such a reduced mantle can happen through redox reaction between diamond, pyroxene, and olivine, in which the initial liquid is a carbonated silicate melt. However, the theoretical oxygen fugacity (fO₂) in the asthenosphere is much lower than that predicted by the reaction and requires elevated fO₂, which can be caused by the addition of relatively oxidized materials from the lower mantle, deep asthenospheric material, and various recycled components. A combination of these processes generates locally oxidized domains in the deep mantle.     

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.     

云南宾川地区峨眉山玄武岩地球化学特征:岩石类型及随时间演化规律

分布于云南宾川地区、厚逾5000m的晚二叠世峨眉山玄武岩以拉斑玄武岩为主,少量碱性玄武岩。根据它们的岩相学和主量元素、微量元素特征,将其划分为两个地球化学类型:低钛玄武岩(LT)和高钛玄武岩(HT)。低钛玄武岩主要分布于岩石剖面的中下部,其主要地球化学标志为低Ti/Y(<500),高Mg^#(048～067),低的∑REE(<120μg/g)、(La/Yb)_N(<6)和(Sm/Yb)_N(<25)。高钛玄武岩分布于剖面的顶部,以高Ti/Y(>500),低Mg^#(039～053)、高的∑REE(>150)、(La/Yb)_N(>9)和(Sm/Yb)_N(>3)为特征。根据LT的分异指数和Th、U异常特征等,将其分为LT1和LT2两个亚类。其中LT1位于岩石剖面下部,表现为明显的Th和U正异常,高Mg^#(061～067)和低Nb/U比值等。LT2位于岩石剖面中部,具Rb、Ba正异常,无Th和U异常和较低的分异指数(Mg^#=048～054)。LT玄武岩可能是峨眉地幔热柱主活动期地幔柱头部熔融的产物,后经较弱的橄榄石+单斜辉石±斜长石结晶分异形成。剖面最底部的LT1玄武岩显示较强的“壳源”印记,可能与岩石圈地幔中富集组份的活化有关。这些壳源物质的参与程度自底部向上有降低的趋势。晚期HT玄武岩为地幔热柱消亡期的产物,其在地壳浅部经历强烈的以斜长石为主的结晶分异,壳源物质的混染不明显。     

Petrogenesis of the Apollo 14 high-alumina basalts: Implications from ion microprobe analyses

In this study, ion microprobe analyses of individual minerals are used to investigate the petrogenesis of the Apollo 14 high-Al basalts. We use trace element concentrations from individual minerals in the Apollo 14 high-Al basalts to evaluate both endogenic and exogenic models. The data show that if the Apollo 14 high-Al basalts were produced by melting within the lunar mantle, these basalts cannot be related to one another by closed-system fractional crystallization of a single basaltic melt. Rather, the trace element data show that variable amounts of a KREEP component were added to the basalts by either assimilation, mixing into mantle sources, or impact melting. Single-stage assimilation-fractional crystallization models can only explain the data from this study if an excessively large mass of urKREEP is assimilated into the parent magma before olivine crystallization. Alternatively, the trace element data can be explained if the Apollo 14 high-Al basalts were produced by melting multiple Al-rich mantle sources that contain different amounts of urKREEP. Finally, for impact melting to be a relevant process, the data require that multiple large impact melts be formed from mixed KREEP-rich target lithologies. The resulting impact melts must then crystallize to produce basalts with igneous textures, high Al₂O₃ concentrations, uniform major element compositions, and a wide range of incompatible trace element concentrations.