Long-term Changes in Dwarf Pine (Pinus mugo) Cover in the High Tatra Mountains, Slovakia

The present article focuses on the distribution of Pinus mugo under conditions in the central Tatra Mountains, the main mountain range of the Western Carpathians. We analyze the response of P. mugo distribution to selected abiotic habitat conditions in the eastern Tatra Mountains. The study also compares data on the distribution of P. mugo in the higher central Tatras and in the hills of the western Tatras published in previous studies. The source data for this study were aerial photographs from 3 periods (1955, 1986, and 2002). Mountain areas covered by mountain pine were identified and analyzed by ArcGIS 10, and pine fields were classified with the help of the gray scale mode. A strip of mountain pine above the upper limit of the forest represents an easily identifiable boundary on the aerial photographs: 25 well-recognized localities were selected to examine the changes in the tree line in the eastern Tatras. The distribution of mountain pine increased in the central granite and eastern limestone Tatra Mountains from 1955 to 2002 at all monitored sites. The percentage of total surface area covered in P. mugo increased from 28.11% in 1955 to 34.74% in 1986 and to 39.01% in 2002. The study also analyzes the dispersal of mountain pine over 40 years in relation to elevation, slope, radiation aspect, flow accumulation, and vertical and horizontal curvature. The results of this study explain ongoing vegetation changes and are of importance as a contribution to monitoring of climate change in the mid-European mountain areas.


