Articles | Volume 17, issue 1
Clim. Past, 17, 229–239, 2021
Clim. Past, 17, 229–239, 2021

Research article 18 Jan 2021

Research article | 18 Jan 2021

Exploring a link between the Middle Eocene Climatic Optimum and Neotethys continental arc flare-up

Exploring a link between the Middle Eocene Climatic Optimum and Neotethys continental arc flare-up
Annique van der Boon1,a, Klaudia F. Kuiper2, Robin van der Ploeg1,b, Margot J. Cramwinckel1,c, Maryam Honarmand3, Appy Sluijs1, and Wout Krijgsman1 Annique van der Boon et al.
  • 1Department of Earth Sciences, Utrecht University, Princetonlaan 8a, 3584 CB Utrecht, the Netherlands
  • 2Department of Earth Sciences, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1085, 1081 HV Amsterdam, the Netherlands
  • 3Department of Earth Sciences, Institute for Advanced Studies in Basic Sciences (IASBS), P.O. Box 45195-1159, Zanjan, Iran
  • anow at: Geomagnetic Laboratory, Oliver Lodge Building, Department of Physics, Oxford Street, University of Liverpool, Liverpool, L69 7ZE, UK
  • bnow at: Shell Global Solutions International B.V., Grasweg 31, 1031 HW Amsterdam, the Netherlands
  • cnow at: National Oceanography Centre Southampton, University of Southampton Waterfront Campus, European Way, Southampton, SO14 3ZH, UK

