Zircon as a recorder of contrasting magma recharge and eruptive recurrence patterns
Graphical abstract
Introduction
Understanding magma evolution and storage conditions is an indispensable prerequisite to elucidate eruptive behavior and thus, ultimately, volcanic hazards. Crystal mush-dominated transcrustal magmatic plumbing systems are widely regarded to form by prolonged, pulsed magma injection (e.g., Cashman et al., 2017, and references therein). The current state of these magmatic systems can, although with limited spatial resolution, be imaged by geophysical methods (e.g., Magee et al., 2018). However, their evolution through time, in particular with respect to the debate of protracted “cold” magma storage below the rheological lockup temperature of a crystal mush (Cooper and Kent, 2014) vs. “warm” magma storage under eruptible conditions (Barboni et al., 2016), can only be inferred from geological and geochemical records. For intermediate to evolved magmas, zircon is uniquely suited to record their thermochemical conditions as its crystallization is protracted, and diffusive re-equilibration is insignificant for zircon even when stored at magmatic temperatures for extended timescales (e.g., Schmitt, 2011). To interpret this record in terms of magma evolution and storage conditions, empirical zircon crystallization age spectra can be compared to synthetic ones generated by varying parameters such as magma volume and recharge rates in model scenarios of magma accumulation in upper crustal settings that are open for recharge of hot, zircon-undersaturated magma (Bindeman and Melnik, 2016; Caricchi et al., 2014; Tierney et al., 2016). Geothermometric data have been employed only very recently to refine such models (Lukács et al., 2021; Weber et al., 2020).
Here, we investigate two neighboring stratovolcanoes, Mt. Hasan and Mt. Erciyes (Turkey), formed in the post-collisional setting of late Quaternary volcanism within the Anatolian continental microplate, and complement published secondary ion mass spectrometry (SIMS) U–Th–Pb zircon rim and interior crystallization ages for >1000 crystals (Friedrichs et al., 2020c) with new trace element analyses on the same spots to unravel their magmatic histories, including Ti-in-zircon geothermometry. Applying the same methods and techniques to zircon from these two volcanoes mitigates bias and can uniquely reveal different styles of zircon crystallization that reflect on pre-eruptive crystallization and storage conditions. These data are the basis for comparative thermochemical modeling of zircon crystallization to reconstruct different magma recharge histories, indicating that different eruptive recurrence and styles are controlled by the thermal inertia of subvolcanic magma reservoirs where “cold” and “warm” storage are end-members related to waxing, peak, and waning stages in a volcano's evolution.
Section snippets
Regional setting and volcanostratigraphy
Mt. Hasan and Mt. Erciyes, separated by ∼120 km, are the two largest stratovolcanic complexes of the post-collisional Central Anatolian Volcanic Province (CAVP; Fig. 1) where widespread ignimbrites erupted between ca. 10 and ca. 2.5 Ma (Aydar et al., 2012). Quaternary activity within the CAVP mainly produced basaltic and compositionally bimodal volcanic fields as well as polygenetic cones. Both Mt. Hasan and Mt. Erciyes are located on major strike-slip faults within the trans-tensional tectonic
Samples and study design
Thirty-nine volcanic rock samples, predominantly lava and composite pumice, were collected at Mt. Hasan and Mt. Erciyes (Friedrichs et al., 2020c, 37 samples; Schmitt et al., 2014, two samples), targeting the stratigraphically and morphologically youngest units (Supplementary Table A.1). After crushing, disc-milling, sieving (<125 μm), and hydrodynamic enrichment, zircon crystals were hand-picked under a binocular microscope, and Fe–Ti oxide minerals extracted by magnetic separation. Zircon
Zircon geochronology
Individual sub- to euhedral zircon crystals of both Mt. Hasan and Mt. Erciyes display typical magmatic features of oscillatory and sector zonation in CL (Figs. 2, B.2–3). SIMS depth profiling of zircon rims (Table A.2) generally indicates continuous crystallization at linear growth rates of ∼0.05–2 μm/ka, with significant intra-crystal variability observed only for Mt. Erciyes. Intra- and inter-crystal variability in zircon U concentrations is minor for Mt. Hasan (∼100–600 ppm) compared to Mt.
Contrasting zircon thermochemical records in two adjacent volcanoes
Zircon age patterns for rims and interiors as well as corresponding melt differentiation indices are clearly distinct for Mt. Hasan and Mt. Erciyes (Fig. 3, Fig. 4, Fig. 5), despite similarities in tectonic setting (e.g., Dilek and Sandvol, 2009) and composition (e.g., Aydar and Gourgaud, 1998; Şen et al., 2003). Moreover, both data sets were acquired under identical conditions, so that sampling or instrumental bias is largely mitigated. There are, however, clear differences in eruptive
Conclusions
A combination of zircon petrochronology data and thermochemical modeling was applied to unravel contrasting magma storage conditions governing distinct eruptive behaviors of two neighboring, long-lived, and active stratovolcanoes in very similar tectonic settings. Mt. Hasan, revealed to be still in its peak stage, is underlain by a large, thermally inert magma system as evident in time-invariant zircon crystallization temperature and trace element concentration ranges, and characterized by
CRediT authorship contribution statement
Bjarne Friedrichs: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft. Axel K. Schmitt: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing. Oscar M. Lovera: Formal analysis, Investigation, Methodology, Software, Validation, Writing – review & editing. Gokhan Atıcı: Funding
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
We thank Anne Sturm for help with zircon separation, Andreas Kronz for supervising EMP analysis, John Hora for guidance through Fe–Ti oxide data evaluation, and Alejandro Cisneros de León for comments on an early version of the manuscript. This work was supported by DFG (German Research Foundation) grant SCHM2521/3-1, MTA (General Directorate of Mineral Research and Exploration of Turkey), and a DAAD (German Academic Exchange Service) doctoral scholarship to BF. Maps were created using ArcGIS®
References (50)
- et al.
The geology of Mount Hasan stratovolcano, central Anatolia, Turkey
J. Volcanol. Geotherm. Res.
(1998) - et al.
Correlation of ignimbrites in the central Anatolian volcanic province using zircon and plagioclase ages and zircon compositions
J. Volcanol. Geotherm. Res.
(2012) - et al.
Zircon saturation re-revisited
Chem. Geol.
(2013) - et al.
Millennial to decadal magma evolution in an arc volcano from zircon and tephra of the 2016 Santiaguito eruption (Guatemala)
Lithos
(2019) - et al.
Evaluating the relationship between climate change and volcanism
Earth-Sci. Rev.
(2018) - et al.
Application of combined U-Th-disequilibrium/U-Pb and (U-Th)/He zircon dating to tephrochronology
Quat. Geochronol.
(2017) - et al.
Petrogenesis of mafic-silicic lavas at Mt. Erciyes, central Anatolia, Turkey
J. Volcanol. Geotherm. Res.
(2013) - et al.
Late Pleistocene eruptive recurrence in the post-collisional Mt. Hasan stratovolcanic complex (Central Anatolia) revealed by zircon double-dating
J. Volcanol. Geotherm. Res.
(2020) - et al.
Sequence modeling in zircon double-dating of early Holocene Mt. Erciyes domes (Central Anatolia)
Quat. Geochronol.
(2021) - et al.
New insights into source and dispersal of Mediterranean S1 tephra, an early Holocene marker horizon erupted at Mt. Erciyes (Turkey)
Quat. Sci. Rev.
(2020)
U–Th whole rock data and high spatial resolution U–Th disequilibrium and U–Pb zircon ages of Mt. Erciyes and Mt. Hasan Quaternary stratovolcanic complexes (Central Anatolia)
Data Brief
An Excel-based tool for evaluating and visualizing geothermobarometry data
Comput. Geosci.
Parameters controlling the density of welded ignimbrites—a case study on the Incesu Ignimbrite, Cappadocia, Central Anatolia
Geochemistry
Zircon geochronology suggests a long-living and active magmatic system beneath the Ciomadul volcanic dome field (eastern-central Europe)
Earth Planet. Sci. Lett.
Chlorine degassing constrained by cosmogenic 36Cl and radiocarbon dating of early Holocene rhyodacitic lava domes on Erciyes stratovolcano, central Turkey
J. Volcanol. Geotherm. Res.
Glaciations and paleoclimate of Mount Erciyes, central Turkey, since the Last Glacial Maximum, inferred from 36Cl cosmogenic dating and glacier modeling
Quat. Sci. Rev.
Impact of volcanism on the evolution of Lake Van (eastern Anatolia) III: periodic (Nemrut) vs. episodic (Süphan) explosive eruptions and climate forcing reflected in a tephra gap between ca. 14 ka and ca. 30 ka
J. Volcanol. Geotherm. Res.
Volcanological evolution of Mount Erciyes stratovolcano and origin of the Valibaba Tepe ignimbrite (Central Anatolia, Turkey)
J. Volcanol. Geotherm. Res.
Unraveling the electrical conductivity structure to decipher the hydrothermal system beneath the Mt. Hasan composite volcano and its vicinity, SW Cappadocia, Turkey
J. Volcanol. Geotherm. Res.
Insight into zircon REE oxy-barometers: a lattice strain model perspective
Earth Planet. Sci. Lett.
Repeated rhyolite eruption from heterogeneous hot zones embedded within a cool, shallow magma reservoir
J. Geophys. Res., Solid Earth
Curie point depth investigation of central Anatolia, Turkey
Pure Appl. Geophys.
Mg/Mn partitioning as a test for equilibrium between coexisting Fe-Ti oxides
Am. Mineral.
Warm storage for arc magmas
Proc. Natl. Acad. Sci. USA
Zircon survival, rebirth and recycling during crustal melting, magma crystallization, and mixing based on numerical modelling
J. Petrol.
Cited by (7)
Hydrothermal fluid circulation within the restless structural frame of Hasandağ volcanic system (Central Anatolia, Türkiye) inferred from Self-potential, CO<inf>2</inf>, and temperature measurements
2024, Journal of Volcanology and Geothermal ResearchExotic REE behaviors of zircon in the Koktokay No. 3 granitic pegmatite, Xinjiang, northwest China
2023, Ore Geology ReviewsTiming and recurrence intervals for voluminous silicic eruptions from Amatitlán caldera (Guatemala)
2023, Quaternary Science ReviewsZirconium and its stable isotopes in igneous systems
2023, Earth-Science ReviewsCitation Excerpt :This implies that large Zr isotope fractionation between major-element mineral and melt may be uncommon during magma differentiation. In comparison, the estimated growth rates of zircon are much lower, mostly ranging from 10−17 to 10−13 cm/s (Friedrichs et al., 2021; Schmitt et al., 2011; Sorokin et al., 2022; Watson, 1996; Zhang and Xu, 2016). The calculated R/D ratio for zircon is still controversial, depending on the Zr diffusivity data in use (Guo et al., 2020; Méheut et al., 2021).