Introduction
High mountain ecosystems are unique areas for detection of climate change and assessment of climate-related impacts (Beniston 2003). The main ecological driving force is climate, with temperature and the duration of the snow cover as key factors (Gottfried et al 1999). Climate change associated with global warming is more pronounced at higher elevations than at low elevations (Beniston and Rebetez 1996;Diaz and Bradley 1997;Giorgi et al 1997). The effect of elevation on surface warming is especially marked in the winter and spring seasons because it is mostly associated with a decrease in snowpack and is thus enhanced by the snow-albedo feedback (Giorgi et al 1997). Changes in air temperature can extend the length of the average annual growing season (Menzel and Fabian 1999) and can also cause a shift in phenology (Parmesan and Yohe 2003;Visser and Both 2005). Many studies show that high mountain ecosystems are vulnerable to climate change (eg Theurillat and Guisan 2001;Dullinger et al 2003aDullinger et al , 2003bDirnbö ck et al 2011). Climate change mainly affects the distribution of plant and animal communities (Beckage et al 2008), and, under expected climate scenarios, in the final perspective, results in the loss of rare species of alpine habitats (Dirnbö ck et al 2011). In global meta-analyses by Parmesan and Yohe (2003) and Root et al (2003), significant range shifts toward the poles or toward higher altitudes were documented for many organisms. Large parts of these changes may be attributed to increased global temperatures. In general, we expected that climaterelated changes in mountain ecosystems would be most pronounced in the ''ecoclines'' (boundary ecosystem) or the ecotones (Theurillat and Guisan 2001).
Scenarios of upward plant species and vegetation shifts have been widely discussed in many research articles. Theurillat and Guisan (2001) published a review that discusses this matter and concludes that, although the alpine vegetation can tolerate an increase of 1-2uC of average air temperature, in the case of a sharper increase, we can expect major changes. Dullinger et al (2012) showed that projected extinctions of high-mountain plants under 21st century climate change have a significant long-term component, and they also believe that endemic Alpine species seem to face the highest range losses. Loss of diversity in alpine communities and fragmentation of plant populations caused by climate warming is expected for comparable high mountain areas around the world (Grabherr et al 1995;Saetersdal et al 1998). Results  hypothesis that alpine plant species above the forest line will be affected by heavy fragmentation and habitat loss but only if the average annual temperature increases by 2uC or more. The loss of most of these Alpine plant species habitats is expected to be caused by the expansion of Pinus mugo in the Alpine zone. The same species of alpine heaths are declining due to expanded dwarf pine, for example, Carex sempervirens, Agrostis rupestris, Festuca supina, and others. Growth and fertility of P. mugo is mostly controlled by temperature (Dullinger et al 2004). Thus, the main limiting factor of P. mugo growth at high altitudes could be the soil temperature (Smith et al 2003), although Rossi et al (2007) refer to varying soil temperature thresholds at different sites, which indicates that soil temperature may not be the main factor that limits xylogenesis of conifers and provides strong evidence that air temperature is a critical factor that limits xylem cell production and differentiation at high altitudes. However, the air temperature alone may not be the dominant factor that determines the tree-line position, because the direct influence of temperature may be masked by interactions with other factors, such as precipitation, cold-induced photoinhibition, disturbance, or plant-plant interactions (Harsch et al 2009). This evidence, therefore, is inconclusive. Differences in expert opinions on this matter have led Smith et al (2009) to formulate 6 current hypotheses about the causes of upper tree-limit movement: climatic stress, mechanical disturbance, insufficient carbon balance, limitations in cell growth and tissue formation, limited nutrient supply, and limited regeneration. In the global meta-analyses by Parmesan and Yohe (2003) and Root et al (2003), significant range shifts toward the poles or toward higher altitudes were documented for many organisms. Large parts of these changes may be attributed to increased global temperatures. The expansion of tree-line-forming species (subalpine zone) to higher altitudes is evident in the Pyrenees (Camarero and Gutié rrez 2004;Peñ uelas et al 2007), in the Alps (Dullinger et al 2003a(Dullinger et al , 2003bGehrig-Fasel et al 2007;Vittoz et al 2008), in the Carpathians (Mihai et al 2007;Martazinova et al 2009;Š vajda et al 2011), and also in Sweden (Kullman 2002) and Caucasus (Akatov 2009).
Generally, climate conditions and land use in high mountain areas have been shown to influence the distribution of mountain pine. The potential model of timberline is based on the assumption that climate change as a factor in forest regeneration is primarily responsible for moving the upper limit of the natural forest above the original climatically determined timberline, whereas the abandonment of farming in the country is assumed to be the dominant factor that determines forest regeneration below this line. The aim of this study was analysis of mountain pine cover in the central and eastern Tatra Mountains and of the abiotic conditions significantly related to P. mugo increments.

Study area
The Tatra Mountains are situated at the Slovak-Polish border (20u119E; 49u119N) and are the highest mountain massif within the Carpathian Range of Central Europe. The highest summit reaches 2656 m; the massif is classified as a high-mountain landscape covered by subalpine and alpine zones. The study area (Figure 1) is situated in the central and eastern Tatra Mountains, composed of the high Tatras and Belianské Tatras. The geology of the study area is based on crystalline bedrock. The high Tatra Mountains contain from significant biotite tonalities to granodiorites, and locally porphyritic and also porphyritic granodiorites to granites. The Belianské Tatras contain various shales, sandstones and dolomites, dark limestones, and quartzites (Biely et al 2002).
The vegetation of the Alpine zone is dominated by Alpine meadows (dry tundra with mostly Festuca picturata, Luzula alpino-pilosa, Calamagrostis villosa, and Juncus trifidus), with patches of dwarf pine (P. mugo) and an increasing percentage of rocks (bare or covered with lichens, commonly Rhizocarpon, Acarospora oxytona, and Dermatocarpon luridum) above the upper tree line of 1800 m (Vološč uk 1994).
The average annual air temperature decreases with elevation by 0.6uC per 100 m, being 1.6 and 23.8uC at elevations of 1778 and 2635 m, respectively (Konč ek and Orlicz 1974). The amount of precipitation increases as elevation rises, varying from approximately 1.0 to approximately 1.6 m/y between 1330 and 2635 m but reaching .2.00 m/y in some valleys (Chomitz and Š amaj 1974). Precipitation is generally higher in the northern part than in the southern part of the mountains, as is runoff, which averages 1.42 and 1.57 m/y for the south and north, respectively (Lajczak1996). Snow cover usually lasts from October to June at elevations .2000 m. The climatedriven tree line in the Tatra Mountains is located at approximately 1550 m and partly includes natural ecotones, with individual conifers that reach ages of 350-450 years (Bü ntgen et al 2007). P. mugo is an obligatory prostrate pine with adult canopy height that varies between 0.3 and 2.5 m in the study area. The typical dwarf pine altitudinal (subalpine) zone extends from 1500 to between 1850 and 1900 m. Mountain pine zone developed especially in the western Tatras with glacial-meadow relief, with great antierosion and water retention potential. Closed mountain-pine thickets stretch up to 300 m above the timberline, which reaches approximately 1600-1750 m in the Tatras and encompasses the upper part of the forest alpine tundra ecotone. Mountain pine plays a significant role in the natural environment: it protects the soil and stabilizes the snow cover, thus restricting the release of avalanches, and it provides habitats for many species of flora and fauna (Jodowski 2006).

Material and methods
The present analysis was carried out by using geographic information system (GIS) (ArcGIS 10), based on aerial photographs from 1955, 1986, and 2003. The aim of the analysis was to verify temporal trends in the distribution of P. mugo and to investigate which environmental variables best explain the changes in the growth and distribution of the mountain pine. The applied modeling approach is based on 3 major assumptions: (1) the abiotic factors are assumed to be the major driving force of species distribution changes as well as of postgrazing succession, (2) rhe models are calibrated by using field data and thus comprise any competitive constraint a species may force upon or experience from its neighbor, (3) the speed of plant migration is consistent with that of climate change so that plant communities are in a permanent equilibrium with their environment (Dirnbö ck et al 2003).
Aerial images from 2002 had a reference system (sjtsk); aerial photographs from 1986 and 1955 (without reference system) were georeferenced on the basis of orthophotos (Table 1). For our analysis, we selected 25 localities from the study area. The main factor in selection of study localities was spatial distribution of P. mugo with respect to a variety of abiotic conditions and shadowed areas on the aerial photographs. Mountain pine fields were extracted according to Š vajda et al (2011) from the aerial photos in gray scale and then reclassified into the range that represents mountain pine occurrence in the study area. Each photo was examined individually. If mountain pine on the slide was gray, with a value from 75 to 110, all such values in the range were reclassified as 1 (pine cover). The remaining values from 0 to 75 and 110 to 256 were reclassified as 0 (no pine cover). Intervals for segmentation were separately fixed for each image on the previous identification of pixel values that represented patches of mountain pine. We created a grid in which each pixel contained either the value 1 or the value 0. Then, the grid was automatically vectorized on the basis of the 2 values ( Figure 2).
Habitat conditions were spatially simulated by using GIS, a digital elevation model, meteorological data, and existing maps. The most significant factors that explain   the presence of P. mugo according to Dirnbö ck et al (2003) are the daily temperature, followed by slope, geology, solar radiation in September, and duration of snow cover.
To test this hypothesis, it was necessary to create an explicit temporal and spatial model of the spread of mountain pine and to analyze their sensitivity to predicted climate change trends. The size of the pixels' (cell) grid was equivalent in all images, because each image was adjusted to the same cell size through the transformation of the grid as well as during georeferencing. Thus, all images and the grids had the same pixels. A digital elevation model of the study area was used for representation of selected abiotic habitat conditions. Finally, we created a grid (25 3 25 m) from raster maps of abiotic conditions overlaid with maps of mountain pine cover. Each square represented a row in the database, with unique values to identify number, site number, cover of mountain pine (1955,1986,2003), elevation, slope, orientation, radiation, water flow, and horizontal and vertical curvature. GIS intersection of study sites divided the sites studied into 23,526 smaller areas with unique characteristics related to pine increase. Two sites (nos. 3 and 11) ( Table 2) were excluded due to the lack of data from 1986. A single matrix was analyzed. Independent variables were as follows: elevation (m); slope (degree); radiation as the whole energy of solar radiation in the vegetation period from May to November (watt hours per square meter); aspect (degree); flow accumulation, accumulated flow to each cell as determined by accumulating the weight for all cells that flow into each downslope cell; vertical curvature, curvature of the surface perpendicular to the slope direction (a positive value indicates that the surface is upwardly convex at that cell; a negative plan indicates that the surface is upwardly concave at that cell), horizontal curvature, and curvature of the surface in the direction of slope (a negative value indicates the surface is upwardly convex at that cell; a positive profile indicates the surface is upwardly concave at that cell). Dependent variables were expressed as relative increments of mountain pine during observed periods 1955-1986 (relative increments of mountain pine from 1955 to 1986 in m 2 ) and 1986-2002 (relative increments of mountain pine from 1986 to 2002 in m 2 ). The increments in mountain pine were reported as means and standard deviations for potential comparison with other studies, but the values showed a highly skewed distribution in most sample groups. Therefore, a nonparametric approach to analysis of the data was necessary. The significance of differences between groups was tested by using the Kruskal-Wallis nonparametric test. When P , 0.05, the data were considered as significantly different. The principal component (PC) analysis, correlation matrix, a multivariate technique, was used to extract the potential relationships between the variables studied. Principal components are linear combinations of original variables (slope, elevation, radiation, flow accumulation, vertical curvature, horizontal curvature, and relative increase of pine during observed periods), each axis being statistically orthogonal to the others. Integration of the variables slope and elevation (m) in different periods enabled us to follow different processes of mountain colonization by mountain pine during the respective periods. Because this statistical technique produces statistically orthogonal axes, we were able to examine potentially independent biological phenomena. We used 8 variables; consequently, we evaluated 8 principal components. The proportions of the total variance accounted for by each component are shown in Table 3 (see results).  Table 2). The results also indicate that the mean increase in mountain pine surface cover for all periods was approximately 0.2% per year (0.21% first period, 0.27% second period) for the total surface area, but results in relation to selected abiotic conditions still showed some differences.

Mountain pine cover in the Tatra
From 1955 to 1986, mountain pine cover probably increased in areas that were subject to grazing. It is likely that this is determined by the PC1 component (Table 3). Mountain pine in the earlier period expanded mainly at lower altitudes with lower slope and expanded in the ridge portions of these areas (vertical versus horizontal curvature). The effect of intensive growth of mountain pine occurred in the period 1955-1986. In the later period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002), this effect decreased (PC1) ( Table 3). PC1 can be designated as a factor in the ingrowing of mountain pine areas in the lower parts of the mountains, mainly due to abandoned grazing ( Figure 4A, B). PC2 (Table 3) is a factor that mainly describes the variability of dwarf pine growth rates in both periods. This happened in the ridge portions of the mountains. The phenomenon of scattering of mountain pine on ridges is especially significant on the western slopes of hills ( Figure 5A). Therefore, we can conclude that mountain pine tended to spread very intensively, mainly on the western and southwestern ridges of the central and eastern Tatra Mountains. The spread of mountain pine on the western ridges is significantly different from that on the eastern ridges ( Figure 5A). PC3 and PC4 are the factors that primarily reflect the mutual variability of the abiotic variables. PC5 ( Figure 5B; Table 3) mainly describes the trend in the growth of mountain pine in the later period, where the effect of mountain pine expansion is visible due to better conditions with enough water. This relates to expansion of mountain pine to territory with a higher altitude. Component PC6 ( Figure 5C) is relatively important for understanding of differential scattering of mountain pine in different periods. It shows the spread of mountain pine by slope, independent of altitude. For mountain pine, steeper places were more suitable in the earlier period, and, therefore, it expanded from more level areas onto steeper positions. In the later period, 1986-2002, mountain pine colonized more level places. This effect also differs according to the aspect ( Figure 5C). The ''steeper slope'' strategy of pine growth known in the earlier period was southerly oriented (compare Figure 4C), whereas the ''level strategy'' from the years 1986-2002 was westerly oriented (compare Figure 4D). A combination of the necessary amounts of hours of daylight and proper temperatures is probably very important for the spread of pine. In the second period, when mountain pine extended to the north, northwest, and western areas, it was distributed mainly in the more level places. This effect is also independent of altitude. Obviously, north-west and westerly orientations once again play an important role, which allows the species to spread to less favorable habitats. The PC7 and PC8 components describe the intervariability among climatic factors, which is not significantly related to the growth of pine. We also tested regression analysis-related components with a dependent variable (P. mugo increments), but this regression analysis is not shown because it led to the same results.

