*Article* **Isotopic 'Altitude' and 'Continental' E**ff**ects in Modern Precipitation across the Adriatic–Pannonian Region**

#### **Zoltán Kern 1,\*, István Gábor Hatvani 1, György Czuppon 1,2, István Fórizs 1, Dániel Erdélyi 1,3, Tjaša Kanduˇc 4, László Palcsu <sup>2</sup> and Polona Vreˇca <sup>4</sup>**


Received: 17 April 2020; Accepted: 19 June 2020; Published: 24 June 2020

**Abstract:** It is generally observed that precipitation is gradually depleted in 18O and 2H isotopes as elevation increases ('altitude' effect) or when moving inland from seacoasts ('continental' effect); the regionally accurate estimation of these large-scale effects is important in isotope hydrological or paleoclimatological applications. Nevertheless, seasonal and spatial differences should be considered. Stable isotope composition of monthly precipitation fallen between January 2016 and December 2018 was studied for selected stations situated along an elevation transect and a continental transect in order to assess the isotopic 'altitude' and 'continental' effects in modern precipitation across the Adriatic–Pannonian region. Isotopic characteristics argue that the main driver of the apparent vertical depletion of precipitation in heavy stable isotopes is different in summer (raindrop evaporation) and winter (condensation), although, there is no significant difference in the resulting 'altitude' effect. Specifically, an 'altitude' effect of <sup>−</sup>1.2%/km for <sup>δ</sup>18O and <sup>−</sup>7.9%/km for <sup>δ</sup>2H can be used in modern precipitation across the Adriatic–Pannonian region. Isotopic characteristics of monthly precipitation showed seasonally different patterns and suggest different isotope hydrometeorological regimes along the continental transect. While no significant decrease was found in δ18O data moving inland from the Adriatic from May to August of the year, a clear decreasing trend was found in precipitation fallen during the colder season of the year (October to March) up to a break at ~400 km inland from the Adriatic coast. The estimated mean isotopic 'continental' effect for the colder season precipitation is <sup>−</sup>2.4%/100 km in <sup>δ</sup>18O and <sup>−</sup>20%/100 km in <sup>δ</sup>2H. A prevailing influence of the Mediterranean moisture in the colder season is detected up to this breakpoint, while the break in the δ18O data probably reflects the mixture of moisture sources with different isotopic characteristics. A sharp drop in the d-excess (>3%) at the break in precipitation δ18O trend likely indicates a sudden switch from the Mediterranean moisture domain to additional (mainly Atlantic) influence, while a gradual change in the d-excess values might suggest a gradual increase of the non-Mediterranean moisture contribution along the transect.

**Keywords:** stable isotopes; oxygen; hydrogen; d-excess; elevation effect; altitude effect; continental effect; Slovenia; Hungary

#### **1. Introduction**

Precipitation is a key element of the atmospheric branch of the water cycle, and the study of precipitation is essential for improving the understanding of the Earth's water balance on both regional and global scales [1]. An important method for tracing the water cycle is that of comparing the ratios of heavy to light stable isotopes of hydrogen and oxygen (2H/ 1H; 18O/ 16O) in the water molecule [2–4]. The ratio of heavy to light stable isotopes is traditionally expressed as per mil (%) deviation relative to Vienna Standard Mean Ocean Water (VSMOW), and delta notation (δ18O and δ2H) is commonly used to report the measured isotope variations [5].

