Combined 238U–230Th and 235U–231Pa constraints on the transport of slab-derived material beneath the Mariana Islands

We present 235U–231Pa measurements for a suite of well-characterised lavas from the Mariana Islands. There is limited variability in (231Pa/235U), 1.18–1.55, with the lowest values typically found in the incompatible element depleted samples. We also report the complementary 238U–234U–230Th–232Th analyses, which are in good agreement previous measurements on these same samples. In keeping with this earlier work, we interpret the range of (230Th/238U) from 0.65 to 1.06 as the result of variable contributions of an enriched, sediment melt component coupled with a near constant flux of a relatively U-rich aqueous fluid, derived from the subducted mafic oceanic crust to the sub-arc mantle. We reassess the composition of the subducting sedimentary assemblage and find that its (230Th/232Th) was previously under-estimated. We performed quantitative models to reproduce the coupled 238U–230Th–232Th and 235U–231Pa data starting from a mantle source modified by subduction-related components according to traditional assumptions of a Pa (and Th)-free fluid and a sedimentary component in 231Pa–235U and 238U–230Th equilibrium. This scenario does not allow the U-series isotope compositions of the Mariana lavas to be reproduced. Successful solutions are instead possible with a sediment melt component with relatively high (230Th/232Th) 0.9 and a modest 230Th- and 231Pa excess, (30%). The need for a 230Th- and 231Pa-excess in the sedimentary component implies that the sediment melted <150 ka before eruption of the lavas that carry this signature and that it melted to sufficiently high degree to exhaust monazite. In contrast, the low (230Th/238U) of the aqueous fluid component is explained by its derivation from lower temperature, mafic crust that contained residual allanite. The maximum (230Th/232Th) of this fluid is set by the well-defined array of lavas on an equiline plot and this argues for a fresh, rather than significantly altered MORB protolith. Moreover, in-growth is insufficient to account for the 231Pa excesses observed in the Marianas lavas given a Pa-free fluid. Thus we infer that the combined U–Th–Pa signatures of the Marianas are more strongly controlled by the initial compositions of the slab components than by mantle melting. Creation of U-series disequilibria is strongly favoured by the presence of key accessory phases, such as allanite, monazite and zircon, that are likely to be stable in the slab and we suggest that the variable U–Th–Pa systematics of arc lavas globally may well reflect the contrasting stabilities of accessory phases in subduction zones with different thermal and compositional controls.