Stable isotopes in cave ice suggest summer temperatures in east-central Europe are linked to Atlantic Multidecadal Oscillation variability

The climate of east-central Europe (ECE) is the result of a combination of influences originating in the wider North Atlantic realm, the Mediterranean Sea, and the western Asian and Siberian regions. Previous studies have shown that the complex interplay between the large-scale atmospheric patterns across the region results in strongly dissimilar summer and winter conditions on timescales ranging from decades to millennia. To put these into a wider context, longterm climate reconstructions are required, but, largely due to historical reasons, these are lacking in ECE. We address these issues by presenting a high-resolution, radiocarbondated record of summer temperature variations during the last millennium in ECE, based on stable isotope analysis of a 4.84 m long ice core extracted from Focul Viu Ice Cave (Western Carpathians, Romania). Comparisons with both instrumental and proxy-based data indicate that the stable isotope composition of cave ice records the changes in summer air temperature and has a similar temporal evolution to that of the Atlantic Multidecadal Oscillation on decadal to multidecadal timescales, suggesting that changes in the North Atlantic are transferred, likely via atmospheric processes towards the wider Northern Hemisphere. On centennial timescales, the data show little summer temperature differences between the Medieval Warm Period (MWP) and the Little Ice Age (LIA) in eastern Europe. These findings are contrary to those that show a marked contrast between the two periods in terms of both winter and annual air temperatures, suggesting that cooling during the LIA was primarily the result of wintertime climatic changes.


