226Ra Excesses Research Articles

Timescales of magma ascent and degassing and the role of crustal assimilation at Merapi volcano (2006–2010), Indonesia: Constraints from uranium-series and radiogenic isotopic compositions

We present new 238U-230Th-226Ra-210Pb-210Po, 87Sr/86Sr and 143Nd/144Nd isotopic data of whole-rock samples and plagioclase separates from volcanic deposits of the 2006 and 2010 eruptions at Merapi volcano, Java, Indonesia. These data are combined with available eruption monitoring, petrographic, mineralogical and Pb isotopic data to assess current theories on the cause of a recent transition from effusive dome-building (2006) to explosive (2010) activity at the volcano, as well as to further investigate the petrogenetic components involved in magma genesis and evolution. Despite the significant difference in eruption style, the 2006 and 2010 volcanic rocks show no significant difference in (238U/232Th), (230Th/232Th) and (226Ra/230Th) activity ratios, with all samples displaying U and Ra excesses. The 226Ra and 210Pb excesses observed in plagioclase separates from the 2006 and 2010 eruptions indicate that a proportion of the plagioclase grew within the decades preceding eruption. The 2006 and 2010 samples were depleted in 210Po relative to 210Pb ((210Po/210Pb)i < 1) at the time of eruption but were variably degassed (69%–100%), with the degree of 210Pb degassing strongly related to sample texture and eruption phase. In good agreement with several activity monitoring parameters, 210Po ingrowth calculations suggest that initial intrusion into the shallow magma plumbing system occurred several weeks to a few months prior to the initial 2010 eruption. The 2006 and 2010 samples show a wide range in (210Pb/226Ra) activity ratio within a single eruption at Merapi and are largely characterised by 210Pb deficits ((210Pb/226Ra) < 1). Assuming a model of complete radon degassing, the 210Pb deficits in the 2006 volcanic rocks indicate relatively longer degassing timescales of ∼2–4 years than those given by the 2010 samples of ∼0–3 years. The uranium-series and radiogenic isotopic data do not support greater crustal assimilation of carbonate material as the explanation for the more explosive behaviour of Merapi in 2010 (as has been previously suggested) and instead indicate that relatively rapid ascent of a more undegassed magma was the primary difference responsible for the transition in explosive behaviour. This interpretation is in good agreement with gas monitoring data, previous petrological studies (mineral, microlite and melt inclusion work) and maximum calculated timescale estimates using Fe-Mg compositional gradients in clinopyroxene, that also suggest more rapid movement of relatively undegassed magma in 2010 relative to 2006.

U-series disequilibria of trachyandesites from minor volcanic centers in the Central Andes

Young trachyandesite lavas from minor volcanic centers in the Central Andes record the magma differentiation processes at the base of the lower continental crust. Here we report U-series disequilibrium data for the historical lavas from the Andagua Valley in Southern Peru to define the time-scale and processes of magmatism from melting in the mantle wedge to differentiation in the crust. The Andagua lavas show (230Th)/(238U), (231Pa)/(235U), and (226Ra)/(230Th) above unity except for one more evolved lava with 230Th depletion likely owing to fractional crystallization of accessory minerals. The 226Ra excess indicates that the time elapsed since magma emplacement and differentiation in the deep crust is within 8000years. Based on the correlations of U-series disequilibria with SiO2 content and ratios of incompatible elements, we argue that the Andagua lavas were produced by mixing of fresh mantle-derived magma with felsic melt of earlier emplaced basalts in the deep crust. Because of the lack of sediment in the Chile-Peru trench, there is no direct link of recycled slabs with 230Th and 231Pa excesses in the Andagua lavas. Instead, 230Th and 231Pa excesses are better explained by in-growth melting in the upper mantle followed by magma differentiation in the crust. Such processes also produced the 226Ra excess and the positive correlations among (226Ra)/(230Th), Sr/Th, and Ba/Th in the Andagua lavas. The time-scale of mantle wedge melting should be close to the half-life of 231Pa (ca. 33ka), while it takes less than a few thousand years for magma differentiation to form intermediate volcanic rocks at a convergent margin.

210Pb-226Ra disequilibria in young gas-laden magmas