It is generally observed that precipitation is gradually depleted in 18O and 2H isotopes as elevation increases [6]. The phenomenon is commonly called the 'altitude' effect and primarily results from the cooling of the air masses as they ascend on a mountain, accompanied by the rainout of the excess moisture [7,8]. In case of precipitation <sup>δ</sup>18O a gradient with elevation has a global average of <sup>−</sup>2.8%/km, ranging from −1.7 to −5.0%/km; the European average is −2.1%/km [9]. A similar general observation is that precipitation is gradually depleted in 18O and 2H isotopes when moving inland from seacoasts. This phenomenon is called the 'continental' effect, and it reflects the gradual depletion of the residual marine moisture in heavy isotopes during the sequential rainout as air masses travel from the marine moisture source to the center of the continents [6]. In the case of Europe, for example, a gradient in the long-term mean precipitation δ18O shows an ~8% decrease over a distance of approximately 4000 km from the Irish coasts to the foothills of the Ural Mountains [10]. Besides the above-mentioned ones, there are a few other empirical relationships observed between environmental parameters and the stable hydrogen and oxygen isotope composition of precipitation, such as 'temperature', 'latitude' and 'amount' effects [6,10].

These large-scale effects are also important from a practical point of view because they have enabled useful applications in isotope hydrology (e.g., identification of the mean elevation of groundwater recharge [11,12]), paleoclimatology [13], paleoaltimetry [14], food authenticity [15], etc. However, numerous studies have reported seasonal and spatial differences in both 'altitude' [16–19] and 'continental' effects [10,20,21]. Therefore, the regional estimation of these empirical relations in modern precipitation is critically important for underpinning the application of appropriate gradient values to subcontinental scale (e.g., river basin) practical problems (e.g., estimation of mean recharge area [22]) or paleoclimatological comparisons [13].

The aim of the study is to assess the isotopic 'altitude' and 'continental' effects in modern precipitation across the Adriatic–Pannonian region, specifically to (i) check the potential seasonal changes and (ii) determine the characteristic gradients specifically for this region which can be applied in isotope hydrological/hydrogeological studies and other research in the future.

#### **2. Study Area and Main Moisture Sources**

The Adriatic–Pannonian region encompasses the northern part of Southeast Europe and the eastern part of Central Europe. The Mediterranean Sea was identified as the main marine moisture source for Central and Southeastern Europe [23–25], and changes in its moisture supply have a decisive role on dryness conditions in this region [26]. More specifically, the Western and Central Mediterranean are the dominant marine moisture source regions throughout the year, although locally recycled continental moisture is the dominant source of the atmospheric moisture in the summer [24]. At the northern part of the study region, the westerlies transport Atlantic marine moisture mixed with recycled terrestrial moisture from Western and West-Central Europe, in addition to the Mediterranean moisture contribution [27,28]. Studies dealing with single precipitation events in Eastern Hungarian localities found occasional moisture contributions from eastern and northern directions mostly in the October to March period [27,28].

The surface water of the Mediterranean Sea is more enriched in heavy isotopes than Atlantic surface water [29]. This feature and the distinct evaporative conditions in these regions imprint the resultant atmospheric moisture with a source-specific isotopic signature (e.g., see [30,31]; more details in Section 4.2).

#### **3. Data and Methods**

Stable hydrogen and oxygen isotope compositions (δ2H and δ18O) of monthly aggregated precipitation were gathered from the national precipitation monitoring networks of Austria [32], Slovenia [33,34] and Hungary [35] and from additional individual stations with accessible data, namely Debrecen [36] and Hraš´cica [37], from three consecutive years (2016, 2017, 2018). Subsets suitable for assessing the 'altitude' and 'continental' effects have been selected from the gathered data set (Figure 1). The sampling methods might be different among the stations, but they always employed a collector configuration recommended for precipitation isotope analysis [38]. The analytical techniques may also differ between data sources (e.g., mass spectrometry for Slovenian samples, laser spectroscopy for the Hungarian and Croatian samples); however, the raw values were converted by all laboratories to the same international reference scale (VSMOW/SLAP), ensuring comparability and joint assessment of the data. In addition, measurements were carried out together with laboratory reference materials that are calibrated periodically against primary IAEA calibration standards to VSMOW/SLAP scale. The highest estimated uncertainties of the data are <sup>±</sup>1% and <sup>±</sup>0.3% for <sup>δ</sup>2H and <sup>δ</sup>18O, respectively.