Introduction
Rapid global warming (IPCC, 2018) and the ensuing suite of climatic changes that it triggers (Coumou and Rahmstorf, 2012) demands a clear understanding of the background mechanisms in order to be able to disentangle natural and anthropogenic processes (Haustein et al., 2017;IPCC, 2018). Especially important are high-resolution reconstructions of the past variability of different climatic variables, such as seasonal air temperatures, precipitation amounts and moisture sources, that allow for direct comparisons with the dynamics of natural forcing and further decipherment of the mechanisms of past and present climate dynamics. The last 1000 years are particularly significant, as the European climate has generally changed from warm to cold (the Medieval Warm Period-Little Ice Age transition; Jones et al., 2009) and back to warm (the present-day warming; Neukom et al., 2019). These transitions allow for the links between forcing and climatic response to be tested. While several global (Jones and Mann, 2004;Mann et al., 2009) and hemispheric (Moberg et al., 2005;Neukom et al., 2019;PAGES 2k Consortium, 2019;Ljungqvist et al., 2019) climatic reconstructions have been published, these made no seasonal differentiation -a task that became recently increasingly necessary to constrain seasonally distinctive climatic changes (e.g., Ljungqvist et al., 2019), as these respond to different forcing mechanisms (e.g., Pers , oiu et al., 2019). On multidecadal timescales, the summer climate over Europe is mainly influenced by the Atlantic Multidecadal Oscillation, or AMO (Schlesinger et al., 1994;Kaplan et al., 1998;Kerr, 2000;Knudsen et al., 2011Knudsen et al., , 2014. The AMO is a climate mode of variability associated with periodic anomalies in sea surface temperatures (SSTs) in northern, extra-tropical latitudes. The positive phase is characterized by positive SST anomalies spanning the whole North Atlantic Ocean and is associated with above normal temperatures over the central and eastern parts of Europe, whereas the negative phase is characterized by negative SST anomalies over the North Atlantic Ocean and is associated with below normal temperatures over the central and eastern parts of Europe. Over Europe, the influence of the AMO is clearest during summer (Sutton and Dong, 2012;Ionita et al., 2012Ionita et al., , 2017O'Reilly et al., 2017).
In temperate climatic regions, one of the most sensitive environmental archives are ice caves (Holmlund et al., 2005;Kern and Pers , oiu, 2013), i.e., rock caves hosting perennial accumulations of ice. In such caves, ice forms either by freezing of water or direct snow deposition in the entrance shafts (e.g., Mavlyudov, 2018). Several studies have shown that these deposits host a wealth of information on past climate variability. Stoffel et al. (2009) , Pers , oiu et al. (2017) and Sancho et al. (2018) have shown that proxies in cave ice forming during winter months record the changes in temperature and moisture sources, which are likely influenced by the dynamics of the North Atlantic Oscillation. Other studies have used pollen and plant macrofossils recovered from cave ice to reconstruct past vegetation dynamics (Feurdean et al., 2011;Leunda et al., 2019), while some works have used the accumulation rate of ice as indicators of past climatic variability (e.g., Kern et al., 2018) or atmospheric processes . Studies of ice caves in southern Europe have also highlighted the sensitivity of cave glaciers to summer climatic conditions (Colucci et al., 2016;Colucci and Guglielmin, 2019;Pers , oiu et al., 2020). Regardless of the deposition style, the ice records the original stable isotope composition of precipitation that further reflects changes in air temperature; thus, it is an important archive of past temperature and moisture source variability (Pers , oiu et al., 2011a, b). The Carpathian Mountains host ice caves (Brad et al., 2018) that preserve a large variety of geochemical information on past climate and environmental changes Kern et al., 2004;Citterio et al., 2005;Pers , oiu et al., 2017). Here, we present a reconstruction of summer climate variability and large-scale circulation drivers during the last 1000 years in east-central Europe based on the δ 18 O and δ 2 H values measured along an ice core drilled in Focul Viu Ice Cave (Western Carpathian Mountains, Romania).

Site information
Focul Viu Ice Cave (FV; 107 m long and ∼ 30 m deep) is located in the central Bihor Mountains, Romania (46 • 34 N,22 • 40 E; 1165 m above sea level; Fig. 1a ;Pers , oiu and Onac, 2019). The cave has a simple morphology (Fig. 1b, c) with a small entrance that opens into the Great Hall (68 m × 46 m), which, in turn, is followed by a narrow gallery (Little Hall, 20 m × 5 m). The ceiling of the Great Hall opens to the surface (Fig. 1c) which allows precipitation to reach the cave. Below the opening, covering the entire surface of the Great Hall, a layered ice block has developed, with an estimated thickness of 20 m and minimum volume of 30 000 m 3 (Orghidan et al., 1984;Brad et al., 2018). The descendent morphology of the cave and the presence of the two openings determine a specific type of air circulation (Pers , oiu and Onac, 2019): cold air inflow through the lower entrance and warm air outflow through the upper entrance in winter, and slow convective circulation within the cave (with no air mass exchange with the outside) during summer. As a result of this air circulation, the dynamics of air temperature inside and outside of the cave follow a similar pattern between October and April, whereas cave temperatures are stable at 0 • C between May and September, regardless of changes outside (Pers , oiu et al., 2007).
A direct consequence of the predominantly negative air temperatures in the cave is the genesis, accumulation and preservation of ice (Fig. 1b, c). During summer, infiltrating rainwater accumulates on top of the existing ice block to form a layer of water that is approximately 0-20 cm deep. Monitoring of air temperature in the cave has shown that air temperature drop below 0 • C outside the cave in early autumn (September), leading to cold air avalanches reaching the inner parts of the cave (Pers , oiu et al., 2007). The inflow of cold air leads to the freezing of the lake water from top to bottom, forming a 1-20 cm thick layer of ice ("summer ice"). Although the air temperature might briefly rise above 0 • C in autumn, temperatures inside the cave do not exceed 0 • C. Thus, the layer of ice formed on top of the lake will prevent the subsequent addition of water to the lake, thereby preserving the original isotopic composition of water reflecting that of summer precipitation. Infiltration and subsequent freezing of water during warm periods in winter result in additional layers of ice on top of the ice block ("winter ice"). However, at the onset of melting, this winter ice melts (Pers , oiu et al., 2011b). The result of these processes is a multiannual, layered, ice block, consisting of annual couplets of clear ice (on top) and a sediment-rich layer beneath. Inflow of warm water in wet summers leads to rapid ablation of the ice at the top of the ice block, partly altering the annual layering. The processes of cave ice formation by water freezing and is the digital elevation model (DEM), which was downloaded from http://geospatial.org (last access: 16 October 2020) and processed in ArcGIS by CAB. The map in panel (b) is our own work (Aurel Pers , oiu), and the map in panel (c) is modified from our previous work (Brad et al., 2018). the registration of environmental signals by various proxies (e.g., stable isotope composition of ice, pollen content) have been described from the nearby Scȃris , oara Ice Cave (Pers , oiu and Pazdur, 2011;Feurdean et al., 2011) and, given the similarities between the two caves, are also pertinent to Focul Viu Ice Cave. The one notable difference is the timing of the onset of freezing: in Scȃris , oara Ice Cave, the onset of freezing is delayed until late autumn and early winter (Pers , oiu et al., 2017), whereas in Focul Viu Ice Cave it starts in early autumn.

Drilling and stable isotope analyses
The FV ice core (4.87 m long and 10 cm diameter) was drilled in May 2016 from the Great Hall of FV (Fig. 1d) using a modified PICO electric drill (Koci and Kuivine, 1984) manufactured by Heavy Duties S.R.L, Cluj-Napoca, Romania. A rock embedded in the ice at 4.87 m below the surface stopped the drilling effort, but previous work in the cave has shown that the thickness of the ice block exceeds 15 m (Orghidan et al., 1984;Kern et al., 2004;Pers , oiu and Onac, 2019). The ice core was cut into 1 cm long pieces (taking the annual lay-ering into account), and each piece was subsequently sealed in plastic bags, allowed to melt at room temperature, transferred to 20 mL HDPE (high-density poly ethylene) scintillation vials and stored at 4 • C prior to analysis.
Water samples were analyzed for stable isotope composition at the Stable Isotope Laboratory,Ştefan cel Mare University (Suceava, Romania), using a Picarro L2130-i CRDS (cavity ring-down spectroscopy) analyzer connected to a high-precision vaporizing module. All samples were filtered through 0.45 µm nylon membranes before analysis and were manually injected into the vaporization module multiple times until the standard deviation of the last four injections was less than 0.03 for δ 18 O and less than 0.3 for δ 2 H. The average of these last four injections was normalized on the SMOW-SLAP scale using two internal standards calibrated against VSMOW2 and SLAP2 standards provided by the IAEA and further used in our interpretation. A third standard was used to check the long-term stability of the analyzer. The stable isotope values are reported using standard 2448 C.-A. Bȃdȃluţȃ et al.: Summer temperatures in east-central Europe are linked to AMO variability δ notation, and the precision is estimated to be better than 0.1 ‰ for δ 18 O and better than 0.5 ‰ for δ 2 H based on repeated measurements of an internal standard.

Radiocarbon dating and age-depth model construction
The wide opening to the surface in the ceiling of the Great Hall ( Fig. 1c) allows for a large volume of organic matter to fall into the cave and subsequently become trapped in the ice, including large pieces of wood that tend to cut through layers of ice, with the latter possibly encompassing decades of ice accumulation. Furthermore, ice melting and water freezing processes usually result in inclined ice surfaces and slow tipping of any heavy materials sitting on top of the ice (see the position of tree trunks in Fig. 1). As a consequence, wood with an age older than that of the newly forming ice can be incorporated in the ice block, resulting in sample ages much older than the ice layers ("old wood effect"). However, the challenge of identifying such organic material contrasts with the desire for precise chronologies, which rely on a high number of data points. Thus, all possible organic samples were recovered and dated, and unreliable ages that were identified at the stage of age-depth modeling were removed. Out of 14 samples recovered from the ice core and potentially suitable for radiocarbon dating, two were not datable due to their extremely small carbon yield. Accelerator mass spectrometry (AMS) radiocarbon analyses were performed at the Institute of Physics, Silesian University of Technology, Poland (Piotrowska, 2013). All samples were precleaned with standard acid-alkali-acid treatment, dried and subjected to graphite preparation using an AGE-3 system (Ionplus, Switzerland) equipped with an Elementar vario MICRO cube elemental analyzer and automated graphitization unit Nemec et al., 2010). The 14 C concentrations in graphite produced from unknown samples, Oxalic Acid II standards and coal blanks of comparable carbon masses were measured by the DirectAMS laboratory, Bothell, USA (Zoppi et al., 2007). The results are reported in Table 1. The radiocarbon dates were calibrated using Ox-Cal v4.3 (Bronk Ramsey, 2009) and the IntCal13 calibration curve (Reimer et al., 2013). The NH1 curve (Hua et al., 2013) was used for one post-bomb date. Because organic material can fall into the cave decades to centuries before being trapped in the ice (see Fig. 1b), we have carefully screened the radiocarbon results prior to agedepth modeling with the aim of selecting the most reliable dates forming a chronological sequence. In total, four dates were selected for age-depth modeling. For the top of the ice core, a uniform age distribution from 1991 to 2016 CE was assigned, allowing for the possibility of surface ice melting. The model was constructed using the OxCal P_Sequence algorithm (Bronk Ramsey, 2008) with a variable prior k parameter (k = 1, U (−2, 2); Bronk Ramsey and Lee, 2013) and extrapolated to a depth of 4.86 m. The agreement index of the model was 85 %, confirming a good statistical performance when the threshold of 60 % is surpassed. The sections between the dated depths were assumed to have a constant deposition rate. The complete age-depth model is shown in Fig. 2. The mean age derived from the model was used for further analysis and is also reported in Table 1. The constructed age-depth model was compared with that of Maggi et al. (2008), plotted in green in Fig. 2, and a broad agreement was found between both chronologies. All of our rejected ages are older than those of Maggi et al. (2008), and we suspect that these were based on dating old wood that had already been in the cave for decades before being incorporated in the ice (see the "old wood effect" discussion above and Fig. 1c). Further, the same authors identified several potential markers of volcanic eruptions, and their ages agree within ±40 years with those of our model. Another evaluation was obtained by comparing our stable isotope record with those of Kern et al. (2004) andForizs et al. (2004); although the latter lack a precise chronology, a simple visual correlation between the records indicates a satisfactory match. In order to avoid any circular reasoning, we decided against the use of other records, e.g., the regional temperature record, to further anchor our chronology.

