Climate warming has arrived in caves of the European Alps

Climate warming has arrived in caves of the European Alps | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Climate warming has arrived in caves of the European Alps Friedrich Obleitner, Christoph Spötl, Martin Trüssel This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4801555/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 10 Nov, 2024 Read the published version in Scientific Reports → Version 1 posted 12 You are reading this latest preprint version Abstract Cave air temperatures in four caves in the European Alps show statistically significant warming trends of about 0.2 °C per decade over the last two decades (2000-2020). These trends are about half as large as those observed outside and are characterized by a remarkable spatial and temporal consistency. The investigated caves represent different types in terms of their ventilation regime and one of them also hosts perennial ice. Key observation sites are located in cave sections where the temporal variability of air temperature is strongly attenuated compared to outside conditions and data from different cave sections show that the main results are valid for large parts of the investigated caves. Continued warming will lead to broad changes in Alpine cave environments, including changes in the ventilation regime, karst hydrology and subsurface ecosystems. The observed subsurface warming has a particular strong effect on the long-term preservation of perennial ice present in some of these caves. This is shown for an ice cave in the Austrian Alps, where enhanced melt of ice correlates with the observed warming. This cave (and similar ones) will not be able to hold perennial ice beyond the next decade. Earth and environmental sciences/Climate sciences Earth and environmental sciences/Environmental sciences Caves temperature trends Europe regional climate change Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Introduction Atmospheric warming is a global phenomenon of the Anthropocene. Global mean temperatures and temperatures averaged over land areas have warmed by 1.1°C and 1.6°C from 1850–1900 to 2011–2020, respectively [1] . Mountain regions warmed at an even higher rate, which is due to reduced snow and ice feedbacks, reduced local cooling due to shrinking glacier extent, and local land use changes [2] . The European Alps (hereafter Alps for short) have warmed twice as much as the global average between the late 19th and the early twenty-first century [3, 4] . In recent decades, the average warming rate in the Alps has been 0.3°C per decade, outpacing the global rate of 0.2 ° C per decade [5] . The environmental impact of this strong atmospheric warming of the Alps is most prominently manifested by the loss of glacier extent and volume as a result of the rise of the equilibrium line altitude [6, 7] . Parallel to this reduction in surface ice masses, high-Alpine regions have also witnessed a degradation of permafrost [8,9] . Rising temperatures have also affected surface waters in the Alps [10] . A recent study of two catchments in the Tyrolian Alps identified mean warming trends of local streams of 0.24 and 0.45°C per decade for the time period 1977 to 2020 [11] . This signal has also propagated into the subsurface and led to an increase in groundwater temperatures. Shallow groundwaters in Austria (average depth 7 ± 4 m) increased on average by 0.37°C per decade between 1994 and 2013 [12,13] . A study in the northern Alpine foreland (Bavaria, Germany) reported a decrease of the warming amplitude with depth, from 0.28°C per decade in 20 m to 0.09°C per decade in 60 m depth between the early 1990s and 2019 [14] . Rising temperatures have also been detected inside High-Alpine glaciers, e.g., in the Mont Blanc area, revealing a mean warming rate of 0.14°C per decade between 1900 and 2004 [15] . Perrier et al. [16] provided evidence that long-term changes in surface climate have also propagated into urban underground. Surprisingly little is known, however, if and to what extent this atmospheric warming signal has already “arrived” in the interior of cave systems. Caves are abundant in many part of the Alps where karstifiable rocks (chiefly limestone, dolomite and/or marble) are present. The Alps not only host some of the largest caves in Europe (e.g., Hölloch, Siebenhengste-Hohgant Höhlensystem, Schönberg-Höhlensystem), but also caves with ice deposits many of which showing strong signs of retreat, e.g., Dachstein Mammuthöhle [17] , Hundsalm Eis- und Tropfsteinhöhle [18] and Obstanser Eishöhle [19] . Warming of underground cavities has been shown to have also significant consequences for subterranean biota and ecosystems [20, 21] . Here, we document for the first time decadal-scale warming in the interior of caves across the Alps. Although local time series of temperature measurements in different caves (and parts thereof) exist since many years [22, 23, 24] , this is the first study that systematically evaluates temperature measurements obtained in four caves in Switzerland and Austria. The aim of this study is to document and quantify the ongoing warming in the interior of these Alpine caves and to encourage follow-up research, in order to take full advantage of the wealth of temperature data available for other caves in this mountain belt. The microclimate of caves Although hidden in the underground, caves are part of the Earth’s surface and mass and energy are exchanged between them in a variety of ways. The most important processes are advection of air and water. Diurnal and seasonal temperature variations that characterize the Earth surface are attenuated in the cave environment. Ignoring latent heat and humidity effects, the attenuation of the diurnal and seasonal temperature amplitudes follows an exponential function from the entrance to the interior of caves [25] . In the near-entrance part of caves – the heterothermic zone – diurnal and seasonal temperature cycles are well expressed. Deeper into the cave, the attenuation approaches an asymptotic value. This is referred to as the homothermic zone, where diurnal cycles are below detection and also the seasonal temperature amplitude typically does not exceed a few tenths of a degree. This stable temperature is commonly close to the mean annual surface temperature, except where local effects play a role. A compilation of data from 40 caves in Austria located between 403 and 2347 m a.s.l. revealed that most caves show interior cave air temperatures close to the multi-annual mean of the outside temperature [26] . Anomalously low temperatures were measured in caves that are characterized by restricted ventilation and act as cold traps in winter. This study also reported preliminary results of long-term measurements in one Alpine cave (Spannagelhöhle), where temperatures had increased since the year 2000 in a near-linear manner. Preliminary data of a similar monitoring work in caves of the Swiss Alps documented rising temperatures since the year 1990 [27] . In contrast to these well-ventilated caves, a study in a side passage of Postojna Cave (Slovenia) with a dead-end gallery with limited air advection (in contrast to the rest of the system) concluded that there the external warming signal was largely transferred into the cave by thermal conduction, resulting in a delay of 20–25 years [28] . An extreme endmember is a monitoring site in the interior of Bärenschacht, a 60 km-long cave network spanning 946 m of vertical extent. The measurement site is 700 m below the surface with limited air exchange and showed a constant air temperature within 0.1°C between 1999 and 2017 (W. Janz, pers. comm.). Study sites Four caves in different climatological regions across the Alps were selected showing a range in size and elevation (Table 1 ). Rasslsystem (RS), Spannagelhöhle (SP) and Hundsalm Eis- und Tropfsteinhöhle (HA) are located in different parts of the Austrian Alps, while Schrattenhöhle (SH) is located in the Swiss Alps (Fig. 1 ). Table 1 Key parameters of the studied caves. Spannagelhöhle Schrattenhöhle Hundsalm Eis- und Tropfsteinhöhle Rasslsystem Abbreviation SP SH HA RS Cave register number 2515/1 OW 4/61 (22363) 1266/1 3925/9 Coordinates (lat./long.) 47.08°/11.67° 46.77°/8.27° 47.55°/12.03° 46.51°/14.54° Length of cave (m) 12000 19718 264 1085 Vertical extension (m) 329 573 51 161 Classification of cave system dynamic dynamic static with ice dynamic Number of entrances 3 8 1 1 Elevation of (upper) entrance (m a.s.l.) 2523 2090 1495 1085 Name of principle cave monitoring site Porzellanladen Windkluft Jugendgang Perlenhalle Abbreviation SP-PL SH-WK HA-JG RS-PH Elevation of monitoring site (m a.s.l.) 2408 1743 1480 1070 Rock overburden of monitoring site (m) 17 60 50 95 Distance of monitoring site to nearest entrance 840 250 85 270 Relevant outside weather stations JE, PK, PL PI, EN, BO JE, PK BA, SE, ZE Spannagelhöhle is a 12 km-long system in the western Zillertal Alps (Tyrol; Fig. 1 ) with three upper entrances and unknown lower entrance(s) extending over 329 m in the vertical dimension (Supplementary Fig. S1 ). The main cave parts are well ventilated and characterized by upward directed air flow in winter and a reversed flow in summer. The main monitoring site, Porzellanladen (SP-PL), is located in the interior of the cave very close to the main conduit where the air is flowing. Additionally, a site in a dead-end side gallery (Tropfsteingang, SP-TG) and one located in the main air flow but closer to the main upper entrance (Gaunhalle, SP-GH) were included in this study (Supplementary Fig. S1 ). We also present data from a small cave that opens 27 m beneath the main entrance of Spannagel Cave and contained an ice body when discovered in 1998. This so-called Spannagel-Eishöhle (SP-EH) shows a descending geometry and has lost its ice by the year 2002. Hundsalm Eis- und Tropfsteinhöhle is located in the Northern Calcareous Alps of Tyrol (Fig. 1 ) and classifies as a sag-type cave given its subvertical geometry and the lack of a lower entrance (Supplementary Fig. S2). HA is also exceptional because the cave consists of two compartments, a higher one that also contains perennial ice, and a lower section, which is ice-free. The two cave parts were artificially connected in 1996 and then subsequently thermally re-separated using an airlock. Our study concerns primarily the lower section (location Jugendgang, HA-JG) whose air temperature is close to the mean air temperature outside (in contrast to the cold trap in the upper cave part). Temperature data obtained in the main shaft and in a short borehole in the upper part of the cave are also considered in this study (Supplementary Fig. S2). Rasslsystem is one of several caves in the Obir Massif of the Northern Karawanks in SE Austria (Fig. 1 ). Its air flow regime was modified as a result of mining activities which came to an end at the beginning of the 20th century. An abandoned mining adit acts as the lower entrance of cave system (Supplementary Fig. S3) with strong inward and outward air flow during the cold and warm season, respectively, while the pathways of the air in the upper inaccessible parts of the cave are largely unknown [22] . The monitoring site is Perlenhalle (RS-PH), a high and elongated terminal chamber in the innermost part of Rasslsystem with no measurable air flow. Schrattenhöhle is located in the Melchsee-Frutt region of Switzerland (Obwalden; Fig. 1 ) and comprises an extensive network of mostly well-ventilated galleries extending over 573 m and showing multiple entrances (Supplementary Fig. S4). The monitoring site, Windkluft (SH-WK), is located in the main air flow route 250 m from the next entrance, while the subsidiary site SH-WG (Wermutgang) is located close to the uppermost entrance of the Schrattenhöhle cave system. RESULTS We examined daily averages of data measured in the four caves and in the free atmosphere of nearby meteorological stations. The data were subjected to quality control and restricted to the period providing the most consistent records (2000–2020). The main analyses concerned climatological features and relationships with external conditions, trends using different methods including seasonal aspects, spatial and temporal representativeness and the relationship with cave ice developments. Climatological characteristics The analyzed records are characterized by a) different mean values, b) regular annual cycles of comparable magnitude, and c) continuously increasing temperatures at all observation sites (Fig. 2 ). The average conditions in the caves as well as in their outside atmosphere are summarized in Table 2 , which refers to the central range of the overall investigation period (2008–2010). Table 2 Mean and standard deviation of air temperature inside and outside of the investigated caves and respective differences (2008–2010). Bracketed numbers in the lowest panel denote Pearson correlation coefficients between air temperature in the caves and outside. SP-PL SH-WK HA-JG RS-PH Elevation of cave sites [m a.s.l.] 2408 1743 1480 1070 Cave temperature [°C] mean annual winter (open period) summer (closed period) + 1.82 ± 0.066 + 1.77 ± 0.044 + 1.86 ± 0.046 + 3.17 ± 0.064 + 3.12 ± 0.046 + 3.21 ± 0.044 + 4.15 ± 0.054 + 4.13 ± 0.055 + 4.17 ± 0.043 + 5.47 ± 0.022 + 5.46 ± 0.023 + 5.47 ± 0.021 Outside temperature [°C] * mean annual winter (open period) summer (closed period) + 0.40 ± 7.55 1 -6.07 ± 5.18 + 5.66 ± 5.09 + 3.75 ± 7.05 2 -1.42 ± 4.96 + 8.79 ± 4.75 + 4.18 ± 7.60 3 -1.65 ± 4.99 + 9.93 ± 4.83 + 4.94 ± 8.21 4 -1.39 ± 5.37 + 11.2 ± 5.20 Temperature difference in-out [°C] (correlation coefficient) mean annual winter (open period) summer (closed period) +1.42 (+ 0.72) +6.96 (+ 0.50) -4.04 (+ 0.31) -0.58 (+ 0.66) + 4.54 (+ 0.35) -5.58 (+ 0.25) -0.03 (+ 0.46) ´+5.78 (+ 0.44) -5.76 (-0.03) + 0.53 (-0.11) + 6.85 (-0.06) -5.73 (-0.23) * at the elevation of measurements in the cave 1 calculated from AWS data JE and PA, bias from comparison with PA 2 calculated from AWS data EN and PI, bias from comparison with BO 3 calculated from AWS data JE and PA, bias from comparison with HA-entrance 4 calculated from AWS data BA and SE, bias from comparison with ZE The average cave temperatures correspond to the elevations of the cave sites. Lowest values are observed at SP-PL (+ 1.8°C) being located at an elevation of 2408 m a.s.l. SH-WK (1743 m a.s.l.) is characterized by an average temperature of + 3.2°C, HA-JG (1480 m a.s.l.) by + 4.2°C and RS-PH by + 5.5°C (1070 m a.s.l.). The mean temperatures are positive at all cave sites, as well as those in the atmosphere outside of the caves, which range from + 0.4°C at the elevation of SP-PL to + 4.9°C (for RS-PH), respectively. The variability of cave temperatures is almost two orders of magnitude smaller compared to that outside ( < ± 0.1°C vs. ca. ±8°C) and increases with elevation (from ± 0.02 to ± 0.07°C). Thus, RS-PH stands out by an almost homothermic regime, i.e., quasi-constant temperatures despite of the somewhat larger outside variability. Figure 3 graphically summarizes these basic characteristics in terms of boxplots. Table 2 includes information about the average conditions during winter and summer, which in the context of cave microclimate are also referred to as the so-called open and closed periods [18] . In this respect it is remarkable that the average cave temperatures are positive for both periods. Moreover, the cave sites experienced strongly attenuated seasonal temperature differences. Winters are just about 0.1°C colder than summer while in the free atmosphere at the same elevation where winters are about 10°C colder than summer and have negative mean temperatures. During winter the cave temperatures are about 5°C higher than outside, while caves remain colder than outside during summer. Table 2 also shows that the average cave temperatures observed in HA and SH are slightly lower than the outside temperatures at the same elevation, while SP and RS are warmer than outside. This is corroborated by consideration of the observed AWS temperatures (Fig. 4 ), although the latter represent a less reliable measure of local outside conditions due to being located at different elevations. The observed cave temperatures follow a quasi-vertical gradient of about − 0.4°C 100 m − 1 (red dots in Fig. 4 ), which is somewhat smaller than the lapse rate of outside temperatures in the relevant elevation range, i.e., 1000–2500 m a.s.l. and the moist adiabatic lapse rate (blue triangles and dotted line), respectively. However, it should be noted that the cave sites are located in different climatological regimes characterizing the Swiss Alps (SH), the Northern Calcareous Alps (HA), the Central Alps (SP) and the southeastern part of the Alps (RS). Judged by AWS data, temperature gradients are different along valleys in the vicinity of the investigated caves (-0.57°C 100 m − 1 in the Central Alps compared to -0.43 and − 0.41°C 100 m − 1 in the western and eastern and western fringes of the Alps, respectively). Figure 4 also suggests that cave temperatures are correlated with the outside temperatures at the same elevation, which can be considered as an indicator of a physically based response of cave temperatures to changes of the outside atmosphere. This is supported by Table 2 (lowest panel) which also shows that the correlation of daily cave temperatures with concurrent outside values is stronger at the high-elevation sites and during winter, RS being outstanding in this respect. These results are related to cave-specific ventilation patterns which themselves strongly depend on their morphological features (concerning entrances essentially) as will be discussed below. Trend analysis Long-term temperature trends were investigated by calculation of linear trends based on ordinary least squares regression and the Theil–Sen trend estimator as well as using Mann-Kendall trend tests. The methods are applied to quality-controlled cave temperatures covering two decades as well as to AWS data representing temperatures at elevations below and above the cave sites. Visual inspection of the observed records suggests that cave temperatures have systematically increased at all observation sites (Fig. 2 ). This is confirmed by quantitative analyses for the period 2000–2020 (Table 3 ), which reveals positive trends for all investigated sites, i.e., in the caves themselves as well as in the outside atmosphere. Cave temperatures increased by about + 0.2°C per decade, with RS and SP showing higher values than SH or HA. Table 3 Decadal temperature trends for the period 2000–2020 (linear and Theil-Sen, the latter in bold). Outside temperatures refer to the elevation of the cave sites. MK + qualifies the calculated trends to be significant according to the Mann-Kendall test. Spannagelhöhle Schrattenhöhle Hundsalm Eis- und Tropfsteinhöhle Rasslsystem Cave [°C per decade] total period open periods closed periods SP-PL + 0.24 l + 0.26 (MK+) + 0.22 | +0.25 (MK+) + 0.25 | +0.29 (MK-) SH-WK + 0.14 l + 0.15 (MK+) + 0.14 | +0.15 (MK+) + 0.13 | +0.14 (MK+) HA-JG + 0.20 l + 0.12 (MK+) + 0.21 | +0.21 (MK+) + 0.20 | +0.19 (MK+) RS-PH + 0.24 l + 0.28 (MK+) + 0.25 |+0.28 (MK+) + 0.24 |+0.26 (MK-) Outside [°C per decade] total period open periods closed periods at elevation of SP-PL + 0.54 l + 0.53 (MK+) + 0.61 | +0.59 (MK+) + 0.22 | +0.23 (MK-) at elevation of SH-WK + 0.66 l + 0.72 (MK+) + 0.63 | +0.94 (MK+) + 0.48 | +0.53 (MK+) at elevation of HA-JG + 0.59 l + 0.57 (MK+) + 0.62 | +0.62 (MK+) + 0.29 | +0.28 (MK+) at elevation of RS-PH + 0.58 l + 0.53 (MK+) + 0.76 | +0.74 (MK+) + 0.10 | +0.03 (MK+) AWS low [°C per decade] total period open periods closed periods JE + 0.68 l + 0.63 (MK+) + 0.69 | +0.68 (MK+) + 0.36 | +0.32 (MK+) EN + 0.65 l + 0.61 (MK+) + 0.68 | +0.65 (MK+) + 0.36 | +0.32 (MK+) JE + 0.68 l + 0.63 (MK+) + 0.69 | +0.68 (MK+) + 0.36 | +0.32 (MK+) BA + 0.63 l + 0.59 (MK+) + 0.64 | +0.68 (MK+) + 0.30 | +0.27 (MK+) AWS high [°C per decade] total period open periods closed periods PA + 0.54 l + 0.53 (MK+) + 0.61 | +0.60 (MK+) + 0.22 | +0.23 (MK-) PI + 0.70 l + 0.77 (MK+) + 0.68 | +0.76 (MK+) + 0.36 | +0.31 (MK+) PA + 0.55 l + 0.53 (MK+) + 0.61 | +0.60 (MK+) + 0.22 | +0.23 (MK-) SE + 1.11 l + 1.17 (MK+) + 0.96 | +0.99 (MK+) + 0.88 | +0.93 (MK+) The outside temperatures at the elevation of the cave sites increased at higher rates of about 0.5°C per decade. This value is independent of the applied calculation method and there is no clear evidence regarding elevation dependencies or regional patterns. It is notable, however, that temperatures at low-lying weather stations are characterized by rather consistent trends in the order of + 0.6°C per decade, while some distinction can be observed regarding the investigated high-elevation sites. Hence there is an indication that the eastern and western fringes of the Alps (i.e., the RS and SH regions) experienced enhanced warming compared to the Central (SP) and the Northern Calcareous Alps (HA), which however, is not reflected by the cave data. We also examined seasonal trends representing open (Nov-May) and closed periods (Jun-Oct), which provide evidence that the outside temperatures increased stronger during winter compared to summer (Table 3 , open vs. closed periods). Naturally, this also concerns the underlying AWS data and contrasts the cave conditions, where significant differences between winter and summer temperature trends can hardly be discerned. Table 3 also shows the results of the Mann-Kendall tests to explore the temporal consistency of the observed cave temperature trends. There is evidence that all investigated cave records are characterized by statistically significant trends, which also holds true for the investigated sub-periods (open and closed periods) and the calculated outside temperatures. The significance of positive trends can also be stated for the weather stations being used for the calculation of the outside temperatures, with the exception of summer temperatures at Patscherkofel. These results are based on observations representing the two decades providing the most consistent data (2000–2020). However, the underlying time series have different length and episodic gaps (Fig. 2 ), which in principle can affect the robustness of the so far presented analysis. Moreover, the spatial representativeness of the presented cave records may be influenced by local effects (e.g., ventilation conditions, episodic water influx or icing) and diverse observation problems (such as unrepresentative sites, different instruments and/or failure due to e.g. icing). A rigorous analysis of respective effects on our trend analysis is hampered by the scarceness of appropriate cave records over longer time spans. To at least qualitatively strengthen the respective validity of the presented results we considered some additional temperature records from other positions in these caves. Corresponding results are addressed in the following sections and we here just mention strong evidence that all investigated cave sites experienced a warming since at least 1990 and beyond 2020 until the recent days. In this context, we also refer to temperature data across the so-called Greater Alpine Region (GAR), which represents a gridded data set across the Alps derived from long-term homogenized station data [3] which allows putting the presented results in a broader temporal and spatial context. We here consider a limited ensemble of these data, i.e., back to 1970 as well as low/high-elevation, West and Southeast regions, respectively. Figure 5 demonstrates that these records are characterized by consistent and fairly linear trends in the order of + 0.5°C per decade, with indication of enhanced warming in the lowlands compared to high-lying regions (+ 0.56 vs. +0.48°C per decade). The analysis further shows that the observed cave records as well as the calculated outside temperatures fit these regional developments, which essentially confirms their spatial and temporal representativeness. Note that the Austrian caves are located in HISTALP-regions West and Southeast. It may be added that Jungfraujoch, i.e., a Swiss mountain station some 50 km south of SH, is included in the GAR-High-Level data, where a warming rate of + 0.45°C per decade was observed during the period 1970–2020. Case study SP SP has several entrances at different elevations, which has implications for the ventilation regime of the cave. Figure 6 (red line) again shows the previously presented SP-PL data and outside temperatures which are both characterized by consistent warming (Table 3 ). Addressing the question of respective spatial representativeness, we explored whether trends can also be detected in two other parts of this cave. Site SP-GH is located closer to the (upper) entrance while site SP-TG is located deeper into the cave, in a dead-end side passage (Supplementary Fig. S1 ). SP-GH largely mirrors the SP-PL data concerning the overall level (+ 1.68 vs. +1.82°C), the seasonal variability and the trend (+ 0.36 vs. +0.24°C per decade). The slightly higher temperatures at SP-PL can be attributed to the larger distance from the main (upper) entrance (840 m compared to 470 m for SP-GH). SP-TG is characterized by a different thermal regime with higher temperatures (+ 2.16°C on average) and a strikingly reduced seasonal variability. This can be related to its location in kind of chamber deeper down in the cave and off the main air flow route. It is remarkable that also this remote cave section experienced a warming (+ 0.27°C per decade) which according to the MK-test is statistically significant. We further studied temperature records from a small separate cave close to the entrance of SP (Spannagel-Eishöhle, SP-EH) which contained perennial ice until the year 2002. Figure 6 (blue line) shows that also this cave experienced a warming during the last two decades. The overall lower air temperature as well as the tenuous trend at the beginning of the record may have been related to the presence of ice. The seasonal amplitude of air temperature is comparatively large due the small size and ventilation of this cave. Summarizing, there is strong evidence that Spannagel cave as a whole and nearby caves were affected by a consistent warming over the last two decades at least. Case study HA HA is special in the context of the four studied Alpine caves in that its ventilation regime is much more restricted given its sag-type geometry and the presence of perennial ice (Supplementary, Fig. S1 ). This leads to a cold-trap setting in the upper part of the cave, which will be examined in more detail in this chapter thereby also considering wall-rock temperature data. First, we again note the progressive damping of outside temperature changes as they progress into the cave, which is a key feature of cave microclimate and strongly depends on ventilation characteristics. In sag-type caves (as HA) the damping is most pronounced during summer when thermal stratification prohibits effective exchange of air between the cave interior and the outside. Conversely, cold air can relatively easily penetrate into the deeper parts of the cave. We observed a seasonal temperature variation of about 30°C outside, while at a depth of 29 m below the main entrance (site HA-29) the air temperature varies by just about 5°C. We note that the conditions in this cave section are not only preconditioned by effective cold air pooling but also by the existence of perennial ice (Supplementary Fig. S1 ). The lower part of the cave is thermally separated from the upper one (site HA-JG), contains no ice and is characterized by a homothermic regime (air temperature hardly varies on an annual time scale) at comparatively high temperatures (average + 4.3 compared to + 0.1°C at site HA-29). Figure 7 also shows long-term observations of rock temperatures at a distance of 126 cm behind the cave wall next to site HA-29. The average rock temperature is + 0.4°C, thus being slightly higher compared to the cave air at the same depth (+ 0.1°C). During winter the rock does not cool as much as the air, while air and rock are thermally adjusted during summer. Hence there is a systematic difference in seasonal variability. Closer inspection also reveals that changes in rock temperature systematically lag behind those in the air [18] . Figure 7 finally supports that each of the considered HA sites experienced a systematic warming throughout the investigation period 2010–2020. This is corroborated by the calculated trends (0.24 / 0.79 / 0.99°C per decade for HA-JG / HA-29 / HA-126) and the according MK-trend tests (all attesting statistical significance). It can thus be concluded that the results, which refer to HA-JG, can serve as a reference for the long-term evolution of air temperature in this cave. Moreover, there is indication that trends are stronger in cave sections containing ice (HA-29) and rock temperatures (HA-126) compared to relatively remote cave sections without ice (HA-JG). Surface height changes of ice in HA (Supplementary Fig. S2) were monitored by annual stake readings close to site HA-29. This effort was motivated by the expectation that long-term warming traces in enhanced melt of ice, as is generally known for snow and glaciers outside [5] . Indeed, Fig. 8 (blue line) shows a continuous loss of ice, which was accelerated during the second half of the investigation period (2008–2018). The values of air temperature averaged over time between the stake readings (black line) show a transition from initially negative values towards mostly positive values. This analysis confirms a strong physical relationship between long-term temperature evolution and the development of ice in caves, as is known from ice caves elsewhere [29, 18, 19, 30] . Discussion Ongoing climate change is leading to a regional and global rise in the temperature of the atmosphere. The evidence is essentially based on a sophisticated synthesis of observations from weather stations, radiosondes and satellites, as well as numerical models [1] . Equally important evidence comes from indirect sources relating to corresponding changes in the ocean (e.g., rising sea-surface temperatures), in the cryosphere (e.g., vanishing glaciers and ice sheets and thawing permafrost) and in the biosphere (e.g., shifting vegetation belts). This raises the question of whether climate change has also “arrived” in the caves, which is investigated in this study on the basis of long-term observations from caves in the Alps. The data from four caves in different parts of this mountain range provide ample evidence that this is indeed the case. Principal data were selected considering their quality as well as their spatial and temporal representativeness. The main analysis is based on data covering two decades (2000–2020) and linear trends were calculated and tested for statistical significance. Extended temporal representativeness of key results is shown back to the 1990s, however, as well as correspondence to developments in other cave sections and regional climate changes outside the caves. Context to previous work The analysis of the studied cave temperature records reveals a consistent picture in terms of linear trends ranging from + 0.1 to + 0.2°C per decade (Table 3 ). These figures are backed by the application of two calculation methods (least squares and Theil-Sen) and the statistical significance is confirmed by the Mann-Kendall test. Additional temperature records from each cave show that the principal records are representative for major parts of the investigated caves and not just for individual sites. These results overall agree with previous findings which are, however, mostly based on shorter time series and less elaborated analyses. Spötl and Pavuza [26] reported temperature trends for Spannagel cave of + 0.5 and + 0.3°C per decade (2000–2015) obtained from sites GH and TG, respectively. These sites are located above and below site SP-PL, which is considered in this work. Although these figures refer to different periods, they also indicate that the sections closer to the upper entrance have experienced greater warming. Trüssel [27] analyzed temperature records from Schrattenhöhle and reported a warming of 0.7°C for Windkluft (SH-WK in this study) for the period 1990–2018, which corresponds to a trend in order of + 0.3°C per decade (decreasing through time). Similar warming rates were observed for a shorter period of time (1990–2010) in cave sections above Windkluft (Fünferkluft and Wermutsgang). The author critically pointed out that the cave sections for climate studies must be carefully selected in order to avoid misinterpretation that can arise, for example, from an inappropriate selection of observation sites. Filipponi [31] proposed criteria suggesting that SH-WK is representative of large parts of the cave. Wind et al. [18] addressed long-term temperature changes in Hundsalm Eis- und Tropfsteinhöhle. Significant warming is apparent in the ice-free part of the cave throughout the period 2008–2020 (+ 0.24°C per decade at site HA-JG, which is a principal site in this study). This value agrees with the findings of this study. Ongoing warming is also observed in the higher ice-bearing section of the cave, where trends depend on the distance to the ice (+ 0.27 for T36 vs. +0.54 for T29; the latter being located closer to the ice). Some related information is available from neighboring countries. Domínguez-Villar et al. [28] observed that temperatures in Postojnia cave (Slovenia) increased at a rate of 0.04°C per decade (2009–2013). The application of a model indicates that the cave was prone to warming since the 1980’s at least and enhanced warming is predicted for the future (+ 0.15°C per decade). Analyses of historical data since 1935 indeed show that in this cave the average temperature has increased up to 2°C [32] , corresponding to a trend of + 0.25°C per decade. Significant warming was also reported for another cave in Slovenia (Predjama Cave) as well as for the average temperature across Slovenia (+ 0.34°C per decade during period 1961 and 2011). The calculated temperature trends for the free atmosphere at the elevation of our studied cave sites are about twice as large as those observed inside the caves (+ 0.5 to + 0.7°C per decade). These figures are in agreement with site-specific observations for Patscherkofel [26] , Bonistock [27] and Hahnenkamm [18] . Since a straightforward comparison is hampered by the different elevations and periods, we also placed our results in an extended temporal and regional context provided by regionally gridded climate data i.e., the HISTALP data set [3, 33, 34] . These data show the generally strong warming of the Alpine region, which is about twice as large as the Northern Hemispheric average and affected all subregions of the Alps. Enhanced and fairly homogeneous warming occurred since the 1980, with rates matching those from our study (ca. +0.5°C per decade). Our analysis indicates that trends in the low-lying regions are relatively uniform, while mountain stations show some regional distinction, i.e., smaller trends in the Central Alps compared to the Swiss Alps. Referring to the vicinity of Monlesi Ice Cave, Luetscher et al. [35] note that the Jura Mountains of northern Switzerland experienced a warming of ca. 0.5°C per decade (1980–2000). Based on a gridded data set across Switzerland, Ceppi et al. [4] found an average warming rate of + 0.35°C per decade (1959–2008). The Obir mountain area in South-Eastern Austria is characterized by outstanding warming rates (ca. +1°C per decade) which, however, is not fully reflected in corresponding HISTALP data. Elevation-dependent and thus also regional warming in the Alps is controversially discussed in literature. Recent reviews conclude that a general picture is difficult to draw due to the variety of influencing factors and processes being strongly related to small-scale topographic effects [36, 2, 34] . It is therefore also difficult to interpret the result that trends are significantly larger during the open periods (winter) compared to the closed periods (ca. +0.6 vs. +0.3°C per decade). A thorough investigation of local characteristics (such as cold air pooling, mesoscale wind systems, snow cover, radiation, cloudiness) that influence outside temperature is beyond the scope of this work. Using data from more than 1000 mountain sites across the globe, Pepin and Lundquist [37] showed that high-elevation temperature trends are highest near the annual 0°C isotherm due to snow-ice feedback, which was also observed in Switzerland [4] . Context for cave-specific conditions and processes A comparison of the average cave temperatures with those outside reveals a small difference and a clear dependence on elevation (Table 2 , Fig. 4 ). Both suggest that the cave temperatures overall correspond to the external conditions, as is corroborated by an independent analysis considering more than 60 caves across the Austrian Alps [26] . This result is also in agreement with theoretical considerations [38, 25] , which show that cave temperatures are characterized by a damped and delayed response to external forcings, both increasing with distance from the entrance(s). Observations from caves across the world support the related hypothesis that temperature in deeper cave zones asymptotically approaches an equilibrium rock temperature. Following Badino [39] , the latter roughly corresponds to the long-term outside temperature. Short-term fluctuations are effectively filtered out due to the large thermal capacity of the rock in which the caves formed. Application of this conceptual model to karstic rock with a thickness of 100 m yields equilibration time scales (i.e., the time to attain thermal equilibrium) in the order of 5–50 decades, depending on the relative importance of energy exchange involving flow of air and water. Since our studied caves have a rock overburden of less than 100 m it can be argued that the analyzed data were obtained at sites where the seasonal variability is strongly damped, while the climatically relevant variability is still captured, albeit with a delay. These theoretical considerations support the observation that warming trends inside caves are smaller that outside (Table 3 ). It also helps to understand the nature of the time series mainly considered in this study, which essentially show long-term changes (trends) superimposed by weak seasonal variability (Fig. 2 ). Systematically delayed responses of cave air temperatures to external forcings cause thermal imbalances, which are reflected in the differences between mean temperatures inside and outside of the caves, as our results also show (Fig. 4 ). The thermal cyclicity observed in cave records has often been used to classify caves and zones therein. A basic distinction concerns the heterothermic and homothermic zones [40] . Considering that the analyzed cave data are characterized by a seasonal variability of the order of a few tenths of a degree (compared to an external variability of about +/-15°C), they cannot represent truly homothermic conditions, but remote sections of heterothermic cave sections. Sedaghatkish et al. [41] used a numerical model to simulate heat transfer in a karst conduit due to air flow and thermal conduction, which suggests a more physically based distinction of a “diffusive region” (where external forcings mainly propagate through thermal conduction in the rock), a “convective region” (where heat is mainly transported by heat transfer in air) and a “deep region” (where temperatures are quasi-constant). This distinction basically matches the observationally based approach by Cropley [42] and suggests that three of the principal sites being considered in this work can be assigned to the second class (SP, SH and RS; i.e., heterothermic), while HA-JG is similar to the third one (which essentially conforms to the homothermal regime). Medina et al. [21] analyzed records from twelve caves across the globe and found three patterns which represent caves sites retracing external variations with little modification except for reduced amplitudes, sites featuring low correlation to surface temperature and a slight thermal delay, and cave sections showing extreme delays or even no relationship to surface variations. Our cave data match the first category, which the authors expect to represent caves with multiple entrances, where air flow plays a dominant role in controlling the cave microclimate (supported by our observations). Seasonally changing airflow patterns are known for three of the investigated caves, i.e. SP, SH and RS [26, 27, 31] , which have two or more entrances at different elevations (Table 1 ). This, as well as the seasonally changing signs of the thermal gradient (inside vs. outside) are strong indications that processes related to advection and/or convection are important for the exchange of matter and energy with the outside. The so-called chimney effect has been recognized as a key process in this context, which is driven by seasonally changing density differences between the subsurface air and the outside air [38, 43] and can be significantly modified by synoptically induced pressure differences [44] . The effect is known to be important in Spannagelhöhle, where it induces a positive deviation of the cave temperature in sections close to the upper entrances compared to outside (Fig. 2 ; Table 2 ) [26] . This is most evident during winter when warm air is forced towards the upper entrance to which sites SP-PL and SP-GH are relatively close. The direction of the air flow switches at a threshold temperature of about 2°C [45] . Trüssel [27] showed that air flow in Schrattenhöhle is also driven by the chimney effect, but it is interesting to see that the temperature at the investigated site (SH-WK) is lower than outside. This can be related to the position of the investigated site closer to lower entrances, which in the context of the chimney effect are prone to be colder than regions closer to the upper entrances. However, air flow in SH is generally more complex due to the many entrances, which favor differential air flow in different branches of the cave. We note that wind-driven pressure differences can additionally complicate cave ventilation [44] , which could play a role in Schrattenhöhle, too. The RS is characterized by relatively high temperatures compared to outside. Also this cave is known for air flow being driven by the chimney effect which, however, is relatively weak due to the remoteness of this site relative to the main pathway of the air (Supplementary Fig. S3). Still, the direction of air flow seasonally reverses at an outside temperature threshold of 5.5°C [22] . Hundsalm Eis- und Tropfsteinhöhle is exceptional in several respects. First, it is the smallest cave and its sag-type geometry effectively traps cold air during winter, which enters the cave by density-driven convection. During summer, the air pools without significant exchange with the surface due to a stable thermal stratification, which favors the existence of perennial ice in the upper cave section [18] . The cave is therefore classified as “static with ice” [46] . These characteristics give rise to an asymmetric seasonal temperature cycle, which is not observed in the other caves. However, the mainly considered data (HA-JG) are much less affected by these processes, since they were collected in a relatively isolated chamber in the deepest section of the cave and therefore represent conditions corresponding to a homothermic regime. Having thus addressed some individual characteristics of the studied caves, it is even more remarkable that the magnitudes of the calculated temperature trends are comparable, i.e., ranging between + 0.14 to + 0.24°C per decade for the principally investigated sites. It should also be remembered that these sites are located in different regions across the Alps and at different elevations. The rather uniform warming rates can be explained by the fact that records with comparable seasonal variability were considered, which in principle can lead to similar long-term external signals. The microclimates of different cave sections, however, retain some spatial variability due to individual physical characteristics such as the relative distance to entrances which influences related processes (e.g., heat conduction, ventilation and flow of water). Our case studies show a correspondingly different seasonality, which is reflected by a progressive attenuation (Figs. 6 and 7 ), which was already discussed. As far as the trends in the different cave sections are concerned, the results indicate that the decadal warming was stronger near the entrance(s) than in more remote sections (SP-GH vs. SP-TG and HA-29 vs. HA-JG). This is also in agreement with theoretical concepts that predict that long-term, i.e. low-frequency, external variability is delayed with increasing distance from the entrances due to the large thermal capacity of the rock [39] . Indirect evidence of warming in the investigated caves Temperature strongly controls the formation and existence of ice, which is present in two of the investigated caves. Hundsalm Eis- und Tropfsteinhöhle hosts perennial firn and ice in its upper section, which reaches a maximum thickness of about 5 m. The ice formed by metamorphism of snow falling though one of the entrances, as well as by refreezing of seepage water [47, 18] . Our results confirm a development towards increasingly negative mass balances during the recent decade (Fig. 8 ), in parallel with the concurrent atmospheric warming in the Alps. Wind et al. [18] developed a model, which allows to calculate annual melt rates as a function of the external temperature. Given the current warming rates, the model predicts that within the next decade, this cave will no longer preserve perennial ice (only seasonal ice). Note, however, that the mass balance of an ice body not only includes surface melt (ablation) but also accumulation and basal melting, both of which also depend on temperature. A small cave next to Spannagel cave (Spannagel Eishöhle, SP-EH) contained perennial ice until about 2002 and has been ice-free since then. This adds evidence that caves are sensitive indicators of regional climate change, similar to snow cover and glaciers outside. Significant reduction of ice in caves is a wide-spread phenomenon across Europe [40, 48, 17, 49, 50, 18] and there is no indication that this trend will not continue (if not accelerate). Potential influences of anthropogenic cave manipulations Certain sections in the investigated caves are open to guided tours, which necessitates some related discussion. Touristic use in Spannagelhöhle began in 1994, involving ca. 500 m of galleries behind the upper entrance (open all year round). HA was opened as a show cave in 1967, and the guided tours cover most of the upper, i.e., ice-bearing part (open from May until October). The lower section (in which site HA-JG is located) was discovered in 1984 and is not accessible to visitors (locked). Rasslsystem (RS) is part of the Obir Cave system. Guided tours are only offered for part of this system that are far away from RS (since 1991; open April to October). Schrattenhöhle is also can also be visited within guided tours, but this s only possible since the year 2000 i.e., after the investigation period and in cave sections far off site SH-WK. Touristic use can have severe impacts on the thermal regime of caves [51] . The main influencing factors are the heat emitted by people and the lighting. Evidence is provided by comparing historical temperature data with recent ones and by comparing cave sections or neighboring caves with and without touristic influence [32] . Microclimatic studies usually show increased daily variability which related to the number of visitors or concurrent changes in, e.g., humidity and CO 2 concentrations [52, 53] . In general, anthropogenic influences are difficult to disentangle from natural changes, as they are cave-specific and less pronounced in well-ventilated caves or sections. Regarding the potential impact on the investigated caves, it should first be noted that in all three caves the main study sites (SP-PL, HA-JG and RS-PH) are not part of guided tours, i.e. the sites are relatively far from the direct impact of visitors. In terms of timing, tourist activities are rather confined to summer and autumn, which means that excess heat in outer cave sections (where guided tours are offered) is likely to be pushed out through the entrances. This may be relevant for Spannagel and Rasslsystem, but not for Hundsalm Eis- und Tropfsteinhöhle. This sag-type cave is special because its microclimatic summer regime is characterized by stable air stratification. In principle, this has meant that the excess energy input from the tourists has slowly warmed the cave, which in turn will likely have been compensated for by the increased melting of snow and ice (Fig. 8 ). Wind et al. [18] noticed small transient disturbances of the cave air temperature being related to guided tours, but these were equilibrated shortly after. Potential long-term effects have not been observed. Regarding the principally considered site in this cave (HA-JG, i.e. the lower, ice-free part), warming by tourists cannot play a role, since it is located in a thermally isolated chamber underneath the ice-bearing part (Supplementary Fig. S1 ). Finally, as far as the decadal trends are concerned, it should be noted that touristic use in all these caves began before the main period under consideration (2000–2020). It can therefore be assumed that the trends analyzed are not distorted by tourist activities. Data from Schrattenhöhle are not at all affected by such effects because guided tours started after the investigation period. Consequences of warming caves The demise of cave ice as a consequence of warming of cave environments was already been addressed in the context of the Hundsalm Eis- und Tropfsteinhöhle and the (former) Spannagel Eishöhle (Figs. 7 and 8 ). Related processes have been elucidated by e.g. Luetscher et al. [40] or Obleitner and Spötl [54] . These studies also showed the importance of refreezing seepage water, whose availability strongly depends on the permeability of the host rock, which at high-elevation Alpine caves may also contains sporadic permafrost. Systematic warming is therefore expected to significantly change the availability (amount and timing) of water in these caves, which in turn will also affect the thermal environment [39] . The degradation of sporadic permafrost can also lead to the destabilization of cave passages, which happened in Hundsalm cave in 2020 (rockfall). Warming can affect the ventilation regime of high-elevation caves by changing the upper entrances (which otherwise may be plugged by snow or ice). Such effects have already been observed, e.g., in Schrattenhöhle [27] . Finally, warming of subterranean environments will potentially have a large impact on biota. For example, Mammola et al. [20] argued that caves are largely unexplored habitats for species (bats, spiders, bacteria) that are likely to be sensitive to even small changes in the microclimate [55] . Conclusions Statistically significant changes of air temperature were observed in four caves across the Alps. Our analysis covers the recent two decades and different cave types (static vs. dynamic and also with ice) as well as different sections therein. Key results, interpretations and conclusions are summarized as follows: Mean air temperatures inside these caves are close to those outside and are characterized by much reduced seasonal variability. Caves are systematically warmer in winter (colder in summer) and quasi-vertical temperature gradients in caves are smaller than outside. Trends of cave temperature increase are about half of that outside (0.2 vs. 0.5 °C per decade). These warming rates are characterized by a remarkable linearity and a pronounced spatial and temporal consistency compared to external conditions. These figures represent cave sections characterized by an almost homothermic regime, although there are indications of increased warming rates near the entrances. This and other observations (reduced and lagged seasonal variability) are in line with basic theoretical expectations. The consistency of these trends is remarkable given the variable geometry of the investigated caves, and the results are considered representative of cave environments in the Alps. The synthesis of the data confirms that global warming has reached the interior of these cave systems. It is not expected that the touristic use of some of these Alpine caves has significantly influenced the results. Continued warming may cause feedbacks and broader changes in cave environments, e.g. changes in the ventilation regime and increased ice melt, which may also affect water resources and subsurface ecosystems. This study shall stimulate sustained long-term monitoring of temperature and other physical or chemical parameters, as well as related modelling, to better understand the current and future responses of caves in a warming world. METHODS Cave air temperatures were logged using different instruments. Most data from the Austrian caves were obtained using HOBO TempPro v2 logger (Onset Computer Corp.; resolution 0.02 °C, accuracy ±0.21 °C); the first few years of some of these time series were recorded by an earlier model, HOBO Optic StowAway Temp (Onset Computer Corp.; resolution 0.16 °C, accuracy ±0.2 °C). In one cave, Hundsalm Eis- und Tropfsteinhöhle, we also logged the rock temperature in a borehole drilled horizontally into the cave wall (diameter 2.5 cm) using a HOBO TempPro v2 logger with an external sensor (Onset Computer Corp.; resolution 0.02 °C, accuracy ±0.21 °C). The calibration of both these logger types was checked using an ice-water bath in a dewar prior to the installation in the cave and after the loggers were brought back to the laboratory. No drift was observed in any of the time series. The final data were corrected for the offset at 0°C, which for most HOBO TempPro v2 logger was less than 0.1 °C. The logging interval in the Austrian caves varied between 1 and 3 hours. Temperature data in the Schrattenhöhle were collected using a different logger model, i.e., Hamster and Hamster-A (Elpro; resolution 0.05 °C, accuracy ±0.2 °C) and the calibration provided by the manufacturer was used. Occasional cross-checks using other logger types suggest an overall accuracy of ±0.1 °C. The logging interval in Schrattenhöhle was 1.5 hours (1990-1999) and 1 hour (since 2000). These instruments proved reliable to recover cave temperatures in a methodically consistent manner and over long periods of time. Minor gaps in the records were mostly due maintenance (downloading the data or change of sensors) which in caves can need longer time. Major gaps evolved in the context of post-processing raw data focused on detection of physically unplausible data and comparison to parallel series as far as available. It is helpful in this context, that cave temperatures are less variable compared to outside temperatures. The reason(s) for physically unplausible data can be related to various environmental issues including, e.g. intermittent influences by water, rime or ice on sensors. About 15% of the original data from Hundsalm and Spannagel caves were finally excluded from further analysis, significantly more from SH and RS (62 and 42%). The postprocessed data were then merged to daily averages, which form the basis of this study. The mostly complex geometrical structure of caves and related impacts on their microclimate questions the spatial representativeness of analysis referring to a certain cave section. Practically, the siting of a sensor is often determined by on-site compromises considering scientific arguments, accessibility and constraints to safely mount sensors, which usually can hardly be visited several times a year on the other hand. Related uncertainties are difficult to quantify and would require an extensive measurement network or detailed numerical simulations of air flow along the investigated caves (both being beyond the scope of this work). Nevertheless, in order to address the topic of this investigation (temperature trends) we also analyzed observations in neighboring sections of the studied caves to study the spatial representativity of the best possible chosen reference time series. Daily values of air temperature measured at automatic weather stations (AWS) in the vicinity of the four caves were used to calculate free-atmosphere temperatures at the elevation of the cave sites. The latter were derived by linear interpolation based on daily lapse rates calculated from station pairs located at the valley bottom and at crest height in representative valleys. Correspondingly, we used the station pair JE (530 m a.s.l.) and PA (2251 m a.s.l.) to derive the lapse rate for the SP and HA environment, while BA (623 m a.s.l.) and SE (1040 m a.s.l.) are used for RS and EN (1000 m a.s.l.) and PI (2100 m a.s.l.) for SR, respectively. The calculated temperatures were validated using additional AWS data (Plattkopf (2260 m a.s.l.), HA-entrance (1520 m a.s.l.), Bonistock (2160 m a.s.l.) and Zell-Pfarre (900 m a.s.l.; Fig. 1) and corresponding offsets were applied to refine the calculated temperatures at the elevation of the cave sites. Finally, the ice development in HA cave was monitored using stakes drilled into the ice yielding annual measures of surface height changes with an estimated accuracy of ±1 cm [18] . Characterization of the typical meteorological conditions in the studied caves and their outside environment is based on standard statistics for a reference period (2008-2010) with most consistent data. Shorter and longer periods are considered in the context of case studies addressing specific problems, e.g. the spatial representativity. Trend analysis employs the calculation of linear trends and Mann-Kendall tests [56, 57] . The former is firstly based on ordinary least squares regression of daily values covering the period 2000-2020 and consideration of respective slopes. Since this method can suffer from its sensitivity to outliers, we also employed the Theil–Sen trend estimator [58, 59] . This is a comparatively robust method to fit a line to sample points by considering the median of slopes determined by all pairs of successive data. Calculations were performed using Matlab® [60] and were also applied to seasonal subsets covering the so-called open and closed periods. Following Wind et al. [18] , the latter comprise the months November to April and May until October, respectively. Declarations Acknowledgements Long-term monitoring in the Austrian caves was partially supported by the Austrian Science Fund (FWF) grants Y122-GEO, P18637-N10, P20618-N10 and P31874-N37 to CS. Data from automatic weather stations were provided by national sevices i.e., Geosphere Austria, Slovenian Environment Agency, Hydrographischer Dienst Tirol, Hydrographischer Dienst Kärnten and VERBUND AG, Austria. For open access purposes, the authors have applied a CC BY public copyright license to any author accepted manuscript version arising from this submission. Funding for research in the Swiss caves was partially provided by the Foundation Naturerbe Karst und Höhlen Obwalden (NeKO), Alpnach, Switzerland. Several members of caving associations provided access to caves and on-site support. Author contributions The idea of the manuscript was developed by CS and MT, who also provided the cave data. Main analysis was provided by FO, who also drafted the manuscript, tables and figures except of Fig. 1. Figures in the supplement were contributed by CS, where Fig. 4 is adapted from Filipponi (2003). The manuscript was reviewed by all authors, thereby also considering comments by Beat Niederberger. The authors declare no competing interests and underlying data can be provided on request. Funding Long-term monitoring in the Austrian caves was partially supported by the Austrian Science Fund (FWF) grants Y122-GEO, P18637-N10, P20618-N10 and P31874-N37 to CS. Competing interests The authors declare no competing interests. Data availability Data being used in this study are available upon request. Please contact [email protected] or [email protected] . References Gulev, S.K. et al.(Changing State of the Climate System. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (ed. by Masson-Delmotte, V. et al.). Cambridge University Press, Cambridge, 287-422. https://doi.org/10.1017/9781009157896.004 (2021). Mountain Research Initiative EDW Working Group. Elevation-dependent warming in mountain regions of the world. Nature Climate Change 5 , 424-430 (2015). Auer I., R. et al.. HISTALP – historical instrumental climatological surface time series of the Greater Alpine Region. Int. J. Climatol. 27 , 17-46 (2007). Ceppi, P. et al.. Revisiting Swiss temperature trends 1959-2008. Int. J. Climatol. 32 , 203-213. https://doi.org/10.1002/joc.2260 (2012). Hock, R.G.et al. . High Mountain Areas. In: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (ed. by H.-O. Pörtner, H.-O. et al.). Cambridge University Press, Cambridge, 131-202. https://doi.org/10.1017/9781009157964.004 (2019). Sommer, C. et al.. Rapid glacier retreat and downwasting throughout the European Alps in the early 21st century. Nature Commun. 11 :3209. https://doi.org/10.1038/s41467-020-16818-0 (2020). Žebre, M. et al.. 200 years of equilibrium‑line altitude variability across the European Alps (1901−2100). Climate Dyn. 56 , 1183-1201 (2021). Beniston, M. et al.(2018). The European mountain cryosphere: a review of its current state, trends, and future challenges. The Cryosphere 12 , 759-794 (2018). Haberkorn, A., Kenner, R., Noetzli, J. & Phillips, M. Changes in ground temperature and dynamics in mountain permafrost in the Swiss Alps. Frontiers in Earth Science 9 :626686, https://doi.org/10.3389/feart.2021.626686 (2021). Webb, B.W. & Nobilis, F. Long term water temperature trends in Austrian rivers. Hydrol. Sci. J. 40 , 83-96 (1995). Niedrist, G.H. Substantial warming of Central European mountain rivers under climate change. Regional Environ. Change 23 :43 https://doi.org/10.1007/s10113-023-02037-y (2023). Kormann, C., Francke, T. & Bronstert, A. Detection of regional climate change effects on alpine hydrology by daily resolution trend analysis in Tyrol, Austria. J. Water Climate Change 6 (1): 124-143. https://doi.org/10.2166/wcc.2014.099 (2015). Benz, S.A., Bayer, P., Winkler, G. & Blum, P. Recent trends of groundwater temperatures in Austria. Hydrol. Earth Syst. Sci. 22 , 3143-3154 (2018). Hemmerle, H. & Bayer, P. Climate change yields groundwater warming in Bavaria, Germany. Frontiers in Earth Sci. 8 : 575894. doi: 10.3389/feart.2020.575894 (2020). Gilbert, A. & Vincent, C. Atmospheric temperature changes over the 20th century at very high elevations in the European Alps from englacial temperatures. Geophys. Res. Lett. 40 , 2102. doi:10.1002/grl.50401(2023). Perrier, F., Le Mouël, J. L., Poirier, J. & Shnirman, M.G. Long-term climate change and surface versus underground temperature measurements in Paris. Int. J. Climatol. 25, 1619-1631. https://doi.org/10.1002/joc.1211(2005). Spötl, C., Wimmer, M., Pavuza, R. & Plan, L. Ice caves in Austria. in: Ice Caves (Perşoiu, A. & Lauritzen, S.E.), 237-262, Amsterdam, Elsevier (2018). Wind, M., Obleitner, F., Racine, T. & Spötl, C. Multi-annual temperature evolution and implications for cave ice development in a sag-type ice cave in the Austrian Alps. The Cryosphere 16 , 3163-3179. https://doi.org/10.5194/tc-16-3163-2022 (2022). Spötl, C., Egger, M., Mangini, A. & Pavuza, R. Das Karstsystem am Obstanser See (Kartitsch, Osttirol) – geowissenschaftliche Ergebnisse. Die Höhle 68 , 59-78 (2017). Mammola, S. et al.. Climate change going deep: The effects of global climatic alterations on cave ecosystems, Anthropocene Rev. 6 , 98-116 (2019). Medina, M.J. et al.. Temperature variation in caves and its significance for subterranean ecosystems. Sci. Rep. 13 : 20735. https://doi.org/10.1038/s41598-023-48014-7 (2023). Spötl, C., Fairchild, I.J. & Tooth, A.F. Cave air control on dripwater geochemistry, Obir Caves (Austria): Implications for speleothem deposition in dynamically ventilated caves. Geochim. Cosmochim. Acta 69 , 2451-2468 (2005). Borreguero, M. & Pahud, A. Climatologie de la grotte des Pingouins. In: Lapi di Bou: recherches et explorations spéléologiques,1974-2009, 41-54 (2009). Wimmer, M. 25 Jahre Eis- und Höhlenklimaforschung im Schönberg-Höhlensystem. Mitt. Landesvereins. Höhlenk. Oberösterreich 65/124 , 17-25 (2019). Wigley, T.M.L. & Brown, M.C. The physics of caves. In: The Science of Speleology (ed. by Ford, T.D. and Cullingford, C.H.D.), 329-358, Academic Press (1976). Spötl, C. & Pavuza, R. Höhlenatmosphäre. in: Höhlen und Karst in Österreich (ed. by C. Spötl, L. Plan, E. Christian), Linz, Oberösterreichisches Landesmuseum, 123-138 (2016). Trüssel, M. Höhlenklima: Langzeit-Temperaturmessungen in der Schrattenhöhle, Melchsee-Frutt (CH-6064 Kerns OW). Akten des 14. Nationalen Kongresses für Höhlenforschung, 227-232, Höhlenforschergemeinschaft Region Hohgant (2019). Domínguez-Villar, D. et al.Is global warming affecting cave temperatures? Experimental and model data from a paradigmatic case study. Climate Dyn. 45 , 569-581 (2015). Colucci, R., Fontana, D., Forte, E. & Potleca, M., Guglielmin, M. Response of ice caves to weather extremes in the southeastern Alps, Europe, Geomorphology 261 , 1-11, http://dx.doi.org/10.1016/j.geomorph.2016.02.017 (2015). Racine, T., Reimer, P. J. & Spötl, C. Multi‐centennial mass balance of perennial ice deposits in Alpine caves mirrors the evolution of glaciers during the Late Holocene. Sci. Rep., 12 : 11374. https://doi.org/10.1038/s41598-022-15516-9 (2022). Filipponi, M. Das Klima der Schrattenhöhle: Auswertung einer 10-jährigen Messreihe, ETH Zürich, Institut für Atmosphäre und Klima, https://doi.org/10.3929/ethz-a-004662252 (2003). Šebela, S. & Turk, J. Comparison of historical and current temperatures in show caves (Slovenia). SN Appl. Sci. 4 , 1. https://doi.org/10.1007/s42452-021-04881-1(2022). Brunetti, M. et al. Climate variability and change in the Greater Alpine Region over the last two centuries based on multi-variable analysis. Int. J. Climatol. 29 , 2197–225.(2009). Kuhn M. & Olefs M. Elevation-Dependent Climate Change in the European Alps, Oxford Research Encyclopedias, Climate Science, 1-26. https://doi.org/10.1093/acrefore/9780190228620.013.762 (2020) Luetscher, M. & Jeannin, P.Y., Haeberli W. Ice caves as an indicator of winter climate evolution: a case study from the Jura Mountains, The Holocene 15 , 982-/993. https://doi.org/10.1191/0959683605hl872ra (2005). Gobiet, A. et al.21st century climate change in the European Alps - a review. Sci. Total Environm. 493 , 1138-1151 (2014). Pepin N.C., & Lundquist J.D. Temperature trends at high elevations: Patterns across the globe. Geophysical Research Letters 35 :, L14701. https://doi.org/10.1029/2008GL034026 (2008). Cigna, A.A. An analytical study of air circulation in caves. Intern. J. Speleol. 3, 41-54 (1968). Badino, G. Cave temperatures and global climatic change. Int. J. Speleol. 33 , 103-114 (2004). Luetscher, M., Lismonde, B. & Jeannin, P.Y. Heat exchanges in the heterothermic zone of a karst system: Monlesi cave, Swiss Jura Mountains, J. Geophys. Res. 113 , F02025, doi:10.1029/2007JF000892 (2008). Sedaghatkish, A. et al. Modeling heat transfer for assessing the convection length in ventilated caves. J. Geophys. Res.: Earth Surface 129 , e2024JF007646. https://doi.org/10.1029/2024JF007646 (2024). Cropley J. B. Influence of surface conditions on temperatures in large cave systems. Bull. Nat. Speleol. Soc. 27 , 1–10 (1965). Gabrovsek, F. How do caves breathe: The airflow patterns inckarst underground. PLoS ONE 18(4) : e0283767, https://doi.org/10.1371/journal.pone.0283767 (2023). Kukuljan, L., Gabrovšek, F. & Covington, M.D. The relative importance of wind-driven and chimney effect cave ventilation: Observations in Postojna Cave (Slovenia). Int. J. Speleol. 50 , 275-288. https://doi.org/10.5038/1827-806X.50.3.2392 (2021). Spötl, C. & Pavuza, R. Zur Bewetterung der Spannagelhöhle (Kat. Nr. 2411/1). Höhlenk. Mitt. Landesvereins für Höhlenkunde in Tirol 39 , 5-9 (2001). Luetscher, M. & Jeannin, P. Y. A process-based classification of alpine ice caves. Theor. Appl. Karst. 17 , 5–10 (2004). Spötl, C., Reimer, P.J. & Luetscher, M. Long-term mass balance of perennial firn and ice in an Alpine cave (Austria): Constraints from radiocarbon-dated wood fragments. The Holocene 24 , 165-175 (2014). Kern, Z. & Perșoiu, A. Cave ice - the imminent loss of untapped mid-latitude cryospheric palaeoenvironmental archives. Quat. Sci. Rev. 67 , 1-7 (2013). Colucci, R.R. & Guglielmin, M. Climate change and rapid ice melt: Suggestions from abrupt permafrost degradation and ice melting in an alpine ice cave. Progress Phys. Geogr. 43 , 561-573 (2019). Perșoiu, A. et al.. Record summer rains in 2019 led to massive loss of surface and cave ice in SE Europe. The Cryosphere 15 , 2383-2399 (2021). Villar, E. et al. Ambient temperature variations in the hall of paintings of Altamira cave due to the presence of visitors. Trans. British Cave Res. Ass. 11 , 99-103 (1984). Prelovšek, M., Šebela, S. & Turk, J. Carbon dioxide in Postojna Cave (Slovenia): spatial distribution, seasonal dynamics and evaluation of plausible sources and sinks. Environ. Earth. Sci. 77 , 289. https://doi.org/10.1007/s12665-018-7459-6 (2018). Fernandez-Cortes, A. et al. Detection of human-induced environmental disturbances in a show cave. Environ. Sci. Pollut. Res. 18 , 1037–1045. https://doi./10.1007/s11356-011-0513-5 (2011). Obleitner, F. & Spötl, C. The mass and energy balance of ice within the Eisriesenwelt cave, Austria. The Cryosphere 5 , 245–257, 2011. https://doi.org/10.5194/tc-5-245-2011 (2011). Perry, R.W. A review of factors affecting cave climates for hibernating bats in temperate North America. Environ. Reviews 21 , 28-39, https://doi.org/10.1139/er-2012-0042 (2012)66. Mann, H.B. Non-parametric tests against trend. Econometrica 13 , 163-171 (1945). Kendall, M.G. Rank Correlation Methods. 4th ed., Charles Griffin, London (1975). Theil, H. A rank-invariant method of linear and polynomial regression analysis. Nederl. Akad. Wetensch. Proc. 53 , 386–392 (1950). Sen, P. K. Estimates of the regression coefficient based on Kendall's tau. J. Amer. Stat. Ass. 63 (324), 1379–1389 (1968). The MathWorks Inc. MATLAB version: 9.13.0 (R2022b), Natick, Massachusetts: The MathWorks Inc. https://www.mathworks.com (2022). Additional Declarations No competing interests reported. Supplementary Files Obleitneretal2024supplement.pdf Cite Share Download PDF Status: Published Journal Publication published 10 Nov, 2024 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 03 Oct, 2024 Reviews received at journal 26 Sep, 2024 Reviews received at journal 17 Sep, 2024 Reviewers agreed at journal 16 Sep, 2024 Reviews received at journal 20 Aug, 2024 Reviewers agreed at journal 06 Aug, 2024 Reviewers agreed at journal 06 Aug, 2024 Reviewers invited by journal 06 Aug, 2024 Editor assigned by journal 06 Aug, 2024 Editor invited by journal 06 Aug, 2024 Submission checks completed at journal 03 Aug, 2024 First submitted to journal 25 Jul, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4801555","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":346431602,"identity":"18e45fe7-c62f-4927-8474-dbea151da059","order_by":0,"name":"Friedrich Obleitner","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABKUlEQVRIie3PsUrDQBjA8e84uCyHXb9Soa/whUJaoeCrpAhx6SC4CEqIFOJS8QX6EF0Fh4NAsvQBBAsmy00ROioW8VKEiEk7O9wf7ggHv/suADbbP8yJWAQSEAhAQG4+TCwHUrSPSPWb+ICV4VQTv4VUm1k/ZDdGmO0A6SW3+dvTqD8Enub+1SgcOvf6Wl6sB+AkOXy8N8nRZOYuNLqPkTgnf4V4Ms+8F0naAxkQmzennEoW96RCn5T0cBIj0nMgDEnGAFNzZcvDJLv7rMmXIa9aXO5IpwS2bSUxr0lUTRGCG+IBToG3T5l1FwrdZVL9S4rd5Srg3QXpgUBNyXHQQpxiU6qwT9kszTc3YYeylG3K7dp96JwVRTlukDr+90CYpQ4Am81ms+3vGz7JYehKBLzeAAAAAElFTkSuQmCC","orcid":"","institution":"Universität Innsbruck","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Friedrich","middleName":"","lastName":"Obleitner","suffix":""},{"id":346431604,"identity":"24f64ede-7459-40e0-a2c9-e649eb0afbb0","order_by":1,"name":"Christoph Spötl","email":"","orcid":"","institution":"Universität Innsbruck","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Christoph","middleName":"","lastName":"Spötl","suffix":""},{"id":346431606,"identity":"10e29fac-9b4d-462d-a574-787abd38e812","order_by":2,"name":"Martin Trüssel","email":"","orcid":"","institution":"Foundation Naturerbe Karst und Höhlen Obwalden","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Martin","middleName":"","lastName":"Trüssel","suffix":""}],"badges":[],"createdAt":"2024-07-25 11:36:15","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4801555/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4801555/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-024-78658-y","type":"published","date":"2024-11-10T15:57:45+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":63782681,"identity":"68410354-5aa0-4e87-a5a7-df80b9ff1889","added_by":"auto","created_at":"2024-09-02 09:58:56","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":785024,"visible":true,"origin":"","legend":"\u003cp\u003eLocations of the studied Alpine caves and the weather stations used for comparison.\u003c/p\u003e\n\u003cp\u003eAbbreviations of cave names: Hundsalm Eis- und Tropfsteinhöhle (HA), Rasslsystem (RS). Schrattenhöhle (SH), Spannagelhöhle (SP). Abbreviations of weather stations: Bad Eisenkappel (BA), Bonistock (BO), Engelberg (EN), Innsbruck (IN), Jenbach (JE), Mayrhofen (MA), Patscherkofel (PA), Pilatus (PI), Plattkopf (PA), Seeberg (SE), Zell-Pfarre (ZE). Base map: Google Earth.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/973ff65dbac0ad20ac974c1f.png"},{"id":63782683,"identity":"77e89c96-a01c-4809-b2d3-b98f9f89382b","added_by":"auto","created_at":"2024-09-02 09:58:56","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":195801,"visible":true,"origin":"","legend":"\u003cp\u003eTime series of cave temperatures observed in the four caves during 2000-2020. The partly lower resolution of some records before about 2005 is due to the use of different instruments.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/bdae12c9bb5062e1a4dd6023.png"},{"id":63782685,"identity":"cba34b3d-bd88-4585-a82e-8fae1d0b8e1b","added_by":"auto","created_at":"2024-09-02 09:58:57","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":72969,"visible":true,"origin":"","legend":"\u003cp\u003eBoxplots of observed cave temperatures (left) and outside temperatures calculated for the same elevation (right). Red lines denote medians, blue boxes indicate the 25th and 75th percentiles, whiskers represent extremes and red crosses outliers.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/439f2d32235ec4aa7902e16c.png"},{"id":63783275,"identity":"7a9d1a99-8f11-4afb-9dee-a222e62b0eed","added_by":"auto","created_at":"2024-09-02 10:06:57","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":65812,"visible":true,"origin":"","legend":"\u003cp\u003eQuasi-vertical gradients of mean temperatures measured in the investigated caves (2008-2010; red dots) compared to outside conditions, i.e., locally relevant AWS data (open triangles) and values interpolated to the elevation of cave observations (filled triangles). The dashed line shows a moist-adiabatic lapse rate (-0.5 °C 100 m\u003csup\u003e-1\u003c/sup\u003e).\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/dda46c3acfd6a8598b04dc1f.png"},{"id":63782687,"identity":"be647190-2f66-485b-a69e-a59a4690ee6e","added_by":"auto","created_at":"2024-09-02 09:58:57","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":190114,"visible":true,"origin":"","legend":"\u003cp\u003eMean annual cave temperatures (dots) compared to calculated temperatures outside of the caves (thick solid lines) and HISTALP data for selected Greater Alpine Regions (thin lines).\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/9a4e73644f06a32c68639b3c.png"},{"id":63782689,"identity":"293a8cf3-733f-4114-8f62-88dadfb29c3a","added_by":"auto","created_at":"2024-09-02 09:58:57","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":291984,"visible":true,"origin":"","legend":"\u003cp\u003eTime series of air temperatures in Spannagel cave. SP-PL and SP-outside denote the so far discussed data (red and grey lines) compared to SP-GH (green), SP-TG (yellow) and SP-EH (blue). The lower resolution of the measurements at the beginning of the SP-GH and SP-PL records is due to use of different instruments that were replaced at about 2005.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/54cc20022a5db53ff64ab59d.png"},{"id":63783273,"identity":"e0f0e54a-df0d-43af-b47b-5ef6216de746","added_by":"auto","created_at":"2024-09-02 10:06:56","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":195114,"visible":true,"origin":"","legend":"\u003cp\u003eTime series of air and rock temperatures in Hundsalm Eis- and Tropfsteinhöhle. HA-JG and HA-outside denote the previously presented data (red and grey lines), while HA-29 (blue) and HA-126 (yellow) refer to additional data measured in air and within rock in the upper, ice-bearing part of the cave.\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/6938ecda0c9637799e228a76.png"},{"id":63783274,"identity":"2737cebb-1ed8-49f2-97ce-72ac3772fbd6","added_by":"auto","created_at":"2024-09-02 10:06:57","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":70878,"visible":true,"origin":"","legend":"\u003cp\u003eMean air temperature observed at site HA-29 during subsequent stake measurements (red line) and concurrent surface height changes of ice located at the same location in HA (blue line).\u003c/p\u003e","description":"","filename":"8.png","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/1129719f18fe59dcb47465a4.png"},{"id":68749956,"identity":"224eaa81-7a86-4e56-904b-3f3f535b228c","added_by":"auto","created_at":"2024-11-11 16:08:02","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2556867,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/911e2be2-9192-4bff-8cfe-c98b079c13ee.pdf"},{"id":63782680,"identity":"413cebcf-cc12-43b2-b755-3d3cf026112d","added_by":"auto","created_at":"2024-09-02 09:58:56","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":717984,"visible":true,"origin":"","legend":"","description":"","filename":"Obleitneretal2024supplement.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4801555/v1/4e0dd144351fbf46bcb77f9a.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":" Climate warming has arrived in caves of the European Alps","fulltext":[{"header":"Introduction","content":"\u003cp\u003eAtmospheric warming is a global phenomenon of the Anthropocene. Global mean temperatures and temperatures averaged over land areas have warmed by 1.1\u0026deg;C and 1.6\u0026deg;C from 1850\u0026ndash;1900 to 2011\u0026ndash;2020, respectively \u003csup\u003e[1]\u003c/sup\u003e. Mountain regions warmed at an even higher rate, which is due to reduced snow and ice feedbacks, reduced local cooling due to shrinking glacier extent, and local land use changes \u003csup\u003e[2]\u003c/sup\u003e. The European Alps (hereafter Alps for short) have warmed twice as much as the global average between the late 19th and the early twenty-first century \u003csup\u003e[3, 4]\u003c/sup\u003e. In recent decades, the average warming rate in the Alps has been 0.3\u0026deg;C per decade, outpacing the global rate of 0.2 \u0026deg;\u003c/p\u003e \u003cp\u003eC per decade \u003csup\u003e[5]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe environmental impact of this strong atmospheric warming of the Alps is most prominently manifested by the loss of glacier extent and volume as a result of the rise of the equilibrium line altitude \u003csup\u003e[6, 7]\u003c/sup\u003e. Parallel to this reduction in surface ice masses, high-Alpine regions have also witnessed a degradation of permafrost \u003csup\u003e[8,9]\u003c/sup\u003e. Rising temperatures have also affected surface waters in the Alps \u003csup\u003e[10]\u003c/sup\u003e. A recent study of two catchments in the Tyrolian Alps identified mean warming trends of local streams of 0.24 and 0.45\u0026deg;C per decade for the time period 1977 to 2020 \u003csup\u003e[11]\u003c/sup\u003e. This signal has also propagated into the subsurface and led to an increase in groundwater temperatures. Shallow groundwaters in Austria (average depth 7\u0026thinsp;\u0026plusmn;\u0026thinsp;4 m) increased on average by 0.37\u0026deg;C per decade between 1994 and 2013 \u003csup\u003e[12,13]\u003c/sup\u003e. A study in the northern Alpine foreland (Bavaria, Germany) reported a decrease of the warming amplitude with depth, from 0.28\u0026deg;C per decade in 20 m to 0.09\u0026deg;C per decade in 60 m depth between the early 1990s and 2019 \u003csup\u003e[14]\u003c/sup\u003e. Rising temperatures have also been detected inside High-Alpine glaciers, e.g., in the Mont Blanc area, revealing a mean warming rate of 0.14\u0026deg;C per decade between 1900 and 2004 \u003csup\u003e[15]\u003c/sup\u003e. Perrier et al.\u003csup\u003e[16]\u003c/sup\u003e provided evidence that long-term changes in surface climate have also propagated into urban underground.\u003c/p\u003e \u003cp\u003eSurprisingly little is known, however, if and to what extent this atmospheric warming signal has already \u0026ldquo;arrived\u0026rdquo; in the interior of cave systems. Caves are abundant in many part of the Alps where karstifiable rocks (chiefly limestone, dolomite and/or marble) are present. The Alps not only host some of the largest caves in Europe (e.g., H\u0026ouml;lloch, Siebenhengste-Hohgant H\u0026ouml;hlensystem, Sch\u0026ouml;nberg-H\u0026ouml;hlensystem), but also caves with ice deposits many of which showing strong signs of retreat, e.g., Dachstein Mammuth\u0026ouml;hle \u003csup\u003e[17]\u003c/sup\u003e, Hundsalm Eis- und Tropfsteinh\u0026ouml;hle \u003csup\u003e[18]\u003c/sup\u003e and Obstanser Eish\u0026ouml;hle \u003csup\u003e[19]\u003c/sup\u003e. Warming of underground cavities has been shown to have also significant consequences for subterranean biota and ecosystems \u003csup\u003e[20, 21]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eHere, we document for the first time decadal-scale warming in the interior of caves across the Alps. Although local time series of temperature measurements in different caves (and parts thereof) exist since many years \u003csup\u003e[22, 23, 24]\u003c/sup\u003e, this is the first study that systematically evaluates temperature measurements obtained in four caves in Switzerland and Austria. The aim of this study is to document and quantify the ongoing warming in the interior of these Alpine caves and to encourage follow-up research, in order to take full advantage of the wealth of temperature data available for other caves in this mountain belt.\u003c/p\u003e\n\u003ch3\u003eThe microclimate of caves\u003c/h3\u003e\n\u003cp\u003eAlthough hidden in the underground, caves are part of the Earth\u0026rsquo;s surface and mass and energy are exchanged between them in a variety of ways. The most important processes are advection of air and water. Diurnal and seasonal temperature variations that characterize the Earth surface are attenuated in the cave environment. Ignoring latent heat and humidity effects, the attenuation of the diurnal and seasonal temperature amplitudes follows an exponential function from the entrance to the interior of caves \u003csup\u003e[25]\u003c/sup\u003e. In the near-entrance part of caves \u0026ndash; the heterothermic zone \u0026ndash; diurnal and seasonal temperature cycles are well expressed. Deeper into the cave, the attenuation approaches an asymptotic value. This is referred to as the homothermic zone, where diurnal cycles are below detection and also the seasonal temperature amplitude typically does not exceed a few tenths of a degree. This stable temperature is commonly close to the mean annual surface temperature, except where local effects play a role.\u003c/p\u003e \u003cp\u003eA compilation of data from 40 caves in Austria located between 403 and 2347 m a.s.l. revealed that most caves show interior cave air temperatures close to the multi-annual mean of the outside temperature \u003csup\u003e[26]\u003c/sup\u003e. Anomalously low temperatures were measured in caves that are characterized by restricted ventilation and act as cold traps in winter. This study also reported preliminary results of long-term measurements in one Alpine cave (Spannagelh\u0026ouml;hle), where temperatures had increased since the year 2000 in a near-linear manner. Preliminary data of a similar monitoring work in caves of the Swiss Alps documented rising temperatures since the year 1990 \u003csup\u003e[27]\u003c/sup\u003e. In contrast to these well-ventilated caves, a study in a side passage of Postojna Cave (Slovenia) with a dead-end gallery with limited air advection (in contrast to the rest of the system) concluded that there the external warming signal was largely transferred into the cave by thermal conduction, resulting in a delay of 20\u0026ndash;25 years \u003csup\u003e[28]\u003c/sup\u003e. An extreme endmember is a monitoring site in the interior of B\u0026auml;renschacht, a 60 km-long cave network spanning 946 m of vertical extent. The measurement site is 700 m below the surface with limited air exchange and showed a constant air temperature within 0.1\u0026deg;C between 1999 and 2017 (W. Janz, pers. comm.).\u003c/p\u003e\n\u003ch3\u003eStudy sites\u003c/h3\u003e\n\u003cp\u003eFour caves in different climatological regions across the Alps were selected showing a range in size and elevation (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Rasslsystem (RS), Spannagelh\u0026ouml;hle (SP) and Hundsalm Eis- und Tropfsteinh\u0026ouml;hle (HA) are located in different parts of the Austrian Alps, while Schrattenh\u0026ouml;hle (SH) is located in the Swiss Alps (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eKey parameters of the studied caves.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSpannagelh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSchrattenh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHundsalm Eis- und Tropfsteinh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRasslsystem\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAbbreviation\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSP\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSH\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHA\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRS\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eCave register number\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2515/1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eOW 4/61 (22363)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1266/1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e3925/9\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eCoordinates (lat./long.)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e47.08\u0026deg;/11.67\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e46.77\u0026deg;/8.27\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e47.55\u0026deg;/12.03\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e46.51\u0026deg;/14.54\u0026deg;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLength of cave (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e12000\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e19718\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e264\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1085\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eVertical extension (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e329\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e573\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e51\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e161\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eClassification of cave system\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003edynamic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003edynamic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003estatic with ice\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003edynamic\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eNumber of entrances\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eElevation of (upper) entrance (m a.s.l.)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2523\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e2090\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1495\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1085\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eName of principle cave monitoring site\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePorzellanladen\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eWindkluft\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eJugendgang\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePerlenhalle\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAbbreviation\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSP-PL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSH-WK\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHA-JG\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRS-PH\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eElevation of monitoring site (m a.s.l.)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2408\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1743\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e1480\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e1070\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRock overburden of monitoring site (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e17\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e60\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e95\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eDistance of monitoring site to nearest entrance\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e840\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e250\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e85\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e270\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eRelevant outside weather stations\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eJE, PK, PL\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePI, EN, BO\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eJE, PK\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eBA, SE, ZE\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eSpannagelh\u0026ouml;hle is a 12 km-long system in the western Zillertal Alps (Tyrol; Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) with three upper entrances and unknown lower entrance(s) extending over 329 m in the vertical dimension (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). The main cave parts are well ventilated and characterized by upward directed air flow in winter and a reversed flow in summer. The main monitoring site, Porzellanladen (SP-PL), is located in the interior of the cave very close to the main conduit where the air is flowing. Additionally, a site in a dead-end side gallery (Tropfsteingang, SP-TG) and one located in the main air flow but closer to the main upper entrance (Gaunhalle, SP-GH) were included in this study (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). We also present data from a small cave that opens 27 m beneath the main entrance of Spannagel Cave and contained an ice body when discovered in 1998. This so-called Spannagel-Eish\u0026ouml;hle (SP-EH) shows a descending geometry and has lost its ice by the year 2002.\u003c/p\u003e \u003cp\u003eHundsalm Eis- und Tropfsteinh\u0026ouml;hle is located in the Northern Calcareous Alps of Tyrol (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) and classifies as a sag-type cave given its subvertical geometry and the lack of a lower entrance (Supplementary Fig. S2). HA is also exceptional because the cave consists of two compartments, a higher one that also contains perennial ice, and a lower section, which is ice-free. The two cave parts were artificially connected in 1996 and then subsequently thermally re-separated using an airlock. Our study concerns primarily the lower section (location Jugendgang, HA-JG) whose air temperature is close to the mean air temperature outside (in contrast to the cold trap in the upper cave part). Temperature data obtained in the main shaft and in a short borehole in the upper part of the cave are also considered in this study (Supplementary Fig. S2).\u003c/p\u003e \u003cp\u003eRasslsystem is one of several caves in the Obir Massif of the Northern Karawanks in SE Austria (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Its air flow regime was modified as a result of mining activities which came to an end at the beginning of the 20th century. An abandoned mining adit acts as the lower entrance of cave system (Supplementary Fig. S3) with strong inward and outward air flow during the cold and warm season, respectively, while the pathways of the air in the upper inaccessible parts of the cave are largely unknown \u003csup\u003e[22]\u003c/sup\u003e. The monitoring site is Perlenhalle (RS-PH), a high and elongated terminal chamber in the innermost part of Rasslsystem with no measurable air flow.\u003c/p\u003e \u003cp\u003eSchrattenh\u0026ouml;hle is located in the Melchsee-Frutt region of Switzerland (Obwalden; Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) and comprises an extensive network of mostly well-ventilated galleries extending over 573 m and showing multiple entrances (Supplementary Fig. S4). The monitoring site, Windkluft (SH-WK), is located in the main air flow route 250 m from the next entrance, while the subsidiary site SH-WG (Wermutgang) is located close to the uppermost entrance of the Schrattenh\u0026ouml;hle cave system.\u003c/p\u003e"},{"header":"RESULTS","content":"\u003cp\u003eWe examined daily averages of data measured in the four caves and in the free atmosphere of nearby meteorological stations. The data were subjected to quality control and restricted to the period providing the most consistent records (2000\u0026ndash;2020). The main analyses concerned climatological features and relationships with external conditions, trends using different methods including seasonal aspects, spatial and temporal representativeness and the relationship with cave ice developments.\u003c/p\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eClimatological characteristics\u003c/h2\u003e \u003cp\u003eThe analyzed records are characterized by a) different mean values, b) regular annual cycles of comparable magnitude, and c) continuously increasing temperatures at all observation sites (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe average conditions in the caves as well as in their outside atmosphere are summarized in Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, which refers to the central range of the overall investigation period (2008\u0026ndash;2010).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eMean and standard deviation of air temperature inside and outside of the investigated caves and respective differences (2008\u0026ndash;2010). Bracketed numbers in the lowest panel denote Pearson correlation coefficients between air temperature in the caves and outside.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSP-PL\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSH-WK\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHA-JG\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRS-PH\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eElevation of cave sites [m a.s.l.]\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2408\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e1743\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1480\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1070\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eCave temperature [\u0026deg;C]\u003c/em\u003e\u003c/p\u003e \u003cp\u003emean annual\u003c/p\u003e \u003cp\u003ewinter (open period)\u003c/p\u003e \u003cp\u003esummer (closed period)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;1.82\u0026thinsp;\u0026plusmn;\u0026thinsp;0.066\u003c/p\u003e \u003cp\u003e+\u0026thinsp;1.77\u0026thinsp;\u0026plusmn;\u0026thinsp;0.044\u003c/p\u003e \u003cp\u003e+\u0026thinsp;1.86\u0026thinsp;\u0026plusmn;\u0026thinsp;0.046\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e+\u0026thinsp;3.17\u0026thinsp;\u0026plusmn;\u0026thinsp;0.064\u003c/p\u003e \u003cp\u003e+\u0026thinsp;3.12\u0026thinsp;\u0026plusmn;\u0026thinsp;0.046\u003c/p\u003e \u003cp\u003e+\u0026thinsp;3.21\u0026thinsp;\u0026plusmn;\u0026thinsp;0.044\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e+\u0026thinsp;4.15\u0026thinsp;\u0026plusmn;\u0026thinsp;0.054\u003c/p\u003e \u003cp\u003e+\u0026thinsp;4.13\u0026thinsp;\u0026plusmn;\u0026thinsp;0.055\u003c/p\u003e \u003cp\u003e+\u0026thinsp;4.17\u0026thinsp;\u0026plusmn;\u0026thinsp;0.043\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e+\u0026thinsp;5.47\u0026thinsp;\u0026plusmn;\u0026thinsp;0.022\u003c/p\u003e \u003cp\u003e+\u0026thinsp;5.46\u0026thinsp;\u0026plusmn;\u0026thinsp;0.023\u003c/p\u003e \u003cp\u003e+\u0026thinsp;5.47\u0026thinsp;\u0026plusmn;\u0026thinsp;0.021\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eOutside temperature [\u0026deg;C] *\u003c/em\u003e\u003c/p\u003e \u003cp\u003emean annual\u003c/p\u003e \u003cp\u003ewinter (open period)\u003c/p\u003e \u003cp\u003esummer (closed period)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e+\u0026thinsp;0.40\u0026thinsp;\u0026plusmn;\u0026thinsp;7.55\u003csup\u003e1\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e-6.07\u0026thinsp;\u0026plusmn;\u0026thinsp;5.18\u003c/p\u003e \u003cp\u003e+\u0026thinsp;5.66\u0026thinsp;\u0026plusmn;\u0026thinsp;5.09\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e+\u0026thinsp;3.75\u0026thinsp;\u0026plusmn;\u0026thinsp;7.05\u003csup\u003e2\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e-1.42\u0026thinsp;\u0026plusmn;\u0026thinsp;4.96\u003c/p\u003e \u003cp\u003e+\u0026thinsp;8.79\u0026thinsp;\u0026plusmn;\u0026thinsp;4.75\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e+\u0026thinsp;4.18\u0026thinsp;\u0026plusmn;\u0026thinsp;7.60\u003csup\u003e3\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e-1.65\u0026thinsp;\u0026plusmn;\u0026thinsp;4.99\u003c/p\u003e \u003cp\u003e+\u0026thinsp;9.93\u0026thinsp;\u0026plusmn;\u0026thinsp;4.83\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e+\u0026thinsp;4.94\u0026thinsp;\u0026plusmn;\u0026thinsp;8.21\u003csup\u003e4\u003c/sup\u003e\u003c/p\u003e \u003cp\u003e-1.39\u0026thinsp;\u0026plusmn;\u0026thinsp;5.37\u003c/p\u003e \u003cp\u003e+\u0026thinsp;11.2\u0026thinsp;\u0026plusmn;\u0026thinsp;5.20\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eTemperature difference in-out [\u0026deg;C]\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e(correlation coefficient)\u003c/em\u003e\u003c/p\u003e \u003cp\u003emean annual\u003c/p\u003e \u003cp\u003ewinter (open period)\u003c/p\u003e \u003cp\u003esummer (closed period)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e+1.42 (+\u0026thinsp;0.72)\u003c/p\u003e \u003cp\u003e+6.96 (+\u0026thinsp;0.50)\u003c/p\u003e \u003cp\u003e-4.04 (+\u0026thinsp;0.31)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e-0.58 (+\u0026thinsp;0.66)\u003c/p\u003e \u003cp\u003e+\u0026thinsp;4.54 (+\u0026thinsp;0.35)\u003c/p\u003e \u003cp\u003e-5.58 (+\u0026thinsp;0.25)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e-0.03 (+\u0026thinsp;0.46)\u003c/p\u003e \u003cp\u003e\u0026acute;+5.78 (+\u0026thinsp;0.44)\u003c/p\u003e \u003cp\u003e-5.76 (-0.03)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e+\u0026thinsp;0.53 (-0.11)\u003c/p\u003e \u003cp\u003e+\u0026thinsp;6.85 (-0.06)\u003c/p\u003e \u003cp\u003e-5.73 (-0.23)\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e* at the elevation of measurements in the cave\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e1\u003c/sup\u003e calculated from AWS data JE and PA, bias from comparison with PA\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e2\u003c/sup\u003e calculated from AWS data EN and PI, bias from comparison with BO\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e3\u003c/sup\u003e calculated from AWS data JE and PA, bias from comparison with HA-entrance\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003e4\u003c/sup\u003e calculated from AWS data BA and SE, bias from comparison with ZE\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe average cave temperatures correspond to the elevations of the cave sites. Lowest values are observed at SP-PL (+\u0026thinsp;1.8\u0026deg;C) being located at an elevation of 2408 m a.s.l. SH-WK (1743 m a.s.l.) is characterized by an average temperature of +\u0026thinsp;3.2\u0026deg;C, HA-JG (1480 m a.s.l.) by +\u0026thinsp;4.2\u0026deg;C and RS-PH by +\u0026thinsp;5.5\u0026deg;C (1070 m a.s.l.). The mean temperatures are positive at all cave sites, as well as those in the atmosphere outside of the caves, which range from +\u0026thinsp;0.4\u0026deg;C at the elevation of SP-PL to +\u0026thinsp;4.9\u0026deg;C (for RS-PH), respectively. The variability of cave temperatures is almost two orders of magnitude smaller compared to that outside (\u0026thinsp;\u0026lt;\u0026thinsp;\u0026plusmn;\u0026thinsp;0.1\u0026deg;C vs. ca. \u0026plusmn;8\u0026deg;C) and increases with elevation (from \u0026plusmn;\u0026thinsp;0.02 to \u0026plusmn;\u0026thinsp;0.07\u0026deg;C). Thus, RS-PH stands out by an almost homothermic regime, i.e., quasi-constant temperatures despite of the somewhat larger outside variability. Figure\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e graphically summarizes these basic characteristics in terms of boxplots.\u003c/p\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e includes information about the average conditions during winter and summer, which in the context of cave microclimate are also referred to as the so-called open and closed periods \u003csup\u003e[18]\u003c/sup\u003e. In this respect it is remarkable that the average cave temperatures are positive for both periods. Moreover, the cave sites experienced strongly attenuated seasonal temperature differences. Winters are just about 0.1\u0026deg;C colder than summer while in the free atmosphere at the same elevation where winters are about 10\u0026deg;C colder than summer and have negative mean temperatures. During winter the cave temperatures are about 5\u0026deg;C higher than outside, while caves remain colder than outside during summer.\u003c/p\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e also shows that the average cave temperatures observed in HA and SH are slightly lower than the outside temperatures at the same elevation, while SP and RS are warmer than outside. This is corroborated by consideration of the observed AWS temperatures (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e), although the latter represent a less reliable measure of local outside conditions due to being located at different elevations. The observed cave temperatures follow a quasi-vertical gradient of about \u0026minus;\u0026thinsp;0.4\u0026deg;C 100 m\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e (red dots in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e), which is somewhat smaller than the lapse rate of outside temperatures in the relevant elevation range, i.e., 1000\u0026ndash;2500 m a.s.l. and the moist adiabatic lapse rate (blue triangles and dotted line), respectively. However, it should be noted that the cave sites are located in different climatological regimes characterizing the Swiss Alps (SH), the Northern Calcareous Alps (HA), the Central Alps (SP) and the southeastern part of the Alps (RS). Judged by AWS data, temperature gradients are different along valleys in the vicinity of the investigated caves (-0.57\u0026deg;C 100 m\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in the Central Alps compared to -0.43 and \u0026minus;\u0026thinsp;0.41\u0026deg;C 100 m\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in the western and eastern and western fringes of the Alps, respectively).\u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e also suggests that cave temperatures are correlated with the outside temperatures at the same elevation, which can be considered as an indicator of a physically based response of cave temperatures to changes of the outside atmosphere. This is supported by Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e (lowest panel) which also shows that the correlation of daily cave temperatures with concurrent outside values is stronger at the high-elevation sites and during winter, RS being outstanding in this respect. These results are related to cave-specific ventilation patterns which themselves strongly depend on their morphological features (concerning entrances essentially) as will be discussed below.\u003c/p\u003e \u003cp\u003eTrend analysis\u003c/p\u003e \u003cp\u003eLong-term temperature trends were investigated by calculation of linear trends based on ordinary least squares regression and the Theil\u0026ndash;Sen trend estimator as well as using Mann-Kendall trend tests. The methods are applied to quality-controlled cave temperatures covering two decades as well as to AWS data representing temperatures at elevations below and above the cave sites.\u003c/p\u003e \u003cp\u003eVisual inspection of the observed records suggests that cave temperatures have systematically increased at all observation sites (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). This is confirmed by quantitative analyses for the period 2000\u0026ndash;2020 (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e), which reveals positive trends for all investigated sites, i.e., in the caves themselves as well as in the outside atmosphere. Cave temperatures increased by about\u0026thinsp;+\u0026thinsp;0.2\u0026deg;C per decade, with RS and SP showing higher values than SH or HA.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eDecadal temperature trends for the period 2000\u0026ndash;2020 (linear and Theil-Sen, the latter in bold). Outside temperatures refer to the elevation of the cave sites. MK\u0026thinsp;+\u0026thinsp;qualifies the calculated trends to be significant according to the Mann-Kendall test.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSpannagelh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSchrattenh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHundsalm Eis- und Tropfsteinh\u0026ouml;hle\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRasslsystem\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003eCave\u003c/b\u003e [\u0026deg;C per decade]\u003c/p\u003e \u003cp\u003etotal period\u003c/p\u003e \u003cp\u003e\u003cem\u003eopen periods\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003eclosed periods\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSP-PL\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.24 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.26 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.22 | +0.25 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.25 | +0.29 (MK-)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSH-WK\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.14 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.15 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.14 | +0.15 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.13 | +0.14 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eHA-JG\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.20 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.12 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.21 | +0.21 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.20 | +0.19 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRS-PH\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.24 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.28 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.25 |+0.28 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.24 |+0.26 (MK-)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003eOutside\u003c/b\u003e [\u0026deg;C per decade]\u003c/p\u003e \u003cp\u003etotal period\u003c/p\u003e \u003cp\u003e\u003cem\u003eopen periods\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003eclosed periods\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eat elevation of SP-PL\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.54 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.53 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.61 | +0.59 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.22 | +0.23 (MK-)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eat elevation of SH-WK\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.66 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.72 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.63 | +0.94 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.48 | +0.53 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eat elevation of HA-JG\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.59 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.57 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.62 | +0.62 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.29 | +0.28 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eat elevation of RS-PH\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.58 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.53 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.76 | +0.74 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.10 | +0.03 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003eAWS low\u003c/b\u003e [\u0026deg;C per decade]\u003c/p\u003e \u003cp\u003etotal period\u003c/p\u003e \u003cp\u003e\u003cem\u003eopen periods\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003eclosed periods\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eJE\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.68 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.63 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.69 | +0.68 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.36 | +0.32 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003eEN\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.65 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.61 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.68 | +0.65 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.36 | +0.32 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eJE\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.68 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.63 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.69 | +0.68 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.36 | +0.32 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eBA\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.63 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.59 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.64 | +0.68 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.30 | +0.27 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cb\u003eAWS high\u003c/b\u003e [\u0026deg;C per decade]\u003c/p\u003e \u003cp\u003etotal period\u003c/p\u003e \u003cp\u003e\u003cem\u003eopen periods\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003eclosed periods\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ePA\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.54 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.53 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.61 | +0.60 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.22 | +0.23 (MK-)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePI\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.70 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.77 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.68 | +0.76 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.36 | +0.31 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003ePA\u003c/p\u003e \u003cp\u003e+\u0026thinsp;0.55 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;0.53 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.61 | +0.60 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.22 | +0.23 (MK-)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSE\u003c/p\u003e \u003cp\u003e+\u0026thinsp;1.11 l\u0026thinsp;\u003cb\u003e+\u0026thinsp;1.17 (MK+)\u003c/b\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.96 | +0.99 (MK+)\u003c/em\u003e\u003c/p\u003e \u003cp\u003e\u003cem\u003e+\u0026thinsp;0.88 | +0.93 (MK+)\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe outside temperatures at the elevation of the cave sites increased at higher rates of about 0.5\u0026deg;C per decade. This value is independent of the applied calculation method and there is no clear evidence regarding elevation dependencies or regional patterns. It is notable, however, that temperatures at low-lying weather stations are characterized by rather consistent trends in the order of +\u0026thinsp;0.6\u0026deg;C per decade, while some distinction can be observed regarding the investigated high-elevation sites. Hence there is an indication that the eastern and western fringes of the Alps (i.e., the RS and SH regions) experienced enhanced warming compared to the Central (SP) and the Northern Calcareous Alps (HA), which however, is not reflected by the cave data.\u003c/p\u003e \u003cp\u003eWe also examined seasonal trends representing open (Nov-May) and closed periods (Jun-Oct), which provide evidence that the outside temperatures increased stronger during winter compared to summer (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, open vs. closed periods). Naturally, this also concerns the underlying AWS data and contrasts the cave conditions, where significant differences between winter and summer temperature trends can hardly be discerned.\u003c/p\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e also shows the results of the Mann-Kendall tests to explore the temporal consistency of the observed cave temperature trends. There is evidence that all investigated cave records are characterized by statistically significant trends, which also holds true for the investigated sub-periods (open and closed periods) and the calculated outside temperatures. The significance of positive trends can also be stated for the weather stations being used for the calculation of the outside temperatures, with the exception of summer temperatures at Patscherkofel.\u003c/p\u003e \u003cp\u003eThese results are based on observations representing the two decades providing the most consistent data (2000\u0026ndash;2020). However, the underlying time series have different length and episodic gaps (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), which in principle can affect the robustness of the so far presented analysis. Moreover, the spatial representativeness of the presented cave records may be influenced by local effects (e.g., ventilation conditions, episodic water influx or icing) and diverse observation problems (such as unrepresentative sites, different instruments and/or failure due to e.g. icing). A rigorous analysis of respective effects on our trend analysis is hampered by the scarceness of appropriate cave records over longer time spans. To at least qualitatively strengthen the respective validity of the presented results we considered some additional temperature records from other positions in these caves. Corresponding results are addressed in the following sections and we here just mention strong evidence that all investigated cave sites experienced a warming since at least 1990 and beyond 2020 until the recent days.\u003c/p\u003e \u003cp\u003eIn this context, we also refer to temperature data across the so-called Greater Alpine Region (GAR), which represents a gridded data set across the Alps derived from long-term homogenized station data \u003csup\u003e[3]\u003c/sup\u003e which allows putting the presented results in a broader temporal and spatial context. We here consider a limited ensemble of these data, i.e., back to 1970 as well as low/high-elevation, West and Southeast regions, respectively. Figure\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e demonstrates that these records are characterized by consistent and fairly linear trends in the order of +\u0026thinsp;0.5\u0026deg;C per decade, with indication of enhanced warming in the lowlands compared to high-lying regions (+\u0026thinsp;0.56 vs. +0.48\u0026deg;C per decade). The analysis further shows that the observed cave records as well as the calculated outside temperatures fit these regional developments, which essentially confirms their spatial and temporal representativeness. Note that the Austrian caves are located in HISTALP-regions West and Southeast. It may be added that Jungfraujoch, i.e., a Swiss mountain station some 50 km south of SH, is included in the GAR-High-Level data, where a warming rate of +\u0026thinsp;0.45\u0026deg;C per decade was observed during the period 1970\u0026ndash;2020.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eCase study SP\u003c/h2\u003e \u003cp\u003eSP has several entrances at different elevations, which has implications for the ventilation regime of the cave. Figure\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e (red line) again shows the previously presented SP-PL data and outside temperatures which are both characterized by consistent warming (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Addressing the question of respective spatial representativeness, we explored whether trends can also be detected in two other parts of this cave. Site SP-GH is located closer to the (upper) entrance while site SP-TG is located deeper into the cave, in a dead-end side passage (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). SP-GH largely mirrors the SP-PL data concerning the overall level (+\u0026thinsp;1.68 vs. +1.82\u0026deg;C), the seasonal variability and the trend (+\u0026thinsp;0.36 vs. +0.24\u0026deg;C per decade). The slightly higher temperatures at SP-PL can be attributed to the larger distance from the main (upper) entrance (840 m compared to 470 m for SP-GH). SP-TG is characterized by a different thermal regime with higher temperatures (+\u0026thinsp;2.16\u0026deg;C on average) and a strikingly reduced seasonal variability. This can be related to its location in kind of chamber deeper down in the cave and off the main air flow route. It is remarkable that also this remote cave section experienced a warming (+\u0026thinsp;0.27\u0026deg;C per decade) which according to the MK-test is statistically significant.\u003c/p\u003e \u003cp\u003eWe further studied temperature records from a small separate cave close to the entrance of SP (Spannagel-Eish\u0026ouml;hle, SP-EH) which contained perennial ice until the year 2002. Figure\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e (blue line) shows that also this cave experienced a warming during the last two decades. The overall lower air temperature as well as the tenuous trend at the beginning of the record may have been related to the presence of ice. The seasonal amplitude of air temperature is comparatively large due the small size and ventilation of this cave. Summarizing, there is strong evidence that Spannagel cave as a whole and nearby caves were affected by a consistent warming over the last two decades at least.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eCase study HA\u003c/h2\u003e \u003cp\u003eHA is special in the context of the four studied Alpine caves in that its ventilation regime is much more restricted given its sag-type geometry and the presence of perennial ice (Supplementary, Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). This leads to a cold-trap setting in the upper part of the cave, which will be examined in more detail in this chapter thereby also considering wall-rock temperature data.\u003c/p\u003e \u003cp\u003eFirst, we again note the progressive damping of outside temperature changes as they progress into the cave, which is a key feature of cave microclimate and strongly depends on ventilation characteristics. In sag-type caves (as HA) the damping is most pronounced during summer when thermal stratification prohibits effective exchange of air between the cave interior and the outside. Conversely, cold air can relatively easily penetrate into the deeper parts of the cave. We observed a seasonal temperature variation of about 30\u0026deg;C outside, while at a depth of 29 m below the main entrance (site HA-29) the air temperature varies by just about 5\u0026deg;C. We note that the conditions in this cave section are not only preconditioned by effective cold air pooling but also by the existence of perennial ice (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). The lower part of the cave is thermally separated from the upper one (site HA-JG), contains no ice and is characterized by a homothermic regime (air temperature hardly varies on an annual time scale) at comparatively high temperatures (average\u0026thinsp;+\u0026thinsp;4.3 compared to +\u0026thinsp;0.1\u0026deg;C at site HA-29). Figure\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e also shows long-term observations of rock temperatures at a distance of 126 cm behind the cave wall next to site HA-29. The average rock temperature is +\u0026thinsp;0.4\u0026deg;C, thus being slightly higher compared to the cave air at the same depth (+\u0026thinsp;0.1\u0026deg;C). During winter the rock does not cool as much as the air, while air and rock are thermally adjusted during summer. Hence there is a systematic difference in seasonal variability. Closer inspection also reveals that changes in rock temperature systematically lag behind those in the air \u003csup\u003e[18]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eFigure \u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e finally supports that each of the considered HA sites experienced a systematic warming throughout the investigation period 2010\u0026ndash;2020. This is corroborated by the calculated trends (0.24 / 0.79 / 0.99\u0026deg;C per decade for HA-JG / HA-29 / HA-126) and the according MK-trend tests (all attesting statistical significance). It can thus be concluded that the results, which refer to HA-JG, can serve as a reference for the long-term evolution of air temperature in this cave. Moreover, there is indication that trends are stronger in cave sections containing ice (HA-29) and rock temperatures (HA-126) compared to relatively remote cave sections without ice (HA-JG).\u003c/p\u003e \u003cp\u003eSurface height changes of ice in HA (Supplementary Fig. S2) were monitored by annual stake readings close to site HA-29. This effort was motivated by the expectation that long-term warming traces in enhanced melt of ice, as is generally known for snow and glaciers outside \u003csup\u003e[5]\u003c/sup\u003e. Indeed, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e (blue line) shows a continuous loss of ice, which was accelerated during the second half of the investigation period (2008\u0026ndash;2018). The values of air temperature averaged over time between the stake readings (black line) show a transition from initially negative values towards mostly positive values. This analysis confirms a strong physical relationship between long-term temperature evolution and the development of ice in caves, as is known from ice caves elsewhere \u003csup\u003e[29, 18, 19, 30]\u003c/sup\u003e .\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eOngoing climate change is leading to a regional and global rise in the temperature of the atmosphere. The evidence is essentially based on a sophisticated synthesis of observations from weather stations, radiosondes and satellites, as well as numerical models \u003csup\u003e[1]\u003c/sup\u003e. Equally important evidence comes from indirect sources relating to corresponding changes in the ocean (e.g., rising sea-surface temperatures), in the cryosphere (e.g., vanishing glaciers and ice sheets and thawing permafrost) and in the biosphere (e.g., shifting vegetation belts). This raises the question of whether climate change has also \u0026ldquo;arrived\u0026rdquo; in the caves, which is investigated in this study on the basis of long-term observations from caves in the Alps. The data from four caves in different parts of this mountain range provide ample evidence that this is indeed the case. Principal data were selected considering their quality as well as their spatial and temporal representativeness. The main analysis is based on data covering two decades (2000\u0026ndash;2020) and linear trends were calculated and tested for statistical significance. Extended temporal representativeness of key results is shown back to the 1990s, however, as well as correspondence to developments in other cave sections and regional climate changes outside the caves.\u003c/p\u003e\n\u003ch3\u003eContext to previous work\u003c/h3\u003e\n\u003cp\u003eThe analysis of the studied cave temperature records reveals a consistent picture in terms of linear trends ranging from +\u0026thinsp;0.1 to +\u0026thinsp;0.2\u0026deg;C per decade (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). These figures are backed by the application of two calculation methods (least squares and Theil-Sen) and the statistical significance is confirmed by the Mann-Kendall test. Additional temperature records from each cave show that the principal records are representative for major parts of the investigated caves and not just for individual sites. These results overall agree with previous findings which are, however, mostly based on shorter time series and less elaborated analyses. Sp\u0026ouml;tl and Pavuza \u003csup\u003e[26]\u003c/sup\u003e reported temperature trends for Spannagel cave of +\u0026thinsp;0.5 and +\u0026thinsp;0.3\u0026deg;C per decade (2000\u0026ndash;2015) obtained from sites GH and TG, respectively. These sites are located above and below site SP-PL, which is considered in this work. Although these figures refer to different periods, they also indicate that the sections closer to the upper entrance have experienced greater warming.\u003c/p\u003e \u003cp\u003eTr\u0026uuml;ssel \u003csup\u003e[27]\u003c/sup\u003e analyzed temperature records from Schrattenh\u0026ouml;hle and reported a warming of 0.7\u0026deg;C for Windkluft (SH-WK in this study) for the period 1990\u0026ndash;2018, which corresponds to a trend in order of +\u0026thinsp;0.3\u0026deg;C per decade (decreasing through time). Similar warming rates were observed for a shorter period of time (1990\u0026ndash;2010) in cave sections above Windkluft (F\u0026uuml;nferkluft and Wermutsgang). The author critically pointed out that the cave sections for climate studies must be carefully selected in order to avoid misinterpretation that can arise, for example, from an inappropriate selection of observation sites. Filipponi \u003csup\u003e[31]\u003c/sup\u003e proposed criteria suggesting that SH-WK is representative of large parts of the cave.\u003c/p\u003e \u003cp\u003eWind et al. \u003csup\u003e[18]\u003c/sup\u003e addressed long-term temperature changes in Hundsalm Eis- und Tropfsteinh\u0026ouml;hle. Significant warming is apparent in the ice-free part of the cave throughout the period 2008\u0026ndash;2020 (+\u0026thinsp;0.24\u0026deg;C per decade at site HA-JG, which is a principal site in this study). This value agrees with the findings of this study. Ongoing warming is also observed in the higher ice-bearing section of the cave, where trends depend on the distance to the ice (+\u0026thinsp;0.27 for T36 vs. +0.54 for T29; the latter being located closer to the ice).\u003c/p\u003e \u003cp\u003eSome related information is available from neighboring countries. Dom\u0026iacute;nguez-Villar et al. \u003csup\u003e[28]\u003c/sup\u003e observed that temperatures in Postojnia cave (Slovenia) increased at a rate of 0.04\u0026deg;C per decade (2009\u0026ndash;2013). The application of a model indicates that the cave was prone to warming since the 1980\u0026rsquo;s at least and enhanced warming is predicted for the future (+\u0026thinsp;0.15\u0026deg;C per decade). Analyses of historical data since 1935 indeed show that in this cave the average temperature has increased up to 2\u0026deg;C \u003csup\u003e[32]\u003c/sup\u003e, corresponding to a trend of +\u0026thinsp;0.25\u0026deg;C per decade. Significant warming was also reported for another cave in Slovenia (Predjama Cave) as well as for the average temperature across Slovenia (+\u0026thinsp;0.34\u0026deg;C per decade during period 1961 and 2011).\u003c/p\u003e \u003cp\u003eThe calculated temperature trends for the free atmosphere at the elevation of our studied cave sites are about twice as large as those observed inside the caves (+\u0026thinsp;0.5 to +\u0026thinsp;0.7\u0026deg;C per decade). These figures are in agreement with site-specific observations for Patscherkofel \u003csup\u003e[26]\u003c/sup\u003e, Bonistock \u003csup\u003e[27]\u003c/sup\u003e and Hahnenkamm \u003csup\u003e[18]\u003c/sup\u003e. Since a straightforward comparison is hampered by the different elevations and periods, we also placed our results in an extended temporal and regional context provided by regionally gridded climate data i.e., the HISTALP data set \u003csup\u003e[3, 33, 34]\u003c/sup\u003e. These data show the generally strong warming of the Alpine region, which is about twice as large as the Northern Hemispheric average and affected all subregions of the Alps. Enhanced and fairly homogeneous warming occurred since the 1980, with rates matching those from our study (ca. +0.5\u0026deg;C per decade). Our analysis indicates that trends in the low-lying regions are relatively uniform, while mountain stations show some regional distinction, i.e., smaller trends in the Central Alps compared to the Swiss Alps. Referring to the vicinity of Monlesi Ice Cave, Luetscher et al. \u003csup\u003e[35]\u003c/sup\u003e note that the Jura Mountains of northern Switzerland experienced a warming of ca. 0.5\u0026deg;C per decade (1980\u0026ndash;2000). Based on a gridded data set across Switzerland, Ceppi et al. \u003csup\u003e[4]\u003c/sup\u003e found an average warming rate of +\u0026thinsp;0.35\u0026deg;C per decade (1959\u0026ndash;2008). The Obir mountain area in South-Eastern Austria is characterized by outstanding warming rates (ca. +1\u0026deg;C per decade) which, however, is not fully reflected in corresponding HISTALP data.\u003c/p\u003e \u003cp\u003eElevation-dependent and thus also regional warming in the Alps is controversially discussed in literature. Recent reviews conclude that a general picture is difficult to draw due to the variety of influencing factors and processes being strongly related to small-scale topographic effects \u003csup\u003e[36, 2, 34]\u003c/sup\u003e. It is therefore also difficult to interpret the result that trends are significantly larger during the open periods (winter) compared to the closed periods (ca. +0.6 vs. +0.3\u0026deg;C per decade). A thorough investigation of local characteristics (such as cold air pooling, mesoscale wind systems, snow cover, radiation, cloudiness) that influence outside temperature is beyond the scope of this work. Using data from more than 1000 mountain sites across the globe, Pepin and Lundquist \u003csup\u003e[37]\u003c/sup\u003e showed that high-elevation temperature trends are highest near the annual 0\u0026deg;C isotherm due to snow-ice feedback, which was also observed in Switzerland \u003csup\u003e[4]\u003c/sup\u003e.\u003c/p\u003e\n\u003ch3\u003eContext for cave-specific conditions and processes\u003c/h3\u003e\n\u003cp\u003eA comparison of the average cave temperatures with those outside reveals a small difference and a clear dependence on elevation (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Both suggest that the cave temperatures overall correspond to the external conditions, as is corroborated by an independent analysis considering more than 60 caves across the Austrian Alps \u003csup\u003e[26]\u003c/sup\u003e. This result is also in agreement with theoretical considerations \u003csup\u003e[38, 25]\u003c/sup\u003e, which show that cave temperatures are characterized by a damped and delayed response to external forcings, both increasing with distance from the entrance(s). Observations from caves across the world support the related hypothesis that temperature in deeper cave zones asymptotically approaches an equilibrium rock temperature. Following Badino \u003csup\u003e[39]\u003c/sup\u003e, the latter roughly corresponds to the long-term outside temperature. Short-term fluctuations are effectively filtered out due to the large thermal capacity of the rock in which the caves formed. Application of this conceptual model to karstic rock with a thickness of 100 m yields equilibration time scales (i.e., the time to attain thermal equilibrium) in the order of 5\u0026ndash;50 decades, depending on the relative importance of energy exchange involving flow of air and water. Since our studied caves have a rock overburden of less than 100 m it can be argued that the analyzed data were obtained at sites where the seasonal variability is strongly damped, while the climatically relevant variability is still captured, albeit with a delay. These theoretical considerations support the observation that warming trends inside caves are smaller that outside (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). It also helps to understand the nature of the time series mainly considered in this study, which essentially show long-term changes (trends) superimposed by weak seasonal variability (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Systematically delayed responses of cave air temperatures to external forcings cause thermal imbalances, which are reflected in the differences between mean temperatures inside and outside of the caves, as our results also show (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe thermal cyclicity observed in cave records has often been used to classify caves and zones therein. A basic distinction concerns the heterothermic and homothermic zones \u003csup\u003e[40]\u003c/sup\u003e. Considering that the analyzed cave data are characterized by a seasonal variability of the order of a few tenths of a degree (compared to an external variability of about +/-15\u0026deg;C), they cannot represent truly homothermic conditions, but remote sections of heterothermic cave sections. Sedaghatkish et al. \u003csup\u003e[41]\u003c/sup\u003e used a numerical model to simulate heat transfer in a karst conduit due to air flow and thermal conduction, which suggests a more physically based distinction of a \u0026ldquo;diffusive region\u0026rdquo; (where external forcings mainly propagate through thermal conduction in the rock), a \u0026ldquo;convective region\u0026rdquo; (where heat is mainly transported by heat transfer in air) and a \u0026ldquo;deep region\u0026rdquo; (where temperatures are quasi-constant). This distinction basically matches the observationally based approach by Cropley \u003csup\u003e[42]\u003c/sup\u003e and suggests that three of the principal sites being considered in this work can be assigned to the second class (SP, SH and RS; i.e., heterothermic), while HA-JG is similar to the third one (which essentially conforms to the homothermal regime). Medina et al. \u003csup\u003e[21]\u003c/sup\u003e analyzed records from twelve caves across the globe and found three patterns which represent caves sites retracing external variations with little modification except for reduced amplitudes, sites featuring low correlation to surface temperature and a slight thermal delay, and cave sections showing extreme delays or even no relationship to surface variations. Our cave data match the first category, which the authors expect to represent caves with multiple entrances, where air flow plays a dominant role in controlling the cave microclimate (supported by our observations).\u003c/p\u003e \u003cp\u003eSeasonally changing airflow patterns are known for three of the investigated caves, i.e. SP, SH and RS \u003csup\u003e[26, 27, 31]\u003c/sup\u003e, which have two or more entrances at different elevations (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). This, as well as the seasonally changing signs of the thermal gradient (inside vs. outside) are strong indications that processes related to advection and/or convection are important for the exchange of matter and energy with the outside. The so-called chimney effect has been recognized as a key process in this context, which is driven by seasonally changing density differences between the subsurface air and the outside air \u003csup\u003e[38, 43]\u003c/sup\u003e and can be significantly modified by synoptically induced pressure differences \u003csup\u003e[44]\u003c/sup\u003e. The effect is known to be important in Spannagelh\u0026ouml;hle, where it induces a positive deviation of the cave temperature in sections close to the upper entrances compared to outside (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e; Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e) \u003csup\u003e[26]\u003c/sup\u003e. This is most evident during winter when warm air is forced towards the upper entrance to which sites SP-PL and SP-GH are relatively close. The direction of the air flow switches at a threshold temperature of about 2\u0026deg;C \u003csup\u003e[45]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eTr\u0026uuml;ssel \u003csup\u003e[27]\u003c/sup\u003e showed that air flow in Schrattenh\u0026ouml;hle is also driven by the chimney effect, but it is interesting to see that the temperature at the investigated site (SH-WK) is lower than outside. This can be related to the position of the investigated site closer to lower entrances, which in the context of the chimney effect are prone to be colder than regions closer to the upper entrances. However, air flow in SH is generally more complex due to the many entrances, which favor differential air flow in different branches of the cave. We note that wind-driven pressure differences can additionally complicate cave ventilation \u003csup\u003e[44]\u003c/sup\u003e, which could play a role in Schrattenh\u0026ouml;hle, too.\u003c/p\u003e \u003cp\u003eThe RS is characterized by relatively high temperatures compared to outside. Also this cave is known for air flow being driven by the chimney effect which, however, is relatively weak due to the remoteness of this site relative to the main pathway of the air (Supplementary Fig. S3). Still, the direction of air flow seasonally reverses at an outside temperature threshold of 5.5\u0026deg;C \u003csup\u003e[22]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eHundsalm Eis- und Tropfsteinh\u0026ouml;hle is exceptional in several respects. First, it is the smallest cave and its sag-type geometry effectively traps cold air during winter, which enters the cave by density-driven convection. During summer, the air pools without significant exchange with the surface due to a stable thermal stratification, which favors the existence of perennial ice in the upper cave section \u003csup\u003e[18]\u003c/sup\u003e. The cave is therefore classified as \u0026ldquo;static with ice\u0026rdquo; \u003csup\u003e[46]\u003c/sup\u003e. These characteristics give rise to an asymmetric seasonal temperature cycle, which is not observed in the other caves. However, the mainly considered data (HA-JG) are much less affected by these processes, since they were collected in a relatively isolated chamber in the deepest section of the cave and therefore represent conditions corresponding to a homothermic regime.\u003c/p\u003e \u003cp\u003eHaving thus addressed some individual characteristics of the studied caves, it is even more remarkable that the magnitudes of the calculated temperature trends are comparable, i.e., ranging between +\u0026thinsp;0.14 to +\u0026thinsp;0.24\u0026deg;C per decade for the principally investigated sites. It should also be remembered that these sites are located in different regions across the Alps and at different elevations. The rather uniform warming rates can be explained by the fact that records with comparable seasonal variability were considered, which in principle can lead to similar long-term external signals.\u003c/p\u003e \u003cp\u003eThe microclimates of different cave sections, however, retain some spatial variability due to individual physical characteristics such as the relative distance to entrances which influences related processes (e.g., heat conduction, ventilation and flow of water). Our case studies show a correspondingly different seasonality, which is reflected by a progressive attenuation (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e and \u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e), which was already discussed. As far as the trends in the different cave sections are concerned, the results indicate that the decadal warming was stronger near the entrance(s) than in more remote sections (SP-GH vs. SP-TG and HA-29 vs. HA-JG). This is also in agreement with theoretical concepts that predict that long-term, i.e. low-frequency, external variability is delayed with increasing distance from the entrances due to the large thermal capacity of the rock \u003csup\u003e[39]\u003c/sup\u003e.\u003c/p\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eIndirect evidence of warming in the investigated caves\u003c/h2\u003e \u003cp\u003eTemperature strongly controls the formation and existence of ice, which is present in two of the investigated caves. Hundsalm Eis- und Tropfsteinh\u0026ouml;hle hosts perennial firn and ice in its upper section, which reaches a maximum thickness of about 5 m. The ice formed by metamorphism of snow falling though one of the entrances, as well as by refreezing of seepage water \u003csup\u003e[47, 18]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eOur results confirm a development towards increasingly negative mass balances during the recent decade (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e), in parallel with the concurrent atmospheric warming in the Alps. Wind et al. \u003csup\u003e[18]\u003c/sup\u003e developed a model, which allows to calculate annual melt rates as a function of the external temperature. Given the current warming rates, the model predicts that within the next decade, this cave will no longer preserve perennial ice (only seasonal ice). Note, however, that the mass balance of an ice body not only includes surface melt (ablation) but also accumulation and basal melting, both of which also depend on temperature.\u003c/p\u003e \u003cp\u003eA small cave next to Spannagel cave (Spannagel Eish\u0026ouml;hle, SP-EH) contained perennial ice until about 2002 and has been ice-free since then. This adds evidence that caves are sensitive indicators of regional climate change, similar to snow cover and glaciers outside. Significant reduction of ice in caves is a wide-spread phenomenon across Europe \u003csup\u003e[40, 48, 17, 49, 50, 18]\u003c/sup\u003e and there is no indication that this trend will not continue (if not accelerate).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003ePotential influences of anthropogenic cave manipulations\u003c/h2\u003e \u003cp\u003eCertain sections in the investigated caves are open to guided tours, which necessitates some related discussion. Touristic use in Spannagelh\u0026ouml;hle began in 1994, involving ca. 500 m of galleries behind the upper entrance (open all year round). HA was opened as a show cave in 1967, and the guided tours cover most of the upper, i.e., ice-bearing part (open from May until October). The lower section (in which site HA-JG is located) was discovered in 1984 and is not accessible to visitors (locked). Rasslsystem (RS) is part of the Obir Cave system. Guided tours are only offered for part of this system that are far away from RS (since 1991; open April to October). Schrattenh\u0026ouml;hle is also can also be visited within guided tours, but this s only possible since the year 2000 i.e., after the investigation period and in cave sections far off site SH-WK.\u003c/p\u003e \u003cp\u003eTouristic use can have severe impacts on the thermal regime of caves \u003csup\u003e[51]\u003c/sup\u003e. The main influencing factors are the heat emitted by people and the lighting. Evidence is provided by comparing historical temperature data with recent ones and by comparing cave sections or neighboring caves with and without touristic influence \u003csup\u003e[32]\u003c/sup\u003e. Microclimatic studies usually show increased daily variability which related to the number of visitors or concurrent changes in, e.g., humidity and CO\u003csub\u003e2\u003c/sub\u003e concentrations \u003csup\u003e[52, 53]\u003c/sup\u003e. In general, anthropogenic influences are difficult to disentangle from natural changes, as they are cave-specific and less pronounced in well-ventilated caves or sections.\u003c/p\u003e \u003cp\u003eRegarding the potential impact on the investigated caves, it should first be noted that in all three caves the main study sites (SP-PL, HA-JG and RS-PH) are not part of guided tours, i.e. the sites are relatively far from the direct impact of visitors. In terms of timing, tourist activities are rather confined to summer and autumn, which means that excess heat in outer cave sections (where guided tours are offered) is likely to be pushed out through the entrances. This may be relevant for Spannagel and Rasslsystem, but not for Hundsalm Eis- und Tropfsteinh\u0026ouml;hle. This sag-type cave is special because its microclimatic summer regime is characterized by stable air stratification. In principle, this has meant that the excess energy input from the tourists has slowly warmed the cave, which in turn will likely have been compensated for by the increased melting of snow and ice (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). Wind et al. \u003csup\u003e[18]\u003c/sup\u003e noticed small transient disturbances of the cave air temperature being related to guided tours, but these were equilibrated shortly after. Potential long-term effects have not been observed. Regarding the principally considered site in this cave (HA-JG, i.e. the lower, ice-free part), warming by tourists cannot play a role, since it is located in a thermally isolated chamber underneath the ice-bearing part (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). Finally, as far as the decadal trends are concerned, it should be noted that touristic use in all these caves began before the main period under consideration (2000\u0026ndash;2020). It can therefore be assumed that the trends analyzed are not distorted by tourist activities. Data from Schrattenh\u0026ouml;hle are not at all affected by such effects because guided tours started after the investigation period.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eConsequences of warming caves\u003c/h2\u003e \u003cp\u003eThe demise of cave ice as a consequence of warming of cave environments was already been addressed in the context of the Hundsalm Eis- und Tropfsteinh\u0026ouml;hle and the (former) Spannagel Eish\u0026ouml;hle (Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e and \u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). Related processes have been elucidated by e.g. Luetscher et al. \u003csup\u003e[40]\u003c/sup\u003e or Obleitner and Sp\u0026ouml;tl \u003csup\u003e[54]\u003c/sup\u003e. These studies also showed the importance of refreezing seepage water, whose availability strongly depends on the permeability of the host rock, which at high-elevation Alpine caves may also contains sporadic permafrost. Systematic warming is therefore expected to significantly change the availability (amount and timing) of water in these caves, which in turn will also affect the thermal environment \u003csup\u003e[39]\u003c/sup\u003e. The degradation of sporadic permafrost can also lead to the destabilization of cave passages, which happened in Hundsalm cave in 2020 (rockfall). Warming can affect the ventilation regime of high-elevation caves by changing the upper entrances (which otherwise may be plugged by snow or ice). Such effects have already been observed, e.g., in Schrattenh\u0026ouml;hle \u003csup\u003e[27]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eFinally, warming of subterranean environments will potentially have a large impact on biota. For example, Mammola et al. \u003csup\u003e[20]\u003c/sup\u003e argued that caves are largely unexplored habitats for species (bats, spiders, bacteria) that are likely to be sensitive to even small changes in the microclimate \u003csup\u003e[55]\u003c/sup\u003e.\u003c/p\u003e \u003c/div\u003e"},{"header":"Conclusions","content":"\u003cp\u003eStatistically significant changes of air temperature were observed in four caves across the Alps. Our analysis covers the recent two decades and different cave types (static vs. dynamic and also with ice) as well as different sections therein. Key results, interpretations and conclusions are summarized as follows:\u003c/p\u003e\n\u003cul type=\"disc\"\u003e\n \u003cli\u003eMean air temperatures inside these caves are close to those outside and are characterized by much reduced seasonal variability. Caves are systematically warmer in winter (colder in summer) and quasi-vertical temperature gradients in caves are smaller than outside.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eTrends of cave temperature increase are about half of that outside (0.2 vs. 0.5 °C per decade).\u003c/li\u003e\n \u003cli\u003eThese warming rates are characterized by a remarkable linearity and a pronounced spatial and temporal consistency compared to external conditions.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eThese figures represent cave sections characterized by an almost homothermic regime, although there are indications of increased warming rates near the entrances. This and other observations (reduced and lagged seasonal variability) are in line with basic theoretical expectations.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eThe consistency of these trends is remarkable given the variable geometry of the investigated caves, and the results are considered representative of cave environments in the Alps. The synthesis of the data confirms that global warming has reached the interior of these cave systems. \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eIt is not expected that the touristic use of some of these Alpine caves has significantly influenced the results.\u003c/li\u003e\n \u003cli\u003eContinued warming may cause feedbacks and broader changes in cave environments, e.g. changes in the ventilation regime and increased ice melt, which may also affect water resources and subsurface ecosystems.\u003c/li\u003e\n \u003cli\u003eThis study shall stimulate sustained long-term monitoring of temperature and other physical or chemical parameters, as well as related modelling, to better understand the current and future responses of caves in a warming world.\u003c/li\u003e\n\u003c/ul\u003e"},{"header":"METHODS","content":"\u003cp\u003eCave air temperatures were logged using different instruments. Most data from the Austrian caves were obtained using HOBO TempPro v2 logger (Onset Computer Corp.; resolution 0.02 °C, accuracy ±0.21 °C); the first few years of some of these time series were recorded by an earlier model, HOBO Optic StowAway Temp (Onset Computer Corp.; resolution 0.16 °C, accuracy ±0.2 °C). In one cave, Hundsalm Eis- und Tropfsteinhöhle, we also logged the rock temperature in a borehole drilled horizontally into the cave wall (diameter 2.5 cm) using a HOBO TempPro v2 logger with an external sensor (Onset Computer Corp.; resolution 0.02 °C, accuracy ±0.21 °C). The calibration of both these logger types was checked using an ice-water bath in a dewar prior to the installation in the cave and after the loggers were brought back to the laboratory. No drift was observed in any of the time series. The final data were corrected for the offset at 0°C, which for most HOBO TempPro v2 logger was less than 0.1 °C. The logging interval in the Austrian caves varied between 1 and 3 hours.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTemperature data in the Schrattenhöhle were collected using a different logger model, i.e., Hamster and Hamster-A (Elpro; resolution 0.05 °C, accuracy ±0.2 °C) and the calibration provided by the manufacturer was used. Occasional cross-checks using other logger types suggest an overall accuracy of ±0.1 °C. The logging interval in Schrattenhöhle was 1.5 hours (1990-1999) and 1 hour (since 2000). \u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThese instruments proved reliable to recover cave temperatures in a methodically consistent manner and over long periods of time. Minor gaps in the records were mostly due maintenance (downloading the data or change of sensors) which in caves can need longer time. Major gaps evolved in the context of post-processing raw data focused on detection of physically unplausible data and comparison to parallel series as far as available. It is helpful in this context, that cave temperatures are less variable compared to outside temperatures. The reason(s) for physically unplausible data can be related to various environmental issues including, e.g. intermittent influences by water, rime or ice on sensors. About 15% of the original data from Hundsalm and Spannagel caves were finally excluded from further analysis, significantly more from SH and RS (62 and 42%). The postprocessed data were then merged to daily averages, which form the basis of this study.\u003c/p\u003e\n\u003cp\u003eThe mostly complex geometrical structure of caves and related impacts on their microclimate \u0026nbsp; questions the spatial representativeness of analysis referring to a certain cave section. Practically, the siting of a sensor is often determined by on-site compromises considering scientific arguments, accessibility and constraints to safely mount sensors, which usually can hardly be visited several times a year on the other hand. Related uncertainties are difficult to quantify and would require an extensive measurement network or detailed numerical simulations of air flow along the investigated caves (both being beyond the scope of this work). Nevertheless, in order to address the topic of this investigation (temperature trends) we also analyzed observations in neighboring sections of the studied caves to study the spatial representativity of the best possible chosen reference time series. \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u003c/p\u003e\n\u003cp\u003eDaily values of air temperature measured at automatic weather stations (AWS) in the vicinity of the four caves were used to calculate free-atmosphere temperatures at the elevation of the cave sites. The latter were derived by linear interpolation based on daily lapse rates calculated from station pairs located at the valley bottom and at crest height in representative valleys. Correspondingly, we used the station pair JE (530 m a.s.l.) and PA (2251 m a.s.l.) to derive the lapse rate for the SP and HA environment, while BA (623 m a.s.l.) and SE (1040 m a.s.l.) are used for RS and EN (1000 m a.s.l.) and PI (2100 m a.s.l.) for SR, respectively. The calculated temperatures were validated using additional AWS data (Plattkopf (2260 m a.s.l.), HA-entrance (1520 m a.s.l.), Bonistock (2160 m a.s.l.) and Zell-Pfarre (900 m a.s.l.; Fig. 1) and corresponding offsets were applied to refine the calculated temperatures at the elevation of the cave sites. Finally, the ice development in HA cave was monitored using stakes drilled into the ice yielding annual measures of surface height changes with an estimated accuracy of ±1 cm\u0026nbsp;\u003csup\u003e[18]\u003c/sup\u003e.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eCharacterization of the typical meteorological conditions in the studied caves and their outside environment is based on standard statistics for a reference period (2008-2010) with most consistent data. Shorter and longer periods are considered in the context of case studies addressing specific problems, e.g. the spatial representativity. Trend analysis employs the calculation of linear trends and Mann-Kendall tests \u003csup\u003e[56, 57]\u003c/sup\u003e. The former is firstly based on ordinary least squares regression of daily values covering the period 2000-2020 and consideration of respective slopes. Since this method can suffer from its sensitivity to outliers, we also employed the Theil–Sen trend estimator\u0026nbsp;\u003csup\u003e[58, 59]\u003c/sup\u003e. This is a comparatively robust method to fit a line to sample points by considering the median of slopes determined by all pairs of successive\u0026nbsp;data.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eCalculations were performed using Matlab®\u0026nbsp;\u003csup\u003e[60]\u003c/sup\u003e and were also applied to seasonal subsets covering the so-called open and closed periods. Following Wind et al.\u0026nbsp;\u003csup\u003e[18]\u003c/sup\u003e, the latter comprise the months November to April and May until October, respectively.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLong-term monitoring in the Austrian caves was partially supported by the Austrian Science Fund (FWF) grants Y122-GEO, P18637-N10, P20618-N10 and P31874-N37 to CS. Data from automatic weather stations were provided by national sevices i.e., Geosphere Austria, Slovenian Environment Agency, Hydrographischer Dienst Tirol, Hydrographischer Dienst K\u0026auml;rnten and \u0026nbsp; VERBUND AG, Austria.\u003c/p\u003e\n\u003cp\u003eFor open access purposes, the authors have applied a CC BY public copyright license to any author accepted manuscript version arising from this submission. Funding for research in the Swiss caves was partially provided by the\u0026nbsp;Foundation Naturerbe Karst und H\u0026ouml;hlen Obwalden (NeKO), Alpnach, Switzerland.\u0026nbsp;Several members of caving associations provided access to caves and on-site support.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe idea of the manuscript was developed by CS and MT, who also provided the cave data. Main analysis was\u0026nbsp;provided by FO, who also drafted the manuscript, tables and figures except of Fig. 1. Figures in the supplement were contributed \u0026nbsp;by CS, where Fig. 4 is adapted from\u0026nbsp;Filipponi (2003).\u0026nbsp;The manuscript was reviewed by all authors, thereby also considering comments by\u0026nbsp;Beat Niederberger. The\u0026nbsp;authors declare no competing interests and underlying data can be provided on request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLong-term monitoring in the Austrian caves was partially supported by the Austrian Science Fund (FWF) grants Y122-GEO, P18637-N10, P20618-N10 and P31874-N37 to CS.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eData being used in this study are available upon request. Please contact [email protected] or [email protected].\u003c/p\u003e"},{"header":"References","content":"\u003cp\u003eGulev, S.K. et al.(Changing State of the Climate System. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (ed. by Masson-Delmotte, V. et al.). Cambridge University Press, Cambridge, 287-422. https://doi.org/10.1017/9781009157896.004 (2021).\u003c/p\u003e\n\u003col\u003e\n \u003cli\u003eMountain Research Initiative EDW Working Group. Elevation-dependent warming in mountain regions of the world. \u003cem\u003eNature Climate Change\u003c/em\u003e\u003cstrong\u003e5\u003c/strong\u003e, 424-430 (2015).\u003c/li\u003e\n \u003cli\u003eAuer I., R. et al.. HISTALP \u0026ndash; historical instrumental climatological surface time series of the Greater Alpine Region. \u003cem\u003eInt. J. Climatol.\u003c/em\u003e\u003cstrong\u003e27\u003c/strong\u003e, 17-46 (2007).\u003c/li\u003e\n \u003cli\u003eCeppi, P. et al.. Revisiting Swiss temperature trends 1959-2008. \u003cem\u003eInt. J. Climatol.\u003c/em\u003e\u003cstrong\u003e32\u003c/strong\u003e, 203-213. https://doi.org/10.1002/joc.2260 (2012).\u003c/li\u003e\n \u003cli\u003eHock, R.G.et al. . High Mountain Areas. In: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (ed. by H.-O. P\u0026ouml;rtner, H.-O. et al.). Cambridge University Press, Cambridge, 131-202. https://doi.org/10.1017/9781009157964.004 (2019).\u003c/li\u003e\n \u003cli\u003eSommer, C. et al.. Rapid glacier retreat and downwasting throughout the European Alps in the early 21st century. \u003cem\u003eNature Commun.\u003c/em\u003e\u003cstrong\u003e11\u003c/strong\u003e:3209. https://doi.org/10.1038/s41467-020-16818-0 (2020).\u003c/li\u003e\n \u003cli\u003eŽebre, M. et al.. 200 years of equilibrium‑line altitude variability across the European Alps (1901\u0026minus;2100). \u003cem\u003eClimate Dyn.\u003c/em\u003e\u003cstrong\u003e56\u003c/strong\u003e, 1183-1201 (2021).\u003c/li\u003e\n \u003cli\u003eBeniston, M. et al.(2018). The European mountain cryosphere: a review of its current state, trends, and future challenges. \u003cem\u003eThe Cryosphere\u003c/em\u003e\u003cstrong\u003e12\u003c/strong\u003e, 759-794 (2018).\u003c/li\u003e\n \u003cli\u003eHaberkorn, A., Kenner, R., Noetzli, J. \u0026amp; Phillips, M. Changes in ground temperature and dynamics in mountain permafrost in the Swiss Alps. \u003cem\u003eFrontiers in Earth Science\u003c/em\u003e\u003cstrong\u003e9\u003c/strong\u003e:626686, https://doi.org/10.3389/feart.2021.626686 (2021).\u003c/li\u003e\n \u003cli\u003eWebb, B.W. \u0026amp; Nobilis, F. Long term water temperature trends in Austrian rivers. \u003cem\u003eHydrol. Sci. J.\u003c/em\u003e\u003cstrong\u003e40\u003c/strong\u003e, 83-96 (1995).\u003c/li\u003e\n \u003cli\u003eNiedrist, G.H. Substantial warming of Central European mountain rivers under climate change. \u003cem\u003eRegional Environ. Change\u003c/em\u003e\u003cstrong\u003e23\u003c/strong\u003e:43 https://doi.org/10.1007/s10113-023-02037-y (2023).\u003c/li\u003e\n \u003cli\u003eKormann, C., Francke, T. \u0026amp; Bronstert, A. Detection of regional climate change effects on alpine hydrology by daily resolution trend analysis in Tyrol, Austria. \u003cem\u003eJ. Water Climate Change\u003c/em\u003e\u003cstrong\u003e6\u003c/strong\u003e (1): 124-143. https://doi.org/10.2166/wcc.2014.099 (2015).\u003c/li\u003e\n \u003cli\u003eBenz, S.A., Bayer, P., Winkler, G. \u0026amp; Blum, P. Recent trends of groundwater temperatures in Austria. Hydrol. \u003cem\u003eEarth Syst. Sci.\u003c/em\u003e\u003cstrong\u003e22\u003c/strong\u003e, 3143-3154 (2018).\u003c/li\u003e\n \u003cli\u003eHemmerle, H. \u0026amp; Bayer, P. Climate change yields groundwater warming in Bavaria, Germany. \u003cem\u003eFrontiers in Earth Sci.\u003c/em\u003e\u003cstrong\u003e8\u003c/strong\u003e: 575894. doi: 10.3389/feart.2020.575894 (2020).\u003c/li\u003e\n \u003cli\u003eGilbert, A. \u0026amp; Vincent, C. Atmospheric temperature changes over the 20th century at very high elevations in the European Alps from englacial temperatures. \u003cem\u003eGeophys. Res. Lett.\u003c/em\u003e \u003cstrong\u003e40\u003c/strong\u003e, 2102. doi:10.1002/grl.50401(2023).\u003c/li\u003e\n \u003cli\u003ePerrier, F., Le Mou\u0026euml;l, J. L., Poirier, J. \u0026amp; Shnirman, M.G. Long-term climate change and surface versus underground temperature measurements in Paris. \u003cem\u003eInt. J. Climatol.\u0026nbsp;\u003c/em\u003e\u003cstrong\u003e25,\u003c/strong\u003e 1619-1631. https://doi.org/10.1002/joc.1211(2005).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C., Wimmer, M., Pavuza, R. \u0026amp; Plan, L. Ice caves in Austria. in: Ice Caves (Perşoiu, A. \u0026amp; Lauritzen, S.E.), 237-262, Amsterdam, Elsevier (2018).\u003c/li\u003e\n \u003cli\u003eWind, M., Obleitner, F., Racine, T. \u0026amp; Sp\u0026ouml;tl, C. Multi-annual temperature evolution and implications for cave ice development in a sag-type ice cave in the Austrian Alps. \u003cem\u003eThe Cryosphere\u003c/em\u003e\u003cstrong\u003e16\u003c/strong\u003e, 3163-3179. https://doi.org/10.5194/tc-16-3163-2022 (2022).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C., Egger, M., Mangini, A. \u0026amp; Pavuza, R. Das Karstsystem am Obstanser See (Kartitsch, Osttirol) \u0026ndash; geowissenschaftliche Ergebnisse. \u003cem\u003eDie H\u0026ouml;hle\u003c/em\u003e\u003cstrong\u003e68\u003c/strong\u003e, 59-78 (2017).\u003c/li\u003e\n \u003cli\u003eMammola, S. et al.. Climate change going deep: The effects of global climatic alterations on cave ecosystems, \u003cem\u003eAnthropocene Rev.\u003c/em\u003e\u003cstrong\u003e6\u003c/strong\u003e, 98-116 (2019).\u003c/li\u003e\n \u003cli\u003eMedina, M.J. et al.. Temperature variation in caves and its significance for subterranean ecosystems. \u003cem\u003eSci. Rep.\u003c/em\u003e\u003cstrong\u003e13\u003c/strong\u003e: 20735. https://doi.org/10.1038/s41598-023-48014-7 (2023).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C., Fairchild, I.J. \u0026amp; Tooth, A.F. Cave air control on dripwater geochemistry, Obir Caves (Austria): Implications for speleothem deposition in dynamically ventilated caves. \u003cem\u003eGeochim.\u0026nbsp;\u003c/em\u003e\u003cem\u003eCosmochim. Acta\u003c/em\u003e\u003cstrong\u003e69\u003c/strong\u003e, 2451-2468 (2005).\u003c/li\u003e\n \u003cli\u003eBorreguero, M. \u0026amp; Pahud, A. Climatologie de la grotte des Pingouins. In: Lapi di Bou: recherches et explorations sp\u0026eacute;l\u0026eacute;ologiques,1974-2009, 41-54 (2009).\u003c/li\u003e\n \u003cli\u003eWimmer, M. 25 Jahre Eis- und H\u0026ouml;hlenklimaforschung im Sch\u0026ouml;nberg-H\u0026ouml;hlensystem. \u003cem\u003eMitt. Landesvereins. H\u0026ouml;hlenk. Ober\u0026ouml;sterreich\u003c/em\u003e\u003cstrong\u003e65/124\u003c/strong\u003e, 17-25 (2019).\u003c/li\u003e\n \u003cli\u003eWigley, T.M.L. \u0026amp; Brown, M.C. The physics of caves. In: The Science of Speleology (ed. by Ford, T.D. and Cullingford, C.H.D.), 329-358, Academic Press (1976).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C. \u0026amp; Pavuza, R. H\u0026ouml;hlenatmosph\u0026auml;re. in: H\u0026ouml;hlen und Karst in \u0026Ouml;sterreich (ed. by C. Sp\u0026ouml;tl, L. Plan, E. Christian), Linz, Ober\u0026ouml;sterreichisches Landesmuseum, 123-138 (2016).\u003c/li\u003e\n \u003cli\u003eTr\u0026uuml;ssel, M. H\u0026ouml;hlenklima: Langzeit-Temperaturmessungen in der Schrattenh\u0026ouml;hle, Melchsee-Frutt (CH-6064 Kerns OW). Akten des 14. Nationalen Kongresses für H\u0026ouml;hlenforschung, 227-232, H\u0026ouml;hlenforschergemeinschaft Region Hohgant (2019).\u003c/li\u003e\n \u003cli\u003eDom\u0026iacute;nguez-Villar, D. et al.Is global warming affecting cave temperatures? Experimental and model data from a paradigmatic case study. \u003cem\u003eClimate Dyn.\u003c/em\u003e\u003cstrong\u003e45\u003c/strong\u003e, 569-581 (2015).\u003c/li\u003e\n \u003cli\u003eColucci, R., Fontana, D., Forte, E. \u0026amp; Potleca, M., Guglielmin, M. Response of ice caves to weather extremes in the southeastern Alps, Europe, \u003cem\u003eGeomorphology\u003c/em\u003e\u003cstrong\u003e261\u003c/strong\u003e, 1-11, http://dx.doi.org/10.1016/j.geomorph.2016.02.017 (2015).\u003c/li\u003e\n \u003cli\u003eRacine, T., Reimer, P. J. \u0026amp; Sp\u0026ouml;tl, C. Multi‐centennial mass balance of perennial ice deposits in Alpine caves mirrors the evolution of glaciers during the Late Holocene. \u003cem\u003eSci. Rep.,\u003c/em\u003e\u003cstrong\u003e12\u003c/strong\u003e: 11374. https://doi.org/10.1038/s41598-022-15516-9 (2022).\u003c/li\u003e\n \u003cli\u003eFilipponi, M. Das Klima der Schrattenh\u0026ouml;hle: Auswertung einer 10-j\u0026auml;hrigen Messreihe, ETH Z\u0026uuml;rich, Institut f\u0026uuml;r Atmosph\u0026auml;re und Klima, https://doi.org/10.3929/ethz-a-004662252 (2003).\u003c/li\u003e\n \u003cli\u003e\u0026Scaron;ebela, S. \u0026amp; Turk, J. Comparison of historical and current temperatures in show caves (Slovenia). \u003cem\u003eSN Appl. Sci.\u003c/em\u003e\u003cstrong\u003e4\u003c/strong\u003e, 1. https://doi.org/10.1007/s42452-021-04881-1(2022).\u003c/li\u003e\n \u003cli\u003eBrunetti, M. et al. Climate variability and change in the Greater Alpine Region over the last two centuries based on multi-variable analysis. \u003cem\u003eInt. J. Climatol.\u003c/em\u003e\u003cstrong\u003e29\u003c/strong\u003e, 2197\u0026ndash;225.(2009).\u003c/li\u003e\n \u003cli\u003eKuhn M. \u0026amp; Olefs M. Elevation-Dependent Climate Change in the European Alps, Oxford Research Encyclopedias, Climate Science, 1-26. https://doi.org/10.1093/acrefore/9780190228620.013.762 (2020)\u003c/li\u003e\n \u003cli\u003eLuetscher, M. \u0026amp; Jeannin, P.Y., Haeberli W. Ice caves as an indicator of winter climate evolution: a case study from the Jura Mountains, \u003cem\u003eThe Holocene\u003c/em\u003e\u003cstrong\u003e15\u003c/strong\u003e, 982-/993. https://doi.org/10.1191/0959683605hl872ra (2005).\u003c/li\u003e\n \u003cli\u003eGobiet, A. et al.21st century climate change in the European Alps - a review. \u003cem\u003eSci. Total Environm.\u003c/em\u003e\u003cstrong\u003e493\u003c/strong\u003e, 1138-1151 (2014).\u003c/li\u003e\n \u003cli\u003ePepin N.C., \u0026amp; Lundquist J.D. Temperature trends at high elevations: Patterns across the globe. \u003cem\u003eGeophysical Research Letters\u003c/em\u003e\u003cstrong\u003e35\u003c/strong\u003e:, L14701. https://doi.org/10.1029/2008GL034026 (2008).\u003c/li\u003e\n \u003cli\u003eCigna, A.A. An analytical study of air circulation in caves.\u003cem\u003e\u0026nbsp;Intern. J. Speleol.\u003c/em\u003e 3, 41-54 (1968).\u003c/li\u003e\n \u003cli\u003eBadino, G. Cave temperatures and global climatic change. \u003cem\u003eInt. J. Speleol.\u003c/em\u003e\u003cstrong\u003e33\u003c/strong\u003e, 103-114 (2004).\u003c/li\u003e\n \u003cli\u003eLuetscher, M., Lismonde, B. \u0026amp; Jeannin, P.Y. Heat exchanges in the heterothermic zone of a karst system: Monlesi cave, Swiss Jura Mountains, \u003cem\u003eJ. Geophys. Res.\u003c/em\u003e\u003cstrong\u003e113\u003c/strong\u003e, F02025, doi:10.1029/2007JF000892 (2008).\u003c/li\u003e\n \u003cli\u003eSedaghatkish, A. et al. Modeling heat transfer for assessing the convection length in ventilated caves. \u003cem\u003eJ. Geophys. Res.: Earth Surface\u003c/em\u003e\u003cstrong\u003e129\u003c/strong\u003e, e2024JF007646. https://doi.org/10.1029/2024JF007646 (2024).\u003c/li\u003e\n \u003cli\u003eCropley J. B. Influence of surface conditions on temperatures in large cave systems. \u003cem\u003eBull. Nat. Speleol. Soc.\u003c/em\u003e\u003cstrong\u003e27\u003c/strong\u003e, 1\u0026ndash;10 (1965).\u003c/li\u003e\n \u003cli\u003eGabrovsek, F. How do caves breathe: The airflow patterns inckarst underground. \u003cem\u003ePLoS ONE\u003c/em\u003e\u003cstrong\u003e18(4)\u003c/strong\u003e: e0283767, https://doi.org/10.1371/journal.pone.0283767 (2023).\u003c/li\u003e\n \u003cli\u003eKukuljan, L., Gabrov\u0026scaron;ek, F. \u0026amp; Covington, M.D. The relative importance of wind-driven and chimney effect cave ventilation: Observations in Postojna Cave (Slovenia). \u003cem\u003eInt. J. Speleol.\u003c/em\u003e\u003cstrong\u003e50\u003c/strong\u003e, 275-288. https://doi.org/10.5038/1827-806X.50.3.2392 (2021).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C. \u0026amp; Pavuza, R. Zur Bewetterung der Spannagelh\u0026ouml;hle (Kat. Nr. 2411/1). \u003cem\u003eH\u0026ouml;hlenk. Mitt. Landesvereins f\u0026uuml;r H\u0026ouml;hlenkunde in Tirol\u003c/em\u003e\u003cstrong\u003e39\u003c/strong\u003e, 5-9 (2001).\u003c/li\u003e\n \u003cli\u003eLuetscher, M. \u0026amp; Jeannin, P. Y. A process-based classification of alpine ice caves. \u003cem\u003eTheor. Appl. Karst.\u003c/em\u003e\u003cstrong\u003e17\u003c/strong\u003e, 5\u0026ndash;10 (2004).\u003c/li\u003e\n \u003cli\u003eSp\u0026ouml;tl, C., Reimer, P.J. \u0026amp; Luetscher, M. Long-term mass balance of perennial firn and ice in an Alpine cave (Austria): Constraints from radiocarbon-dated wood fragments. \u003cem\u003eThe Holocene\u003c/em\u003e\u003cstrong\u003e24\u003c/strong\u003e, 165-175 (2014).\u003c/li\u003e\n \u003cli\u003eKern, Z. \u0026amp; Perșoiu, A. Cave ice - the imminent loss of untapped mid-latitude cryospheric palaeoenvironmental archives. \u003cem\u003eQuat. Sci. Rev.\u003c/em\u003e\u003cstrong\u003e67\u003c/strong\u003e, 1-7 (2013).\u003c/li\u003e\n \u003cli\u003eColucci, R.R. \u0026amp; Guglielmin, M. Climate change and rapid ice melt: Suggestions from abrupt permafrost degradation and ice melting in an alpine ice cave. \u003cem\u003eProgress Phys. Geogr.\u003c/em\u003e\u003cstrong\u003e43\u003c/strong\u003e, 561-573 (2019).\u003c/li\u003e\n \u003cli\u003ePerșoiu, A. et al.. Record summer rains in 2019 led to massive loss of surface and cave ice in SE Europe. \u003cem\u003eThe Cryosphere\u003c/em\u003e\u003cstrong\u003e15\u003c/strong\u003e, 2383-2399 (2021).\u003c/li\u003e\n \u003cli\u003eVillar, E. et al. Ambient temperature variations in the hall of paintings of Altamira cave due to the presence of visitors. \u003cem\u003eTrans. British Cave Res. Ass.\u003c/em\u003e\u003cstrong\u003e11\u003c/strong\u003e, 99-103 (1984).\u003c/li\u003e\n \u003cli\u003ePrelov\u0026scaron;ek, M., \u0026Scaron;ebela, S. \u0026amp; Turk, J. Carbon dioxide in Postojna Cave (Slovenia): spatial distribution, seasonal dynamics and evaluation of plausible sources and sinks. \u003cem\u003eEnviron. Earth. Sci.\u0026nbsp;\u003c/em\u003e\u003cstrong\u003e77\u003c/strong\u003e, 289. https://doi.org/10.1007/s12665-018-7459-6 (2018).\u003c/li\u003e\n \u003cli\u003eFernandez-Cortes, A. et al. Detection of human-induced environmental disturbances in a show cave. \u003cem\u003eEnviron. Sci. Pollut. Res.\u003c/em\u003e\u003cstrong\u003e18\u003c/strong\u003e, 1037\u0026ndash;1045. https://doi./10.1007/s11356-011-0513-5 (2011).\u003c/li\u003e\n \u003cli\u003eObleitner, F. \u0026amp; Sp\u0026ouml;tl, C. The mass and energy balance of ice within the Eisriesenwelt cave, Austria. \u003cem\u003eThe Cryosphere\u003c/em\u003e\u003cstrong\u003e5\u003c/strong\u003e, 245\u0026ndash;257, 2011. https://doi.org/10.5194/tc-5-245-2011 (2011).\u003c/li\u003e\n \u003cli\u003ePerry, R.W. A review of factors affecting cave climates for hibernating bats in temperate North \u003cem\u003eAmerica. Environ. Reviews\u003c/em\u003e\u003cstrong\u003e21\u003c/strong\u003e, 28-39, https://doi.org/10.1139/er-2012-0042 (2012)66.\u003c/li\u003e\n \u003cli\u003eMann, H.B. Non-parametric tests against trend. \u003cem\u003eEconometrica\u003c/em\u003e\u003cstrong\u003e13\u003c/strong\u003e, 163-171 (1945).\u003c/li\u003e\n \u003cli\u003eKendall, M.G. Rank Correlation Methods. 4th ed., Charles Griffin, London (1975).\u003c/li\u003e\n \u003cli\u003eTheil, H. A rank-invariant method of linear and polynomial regression analysis. \u003cem\u003eNederl. Akad. Wetensch. Proc.\u003c/em\u003e\u003cstrong\u003e53\u003c/strong\u003e, 386\u0026ndash;392 (1950).\u003c/li\u003e\n \u003cli\u003eSen, P. K. Estimates of the regression coefficient based on Kendall\u0026apos;s tau. \u003cem\u003eJ. Amer. Stat. Ass.\u003c/em\u003e\u003cstrong\u003e63\u003c/strong\u003e (324), 1379\u0026ndash;1389 (1968).\u003c/li\u003e\n \u003cli\u003eThe MathWorks Inc. MATLAB version: 9.13.0 (R2022b), Natick, Massachusetts: The MathWorks Inc. https://www.mathworks.com (2022).\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Caves, temperature trends, Europe, regional climate change","lastPublishedDoi":"10.21203/rs.3.rs-4801555/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4801555/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eCave air temperatures in four caves in the European Alps show statistically significant warming trends of about 0.2 °C per decade over the last two decades (2000-2020). These trends are about half as large as those observed outside and are characterized by a remarkable spatial and temporal consistency. The investigated caves represent different types in terms of their ventilation regime and one of them also hosts perennial ice. Key observation sites are located in cave sections where the temporal variability of air temperature is strongly attenuated compared to outside conditions and data from different cave sections show that the main results are valid for large parts of the investigated caves.\u003c/p\u003e\n\u003cp\u003eContinued warming will lead to broad changes in Alpine cave environments, including changes in the ventilation regime, karst hydrology and subsurface ecosystems. The observed subsurface warming has a particular strong effect on the long-term preservation of perennial ice present in some of these caves. This is shown for an ice cave in the Austrian Alps, where enhanced melt of ice correlates with the observed warming. This cave (and similar ones) will not be able to hold perennial ice beyond the next decade.\u003c/p\u003e","manuscriptTitle":" Climate warming has arrived in caves of the European Alps","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-02 09:58:52","doi":"10.21203/rs.3.rs-4801555/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-10-03T15:04:03+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-26T07:26:16+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-17T05:50:28+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"157685675965073356972667877836083389781","date":"2024-09-16T07:21:42+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-08-20T16:57:55+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"168128677513096880606789062334680085433","date":"2024-08-06T16:10:22+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"221288058576204865314310882578587256535","date":"2024-08-06T16:01:39+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-08-06T15:58:23+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-08-06T15:53:22+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2024-08-06T15:44:13+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-08-03T06:36:52+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2024-07-25T11:34:49+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"19a26147-8e7d-4b9e-a444-8d1a65f76741","owner":[],"postedDate":"September 2nd, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":36734767,"name":"Earth and environmental sciences/Climate sciences"},{"id":36734768,"name":"Earth and environmental sciences/Environmental sciences"}],"tags":[],"updatedAt":"2024-11-11T16:01:46+00:00","versionOfRecord":{"articleIdentity":"rs-4801555","link":"https://doi.org/10.1038/s41598-024-78658-y","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2024-11-10 15:57:45","publishedOnDateReadable":"November 10th, 2024"},"versionCreatedAt":"2024-09-02 09:58:52","video":"","vorDoi":"10.1038/s41598-024-78658-y","vorDoiUrl":"https://doi.org/10.1038/s41598-024-78658-y","workflowStages":[]},"version":"v1","identity":"rs-4801555","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4801555","identity":"rs-4801555","version":["v1"]},"buildId":"ApUGefWb6u5IBVtyqm6d5","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}