**Figure 1.** The study area and the location of the precipitation stable isotope monitoring stations. The sites used in the evaluation of the 'continental' effect are indicated by white dots, while the ones used in studying the 'altitude' effect are indicated by black dots; Ljubljana was used in both. Site details can be found in Table 1. The black rectangle in the inset map of Europe shows the position of the studied region. Country map of Europe by FreeVectorMaps.com. Digital elevation model taken from http://srtm.csi.cgiar.org/.

Stations situated within a relatively close distance but spanning a pronounced elevation range (from 282 to 2515 m a.s.l.) around the Slovenian–Austrian border (Figure 1 and Table 1) offer the opportunity to estimate the regional 'altitude' effect in modern precipitation in the Adriatic–Pannonian region. Due to the restricted latitudinal extension of this elevation transect, latitudinal or continental dependence are assumed to have only a negligible effect on the stable isotope composition of precipitation and a similarly uniform moisture source can be plausibly assumed for these closely located stations. The elevation range of this transect practically comprises the elevation extent of

the entire Adriatic–Pannonian region, so the derived empirical relation should be applicable to the entire region.

To study the continental effect, stations situated along a ~750 km transect, stretching in a SW–NE direction from the Adriatic coast (Istria) to the NE corner of the Great Hungarian Plain, were selected (Figure 1 and Table 1). To eliminate potential interference with 'altitude' effect, only stations below 350 m a.s.l. were considered. This transect was chosen due to the predominance of the Mediterranean as the moisture source over the studied region [25,26].

**Table 1.** Basic geographical information of precipitation stable isotope monitoring stations operating in East-Central Europe between January 2016 and December 2018. Country codes follow the ISO-3166-1 ALPHA-2 standard. Latitude and longitude are in WGS84 projection; EPSG: 4326.


To statistically investigate the 'altitude' and 'continental' effects, the monthly precipitation stable isotope values were plotted against station elevation and distance from the Adriatic coast, respectively. Next, linear trendlines were fitted with least squares regression and the fits' statistics were investigated (slope (coefficients), *p*, adjusted *r*2). Only the slope coefficients with their significance at the usually applied significance thresholds are presented. Besides the overall trends, in the seasonal (winter and summer half-years) difference in 'altitude' effect was explored with an independent sample Mann–Whitney U-test (MW U-test) [39].

Deuterium excess (d-excess), defined as d <sup>=</sup> <sup>δ</sup>2H <sup>−</sup> <sup>8</sup> <sup>×</sup> <sup>δ</sup>18O [6] was also calculated. The uncertainty of d-excess is ≤1.3%, as propagated from the highest uncertainties of the primary isotopic parameters. This second-order isotopic parameter is a useful tracer of the origin of moisture because it is relatively invariant during transport and during the formation of condensate accompanied by equilibrium isotopic fractionation [40] but is an informative indicator of non-equilibrium isotopic effects [6] and the mixing of moistures of different origin [41].

#### **4. Results and Discussion**

In δ*-*δ space, the global meteoric water line [42] fits nicely to the used data (Figure S1), certifying that the isotopic compositions correspond to precipitation origin. Due to the strong linear relationship between δ2H and δ18O in precipitation (Figure S1), the results for 'altitude' and 'continental' effects on the monthly scale are only presented graphically for δ18O for the studied 36-month period between 2016 January and 2018 December.

#### *4.1. 'Altitude' E*ff*ect*

The analysis of the precipitation δ18O gradient indicated a decreasing trend in most of the months, ranging from 0.1 to <sup>−</sup>3.5%/km (Figure 2a). The distribution of these monthly <sup>δ</sup>18O lapse rates tends to be higher than the reported global range (from −1.7 to −5%/km [9]), although overlapping with the upper portion of the global range, and corresponds well with the range of other large-scale

