The changes in precipitation and temperature regimes brought on by the current climate change have influenced ecosystems globally. The consequences of climate change on plant phenology have been widely investigated during the last few years. However, the underlying causes of the timing of autumn phenology have not been fully clarified yet. Here, we focused on the onset (10%) of leaf colouring—LCO—(Biologische Bundesanstalt, Bundessortenamt und Chemische Industrie (BBCH) 92) of European beech (Fagus sylvatica, L.) as an important native tree species growing throughout Europe. Studied beech stands are located along the natural distribution range of the European beech in Western Carpathians (Slovakia) at different altitudes from lowlands (300 m a.s.l.) to uplands (1050 m a.s.l.) and climatic regions from warm to cold. To define limiting climate conditions for LCO, we established several bioclimatic indices as indicators of meteorological drought: climatic water balance (CWB), standardized precipitation index (SPI), standardized precipitation-evapotranspiration index (SPEI), dry period index (DPI), and heat waves (HW). In addition, meteorological variables such as monthly mean temperatures and precipitation totals were taken into account. Throughout the 23-year period (1996–2018) of ground-based phenological observations of temperate beech forests, the timing of LCO was significantly delayed (p ≤ 0.05) in the middle to high altitudes, while in the lowest altitude, it remained unchanged. Over the last decade, 2009–2018, LCO in middle altitudes started at comparable to low altitudes and, at several years, even later. This resulted mainly from the significant negative effect of drought prior to this phenological phase (p ≤ 0.01) expressed through a 1-month SPI in September (SPI_IX) at the stand at the low-altitude and warm-climatic region. Our results indicate that the meteorological drought conditioned by lower total precipitation and higher evapotranspirative demands in the warmer climate advance leaf senescence. However, at present time, growth in rising temperature and precipitation is acceptable for most beech stands at middle to high altitudes. Beech utilizes these conditions and postpones the LCO by 0.3–0.5 and 0.6–1.2 day per year at high and middle altitudes, respectively. Although we show the commencing negative effect of drought at mid-altitudes with lower (below 700 mm) total annual precipitation, the trend of LCO in favourable warm climates is still significantly delayed. The ongoing warming trend of summer months suggests further intensification of drought as has started to occur in middle altitudes, spreading from the continual increase of evapotranspiration over the next decades.

In the current alteration of temperature and snow cover regimes, the impacts of winter climate have received considerably less attention than those of the vegetation period. In this study, we present the results demonstrating the influence of the winter climate conditions on the Mountain pine (Pinus mugo Turra) communities in High Tatra Mts (Western Carpathians). The changes in greenness in 2000–2020 were represented by the inter-annual differences of satellite-derived Normalized Difference Vegetation Index (NDVI). The winter climate conditions were characterized by climate indices calculated from the temperature and snow cover data measured at Skalnaté Pleso Observatory (1778 m a.s.l.) over the period between 1941–2020. Areas with P. mugo were classified into two density classes and five altitudinal zones of occurrence. The partial correlation analyses, which controlled the influence of summer climate, indicated that winter warm spells (WWS) caused a significant decrease in the greenness of the P. mugo thickets growing in the dense class D2 (R = −0.47) and in the altitudinal zones A2 (1600–1700 m a.s.l.) and A3 (1700–1800 m a.s.l.) with R = −0.54 for each zone. The changes in greenness were related to the average snow depth (ASD) as well, particularly in the dense class D2 (R = 0.45) and in the altitudinal zone A2 (R = 0.50). Here, in the summers following winters with the incidence of WWS or low ASD, we found decreased greenness following the injury of P. mugo shrubs, but NDVI after winters with higher ASD indicated more greenness. At lower altitudes, injuries may result in the loss of competition capacity of P. mugo near the timberline, where taller mountain tree species can utilize the conditions of warmer climate for expansion. We also found a significant positive effect of warmer winter seasons in the sparse P. mugo thickets (D1) with R = 0.50 and at higher altitudes (R = 0.49 in A4—1800–1900 m a.s.l.; R = 0.53 in A5—1900–2000 m a.s.l.). The increased temperatures in December correlated significantly with the increase of the greenness in all P. mugo pixels (R = 0.47), with the most pronounced effect in the sparse class D1 (R = 0.57) and in altitudinal zones A4 (R = 0.63) and A5 (R = 0.44), creating advantageous conditions for the thermophilisation of the alpine zone by P. mugo.

Variability of precipitation and runoff in the entire high Tatra Mountains in the period 1961–2010

This study is focused on the research of selected Pinus species exposed to high ozone concentrations in the mountain environment. We noticed different values of modelled ozone doses (MOD) up-taken by Mountain pine ( Pinus mugo Turra) in the High Tatra Mts (SK–HTMts) and Swiss stone pine ( Pinus cembra L.) in the Alpes-Mercantour (FR–AlpMar) during the growing season 2019. The MOD values were obtained by multiplicative DO3SE model, while we also tested a new approach based on modification of input ozone data. The MOD values were obtained by multiplicative DO3SE model, while we also tested a new approach based on modification of input ozone data. Testing has shown that ozone input based on passive sampling may be used in MOD modelling for sites situated in the subalpine zone where the operation of active monitors is limited.. Presented results confirmed the assumption regarding stomatal ozone flux reduction due to the occurrence of soil drought in hot and dry summer weather typical for the Mediterranean climate region. Despite the limitation of stomatal flux, foliar ozone specific injury on two years needles of P. cembra was substantially higher in comparison to the incidence of ozone injury symptoms observed on two years needles of P. mugo in SK–HTMts. It may suggest low phytotoxicity of given MOD or efficient resistance of P. mugo against oxidative stress. In addition, the visible injury index (VINX) covering the broad effect of biotic and abiotic harmful agents was appraised on P. mugo . Percentage of affected surface indicated moderate deterioration of needle injury at the end of the growing season, particularly due to traces of mechanical damage.

The high ambient ozone concentrations cause impairing effects on vegetation leading to plant injuries. The potential ozone uptake to vegetation through open stomata can be quantified using stomatal conductance measurements under the local environmental conditions. This study compares the ozone stomatal conductance to vegetation obtained with a modified Jarvis formula adopted from the Vegetation Manual of United Nations Economic Commission for Europe, and experimental field measurements’ data. The stomatal conductance was measured by a portable photosynthesis and gas exchange analyzer system LiCOR6400. The measurements were performed in the submontane environment of the High Tatra Mountains in Slovakia on Swiss pine (Pinus cembra), as a native species of the local flora. According to previous studies, Swiss pine is considered as an ozone-sensitive species. The modified Jarvis model for the ozone stomatal conductance is compared with the field measurements. The suitable parameterization of the modified Jarvis model for Swiss pine is obtained. The parameterization of stomatal conductance for Swiss pine in the local environment would help understand its specificity and similarity to other conifer species. In the case of using parameterization for a boreal coniferous from the Vegetation Manual of the International Cooperative Programme on Effects of Air Pollution on Natural Vegetation and Crops, validation of the model with the measurements without temperature adjustment of the conifer chamber achieved a coefficient of determination of $R^2=0.75$. This result is not in contradiction with the previous researches. With the optimal set of parameters, obtained in this paper, the Jarvis model reaches $R^2=0.85$. The data suggest that Jarvis-type models with appropriate parameterization are applicable for stomatal conductance estimation for Pinus cembra when the measurements do not modify the temperature regime.