New paleohydroclimate record of the MIS 5e/5d transition from Yelini Cave, central Anatolian region of Türkiye

Abstract This study presents the results of a detailed paleoclimate investigation on stalagmite YL-1 (Yelini Cave, Türkiye). YL-1 grew between 117.13 (+0.57/−0.44) ka and 114.87 (+1.63/−2.89) ka within Greenland Stadial 26, indicating a positive moisture balance during the stadial conditions in this semi-arid region. Rainfall is significantly affected by sub-cloud and surface evaporation and decreasing net effective winter precipitation is recorded by high isotope values. Enriched δ18O and δ13C at 116.65 (+0.51/−0.39) ka are interpreted as a drought event that took place ca. 400 years before the end of the MIS 5e. This event, which was reported simultaneously in marine and terrestrial archives in the Northern Hemisphere, is a result of decreased cyclone activity linked to weakening of the Atlantic meridional overturning circulation. 87Sr/86Sr values of YL-1 are close to the host-rock values. Decreased 87Sr/86Sr ratio at 116.67 (+0.54/−0.38) ka reflects the intensified water–rock interaction due to lower precipitation. Along with prior calcite precipitation effect, this is also observed by increased Mg/Ca and Sr/Ca, while low P, Cu, Be, Y, and Zr concentrations indicate a lowered amount of soil-derived colloidal material. The MIS 5e/5d transition is marked by reduced insolation and enriched δ18O at 116.24 (+0.53/−0.86) ka. The Greenland Interstadial 25 phase at 115.87 (+0.83/−1.71) ka is represented by more negative δ18O and δ13C.