compilations [8]. The fitted linear slope was significant in roughly half of the cases, and most of the insignificant (*p* > 0.1) trends were observed in the October–March half-year (59% in the case of δ18O, 62.5% in the case of δ2H).

**Figure 2.** Stable isotope characteristics of monthly precipitation along a representative elevation transect in the Adriatic–Pannonian region: δ18O (**a**) and deuterium excess (d-excess) (**b**) vs. elevation. Best fit lines (blue) for precipitation δ18O data and elevation are plotted. The slope and significance of the linear fit are coded as α-values of 0.01, 0.05 or 0.1 and are marked with (\*\*\*), (\*\*) or (\*), respectively.

Corresponding d-excess values along the elevation transect show a different vertical pattern when comparing colder (October to March) and warmer (April to September) months (Figure 2b). While d-excess values are scattered in a relatively narrow range without any systematic shift along the elevation gradient in the colder season, steep gradients with low d-excess values for low elevation stations can be seen in the warmer season (Figure 2b). This suggests distinct isotopic precipitation regimes for the warmer and colder parts of the year.

According to the classical concept of Rayleigh rainout, process whatever isotopic properties acquired by the cloud before, precipitation is removed under isotopic equilibrium conditions [6]. Under this ideal Rayleigh rainout scenario, the atmospheric moisture is depleted in heavy isotope content commensurate with the loss of moisture by precipitation; the d-excess during this process remains invariant [43], even though there is a tendency for an increase in this value at higher elevations [8].

The remarkably constant d-excess data along the elevation transect suggests a single dominant moisture source during the colder season [19]. Thus, in winter conditions, Rayleigh rainout effect can be taken as the predominant one, with a single dominating moisture source and with the moisture mass ascending to form orographic precipitation which is continuously depleted in heavy isotopes with increasing elevation [8,19]. However, the change of d-excess values along the elevation transect in the warmer season suggests that other factors which change the isotope composition besides the basic Rayleigh effect need to be considered. Evaporative enrichment of isotopes in raindrops during the fall beneath the cloud base reduces the d-excess [6,44]. The amplification of this effect is indeed expected in the warmer season; and it is expected to be larger at low elevations, conforming to the observed vertical d-excess trend (Figure 2b). Therefore, this appears to be a 'pseudo-altitude effect' [45] when vertical change of stable isotopes in precipitation is not due to the rainout of vapor of a single origin [8].

It should be noted that precipitation is dominantly rain in the warmer season over the studied elevation transect; while in the colder season, only snow is expected above ~1000 m, and a mixture of snow and rain is expected in the lower part of the transect. Since vapor → ice fractionation is slightly larger than vapor → water fractionation around the altitude where the regime switches from rain to snow (roughly corresponding to the 0 ◦C isotherm), the condensed hydrometeor is expected to be less depleted in heavy isotopes than it would be if the vapor → water fractionation was still in action. This phenomenon might be an additional effect that influences the monotony of the vertical upward decrease in precipitation δ18O and causes the fitted linear regression to be frequently non-significant in the colder season.

Comparing the frequency distribution of monthly estimated 'altitude' effect for the colder and warmer seasons for the 2016–2018 period (Figure 3), no significant difference is found in the seasonal distributions (*p* = 0.351 for δ18O and *p* = 0.832 for δ2H according to the MW U-Test) despite the undeniably different responsible processes dictating the vertical trend of precipitation stable isotopes in the warmer and colder seasons. Hence, from a practical point of view, the mean of the monthly gradients calculated for the studied 3-year-long period (δ18O: <sup>−</sup>1.2%/km; <sup>δ</sup>2H: <sup>−</sup>7.9%/km) can be considered as empirical isotopic 'altitude' effect in modern precipitation across the Adriatic–Pannonian region.

