Middle Miocene climate of southwestern Anatolia from multiple botanical proxies

The middle Miocene climate transition (MMCT) was a phase of global cooling possibly linked to decreasing levels of atmospheric CO₂. The MMCT coincided with the European Mammal Faunal Zone MN6. From this time, important biogeographic links between Anatolia and eastern Africa include the hominid Kenyapithecus. Vertebrate fossils suggested mixed open and forested landscapes under (sub)tropical seasonal climates for Anatolia. Here, we infer the palaeoclimate during the MMCT and the succeeding cooling phase for a middle Miocene (14.8–13.2 Ma) intramontane basin in southwestern Anatolia using three palaeobotanical proxies: (i) Köppen signatures based on the nearest living-relative principle; (ii) leaf physiognomy analysed with the Climate Leaf Analysis Multivariate Program (CLAMP); (iii) genus-level biogeographic affinities of fossil flora with modern regions. The three proxies reject tropical and hot subtropical climates for the MMCT of southwestern Anatolia and instead infer mild warm temperate C climates. Köppen signatures reject summer-dry Cs climates but cannot discriminate between fully humid Cf and winter-dry Cw; CLAMP reconstructs Cf climate based on the low X3.wet∕X3.dry ratio. Additionally, we assess whether the palaeobotanical record resolves transitions from the warm Miocene Climatic Optimum (MCO, 16.8–14.7 Ma) to the MMCT (14.7–13.9 Ma), and a more pronounced cooling at 13.9–13.8 Ma, as reconstructed from benthic stable isotope data. For southwestern Anatolia, we find that arboreal taxa predominate in MCO flora (MN5), whereas in MMCT flora (MN6) abundances of arboreal and non-arboreal elements strongly fluctuate, indicating higher structural complexity of the vegetation. Our data show a distinct pollen zone between MN6 and MN7+8 dominated by herbaceous taxa. The boundary between MN6 and MN7+8, roughly corresponding to a first abrupt cooling at 13.9–13.8 Ma, might be associated with this herb-rich pollen zone.

1 Introduction

The middle Miocene (15.97–11.63 Ma, ICS-chart 2017/02, Cohen, 2013) is characterized by a warm phase lasting until ca. 15 Ma that was followed by a gradual cooling and the restoration of a major Antarctic ice sheet and the first northern hemispheric glaciations (Holbourn et al., 2014). It has been suggested that the final closure of the Mediterranean gateway connecting the Mediterranean with the Indian Ocean and the resulting changes in ocean circulation might have been one of the reasons for the final expansion of the East Antarctic Ice Sheet around 14.8 Ma (Flower and Kennett, 1995). During the middle Miocene climate transition (MMCT) at 14.7 to 13.8 Ma, a drop in sea surface temperatures of 6–7 ^∘C occurred (Shevenell et al., 2004). At the same time, different proxies to reconstruct atmospheric CO₂ levels for the Miocene Climatic Optimum (MCO), MMCT, and the succeeding more pronounced cooling do not concur (Beerling and Royer, 2011). Specifically, stable isotope data from phytoplankton infer stable CO₂ levels for the Neogene, with minor fluctuations (MCO, 227–327 ppm; MMCT, 265–300 ppm; see Table S1 of Beerling and Royer, 2011), while stomata densities from fossil leaves suggest a pronounced drop in CO₂ after the MCO (see Table S1 of Beerling and Royer, 2011).

The European Mammal Faunal Zone MN6 (14.8–13.8 Ma; Neubauer et al., 2015) coincides with the MMCT. From this period world-famous vertebrate faunas are known from western Anatolia (e.g. Andrews and Tobien, 1977; Mayda et al., 2015) including the hominoids Griphopithecus alpani in Çandır and Paşalar and Kenyapithecus kizili in Paşalar (Stringer and Andrews, 2011). Geraads et al. (2003) investigated the depositional environment and large mammal fauna of Çandır close to Ankara and inferred open landscapes for this locality. Bernor et al. (1979, p. 86) analysed the community structure of Turkish and European middle Miocene fauna and suggested that “faunas adapted to woodland conditions were present … at localities such as Paşalar and Yeni Eskihisar [MN7 + 8]” while the “Çandır fauna has a community structure more suggestive of closed woodland conditions”. This interpretation is the exact opposite of that by Geraads et al. (2003). Recent investigations using carnivore guild structure suggest a “mixed environment between tropical forest and open savannah landscapes” for Çandır and Paşalar (Mayda et al., 2015). Strömberg et al. (2007) investigated phytoliths (plant silica bodies) from early to late Miocene deposits of Turkey and suggested that open, grass-dominated habitats had become common in Turkey and adjacent areas by the early Miocene (ca. 20 Ma). In contrast, Kayseri-Özer (2017) using “integrated plant record” (IPR) analysis (Kovar-Eder et al., 2008) suggested that various forest types covered most of western and Central Anatolia during the middle Miocene (broad-leaved evergreen and mixed mesophytic forests and ecotones between these forests).

