Time-scales for magmatic differentiation at the Snaefellsjökull central volcano, western Iceland: Constraints from U-Th-Pa-Ra disequilibria in post-glacial lavas

We present major and trace element and Sr-Nd-Pb and U-Th-Pa-Ra isotope data for a small sample suite of primarily post-glacial, mildly alkalic volcanic rocks from the Snaefellsjökull central volcano situated off the main rift systems in western Iceland. The volcanic rocks are crystal-poor and range from olivine alkali basalt to trachyte and show tight correlations of major and trace elements that are explained by fractional crystallization involving removal of olivine, plagioclase, clinopyroxene, Fe-Ti oxide and apatite. Sr-Nd-Pb isotopes are practically invariant, consistent with derivation from the same source region. During fractionation from primitive basalt to evolved trachyte, (²³⁰Th/²³²Th), (²³⁰Th/²³⁸U) and (²³¹Pa/²³⁵U) decrease progressively at broadly constant (²³⁸U/²³²Th). A continuous closed-system fractionation model that assumes constant initial (²³⁰Th/²³²Th) in the basaltic precursor melt indicates that hawaiite was derived from olivine basalt by ∼50% fractional crystallization within 12.1_{- 5.2} ^{+ 5.1} kyr and trachyte by ∼80% fractionation within 53.9_{- 6.7} ^{+ 7.0} kyr. An overrepresentation of evolved basalts and hawaiites with young inferred magma ages in the dataset is consistent with the parental precursor to these magmas intruded into the sub-volcanic magma plumbing system as a consequence of lithospheric rebound caused by deglaciation. Lavas affected by apatite removal have higher (²³¹Pa/²³⁵U) than predicted for simple radioactive decay, suggesting apatite significantly fractionates U from Pa. The proposed fractionation model consistently explains our U-series data assuming D_Pa ^apa ≤ 0.5 and D_U ^apa ∼ 3.2 and D_Th ^apa ∼ 3.0. If applicable, these D values would indicate that the effect of apatite fractionation must be adequately considered when assessing differentiation time scales using (²³¹Pa/²³⁵U) disequilibria data.