Correspondence: Annique van der Boon (


The Middle Eocene Climatic Optimum (MECO), a ∼500kyr episode of global warming that initiated at  40.5 Ma, is postulated to be driven by a net increase in volcanic carbon input, but a direct source has not been identified. Here we show, based on new and previously published radiometric ages of volcanic rocks, that the interval spanning the MECO corresponds to a massive increase in continental arc volcanism in Iran and Azerbaijan. Ages of Eocene igneous rocks in all volcanic provinces of Iran cluster around 40 Ma, very close to the peak warming phase of the MECO. Based on the spatial extent and volume of the volcanic rocks as well as the carbonaceous lithology in which they are emplaced, we estimate the total amount of CO2 that could have been released at this time corresponds to between 1052 and 12 565 Pg carbon. This is compatible with the estimated carbon release during the MECO. Although the uncertainty in both individual ages, and the spread in the compilation of ages, is larger than the duration of the MECO, a flare-up in Neotethys subduction zone volcanism represents a plausible excess carbon source responsible for MECO warming.

1 Introduction

The Middle Eocene Climatic Optimum (MECO) is characterized by surface and deep ocean warming, both of approximately 2–6 C. MECO warming initiated at ∼40.5Ma, culminating in a short peak warming phase at ∼40.0Ma and terminating at ∼39.9Ma with a comparatively rapid cooling (Bijl et al., 2010; Bohaty et al., 2009; Bohaty and Zachos, 2003; Boscolo Galazzo et al., 2013, 2014; Cramwinckel et al., 2018). The MECO is associated with a rise in atmospheric CO2 concentrations (Bijl et al., 2010; Henehan et al., 2020), extensive deep sea carbonate dissolution (Bohaty et al., 2009) and marine biotic change (Bijl et al., 2010; Cramwinckel et al., 2019; Edgar et al., 2013; Witkowski et al., 2012). The MECO inherently differs from the early Paleogene transient warming events such as the Paleocene–Eocene Thermal Maximum (PETM; ∼56Ma) primarily in its longer duration (∼500kyr) of warming, precluding a sudden trigger but rather suggesting a continued driver (Bohaty and Zachos, 2003; Sluijs et al., 2013). Furthermore, unlike the PETM and similar transients, the MECO is not characterized by a negative δ13C excursion of the exogenic carbon pool, ruling out the input of 13C-depleted organic-sourced carbon as a driver, but suggesting a volcanic source (Bohaty and Zachos, 2003). Reconstructions and simulations of the carbon cycle indeed point to an imbalance in the long-term inorganic carbon cycle during the MECO (Sluijs et al., 2013), caused by enhanced volcanism and sustained by diminished continental silicate weathering (van der Ploeg et al., 2018). However, this scenario is quantitatively far from settled, partly because recent analyses based on foraminifer boron isotope ratios suggest that atmospheric CO2 concentrations rose by significantly less than a doubling and did not rise substantially during the onset of the MECO (Henehan et al., 2020). In addition, a plausible source of excess volcanic CO2 remains to be identified.

Here, we explore a volcanic arc flare-up in the Neotethys subduction zone as a potential source. Arc flare-ups can generate 80 %–90 % of the total volume of igneous rocks in arc systems in periods of a few million years (Ducea and Barton, 2007). During the Eocene, a large flare-up took place in vast areas of present-day Iran (see Fig. 1a), and these volcanic rocks show subduction-related geochemical signatures, representative of continental arc volcanism (e.g., Moghadam et al., 2015; Pang et al., 2013; Verdel et al., 2011). Geologic settings of the Eocene volcanic regions in Iran differ. Extensive magmatism in the Lut block is regarded by Pang et al. (2013) to be the result of post-collisional convective removal of the lithosphere and not directly related to subduction. Volcanism in the Sabzevar zone is linked by Moghadam et al. (2016) to lithospheric delamination, possibly assisted by slab breakoff. In the Talesh–Alborz region, there are conflicting theories on the formation of the volcanic rocks. Asiabanha and Foden (2012) mention a post-collisional transition to a continental arc in their title, but then they describe the volcanism as back-arc volcanism. Van der Boon (2017) gives an overview of proposed conflicting settings for volcanism in the Alborz. It is striking that in most of the areas in Iran, the flare-up is linked to an extensional setting (e.g., Verdel et al., 2011), which makes it different from other flare-ups (e.g., Ducea et al., 2015; Ducea and Barton, 2007).

Figure 1(a) Map showing the outcrop of Eocene volcanic rocks in Iran (modified after Agard et al., 2011, and shapefiles of Sahandi et al., 2014). The five largest areas are shown with white outlines. (b) Sample locations of newly acquired 40Ar∕39Ar ages.

The main volcanic arc associated with the Neotethys subduction zone stretches from Bazman in southeastern Iran towards Azerbaijan in the northwest, where it continues westwards into Armenia, Georgia and Turkey (Van Der Boon et al., 2017). North of the volcanic arc, in the Peri-Tethys basin of Azerbaijan and Russia, thick bentonites and ash layers are found within middle Eocene marine sediments (Beniamovski et al., 2003; Seidov and Alizade, 1966). In the past, it has been hypothesized based on field studies that the middle Eocene part makes up the bulk of the Eocene volcanic succession (e.g., Glaus, 1965) because of the presence of middle Eocene nummulites within the Karaj formation (e.g., Sieber, 1970) and volcanism climaxes during middle Eocene time (Berberian and King, 1981; Davoudzadeh et al., 1997; Verdel, 2009). Sahandi et al. (2014) produced a compilation of geological maps of Iran, which shows that more than half of the outcrop area of igneous rocks in Iran is of Eocene age (see Fig. 1a). The total surface area that is covered by Eocene igneous rocks is almost 70 000 km2 (including units mapped as middle Eocene, Eocene–Oligocene, etc.). A causal relationship between peak volcanism in this region and the MECO has been suggested (Allen and Armstrong, 2008; Kargaranbafghi and Neubauer, 2018), but radioisotopic age constraints to test this hypothesis are insufficient. To quantitatively assess whether volcanism in the Iran–Azerbaijan region could have been a contributor to global warming during the MECO, we present a compilation of new and previously published radiometric ages for volcanic rocks and estimate eruptive volumes of the flare-up in Iran to evaluate how much CO2 could have been released during this continental arc flare-up.

2 Dating the continental arc flare-up of the Neotethys subduction zone

2.1 New 40Ar∕39Ar data

We analyzed 48 samples of Eocene volcanic rocks of the Azerbaijan–Bazman Arc in Iran and Azerbaijan. Lava flows of the Peshtasar Formation were dated by Vincent et al. (2005) and van der Boon et al. (2017), but ages suffered from severe excess argon. Here, we re-dated lava flows from the lower and middle part of the Peshtasar Formation using new instrumentation to check for potential age bias caused by hydrocarbon interferences in previous data. We further dated samples of two ash layers in the Kura basin in Azerbaijan, as well as four volcanic rocks from the Talesh and western Alborz in Iran (see Fig. 1b). Depending on the rock type, groundmass, plagioclase, sanidine, biotite and/or glass was measured (see Table 1). Thin section analysis showed pervasive alteration of volcanic rocks, disqualifying many sampled units for radioisotope dating (see Supplement S1 in the for a comparison of some thin sections). However, eight samples showed no significant alteration and were prepared for 40Ar∕39Ar dating using standard mineral separation techniques including heavy liquid and magnetic separation and handpicking. In general, fractions between 250–500 µm size were taken. For some minerals, both groundmass or glass and plagioclase or biotite could be separated.

Table 1Details of samples for Ar−Ar dating.

NA: not available.

Download Print Version | Download XLSX

Samples were leached with diluted HNO3 and/or HF. Samples were irradiated for 12 and 18 h, respectively, in two irradiations (VU101 in 2014 and VU107 in 2016) at the Oregon State University Triga CLICIT facility, together with Fish Canyon Tuff sanidine as a standard (FCs; 28.201±0.023Ma; Kuiper et al., 2008). After irradiation samples were loaded on Cu trays and run on a 10-collector Helix-MC mass spectrometer with an in-house built extraction with SEAS NP10, St172 and Ti sponge getters and a Lauda cooler run at 70 C, at the Vrije Universiteit Amsterdam. The used cup configuration was either 40Ar on the H2 Faraday cup and 39–36 argon isotopes on compact discrete dynodes or both 40Ar and 39Ar on respectively H2 and H1 Faraday. Gain calibration was done by peak-jumping CO2 in dynamic mode on the different cups (see Monster, 2016, for details). Samples were analyzed using step-heating experiments, while for the ash layers usually single or a few grains were fused in one step and analyzed. Initial measurements were on single or a small number of grains, leading in some samples to very low intensities of 40Ar (3–4 times higher than blanks). In those cases, more grains were loaded in the next experiment. Ages are calculated relative to the age of FCs reported in Kuiper et al. (2008; 28.201±0.023Ma) with decay constants of Min et al. (2000).

Figure 2(a) 40Ar∕39Ar ages of rocks from northwestern Iran, southern Azerbaijan and the Kura basin with uncertainties (1σ). (b) Kernel density plot in blue (Vermeesch, 2012) of combined 40Ar∕39Ar ages. The duration and timing of the MECO event is indicated by the pink band. (c) Timescale (created with TS-Creator) with scaled eruptive areas (from Sahandi et al., 2014), color legend for areas is the same as in Fig. 1a. Also plotted are radiometric ages from the literature with associated 1σ errors, sorted by age. y axis is an arbitrary unit. Red markers represent extrusive ages; blue markers represent intrusive ages.


Out of the eight prepared samples, seven gave results. Our new 40Ar∕39Ar ages from igneous rocks and ash layers fall within a range of ∼36–45 Ma (Fig. 2a), with weighted mean ages per sample between 39.3–43.1 Ma (Fig. 2b). Detailed results per sample are described in Supplement S5, and detailed results per experiment can be found in Supplement S6–S32. Multiple aliquots of the same samples were measured. The integrated density distribution of these data reveals a peak at around 40.0 Ma (Fig. 2b). All compiled ages are shown together with the scaled areal extent of mapped units of Sahandi et al. (2014) (see Fig. 2c).

2.2 Compilation of literature data

We combined our newly acquired data with more than 420 ages from 72 published studies, including K–Ar, Ar−Ar, U−Pb, Rb−Sr and Re−Os ages (but mainly Ar−Ar and U−Pb; see Supplement S2 and S3). Our age compilation aimed at pre-Quaternary rocks and is incomplete with respect to Quaternary and pre-Paleogene igneous rocks in Iran. We then used a kernel density plot (Vermeesch, 2012) to integrate all ages from 60–0 Ma, together with our newly acquired data. Ages and their 1σ uncertainties are used as input in the calculation of these distributions. Optimal bandwidth is calculated automatically, and we have set the bin width to 1 Myr. When studies did not report the significance level of their uncertainties, we assumed a 1σ uncertainty. Where possible, Ar−Ar ages were recalibrated to the standard of the Fish Canyon Tuff according to the Kuiper et al. (2008) calibration model. In some cases, original studies did not provide sufficient information for recalibration and then the original ages were used. All details of literature ages and associated references are added in Supplement S2 and S3.

Figure 3Radioisotope ages from 0–60 Ma, compiled from the literature, combined with our newly obtained ages. Black thin line represents the probability density plot (PDP); blue filled line represents the kernel density estimate (KDE); boxes represent histograms (numbers on y axis). The timing and duration of the MECO is indicated by the pink box. (a) 40Ar∕39Ar ages from the literature combined with newly obtained ages for extrusive rocks only. (b) U−Pb ages from the literature for extrusive rocks only. (c) 40Ar∕39Ar data and U−Pb data from extrusive rocks only. (d) Combined ages from the literature for intrusive rocks only.


The compilation of 40Ar∕39Ar ages from the literature, mostly from extrusive rocks (only five Ar−Ar ages are from intrusive rocks), yields a highly similar age density distribution to our dated samples (see Fig. 3a), showing a peak at 39.7 Ma. Published U−Pb ages are typically obtained from zircons, which provide less accuracy for eruption ages than 40Ar∕39Ar ages from groundmass, plagioclase, sanidine or biotite (Simon et al., 2008), which is reflected in the greater width of the peaks from extrusive U−Pb ages (see Fig. 3b). Combined, the Ar−Ar and U−Pb ages obtained from extrusive rocks record a wide peak around 42 Ma, with two sub-peaks at 43.4 and 39.4 Ma. Two smaller peaks at 29.8 and 17.1 Ma are apparent (see Fig. 3c). Intrusive activity also peaks around the same time, with radiometric ages from intrusive rocks (n=201) showing a peak at 40.5 Ma, with another sub-peak at 36.6 Ma (Fig. 3d). Smaller peaks in intrusive activity are present at 29.7 and 19.9 Ma.

3 Neotethys volcanism and the MECO

Considering that the Neotethys subduction zone has been active since the Late Triassic (Arvin et al., 2007), our compilation shows a remarkable clustering of ages during the middle Eocene at ∼40Ma. Estimation of the areal extent of middle Eocene volcanic rocks is done using the shapefiles of Sahandi et al. (2014), who made a compilation of geologic maps. According to the geologic maps, 54 % of all area covered by volcanic rocks in Iran is of Eocene age. For the Eocene, shapefiles are classified as “Eocene”, “Eocene–Oligocene”, “Late Eocene–Oligocene”, “Middle Eocene” and “Middle–Late Eocene”. More than half are marked as “Eocene” and not specified further, but of the rest that are specified, almost half are “Middle Eocene”. Assuming that the unspecified Eocene rocks have approximately the same age distribution as the specified Eocene rocks, we estimate that roughly half of the Eocene volcanic rocks in Iran and a quarter of the total area covered by volcanic rocks in Iran are of middle Eocene age. We use these areas to estimate the volumes of volcanic rocks formed in the middle Eocene. We thus assumed that shapefiles specified as “Eocene” had the same proportion of middle Eocene igneous rocks and calculated an areal extent of 38 223 km2 of middle Eocene igneous rocks.

Our compilation indicates that many volcanic provinces in Iran were active simultaneously around 40 Ma (see Fig. 2c), including the Azerbaijan–Bazman magmatic arc in the west, the Sabzevar zone in northeastern Iran (Shafaii Moghadam et al., 2015) and the Lut block in the east (Pang et al., 2013). Some of the largest volumes of middle Eocene volcanic rocks are located in the Talesh Mountains, where four out of five exposures with the largest areal extent are mapped (marked in white in Fig. 1a). Almost three-quarters of U−Pb ages (ntotal=329) in Iran are derived from intrusive rocks (nintrusive=239). All ages of the intrusive rocks together reveal a peak at ∼40.5Ma (Fig. 3d), indicating that the peak of middle Eocene volcanism is also close in time to peak intrusive activity.

It is thus clear that the MECO corresponds to a phase of intense volcanism in the studied area. However, the average error (1σ) of the literature-based ages from 20–60 Ma is 585 kyr and thus exceeds the duration of the MECO (500 kyr). Furthermore, the exact ages of the peaks in volcanic activity in Fig. 2 are sensitive to the number of data points included and are thus not particularly robust – the addition of a few new data points may shift the peaks by thousands of years.

4 Volcanic CO2 emissions in Iran and the MECO

The surface area of Iran covered by middle Eocene volcanic rocks is almost 40 000 km2 (Sahandi et al., 2014; Table 2). These volcanic rocks were produced by numerous eruptions throughout the middle Eocene. In the Alborz and central Iran, middle Eocene volcanic formations are reported to be very thick, with estimates ranging from 3–5 km in the Alborz Mountains (Stöcklin, 1974) to 6–12 km locally throughout nearly all of Iran (Berberian and King, 1981). More recent estimates of the thickness are 3–9 km (e.g., Morley et al., 2009; Verdel et al., 2011). These estimates are supported by geologic maps and their descriptions that are based on extensive fieldwork. Estimates from maps range mostly between 2 and 7 km. On the lower side are for example Saein Qaleh (Kholghi Khasraghi, 1994), Saveh (Ghalamghash et al., 1998a) and Kuhpayeh (Radfar et al., 2002), with thicknesses of ∼2km, Tafresh (Hadjian et al., 1999) with ∼3km, then Meyamey (Amini Chehragh and Ghalamghash, 2021), Tarom (Hirayama et al., 1966) and Kalateh (Jafarian, 2021) with around ∼4km, while Kajan (Amini and Amini Chehragh, 2001), Kahak (Ghalamghash et al., 1998b) and Lahrud (Babakhani et al., 1991) display thicknesses of the Eocene volcanic succession of approximately 6 km, and Bardsir (Mohajjel Kafshdouz and Khodabandeh, 1992) of around 7 km. On the other hand, Iwao and Hushmand-Zadeh (1971) show a generalized lithostratigraphic log of the Karaj formation and mention that the succession reaches a thickness of more than 10 km in the Alborz Mountains. In Table 2, we calculate how much CO2 could have been released through formation of different volumes of volcanic rocks. We calculate this for a range of thicknesses between 2 and 10 km. Extrapolating these thicknesses, this implies a total volume of middle Eocene volcanic rocks between 7.6×104 and 3.8×105km3 (see Table 2) that potentially produced significant amounts of CO2. Our estimates of CO2 released due to middle Eocene volcanism in Iran are likely underestimates, as there is volcanism in other regions along the Neotethys subduction zone. Unfortunately, the lack of shapefiles of Eocene volcanic and intrusive rocks in Armenia and Azerbaijan, along the Lesser Caucasus Mountains (e.g., Allen and Armstrong, 2008), and plutons and volcanic rocks in Armenia (e.g., Moritz et al., 2016; Sahakyan et al., 2016), hampers calculations on additional CO2 emissions within these regions.

Table 2Estimates for the amount of released carbon for different volumes of middle Eocene volcanic rocks in Iran. Estimates marked in bold are within estimates for the MECO (2000–4000 Pg C; Henehan et al., 2020; Sluijs et al., 2013; van der Ploeg et al., 2018).

Download Print Version | Download XLSX

Due to the absence of quantifications of the relation between the erupted volumes of volcanic rocks and emission of CO2 in continental arcs, we make a comparison with the Deccan traps, for which this relation has been calculated. The Deccan traps have an estimated eruptive volume of volcanic and volcaniclastic rocks of 1.3×106km3 (Jay and Widdowson, 2008), with an associated emission 4.14×1017mol CO2 (Tobin et al., 2017). From different estimates of volume and related CO2 emissions of Tobin et al. (2017), we obtain a linear relation of lava volume (in 106km3)/total CO2 (in 1017 mol) ≈0.31 for the Deccan traps.

CO2 degassing rates for continental arcs may be similar to (Marty and Tolstikhin, 1998) or larger than for continental flood basalts (McKenzie et al., 2016; Wignall et al., 2009). As a conservative starting point, we assume a similar volume vs. emission relationship as the Deccan traps, which implies a minimum estimate for CO2 released from middle Eocene volcanism in Iran between 2.34×1016 and 1.22×1017 mol (see Table 2), which corresponds to 292–1461 Pg C. Moreover, the amount of CO2 released during volcanic episodes has been shown to increase substantially if eruptions occur among carbonate-rich sediments (Lee et al., 2013; Lee and Lackey, 2015). For example, CO2 released from carbonate sediments during the emplacement of the Emeishan large igneous province in the end-Guadalupian was estimated to be 3.6–8.6 times higher than the amount of CO2 released by volcanic outgassing alone (Ganino and Arndt, 2009). Indeed, the Eocene volcanism in Iran erupted in shallow marine basins and through significant amounts of carbonate-rich rocks of Jurassic, Cretaceous, and Paleogene age (e.g., Berberian and King, 1981). Glaus (1965) mentions that middle Eocene limestones occur as lenticular masses within the basaltic flows or as consistent horizons associated with tuffs. Verdel (2009) shows that Eocene volcanic rocks are formed in close association with Eocene limestones in northern, western and eastern Iran. This is also the case in central Iran, which can be seen from geologic maps, such as the one from Qom (Emami, 1981). As a result, carbon release associated with the production of volcanic rocks in Iran could be much larger, potentially ranging from 1052 to 12 565 Pg C (see Table 2). This range of CO2 emissions is compatible with the carbon cycle imbalance that drives the MECO in simple carbon cycle simulations constrained by available proxy data (roughly 2000–4000 Pg C; Henehan et al., 2020; Sluijs et al., 2013; van der Ploeg et al., 2018). Table 2 shows that middle Eocene volcanic rocks with thicknesses between 2 and 7 km, and a contribution from limestones, give estimates that lie within the range expected for the MECO (marked in bold in Table 2). There could have been a contribution to CO2 through skarn formation by intrusive activity, which clusters around 40.5 Ma (see Fig. 3d), although we currently lack the constraints to quantitatively assess this. Erosion has affected the entire Iranian plateau and could have eroded away significant volumes of Eocene volcanic rocks. Morley et al. (2009) and Ballato et al. (2011) note that clasts in the Lower and Upper Red Formation (Oligocene–Miocene age), which in many places overlie Eocene volcanics, are for a large part made up of eroded Eocene volcanic rocks. Original thicknesses of Eocene volcanic rocks in Iran could thus have been larger, making our CO2 output estimate a minimum estimate. Despite the fact that sampling biases (i.e., sampling is often focused on easily accessible sites and certain time periods) can never be avoided, our compilation of radiometric ages shows a good correlation to the geologic maps, in the sense that the radiometric ages confirm that the flare-up took place during the middle Eocene. We note that the Miocene peak (Fig. 3c) is relatively high compared to the Eocene, which could be caused by a sampling bias, as the geologic maps (Sahandi et al., 2014) indicate that only 2 %–4 % of Iran is covered by Miocene volcanic rocks.

5 Future perspectives

There are several obstacles in solidifying the link between warming during the MECO and volcanism in the Neotethys subduction zone. First of all, continental arcs are generally active for (tens of) millions of years, while the MECO has a duration of 500 kyr. Moreover, this duration is shorter than common uncertainties for radiometric ages in the Eocene, complicating the establishment of a causal relationship. This is important because a driver for the MECO requires excess CO2 input only during the ∼500kyr spanning the MECO, and not during the time surrounding it (Sluijs et al., 2013). This is also supported by the drop in global ocean osmium isotope ratios, which is specifically associated with the MECO interval (van der Ploeg et al., 2018). Secondly, Iran is a relatively understudied area compared to other (continental) arcs. As a result of this, the number of radiometric ages is low, with on average about one radiometric age for every several hundred km2 of outcrop.

Therefore, the relation in time between the MECO and Neotethys arc flare-up calls for the development of much better age constraints of the volcanic deposits in Iran, and this is certainly feasible. While most flare-ups have to be studied via their intrusive roots, as the extrusive record is removed through erosion (Ducea and Barton, 2007; De Silva et al., 2015), the extrusive record in Iran is extensive so that the ages can be mapped in high detail. Acquisition of radiometric ages throughout sections that cover the entire Eocene volcanic succession could aid in quantification of magmatic flux over time. Moreover, the respective roles of intrusive and extrusive rocks can be assessed to estimate the amount of volatiles of the igneous rocks, and sedimentological studies can provide minimum estimates on how much extrusive rock has been lost through erosion. Studies that focus on the interaction between carbonates and magma chambers could aid in quantifying the carbonate contribution to CO2 release. This would help constrain CO2 input rates across from the Neotethys flare-up to a narrower interval around the MECO.

6 Conclusions

We provide new Ar−Ar ages from volcanic rocks of the Azerbaijan–Bazman Arc in Iran and combine these with literature data to show that a flare-up of continental arc volcanism in Iran peaked about 40 Ma ago, conspicuously close to the Middle Eocene Climatic Optimum. We estimated volumes of middle Eocene volcanism in Iran to be between 7.6×104 and 3.8×105km3. We compared the volume of middle Eocene volcanics in Iran to that of the Deccan traps and estimate that between 292 and 1461 Pg of carbon in the shape of CO2 was released during deposition. Taking into account the fact that all volcanism occurred in shallow marine basins and erupted in and through pre-existing carbonate-rich rocks, CO2 release might have been between 1052 and 12 565 Pg. Thicknesses of the middle Eocene volcanic succession between 2 and 7 km, with a contribution from carbonate-rich rocks, result in estimates of released carbon that are in line with estimates for the MECO. Although the flare-up must be dated much better to establish its chronological relation with the MECO in more detail, we consider it a plausible major contributor to greenhouse warming during the MECO.

Data availability

All data are available in the Supplement of this paper.


Examples of scans of thin sections are supplied in Supplement S1. All details of literature ages and associated references are added in Supplement S2 and S3. S4 is a KMZ file that contains the GPS locations of the literature ages (except for Shafaii Moghadam et al., 2020, who did not provide GPS locations) and can be opened in Google Earth. A detailed description of Ar−Ar results per sample is provided in Supplement S5. Supplement S6–S32 show the results of the 40Ar∕39Ar geochronology per experiment. Supplement S33 shows an extended version of the literature age plot of Fig. 2c. The supplement related to this article is available online at:

Author contributions

Fieldwork was undertaken by AvdB, MH and WK. AvdB, KFK and MH performed Ar−Ar dating. Data analysis was performed by AvdB, KFK, RvdP, MJC and AS. All authors contributed to scientific discussions and were involved in writing the manuscript.

Competing interests

The authors declare no competing interests.


We thank Roel van Elsas for help with Ar−Ar sample preparation. We thank two anonymous reviewers for their comments on an earlier version of this article and three anonymous reviewers for their comments that have improved this article.

Financial support

This research has been supported by the Netherlands Organization for Scientific Research (grant no. 865.10.011) and the European Research Council (SPANC (grant no. 771497)).

Review statement

This paper was edited by Zhengtang Guo and reviewed by three anonymous referees.


Allen, M. B. and Armstrong, H. A.: Arabia – Eurasia collision and the forcing of mid-Cenozoic global cooling, Palaeogeogr. Palaeocl., 265, 52–58,, 2008. 

Amini, B. and Amini Chehragh, M. R.: Geological Map of Iran 1:100.000 Series Sheet 2555 – Kajan, Geological Survey of Iran, Tehran, 2001. 

Amini Chehragh, M. R. and Ghalamghash, J.: Geological Map of Iran 1:100.000 Series Sheet 7162 – Meyamey, Geological Survey of Iran, Tehran, available at:, last access: 11 January 2021. 

Arvin, M., Pan, Y., Dargahi, S., Malekizadeh, A., and Babaei, A.: Petrochemistry of the Siah-Kuh granitoid stock southwest of Kerman, Iran: Implications for initiation of Neotethys subduction, J. Asian Earth Sci., 30, 474–489,, 2007. 

Asiabanha, A. and Foden, J.: Post-collisional transition from an extensional volcano-sedimentary basin to a continental arc in the Alborz Ranges, N-Iran, Lithos, 148, 98–111,, 2012. 

Babakhani, A. R., Nazer, N. H., and Amidi, M.: Geological Map of Iran 1:100.000 Series Sheet 5567 – Lahrud, (5567), Geological Survey of Iran, Tehran, 1991. 

Ballato, P., Uba, C. E., Landgraf, A., Strecker, M. R., Sudo, M., Stockli, D. F., Friedrich, A., and Tabatabaei, S. H.: Arabia-Eurasia continental collision: Insights from late Tertiary foreland-basin evolution in the Alborz Mountains, northern Iran, Geol. Soc. Am. Bull., 123, 106–131,, 2011. 

Beniamovski, V. N., Alekseev, A. S., Ovechkina, M. N., and Oberhänsli, H.: Middle to upper Eocene dysoxic-anoxic Kuma Formation (northeast Peri-Tethys): Biostratigraphy and paleoenvironments, Geol. Soc. Am. Spec. Pap., 369, 95–112, 2003. 

Berberian, M. and King, G. C. P. C. P.: Towards a paleogeography and tectonic evolution of Iran, Can. J. Earth Sci., 18, 1764–1766,, 1981. 

Bijl, P. K., Houben, A. J. P. P., Schouten, S., Bohaty, S. M., Sluijs, A., Reichart, G.-J. J., Sinninghe Damsté, J. S., Brinkhuis, H., Damsté, J. S. S., and Brinkhuis, H.: Transient middle eocene atmospheric CO2 and temperature variations, Science, 330, 819–821,, 2010. 

Bohaty, S. M. and Zachos, J. C.: Significant Southern Ocean warming event in the late middle Eocene, Geology, 31, 1017–1020, 2003. 

Bohaty, S. M., Zachos, J. C., Florindo, F., and Delaney, M. L.: Coupled greenhouse warming and deep-sea acidification in the middle Eocene, Paleoceanography, 24, 1–16,, 2009. 

Boscolo Galazzo, F., Giusberti, L., Luciani, V., and Thomas, E.: Paleoenvironmental changes during the Middle Eocene Climatic Optimum (MECO) and its aftermath: The benthic foraminiferal record from the Alano section (NE Italy), Palaeogeogr. Palaeocl., 378, 22–35,, 2013. 

Boscolo Galazzo, F., Thomas, E., Pagani, M., Warren, C., Luciani, V., and Giusberti, L.: The middle Eocene climatic optimum (MECO): A multiproxy record of paleoceanographic changes in the southeast Atlantic (ODP Site 1263, Walvis Ridge), Paleoceanography, 29, 1143–1161,, 2014. 

Cramwinckel, M. J., Huber, M., Kocken, I. J., Agnini, C., Bijl, P. K., Bohaty, S. M., Frieling, J., Goldner, A., Hilgen, F. J., Kip, E. L., Peterse, F., Van Der Ploeg, R., Röhl, U., Schouten, S., and Sluijs, A.: Synchronous tropical and polar temperature evolution in the Eocene, Nature, 559, 382–386,, 2018. 

Cramwinckel, M. J., van der Ploeg, R., Bijl, P. K., Peterse, F., Bohaty, S. M., Röhl, U., Schouten, S., Middelburg, J. J., and Sluijs, A.: Harmful algae and export production collapse in the equatorial Atlantic during the zenith of Middle Eocene Climatic Optimum warmth, Geology, 47, 247–250, 2019. 

Davoudzadeh, M., Lammerer, B., and Weber-Diefenbach, K.: Paleogeography, Stratigraphy and tectonics of the Tertiary of Iran, Neues Jahrb. fur Geol. und Palaontologie, Abhandlungen 205, 33–67, 1997. 

De Silva, S. L., Riggs, N. R., and Barth, A. P.: Quickening the pulse: Fractal tempos in continental arc magmatism, Elements, 11, 113–118,, 2015. 

Ducea, M. N. and Barton, M. D.: Igniting flare-up events in Cordilleran arcs, Geology, 35, 1047–1050,, 2007. 

Ducea, M. N., Paterson, S. R., and DeCelles, P. G.: High-volume magmatic events in subduction systems, Elements, 11, 99–104,, 2015. 

Edgar, K. M., Bohaty, S. M., Gibbs, S. J., Sexton, P. F., Norris, R. D., and Wilson, P. A.: Symbiont “bleaching” in planktic foraminifera during the Middle Eocene Climatic Optimum, Geology, 41, 15–18,, 2013. 

Emami, M. H.: Geological Quadrangle Map of Iran, 1:250.000 scale, Sheet E6, Qom, Geological Survey of Iran, Tehran, 1981. 

Ganino, C. and Arndt, N. T.: Climate changes caused by degassing of sediments during the emplacement of large igneous provinces, Geology, 37, 323–326,, 2009. 

Ghalamghash, J., Fonoudi, M., and Mehrpartou, M.: Geological Map of Iran 1:100.000 Series Sheet 6060 – Saveh, Geological Survey of Iran, Tehran, 1998a. 

Ghalamghash, J., Babakhani, A. R., Bahroudi, A., and Fonoudi, M.: Geological Map of Iran 1:100.000 Series Sheet 6158 – Kahak, Geological Survey of Iran, Tehran, 1998b. 

Glaus, M.: Die Geologie des Gebietes nordlich des Kandevan-Passes (Zentral-Elburz), Iran, PhD thesis, Nr. 3768, ETH, Zurich, 1965. 

Hadjian, J., Amini, B., and Amini Chehragh, M. R.: Geological Map of Iran 1:100.000 Series Sheet 6059 – Tafreshm, Geological Survey of Iran, Tehran, 1999. 

Henehan, M. J., Edgar, K. M., Foster, G. L., Penman, D. E., Hull, P. M., Greenop, R., Anagnostou, E., and Pearson, P. N.: Revisiting the Middle Eocene Climatic Optimum “Carbon Cycle Conundrum” with new estimates of atmospheric pCO2 from boron isotopes, Paleoceanogr. Paleocl., 35: e2019PA003713,, 2020. 

Hirayama, K., Samimi, M., Zahedi, M., and Hushmand-Zadeh, A.: Geology of the Tarom district, western part (Zanjan area, northwest Iran), Geological Survey of Iran, Tehran, 1966. 

Iwao, S. and Hushmand-Zadeh, A.: Stratigraphy and Petrology of the Low-Grade Regionally Metamorphosed rocks of the Eocene Formation in the Alborz Range, North of Tehran, Iran, J. Japanese Assoc. Mineral. Petrol. Econ. Geol., 65, 265–285, 1971. 

Jafarian, M. B.: Geological Map of Iran 1:100,000 Series Sheet 6860 – Kalateh, Geological Survey of Iran, Tehran, avalable at:, last access: 11 January 2021. 

Jay, A. E. and Widdowson, M.: Stratigraphy, structure and volcanology of the SE Deccan continental flood basalt province: implications for eruptive extent and volumes, J. Geol. Soc. London., 165, 177–188,, 2008. 

Kargaranbafghi, F. and Neubauer, F.: Tectonic forcing to global cooling and aridification at the Eocene-Oligocene transition in the Iranian plateau, Global Planet. Change, 171, 248–254,, 2018. 

Kholghi Khasraghi, M. H.: Geological Map of Iran 1:100.000 Series Sheet 5365, Saein Qaleh, Geological Survey of Iran, Tehran, 1994. 

Kuiper, K. F., Deino, A., Hilgen, F. J., Krijgsman, W., Renne, P. R., and Wijbrans, J. R.: Synchronizing rock clocks of earth history, Science, 320, 500–504,, 2008. 

Lee, C. A. and Lackey, J. S.: Global Continental Arc Flare-ups and Their Relation to Long-Term Greenhouse Conditions, Elements, 11, 125–130,, 2015. 

Lee, C.-T. A., Shen, B., Slotnick, B. S., Liao, K., Dickens, G. R., Yokoyama, Y., Lenardic, A., Dasgupta, R., Jellinek, M., Lackey, J. S., Schneider, T., and Tice, M. M.: Continental arc – island arc fluctuations, growth of crustal carbonates, and long-term climate change, Geosphere, 9, 21–36,, 2013. 

Marty, B. and Tolstikhin, I. N.: CO2 fluxes from mid-ocean ridges, arcs and plumes, Chem. Geol., 145, 233–248, 1998. 

McKenzie, N. R., Horton, B. K., Loomis, S. E., Stockli, D. F., Planavsky, N. J., and Lee, C. A. C.-T. A. C. A. C.-T. A. C. A.: Continental arc volcanism as the principal driver of icehouse-greenhouse variability, Science, 352, 444–447,, 2016. 

Min, K., Mundil, R., Renne, P. R., and Ludwig, K. R.: A test for systematic errors in 40Ar39Ar geochronology through comparison with U/Pb analysis of a 1.1-Ga rhyolite, Geochim. Cosmochim. Acta, 64, 73–98,, 2000. 

Moghadam, H. S., Li, X.-H., Ling, X.-X., Santos, J. F., Stern, R. J., Li, Q.-L., and Ghorbani, G.: Eocene Kashmar granitoids (NE Iran): Petrogenetic constraints from U-Pb zircon geochronology and isotope geochemistry, Lithos, 216–217, 118–135,, 2015. 

Moghadam, H. S., Rossetti, F., Lucci, F., Chiaradia, M., Gerdes, A., Martinez, M. L., Ghorbani, G., and Nasrabady, M.: The calc-alkaline and adakitic volcanism of the Sabzevar structural zone (NE Iran): Implications for the Eocene magmatic flare-up in Central Iran, Lithos, 248–251, 517–535,, 2016. 

Mohajjel Kafshdouz, M. and Khodabandeh, A. A.: Geological Map of Iran 1:100.000 Series Sheet 7335 – Bardsir, Geological Survey of Iran, Tehran, 1992. 

Monster, M.: Multi-method palaeointensity data of the geomagnetic field during the past 500 kyrs from European volcanoes, UU Dept. of Earth Sciences, Utrecht, 2016. 

Moritz, R., Rezeau, H., Ovtcharova, M., Tayan, R., Melkonyan, R., Hovakimyan, S., Ramazanov, V., Selby, D., Ulianov, A., Chiaradia, M., and Putlitz, B.: Long-lived, stationary magmatism and pulsed porphyry systems during Tethyan subduction to post-collision evolution in the southernmost Lesser Caucasus, Armenia and Nakhitchevan, Gondwana Res., 37, 465–503,, 2016. 

Morley, C. K., Kongwung, B., Julapour, A. A. A., Abdolghafourian, M., Hajian, M., Waples, D., Warren, J., Otterdoom, H., Srisuriyon, K., and Kazemi, H.: Structural development of a major late Cenozoic basin and transpressional belt in central Iran: The Central Basin in the Qom-Saveh area, Geosphere, 5, 325–362,, 2009. 

Pang, K.-N., Chung, S.-L., Zarrinkoub, M. H., Khatib, M. M., Mohammadi, S. S., Chiu, H.-Y., Chu, C.-H., Lee, H.-Y., and Lo, C.-H.: Eocene–Oligocene post-collisional magmatism in the Lut–Sistan region, eastern Iran: Magma genesis and tectonic implications, Lithos, 180–181, 234–251,, 2013. 

Radfar, J., Kohansal, R., and S. Zolfagh: Geological Map of Iran 1:100.000 Series Sheet 6455 – Kuhpayeh, Geological Survey of Iran, Tehran, 2002. 

Sahakyan, L., Bosch, D., Sosson, M., Avagyan, A., Galoyan, G. H., Rolland, Y., Bruguier, O., Stepanyan, Z. H., Galland, B., Vardanyan, S., and Bataillon, P. E.: Geochemistry of the Eocene magmatic rocks from the Lesser Caucasus area (Armenia): evidence of a subduction geodynamic environment, Geol. Soc. Spec. Publ., 428, 73–98,, 2016. 

Sahandi, R., Soheili, M., Sadeghi, M., Delavar, T., and Jafari Rad, A.: Compiled geological map of Iran, scale 1:1 000 000, digitally published by the Geological Survey of Iran, 2014. 

Seidov, A. G. and Alizade, K. A.: The formation and mineralogy of bentonites in Azerbaijan, Clay Miner., 6, 157–166, 1966. 

Shafaii Moghadam, H., Li, X.-H., Ling, X.-X., Santos, J. F., Stern, R. J., Li, Q.-L., and Ghorbani, G.: Eocene Kashmar granitoids (NE Iran): Petrogenetic constraints from U−Pb zircon geochronology and isotope geochemistry, Lithos, 216–217, 118–135,, 2015. 

Shafaii Moghadam, H., Li, Q. L., Li, X. H., Stern, R. J., Levresse, G., Santos, J. F., Lopez Martinez, M., Ducea, M. N., Ghorbani, G., and Hassannezhad, A.: Neotethyan Subduction Ignited the Iran Arc and Backarc Differently, J. Geophys. Res.-Sol. Ea., 125, 1–30,, 2020. 

Sieber, N.: Zur Geologie des Gebietes südlich des Taleghan-Tales Zentral-Elburz (Iran), PhD thesis, Nr. 4482, ETH Zurich,, 1970. 

Simon, J. I., Renne, P. R., and Mundil, R.: Implications of pre-eruptive magmatic histories of zircons for U−Pb geochronology of silicic extrusions, Earth Planet. Sc. Lett., 266, 182–194,, 2008. 

Sluijs, A., Zeebe, R. E., Bijl, P. K., and Bohaty, S. M.: A middle Eocene carbon cycle conundrum, Nat. Geosci., 6, 429–434,, 2013. 

Stöcklin, J.: Northern Iran: Alborz Mountains, Geol. Soc. Spec. Publ., 4, 213–234,, 1974. 

Tobin, T. S., Bitz, C. M., and Archer, D.: Modeling climatic effects of carbon dioxide emissions from Deccan Traps volcanic eruptions around the Cretaceous – Paleogene boundary, Palaeogeogr. Palaeocl., 478, 139–148,, 2017.  

van der Boon, A.: From Peri-Tethys to Paratethys: Basin restriction and anoxia in central Eurasia linked to volcanic belts in Iran, Utrecht University, the Netherlands, 2017. 

van der Boon, A., Kuiper, K. F., Villa, G., Renema, W., Meijers, M. J. M. M., Langereis, C. G., Aliyeva, E., and Krijgsman, W.: Onset of Maikop sedimentation and cessation of Eocene arc volcanism in the Talysh Mountains, Azerbaijan, Geol. Soc. Spec. Publ., 428, 145–169,, 2017. 

van der Ploeg, R., Selby, D., Cramwinckel, M. J., Li, Y., Bohaty, S. M., Middelburg, J. J., and Sluijs, A.: Middle Eocene greenhouse warming facilitated by diminished weathering feedback, Nat. Commun., 9, 2877,, 2018. 

Verdel, C.: I. Cenozoic geology of Iran: an integrated study of extensional tectonics and related volcanism II. Ediacaran stratigraphy of the North American Cordillera: new observations from Eastern California and Northern Utah, Dissertation (PhD), California Institute of Technology,, 2009. 

Verdel, C., Wernicke, B. P., Hassanzadeh, J., and Guest, B.: A Paleogene extensional arc flare-up in Iran, Tectonics, 30, 1–20,, 2011. 

Vermeesch, P.: On the visualisation of detrital age distributions, Chem. Geol., 312–313, 190–194,, 2012. 

Vincent, S. J., Allen, M. B., Ismail-Zadeh, A. D., Flecker, R., Foland, K. A., and Simmons, M. D.: Insights from the Talysh of Azerbaijan into the Paleogene evolution of the South Caspian region, Geol. Soc. Am. Bull., 117, 1513–1533,, 2005. 

Wignall, P. B., Sun, Y., Bond, D. P. G. G., Izon, G., Newton, R. J., Védrine, S., Widdowson, M., Ali, J. R., Lai, X., Jiang, H., Cope, H., and Bottrell, S. H.: Volcanism, Mass Extinction, and Carbon Isotope Fluctuations in the Middle Permian of China, Science, 324, 1179–1182,, 2009. 

Witkowski, J., Bohaty, S. M., McCartney, K., and Harwood, D. M.: Enhanced siliceous plankton productivity in response to middle Eocene warming at Southern Ocean ODP Sites 748 and 749, Palaeogeogr. Palaeocl., 326–328, 78–94,, 2012. 

Short summary
40.5 million years ago, Earth's climate warmed, but it is unknown why. Enhanced volcanism has been suggested, but this has not yet been tied to a specific region. We explore an increase in volcanism in Iran. We dated igneous rocks and compiled ages from the literature. We estimated the volume of igneous rocks in Iran in order to calculate the amount of CO2 that could have been released due to enhanced volcanism. We conclude that an increase in volcanism in Iran is a plausible cause of warming.