The estimated 'altitude' effect based on spring water <sup>δ</sup>18O values (−1.1%/km [22]) investigated between 2005 and 2007 in the vicinity of the studied elevation transect (in the valley of the Radovna River in NW Slovenia) is in good agreement with the current results. Comparing these results to former studies from the region, which typically employed fewer stations and/or a narrower elevation range, one can say that the magnitude of the determined isotopic 'altitude' effect tends to be smaller than the <sup>δ</sup>18O vertical gradient calculated to the south of the study area, e.g., near Zagreb (−2.8%/km) [46] or in the northern Adriatic (−3.0%/km) [47]. However, at the northern border of the study area, in the Western Carpathians, the 'altitude' effect estimated both from spring water <sup>δ</sup>18O values (−1.0 to <sup>−</sup>1.5%/km) [48] and from precipitation <sup>δ</sup>18O values (−2.1%/km) [49] is relatively closer to the empirical δ18O lapse rate obtained here for modern precipitation. At the eastern border of the study area, in the Eastern Carpathians, the estimated long-term 'altitude' effect based on spring water δ18O values also returned comparable results (−1.5%/km [50]).

**Figure 3.** Frequency distribution of 'altitude' effects estimated from monthly precipitation δ18O (**a**) and δ2H (**b**) data for colder (October to March) and warmer (April to September) seasons.

#### *4.2. 'Continental' E*ff*ect*

As an initial observation, one can see that the seasonal change in precipitation stable isotopes is well reflected by the shifting of the whole data set, with more negative δ18O values seen in winter and less negative δ18O values in summer in each studied year. The monthly patterns in precipitation δ18O values along the transect streching from the Adriatic coast to the northeastern corner of the Great Hungarian Plain did not indicate any trend in precipitation δ18O data from May to August (Figure 4a). However, a generally negative gradient is outlined for the dataset from April to September, resembling an expected 'continental' effect.

A similar seasonal contrast has been reported for the long-term mean δ2H gradient along the Atlantic stormtrack in Western Europe and was explained by enhanced raindrop evaporation in summertime at inland locations [20]. Occasionally, a remarkable shift (~10%) can be seen (e.g., June and July 2017) toward lower d-excess values for single inland stations, which might indicate a kinetic isotope effect such as raindrop evaporation (Figure 4b). However, the situation seems to be more complicated here because usually there is no visible shift toward lower d-excess values during the summer months (see e.g., 2018) or, at most, a wide scattering (e.g., May 2016) can be observed (Figure 4b). The absence of a 'continental' effect was reported over the Amazon Basin [51], which was explained by the intense terrestrial moisture recycling [52]. In the studied region, summer is the high season for vegetation activity, accompanied by the highest transpiration flux within the year. Transpiration is essentially a non-fractionation process [53] and thus re-introduces isotopically unchanged water to the atmosphere through evapotranspiration. The large-scale dominance of transpired moisture might explain the coupled lack of trend of δ18O and d-excess values, while the variable importance of evaporation and transpiration might spoil the empirical isotopic 'continental' effect and cause an accompanying wide scattering in d-excess values in precipitation. Supporting this explanation, Lagrangian moisture source modeling clearly indicated the dominant contribution of locally recycled continental moisture in the summer precipitation for Southeastern Europe [24].

A strikingly different pattern is discernible from September to April. Precipitation δ18O values remarkably decrease with increasing distance from the Adriatic coast over the first part of the studied profile, with a break in this trend at ~400 km (Figure 4a). The δ18O values beyond this point are not as depleted as would be expected from the linear trend fitted to the data of the first segment of the profile. The largest gradients were observed from November to March, while the estimated effect is smaller and usually not significant (*p* > 0.1) in April and September (Figure 4a). It is likely that the transition between the winter regime and the summer regime discussed above occurs during these months. An apparent nonlinearity of the 'continental' effect in the winter season was also reported along the distillation route of Atlantic moisture, with a gradient larger by a factor of two from the Irish coast to Central Europe compared to the rest of the transect to the foothills of the Ural Mountains. [10]. The decreased gradient was explained by additional moisture supply (e.g., Black Sea, Caspian Sea) substantially contributing to the water balance over Eastern Europe.