We present new 238U-230Th-226Ra-210Pb and supporting data for young lavas from southwest Pacific island arcs, Eyjafjallajökull, Iceland, and Terceira, Azores. The arc lavas have significant 238U and 226Ra excesses, whereas those from the ocean islands have moderate 230Th and 226Ra excesses, reflecting mantle melting in the presence of a water-rich fluid in the former and mantle melting by decompression in the latter. Differentiation to erupted compositions in both settings appears to have taken no longer than a few millennia. Variations in the (210Pb/226Ra)0 values in all settings largely result from degassing processes rather than mineral-melt partitioning. Like most other ocean island basalts, the Terceira basalt has a 210Pb deficit, which we attribute to ~8.5 years of steady 222Rn loss to a CO2-rich volatile phase while it traversed the crust. Lavas erupted from water-laden magma systems, including those investigated here, commonly have near equilibrium (210Pb/226Ra)0 values. Maintaining these equilibrium values requires minimal persistent loss or accumulation of 222Rn in a gas phase. We infer that degassing during decompression of water-saturated magmas either causes these magmas to crystallize and stall in reservoirs where they reside under conditions of near stasis, or to quickly rise towards the surface and erupt.

Localised magmatic constraints on continental back-arc volcanism in southern Mendoza, Argentina: the Santa Maria Volcano

The Payun Matru Volcanic Field constitutes part of the continental back-arc in Argentina. This volcanic field has been the focus of several regional investigations; however, geochemical analysis of recent volcanoes (<8 ka) at the scale of an individual volcano has not been conducted. We present a morphological description for the Santa Maria Volcano in addition to results from major and trace element analysis and 238U-230Th-226Ra disequilibria. The trace element evidence suggests that the Santa Maria magmatic source has a composition similar to that of the local intraplate end member (resembling an ocean island basalt-like source), with a slight contribution from subduction-related material. The U-series analyses suggest a high 226Ra excess over 230Th for this volcano, which is not derived from a shallow process such as hydrothermal alteration or upper crustal contamination. Furthermore, intermediate-depth processes such as fractional crystallisation have been inferred for the Santa Maria Volcano, but they are not capable of producing the 226Ra excess measured. The 226Ra excess is explained by deep processes like partial melting of mantle lithologies with some influence from subducted Chilean trench sediments. Due to the short half-life of 226Ra (1600 years), we infer that fast magma ascent rates are required to preserve the high 226Ra excess.

U-series disequilibria in subduction zone lavas: Inherited from subducted slabs or produced by mantle in-growth melting?

U-series disequilibria of subduction zone lavas have been used to provide temporal constraints on flux transfer in convergent margins. There are debates on whether the disequilibria are inherited from the source modified by slab materials (fluids or sediments) or produced by in-growth melting of the mantle. After filtering out the effect of crustal modification on U-series data, (226Ra)/(230Th) is not positively correlated with Sr/Th and Ba/Th in arc basalts with SiO2 <52wt%. This invalidates the key evidence for addition of slab-fluids to the mantle producing 226Ra excess in subduction zone lavas. Sediment melts with a zircon-free residue may have 231Pa excess, but it is not consistent with the overwhelming Zr-Hf depletions in subduction zone lavas, which instead require a zircon-saturated residue during sediment melting.With appropriate parameters applied, in-growth melting in the mantle followed by modification in the crust can provide a self-consistent explanation of the U-series disequilibria and the relationships with other geochemical signatures in most subduction zone lavas except the few Tonga-Kermadec lavas with slight 231Pa deficit. On the basis of melting experiments of hydrous peridotites, we further propose that most mantle-derived melts and thus U-series disequilibria should be mainly produced by melting in the hot zone of the mantle wedge, as controlled by thermal structure of the convergent margins. In order to produce enough short-lived nuclides, the timescales of mass transfer and magmatism in the subduction zone should at least range from a few to a few hundred millennia. Future studies on more lavas, numerical modelling, and partition coefficients of U-series nuclides will help for better application of U-series disequilibria into subduction zone magmatism.

Melt transport rates in heterogeneous mantle beneath mid-ocean ridges