Climate data
The sea surface temperature (SST) is extracted from version 5 of the Extended Reconstructed Sea Surface Temperature data (ERSSTv5) of Huang et al. (2018). This dataset covers the period from 1854 to present and has a spatial resolution of 2 • × 2 • . The AMO index used in this study was obtained from https://climexp.knmi.nl/data/iamo_ersst_ ts.dat (last access: 19 October 2020) and is also based on the ERSSTv5 dataset. Station-based meteorological data were provided by the Romanian National Meteorological Administration for three stations (Baia Mare, Sibiu and Timis , oara) that have some of the longest instrumental records in Romania and bracket the location of the study site. To remove the short-term variability and retain only the multidecadal signal in our data, both the temperature time series and the SST data were smoothed with a 21-year running mean filter prior to correlation analysis.

Ice accumulation in Focul Viu Ice Cave
The results of the radiocarbon analyses performed on organic matter recovered from the ice are shown in Table 1, and the age-depth model is displayed in Fig. 2. The maximum age of the ice is 1000 ± 20 cal BP at 4.45 m below surface, based on direct dating of organic remains (Table 1) and extrapolated to 1100 cal BP at 4.86 m below surface High accumulation rates were recorded between 850 and 950 CE (0.39-0.41 cm yr −1 ) and between 1220 and 1970 CE  (0.36-0.44 cm yr −1 ). Between 950 and 1220 CE, the net accumulation rate dropped to between 0.29 and 0.34 cm yr −1 .
The highest net accumulation rates recorded after 1970 CE (0.56 cm yr −1 ) contradict recent findings from other ice caves in the Carpathian Mountains (Kern and Pers , oiu, 2013), which all register record melting. However, this value might be an artifact of the age-depth modeling (see above) as well as of the very short time span considered; thus, it is unreliable for further interpretation. The low accumulation rates spanning the MWP are similar to those recorded in the Scȃris , oara Ice Cave (Pers , oiu et al., 2017, Bȃdȃluţȃ, 2019, the Hundsalm ice cave in Austria (Spötl et al., 2014) and the ice caves in the Velebit Mountains in Croatia (Kern et al., 2018), potentially suggesting a regional signal of climatic conditions unfavorable for ice accumulation. Ice can melt as a result of either warm summers with enhanced conductive heat transfer to the cave or wet summers with rapid ablation resulting from water flowing across the top of the ice block. Subsequently, ice growth is influenced by the amount of water present at the onset of freezing, the timing of this onset and its duration. The low accumulation rates during the MWP were likely the result of enhanced melting during warm (see Sect. 4.2 below) and wet (Feurdean et al., 2015) conditions. After 1450 CE, the climate in the region was dominated by dry summers with frequent storms and cold winters (Pers , oiu, 2017). These conditions led to reduced summer melting and enhanced winter growth, which are conditions favorable for net ice accumulation.

Stable isotopes in Focul Viu cave ice -a proxy for summer air temperatures and AMO variability
The variability of δ 18 O and δ 2 H in precipitation at Gheţar (∼ 13 km south of the cave's location and at the same altitude), assessed for the 2012-2017 CE period, follows that of temperature (Fig. 3a), with the maximum values (−3.6 ‰ and −26 ‰ for δ 18 O and δ 2 H, respectively) in July and August and the minimum values (−19.8 ‰ and −140 ‰ for δ 18 O and δ 2 H, respectively) in January. Similar results were found by Bojar et al. (2009) and Ersek et al. (2018) for the same region, suggesting that the 18 O / 16 O and 2 H / 1 H ratios in precipitation register temperature changes on a regional scale. The Local Meteoric Water Line, defined by the equation δ 18 O = 7.4 × δ 2 H + 6.1 (Fig. 3b), has a slope and intercept very similar to those found by Ersek et al. (2018). Precipitation in the region is mainly delivered by weather systems carrying moisture from the Atlantic Ocean, with the Mediterranean Sea contributing moisture during autumn and winter (Nagavciuc et al., 2019b). The deuterium excess in precipitation (d-excess or d), defined as d = δ 2 H − 8 × δ 18 O (Dansgaard, 1964) allows for a clear separation of the air masses: the Atlantic Ocean (d-excess close to the global average of 10; Craig, 1961) and the Mediterranean Sea (dexcess between 12 and 17, resulting from the high evaporative conditions in the eastern Mediterranean Sea). Similarly high values of d-excess have been measured in precipitation in southwestern Romania (Drȃgus , in et al., 2017) and in precipitation in northeastern Romania  and have been linked to air masses originating in the strongly evaporated Mediterranean and Black seas, respectively. Observations on the dynamics of cave ice during the past 18 years have shown that it starts to grow in early autumn due to the freezing of water accumulated during summer. As the ceiling of the cave is opened to the outside, precipitation directly reaches the site of ice formation; thus, the stable isotope composition of precipitation is not modified in the epikarst above the cave, and the original δ 18 O and δ 2 H values of summer (June-July-August, JJA) precipitation are preserved in the cave water. However, while freezing processes in caves could alter the original δ 18 O (and δ 2 H) values in cave ice, several studies have shown (e.g., in the nearby Scȃris , oara Ice Cave; Pers , oiu et al., 2011b) that the original climatic signal embedded in the stable isotope composition of cave ice is preserved and can be used as a proxy for external climate variability.
Overall, our observations of cave ice genesis and dynamics and the stable isotope monitoring data clearly indicate that summer air temperatures are registered and preserved in the ice block in the FV. In order to test the long-term stability of these connections, we have analyzed the links between the FV δ 18 O ice record and instrumental data from three nearby meteorological stations over the 1851-2016 CE period. On multidecadal timescales, summer air temperature changes in the region are mainly controlled by the dynamics of the Atlantic Multidecadal Oscillation (Ionita et al., 2012). Figure 4 shows the JJA air temperature at the Baia Mare (BM), Sibiu (SB) and Timis , oara (TM) stations, the AMO index and the FV δ 18 O ice . The instrumental temperature data indicate large multidecadal variability, with a cold period between 1890 and 1920 CE, followed by a warm period between 1921 and 1960 CE, a slightly colder period between 1960 and 1980 CE and enhanced warming after 1980 CE, which all follow the AMO variability. The δ 18 O ice values show a similar temporal evolution, with the slight offsets between the observational data and δ 18 O likely being due to the dating uncertainty (20-35 years). Further, we have computed the correlation map between the summer mean air temperature at SB station (with the longest instrumental record) and the summer SST as indicator of AMO variability (Sutton and Dong, 2012). Figure 5 clearly shows that positive (negative) temperature anomalies over the analyzed region are associated with positive (negative) SST anomalies over the North Atlantic Ocean, resembling the SST anomalies associated with the positive (negative) phase of the AMO (Mesta-Nuñez and Enfield 1999; Latif et al., 2004;Knight et al., 2005). The strongest corre- lations (Fig. 5) are found between local air temperature and SSTs in the North Atlantic and the eastern Mediterranean Sea, which are the main sources of moisture feeding local precipitation . These results are also in agreement with the results of Della-Marta et al. (2007), who showed that summer positive temperature anomalies and heat waves over Europe are triggered, at least partially, by the phase of the AMO. A recent study of δ 18 O variability in oak tree rings in northwestern Romania (∼ 75 km north of our site) also indicated the influence of the AMO on summer temperatures and drought conditions (Nagavciuc et al., 2019a). Thus, we suggest that, during summer, strongly meandering Rossby waves (Ionita et al., 2015 result in blocking conditions over central Europe that lead to the persistence of high-pressure systems and the occurrence of regional heat waves. These, in turn, favor regional recycling of moisture, resulting in positive δ 18 O anomalies in precipita-tion that are further recorded by climate proxies: stable isotopes in cave ice (this study) and tree ring parameters (width and cellulose stable isotopes) (Popa and Kern, 2009;Nagavciuc et al., 2019a). The increased contribution of recycled moisture to precipitation has been reported from several locations in central Europe, further supporting our inference (Goìmez-Hernaìndez et al., 2013;Kern et al., 2020).
Combining all of the abovementioned data, we find that, on timescales ranging from years to decades, prolonged periods of positive temperature anomalies throughout the summer months, linked to prolonged warm SSTs in the North Atlantic Ocean (and thus a positive AMO index), could be preserved by the δ 18 O ice in FV.
The FV δ 18 O ice and δ 2 H ice records span the 850-2016 CE period (Fig. 6). Decadal-to multidecadal-scale oscillations occur over the entire record, but no discernable long-term trend has been identified. Several periods of excursions towards low δ 18 O values (defined as δ 18 O ice below the long-term average), indicating low summer temperatures have been observed at 875-930 CE, 1050-1080CE, 1260-1330CE, 1430-1480CE, 1520-1550CE, 1710-1750CE, 1820-1870CE and 1880-1930 Significant maxima in the FV δ 18 O ice record occurred during 850-870 CE, 1000-1050CE, 1080-1260CE, 1350-1390CE, 1480-1520CE, 1625-1710CE and 1950-1970). Both the δ 18 O ice and δ 2 H ice records display a remarkable similarity throughout the entire period; thus, we have relied solely on the δ 18 O record in our discussion. We have compared the FV δ 18 O ice with a tree ring width-based reconstruction of summer (JJA) temperature anomalies from the Eastern Carpathian Mountains (Popa and Kern, 2009). The highest similarities between the FV ice core and summer temperature records were found for the cold periods during 1260-1330CE, 1430-1480CE, 1520-1550CE and 1820-1870CE and the warm periods during 1080-1260CE, 1625-1710CE and 1950-1970 (Fig. 6). Given the very different nature of the two archives (trees vs. cave ice), of the proxies (tree ring width and δ 18 O) and of the chronologies (annual tree ring counting vs. 14 C dating with a ±30-year error), the two records agree remarkably well, further supporting the hypothesis that δ 18 O and δ 2 H values in the FV ice core register both summer air temperature variability during the past ca. 1000 years in eastcentral Europe and, on a broader spatial scale, the variability of the AMO.
Similar to the 1850-2016 CE interval described above, the stable isotope record closely mirrors the AMO variability over the entire studied interval (Fig. 7). The relation-ship between the FV δ 18 O ice and AMO records (Wang et al., 2017) is strongest between 1125 and 1525 CE and 1750 and 2016 CE, with decadal-scale variability in the two records being synchronous. Peak-to-peak matching of the FV δ 18 O ice and AMO records shows that the two records correlate well within the dating uncertainty (±30 years) during the past 1000 years. However, peak-to-peak matching indicates a difference between the two records of up to ∼ 50 years between ∼ 1600 and 1750 CE (Fig. 7), which is likely the result of high (> 50 years) uncertainty in the ice core chronology between 1525 and 1750 CE (Fig. 7). The combined correlations between (1) FV δ 18 O ice and instrumental summer temperature reconstruction over the past 150 years, (2) the instrumental (this study) and proxy-based (Nagavciuc et al., 2019a) summer air temperatures and AMO variability, and (3) the FV δ 18 O ice and AMO records over the past 1000 years suggest that changes in the North Atlantic are transferred, likely via atmospheric processes, towards the wider Northern Hemisphere, resulting in hemisphere-wide climatic responses to these changes.
The FV δ 18 O ice record is in agreement with other summer temperature reconstructions (e.g., Buntgen et al., 2011) at both regional and hemispheric scales (Fig. 6). Further, regional summer temperature (e.g., Popa and Kern, 2009) and summer temperature-sensitive drought (Seim et al., 2012) reconstructions show warm peaks around 1320 CE, 1420 CE, 1560 CE and 1780 CE and cold ones at around 1260 CE, 1450 CE and 1820 CE, similar to reconstructions and models at the global level (Neukom et al., 2019) and the FV temperature reconstruction (this study). Contrary to the summer season temperature reconstructions, a late-autumn through early-winter season temperature reconstruction from the nearby Scȃris , oara Ice Cave (Pers , oiu et al., 2017) shows that the MWP was rather warm and also wet (Feurdean et al., 2011), whereas the LIA was cold and likely dry, with erratically distributed precipitation. Together, these data suggest a complex picture of climate variability in the wider Carpathian region, with much of the yearly temperature variability during the past 1000 years being attributed to the influence of winter conditions.

Conclusions
The analysis of the oxygen and hydrogen stable isotope ratios along a ∼ 5 m long ice core extracted from the Focul Viu Ice Cave (northwest Romania) provided an unprecedented view on the dynamics of summer air temperature and atmospheric circulation changes during the past 1000 years in east-central Europe. The comparison of ice core δ 18 O (and δ 2 H) with instrumental data over the past 150 years indicates that the stable isotope composition of cave ice records the summer air temperatures on multidecadal timescales. Given the apparent relationship between our stable isotope record and both instrumental and proxy-derived records of the At-   (Wang et al., 2017), and the 14 C measurement uncertainty between 850 and 2000 CE. Shading indicates the offset (in years) between the FV δ 18 O and AMO index values: green denotes less than 20 years, yellow denotes between 20 and 50 years and orange denotes above 50 years.
lantic Multidecadal Oscillation (AMO), we further suggest that δ 18 O ice in the Focul Viu Ice Cave can be used to infer past AMO variability. We subsequently hypothesize that changes in summer climatic conditions over the Northern Atlantic are transferred through atmospheric processes across the Northern Hemisphere, influencing summer temperatures across Europe. However, we emphasize that given the ±30year uncertainty in the chronology, detailed studies using better constrained age models (e.g., tree-ring-based proxies) are required to test this hypothesis.
The data show little centennial-scale summer temperature variability since the onset of the Medieval Warm Period and through the Little Ice Age. Nevertheless, well-expressed minima and maxima occurred synchronously with data from other records in the region and through Europe, suggesting that the stable isotope composition from cave ice records a regional climatic signal. In contrast, winter air temperature records from the region indicate colder conditions during the LIA compared with the MWP, pointing towards a seasonally distinct climatic signal during these two periods. This suggests that forcing factors acting seasonally had a strong imprint on temperature variability, overriding long-term, global forcing.
Our results offer a potential hypothesis to be further tested by extending this and similar records back in time and also incorporating other proxy-base reconstructions to investigate the spatial extent of the influence of North Atlantic climate further east. Data availability. The Focul Viu δ 18 O and δ 2 H data as well as the 14 C data and the modeled ages will be made available on the NOAA World Data Service for Paleoclimatology web page (https://www.ncdc.noaa.gov/paleo-search/, last access: 9 December 2020). The meteorological data plotted in Fig. 4 were provided by the Romanian National Meteorological Administration, except for the AMO data that were downloaded from https://climexp.knmi. nl/data/iamo_ersst_ts.dat (last access: 7 December 2020), which are created and maintained by the Royal Netherlands Meteorological Institute (https://www.knmi.nl, last access: 7 December 2020). The paleoclimate data used to plot Fig. 6a were downloaded from the NOAA World Data Service for Paleoclimatology web page (https://www.ncdc.noaa.gov/paleo-search/). Author contributions. CAB and AP designed the project, AP and CAB collected the ice core, and CAB performed the stable isotope analyses. NP performed the radiocarbon analyses and constructed the age-depth model. MI analyzed the climate and large-scale circulation data. CAB and AP wrote the text, with input from MI and NP.