**Figure 4.** Stable isotope characteristics of monthly precipitation from the Adriatic coast to the northeastern corner of the Great Hungarian Plain: δ18O (**a**) and deuterium excess (d-excess) (**b**) vs. distance from the Adriatic coast. Best fit lines (blue) for precipitation δ18O data and distance are plotted over the first half of the continental transect. The significance of the linear fit is coded as α-values of 0.01, 0.05 or 0.1 that are marked with (\*\*\*), (\*\*) or (\*), respectively.

Monthly d-excess patterns along the continetal transect across the Adriatic–Pannonian region also suggest additional moisture supply for the winter precipitation over the most inland section of the transect. Usually a shift can be observed in the d-excess values at the place where the break appears in the decrease of the δ18O values (Figure 4b). Deuterium-excess values are grouped around a relatively higher mean for the first part of the transect accompanying the prevailing δ18O gradient, while the d-excess values of the stations drop and usually group around a lower mean beyond this break. Characteristic examples are November 2016, December 2017 and February 2018, when precipitation d-excess values along the first part of the profile are typically at least 3% higher compared to the last part of the profile (Figure 4b). The average October-to-March d-excess value of the stations with a distance < 350 km to the Adriatic coast is 1.2%, higher compared to the average for those stations situated further than 450 km from the sea.

This isotopic pattern agrees with the results of Lagrangian moisture source modeling for the winter (wet) season in pointing out the Western and Central Mediterranean as the dominant moisture source regions in Southeast Europe [24], while the northern sector of the study area is connected to the Northern Alps where the North Atlantic was found to be the primary marine moisture source [23]. Water vapor originating from the Western and Central Mediterranean has characteristically higher d-excess values [30] compared to the water vapor of the Atlantic atmospheric marine boundary layer [31]. Therefore, it is likely that the lower d-excess values in the northern part of the study area confirm the influence of the Atlantic moisture supply. Indeed, studies on the stable isotope composition of precipitation in Hungary (representing the inland sector of the studied continental transect; see Figure 1) found characteristically higher d-excess values in the case of precipitation arriving from the Mediterranean compared to precipitation arriving from the Atlantic directions [27,28]. The break in the inland-trend of precipitation δ18O values and the accompanying sharp shift in the d-excess values probably indicate a sudden switch from the Mediterranean dominated regime to the Atlantic influence, while the gradual change in d-excess values (e.g., January 2016, February 2017) suggests gradual mixing of moisture delivered from both sources (Figure 4). Occasional contribution from additional moisture sources (e.g., from eastern and northern directions) between October and March can further complicate the isotope hydrometeorological situation in the NE sector of the Great Hungarian Plain [27,28].

Taking all of these results into consideration, a classical isotopic 'continental' effect can be considered only up to ~400 km inland from the Adriatic coast, and the mean gradients estimated from the monthly <sup>δ</sup>18O and <sup>δ</sup>2H data are <sup>−</sup>2.4 and <sup>−</sup>20 %/100 km, respectively. Comparing these results to former studies, we can conclude that this regional gradient is substantially larger than that estimated for Western Europe (−0.38%/100 km) from the Irish coast to Central Europe for <sup>δ</sup>18O [10] and the <sup>−</sup>3.3%/100 km estimated for <sup>δ</sup>2H [20]. This observation of a much steeper gradient along the Adriatic moisture source pathway can be explained by the stronger temperature gradient between the Adriatic coast and the center of the Pannonian Basin (i.e., mild Mediterranean winter vs. cold continental winter) compared to the temperature gradient prevailing over Western Europe in the winter season. The stronger temperature gradient likely forces a higher degree of condensation and a higher degree of fractionation over the same distance.

#### **5. Conclusions**

Assessing the stable isotope composition of monthly precipitation fallen between January 2016 and December 2018 for selected stations situated along an elevation transect and a continental transect provided empirical estimates for the isotopic 'altitude' and 'continental' effects in modern precipitation across the Adriatic–Pannonian region.

