Young Silicic Magmatism of the Greater Caucasus, Russia, with implication for its delamination origin based on zircon petrochronology and thermomechanical modeling

https://doi.org/10.1016/j.jvolgeores.2021.107173Get rights and content

Highlights

  • Greater Caucasus is one of the best-preserved examples of silicic magmatism in a continental collision environments.

  • We report results of a multi-method zircon petrochronologic (U--Pb, O--Hf isotopes) and whole rock geochemistry.

  • Thermomechanical modeling describe tempo and volume of magmatism.

  • Voluminous silicic volcanism at Chegem, Tyrnyauz and Elbrus are due to a relict slab detachment.

Abstract

Several kilometers of rapid uplift in the past 2–3 million years in the Greater Caucasus in Russia has produced world-class exposures of Pliocene-Pleistocene age granites and ignimbrites at the Tyrnyauz and Chegem volcanic centers just to the east of Europe's highest mountain and active volcano, Mount Elbrus. This represents one of the world's best-preserved examples of silicic magmatism in a continental collision environment. We report results of a multi-method zircon petrochronologic (U--Pb, O--Hf isotopes, trace elements) investigation of six ignimbrites, lavas, and associated granodioritic porphyries from these localities. We observe two pulses of magmatism at 2.92 and 1.98 Ma related to Chegem and Tyrnyauz, respectively. High precision CA-ID-TIMS dating of zircons from the top and bottom of the Chegem ignimbrite and its associated porphyry yield indistinguishable age spectra recording 160 kyr of magma assembly and differentiation, with 2.9181 ± 0.0014 Ma eruption age as is constrained by the youngest dated zircon population. Together with overlapping O and Hf isotopic values, this suggests that they represent a large volume of pre-eruptively homogenized magma with ~20% crystals. Zircons have isotopically diverse cores (δ18O = +3.4–6.7‰, εHf = +0.8 to +5.6) while the rims (+5.75 ± 0.20‰, εHf = 3.3 ± 0.7) are in isotopic equilibrium with associated unaltered glass and major phenocryst phases. The neighboring and younger ignimbrites and granites of the nearby Elbrus and Tyrnyauz centers have overlapping higher δ18O and lower εHf zircon values and cluster narrowly in age around 1.98 Ma, suggesting a common source for all these ignimbrites, likely in the Tyrnyauz area. Isotopic data for Chegem, Elbrus and Tyrnyauz zircons and rocks require the contribution of a high proportion (25–55%) of Paleozoic crust to magma petrogenesis.

The origin of sudden spikes of voluminous silicic magmatism is further investigated using thermomechanical modeling of the collisional environment. We show that the voluminous silicic volcanism at Chegem, Tyrnyauz and Elbrus is best explained by heating of the lower crust by asthenospheric upwelling after a relict slab detached and sank into the mantle at approximately 5 Ma. This timing is coincident with both the beginning of rapid uplift of the Greater Caucasus and initiation of major volcanism. After this initial melting, continuing delamination of the lithosphere and lower crust from the subducting plate produced broad mantle melting and basaltic volcanism, potentially explaining continued regional-scale magmatism in both the Greater and Lesser Caucasus.

Introduction

Magmatism during the continent-continent collisions is perhaps the least understood source of felsic magmas on Earth today, with no one melting mechanism able to explain all magma production in these environments. In Tibet, the most famous continental collision zone, the production of silicic magmas has been attributed both to radiogenic heating of thickened crust (e.g., Bea, 2012) and to partial melting of hydrated oceanic crust residing in the upper mantle after being subducted prior to collision (Mo et al., 2008). In the Greater and Lesser Caucasus mountains, which comprise the active collision zone between the Arabian and Eurasian plates (Fig. 1), there is also significant syn-collisional volcanism. This magmatism has been suggested by a variety of workers using both geochemical and geophysical evidence to be primarily caused by post-collision slab breakoff and lithospheric delamination leading to mantle decompression and melting (Pearce et al., 1990; Keskin, 2003; Şengör et al., 2008; Koulakov et al., 2012; Sugden et al., 2019). However, mantle melting cannot be the sole source of Caucasus magmas, as Sr and Nd isotopes have been used to demonstrate that melting of the Paleozoic crust must also significantly contribute to erupted material (e.g., Lebedev and Vashakidze, 2014).