Recent insights to melt migration beneath ridges suggest that channelized flow is a consequence of melting of a heterogeneous mantle, and that spreading rate modulates the dynamics of the localised flow. A corollary of this finding is that both mantle heterogeneity and spreading rate have implications for the speed and time scale of melt migration. Here, we investigate these implications using numerical models of magma flow in heterogeneous mantle beneath spreading plates. The models predict that a broad distribution of magma flow speeds is characteristic of the sub-ridge mantle. Within the melting region, magmatic flow is fastest in regions of average fusibility; surprisingly, magmas from sources of above-average fusibility travel to the ridge in a longer time. Spreading rate has comparatively simple consequences, mainly resulting in faster segregation speeds at higher spreading rates. The computed time scales are short enough to preserve deep origin 230Th disequilibria and, under favourable parameter regimes, also 226Ra excesses. An important prediction from the models is that mantle heterogeneity induces significant natural variability into U-series disequilibria, complicating the identification of relationships between disequilibria and ridge properties or chemical signatures of heterogeneity.

238U–230Th–226Ra Disequilibria Constraints on the Magmatic Evolution of the Cumbre Vieja Volcanics on La Palma, Canary Islands

A suite of 48 samples, including both historical and prehistoric lavas and some plutonic rocks, have been analysed from the Cumbre Vieja rift, La Palma, Canary Islands. Additionally, mineral–melt partition coefficients have been measured for clinopyroxene, plagioclase, amphibole, titanite and apatite in selected rocks. The lavas range from basanite to phonolite (SiO2 = 41·2–57·5 wt % and MgO = 10–0·8 wt %) in composition and form coherent, curvilinear major and trace element arrays in variation diagrams, irrespective of eruption age. The mafic lavas have typical ocean island incompatible trace element patterns and Sr, Nd and Pb isotope compositions show little variation but have a HIMU-type character. Generation of the parental magmas is inferred to have involved ∼4% dynamic melting of a garnet lherzolite source that may have previously been metasomatized by melts derived from a recycled mafic component containing residual phlogopite. The major process of differentiation to phonotephrite involved fractional crystallization of basanitic magmas that evolved along the same liquid line of descent under similar pressure–temperature conditions. Numerical simulations using the MELTS algorithm suggest that this occurred across a temperature interval from c. 1320 to 950°C at 400 MPa and an oxygen fugacity equivalent to quartz–fayalite–magnetite (QFM), with an initial H2O content of 0·3 wt %. The later stages of differentiation (<5 wt % MgO) were dominated by mixing with partial melts of young syenites formed from earlier magma batches. All of the lavas are characterized by 230Th and 226Ra excesses and (230Th/238U) decreases with decreasing Nb/U and increasing SiO2, with no accompanying change in (226Ra/230Th). To explain the observations, we propose a model in which there was a significant role for amphibole, and more importantly accessory titanite, in decre'asing Nb/U, Ce/Pb and Th/U ratios and increasing or buffering (226Ra/230Th) ratios during the later stages of differentiation and magma mixing. These processes all occurred over a few millennia in small magma batches that were repeatedly emplaced within the mid-crust of the Cumbre Vieja rift system prior to rapid transport to the surface.

Melting versus contamination effects on 238U–230Th–226Ra and 235U–231Pa disequilibria in lavas from São Miguel, Azores

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

238U–230Th equilibrium in arc magmas and implications for the time scales of mantle metasomatism

