AgeSpectraAnalyst: A MATLAB based package to model zircon age distributions in silicic magmatic systems.

In the last decade, improvements in the analytical precision achievable by zircon U-Pb geochronological techniques have allowed to resolve complexities of zircon crystallization histories in magmatic rocks to an unprecedented level. A number of studies have strived to link resolvable dispersion in zircon age spectra of samples from fossil magmatic systems to the physical parameters of their parent magma bodies. However, the methodologies developed have so far been limited to reproduce the effect of simple thermal histories on the final distribution of zircon ages. In this work we take a more nuanced approach, fine-tuning a thermodynamics-based zircon saturation model to predict the relative distribution of zircon ages in samples from silicic magma reservoirs experiencing open-system processes (e.g. heat/mass addition, mechanical mixing). Employing the MATLAB package (AgeSpectraAnalyst) presented in this contribution:•Users can forward model the effect that diverse thermal histories and mechanical mixing processes characteristic of silicic magma bodies have on zircon age distributions as measured by high-precision, chemical abrasion thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb geochronology.•Zircon CA-ID-TIMS datasets from silicic magmatic systems can be easily compared with model output to gain semi-quantitative information on thermo-mechanical history of the system of interest.•We demonstrated (Tavazzani et al., in press) that distribution of high-precision zircon ages in crystallized remnants of shallow (∼ 250 MPa), silicic magma reservoirs can discriminate between systems that experienced catastrophic, caldera-forming eruptions and systems that underwent monotonic cooling histories.

Garnet petrochronology reveals the lifetime and dynamics of phonolitic magma chambers at Somma-Vesuvius.

Somma-Vesuvius is one of the most iconic active volcanoes with historic and archeological records of numerous hazardous eruptions. Petrologic studies of eruptive products provide insights into the evolution of the magma reservoir before eruption. Here, we quantify the duration of shallow crustal storage and document the evolution of phonolitic magmas before major eruptions of Somma-Vesuvius. Garnet uranium-thorium petrochronology suggests progressively shorter pre-eruption residence times throughout the lifetime of the volcano. Residence times mirror the repose intervals between eruptions, implying that distinct phonolite magma batches were present throughout most of the volcano's evolution, thereby controlling the eruption dynamics by preventing the ascent of mafic magmas from longer-lived and deeper reservoirs. Frequent lower-energy eruptions during the recent history sample this deeper reservoir and suggest that future Plinian eruptions are unlikely without centuries of volcanic quiescence. Crystal residence times from other volcanoes reveal that long-lived deep-seated reservoirs and transient upper crustal magma chambers are common features of subvolcanic plumbing systems.

Long-term repeatability and interlaboratory reproducibility of high-precision ID-TIMS U-Pb geochronology.

Age determination of minerals using the U-Pb technique is widely used to quantify time in Earth's history. A number of geochronology laboratories produce the highest precision U-Pb dates employing the EARTHTIME 202Pb-205Pb-233U-235U tracer solution for isotope dilution, and the EARTHTIME ET100 and ET2000 solutions for system calibration and laboratory intercalibration. Here, we report ET100 and ET2000 solution data from the geochronology laboratory of University of Geneva obtained between 2008 and 2021 and compare the most recent data with results from the geochronology laboratories of Princeton University and ETH Zürich. This compilation demonstrates that (i) the choice of the thermal ionization mass spectrometer model has no influence on precision and accuracy of the data; (ii) the often observed excess scatter of apparent ET100 solution 206Pb/238U dates can be mitigated by more careful tracer-sample equilibration; and (iii) natural zircon reference materials are not suitable for evaluating intra-laboratory repeatability and inter-laboratory reproducibility, since they combine several phenomena of natural system complexities (especially domains of different age within the same zircon grain, and residual loss of radiogenic lead in domains of high decay damage after chemical abrasion pre-treatment). We provide our best estimates of apparent dates for the ET100 solution (206Pb/238U date, 100.173 ± 0.003 Ma), for ET2000 solution (207Pb/206Pb date, 1999.935 ± 0.063 Ma), as well as for natural reference zircon Temora-2 (206Pb/238U date, 417.353 ± 0.052 Ma). These data will allow U-Pb laboratories to evaluate their analytical performance and to independently calibrate non-EARTHTIME tracer solutions in use.