Discussion and conclusion
At high altitudes, vegetation is under constant environmental stress and thus abiotic conditions become more important for community development than biotic relationships (Pauli et al 1996). Kö rner (2003)   MountainResearch solar radiation, slope and exposure, and plant stature. Additional components that modulate the interplay of the 3 main drivers are wind velocity; ambient air temperature above the vegetation; and soil properties such as surface structure, moisture, and thermal conductivity. Dwarf pine is a tree species that is well adapted to high variability of climatic conditions in the high mountains. The most significant factors that explain the presence of P. mugo are the daily temperature, followed by slope, geology, solar radiation in September, and duration of snow cover (Dirnbö ck et al 2003). Dispersal of pine in the central and eastern Tatra Mountains is westerly oriented. The air temperature increases gradually during the day, and, by comparison with eastern slopes, when westerly oriented slopes are exposed to the sun, the temperature is usually high, the snow melts, and mountain pine have enough heat, water, and sunlight for photosynthesis. The onset of photosynthetic functions depends mainly on weather conditions in early spring (Lehner and Lü tz 2003), when wood cell production starts in April and May, followed by maximum cell production in June-July (Bü ntgen et al 2007). The autumn effects (September) may also play an important role in pine development (Dirnbö ck et al 2003). A combination of precipitation, temperature, and the timing of snow melt is evidently responsible for early pine wood production. Water supply is important during the first part of the growing season (Bü ntgen et al 2007). In the main chain of the Tatra Mountains, the western aspect and westerly oriented slopes seem to be crucial for effective dispersal of this wood species.
Similar to other high mountains, the Carpathians can show climate change trends and possible shifts of vegetation types with altitude. Climate development can to some extent affect phenological trends (Bauer 2006;Š kvareninová 2008), which reflect the changing climate conditions. The onset of individual phenological stages and their progress is mainly influenced by air temperature as well as temperature and humidity of soil and other meteorological variables (Š kvareninová 2009). Over the past 50 years, summer temperatures in the Tatra Mountains have increased by 0.7uC at higher elevations and by 1.4uC at lower elevations. Winter temperatures have increased by 1.4uC at higher elevations and by 1.9uC at lower elevations (Melo 2005). The temperature limit of the mountain pine zone is determined by the biotemperature threshold in the range of 3.0 to 2.0uC (maximum uC to minimum uC) (Miň dá š et al 2003). When considering the rate of the current temperature change (2-3uC in 100 years), we can expect more changes in  Boltižiar (2007) suggests extension of mountain pine cover, advance of forest, and reduction of grassland areas. Martazinova et al (2009) conducted research on grasslands above the upper forest limit in the Ukrainian Carpathians. Grass cover significantly decreased at the sites with conifers. Spruce stands moved to higher altitudes, mainly on the northern slopes, whereas the beech stands in the same area on the southern slopes did not show any significant movement. Apparently, the greatest changes were recorded at those sites where upper forest limit was marked at higher elevations. In a study of alpine, subalpine, and forest landscapes in the Iezer Mountains (Southern Carpathians), Mihai et al (2007) described how mountain pine-subalpine associations developed and gradually covered subalpine meadows and barren land (between 1986 and 2002, colonization averaged 0.14 km 2 /y). However, the mountain pine area has lost some lower stands because of spruce forests, which increased in elevation. This is largely a feature of southern aspect slopes (sunny), where the natural timberline is higher under some local conditions. Another example from the Austrian Alps (Dullinger et al 2004), which ran a model with projections for 1000 years, predicted that the area covered by pines will increase from 10% to between 24% and 59% of the landscape studied. The shape of the dispersal curve and spatial patterns of competitively controlled recruitment suppression affect range size dynamics at least as much as does variation in assumed future mean annual temperature (between 0 and 2uC above the current mean). Moreover, invasibility and the shape of the dispersal curve interacted with each other due to the spatial patterns of vegetation cover in the region.
Upward pine shifts above the potential regional tree line will be primarily influenced by climate change (Gehrig-Fasel et al 2007), whereas upward shifts below the potential regional tree line are usually interpreted as primarily influenced by land abandonment. Dullinger et al (2003a) indicated that a shift of tree and shrub species caused by land use and expected climate change can be expected in the European Alps. Abandonment of pasture will allow invasive expansion of P. mugo scrubs to new areas. Our study clearly offers the field data for time scheme prognoses in the Alps, because pasturing in the central Tatras was forbidden and abandoned in the 1950s. In peripheral areas, this process will be dependent on the competitive struggle with abandoned grasslands for light after grazing has ceased. In recent decades, dwarf pine forest has generally lost total surface area as the result of pressure from lower vegetation communities and even secondary pastures (Mihai et al 2007).
In addition to climate change, human land use may drive changes in the tree line. Land use in subalpine and alpine areas (grazing and extraction) affects the distribution of flora just as much as climate. Since the 13th and 14th centuries, anthropogenic land cover change has involved clearing mountain-pine thickets to obtain new pastures for sheep and cattle grazing, for extensive charcoal and oil production, and for copper and iron-ore mining, sometimes leading to degradation. Jodowski (2007) described how establishing national parks in the Tatras, Babia Gó ra, and Giant Mountains, enabled secondary succession, which has led to colonization of previously abandoned habitats. However, these processes have been hampered by harsh edaphic and climatic conditions as well as by avalanches and debris flows. Extensive planting of mountain pine in the former Czechoslovakia significantly facilitated the regeneration of mountain pine thickets. After the absolute restriction of grazing in some national parks, we observed progressive long-term trends in secondary succession and patterns of plant establishment driven by climate. Habitats in the peripheral or isolated mountain belts at or above the tree line are generally rich in diversity of endemic species. In these habitats, tree-line expansion disproportionally reduces the habitats of highaltitude species. Such legacies of climate history, which may aggravate extinction risks under future climate change, have to be expected for many temperate mountain ranges (Dirnbö ck et al 2011).
In the Tatras, mountain pine is westerly oriented. Longer growing seasons, milder winters, and shorter duration of snow cover create favorable conditions for the growth of mountain pine. This shift has not only had a FIGURE 4 Increments in mountain pine cover in the Tatra Mountains in the periods 1955-1986 and 1986-2002, according to (A and B) elevation and (C and D) aspect. In both periods, the tempo of growth differed at the different elevations from the period 1955-1986: Kruskal-Wallis analysis of variance (ANOVA) nonparametric test at H (33, N 5 13,594) 5 1606.419, P 5 0. 000, 1986-2002: H (33, N 5 13594) 5 378.0978, P 5 0.000. At lower elevations, P. mugo tended to grow more intensively in the earlier than in the later period (compare Figure 1A and 1B, PC1 in Table 3) (C) In the period 1955-1986, the groups did not differ according to aspect. Kruskal-Wallis ANOVA H (7, N 5 13586) 5 125.91, P 5 0.000. (D) In the period 1986-2002, the westerly aspects were significant, Kruskal-Wallis ANOVA H (7, N 5 13586) 5 970.90, P 5 0.000.

MountainResearch
devastating effect on alpine plant communities due to habitat loss but also due to greater fragmentation, which ultimately may affect the population of different animal species dependent on these habitats. Further understanding of dispersal, persistence, and survival strategies of dwarf pine in the other mountain in Slovakia is also required. We will continue to monitor dispersal of P. mugo and extend our studies to the low Tatras. Continued research on vegetation dynamics in Slovakia's mountain areas is needed in light of the significance of vegetation in the context of global change. This work will help to characterize and evaluate the total tree surface area as a basis for the state Nature Conservancy's management of mountain national parks and protected areas in Slovakia.