Large excesses of 238U and 226Ra relative to 230Th characterize many arc magmas and are commonly interpreted to represent recent addition of slab-derived fluid to the mantle wedge beneath the arc. A significant proportion of arc magmas are, however, in 238U–230Th radioactive equilibrium. This is generally thought to result from “buffering” of the young slab fluid U-series signal by a sediment component in secular equilibrium. Here we present new 238U–230Th–226Ra and 235U–231Pa measurements for historic andesites from Volcán de Colima, Mexico. In all lavas (230Th/238U) are in equilibrium, whereas (231Pa/235U) and (226Ra/230Th) are significantly greater than one. These data demonstrate that arc magmas with (230Th/238U) equilibrium can have significant 231Pa and 226Ra excesses, precluding ageing of the magmas in the crust as the cause of 230Th/238U equilibrium. Quantitative modeling of metasomatic and melting processes further indicates that addition of sediment melts to a depleted mantle wedge produces significant 230Th excesses and that 238U excesses induced by recent addition of fluids derived from the altered oceanic crust are not sufficient to compensate these 230Th excesses. U-series activity ratios in Colima magmas are best explained by models in which the metasomatised mantle returns to secular equilibrium before melting, implying a time lag ⩾350 kyr, with subsequent production of 231Pa and 226Ra excesses by in-growth during melting rather than by addition of slab fluids. Investigation of a global compilation of U-series data in arc magma indicates that our model proposed for Colima applies to most arc magmas in or near (230Th–238U) equilibrium. The time lag between mantle metasomatism and melting appears to vary between hundreds years to more than 350 kyr in subduction zones. We posit that the absence of U/Th elemental fractionation during melting of arc sources in (230Th/238U) equilibrium reflects a higher fO2 compared to MORB sources that yield magmas with 230Th excesses.

Luminescence dating of scoria fall and lahar deposits from Somma–Vesuvius, Italy

Luminescence dating has been applied to scoria and lahar deposits from Somma–Vesuvius, Italy. Samples include scoria from the AD472 and 512 (or 536) eruptions and lahar deposits. In order to find a stable luminescence signal which is less affected by anomalous fading, infrared stimulated luminescence (IRSL) signals at elevated temperatures after bleaching with IR at 50 °C (termed post-IR IRSL; pIRIR) were tested at different preheat and elevated stimulation temperatures. The fading rates of both IRSL and pIRIR signals reduced dramatically with increasing preheat and pIRIR stimulation temperatures. A pIRIR signal measured at 290 °C after a preheat at 320 °C (60 s) and an IR stimulation at 50 °C (100 s) was selected to calculate the equivalent dose (De). The gamma spectrometry results indicate that the U-series nuclides are not in equilibrium and there is a large 226Ra excess. The dose rates and ages were calculated by assuming a 226Ra excess (over its parent 230Th) at deposition, and that this unsupported excess then decayed to the present level. The resulting luminescence ages of the two scoria samples agreed with the expected ages, and the ages of the lahar deposits indicate that they are associated with the AD1631 eruption.

Mantle flow, volatiles, slab‐surface temperatures and melting dynamics in the north Tonga arc–Lau back‐arc basin

The Fonualei Spreading Center affords an excellent opportunity to evaluate geochemical changes with increasing depth to the slab in the Lau back‐arc basin. We present H2O and CO2concentrations and Sr, Nd, Pb, Hf and U‐Th‐Ra isotope data for selected glasses as well as new Hf isotope data from boninites and seamounts to the north of the Tonga arc. The Pb and Hf isotope data are used to show that mantle flow is oriented to the southwest and that the tear in the northern end of the slab may not extend east as far as the boninite locality. Along the Fonualei Spreading Center, key geochemical parameters change smoothly with increasing distance from the arc front and increasing slab surface temperatures. The latter may range from 720 to 866°C, based on decreasing H2O/Ce ratios. Consistent with experimental data, the geochemical trends are interpreted to reflect changes in the amount and composition of wet pelite melts or super‐critical fluids and aqueous fluids derived from the slab. With one exception, all of the lavas preserve both238U excesses and 226Ra excesses. We suggest that lavas from the Fonualei Spreading Center and Valu Fa Ridge are dominated by fluid‐fluxed melting whereas those from the East and Central Lau Spreading Centers, where slab surface temperatures exceed ∼850–900°C, are largely derived through decompression. A similar observation is found for the Manus and East Scotia back‐arc basins and may reflect the expiry of a key phase such as lawsonite in the subducted basaltic crust.

Rapid magmatic processes accompany arc–continent collision: the Western Bismarck arc, Papua New Guinea