Here we use a large dataset from recently published macrofossils and pollen, spores, and cysts from a well-dated middle Miocene basin in western Anatolia to infer palaeoclimate and palaeoenvironments using three palaeobotanical proxies: climate affinity of modern analogues (“nearest living relatives”; taxon-based approach), leaf physiognomy (a-taxonomic), and biogeographic affiliation of plant communities (also taxon-based). The following research questions are addressed. How do the three approaches resolve local climate conditions of Anatolia during a phase of global cooling, ca. 15–13 million years ago? Do different proxies agree on climate inference? Where do modern climates occur that correspond to middle Miocene climates of western Anatolia? Can the palaeobotanical record resolve transitions between MCO, MMCT, and the succeeding more pronounced cooling during the middle Miocene?

Figure 1 Geographic and regional geologic setting of the Yatağan Basin. (a) Map showing the geographical position of the Yatağan Basin (2) and the MN6 vertebrate fossil localities (*) Paşalar (A) and Çandir (B). (b) Simplified regional geological map of the Yatağan Basin based on Becker-Platen (1970) and Atalay (1980); lignite mines Eskihisar (A), Tınaz (B), Salihpaşalar (C); vertebrate fossil locality (*) Yeni Eskihisar MN7/8 (D).

2 Material and methods

2.1 Geological setting

The Yatağan Basin is a southeast-trending graben (50 km long, 15 km wide) in the province of Muğla, southwestern Turkey (Fig. 1). The Neogene basin fill is up to 600 m thick and divided into the Eskihisar Formation (early to middle Miocene), the Yatağan Formation (late Miocene to early Pliocene), and the Milet Formation (middle to late Pliocene; Alçiçek, 2010). The Eskihisar Formation comprises the Turgut Member (reddened alluvial-fan deposits followed by fluviatile deposits and lignites) and the Sekköy Member (fossiliferous limnic marls and limestones); all economically exploited lignite seams of the Yatağan Basin are confined to the transition zone of these two members (Atalay, 1980; Becker-Platen, 1970).

For the present study, we investigated the palaeobotanical content (pollen and plant macrofossils) of the upper Turgut and the Sekköy members exposed at the lignite mines of Eskihisar, Salihpaşalar, and Tınaz (Fig. 1b). The age of the investigated sediments is well constrained by mammal fossils (Eskihisar lignite gallery locality, MN6, Gomphotherium angustidens Cuvier 1817, Percrocuta miocenica Pavlov et Thenius 1965 – Bouchal et al., 2017; Yeni Eskihisar vertebrate locality, MN 7/8 – The NOW Community, 2018) and by radiometric dates from the upper Sekköy Member (13.2 Ma±0.35, Becker-Platen et al., 1977). Hence, the investigated pollen zones (PZs) 1, 2, and 2/3 and the Yeni Eskihisar pollen assemblage represent the Neogene mammal zones MN6 and MN7 + 8, 14.8–13.2 Ma (Neubauer et al., 2015). The layers from which most of the leaf fossils originate correspond to PZ 2. A ∼20 m section comprised of limestone and clayey limestone between PZ 2/3 and the Yeni Eskihisar assemblage is barren of palynological content (Fig. 2).

Figure 2 Generalized lithostratigraphic column for the Eskihisar lignite mine and pollen zones (PZs). The main part of the investigated plant macrofossils originates from ca. 10 m thick deposits overlying the exploited lignite seams (part of the section highlighted in grey corresponding to PZ 2). Yeni Eskihisar 2 (YE2) and Yeni Eskihisar 1 (YE1) vertebrate fossil localities (Becker-Platen et al., 1977). Radiometrically dated tuff layers (*): 1* – 11.2±0.2 Ma; 2* – 13.2±0.35 Ma (Becker-Platen et al., 1977).

2.2 Plant material

The investigated plant material comprises roughly 1800 macrofossils (mainly leaf fossils) collected between 2010 and 2017. Macrofossils represent 77 taxa, of which 5 belong to gymnosperms and 72 to angiosperms. Pollen, spores, and cysts from five pollen zones (Fig. 2) represent 182 taxa, of which 1 is a fungus, 9 are algae, 17 are moss or fern ally spores, 15 are gymnosperms, and 140 are angiosperms (Supplement S1; for taxonomic descriptions of the plant taxa see Yavuz-Işık et al., 2011; Bouchal et al., 2016, 2017; Bouchal, 2018; Güner et al., 2017).

Figure 3 Köppen signal for genus Tilia extracted from 26 extant species. (a) Köppen–Geiger climates in which Tilia is present. (b) Combined Köppen signature of all 26 extant Tilia species. Diagrams based on extracted data from Supplement S2.

2.3 Köppen signatures