Precise initial abundance of Niobium-92 in the Solar System and implications for p-process nucleosynthesis.

The niobium-92-zirconium-92 (92Nb-92Zr) decay system with a half-life of 37 Ma has great potential to date the evolution of planetary materials in the early Solar System. Moreover, the initial abundance of the p-process isotope 92Nb in the Solar System is important for quantifying the contribution of p-process nucleosynthesis in astrophysical models. Current estimates of the initial 92Nb/93Nb ratios have large uncertainties compromising the use of the 92Nb-92Zr cosmochronometer and leaving nucleosynthetic models poorly constrained. Here, the initial 92Nb abundance is determined to high precision by combining the 92Nb-92Zr systematics of cogenetic rutiles and zircons from mesosiderites with U-Pb dating of the same zircons. The mineral pair indicates that the 92Nb/93Nb ratio of the Solar System started with (1.66 ± 0.10) × 10-5, and their 92Zr/90Zr ratios can be explained by a three-stage Nb-Zr evolution on the mesosiderite parent body. Because of the improvement by a factor of 6 of the precision of the initial Solar System 92Nb/93Nb, we can show that the presence of 92Nb in the early Solar System provides further evidence that both type Ia supernovae and core-collapse supernovae contributed to the light p-process nuclei.

Rapid heterogeneous assembly of multiple magma reservoirs prior to Yellowstone supereruptions.

Large-volume caldera-forming eruptions of silicic magmas are an important feature of continental volcanism. The timescales and mechanisms of assembly of the magma reservoirs that feed such eruptions as well as the durations and physical conditions of upper-crustal storage remain highly debated topics in volcanology. Here we explore a comprehensive data set of isotopic (O, Hf) and chemical proxies in precisely U-Pb dated zircon crystals from all caldera-forming eruptions of Yellowstone supervolcano. Analysed zircons record rapid assembly of multiple magma reservoirs by repeated injections of isotopically heterogeneous magma batches and short pre-eruption storage times of 10(3) to 10(4) years. Decoupled oxygen-hafnium isotope systematics suggest a complex source for these magmas involving variable amounts of differentiated mantle-derived melt, Archean crust and hydrothermally altered shallow-crustal rocks. These data demonstrate that complex magma reservoirs with multiple sub-chambers are a common feature of rift- and hotspot related supervolcanoes. The short duration of reservoir assembly documents rapid crustal remelting and two to three orders of magnitude higher magma production rates beneath Yellowstone compared to continental arc volcanoes. The short pre-eruption storage times further suggest that the detection of voluminous reservoirs of eruptible magma beneath active supervolcanoes may only be possible prior to an impending eruption.

Mass spectrometry in Earth sciences: the precise and accurate measurement of time.

Precise determinations of the isotopic compositions of a variety of elements is a widely applied tool in Earth sciences. Isotope ratios are used to quantify rates of geological processes that occurred during the previous 4.5 billion years, and also at the present time. An outstanding application is geochronology, which utilizes the production of radiogenic daughter isotopes by the radioactive decay of parent isotopes. Geochronological tools, involving isotopic analysis of selected elements from smallest volumes of minerals by thermal ionization mass spectrometry, provide precise and accurate measurements of time throughout the geological history of our planet over nine orders of magnitude, from the accretion of the proto-planetary disk, to the timing of the last glaciation. This article summarizes the recent efforts of the Isotope Geochemistry, Geochronology and Thermochronology research group at the University of Geneva to advance the U-Pb geochronological tool to achieve unprecedented precision and accuracy, and presents two examples of its application to two significant open questions in Earth sciences: what are the triggers and timescales of volcanic supereruptions, and what were the causes of mass extinctions in the geological past, driven by global climatic and environmental deterioration?