Isotopic characteristics suggest that the main driver of the apparent vertical depletion of heavy stable isotopes in precipitation differs in warmer (raindrop evaporation) and colder (condensation) seasons, although there is no significant difference in the resulting 'altitude' effect. Therefore, an 'altitude' effect of <sup>−</sup>1.2%/km for <sup>δ</sup>18O and <sup>−</sup>7.9%/km for <sup>δ</sup>2H can be recommended to be used for modern precipitation in the Adriatic–Pannonian region.

The monthly patterns in precipitation δ18O values along a ~750 km long transect stretching from the Adriatic coast to the norteastern corner of the Great Hungarian Plain did not indicate any significant trend in precipitation δ18O data from May to August. This can be explained by the overwhelming influence of terrestrial recycled moisture by evapotranspiration. The recycled moisture negates the effect of the classical Rayleigh distillation process, erasing the expected decreasing trend from the precipitation stable isotope data moving inland from the Adriatic coast. A significant 'continental' effect emerged in the stable isotope composition of winter precipitation along the same transect. However, isotopic characteristics argue for the prevailing influence of the Mediterranean moisture source in winter precipitation only up to ~400 km distance inland from the Adriatic coast. The classical concept of the 'continental' effect is applicable practically from the Adriatic coast to the Danube, and the estimated mean isotopic 'continental' effect for the winter precipitation is <sup>−</sup>2.4%/100 km in <sup>δ</sup>18O and <sup>−</sup>20%/100 km in <sup>δ</sup>2H. A break in the spatial trend of isotopic characteristics in winter precipitation at ~400 km probably reflects the mixture of moisture sources that travel along different distillation routes and delivering moisture with different isotopic composition. A sharp drop in the d-excess (>3%) corresponding in space with the break in precipitation δ18O trend likely indicates a sudden switch from the Mediterranean moisture domain to additional (mainly Atlantic) influence, while a gradual change in the d-excess values might suggest a gradual increase of the non-Mediterranean moisture contribution along the transect.

The regionally determined isotopic 'altitude' and 'continental' effects of modern precipitation can be recommended for use in future isotope hydrological or paleoclimatological applications in the Adriatic–Pannonian region.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4441/12/6/1797/s1, Figure S1: Crossplot of the δ18O vs. δ2H values used in the present study, Table S1: Stable isotopic data used in the study.

**Author Contributions:** Conceptualization, Z.K.; methodology, Z.K. and I.G.H.; formal analysis, I.G.H. and Z.K.; resources, G.C., L.P., P.V. and T.K.; data curation, P.V. and D.E.; writing—original draft preparation, Z.K. and I.G.H.; writing—review and editing, P.V., I.F. and G.C.; visualization, I.G.H., D.E. and Z.K.; project administration, Z.K. and P.V.; funding acquisition, Z.K., P.V. and L.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the National Research, Development and Innovation Office under Grants SNN118205 and PD121387 and by the Slovenian Research Agency ARRS under Grants N1-0054, J4-8216 and P1-0143. The research was partly supported by the European Union and the State of Hungary, co-financed by the European Regional Development Fund in the project of GINOP-2.3.2-15-2016-00009 'ICER'.

**Acknowledgments:** The authors express their gratitude to the staff of the Slovenian Environment Agency at the meteorological stations and I. Kanduˇc for their valuable help in sampling and to S. Žigon for assistance with laboratory work. The results of this study have been discussed within the COST Action: "WATSON" CA19120. The authors applied the FLAE approach for the sequence of authors. See https://doi.org/10.1371/journal.pbio. 0050018 for further details. This is contribution No. 73 of the 2ka Paleoclimatology Research Group.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Water* Editorial Office E-mail: water@mdpi.com www.mdpi.com/journal/water

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18

www.mdpi.com ISBN 978-3-03943-267-7