1 Introduction
Mountain altitudinal vegetation is very sensitive to climate change in the form of global warming (
Steinbauer et al., 2018;
Hagedorn et al., 2019), as exhibited by the reoccurrence of succession in many areas (
Chen et al., 2011;
Guo et al., 2018;
Elsen et al., 2020;
Zu et al., 2021;
Rana et al., 2022). The dynamic changes in mountain altitudinal vegetation belts are considered to be the result of climate change (
Li et al., 2021;
Zhang et al., 2021). In particular, the migration of mountain vegetation to higher altitudes has been studied as a response to global warming (
Parmesan and Yohe, 2003;
Lenoir et al., 2008;
Chen et al., 2011;
Hagedorn et al., 2019;
Zischg et al., 2021;
Zu et al., 2021;
Rana et al., 2022), and plant species have been observed migrating to higher altitude in the European Alps (
Lamprecht et al., 2018) and the Chimborazo volcano in Ecuador (
Morueta-Holme et al., 2015), whereas the migration to lower altitudes can be explained as a response to changes in regional water balance, for example, the altitude of plant species in California decreased significantly (
Crimmins et al., 2011). An increase in water resource utilization by vegetation may offset the impact of warming on distribution. However, recent studies have shown that, at different scales, the direction and magnitude of the spatial migration of mountain vegetation in response to climate change is inconclusive. The majority of previous studies have reported that global warming causes plant species to move upward along the altitudinal gradient (
Gottfried et al., 2012;
Zu et al., 2021), however, studies have also reported the downward movement of plant species along the altitudinal gradient, as observed in the Gongga Mountains of Sichuan where 63.9% of the species moved upward and 22.9% of the species moved downward (
Zu et al., 2021).
As unique geographical units, mountain ecosystems exhibit spatial variation in environmental factors owing to changes in altitude, slope, and aspect associated with the vertical gradient (
Zhang et al., 2021). In the context of climate change, mountain vegetation has undergone local extinction and secondary succession, compelling plants to migrate to more suitable habitats (
Kelly and Goulden, 2008;
Giménez-Benavides et al., 2011) and driving changes in the distribution of vegetation along vertical gradients. In addition, studies have reported that, in different mountain ecosystems, mountain vegetation at different altitude gradients also varies with climate change (
Pauli et al., 2012;
Zu et al., 2021). Climate change affects the migration of vegetation by altering precipitation and temperature-related climatic factors. For example, the global distribution range of
Aletris spicata (Thunb.) Franch. was chiefly explained by changes in precipitation, and its altitude range shifted upward by 370 m in the last decades (
Zu et al., 2021). Therefore, understanding the spatial and temporal distribution of and dynamic changes in mountain vegetation in different regions and their climatic driving factors is crucial for assessing the impact of climate change on mountain ecosystems.
The Tianshan Mountains, known as the “water tower of Central Asia”, provide water sources for major rivers in the arid regions of Central Asia. They are also an important ecological barrier in the arid regions of northwest China (
Li et al., 2020) and provide important ecological services for the social and economic development of arid oases (
Jiao et al., 2019). Over the past half-century, the temperature increase in the Tianshan Mountains has been higher than the rate of temperature change in the Northern Hemisphere as well as the global increase (
Li et al., 2021), and the Tianshan Mountains have shown varying degrees of warming and wetting since the 1980s (
Sánchez-González and López-Mata, 2005;
Zhang et al., 2022). In addition, as a typical representative of mountain ecosystems in arid regions, not only do these mountains provide a location for human activities but they also perform important carbon sink functions, which are of great importance to the local ecological and economic development, human society, and global carbon cycle (
Ni, 2004). However, the response of the Tianshan Mountain ecosystem to global climate change remains unclear. In particular, the response to climate change of vegetation belts at different altitudes on the southern slope of the Tianshan Mountains is still unknown; thus, further research is urgently required. Studies in this field have important scientific implications for environmental protection and monitoring regional and global carbon cycles in the arid region of northwest China. Therefore, this study examined Landsat data from the Google Earth Engine (GEE) cloud platform to analyze altitudinal vegetation belts on the southern slope of the Tianshan Mountains. We aimed to answer the following key questions. 1) What has been the migration trend of the forest belt on the southern slope of the Tianshan Mountains over the past 30 years? 2) How has the forest belt on the southern slope of the Tianshan Mountains responded to climate change?
2 Materials and methods
2.1 Study area
The study area was located in the Tianshan Mountains in Central Asia, near the center of Eurasia (Fig.1). It spans 1700 km from east to west, with an average elevation of 2300 m. The large mountain system influences the atmospheric circulation in the region, resulting in clear differences in the natural climate between the northern and southern Tianshan Mountains. The southern region has a warm temperate climate, whereas the northern region’s climate is temperate. Additionally, this area is affected by westerly winds from the North Atlantic, resulting in a unique humid climate in the Tianshan Mountains, which are also known as the “Central Asian wet island” (
Li et al., 2020;
Yue et al., 2020;
Zhang et al., 2021). The unique climate of the Tianshan Mountains results in well-developed altitudinal vegetation belts, thereby forming the most representative mountain altitudinal vegetation belt in the temperate arid regions of the world (
Zhang et al., 2021). The Tianshan Mountains are composed of a series of east–west parallel mountains and their corresponding basins and valleys. Consideration was given to minimizing the impact of human factors in the selection and division of the study area. Although vegetation change results from the combined effect of climate change and human activities (
Ma et al., 2023), global warming is considered the main driver of changes in terrestrial ecosystems (
Urban, 2015;
Feeley et al., 2020;
Trisos et al., 2020;
Antão et al., 2022), which determines the survival and distribution of surface vegetation. Therefore, based on the first law of geography that “everything is related to everything else, but near things are more related to each other” (
Tobler, 1970), we avoided selecting a study area near the center of human activity in the city. Instead, we focused on three typical mountain ranges (Jengish Chokusu, Bogeda Mountain, and Barikun Mountain, from west to east) on the southern slope of the Tianshan Mountains (Fig.1) to study the spatial dynamics of altitudinal vegetation belts over the past 30 years. Jengish Chokusu is located near the border between China and Kyrgyzstan, at an elevation of 7443.8 m. It is the highest peak in the Tianshan Mountains and is listed as a national nature reserve. Bogeda Mountain, at an elevation of 5445 m, is one of the most prominent mountains in the western part of the East Tianshan Mountains. Barikun Mountain is one of the most prominent mountains in the eastern section of the East Tianshan Mountains, stretching 160 km, with an average elevation of 3300 m. A large area, including the main peak, was selected for analysis based on the direction of the mountains.
2.2 Data resources
Owing to the rapid development of remote sensing technology and the emergence of the Google Earth Engine (GEE) cloud platform, vegetation cover types can be quickly analyzed with long temporal scales, large spatial scales, and at high spatial resolution (
Zurqani et al., 2020). Therefore, from the data provided on the GEE platform, we selected seven-time nodes (1990, 1995, 2000, 2005, 2010, 2015, and 2020) and used the random forest (RF) classification method to divide the vegetation into four types: desert, grassland, forest, and ice–snow. However, five distinct vegetation belts are present in the study region as deserts appear at two altitude ranges simultaneously. According to field investigations, the higher altitude range is considered an alpine desert. Therefore, the spatial distribution of the five altitudinal vegetation belts on the southern slope of the Tianshan Mountains during different periods was obtained.
The remote sensing data used in this study consisted of Landsat Thematic Mapper/Enhanced Thematic Mapper/Operational Land Imager (TM/ETM/OLI) images from 1990 to 2020 in the Landsat Collection 2 Level-2 data set, which was accessed through GEE. To facilitate the accurate extraction of vegetation belts, the median of the image set was used as a single image with a spatial resolution of 30 m, mainly from the vegetation growth season (June–September), for the specific acquisition time (Fig.2). The revisit period of Landsat data is 16 days, but it is discontinuous in time due to frequent cloud pollution and other factors. The data loss and stitching completion are shown in Tab.1. For the years with missing data or large cloud cover, the image mosaic of adjacent years is used to fill in. The specific method is to first take the research year as the center, and then look for the next year to find whether there is a suitable image. If there is, the image of the year will be used to fill in the holes of the research year. On the contrary, the research year is used as the center to look forward to one year, and finally, a suitable and effective image is obtained.
The China Land Cover Data set (CLCD) contains multi-temporal land use classification data from China, with a spatial resolution of 30 m and a high level of classification accuracy, which can be used as potential training samples (
Yang and Huang, 2021). The sample data are mainly based on CLCD, and its spatial distribution is shown in Fig.3, Fig.4, and Fig.5. Among them, there are 520 samples per period in the Jengish Chokusu, a total of 7 periods; there are 200 samples per period in the Bogeda Mountain, a total of 7 periods; there are 160 samples in each period of Barikun Mountain, a total of 7 periods. In addition, to improve the accuracy of the sample points, Google Earth Pro historical images were used for correction. The selection process of training samples is as follows: First, a certain number of random sample points are generated in the study area, and then the land classification attributes of the sample points are extracted through the CLCD product. Finally, the sample points with errors and fuzzy attributes are eliminated by Google Earth Pro, and the training samples are finally obtained. The historical image tool in Google Earth Pro can view images of different periods, but in fact, Google images are not continuous or homogenous but are mosaicked using multiple images from different periods, different spatial resolutions, and from different imagery providers (Lesiv et al.,
2018;
Nwilo et al., 2022). The images are compiled from a wide variety of sources such as SPOT 5, Rapid Eye, Earth Resource Observation Satellites (EROS), Meteosat 2, Geoeye 1, and Digital Globe World View 2 satellite (
Buka et al., 2015). Because Google Earth images come from multiple information sources, they do not have the same spatial resolution (Goudarzi and Landry, 2017). In the correction process, we mainly consider the historical Google Earth images with the highest possible resolution (15 m−10 cm) from 1990 to 2020.
A digital elevation model (DEM) was used to extract the elevation of vegetation, mainly from the Shuttle Radar Topography Mission (SRTM) 30 m spatial resolution DEM data provided in GEE.
Meteorological data included monthly average temperature (
Peng, 2019), monthly precipitation (
Peng, 2020), and monthly potential evapotranspiration (
Peng, 2022). The spatial resolution of the three data is 1 km, and the data format is NetCDF, that is .nc format. These data are mainly used to generate relevant bioclimatic variables. Spatial interpolation data will result in greater uncertainty due to interpolation methods and random errors (
Huang et al., 2019). To this end, we use the temperature and precipitation observation data in the Tianshan Mountains to evaluate the accuracy of this interpolated climate data. The interpolation data of temperature and precipitation are significantly correlated with the observation data, indicating that the study of long-term climate trends in the Tianshan Mountains is accurate and effective (Fig.6). We selected the temperature, precipitation, and potential evapotranspiration data of the National Tibetan Plateau Scientific Data Center from 1990 to 2020.
2.3 Random forest
Random forest (RF) is an ensemble learning method based on a decision tree, which combines several ensemble regression or classification trees (
Breiman, 2001a) and uses Bagging or Bootstrap aggregating algorithms to construct numerous different training subsets. Each decision tree provides classification results for samples that are not selected as training samples, and the final class is determined by the highest score (
Tian et al., 2016). Owing to the RF method’s high prediction accuracy and tolerance to outliers and noise, it is a non-parametric modeling tool with an adaptive function (
Breiman, 1996;
Yu et al., 2018); thus, it was chosen for this study.
First, an RF algorithm was applied to the classification of remote sensing images in GEE. The characteristic parameters involved in the classification mainly comprised 16 variables. In the classification process, we randomly divide the sample points into a training set of 70% and a validation set of 30%, and then obtain the overall classification accuracy and kappa coefficient to evaluate the classification results. The variables and their importance in the classification process are shown in Tables A1–A3, and the classification accuracy is shown in Tables A4–A6. When using the random forest classification method, we should pay attention to its uncertainty. The uncertainty of the random forest classification method mainly comes from data uncertainty and model uncertainty. The random forest generates multiple training subsets from the original data set through self-sampling, but the data distribution of each subset may be different. This may lead to differences between the trained decision trees, which in turn affects the prediction performance of the random forest. Each decision tree in the random forest is constructed independently, and the features are randomly selected for splitting. The randomness of feature selection can increase the diversity of the model, but it may also lead to the inability of some decision trees to capture key features, making each decision tree capture information slightly different, resulting in model uncertainty. In general, in the process of using random forest classification, multiple tests are required, including the depth of the tree, the selection of classification features, etc., to optimize the performance of the model.
Next, RFs were applied to the regression of vegetation belts in response to climate change and environmental factors. As climate change involves multiple dimensions, the relative importance of different climatic variables differs with vegetation type, space, and time. Asymmetric changes in climatic conditions may lead to complex responses among the species in regional communities (
Antão et al., 2022). To evaluate the importance of climate variables, it is important to select biological response variables with known relationships to climate variables (
Greenberg et al., 2015), namely bioclimatic variables. Precipitation, temperature, and potential evapotranspiration are known to have a substantial effect on vegetation change (
Ma et al., 2023). Based on this understanding, with reference to the WorldClim bioclimatic variables, nine bioclimatic variables were selected for RF regression analysis: annual potential evapotranspiration (APET), annual precipitation (AP), mean annual temperature (MAT), temperature seasonality (TS), precipitation seasonality (PS), precipitation of coldest quarter (PCQ), precipitation of warmest quarter (PWQ), mean temperature of coldest quarter (MTCQ), and mean temperature of warmest quarter (MTWQ). MAT, AP, and APET are very important for forest growth and are widely used in the study of forest growth under the background of climate change (
Yu et al., 2014;
Zhang et al., 2016;
Brandt et al., 2017;
Liu et al., 2021). In addition, we used precipitation and temperature in the warmest and coldest seasons (corresponding to growing and non-growing seasons) to characterize the difference between the growing and non-growing seasons. Seasonal temperature and seasonal precipitation are also important factors for the survival and development of organisms (
Shrestha et al., 2018;
Ge et al., 2019). These bioclimatic variables are calculated based on monthly average temperature (
Peng, 2019), monthly precipitation (
Peng, 2020), and monthly potential evapotranspiration (
Peng, 2022). TS is defined as the amount of temperature variation over a given period based on the ratio of the standard deviation of the monthly mean temperatures to the mean monthly temperature (also known as the coefficient of variation (CV)). PS is a measure of the variation in monthly precipitation totals over the course of the year. This index is the ratio of the standard deviation of the monthly total precipitation to the mean monthly total precipitation (also known as the coefficient of variation). PCQ approximates the total precipitation that prevails during the coldest quarter. PWQ approximates the total precipitation that prevails during the warmest quarter. MTCQ approximates mean temperatures that prevail during the coldest quarter. MTWQ approximates mean temperatures that prevail during the warmest quarter. The response of forest belt with dense vegetation coverage in a typical area of the southern slope of the Tianshan Mountains to climate change was analyzed, and the relative importance of each factor was evaluated using an RF regression model. The important parameters of random forest regression are as follows: the number of trees is 100, the maximum depth of the tree is 10, the maximum depth of the leaves is 50, the minimum number of samples required for the re-division of internal nodes is 2, and the minimum number of samples for leaf nodes is 1. This type of regression analysis has certain limitations in explaining the relationships between variables. Although this model cannot explain whether a causal relationship exists between the variables, it is not hindered by the multicollinearity problem faced by general regression analysis. Furthermore, RF does not require variable selection, which is convenient for calculating the nonlinear effects of variables. Instead, this model reflects the interaction between variables and is suitable for evaluating importance (
Breiman, 2001b;
Li, 2013;
Zhang et al., 2018;
Yan et al., 2021). Another factor that we accounted for was the time lag in the response of vegetation dynamics to changes in climate (
Zhao et al., 2020). As the spatial migration of mountain vegetation is a long-term process, we averaged the bioclimatic indicators over time to avoid the effect of the response time lag.
2.4 The lower and upper limits of altitudinal vegetation belts
Changes in the spatial distribution of altitudinal vegetation belts are measured by the differences in the lower and upper limits of altitude (
Menéndez et al., 2014;
Zu et al., 2021). According to our field investigation, altitudes corresponding to 25% and 75% of the maximum number of pixels were set as the lower and upper limits of the vegetation belt, respectively. The upper and lower limits of the forest we extracted refer to the timberline of a closed forest.
3 Results
3.1 Analysis of climate change and spatial distribution of vegetation belts in the study area
Based on meteorological data, we evaluated climate change in three representative mountain regions. The climate records for Jengish Chokusu, Bogeda Mountain, and Barikun Mountain show that over the past 30 years, the three sub-regions experienced discernible climate warming trends (Fig.7(a)–Fig.7(c)), which have passed the significance test and increasing by (0.031 ± 0.005)°C/yr, (0.021 ± 0.005)°C/yr, and (0.030 ± 0.005)°C/yr, respectively. However, the precipitation fluctuations varied by region (Fig.7(d)–Fig.7(f)). The precipitation in Jengish Chokusu showed an upward trend, increasing by 0.776 ± 0.347 mm/yr. In contrast, the precipitation in Bogeda and Barikun Mountains showed a downward trend, decreasing by 0.842 ± 0.328 mm/yr and 0.932 ± 0.280 mm/yr, respectively, and passed the significance test. Moreover, the seasonal climate records of the three regions show that the study area was warm and humid in summer and cold and dry in winter during the study period (Fig.7(g)–Fig.7(h)). Therefore, the warmest season mentioned in this paper refers to summer, and the coldest season refers to winter.
Based on the random forest classification method, the spatial distribution of vegetation on the southern slope of the Tianshan Mountains was extracted (Fig.8–Fig.9). It was found that the spatial distribution of vegetation in the three study areas is basically the same, and different vegetation belts are radial from the center to the periphery. The ice and snow belt was followed by what we believe to be an alpine desert. This was followed by grassland, forest (which is surrounded by grassland), and lastly desert.The north and south slopes are divided by the ridge line, which is basically symmetric. The forest on the north slope appears in pieces, but the forest on the south slope is more fragmented, also known as spotted. Based on the Landsat images, the resolution of these patchy forest areas is at least equal to or greater than 0.009 km2.
Fig.11 shows that the distribution of altitudinal vegetation belts (represented as pixels on the map) in the three study areas on the southern slope of the Tianshan Mountains was consistent. Moving up the altitude gradient, the land types were: desert, grassland, forest, alpine desert, and ice–snow. A considerable decrease in grassland pixels was observed in areas occupied by forest pixels, indicating that the forest belts were spatially distributed within the grassland belt. Forest, alpine desert, and ice–snow belts showed an inverted bell curve similar to the shape of a Gaussian distribution, whereas desert and grassland pixels were more disordered with the altitude distribution curve. In addition, the pixels of each vegetation type clearly coincide with each other, indicating that there was a wide transition area between each vegetation belt. This phenomenon demonstrates the coexistence of different vegetation belts in the same altitude range.
3.2 Variation trend of forest belt in the study area
From the context, it can be seen that the forest belt, as an extremely unique vegetation belt on the southern slope of the Tianshan Mountains, is embedded in the grassland belt in a spot-like manner and shows discontinuous changes or replacements in the vertical gradient, which is different from the previous understanding. From the analysis of the change trend of the upper and lower limits of the typical regional forest belt on the southern slope of the Tianshan Mountains from 1990 to 2020 (Fig.12), the slope of the trend line of the upper and lower limits of the Jengish Chokusu and the Bogeda Mountain regional forest belt and the lower limit of the Barikun Mountain regional forest belt is greater than 0, and the upper limit of the Barikun Mountain regional forest belt is slightly less than 0, which indicates that the upper and lower limits of the typical regional forest belt on the southern slope of the Tianshan Mountains from 1990 to 2020. The trend of migration in the altitude direction is basically the same. Specifically, the upper limit of the Jengish Chokusu and the Bogeda Mountain regional forest belts on the southern slope of the Tianshan Mountains have migrated to high-altitude areas, while the upper limit of the Barikun Mountain regional forest belt has no obvious trend (a slight downward trend). The changing trend of the lower limit of the forest belt on the southern slope of the Tianshan Mountains was the same, and the three regions showed a trend of migration to high-altitude areas. In addition, for the Jengish Chokusu, the slope of the trend line of the upper and lower limits of the forest belt is greater than 0, and the slope of the trend line of the upper limit is greater than the slope of the lower limit of the forest belt, which indicates that the forest belt is expanding during the overall upward migration of the forest belt in the Jengish Chokusu. For the Bogeda Mountain, the slope of the trend line of the upper and lower limits of the forest belt is greater than 0, and the slope of the trend line of the upper limit is less than the slope of the lower limit of the forest belt, which indicates that the forest belt is shrinking during the overall upward migration of the forest belt in the Bogeda Mountain. For the Barikun Mountain, the slope of the upper limit trend line of the forest belt is close to 0, and the slope of the lower limit trend line is greater than 0, which indicates that the forest belt in the Barikun Mountain is shrinking, and the forest at high altitude is almost unchanged. The main change area is the forest at low altitude. Comparing the three study areas horizontally, it can be found that the slope of the trend line of the Bogeda Mountain is much larger than that of the other two areas, whether it is the upper limit of altitude or the lower limit of altitude, which indicates that the Bogeda Mountain area is the most dramatic area of the Tianshan forest belt. For the upper limit of the forest belt, the Barikun Mountain area is undoubtedly the most stable area; for the lower limit of the forest belt, the slope of the trend line in the T Jengish Chokusu is the smallest, which is the least affected.
3.3 Responses of forest belt to climate change in the study area
As an extremely special forest belt on the southern slope of the Tianshan Mountains, the factors affecting the distribution of the upper and lower limits have spatial differences. From the analysis of the importance of climatic factors affecting the distribution of the upper and lower limits of the forest belt on the southern slope of the Tianshan Mountains (Fig.13), it can be seen that in the Jengish Chokusu, the most important climatic factor affecting the upper limit of the forest belt is PWQ, followed by AP; the most important climatic factor affecting the lower limit of forest belt is MAT, followed by MTWQ. In the Bogeda Mountains, the most important climate factor affecting the upper limit of the forest belt is MTWQ; the most important climatic factor affecting the lower limit of the forest is AP, followed by MTWQ. In the Barikun Mountains, PWQ is the most important climate factor affecting the upper limit of the forest belt. The most important climatic factor affecting the lower limit of the forest belt is MTCQ. Horizontal comparison, the climatic factors that mainly affect the upper limit of the forest belt in the Jengish Chokusu and the Barikun Mountain are consistent, which is PWQ. The warmest season here corresponds to the growing season, indicating that the precipitation in the growing season is the main factor affecting the upper limit of the forest belt in these two regions. The Bogeda Mountain is different from the previous two regions. The climate factor that mainly affects the upper limit of the forest belt is MTWQ. Similarly, the warmest season here corresponds to the growing season, but the main factor is the MAT, and PWQ is no longer important. The climatic factors affecting the lower limit of the forest belt in Jengish Chokusu, Bogeda Mountain, and Barikun Mountain are very obvious, followed by MAT, AP, and MTCQ.
4 Discussion
4.1 Comparison of responses in forest belt to climate change
Mountain vegetation is extremely sensitive to climate change and mainly relies on the succession of vegetation types to adapt. Mountains have a characteristic vegetation distribution that transitions from deserts, forests, and grasslands to permanently covered snow belts (
Hagedorn et al., 2019). The response patterns of the upper and lower altitudes of vegetation belts with different vegetation cover to climate change vary, which leads to the complexity of spatial vegetation migration (
Grytnes et al., 2014;
Dolezal et al., 2016). Climatic factors drive vegetation changes, especially in the arid regions of north-west China, where ecosystems are fragile and sensitive to climate change. The forest belt on the southern slope of the Tianshan Mountains displays an unusual distribution pattern. Our investigation revealed that the forest is embedded in the grassland belt in the form of patches, showing discontinuous changes or replacements on the vertical gradient. As a plant community, forests should have suitable climatic conditions. Changes in cold, warm, dry, and wet climates directly or indirectly affect the geographical distribution of forests.
The examination of meteorological data showed that the climate in the Jengish Chokusu area had a warming and wetting trend, whereas the climate in the two sub-areas of Bogeda Mountain and Barikun Mountain showed a warming and drying trend. Climate warming is common in these three regions, yet they differ in moisture level between dry and wet climates. The trend of warming-wetting and warming-drying is significant. Although they have different climatic backgrounds, the response trends of the upper and lower limits of the forest belts in these areas to climate change were similar. The upper and lower limits of the altitude of the Jengish Chokusu and Bogeda Mountain forest belts, and the lower limit of the altitude of the Barikun Mountain forest belt tended to migrate to higher-altitude areas. Notably, the upper limit of altitude in the Barikun Mountain forest belt did not show a distinct trend. We concluded that the migration of the forest belt to the higher-altitude area on the southern slope of the Tianshan Mountains was similar to the current species migration.
4.2 Primary factors affecting changes in forest belts
The migration of vegetation to higher-altitude areas is an adaptation to climate warming (
Parmesan and Yohe, 2003;
Lenoir et al., 2008;
Chen et al., 2011;
Morueta-Holme et al., 2015;
Hagedorn et al., 2019;
Zischg et al., 2021;
Zu et al., 2021;
Rana et al., 2022), but this phenomenon does not imply that temperature is the primary factor leading to the migration of forest belts to higher altitudes. The limiting effects of other variables, such as water availability, are also very important (
Greenberg et al., 2015;
Villén-Peréz et al., 2020), and they may contribute more than temperature in some cases. An increase in water use efficiency may offset the effects of warming. The dry and wet climate in this region does not directly represent the actual use of water resources in the ecosystem. For example, on the southern slope of the Tianshan Mountains, numerous glaciers in the high-altitude ecosystem are present, and climate warming will lead to an increase in glacial-melt water in the region. However, another consideration is the highly abundant precipitation levels typically found at higher altitudes. Therefore, the main factors affecting the migration of the upper and lower limits of altitude in the three typical forest belts in this region should be investigated further. In addition, vegetation type, which is mainly related to the biological characteristics of the dominant vegetation, is also one of the factors affecting the migration of the vegetation belt in response to climate change.
From the random forest regression model, it can be found that the relative importance of the PWQ in the Jengish Chokusu and Barikun Mountains is the highest, while the relative importance of the MTWQ in the Bogeda Mountains is the highest. However, if taken together, the PWQ contributes the most, and the warmest season corresponds to the growing season. Therefore, the upper limit of the forest belt on the southern slope of the Tianshan Mountains is more correlated with the precipitation in the vegetation growing season. At the lower limit of the forest belt, the relative importance of the MAT in the Jengish Chokusu is the highest, and the relative importance of the MTCQ in the Barikun Mountain is the highest, all of which are the largest contribution of temperature-related factors, while the relative importance of the AP in the Bogeda Mountain is the highest, followed by the MTCQ. Similarly, on the whole, the largest contribution is the MTCQ and the MTWQ. Therefore, the lower limit of the forest belt on the southern slope of the Tianshan Mountains is more affected by temperature-related factors. The response of the forest belt to climate varies at different altitudes. In low-altitude areas, temperature-related factors contribute the most, while in high-altitude areas, precipitation in the warmest season contributes more. Attention should be paid to the prerequisites. The warmest season already represents temperature-related factors. The migration of vegetation belts is the result of many factors, not the result of a single condition. By comparing the upper and lower limits of the forest belt, it can be found that the temperature-related factor is the primary condition for the migration of the forest belt. With the increase in altitude, the effect of precipitation gradually becomes prominent, and sufficient precipitation in the growing season is conducive to the growth of trees. In general, the higher the altitude, the lower the temperature. Under the background of climate warming, the migration of forest belts to high altitudes is to resist the impact of climate warming. High-altitude areas are more sensitive to climate change (
Li et al., 2019;
Zheng et al., 2021; Mohd et al., 2022) before trees begin to grow, they also need to reach a certain accumulated temperature (
Wang et al., 2013). When the temperature conditions are met, the precipitation factor is highlighted. On the contrary, the sensitivity of low-altitude areas to climate change is reduced, so it may cover up the role of some factors.
5 Conclusions
The response of mountain altitudinal vegetation belts to climate change is mainly reflected by the succession of the dominant vegetation. The forest belt is an extremely unique vegetation belt on the southern slope of the Tianshan Mountains owing to the fact that it is embedded in the grassland belt in the form of patches. Under different climatic backgrounds in different regions on the southern slope of the Tianshan Mountains, the response of the upper and lower limits of the forest belt to climate change is almost the same, showing migration to high-altitude areas, but the main climatic factors affecting the migration of the upper and lower limits are different and have spatial differences. The upper limit of the forest belt is more affected by the precipitation-related factors in the warmest season (corresponding to the growing season). With lower limits of the forest belt, the contribution of temperature-related factors is greater.
As an important ecological barrier and development support in the arid region of northwest China, especially in the southern slope of the Tianshan Mountains, the forest belt migrating to high altitude is expected to break the original carbon balance, and climate change brings multiple uncertainties to it. Under the background of ‘double carbon’, this will be one of the main research directions in the future.