New U–Th–Ra, major and trace element, and Sr–Nd–Pb isotope data are presented for young lavas from the New Britain and Western Bismarck arcs in Papua New Guinea. New Britain is an oceanic arc, whereas the latter is the site of an arc–continent collision. Building on a recent study of the Manus Basin, contrasts between the two arcs are used to evaluate the processes and timescales of magma generation accompanying arc–continent collision and possible slab detachment. All three suites share many attributes characteristic of arc lavas that can be ascribed to the addition of a regionally uniform subduction component derived from the subducting altered oceanic crust and sediment followed by dynamic melting of the modified mantle. However, the Western Bismarck arc lavas diverge from the Pb isotope mixing array formed by the New Britain and the Manus Basin lavas toward elevated 208Pb/204Pb. We interpret this to reflect a second and subsequent addition of sediment melt at crustal depth during collision. 238U and 226Ra excesses are preserved in all of the lavas and are greatest in the Western Bismarck arc. High-Mg andesites with high Sr/Y ratios in the westernmost arc are attributed to recent shallow mantle flux melting at the slab edge. Data for two historical rhyolites are also presented. Although these rhyolites formed in quite different tectonic settings and display different geochemical and isotopic compositions, both formed from mafic parents within millennia.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52·5–57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene ± orthopyroxene with microlites of Ti-magnetite. Olivine (Fo83–88) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950–1200°C and 0·8–1·8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4·16 wt % H2O, and CO2–H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4–5·5 km and 1·5–2·5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0·2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U–Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (226Ra/230Th) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Insights into the Galápagos plume from uranium-series isotopes of recently erupted basalts