Fossil taxa that are resolved to genus or sectional level were represented by extant members of the genera and sections as modern analogues. First, for accepted taxa (IPNI, http://www.ipni.org/index.html, last access: 13 September 2018; most recent regional floras and monographs), the distribution ranges were determined. Then, 26 Köppen–Geiger climate types (see Table 3 for detailed explanation of Köppen–Geiger climate types, and Kottek et al., 2006; Peel et al., 2007; Rubel et al., 2017; Global_1986-2010_KG_5m.kmz) were mapped onto modern distribution ranges using Google Earth to establish “Köppen signatures” (Denk et al., 2013) for each modern analogue. The representation of different climate types was first scored for each species within a genus as present (1)/absent (0). To summarize preferences for climate types of all modern analogues, we used an implicit weighting scheme to discriminate between modern analogues that are highly decisive (climatically constrained) and those that can be found in many climate zones. The sum of each modern species' Köppen signature is always 1. For example, Tilia chingiana is present in two Köppen–Geiger climate types, Cfa and Cfb, which count as 0.5 for each type, while Tilia americana is present in 10 Köppen–Geiger climate types (As, Aw, Cfa, Cfb, Dfa, Dfb, Cwa, Cwb, BSk, BWh), all counting as 0.1. The Köppen signature of a genus or section, the modern analogue of a fossil taxon, is the sum of its species' Köppen signatures for each climate type divided by the total number of scored species for this genus. By this, the percentage representation of each Köppen–Geiger climate type was determined for a genus/section. In the case of Tilia, the distribution ranges of 26 species resulted in the following genus Köppen signature: Cfa, 22.1 %, Cfb, 14.7 %, Cwa, 19.9 %, Cwb, 9.1 %, Dfb, 5.7 % for the five most common climate types. Figure 3a shows all climate types realized in genus Tilia; Fig. 3b shows that the genus occurs predominantly in Cf and Cw Köppen–Geiger climate types and that tropical and desert climates are nearly absent (see Supplement S3 for genus-level scoring of Köppen–Geiger climate types for all plant taxa encountered in the Yatağan basin fossil assemblages).

Figure 4 Selected pollen grains light microscopy (LM) (a, c, e, g, i, k) and scanning electron microscopy (SEM) (b, d, f, h, j, l) micrographs of the same fossil pollen grain of the Eskihisar (E), Tınaz (T), and Salihpaşalar (S) sections. (a–b) Nitraria sp., EV (E, S153567). (c–d). Sapotaceae gen. indet., EV (T, S143604). (e–f) Decodon sp., EV (S, S153635). (g–h) Fagus sp., PV (T, S143621). (i–j) Cathaya sp., (i) PV, (j) PRV (S, S153632). (k–l) Cedrus sp., EV (E, S153590). EV: equatorial view; PV: polar view; PRV: proximal view. Scale bar: 10 µm (a–l).

For taxa that are resolved to family-level only, mainly pollen taxa of herbaceous and a few woody angiosperm groups, the distributions of extant members of the family were combined into a general family distribution range and the corresponding Köppen–Geiger climate types determined.

Genus-level Köppen–Geiger signals were used to account for possible niche evolution within lineages or species groups of a genus. For example, modern species of Quercus section Ilex are typical members of sclerophyllous, evergreen Mediterranean forest and shrubland vegetation thriving under a Cs (summer-dry warm temperate) climate in western Eurasia and to the south of the eastern Hindu Kush and Karakorum ranges but also occur in humid, mesophytic forests from Afghanistan to East Asia (Cf and Cw climates). To account for this climate niche variability, all species of section Ilex were scored for the general Köppen signature of section Ilex. Hence, the entire section was used as modern analogue, the nearest living relative (NLR), for several fossil species of Quercus section Ilex.

2.4 CLAMP

We inferred quantitative palaeoclimate parameters for the three Yatağan Basin floras using the Climate Leaf Analysis Multivariate Program (CLAMP) (Yang et al., 2011). CLAMP makes use of the relationship between leaf physiognomy of dicotyledonous flowering plants and climate and, hence, is a non-taxonomic approach to palaeoclimate inference (Spicer, 2008). CLAMP calibrates the numerical relations between leaf physiognomy of woody dicots and meteorological parameters in modern terrestrial environments. With this calibration, past climatic data can be determined from leaf fossil assemblages if the sampling of the fossil assemblage represents the characteristics of the living source vegetation well (http://clamp.ibcas.ac.cn, last access: 13 September 2018). Modern and fossil leaf physiognomic data are positioned in multidimensional physiognomic space using canonical correspondence analysis (CANOCO; Ter Braak, 1986). CANOCO orders vegetation sites based on a set of attributes (leaf physiognomic characters).

For modern sites, climate variables are known from long-term observations of climate stations or from high-resolution gridded climate data (New et al., 1999, 2002; Spicer et al., 2009). Vectors for each of the measured climate variables can be positioned in physiognomic space and calibrated. Palaeoclimate variables can then be quantified by scoring a fossil assemblage in the same manner as for the modern vegetation and positioning the fossil site in physiognomic space (http://clamp.ibcas.ac.cn, last access: 13 September 2018).

For the present study, 36 different leaf characters (including leaf shape and size, apex shape, base shape, and leaf margin characteristics) were scored for 61, 63, and 14 dicotyledonous leaf morphotypes from three localities: Tınaz, Eskihisar, and Salihpaşalar (see Supplement S3 for scoring of morphotypes), following the CLAMP protocols (http://clamp.ibcas.ac.cn, last access: 13 September 2018). At genus level, the flora of the Yatağan Basin shows the highest similarity with Eurasian extant woody angiosperms (Table 1), thus the PhysgAsia1 Calibration file dataset of CLAMP was used to position the fossil data.

Table 1 Genus-level biogeographic affinities of fossil taxa of the Yatağan Basin flora.

^a Magri et al. (2017). ^b Corbett and Manchester (2004). ^c Jia et al. (2015). ^d Martinetto (2001). ^e Wang et al. (2007). ^f Including northern Africa. ^g Extinct genus. wEUR: western Eurasia; EA: East Asia; eNA: eastern North America; wNA: western North America; SA: South America; AF: Africa (excluding northern Africa); AUS: Australia.

2.5 Genus-level biogeographic affinities

For all fossil taxa determined to genus level, the present distribution was tabulated indicating the presence or absence of a genus in western Eurasia, East Asia, eastern North America, western North America, or Africa (Table 1).

3 Results

3.1 Climate inference from Köppen signatures (Fig. 5, Supplement S4, S5)

For the fossil plant assemblages warm temperate to temperate C and D climates accounted for almost 80 % of the realized Köppen–Geiger climate types of all taxa in a fossil plant assemblage (using genus-level NLR). The sum of Cf, Df, Cw, and Dw climates amounted to 60 %–70 % in all assemblages (highest scores in macrofossil assemblages).

Figure 5 Köppen signals for the Yatağan Basin flora. (a) Pollen zone (PZ) 1 (MN6; 14.95–13.9 Ma) of the Eskihisar (E), Tınaz (T), and Salihpaşalar (S) localities and the combined signal of all present taxa from PZ 1 of the three Yatağan Basin localities (YB). (b) PZ 2 (MN6) of E, T, S, YB. (c) Macrofossil (MF) assemblages (same level as PZ 2) of E, T, S. (d) PZ 2/3 of T (younger than Yeni Eskihisar vertebrate locality). (e) Yeni Eskihisar vertebrate locality pollen assemblage (MN7/8, younger than radiometric age 13.2 Ma). Diagrams based on extracted data from Supplement S2.

Overall, the best represented Köppen–Geiger climate types when using genus-level NLR were Cfa (warm temperate, fully humid, hot summer), followed by Cfb (warm temperate, fully humid, warm summer; e.g. Fig. 4g, h), Cwa (warm temperate, winter-dry, hot summer), and Cwb (warm temperate, winter-dry, warm summer). Summer-dry Cs climates were represented by 9 %–13 % and arid (generally dry) B climates by 6 %–11 % (e.g. Fig. 4a, b, Table 3, Supplement S4). Tropical (equatorial) climates (A) are represented by 9 %–11 % in older assemblages and 7 %–8 % in the two youngest assemblages (PZ 2/3 and Yeni Eskihisar). Of 1555 modern species used to inform the Köppen signatures of the NLRs for the fossil taxa, 119 show marginal range extensions into Af climate, 168 into Am (heavy monsoon), 85 into As, and 295 into Aw (Supplement S2). Taxa extending into tropical climates are mainly species of Pinus, Celtis, Smilax, and Viburnum, Quercus sections Quercus and Lobatae, Juglandaceae subfamily Engelhardioideae, Oleaceae, and Sapotaceae (Fig. 4c, d). The exclusion of Köppen–Geiger climate signals extracted from cosmopolitan and/or gymnospermous taxa (Fig. 4i–l) did not change the general trends (Supplement S5).

Figure 6 CLAMP climate inference for the macrofossil assemblage of (E) Eskihisar and (T) Tınaz (same level as PZ 2). (a) Mean annual temperature (MAT). (b) Coldest-month mean temperature (CMMT). (c) Precipitation of the 3 wettest months. (d) Precipitation of the 3 driest months.

3.2 CLAMP

Sixty-three morphotypes were scored for Eskihisar (Fig. 6; see Supplement S3 for score sheets and other reconstructed climate parameters). The inferred values for mean annual temperature (MAT) were (11.2–) 12.6 (–14) ^∘C, for coldest-month mean temperature (CMMT) (0.3–) 2.3 (–4.4) ^∘C, and for the three wettest months (X3.wet) (410–) 666 (–936) mm and for the three driest months (X3.dry) (148–) 204 (–262) mm. The ratio X3.wet∕X3.dry was between 2.9 and 3.6. For Tınaz, the reconstructed MAT was (12.3–) 13.8 (–15.2) ^∘C, CMMT (1.5–) 3.6 (–5.6) ^∘C, X3.wet (420–) 700 (–980) mm, and X3.dry (146–) 205 (–260) mm. The ratio X3.wet∕X3.dry was between 2.9 and 3.8. Values for Salihpaşalar are not considered here as they are based on too small a set of morphotypes (see Supplement S3).

3.3 Genus-level biogeography

The genus-level biogeographic analysis of the four Yatağan Basin floras ranging in age from 14.8 to 13.2 Ma (MN6 into MN7 + 8; Table 1) shows that closest biogeographic relationships are with the modern East Asian flora (54 of 59 taxa shared with East Asia); 48 and 44 genera are shared with the modern western Eurasian and eastern North American flora, respectively. Among modern tropical floras, the closest relationships are with South America (21), followed by Africa (16) and northern–northeastern Australia (13). Most taxa extending to tropical regions are cosmopolitan (e.g. Euphorbia, Drosera, Phragmites) and hence of little discriminative power. This is also true for higher taxa such as Polygalaceae and Valerianoideae. The fossil species Smilax miohavanensis belongs to a subtropical–tropical clade of extant species (Denk et al., 2015) and is the only member of this group in Eurasia; it has its last occurrence in the middle Miocene flora of the Yatağan Basin. Overall, the dominating biogeographic signal is a northern hemispheric one.

Table 2 Arboreal to non-arboreal pollen ratios in southwestern Anatolia across the MCO, MMCT, and subsequent cooling phase.

^* Perhaps linked with the 13.9–13.8 Ma cooling event (Holbourn et al., 2014).
AP: arboreal pollen (angiosperms).
NAP: non-arboreal pollen (angiosperms).
Wavy line: profiles separated by tens of metres of sediment barren of pollen.
Red to blue indicates change from warmer to cooler climate.

Table 3 Description of Köppen–Geiger climate symbols and defining criteria (Kottek et al., 2006; Peel et al., 2007). MAP: mean annual precipitation; MAT: mean annual temperature; T_hot: temperature of the hottest month; T_cold: temperature of the coldest month; T_mon10: number of months where the temperature is above 10 ^∘C; P_dry: precipitation of the driest month; P_sdry: precipitation of the driest month in summer; P_wdry: precipitation of the driest month in winter; P_swet: precipitation of the wettest month in summer; P_wwet: precipitation of the wettest month in winter; P_threshold: varies according to the following rules – if 70 % of MAP occurs in winter, then $P_{\mathrm{threshold}}=\mathrm{2}\times \text{MAT}$; if 70 % of MAP occurs in summer, then $P_{\mathrm{threshold}}=\mathrm{2}\times \text{MAT}+\mathrm{28}{}^{\circ }C$; otherwise $P_{\mathrm{threshold}}=\mathrm{2}\times \text{MAT}+\mathrm{14}$. Summer (winter) is defined as the warmer (cooler) 6-month period of October–March and April–September.

3.4 Changes in ratios of arboreal to non-arboreal pollen

Ratios of arboreal pollen (AP) to non-arboreal pollen (NAP) change considerably among and within pollen zones of the Yatağan Basin assemblages (Table 2, Supplement S6). Pollen zone 1 (main lignite seam) consistently has high percentages of AP (94 %–70 %). In contrast, AP percentage values fluctuate throughout pollen zone 2, with values from 89 to 29. Pollen zone 2–3, only covered in the Tınaz section, records AP percentages of 50 to 19. Above, the MN7 + 8 assemblage of Yeni Eskihisar shows again a higher proportion of arboreal taxa (67 %). Similarly, from the vertebrate locality Çatakbağyaka (revised age MN7 + 8, 12 km south of the Yatağan Basin), AP percentages range from ca. 50 % to ca. 80 % (Jiménez-Moreno, 2005; Mayda et al., 2016; Bouchal et al., 2017; Aiglstorfer et al., 2018).

We used the threshold ($\text{AP}/\text{NAP}=\mathrm{3.85}$) proposed by Favre et al. (2008) to separate between tree- and herb-prevalent environments. This ratio translates into AP percentages of close to 80 % to predict reliably tree-prevalent landscapes. As can be seen in Supplement S6, pollen zones 1 and 2 are largely dominated by forested environments. In the upper part of PZ2 (Tınaz, Eskihisar), PZ2/3, and PZ3 (Tınaz) herb-prevalent landscapes are inferred. However, it is noteworthy, that although NAP taxa are more abundant in these pollen zones, AP taxa continue to have fairly high percentages as well (Bouchal et al., 2016, 2017). For example, Fagus and Quercus deciduous and evergreen type are still above the threshold values indicative of local tree presence (Lisitsyna et al., 2011). Hence, the opening of the vegetation in the upper parts of PZ2 and in PZ2/3 and PZ3 may actually represent a coexistence of forest and open vegetation.

4 Discussion

4.1 Climate inference using Köppen signatures and CLAMP

Using Köppen signatures, we made a semi-quantitative reconstruction of the palaeoclimate of the Yatağan Basin during the middle Miocene. All Köppen signatures used here rely on the nearest living-relative principle (Denk et al., 2013). Such approaches are prone to error because niche evolution may have occurred in lineages; the morphologically NLRs of fossil taxa, a species or group of morphologically similar species, may have different niches, and the shift is difficult to quantify (Ackerly, 2004; Grimm and Potts, 2016; Denk et al., 2017). Hence, we decided against applying quantitative NLR methods and determined Köppen signatures for fossil taxa using information from all extant species of a genus used as NLR to avoid bias from undetected niche shifts.

It is important to keep climatic niche shift in mind when using NLR based approaches to palaeoclimate inference and interpreting their results (Grimm and Potts, 2016; Denk et al., 2017). In our dataset of 1555 modern species, 295 also occur in tropical Aw climates. Most of them belong to clades (monophyletic sections, genera, families) that occur in a wide range of climate types (e.g. Amaranthaceae, Celtis, white and red oaks). Others, such as Engelhardia are usually interpreted as tropical–subtropical evergreen elements (Kvaček, 2007) based on the distribution range of the extant genera of the comprising subfamily, the Engelhardoideae. However, Engelhardia of the western Eurasian Cenozoic belongs to its own (extinct) section or genus Palaeocarya (Kvaček, 2007) with a stratigraphic range from the Eocene to the Pliocene. Pollen, foliage, and reproductive structures of fossil material clearly belong to the subfamily Engelhardioideae but cannot be assigned to just a single modern genus Engelhardia (tropical southeast Asia). Instead, the fossil taxon is a mosaic taxon having characteristics of both American and Asian members of the subfamily. Kvaček (2007) noted that the fossil genus/subgenus flourished in subtropical climates during the Eocene but in distinctly temperate climates with coldest-month mean temperatures close to the freezing point in the Neogene, in stark contrast to the surviving four, likely relict genera of the Engelhardioideae. Hence, this extinct lineage of Engelhardioideae is not well represented by a single genus or the combination of all extant genera and their constituent species. Similarly, representatives of Smilax havanensis and allied species are part of a New World clade with most species occurring in tropical climates. However, the single Old World member of the clade, the fossil species S. miohavanensis, is known from early to middle Miocene strata of Anatolia and central Europe (Denk et al., 2015). This fossil species formed part of plant assemblages that rule out tropical climates. In this case, inferring palaeoclimate from extant distribution data will inevitably only produce noise in the climatic signal.

Overall, the most common Köppen–Geiger climate types of NLR taxa of the Yatağan flora were warm temperate C types, and among C types fully humid Cf climates were better represented than more seasonal Cw and Cs types (Fig. 5; Supplement S4). Cs types played only a minor role; however, there was no clear preference of Cf over Cw climates in the representation of Köppen–Geiger climate types. Removing azonal taxa or taxa commonly associated with higher elevations (conifers) did not affect the general signal.

In contrast, CLAMP is not based on NLR and hence not potentially biased by taxonomic error. Its combination with the Köppen signature analysis provides a powerful tool for climate inference and to discern between seasonal Cw (winter dry) and Cs (summer dry) and fully humid Cf climates. Specifically, the ratio of the wettest and the driest month clearly distinguishes strongly seasonal summer rain (monsoon) climates (Cw; precipitation wettest month $>\mathrm{10}\times$ precipitation driest month, $P_{\mathrm{wdry}}<P_{\mathrm{swet}}/\mathrm{10}$; Peel et al., 2007) from weakly seasonal, fully humid climates (Cf; precipitation wettest month $\ll \mathrm{10}\times$ precipitation driest month). Precipitation values for X3.wet and X3.dry inferred by CLAMP and the ratio between these ranges being between 2.9 and 3.8 thus largely rules out a Cw climate (X3.wet and X3.dry are closely correlated to P_dry∕wet). In conjunction with the Köppen signature results ruling out summer-dry conditions, the CLAMP precipitation and temperature estimates point towards cold subtropical to mild temperate Cfa climates at the margin to fully temperate Cfb climates.

4.2 Comparison to palaeoclimate and palaeoenvironment inferences from other proxies

A further refinement of previous climate and vegetation inferences can be made regarding the distinction between tropical ($T_{\min }\equiv \text{CMMT}\ge \mathrm{18}$ ^∘C), subtropical (8–12 months with T≥10 ^∘C; ∼MAT 12–18 ^∘C, and CMMT <18 ^∘C), and temperate climates. CLAMP consistently resolves MAT <18 ^∘C and CMMT <6 ^∘C for the localities Eskihisar and Tınaz, and this agrees with the results from Köppen signatures and a previous qualitative assessment of palaeoenvironments in the Yatağan Basin (Güner et al., 2017). Both these results, strongly rejecting seasonal Cw climates, summer-dry Cs, and tropical A climates (at least for non-coastal areas) for the middle Miocene of western Anatolia, have implications for the reconstruction of palaeoenvironments of famous vertebrate localities in Anatolia that are assigned to MN6. The δ¹³C composition from fossil tooth enamel at Paşalar, western Anatolia, MN6, indicates that animals were feeding on C₃ vegetation (Quade and Cerling, 1995). The palaeoenvironment for this locality was determined as closer to Indian subtropical forests, with seasonal summer rainfalls (i.e. warm Cwa climates), semi-deciduous forest, and dense ground vegetation (Stringer and Andrews, 2011; Mayda et al., 2015). Using carnivore guild structures, Morlo et al. (2010) inferred open (Serengeti type, Aw climate) landscapes for the Central Anatolian MN6 vertebrate locality Çandır. Also, the NOW database (http://www.helsinki.fi/science/now/, last access: 13 September 2018; The NOW Community, 2018) refers to Çandır as more open (“woodland biome”, “open vegetation structure”, “grassland with mosaic of forests”) and to Paşalar as more forested landscapes (“subtropical”, “closed vegetation structure”, “semi-deciduous forests”). Bernor et al. (1979) using the community structure of vertebrate fauna inferred densely wooded environments for Çandır. In a later study based on a taxonomic revision of carnivores, Mayda et al. (2015) proposed a mixed environment between tropical forests and open savannah landscapes for Çandır. It is important to note that these carnivore guild structure studies used only two modern calibration faunas to estimate palaeoenvironments: one tropical rainforest fauna in Guyana and one savannah (tropical) fauna in the Serengeti (Morlo et al., 2010). Thus, using this proxy, only two environments can be reconstructed: tropical savannah or rainforest.

Our plant-proxy-based climate reconstruction unambiguously rejects a tropical climate for the middle Miocene Yatağan Basin, and major biogeographic patterns strongly suggest northern hemispheric affinities. Similar environmental conditions as reconstructed in our study have been inferred for most of western Anatolia during the late early and middle Miocene (Kayseri-Özer, 2017). Most proxies currently used to infer climate and vegetation in western Anatolia during the middle Miocene (carnivore guild structures, vertebrate community structure, plant functional types, plant macrofossils, pollen and spores; Mayda et al., 2015, 2016; Kayseri-Özer, 2017; Güner et al., 2017; Bouchal et al., 2016, 2017; Bouchal, 2018) clearly infer forested vegetation with varying contributions of open vegetation. In contrast, Strömberg et al. (2007) found that “all Miocene phytolith assemblages point to relatively open vegetation, such as savanna or open woodland dominated by open-habitat grasses, or a mixture of grassland and wooded areas”. This result may be biased (see Jokela, 2015, p. 44) and increased diversity of grass types in the phytolith record may not necessarily indicate the presence of widespread open, grass-dominated landscapes.

4.3 Modern climate analogues

The inferred climate for the middle Miocene Yatağan Basin plant assemblages is characterized by a MAT of 11–15 ^∘C, coldest-month mean temperature (CMMT) of 0–6 ^∘C, MAP of ca. 1000–2000 mm, and ratios of X3.wet∕X3.dry of 2.9–3.8. A non-exhaustive search for climate stations with this combination of climate parameters (Supplement S7) identified a single closest match: Pacific central Honshu in Japan. X3.wet∕X3.dry ratios and MAT are similar to the upper limits of the ranges reconstructed for the middle Miocene Yatağan Basin. East Asian Cf climates are generally characterized by distinct summer rain maxima. The modern vegetation of Japan is home to many plant taxa that are currently absent from western Eurasia but were abundant in Neogene plant assemblages of western Eurasia (e.g. Cephalotaxus, Cryptomeria, Torreya, Alangium, Camellia, Castanopsis, Cercidiphyllum, Daphniphyllum, Eurya, Fatsia; Mai, 1995; Miyawaki, 1984; see also Milne, 2004). These taxa require warm and humid equable climates.

A further close match is the area from northern Turkey via Georgia to northern Iran: the Euxinian–Hyrcanian region (Supplement S7). Climates at the transition between Csa and Cfa/b of the region north of Istanbul have up to 1166 mm MAP (Ustaoğlu, 2012) and other climate parameters in this area match the Miocene climate of southwestern Turkey inferred by CLAMP. Towards the humid northeastern part of Turkey, X3.wet∕X3.dry ratios are lower (2.4 for Rize, Hopa, and P'ot'i and K'obulet'i in adjacent western Georgia). Further to the east, south of the Caspian Sea, Rasht and Kiashahr have Cfa and borderline Csa to Cfa climates with slightly more pronounced seasonality than the reconstructed climate for the Miocene of southwestern Turkey (X3.wet∕X3.dry ratios of 4.4 and 4.2). In contrast, X3.wet∕X3.dry ratios in modern Mediterranean western and southwestern Turkey amount to 25 (Izmir) and 21.8 (Muğla, Yatağan Basin). It is noteworthy that modern Cf climates of the Euxinian–Hyrcanian region differ markedly from those of the Pacific part of Honshu by their summer minima in rainfall (Supplement S7). This feature indicates a (weak) Mediterranean influence in this region. According to Biltekin et al. (2015), the Anatolian refugium emerged after the retreat of the Paratethys Sea in the Pliocene and increasing monsoon influence (increased summer rainfall) over the northeastern Mediterranean region (the latter accounting for the much higher summer precipitation in the Euxinian–Hyrcanian than in the Mediterranean region). The Mediterranean climate type in Europe appeared first during the late Pliocene and early Pleistocene (ca. 3.2–2.3 Ma; Suc, 1984) coinciding with first large-scale north hemispheric glaciation in the North Atlantic (Denk et al., 2011).

4.4 Detection of Miocene global climatic changes in the terrestrial fossil record

High-resolution benthic stable isotopic data provide a detailed chronology of (global) climatic changes across the MCO, the MMCT, and the subsequent more pronounced cooling (Holbourn et al., 2014). The terrestrial record usually does not provide the same temporal resolution but allows focussing on regional patterns. The transition from MCO to MMCT has previously been documented in high-resolution palynological analyses. For example, Jiménez-Moreno et al. (2005) investigated a core from the Pannonian Basin and observed a decline in megathermic taxa at the transition from MCO to MMCT. Also, Ivanov and Worobiec (2017) reported a decrease in thermophile taxa for the transition for Bulgaria and Poland. In southwestern Anatolia, Kayseri-Özer et al. (2014) investigated three localities in the Muğla-Ören area south of the Yatağan Basin, which are dated by vertebrate fossils as early and late MN5 and thus correspond to the MCO. These authors report a few warmth-loving elements (palms, Avicennia) that are missing in the younger strata of the Yatağan Basin. This could be due to the deltaic setting of these types of flora as opposed to the intramontane setting of the Yatağan Basin flora. In general, the flora of the Muğla-Ören area is very similar to the flora of the Yatağan Basin (Bouchal et al., 2017). However, a striking difference to the MN6 and MN7 + 8 assemblages of the Yatağan Basin is the almost entire absence of herbaceous taxa (non-arboreal pollen) in the MN5 assemblages of Ören (see Figs. 7–9 in Kayseri-Özer et al., 2014). This may indicate the presence of more closed forest vegetation of the laurisilva type. The extant laurisilva or laurel forest is a type of subtropical forest found in areas with high humidity and relatively stable, mild temperatures. The assemblages of the Yatağan Basin, show fluctuating arboreal to non-arboreal pollen (AP:NAP) ratios with a peak in NAP in the transition zone MN6 to MN7 + 8 (pollen zone PZ 2–3). This peak could possibly correspond to a sharp cooling detected in the benthic stable isotopic data at 13.9–13.8 Ma (Holbourn et al., 2014). In the European mammal stratigraphy (Neubauer et al., 2015), the boundary MN6 to MN7 + 8 is at 13.9 Ma. Above PZ 2–3, the radiometrically dated Yeni Eskihisar pollen assemblage clearly belongs to MN7 + 8. Here, and in the nearby locality Çatakbağyaka woody taxa (including some warmth-loving taxa) are again more prominent. Thus, although the correlation of pollen zones 2–3 with the cooling event at 13.9–13.8 Ma is highly speculative, it is clear that the MCO in southwestern Anatolia was characterized by laurisilva vegetation with little contribution of herbaceous taxa. During the MMCT the main woody taxa did not change much, but herbaceous taxa played a much greater role. This indicates higher structural complexity of the vegetation. The presence of early hominids in western Anatolia during this time might be connected to this more complex vegetation. It is unclear at present, whether these changes were accompanied by changes in concentrations of atmospheric CO₂. The compilation of reconstructed CO₂ values across the Cenozoic from hundreds of proxy data (Beerling and Royer, 2011) shows that there is no agreement between different proxies for the MCO and the subsequent middle Miocene climate cooling. Phytoplankton stable isotopic data suggest nearly stable CO₂ concentrations (MCO, 227–327 ppm; MMCT, 265–300 ppm; see Table S1 of Beerling and Royer, 2011). In contrast, stomata densities from fossil leaves suggest a pronounced decline in CO₂ across this interval.

5 Conclusion

Here we used three proxies to infer climate, palaeoenvironments, and biogeographic affinities of middle Miocene floras of southwestern Anatolia. We showed that the palaeobotanical record resolves transitions from the warm MCO (16.8–14.7 Ma) to the MMCT (14.7–13.9 Ma) and a more pronounced cooling at 13.9–13.8 Ma, mainly expressed in the changing and fluctuating ratios between AP and NAP taxa. Using threshold percentages for main tree taxa, we further show that although NAP values significantly increased during the MMCT, AP taxa remained relatively abundant, signifying the coexistence of forested and open landscapes during this transition. In addition, the biogeographic analysis indicates mainly northern hemispheric biogeographic affinities of the middle Miocene flora of southwestern Anatolia and thus invalidates previous comparisons with tropical environments. Tropical climate conditions are also rejected by the Köppen signatures of the investigated flora and by the CLAMP analysis. Finally, the CLAMP data readily distinguish between strongly seasonal Cs and Cw and fully humid Cf climate types. More combined macrofossil and microfossil studies are needed for the Neogene of Turkey in order to establish a robust framework of terrestrial climate evolution in this important region.