INTRODUCTION
Glacial-interglacial transitions provide precise archives for climate events.This is due to the influence of orbital parameters on the timing of events occurring in these intervals, as well as due to the expansion of the ice sheets and its effect on the oceans and moisture balance.Marine Isotope Stage (MIS) 5 interval (ca.130-74 ka) is represented by the last interglacial period (LIG) and the onset of the last glacial (Capron et al., 2010).The MIS 5e substage covers the interval of 130-116 ka (130 through ca.120 ka in the eastern Mediterranean [EM] and eastern Europe based on terrestrial archives [Demény et al., 2017]), during which the sea level is estimated to have been 6-9 m higher than the current global average (Kukla et al., 2002;Shackleton et al., 2003;Dutton et al., 2015).The LIG and following glacial inception periods are climatically important periods since they have the potential to allow us to reconstruct the possible responses to climatic variations triggered by current global warming (Govin et al., 2015).
The continuous ca.123 ka-long δ 18 O record of the North Greenland Ice Core Project (NGRIP) ice core is very valuable because it also includes a part of the LIG period (NGRIP members, 2004).High temperatures during the LIG are represented by high values of δ 18 O (−32‰) (NGRIP members, 2004).During the LIG, the Mediterranean woody taxa reached its highest abundance (Brauer et al., 2007).Later, this warm interval was interrupted by glacial expansion and some cold-water events in the North Atlantic (Kukla et al., 2002).The NGRIP records indicate a gradual depletion of δ 18 O towards 116 ka with a transition to cold conditions.The end of the cooling phase is defined as the Greenland Stadial 26 (GS 26) interval, which marks the MIS 5e/ 5d transition.
A cooling of 4°C is suggested to have occurred between the C25 cold-water event and the last interglacial optimum (Oppo et al., 2006).In addition, the planktonic δ 18 O record indicates the C26 cold-water event during early ice sheet growth, pointing to a 2-3°C decrease in North Atlantic sea surface temperature (Chapman and Shackleton, 1999) and the warm period of MIS 5e lasting ca.9-10 ka that gradually ended with the C26 and C25 cold-water events (Oppo et al., 2006).The C25 cooling event associated with the end of GS 26 also indicates a decrease in woody taxa, pointing towards a decrease in moisture availability in the eastern Anatolian pollen records (Pickarski et al., 2015).Proxies of the LIG and the beginning of the glacial period have been discussed in detail using glacial (Johnsen et al., 1997;NGRIP members, 2004), marine (McManus et al., 1994;Fronval and Jansen, 1997;Chapman and Shackleton, 1999;Sánchez Goñi et al., 1999, 2012;Oppo et al., 2001Oppo et al., , 2006;;Knudsen et al., 2002;Melki et al., 2010;Galaasen et al., 2014;Irvali et al., 2016), and speleothem records of the North Atlantic and European regions (Drysdale et al., 2005(Drysdale et al., , 2007(Drysdale et al., , 2009;;Boch et al., 2011;Moseley et al., 2015;Vansteenberge et al., 2016Vansteenberge et al., , 2019;;Dumitru et al., 2018;Nicholson et al., 2021).The effects of climatic changes in the eastern Mediterranean occurring in both the last interglacial and glacial inception periods have been reported in several studies (e.g., Bar-Matthews et al., 2003;Nehme et al., 2015;Columbu et al., 2019), but are not fully known in Türkiye.For this reason, it is important to investigate the potential effects of sub-millennial climate changes during the last interglacial-glacial inception on the Anatolian Peninsula.
Due to its geographic setting, Türkiye is affected by various climatic domains (e.g., Atlantic, Mediterranean, and Asian; Rowe et al., 2012).Regional climate is dominated mainly by humid and cold air masses originating from the north and northwest, as well as by mid-latitude subtropical anticyclones (Akçar and Schlüchter, 2005;Türkeş, 2020).The high-pressure systems initiating from Siberia, which are cold and dry in winter, and the extensions of the low-pressure monsoon system originating from Asia in summer are other climate systems affecting the region (Akçar and Schlüchter, 2005;Türkeş, 2020;Jacobson et al., 2022).Precipitation in the region is mainly controlled by the North Atlantic Oscillation (NAO) (Türkeş, 2020).While strong precipitation conditions are observed in winter and autumn months in Türkiye during the NAO negative phase, weaker precipitation conditions prevail during the positive phases (Türkeş andErlat, 2003, 2005).Therefore, considering its effect on winter precipitation, the Anatolian Peninsula is extremely important in order to explore the eastern extensions of climatic changes in the last glacial inception period.
Stalagmites (secondary cave carbonates) can physically and chemically record climatic and environmental events during the time of their formation in isotopic equilibrium with the dripping water (McDermott et al., 2006).Precise dating of these deposits up to 600,000 years using the U-series dating technique (White, 2004) has played an important role in establishing long-term terrestrial climate records.The chronological records of the LIG are generally based on the alignment (correlation) technique of various proxies, which may lead to some discrepancies in ages.High-precision age data obtained from speleothems, however, can provide good age-calibrated climate records (Govin et al., 2015).
Speleothem-based paleoclimate studies in Türkiye are mostly focused on the last ice age or the Holocene period, some with important linkages with early human settlements (Fleitmann et al., 2009;Göktürk et al., 2011;Rowe et al., 2012;Ünal-İmer et al., 2015, 2016a;Koç et al., 2020;Erkan et al., 2021;Jacobson et al., 2021Jacobson et al., , 2022)).The oldest records from the region are from Sofular (Badertscher et al., 2011) in NW Türkiye and Dim Cave (Rowe et al., 2020) from SSW Türkiye.The study on Dim Cave (Alanya) by Rowe et al. (2020) focused on the beginning of the last interglacial period, but the climate record does not extend into the last interglacial-glacial transition.In addition, the Sofular Cave record characterizes different climatic conditions that do not reflect true Mediterranean climate due to significant humidity over the Black Sea.There is a lack of speleothem-based record of the last glacial inception in Türkiye.This study fills this gap by presenting new stable and radiogenic isotope (δ 13 C and δ 18 O, 87 Sr/ 86 Sr) and trace element (Ca, Mg, Sr, P, Cu, Y, Be, and Zr) data from a MIS 5 stalagmite (YL-1) collected from Yelini Cave (Fig. 1) and provides the first multi-proxy climate record of the MIS 5e/5d transition from the central Anatolia.

Yelini Cave
The Yelini Cave (also known as the Yılanlı Cave) is located near the town of Kayakent in the Günyüzü district of Eskişehir Province in central Türkiye.The cave is accessed through a narrow entrance along the southwest side of the Sivrihisar Mountains (Fig. 1).Yelini has archaeological importance because artifacts there include remains and cave paintings suggesting prolonged periods of human use in Neolithic times (Nazik et al., 2001).
Yelini Cave (39°14'23.28"N,31°48'25.33"E)(Fig. 1), located at an altitude of 1105 m above sea level (asl), is a horizontally extending fossil cave with a total length of 271 m (Nazik et al., 2001) (Fig. 1b).Yelini Cave is dry except for the drip waters seeping from the ceilings during wet seasons (Nazik et al., 2001) through a thin cover of limestone in the region.Its deepest point lies ∼26 meters below the cave entrance.Relative humidity within the cave varies from 45% in the entrance to 79% in the deepest gallery where the YL-1 stalagmite was collected in 2019 (Fig. 1).

Climate
The modern climate of the Anatolian Peninsula is heterogeneous.Located in the middle latitude temperate zone, Türkiye experiences four seasons throughout the year due to its close proximity to the sub-tropical zone (Türkeş, 2020).While the coastal regions are represented by a Mediterranean climate, the inland Anatolian Plateau experiences a continental climate (Bozkurt and Şen, 2009).Mountain belts extending parallel to the coast in the southern and northern regions prevent the marine climate from reaching the inner regions (Ünal et al., 2012).Therefore, precipitation received by the coastal regions is significantly higher than in the inner high plateaus of central Anatolia and the transportation of moisture to the interior is provided through the E-W trending valleys between the mountain ranges extending perpendicular to the Aegean coast in the western part of Anatolia (Bozkurt and Şen, 2009).Because of these topographical barriers, the central Anatolian region, where the Yelini Cave is situated, has semi-arid climate conditions (Altın et al., 2012).For this reason, winters are cold and rainy, spring months are cool and rainy, and summers are mildly rainy (Türkeş, 2003).According to data from the Turkish General Directorate of Meteorology for the year 2020, the average annual temperature of the Eskişehir region is 11°C, and the average annual precipitation is 365 mm.The highest precipitation occurs in January, March, May, and June, while lowest precipitation is recorded in July and November (Fig. 2).Potential evapotranspiration exceeds precipitation in average conditions throughout the year (Türkeş, 2020).
The predominance of winter precipitation in the region reflects the annual NAO effect.The dominant influence of NAO on Europe and the Mediterranean basin is particularly important in terms of received precipitation in western and central parts of Anatolia (Türkeş and Erlat, 2003).In the negative phase of NAO, weakening of Azores high pressure and shifting of maritime polar air masses to a southward route and evaporation over the western Mediterranean increase the extent of rainy conditions (Bozkurt and Şen, 2009).

Sample description
The YL-1 stalagmite was collected from the interior of the cave far from the entrance in an area with limited ventilation (Fig. 1b).The total length of the stalagmite is 14.1 cm, and due to the absence of laminae and textural evidence for possible secondary crystallization and deformation, the upper 0.9 cm section was not studied (Fig. 3).
The laminae are quite prominent along the stalagmite growth.Thinly developed laminae (1 mm thickness or less) of the sample are generally light brown, cream, beige, and white in color.A dark-brown layer is observed at 7.7-8.0cm from the apex, which divides the stalagmite into two distinct growth periods (Fig. 3).Changing of growth axes is evident for these growth periods.The lower and upper axes are 5.0 and 7.5 cm long, respectively (Fig. 3).The width of the YL-1 stalagmite is ∼6 cm, but there are minor variations in thickness.

Mineralogy
Thin-section petrography and XRD studies were carried out for mineralogical investigations.Three thin sections were prepared from one-half of the YL-1 stalagmite at the Ankara University  Earth Sciences Application and Research Center (YEBİM) and textural changes of the laminae were examined under a polarizing microscope.
XRD analyzes were performed on laminae from which age data were obtained.The analyses were carried out with the Panalytic-Aeris device in YEBİM by selecting ∼500 mg of powder sample for each.

U-Series dating
The U-series dating method was used for dating of the YL-1 stalagmite.Ten sub-samples, each weighing 100-150 mg, were collected from the stalagmite (Fig. 3).Microsampling was performed at the points of 1, 13, 25, 33, 45, 72, 82, 98, 112, and 131 mm from the top of the stalagmite.Sampling was carried out as much as possible from the impurity-free parts of the calcite layers.Drill bits used for microdrilling were cleaned with distilled water and ethyl alcohol after each sampling.
U-Series dating of the samples was performed at the University of Queensland (Australia) Radiogenic Isotope Laboratory (School of the Environment) using the Nu Plasma HR Multi-collector Inductively Coupled Plasma Mass Spectrometer (MC-ICP-MS).The analysis was carried out using the analytical protocols described by Zhao et al. (2001) and Ünal-İmer et al. (2015), and the results are reported in Table 1.

Stable isotope analysis
For carbon and oxygen isotopic analyses, 127 calcite sub-samples were collected from the stalagmite at 1 mm intervals using an micro hand drill.In addition, one marble wall-rock sample was also analyzed for its carbon and oxygen isotope composition.
Carbon and oxygen isotope analyses were performed at the Environmental Isotope Laboratory of the University of Arizona (USA).δ 18 O and δ 13 C isotope compositions of carbonates were measured with a gas ratio mass spectrometer (KIEL-III) coupled with a carbonate preparation device (Finnigan MAT 252).Powder samples were treated with phosphoric dehydrate under vacuum at 70°C.Isotope ratios were analyzed with an accuracy of ± 0.1‰ for δ 18 O and ± 0.08‰ (1σ) for δ 13 C and calibrated based on repeat measurements of NBS-19 and NBS-18.All oxygen and carbon isotope ratios are shown in standard notation and reported relative to Vienna Pee Dee Belemnite (VPDB) (Supplementary Table 1) with calculated StalAge age results and host rock δ 18 O and δ 13 C isotope ratios.
In addition, for the Hendy Test (Table 2), carbon and oxygen isotope analyses were carried out on 10 sub-samples collected from 2 prominent layers of 27 and 103 mm from the top of the stalagmite.Microsampling was performed at ∼1 mm intervals perpendicular to the growth axis of the stalagmite.Isotope analyses were conducted using the Termofinnigan Delta PlusXP mass spectrometer at the Central Laboratory of the Middle East Technical University (METU), Türkiye.
87 Sr/ 86 Sr isotope and trace element concentrations Nine sub-samples (eight from the stalagmite and one from marble host rock) were selected for trace element and strontium isotope ( 87 Sr/ 86 Sr) analyses.Sampling was performed at the sampling locations selected for dating.Analyses were carried out at the University of Queensland using the Thermo X-series ICP-MS device for trace element analysis and Thermal Ionization Mass Spectrometer (TIMS) for the strontium isotope analysis.Detailed methodology is given in Ünal-İmer et al. (2016b).

Calcite textures
Different calcite textures have been identified in the YL-1 stalagmite (Fig. 4a-c).Calcite crystals exhibit an irregular structure in the interval of ca.116.16 to ca. 116.32 ka.Polycrystals of calcite are dispersedly intertwined (Fig. 4c) along with high porosity.A more-arranged calcite fabric appears at ca. 115.49 to ca. 115.83 ka, showing a distinctive layered structure.Layers are characterized by alternations of light-and dark-colored bands.While the light-colored layers vary between ∼0.05-0.1 mm in thickness, dark-colored layers have thicknesses changing between 0.2 and 0.9 mm (Fig. 4b).A transition to a coarse-grained compact texture is observed in the layers at around 115.27 ka.Anhedral or subhedral mosaic calcite crystals have variable diameter size, ranging from ∼10 μm to 1 mm.There is no significant porosity developed between the crystals (Fig. 4a).
Table 1.MC-ICP-MS U-Series age results (ka) of the YL-1 stalagmite (Yelini Cave).Ages in bold are out of stratigraphic order and were not utilized in StalAge calculations (Scholz and Hoffmann, 2011) XRD results indicate that the main mineral in the stalagmite is calcite.In addition to calcite, the results also show minor aragonite peaks.However, no residual aragonite was detected petrographically in the thin sections.

U-series dating and age-depth model
Results of U-series dating analyses of 10 samples from the YL-1 stalagmite are presented in Table 1.Accordingly, deposition of the stalagmite started at 117.15 ± 0.60 ka and stopped at 109.05 ± 0.74 ka.The analyses at depths of 13, 25, 33, and 72 mm yielded dates beyond the growth period of the stalagmite.The outlier ages may be related to the transport of clastic thorium ( 232 Th) in groundwater and its inclusion in the stalagmite (White 2004;Fairchild and Baker, 2012).Therefore, the age-depth model was constructed through the StalAge code in R (Scholz and Hoffmann, 2011;Scholz et al., 2012) using the ages (n = 6) in stratigraphic order within the margin of error (Fig. 3).
There is a weak positive correlation between the δ 13 C and δ 18 O values of the YL-1 stalagmite (R 2 = 0.31) (Fig. 5b).The sub-  sample from 1 mm (distance from top) has significantly high oxygen and carbon isotope values (−8.41‰ and 0.74‰, respectively).
If the dripping is extremely slow or interrupted during the end of calcite deposition, isotopic enrichment may occur as a result of prolonged separation of CO 2 from the solution (Frisia et al., 2000), resulting from low effective infiltration inside the cave.

Hendy Tests
The kinetic fractionation in stalagmites may mislead the paleoclimate data obtained from the calcite (Hendy, 1971).This causes the light isotopes to be preferentially separated from the solution and concomitant enrichment of heavy isotopes in the stalagmite.
For this reason, we applied the Hendy Test (Hendy, 1971) to the YL-1 stalagmite.According to the Isotopic Equilibrium Test study, the change in δ 18 O values along a single lamina is expected not to exceed 0.5‰ (Hendy 1971;Ünal-İmer et al., 2020).δ 18 O values of YL-1 stalagmite do not show such a variation (Table 2; Fig. 6).In addition, farther from the axis there are no reciprocal changes between δ 13 C and δ 18 O that indicate kinetic fractionation (Fig. 6).Therefore, it is suggested that CO 2 separation from dripping water was slow during stalagmite formation (Ünal-İmer et al., 2020) and isotopic equilibrium conditions were maintained.How effectively the Hendy Test applications in stalagmites reflect the isotopic equilibrium is still controversial since the thickness of stalagmite layers at the growth center and edges may vary.Different kinetic fractionation effects may occur in the same layer along with possible problems with the sampling (Dorale and Liu, 2009).In addition, it is not possible to establish a precise oxygen isotope equilibrium between calcite and dripping water (McDermott et al., 2006).The Hendy Test is widely used to investigate such kinetic effects, although it does not give as precise results as the comparative isotope records from different stalagmites (Dorale and Liu, 2009).
87 Sr/ 86 Sr and trace element compositions 87 Sr/ 86 Sr ratios of YL-1 stalagmite vary in a narrow range from 0.707946 to 0.707987 (Table 3).These values are very close to the Sr isotope ratio of the host rock (0.707851), indicating a close genetic link between the marble host rock and the YL-1 stalagmite.The highest and lowest 87 Sr/ 86 Sr ratios in stalagmite are recorded at 71 mm and 82 mm (at the beginning of a dark layer), respectively.

Speleothem growth and deposition rate
Speleothem formation is an important indicator of climatic conditions since the growth of speleothem depends crucially on water availability (Genty et al., 2006).Although growth rate is an  important proxy, growth mainly responds to changes in temperature, precipitation, and soil activities (Kaufmann and Dreybrodt, 2004).Sufficient Ca 2+ dissolution from the marble wall-rock and the amount of CO 2 in the soil, which controls the solution acidity, are closely correlated with the extent of calcite precipitation in speleothems (Kaufmann, 2003;Vansteenberge et al., 2019).In addition, water loss through evaporation in the unsaturated region, evapotranspiration, and surface flow also can reduce the effect of dripping water (Vaks et al., 2003).Therefore, lowering of the P/E ratio may cause a decrease in growth rate of speleothems or even hiatuses in speleothem growth.Speleothem growth occurs under humid climatic conditions (Genty et al., 2006;Fairchild and Baker, 2012).Based on the age-depth model, the YL-1 stalagmite was deposited at 117.13 (+0.57/−0.44) to 114.87 (+1.63/−2.89)ka and is characterized dominantly by thin layers of fine-grained calcite (Fig. 4).
The growth rate of YL-1 is found to be ∼0.10 mm/yr from the beginning towards ca.116.94 ka (Fig. 8).After this period, the deposition rate increased to 0.12 mm/yr and slightly dropped to 0.10 mm/year at ca. 116.83 ka.The first notable decrease in growth rate to 0.63 mm/yr is seen at ca. 116.67 ka, which corresponds to the beginning of the detrital-rich brown layer at ca. 116.67 ka (77-80 mm depth; Fig. 3).Age results, however, suggest that the hiatus here was not prolonged.A second noticeable growth rate decrease, down to 0.04 mm/yr, appears at ca. 116.08 ka where the calcite fabric displays interrupted crystallization and dendritic-type porous fabric (Figs. 4 and 8).Generally slow growth rate throughout the YL-1 deposition and even slower growth rates in the upper 1 ka section are linked to the lowering P/E (precipitation/evaporation) ratio in the cave.The onset of stalagmite deposition coincides with the GS 26 cold phase represented in the Greenland ice core records (NGRIP members, 2004).This indicates the presence of a positive moisture balance in the region during the cold phase, which is not surprising considering the semi-arid climatic conditions of the central Anatolian region.Under current climatic conditions, a relatively low annual precipitation rate of 300-500 mm/yr and average summer temperatures reaching 23°C cause sub-cloud and/or surface evaporation in the region (Schemmel et al., 2013).Therefore, "effective moisture" (precipitation−evaporation) controlling speleothem growth is of great importance (Hodge et al., 2008).Increased overall precipitation and lower temperatures in the region indicate higher "effective precipitation" or positive P-E.
δ 18 O and δ 13 C variations In order to define the mechanisms controlling changes in δ 18 O cc , it is important to outline controls on the δ 18 O water in that region.There are a number of "isotopic effects" that cause variations in δ 18 O water , including the amount of precipitation, latitude differences, ice volume changes, seasonal changes, temperature, and changes in moisture source (Lachniet, 2009).In addition, other factors such as evaporation at the surface and in the epikarst region, water residence time in the karstic zone, groundwater mixing, and calcite precipitation rate provide important controls on δ 18 O cc (McDermott et al., 2006;Rowe et al., 2012).Hence, the factors acting on δ 18 O cc represent a complex system (Lachniet, 2009;Rowe et al., 2012).
At Yelini, the rock column (marble/recrystallized limestone) above the cave has low thickness (5 m at most) (Nazik et al., 2001).This might have reduced the possibility of mixing of dripping water with groundwater having different isotopic signatures due to lower residence times (Lachniet, 2009).Therefore, it is thought that atmospheric processes exerted a major control on δ 18 O cc .
Previous studies on speleothem δ 18 O cc in the Mediterranean region examined the "amount effect," which is based on the inverse correlation between rainfall amount and δ 18 O ppt , the "source effect" focusing on Atlantic-vs.Mediterranean-sourced humidity that varies during stadial and interstadial transitions, and the "temperature effect" explaining the altitude-based inverse correlation between cold conditions and δ 18 O water (Bar-Matthews et al., 1999;Jex et al., 2010;Boch et al., 2011;Rowe et al., 2012;Dumitru et al., 2018;Columbu et al., 2019).In addition to these mechanisms, there are also some studies in Türkiye that have reported a negative correlation between effective winter rainfall and δ 18 O water .The speleothem record of Akçakale Cave in NE Türkiye showed a strong and consistent correlation between the October-January rainfall records and δ 18 O cc values (Jex et al., 2010), which is attributed to recharge of the karst aquifer by the autumn-winter precipitation with depleted δ 18 O water values.The isotope data of dripping waters, local springs, and streams around the Karaca Cave, which is located in southern Türkiye, also display a strong linkage with contribution from winter rainfall (Rowe et al., 2012).However, in both of these examples, the overburden above the cave is thick with high rainfall amounts during the winter, allowing for longer residence time for infiltrated water, which is less exposed to evaporation and water mixture inside the epikarst.Therefore, the P−E balance is more   (Dilaver et al., 2018).Similarly, stable water isotope data from CAP (Ankara) and Global Network of Isotopes in Precipitation (GNIP) Ankara station data suggest a strongly seasonal (winter dominant) rainfall pattern along the CAP (Schemmel et al., 2013).This is explained by a positive precipitation-evaporation balance during the cold season and winter-dominated groundwater recharge.Consequently, increased or decreased intensity of winter precipitation is the most likely mechanism reflecting the changes in δ 18 O cc values of the YL-1 stalagmite (Fig. 8).Lower δ 18 O cc values of the YL-1 stalagmite are attributed to the increased effect of net  The anomalies in all proxies between 116.67 (+0.54/−0.38)ka and 116.65 (+0.51/ −0.39) ka are noted (green vertical bar).δ 18 O at 116.24 (+0.58/−0.86)ka (MIS 5e/5d transition) increases by 1.79‰ (yellow vertical bar).δ 13 C values record this transition at 116.07 (+0.67/−1.24)ka as a weaker anomaly (yellow vertical bar).The decrease in growth rate is also noticeable during this period.
effective winter precipitation and higher δ 18 O cc values are due to decreasing net effective precipitation arising from the increase in atmospheric evaporation.Adding to this is the plausible evaporation within the epikarst of Yelini Cave, due to the location of the cave close to the surface.
Changes in speleothem δ 13 C cc are directly associated with variations in isotopic contributions of carbon reservoirs including soil and atmospheric CO 2 as well as the carbonate host rock (McDermott et al., 2006;Dumitru et al., 2018).δ 13 C values are affected by fluctuations in isotopically light soil CO 2 (Hodge et al., 2008), therefore, speleothem δ 13 C cc values strongly depend on carbon isotope composition of vegetation (C 3 and C 4 plants) growing above the cave (McDermott et al., 2006).Changes in δ 13 C cc values also may be affected by several factors that include stalactite-stalagmite fractionation (evaporation and CO 2 degassing), CO 2 gas loss of groundwater in the aquifer zone, and short water-residence times in the soil (i.e., inability to establish an equilibrium between soil water and soil CO 2 ) (Baker et al., 1997).
Biogenic CO 2 increase due to increase in rainfall during hot and humid periods has a negative-trending effect on δ 13 C cc (Genty et al., 2003).C 4 -type plants reflecting arid environments and C 3 -type plants reflecting humid environments cause different δ 13 C values in speleothems depending on the climate (δ 13 C for C 4 and C 3 plants are −6 to +2‰ and −14 to −6‰, respectively) (McDermott et al., 2006).On the other hand, decreasing CO 2 contribution of biogenic carbon during cold and dry periods, increasing dissolved carbon content, PCP, primary carbonate host rock, and atmospheric dissolved carbon concentrations result in high δ 13 C cc values in speleothems (Baker et al., 1997;Genty et al., 2006;Dumitru et al., 2018).In general, the δ 13 C values of a speleothem are products of this equilibrium.In the YL-1 stalagmite, δ 13 C values varying from −14 to−6‰ indicate C 3 -type vegetation (McDermott et al., 2006) similar to the other cave records from Anatolia (e.g., Ünal-İmer et al., 2015).Modern vegetation is quite scarce at the studied cave site, but the δ 13 C profile of the YL-1 stalagmite is associated with changes in effective winter precipitation (precipitation−evaporation).The changes towards more negative values can be explained by rises in net effective precipitation and growing seasonal moisture in the soil and carbonate rock/atmospheric components during dry seasons (Fig. 8).
87 Sr/ 86 Sr ratios Sr isotope values ( 87 Sr/ 86 Sr) measured in speleothems reflect the relative contributions of Sr dissolved from the marble host rock and exogenic Sr accumulating in the soil (Goede et al., 1998;Bar-Matthews et al., 1999;Frumkin and Stein, 2004;Li et al., 2005).The negligible extent of isotopic fractionation of 87 Sr/ 86 Sr during calcite formation enables assessment of environmental factors affecting the 87 Sr/ 86 Sr values of dripping water in a cave (Goede et al., 1998).In addition to the Sr isotope composition of the host rock, different radiogenic Sr sources (e.g., dust transport, marine salts, volcanic activity) are reflected in the 87 Sr/ 86 Sr values of the speleothems (Banner et al., 1996;Goede et al., 1998;Bar-Matthews et al., 1999;Verheyden et al., 2000;Frumkin and Stein, 2004;Li et al., 2005).This provides insights into the origin of radiogenic Sr, the effect of varying climatic conditions, and changes in the karst.
Previous studies have shown that weathering of dolomitic host rock with an increase in precipitation might cause low 87 Sr/ 86 Sr values in speleothems (Ayalon et al., 1999).Verheyden et al. (2000) reported the dominant effect of low 87 Sr/ 86 Sr bedrock endmembers and the effect of high 87 Sr/ 86 Sr clastic-phase endmembers with their time-dependent contributions.In addition to the dominant host rock contribution represented by low 87 Sr/ 86 Sr values during humid periods, terrestrial dust can cause high 87 Sr/ 86 Sr ratios in speleothems (Frumkin and Stein, 2004).The dissolution of silica-rich metamorphic rocks (as cave basement) or quartzofeldspathic clastics in bedrock also may cause high 87 Sr/ 86 Sr in speleothems (Ünal-İmer et al., 2016a;Carolin et al., 2019).87 Sr/ 86 Sr ratios of the YL-1 stalagmite mostly indicate genetic signatures that are similar to the host carbonate rock (0.707851) (Fig. 9).Therefore, the contribution of terrestrial dust (0.7082-0.7086;Frumkin and Stein, 2004) input (from the Sahara Desert) is thought to be limited (Ünal-İmer et al., 2016a).87 Sr/ 86 Sr ratios increase from 0.707968 to 0.707985 in the period of 116.94 (+0.49/−0.35)ka to 116.83 (+0.58/−0.42)ka in the YL-1 stalagmite (Fig. 9).A decrease is observed at 116.67 (+0.54/−0.38)ka to the value of 0.707946, which is concurrent with enrichments in the δ 18 O and δ 13 C records (Fig. 8). 87Sr/ 86 Sr values initially show a similar increasing trend with δ 18 O; however, they then follow a different pattern at 116.67 (+0.54/−0.38)ka (Fig. 8).This difference between O-C systematics and 87 Sr/ 86 Sr ratios is likely due to the thin soil cover above the cave, increased groundwater residence time, and increase in the extent of water-host rock interaction possibly due to the decrease in groundwater recharge (Rowe et al., 2020).Lowered levels of 87 Sr/ 86 Sr at 116.67 (+0.54/−0.38)ka can be attributed to the negative P/E balance due to the decrease in winter precipitation (Fig. 9).This indicates higher contributions from the marble host rock as the water-rock interaction intensified.87 Sr/ 86 Sr ratio rises to 0.707987 at 116.51 (+0.38/−0.50)ka, which might reflect rapid return to wetter conditions after drought (Fig. 9).

Trace element concentrations
Variations in trace element compositions of speleothems result from changes in environmental conditions, which include chemical composition of the dripping water, the amount of precipitation, the types of elements in the soil, water-rock interaction, chemical interactions during the entrance of water to the cave, and trace element exchange reactions during calcite deposition (White, 2004).Water with a certain P CO2 leaking from the soil zone dissolves calcium carbonate until it attains equilibrium, and absorbs trace elements (White, 2004).Some elements react rapidly with the infiltrated water while others are leached by slowly flowing seepage waters (Fairchild and Treble, 2009).In the YL-1 stalagmite, both Sr/Ca and Mg/Ca ratios tend to increase from the beginning of deposition up to 116.67 (+0.54/-0.38)ka (Fig. 8), despite their low sampling resolution.Sr/Ca ratio attains its maxima at 116.67(+0.54/-0.38)ka.This might be related to the increased water residence times in the epikarst.The soil CO 2 increased by the vegetation controls the bedrock dissolution and has a significant effect on the chemical composition of the water entering the cave (Hellstrom and McCulloch, 2000).This may cause changes in the amount of calcite or dolomite dissolved, and trace element concentrations in the solution may also change accordingly (Fairchild et al., 2000).
Depending on the decrease in the amount of precipitation, residence time of water in the epikarst may increase, which is accompanied by elevated concentrations of elements such as Mg and Sr as a result of extended water-rock interaction (Fairchild et al., 2000;Hellstrom and McCulloch, 2000).However, linear relationships have been reported between Mg content and increase in the deposition temperature of speleothems (Şenoğlu, 2006;Zhou et al., 2011).Şenoğlu (2006) discussed this increase using a Holocene-aged stalagmite, also from the Yelini Cave, that is located in a semi-arid climate zone.Changes in Mg/Ca ratios in a specific period of the stalagmite are associated with dry (or humid) conditions.Because the stalagmite studied by Şenoğlu (2006) was collected close to the entrance of the cave, it probably was affected by changes in air temperature.However, the YL-1 stalagmite was sampled from the distal part of the cave from the entrance (Fig. 1b) so that it should be represented by negligible airflow (Fohlmeister et al., 2020).
Concurrent increases of δ 13 C and Mg/Ca-Sr/Ca ratios in YL-1 at ca. 116.65 ka are attributed to the Prior Calcite Precipitation (PCP) effect during a probable drought event.In addition to the dissolution process, proportional increases in Mg/Ca and Sr/Ca may also occur due to PCP along the flow path (Fairchild et al., 2000;Fairchild and Treble, 2009).As the P CO2 decreases in the aquifer zone due to drought, CO 2 is separated from the solution, which then becomes supersaturated with respect to calcite (Fairchild et al., 2000).In this case, Mg/Ca and Sr/Ca ratios may increase with the separation of calcium from the solution (Fairchild and Treble, 2009), leading to a simultaneous enrichment of Mg/Ca and Sr/Ca ratios and δ 13 C compositions of speleothems.This is because the change in speleothem δ 13 C values is controlled by changes in plant thickness and soil CO 2 .The simultaneous increases in YL-1 Mg/Ca, Sr/Ca, and δ 13 C at 116.67 (+0.54/−0.38)ka suggest that calcite was deposited along the flow path and PCP is the dominant mechanism, which could be considered as evidence of a decrease in soil CO 2 during that time.A subsequent negative shift in δ 13 C at 116.51 (+0.38/ −0.5) ka is also reflected in variations in Sr/Ca ratios as well as in relatively less-pronounced variations in Mg/Ca ratios (Fig. 8).This can be regarded as a transition to humid conditions during which the soil CO 2 increases, possibly due to elevated vegetation density (Fairchild and Treble, 2009).
We also examined concentrations of P, Cu, Y, Be, and Zr in the YL-1 stalagmite despite low-resolution sampling for these trace elements (Figs. 8 and 10).Studies on P indicate the existence of a linear relationship between rainfall and P concentration (Huang et al., 2001;Borsato et al., 2007).The increase in the concentrations of P, Y, and Cu in an Italian stalagmite (Grotta di Ernesto Cave, ER78 stalagmite) was explained by transport of these elements as colloidal particles from the soil zone during high water seepage in autumn (Borsato et al., 2007).P, Y, Cu, and Be concentrations in YL-1 show low values in the stadial period, whereas Zr contents are relatively high during the same interval (< 116.2 ka) (Fig. 8).Simultaneous decreases in P, Y, Cu, Be, and Zr at 116.67 (+0.54/−0.38)ka might be due to the decrease in the winter precipitation.During GI 25, P, Cu, Be, and Zr, as well as δ 13 C and δ 18 O, follow a similar (contemporaneous) change (Fig. 10).Relatively high P, Y, Cu, and Be concentrations at the onset of GI 25 (115.87ka) are associated with increased net effective precipitation and moisture availability in the soil zone (Figs. 8 and 10).

Comparison of the YL-1 stalagmite record with regional records
The climatic changes recognized in the YL-1 stalagmite record during the MIS 5e/5d transition do not match with other regional records (Fig. 11).The lack of a similar climate signal between the YL-1 stalagmite and the δ 18 O record of the So-17A stalagmite from Sofular Cave (Badertscher et al., 2011) is due to Türkiye's northern mountains (i.e., Pontides) acting as a barrier, preventing moisture from reaching the inland regions of Türkiye via the Black Sea.This can be associated with Sofular Cave in northwestern Türkiye mainly receiving moisture from the Black Sea (Fleitmann et al., 2009;Badertscher et al., 2011).The percentages of arboreal pollen from Tenaghi-Phillippon in northern Greece (Koutsodendris et al., 2023) and from Lake Van in eastern Türkiye (Pickarski et al., 2015) also do not show a similar trend with YL-1 stalagmite profiles (Fig. 11).However, an increase in arboreal pollen percentage at 115.6 ka has been labeled as the C25 cool event in the Lake Van pollen record.This event corresponds to GS 26 in the Greenland records (NGRIP members, 2004).With uncertainties in age, a correlation can be argued between this cold event in the Lake Van pollen record and the δ 18 O increase at 116.2 ka in the YL-1 stalagmite (Fig. 11).
There are no clear synchronous climatic signatures between the YL-1 stalagmite and records from the Levant Basin (Fig. 11).The YL-1 stalagmite does not cover a significant part of the LIG period (130-116 ka); indeed, its growth period (117-114 ka) is considered regionally to fall into the glacial inception period (Nehme et al., 2020).However, the δ 18 O increase at 116.2 ka in YL-1 stalagmite is synchronous with the end of MIS 5e, which is associated with the LIG.Different ages related to the last interglacial period are present in the EM records.The record from Pentadactylos Cave in Cyprus (Nehme et al., 2020) shows that the last interglacial conditions occurred between ca.127 and ca.123 ka, with colder and drier climate conditions prevailing during 123-110 ka (Fig. 11).A decrease in growth rate is also recorded at the end of the LIG.The end of moist conditions during the LIG period is noted at 120.3 ka in the Kanaan Cave speleothem record in Lebanon (Nehme et al., 2015), while the caves in Israel (Soreq and Peqiin) (Bar-Matthews et al., 2003;Grant et al., 2012) recorded the LIG at ca. 120 ka.Since the age interval of the YL-1 stalagmite does not reach back that far, it is not possible to compare these transitions.The isotopic changes in YL-1 during the glacial inception period (ca.117-114 ka) also do not correspond to EM records.This discrepancy could be due to topographical effects similar to those observed in the Sofular record in the north.The southern mountain ranges (Taurides) also significantly hinder moisture transport to central Anatolia.Therefore, it is expected that the YL-1 stalagmite would show differences from records primarily reflecting eastern Mediterranean moisture.Another possibility is that these EM records may have lower resolution than the YL-1 stalagmite has during the glacial inception (ca.114-117 ka) (Fig. 11).

Comparison of YL-1 stalagmite recordglobal relationships
The glacial inception period at Yelini Cave According to the age modeling, YL-1 stalagmite growth covers the time interval of 117.13 (+0.57/−0.44) to 114.87 (+1.63/ −2.89) ka, with slower growth rates after ca.116.2 ka (Figs. 8 and 12).The onset of stalagmite formation at 117.13 (+0.57/ −0.44) ka corresponds to the GS 26 period during which ice sheet growth was recorded in the Greenland ice cores (NGRIP members, 2004).As a result of reduction in insolation, a number of cooling events were reported in the North Atlantic by the end of MIS 5e during the increase of ice volume at 119-116 ka, which was recognized in the Greenland ice core records at the end of the LIG (Oppo et al., 2006;Galaasen et al., 2014;Irvali et al., 2016).An enrichment of δ 18 O values in the YL-1 stalagmite is observed at 116.83 (+0.58/−0.42)and 116.63 (+0.46/ −0.40) ka, which is close to the MIS 5e boundary.A similar enrichment trend is also noticeable for δ 13 C at 116.65 (+0.51/ −0.39) ka (Fig. 8).In addition, increases in Mg/Ca and Sr/Ca and decreases in 87 Sr/ 86 Sr ratios as well as decreases in P, Y, Cu, Be, and Zr concentrations are observed at 116.67 (+0.54/ −0.38) ka (Fig. 8).Further, a substantial decrease in growth rate and a significant reduction in the diameter of the stalagmite are also recorded in the same period (Fig. 3).Therefore, we argue that a multi-centennial drought event occurred in central Anatolia at 116.65 (+0.51/−0.39)ka, which is chiefly deduced from the YL-1 stalagmite δ 13 C profile as well as the anomalies in δ 18 O and Sr-trace element data (Fig. 8).
The YL-1 stalagmite record shows significant correlations with Atlantic and European records.The δ 18 O increase in the BAR-II stalagmite record from Baradla Cave between ca.118.5 and 117 ka is interpreted as an arid pulse (Demény et al., 2017) (Fig. 12m).This event is recognized with simultaneous increases in δ 18 O and δ 13 C values in the BAR-II stalagmite.This period can be compared with the high δ 18 O values at the beginning of the YL-1 growth.In the MD01-2443 marine record obtained from the Iberian margin (Fig. 12f ), a small increase in δ 18 O values of planktonic foraminifer (G.bulloides) at 116.8 ka corresponds to a decrease in SST (de Abreu et al., 2005; Voelker and de Abreu, 2011) (Fig. 12g).The HÖL-10 stalagmite from Hölloch Cave (Moseley et al., 2015) records a small δ 18 O decrease at 117.2 ka (Fig. 12l).A similar change is represented by the δ 18 O increase in the CAM-1 stalagmite from Campanet Cave at 117 ka (Dumitru et al., 2018) (Fig. 12i).
Elevated δ 18 O values at 116.24 (+0.58/−0.86)ka in YL-1 coincides with the depletion of δ 18 O at ca. 115.5 ka, representing the end of the last interglacial period in the Greenland ice record (NGRIP members, 2004) (Fig. 12h).It also matches with the upper limit of the benthic "plateau" of MIS 5e at 116.1 ± 0.9 ka (Shackleton et al., 2003).Neogloboquadrina pachyderma abundance reaching up to 90-95% in the eastern Nordic seas during the same time interval and the decrease in SST at the MIS 5e/ 5d boundary at ca. 116 ka indicate an interglacial-glacial transition in northwest Europe (Fronval and Jansen, 1997;Knudsen et al., 2002).The pollen record of Lake Monticchio in Italy also indicates a decrease in Mediterranean woody taxa at ca.115.8 ka (Brauer et al., 2007).
The MIS 5e/5d transition is also recognized in the δ 18 O records of the Campanet Cave CAM-1 stalagmite in Mallorca (Dumitru et al., 2018;Fig. 12i).The CAM-1 stalagmite records the MIS 5e/5d transition through increasing δ 18 O values towards 115.4 ka.The MIS 5e/5d transition in the YL-1 stalagmite record also overlaps with the transition recognized in the Schneckenloch Cave SCH-7 (Boch et al., 2011) and the Hölloch Cave HÖL-10 stalagmites (Moseley et al., 2015) in the northern Alps (Fig. 12).SST records from the Iberian margin MD95-2040 core do not show significant changes during the MIS 5e/5d transition and at 116.65 ka (de Abreu et al., 2003).In addition, the CC5 stalagmite record from Corchia Cave in Italy (Drysdale et al., 2005) does not provide a correlation with the YL-1 stalagmite record during this period (Fig. 12).
The YL-1 δ 13 C profile at 115.87 (+0.83/−1.71)ka also displays a negative shift during interstadial conditions.Well-laminated, fine-grained crystals in the stalagmite after the transition may reflect regular (or near steady) dripping (and rainfall) conditions (Fig. 4).Mosaic-textured calcite, which replaces fine-grained calcite layers at ca. 115.27 ka, are concurrent with an increase in δ 18 O profile (Figs. 4 and 5).Relatively coarser calcite crystals suggest a decrease in water availability and a decrease in the dripping regime of cave water, which in turn diminishes (or halts) stalagmite growth (Gonzalez et al., 1992;Ayalon et al., 1999).

CONCLUSIONS
This study presents a detailed investigation of the paleohydroclimatic conditions that prevailed in central Anatolia during the time period of 117.13 (+0.57/−0.44)ka to 114.87 (+1.63/−2.89)ka.δ 18 O ratio changes recorded in the Yelini Cave YL-1 stalagmite are affected by increases or decreases in effective winter precipitation regulated by annual changes in NAO.The δ 13 C profile of the YL-1 stalagmite reflects time-dependent changes in soil density and host rock (marble) interactions, which are controlled by variations in net effective precipitation (precipitation−evaporation) in this semi-arid region.In the YL-1 isotope record, the MIS 5e/5d boundary is marked at 116.24 (+0.53/−0.86)ka, which coincides (within the age error margin) with the MIS 5e/5d transition records in several global localities.It is also noted that the YL-1 record shows fewer correlations with regional continental records.We argue that a significant multi-centennial drought event occurred near the end of MIS 5e and ca.400 years before the MIS 5e/5d transition in central Anatolia.This event at 116.65 (+0.51/−0.39)ka (mainly based on swings of the δ 13 C profile) is likely due to the decrease in winter precipitation in the Anatolian Peninsula and correlates well with various global paleoclimatic archives.The YL-1 stalagmite also recorded the transition to the GI 25 wet phase at 115.87 (+0.83/−1.71)ka.As shown by trace element and stable isotope data, the transition to the GI 25 phase began with a significant increase in rainfall amount.Stalagmite growth ceased, with a gradual transition towards a colder period.The Yelini Cave YL-1 stalagmite record constitutes the first paleoclimate archive revealing the climatic conditions and timing of the MIS 5e/5d transition in central Anatolia, Türkiye.
Supplementary material.The supplementary material for this article can be found at https://doi.org/10.1017/qua.2023.56

Figure 3 .
Figure 3. Photograph of the YL-1 stalagmite half showing the growth axes and locations of samples dated by U-series (left).Age-depth model generated on StalAge (Scholz and Hoffmann, 2011; Scholz et al., 2012) using U-series ages of sub-samples at distances of 1, 45, 82, 98, 112, and 131 mm (right).Red lines indicate the age error boundaries.

Figure 5 .
Figure 5. (a) Diagram showing variations in δ 18 O and δ 13 C profiles of YL-1 stalagmite with respect to its depth.(b) Correlation between δ 18 O and δ 13 C values.

Figure 6 .
Figure 6.Hendy Test plots for the δ 18 O and δ 13 C analyses performed on laminae at 27 and 103 mm of YL-1.
biased towards winter season, allowing a negative correlation with δ 18 O cc .In the case of Yelini Cave, the evaporation component in the P−E balance might be important due to atmospheric circulation at the geographical location of the studied cave.The barrier effect of the mountain ranges in the northern and southern parts of Türkiye prevents moisture from moving inland(Erkan et al., 2021), resulting in a prevalence of semi-arid climatic conditions around Yelini Cave.In arid and semi-arid areas, near-surface evaporation processes are reflected as seasonal variations in δ 18 O (McDermott et al., 2006, and references therein).Schemmel et al. (2013) reported a modern meteoric water study suggesting that the general trend of increase in δD and δ 18 O from the sheltered sides of the mountain ranges in the southern and northern parts of Türkiye to the interior Central Anatolian Plateau (CAP) is the result of increased sub-cloud and surface evaporation due to drought conditions.Increased evaporation conditions significantly reduce effective precipitation (P-E) in the region.This also indicates significant seasonal changes on δ 18 O water values.The δ 18 rainfall data for the years between 2012 and 2016 taken from the Turkish State Meteorological Service (MGM) Ankara station, which is the closest meteorological station to Yelini Cave, indicate that winter precipitation is represented by more-negative δ 18 O values compared to other months

Figure 7 .
Figure 7. (a) Weak positive correlation between Sr/Ca and Mg/Ca.(b) Moderate positive correlation of Cu and Be concentrations with P concentration.(c) Weak positive correlation of Zr and Y concentrations with P concentration.

Figure 9 .
Figure 9. Graph showing the variation of 87 Sr/ 86 Sr with respect to stalagmite depth and U-series ages. 87Sr/ 86 Sr ratios of the stalagmite indicate contribution from the marble host rock. 87Sr/ 86 Sr values of the Jerusalem Cave stalagmite (Frumkin and Stein, 2004) are taken as the reference for terrestrial dust from the Sahara Desert.

Table 3 .
ICP-MS trace elements concentrations and TIMS 87 Sr/ 86 Sr ratios of the YL-1 stalagmite and marble host rock.