Uranium-series isotopes (238U-230Th-226Ra-210Pb), major element, trace element, and Sr-Nd isotopic data are presented for recent (<60 years old) Galapagos archipelago basalts. Volcanic rocks from all centers studied (Fernandina, Cerro Azul, Sierra Negra, and Wolf Volcano) display 230Th excesses (4%–15%) and steep rare earth element (REE) patterns indicative of residual garnet during partial melting of their mantle source. Rare earth element modeling suggests that only a few percent of garnet is involved. Correlations between (238U/232Th), radiogenic isotopes and Nb/Zr ratio suggest that the U/Th ratio of these Galapagos volcanic rocks is primarily controlled by geochemical source variations and not fractionation during partial melting. The lowest (230Th/238U) ratio is not observed at Fernandina (the supposed center of the plume) but at the more geochemically “depleted” Wolf Volcano, further to the north. Small radium excesses are observed for all samples with (226Ra/230Th) ranging from 1.107 to 1.614. The 226Ra-230Th disequilibria do not correlate with other uranium-series parent-daughter nuclide pairs or geochemical data, suggesting modification at shallow levels on timescales relevant to the half-life of 226Ra (1600 years). The combination of 226Ra and 210Pb excesses is inconsistent with interaction of magma with cumulate material unless decoupling of 210Pb (or an intermediate daughter, such as 222Rn) occurs prior to modification of Ra-Th disequilibria. An intriguing correlation of (210Pb/226Ra)0 with Nb/Zr and radiogenic isotopes requires further investigation but suggests possible control via magmatic degassing and accumulation that may somehow be related to source heterogeneities.

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 230Th, 226Ra, 231Pa, and 210Po 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 ( 226Ra)/( 210Po) within ±15% of unity. Ten out of 14 samples have significant 226Ra excesses over 230Th, with ( 226Ra)/( 230Th) up to 2.97, while four samples are in 226Ra– 230Th equilibrium within error. All KEJ samples have high ( 231Pa)/( 235U), ranging from 1.56 to 2.64 and high 238U excesses (up to 43%), providing a global end-member of high 238U and high 231Pa 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 ( 231Pa)/( 235U) does not correlate with differentiation indices such as Ho/Sm, ( 230Th)/( 238U) shows a slight negative correlation, likely due to assimilation of materials with slightly higher ( 230Th)/( 238U). Samples with 226Ra excess have higher Sr/Th and Ba/Th than those in 226Ra– 230Th equilibrium, forming rough positive correlations of ( 226Ra)/( 230Th) 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 231Pa excesses require an in-growth melting process operating at low melting rates and small residual porosity; such a model will also produce significant 238U– 230Th and 226Ra– 230Th 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 238U– 230Th, 226Ra– 230Th, and 235U– 231Pa disequilibria in erupted KEJ lavas reflect the in-growth melting process in the mantle wedge (reflecting variations in U/Th, daughter–parent ratios, fO 2, 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 226Ra excess is derived from the slab.

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 238U– 230Th– 226Ra 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 ( 87Sr/ 86Sr = 0.70319–0.70322; 143Nd/ 144Nd = 0.51302–0.51305; 206Pb/ 204Pb = 19.04–19.06; 207Pb/ 204Pb = 15.53–15.54; 208Pb/ 204Pb = 38.61–38.65; 176Hf/ 177Hf = 0.28311–0.28312). All the samples possess ( 230Th/ 238U) disequilibrium with 230Th excesses, and they show systematic variations in ( 230Th/ 232Th) and ( 238U/ 232Th) ratios. The highest 226Ra excesses occur in the basalt and most differentiated andesite lavas, while some basaltic-andesite lavas have ( 226Ra/ 230Th) ratio that are close to equilibrium. The 238U– 230Th– 226Ra disequilibria variations cannot be produced by simple closed-system fractional crystallization with radioactive decay of 230Th and 226Ra 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 230Th– 238U and 226Ra– 230Th 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.

Influence of subducted components on back‐arc melting dynamics in the Manus Basin

Lavas erupted in back‐arc basins afford the opportunity to explore the extent to which decompression and subduction‐related components influence partial melting in this setting. We present U‐Th‐Ra disequilibria data from 24 well‐characterized lavas from the Manus Basin behind the New Britain volcanic arc, supplemented by some additional trace element and Sr‐Nd‐Pb isotope data. The lavas range in composition from 49.6 to 57.7 wt % SiO2 and can be subdivided into those that are broadly like mid‐ocean ridge basalts (MORB) with Ba/Nb &lt; 16 and back‐arc basin basalts (BABB) that are variably influenced by subduction components and have Ba/Nb &gt; 16. Rifts closest to the arc are dominated by BABB, whereas both lavas types erupt further away at the Manus Spreading Center. The MORB have small 230Th excesses (up to 5%) and are displaced below the global correlation of (230Th/238U) with ridge depth. In most respects the BABB closely resemble lavas erupted along the New Britain arc front, including 238U excesses that reach 26%. The Pb isotope data can be explained by mixing of a subduction component into an Indian MORB mantle source. The Pb in the subduction component is derived from both the subducted sediment (5%) and fluids from the subducting altered Solomon Sea oceanic crust (95%), and these were mixed prior to addition to the mantle wedge. U/Th ratios, Fe3+/ΣFe, and H2O contents all increase with increasing 206Pb/204Pb. A model in which addition of the subduction component to the mantle wedge is followed by 230Th in‐growth during decompression and dynamic melting all less than 140 kyr prior to eruption can simulate the data. However, our preferred model is one of dynamic decompression melting in which subduction‐modified, more oxidized mantle had DU ≪ DTh leading to 238U excesses in contrast to unmodified mantle that yields 230Th excess. Large 226Ra excesses in some southern rift samples require addition of a fluid &lt;8 kyr ago but elsewhere reflect melting under low‐porosity conditions.

Generation and evolution of magma beneath the East Pacific Rise: Constraints from U-series disequilibrium and plagioclase-hosted melt inclusions

Mid-ocean ridge basalt (MORB) samples from the East Pacific Rise (EPR 12°50′N) were analyzed for U-series isotopes and compositions of plagioclase-hosted melt inclusions. The 226Ra and 230Th excesses are negatively correlated; the 226Ra excess is positively correlated with Mg# and Sm/Nd, and is negatively correlated with La/Sm and Fe 8; the 230Th excess is positively correlated with Fe 8 and La/Sm and is negatively correlated with Mg# and Sm/Nd. Interpretation of these correlations is critical for understanding the magmatic process. There are two models (the dynamic model and the “two-porosity” model) for interpreting these correlations, however, some crucial parameters used in these models are not ascertained. We propose instead a model to explain the U-series isotopic compositions based on the control of melt density variation. For melting either peridotite or the “marble-cake” mantle, the FeO t content, 230Th excess and La/Sm ratio increases and Sm/Nd decreases with increasing pressure. A deep melt will evolve to a higher density and lower Mg# than a shallow melt, the former corresponds to a long residence time, which lowers the 226Ra excess significantly. This model is supported by the existence of low 226Ra excesses and high 230Th excesses in MORBs having a high Fe 8 content and high density. The positive correlation of 226Ra excess and magma liquidus temperature implies that the shallow melt is cooled less than the deep melt due to its low density and short residence time. The correlations among Fe 8, Ti 8 and Ca 8/Al 8 in plagioclase-hosted melt inclusions further prove that MORBs are formed from melts having a negative correlation in melting depths and degrees. The negative correlation of 226Ra excess vs. chemical diversity index (standard deviation of Fe 8, Ti 8 and Ca 8/Al 8) of the melt inclusions is in accordance with the influence of a density-controlled magma residence time. We conclude that the magma density variation exerts significant control on residence time and U-series isotopic compositions.

The effect of assimilation, fractional crystallization, and ageing on U-series disequilibria in subduction zone lavas

Although most arc lavas have experienced significant magma differentiation, the effect of the differentiation process on U-series disequilibria is still poorly understood. Here we present a numerical model for simulating the effect of time-dependent magma differentiation processes on U-series disequilibria in lavas from convergent margins. Our model shows that, in a closed system with fractional crystallization, the ageing effect can decrease U-series disequilibria via radioactive decay while in an open system, both ageing and bulk assimilation of old crustal material serve to reduce the primary U-series disequilibria. In contrast, with recharge of refresh magma, significant 226Ra excess in erupted lavas can be maintained even if the average residence time is longer than 8000 years. The positive correlations of ( 226Ra/ 230Th) between Sr/Th or Ba/Th in young lavas from convergent margins have been widely used as evidence of fluid addition generating the observed 226Ra excess in subduction zones. We assess to what extent the positive correlations of ( 226Ra/ 230Th) with Sr/Th and Ba/Th observed in the Tonga arc could reflect AFC process. Results of our model show that these positive correlations can be produced during time-dependent magma differentiation at shallow crustal levels. Specifically, fractional crystallization of plagioclase and amphibole coupled with contemporaneous decay of 226Ra can produce positive correlations between ( 226Ra/ 230Th) and Sr/Th or Ba/Th (to a lesser extent). Therefore, the correlations of ( 226Ra/ 230Th) with Sr/Th and Ba/Th cannot be used to unambiguously support the fluid addition model, and the strength of previous conclusions regarding recent fluid addition and ultra-fast ascent rates of arc magmas is significantly lessened.

U-series dating of Lake Nyos maar basalts, Cameroon (West Africa): Implications for potential hazards on the Lake Nyos dam

From previously published 14C and K–Ar data, the age of formation of Lake Nyos maar in Cameroon is still in dispute. Lake Nyos exploded in 1986, releasing CO 2 that killed 1750 people and over 3000 cattle. Here we report results of the first measurements of major elements, trace elements and U-series disequilibria in ten basanites/trachy-basalts and two olivine tholeiites from Lake Nyos. It is the first time tholeiites are described in Lake Nyos. But for the tholeiites which are in 238U– 230Th equilibrium, all the other samples possess 238U– 230Th disequilibrium with 15 to 28% enrichment of 230Th over 238U. The ( 226Ra/ 230Th) activity ratios of these samples indicate small (2 to 4%) but significant 226Ra excesses. U–Th systematics and evidence from oxygen isotopes of the basalts and Lake Nyos granitic quartz separates show that the U-series disequilibria in these samples are source-based and not due to crustal contamination or post-eruptive alteration. Enrichment of 230Th is strong prima facie evidence that Lake Nyos is younger than 350 ka. The 230Th– 226Ra age of Nyos samples calculated with the ( 226Ra/ 230Th) ratio for zero-age Mt. Cameroon samples is 3.7 ± 0.5 ka, although this is a lower limit as the actual age is estimated to be older than 5 ka, based on the measured mean 230Th/ 238U activity ratio. The general stability of the Lake Nyos pyroclastic dam is a cause for concern, but judging from its 230Th– 226Ra formation age, we do not think that in the absence of a big rock fall or landslide into the lake, a big earthquake or volcanic eruption close to the lake, collapse of the dam from erosion alone is as imminent and alarming as has been suggested.