At the northern end of the collision zone in the Greater Caucasus, rapid uplift since the Pleistocene has produced world-class exposures of syn-orogenic silicic intrusions and volcanism near Mount Elbrus (at 5642 m, the highest mountain and active volcano in Europe). These include the 2 Ma, 4 × 5 km Eldjurta Granite stock at the Tyrnyauz mining district, 35 km east of Elbrus (one of the world's youngest granites). The Eldjurta Granite is affiliated with world class Mo-W mineralization and is crosscut by brecciated rhyolitic bodies, vitrophyres, and aplites (Milanovsky and Koronovsky, 1973; Kostitsyn and Kremenetsky, 1995; Grün et al., 1999). About 10 km southeast of Tyrnyauz, the ~11 × 15 km-wide, 2.9 Ma Chegem Caldera has been uplifted and dissected by erosion, exposing a 2 km-thick intracaldera tuff and a central granodioritic porphyry which is texturally and compositionally similar to the Eldjurta granite, fueling speculation that the two may be surface expressions of a larger underlying batholith (Fig. 2, Fig. 3; Milanovsky and Koronovsky, 1973; Lipman et al., 1993; Gazis, 1994; Gazis et al., 1995, Gazis et al., 1996). The quality of these exposures has made them the focus of important early debates about the formation of welded ignimbrites and texturally similar lavas (“tufo-lavas”, Levinson-Lessing, 1913; Milanovsky et al., 1962; Koronovsky et al., 1982; Koronovsky and Demina, 2007, Bindeman et al. 2021 in press, Fig. 1 therein).

Section snippets

Objectives of the present study

This study investigates the origins of recent silicic magmatism in the Greater Caucasus via two interrelated goals:

  • 1)

    The young age and good preservation of Greater Caucasus lavas, tuffs, and plutons presents a unique opportunity to conduct a detailed study of its origins. Newly developed CA-ID-TIMS methods can reduce 2-sigma uncertainties of calculated ages for a single zircon crystal to ~0.1% (e.g., ±2 kyr on a 2 Myr unit). For samples with ages of less than 3 Ma as in the studied region (Fig. 2

Regional background: Silicic magmatism in the Elbrus-Tyrnyauz-Chegem region of the Greater Caucasus

The Greater Caucasus of southern Russia formed as the result of the collision between the Arabian, Anatolian, and Eurasian plates, which has been ongoing since the initial closure of the local Tethys ocean in the Jurassic. The current phase of the collision from 35 Ma to the present has culminated in rapid uplift of the Greater Caucasus which has accelerated in the last 5 million years (Ershov and Nikishin, 2004; Zakariadze et al., 2007; Adamia et al., 2011; Vincent et al., 2020; Fig. 1). The

Sampling

Ignimbrites, volcanic necks, and lavas in the Elbrus (samples ELB-1 to ELB-19) and Tyrnyauz (samples ELB-20, −21) areas are preserved as isolated remnants with thicknesses reaching 50 m. We sampled the most prominent ignimbrites on the northern slopes of modern Elbrus edifice, at Tuzluk Mt. (sample ELB-1), and the Lower (ELB-13, −12) and Upper (ELB-9, −11) Stone Mushroom ignimbrites near the upper Birdjalisu river (Fig. 2). Also at Elbrus, we sampled an additional 10 lavas of varying

U-Pb ages

Zircon ages from Chegem and five other regional ignimbrites, two granites, and two lavas were determined by LA-ICP-MS, SHRIMP and CA-ID-TIMS (Bindeman et al. in press, Tables 2, A1-2 therein). Zircons were generally U-rich, enabling precise dating. The 230Th-corrected 206Pb/238U ignimbrite dates cluster around two age groups: 2.92 Ma for the Chegem units and 1.96–1.98 Ma for the Elbrus and Tyrnyauz/Eldjurta samples (Fig. 4). The Elbrus Upper and Lower ignimbrites, mapped by Gazeev and Gurbanov

Trace elements

The trace element geochemistry of the Chegem rhyolites straddles the boundary between collisional and arc-related magmas (Bindeman et al. in press, Fig. 5 therein; Pearce et al., 1984). Basalts and dacites of the Elbrus, Tyrnyauz, and Kazbek volcanic centers are calc-alkaline with elevated alkalis and show enrichments in large ion lithophile elements such as Rb, Sr, Ba, U, and Pb and depletions in the high field strength elements Ti, Zr, and Nb (Lebedev et al., 2010; Lipman et al., 1993;

Conclusions

We present a broad, multi-disciplinary study of recent volcanism in the Greater Caucasus focusing on the young and well-preserved Chegem, Elbrus and Tyrnyauz volcanic centers. We observe clustering of zircon ages at 2.91 and 1.98 Ma, signifying the presence of two pulses of voluminous silicic magmatism. It is likely that a now-eroded caldera at Tyrnyauz is the source of 1.98 Ma voluminous ignimbrites with the Eldjurta Granite representing its parental magma. Likewise, Chegem and Zayukovo

Author statement

Ilya Bindeman: Fieldwork, O-C-H isotope analyses, zircon extraction, paper writing and editing Dylan Colón: Thermomechanical modeling, paper editing Jörn Wotzlaw: zircon ID TIMS and LA-ICP-MS analysis, data plotting and paper editing. Richard Stern: SIMS analysis of zircons for O isotopes. Massimo Chiaradia: Hf isotope analysis of zircons by LA-ICP-MS.: Marcel Guillong: LA-ICP-MS analysis of zircons].

Declaration of Competing Interest

None.

Acknowledgements

Supported by Russian National Science Fund (RNF) grant #19-17-00241. We thank S. Popov, N. Nekrylov and O.E. Melnik for help with fieldwork. Taras Gerya for access to the EULER computer cluster network and use of his modeling code, Caltech sample repository is thanked for access to sample collection associated with C. C. Gazis and H.P. Taylor et al. 1995 study. Careful and very insightful reviews by Juliana Troch, Kelly Russell, Gerhard Woerner and Shan de Silva and an anonymous reviewer from

References (125)

  • J.C. Hess et al.

    The cooling history of the late Pliocene Eldzhurtinskiy granite (Caucasus, Russia) and the thermochronological potential of grain size/age relationships

    Earth Planet. Sci. Lett.

    (1993)
  • M.R. Hudak et al.

    Conditions of pinnacle formation and glass hydration in cooling ignimbrite sheets from H and O isotope systematics at Crater Lake and the Valley of ten Thousand Smokes

    Earth Planet. Sci. Lett.

    (2018)
  • O. Karakas et al.

    Melt evolution and residence in extending crust: thermal modeling of the crust and crustal magmas

    Earth Planet. Sci. Lett.

    (2015)
  • L. Karlstrom et al.

    The effect of pressurized magma chamber growth on melt migration and pre-caldera vent locations through time at Mount Mazama, Crater Lake, Oregon

    Earth Planet. Sci. Lett.

    (2015)
  • P.W. Lipman et al.

    2.8-Ma ash-flow caldera at Chegem River in the northern Caucasus Mountains (Russia), contemporaneous granites, and associated ore deposits

    J. Volcanol. Geotherm. Res.

    (1993)
  • J.M. Mattinson

    Zircon U–Pb chemical abrasion (“CA-TIMS”) method: combined annealing and multi-step partial dissolution analysis for improved precision and accuracy of zircon ages

    Chem. Geol.

    (2005)
  • X. Mo et al.

    Contribution of syncollisional felsic magmatism to continental crustal growth: a case study of the Paleogene Linzizong volcanic Succession in southern Tibet

    Chem. Geol.

    (2008)
  • T. Mumladze et al.

    Subducted, detached, and torn slabs beneath the Greater Caucasus

    Geo. Res. J.

    (2015)
  • J.A. Pearce et al.

    Genesis of collision volcanism in Eastern Anatolia, Turkey

    J. Volcanol. Geotherm. Res.

    (1990)
  • H. Philip et al.

    The Caucasus: an actual example of the initial stages of continental collision

    Tectonophysics

    (1989)
  • A.W. Rempel et al.

    A model for the development of stable isotopic water signatures of tephra deposited on ice following subglacial caldera collapse

    J. Volcanol. Geotherm. Res.

    (2019)
  • D. Rubatto et al.

    Experimental zircon/melt and zircon/garnet trace element partitioning and implications for the geochronology of crustal rocks

    Chem. Geol.

    (2007)
  • A.M.C. Şengör et al.

    Eastern Turkish high plateau as a small Turkic-type orogen: Implications for post-collisional crust-forming processes in Turkic-type orogens

    Earth Sci. Rev.

    (2008)
  • S. Adamia et al.

    Geology of the Caucasus: A Review

    Turk. J. Earth Sci.

    (2011)
  • G.D. Afanasyev et al.

    Materials to justification of ages for boundaries between several geologic systems and epochs

    Izvestiya USSR Acad. Sci.

    (1963)
  • N.L. Andersen et al.

    Repeated rhyolite eruption from heterogeneous hot zones embedded within a cool, shallow Magma reservoir

    J. Geophys. Res. Solid Earth

    (2019)
  • C. Annen et al.

    The Genesis of intermediate and silicic magmas in deep crustal hot zones

    J. Petrol.

    (2006)
  • A.S. Avdeenko et al.

    Isotopic (Nd, Sr, δ18O) signatures of the late Miocene high Sr-Ba granitoids in the Mineral’nye Vody area, Northern Caucasus: evidence of the protolith nature

    Dokl. Earth Sci.

    (2008)
  • I.N. Bindeman et al.

    Zircon survival, rebirth and recycling during crustal melting, magma crystallization, and mixing based on numerical modelling

    J. Petrol.

    (2016)
  • I.N. Bindeman et al.

    Rhyolites - hard to produce but easy to recycle and sequester: bringing together microgeochemical observations and numerical models

    Geosphere

    (2014)
  • I.N. Bindeman et al.

    Origin and Evolution of Yellowstone silicic magmatism based on the ion microprobe analysis of isotopically-zoned zircons

    J. Petrol.

    (2008)
  • I.N. Bindeman et al.

    Isotopic and petrologic investigation, and a thermomechanical model of genesis of large-volume Rhyolites in arc environments: Karymshina Volcanic complex, Kamchatka, Russia

    Front. Earth Sci.

    (2019)
  • I.N. Bindeman et al.

    Geochronology and Geochemistry Data of the Elbrus, Tyrnyauz, and Chegem Magmatic Centers, Greater Caucasus, Russia

    (2021)
  • O.A. Bogatikov et al.

    The Elbrus caldera in the northern Caucasus

    Dokl. Earth Sci.

    (1998)
  • J.F. Bowring et al.

    Engineering cyber infrastructure for U–Pb geochronology: Tripoli and U–Pb_Redux

    Geochem. Geophys. Geosyst.

    (2011)
  • S.N. Bubnov

    The Chronology of Eruptions and Magma Sources for Youngest Volcanic Centers of the Greater Caucasus

  • C. Cauley

    Digital Imagery and Geologic Map of CHEGEM Caldera, Russia. HSGC Report Number 15-24

    (2015)
  • I.V. Chernyshev et al.

    Isotopic geochronology of Quaternary volcanic eruptions in the Greater Caucasus

    Geochem. Int.

    (2002)
  • I.V. Chernyshev et al.

    Two stages of explosive volcanism of the Elbrus area: Geochronology, petrochemical and isotopic-geochemical characteristics of volcanic rocks, and their role in the Neogene-Quaternary evolution of the Greater Caucasus

    Stratigr. Geol. Correl.

    (2014)
  • M. Chiaradia

    Crustal thickness control on Sr/Y signatures of recent arc magmas: an Earth scale perspective

    Sci. Rep.

    (2015)
  • S.L. Chung et al.

    Adakites from continental collision zones: Melting of thickened lower crust beneath southern Tibet

    Geology

    (2003)
  • D.P. Colón et al.

    Isotopically diverse rhyolites coeval with the Columbia River Flood Basalts: evidence for mantle plume interaction with the continental crust

    Terra Nova

    (2015)
  • D.P. Colón et al.

    Thermomechanical modeling of the formation of a multilevel, crustal-scale magmatic system by the Yellowstone plume

    Geophys. Res. Lett.

    (2018)
  • D.P. Colón et al.

    Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology, oxygen isotope, and hafnium isotope analyses of zircon

    Contrib. Mineral. Petrol.

    (2018)
  • E. Cowgill et al.

    Relict basin closure and crustal shortening budgets during continental collision: an example from Caucasus sediment provenance

    Tectonics

    (2016)
  • F.X. D’Abzac et al.

    Hf isotope analysis of small zircon and baddeleyite grains by conventional multi collector-inductively coupled plasma-mass spectrometry

    Chem. Geol.

    (2015)
  • M.N. Ducea et al.

    High- volume magmatic events in subduction systems

    Elements

    (2015)
  • J. Dufek et al.

    Lower crustal magma genesis and preservation: a stochastic framework for the evaluation of basalt–crust interaction

    J. Petrol.

    (2005)
  • A.V. Ershov et al.

    Recent Geodynamics of the Caucasus-Arabia-East Africa Region

    Geotectonics

    (2004)
  • J.M. Ferry et al.

    New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers

    Contrib. Mineral. Petrol.

    (2007)
  • Cited by (12)

    • Temporal changes in subduction- to collision-related magmatism in the Neotethyan orogen: The Southeast Iran example

      2022, Earth-Science Reviews
      Citation Excerpt :

      Convergent margins continue to deliver arc-like magmas even during and for some time after continental collision as continental lithosphere begins to subduct. In several places, very young or ongoing magmatic activity occurs in collision zone that is not related to oceanic subduction, such as the Carpathians or the Caucasus (e.g., Bindeman et al., 2021; Seghedi and Downes, 2011). Collision-related magmatic products have similar bulk chemistry (in both major-trace elements and isotopes) to subduction magmatism, and are dominated by calc-alkaline andesites, dacites and rhyolites as well as their plutonic equivalents.

    • Geochronology and geochemistry data for the Elbrus, Tyrnyauz, and Chegem magmatic centers, Greater Caucasus, Russia

      2021, Data in Brief
      Citation Excerpt :

      Zircon grains were imaged via cathodoluminescence and scanning electron microscope (Fig. 4) prior to isotopic and chemical analysis. Details of analytical methods are provided in [1]. Three epoxy mounts were prepared.

    View all citing articles on Scopus
    View full text