Assessing future vegetation trends and restoration prospects in the Karst regions of Southwest China

: To alleviate the severe rocky desertiﬁcation and improve the ecological conditions in Southwest China, the national and local Chinese governments have implemented a series of Ecological Restoration Projects since the late 1990s. In this context, remote sensing can be a valuable tool for conservation management by monitoring vegetation dynamics, projecting the persistence of vegetation trends and identifying areas of interest for upcoming restoration measures. In this study, we use MODIS satellite time series (2001–2013) and the Hurst exponent to classify the study area (Guizhou and Guangxi Provinces) according to the persistence of future vegetation trends (positive, anti-persistent positive, negative, anti-persistent negative, stable or uncertain). The persistence of trends is interrelated with terrain conditions (elevation and slope angle) and results in an index providing information on the restoration prospects and associated uncertainty of different terrain classes found in the study area. The results show that 69% of the observed trends are persistent beyond 2013, with 57% being stable, 10% positive, 5% anti-persistent positive, 3% negative, 1% anti-persistent negative and 24% uncertain. Most negative development is found in areas of high anthropogenic inﬂuence (low elevation and slope), as compared to areas of rough terrain. We further show that the uncertainty increases with the elevation and slope angle, and areas characterized by both high elevation and slope angle need special attention to prevent degradation. Whereas areas with a low elevation and slope angle appear to be less susceptible and relevant for restoration efforts (also having a high uncertainty), we identify large areas of medium elevation and slope where positive future trends are likely to happen if adequate measures are utilized. The proposed framework of this analysis has been proven to work well for assessing restoration prospects in the study area, and due to the generic design, the method is expected to be applicable for other areas of complex landscapes in the world to explore future trends of vegetation.


Introduction
Vegetation has considerable impacts on almost all land surface energy exchange processes acting as an interface between land and atmosphere. It affects local and regional climate [1,2] and hydrologic balance of the land surface [3]. The dynamics of vegetation have been recognized as being of primary importance in global change of terrestrial ecosystems [4,5]. Karst regions in Southwest China are some of the largest exposed carbonate rock areas in the world, covering about 540,000 km 2 , within which more than 30,000,000 people under poverty [6]. Under intensive anthropogenic disturbances, these fragile karst regions have been reported to undergo severe rocky desertification. This type of land degradation is a process where a karst area covered by vegetation and soil transforms into a rocky barren landscape, which is considered to be one of the most dangerous ecological and environmental problems in China [7][8][9]. Vegetation plays an important role in ecological conservation and restoration [10]. To alleviate the severe rocky desertification in Southwest China, national and local Chinese governments have implemented a series of Ecological Restoration Projects (ERPs) to improve vegetation and ecosystem conditions since the late 1990s. Examples are the Natural Forest Protection Project, the Grain to Green Program, the Public Welfare Forest Protection and the Karst Rocky Desertification Comprehensive Control and Restoration Project. A rapid and area-wide mapping of the dynamics and spatial patterns of vegetation changes as well as the prediction of future vegetation trends are thus essential for further ecological measures and sustainable development.
Remote sensing data have been increasingly used for ecological studies [11]. The Normalized Difference Vegetation Index (NDVI), a proxy for photosynthetically active vegetation, is strongly correlated with variables like Leaf Area Index, aboveground biomass, and percentage vegetation cover. NDVI is also widely used for phenological change studies [12][13][14][15]. Thus, trends and dynamics in NDVI can be used to monitor changes in vegetation cover, productivity, and health status at both large spatial and long temporal scales [16][17][18][19]. Moderate Resolution Imaging Spectroradiometer (MODIS) data have been considered state-of-the-art since the launch of the Terra platform in 1999 and dense temporal resolution NDVI products are widely used for vegetation dynamic research [20,21]. Several studies have been conducted in Southwest China using NDVI time series [22][23][24][25][26]. However, few studies analyzed the persistence of vegetation trends, which can provide information about the long-term memory of vegetation trends beyond the period of analysis and data availability. The long-term memory effects can be captured by the auto-covariance function that decays exponentially, with a spectral density that tends to infinity [27,28]. The Hurst exponent, based on auto-covariance, estimated via Rescaled Range Analysis (R/S analysis), is a parameter which can be applied to quantitatively detect the persistence of time series data of natural phenomena. The Hurst exponent has been widely used in the fields of hydrology [29], climatology [30], geology [31] and economics [32]. Only recently has it been applied in the time series detection of vegetation variations [33][34][35], though not yet in a complex terrain region like the karst regions of Southwest China.
Vegetation dynamics could be controlled by environmental variables, such as water availability, temperature, and incident radiation [18,36]. Topographic heterogeneity imposes environmental constraints on vegetation dynamics by producing complex substrates with variable structure, hydrology, and chemistry [37]. In particular, elevation determines the altitudinal zonality of soil [38,39] and slope controls the velocity of surface flows [40]. The aspect affects the direction of flows, insolation and the intensity of evaporation [41]. The surface curvature influences water migration and accumulation in landscapes by gravity [42]. A combination of these topographic factors (elevation, slope, aspect and surface curvature) determines the conditions for the growth and the distribution of vegetation [43,44].
To the authors' best knowledge, previous research did not include the combined effects caused by topographic variability in vegetation dynamic studies, and also the spatial pattern of future vegetation dynamics along a topographic gradient has not yet been explored.
In the context of monitoring ecological conservation projects in Southwest China, this paper aims to provide information on the persistence of vegetation trends in karst terrains of varying complexity. The specific objectives of this study are, to (1) analyze the spatio-temporal vegetation trends during 2001-2013 and classify the study area according to predicted future trends; and (2) identify terrain types suited for restoration but also quantify the uncertainty of future trends under different terrain conditions. The findings of this study are relevant for decision-making processes in vegetation restoration programs and the generic design of the framework allows for applications in other regions of complex landscapes around the world.

Materials and Methods
MODIS NDVI data are utilized to generate growing season NDVI (GSN) for each year from 2001 to 2013 to identify regions with the best prospects/high uncertainty for future vegetation restoration in Southwest China (Figure 1). Linear trend analysis and Hurst exponent are used to produce maps of vegetation trends and their persistence. Based on this, a future vegetation trend is generated. Elevation and slope are extracted from a digital elevation model (DEM) and a terrain niche index (TNI) map is calculated, which is classified into 7 separate classes. Two distribution indexes, namely Future Recovery Distribution Index (FRDI) and Future Uncertainty Distribution Index (FUDI), are proposed to explore the distribution of future vegetation dynamics and to identify areas where vegetation has the best prospects for restoration or where future trends are uncertain.

Materials and Methods
MODIS NDVI data are utilized to generate growing season NDVI (GSN) for each year from 2001 to 2013 to identify regions with the best prospects/high uncertainty for future vegetation restoration in Southwest China (Figure 1). Linear trend analysis and Hurst exponent are used to produce maps of vegetation trends and their persistence. Based on this, a future vegetation trend is generated. Elevation and slope are extracted from a digital elevation model (DEM) and a terrain niche index (TNI) map is calculated, which is classified into 7 separate classes. Two distribution indexes, namely Future Recovery Distribution Index (FRDI) and Future Uncertainty Distribution Index (FUDI), are proposed to explore the distribution of future vegetation dynamics and to identify areas where vegetation has the best prospects for restoration or where future trends are uncertain.

Study Area
The study area is the Guizhou and Guangxi provinces in Southwest China (20°54′-29°13′N and 103°36′-112°04′E). The area has a total coverage of 4.13 × 10 5 km 2 and a population of 87.84 million people. The elevation ranges from −8 to 2855 m.a.s.l., with maximum heights in the northwest and the minimum in the southeast (Figure 2b). The climate is a subtropical and tropical monsoon climate, with 1465 mm falling during the growing season from April to November. The precipitation from April to September accounts for 70%-85% of the annual total precipitation [45] and the mean growing season temperature is 25 °C. A Total of 45% and 15% of the bedrocks are pure carbonate and impure carbonate (karst region), respectively. The fragile soils of the karst areas are prone to rocky degradation if over-utilized (e.g., by livestock farming, Figure 2a). The bedrock for the rest of the region consists of clastic rocks (non-karst region) [46]. Governmental rocky desertification control projects implemented in karst regions can be divided into six main project regions [47] shown in Figure 2b. The main vegetation cover types are cultivated vegetation (29%), scrub (28%), grassland (17%), needle-leaf forest (16%), and broad-leaf forest (5%) (Figure 2c). Although the seasonal cycle of the main vegetation types differs, the general pattern is uniform, reflecting a growing season from April to November (Figure 2d).

Study Area
The study area is the Guizhou and Guangxi provinces in Southwest China (20˝54 1 -29˝13 1 N and 103˝36 1 -112˝04 1 E). The area has a total coverage of 4.13ˆ10 5 km 2 and a population of 87.84 million people. The elevation ranges from´8 to 2855 m.a.s.l., with maximum heights in the northwest and the minimum in the southeast (Figure 2b). The climate is a subtropical and tropical monsoon climate, with 1465 mm falling during the growing season from April to November. The precipitation from April to September accounts for 70%-85% of the annual total precipitation [45] and the mean growing season temperature is 25˝C. A Total of 45% and 15% of the bedrocks are pure carbonate and impure carbonate (karst region), respectively. The fragile soils of the karst areas are prone to rocky degradation if over-utilized (e.g., by livestock farming, Figure 2a). The bedrock for the rest of the region consists of clastic rocks (non-karst region) [46]. Governmental rocky desertification control projects implemented in karst regions can be divided into six main project regions [47] shown in Figure 2b. The main vegetation cover types are cultivated vegetation (29%), scrub (28%), grassland (17%), needle-leaf forest (16%), and broad-leaf forest (5%) (Figure 2c). Although the seasonal cycle of the main vegetation types differs, the general pattern is uniform, reflecting a growing season from April to November (Figure 2d

Data and Processing
NDVI data used in this study are MODIS (MOD13Q1) from 2001 to 2013 with a temporal resolution of 16 days. The images have a spatial resolution of 250 m and are retrieved from daily, atmospherically corrected surface reflectance. They have been composited using methods based on product quality assurance metrics to remove low quality pixels and the highest quality pixels were chosen for each 16 day composite [48]. For this study, we used the maximum value composite (MVC) method to select the higher value of two 16-day composite images and obtain the monthly NDVI. This further minimizes atmospheric and scan/angle effects, cloud contamination and high solar zenith angle effects [49]. Then the NDVI values from April to November were averaged to obtain the GSN for each year from 2001 to 2013 (see Figure 2d).

Data and Processing
NDVI data used in this study are MODIS (MOD13Q1) from 2001 to 2013 with a temporal resolution of 16 days. The images have a spatial resolution of 250 m and are retrieved from daily, atmospherically corrected surface reflectance. They have been composited using methods based on product quality assurance metrics to remove low quality pixels and the highest quality pixels were chosen for each 16 day composite [48]. For this study, we used the maximum value composite (MVC) method to select the higher value of two 16-day composite images and obtain the monthly NDVI. This further minimizes atmospheric and scan/angle effects, cloud contamination and high solar zenith angle effects [49]. Then the NDVI values from April to November were averaged to obtain the GSN for each year from 2001 to 2013 (see Figure 2d). A SRTM DEM with 90 m resolution was downloaded from the International Scientific Data Service Platform [50] and used as the source for elevation and slope. The DEM was resampled to 250 m spatial resolution using a nearest neighbor resampling algorithm to match the MODIS NDVI data.
Monthly temperature and rainfall data for 53 weather stations from 2001 to 2013 were obtained from the China Meteorological Data Sharing Service System [51]. The human population data in 2005 and 2010 were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) [52].

Methods
A linear regression between annual GSN and time was applied to derive the slope as an indication of the direction and magnitude of linear trends in the time series. Whereas a positive slope value indicates a positive trend in vegetation (e.g., successful conservation projects), a negative slope value refers to negative vegetation trends which can be an indicator of degradation (e.g., erosion due to anthropogenic overuse). The trends are considered statistically significant if the p value of the regression is smaller than 0.05. If the p value is larger 0.05, no significant change is detected and the pixel is considered to be stable. We thus categorize GSN trends as three types: increasing (slope > 0 and p < 0.05), decreasing (slope < 0 and p < 0.05) and stable (p ě 0.05).
The rescaled range (R/S) analysis developed by Hurst is a method used to estimate auto-correlation properties of time series [53]. We apply the Hurst exponent (H), estimated by the R/S analysis, as a measure of the long-term memory in our GSN time series (annual data). The main calculation steps are [54]: 1. To divide the time series tξ pτqu(τ = 1, 2 . . . , n) into τ sub series x ptq, and for each sub series t = 1, . . . , τ. 2. To define the mean sequence of the time series, x ptq , τ " 1, 2 . . . , n 3. To calculate the cumulative deviation, To create the range sequence, To create the standard deviation sequence, The values of the H range between 0 and 1. A value of 0.5 indicates a true random walk. In a random walk there is no correlation between any element and a future element. An H value between 0.5 and 1 (H, 0.5 < H < 1) indicates "persistent behavior" (i.e., a positive autocorrelation). The closer the H values are to 1, the more persistent is the time series. An H value below 0.5 (H, 0 < H < 0.5) means a time series with "anti-persistent behavior" (or negative autocorrelation). The closer the H values are to 0, the more anti-persistent is the time series [55]. The H value was calculated for each pixel to test the persistence of the GSN time series.
Based on the GSN trends and the H value, we categorized the supposed future vegetation trends into six types (Table 1). When a pixel showed a positive trend and its H value was higher than 0.5 (i.e., the positive trend is likely to continue), the vegetation in this pixel is supposed to have an overall positive development (PD) after the time series ends in 2013. If the trend is positive but its H value was less than 0.5 (the positive trend is not likely to continue), the vegetation in this pixel has an anti-persistent positive development (APD) after 2013. When a pixel showed a negative trend and its H value was higher than 0.5 (the negative trend is likely to continue), the vegetation in this pixel is supposed to experience a negative development (ND) after 2013. If the trend is negative and its H value was less than 0.5, (i.e., the negative trend is not likely to continue), the vegetation in this pixel is characterized by an anti-persistent negative development (AND) after 2013. When a pixel was detected to have no significant trend and its H value was higher than 0.5, the vegetation in this pixel is supposed to have a stable and steady development (SSD) after 2013. If its H value was less than 0.5, the future vegetation trend in this pixel is uncertain, which means the trend after 2013 is an undetermined development (UD).  The terrain niche index (TNI) was utilized to characterize the topographic variation in Southwest China [56]. It was calculated as: where e and E are the elevation of the pixel and the average elevation of the study area respectively, whereas s and S signify the slope of the pixel and the average slope of the study area. Large TNI values correspond to higher elevation and larger slope angles, being typical for karst plateaus and gorges. In contrast, smaller TNI values indicate lower elevation and smaller slope angles. Medium TNI values are found in areas of a higher elevation but small slope angle, or low elevation but with larger slope angles, or moderate elevation and slope angle. In this study the TNI values were categorized into 7 classes with an interval of 0.2.
To characterize the status of the vegetation for each terrain interval beyond 2013, the future vegetation trend classes shown in Table 1 were used as reference. All pixels within each of the 7 terrain intervals were classified according to one of the future trend classes and the overall sum of the classes was calculated (where ND has´1, others have the value 1) for each terrain interval. To account for different area sizes of the terrain intervals and capture supposed future change for each specific interval, we propose two distribution indexes, FRDI (identifying areas where vegetation has the best prospects for restoration) and FUDI (showing the distribution of the uncertainty). These two comprehensive indexes are calculated as follows: where i refers to the 7 terrain intervals; e ij and E j (j = 1 (PD), 2 (SSD), 3 (ND), 4 (UT), 5 (APD) and 6 (AND)) are the area of the six future vegetation trends (Table 1) in each terrain interval i and the study area; and a i and A refer to the area of the terrain interval i and the study area, respectively. FRDI i and FUDI i are standardized and dimensionless indexes. If FRDI i > 1, the corresponding terrain interval i is supposed to be more favorable for vegetation restoration in the future than other terrain intervals. The larger the FRDI i is, the better the future prospect changes. If the FRDI i is lower than 1, the terrain interval i is supposed to have less favorable prospects. If FUDI i > 1, trend predictions in the terrain interval i have a higher uncertainty than other terrain intervals. A lower FUDI i reduces prediction uncertainty and a FUDI i below 1 implies more certain prospects.

Vegetation Patterns and Dynamics during 2001-2013
The growing season vegetation in the study area has an average GSN of 0.72 for the period 2001-2013 with a high spatial heterogeneity between the land-cover classes (Figure 3a). For the Guangxi Province, the average GSN is 0.74, which is slightly higher than in the Guizhou Province (0.69). The vegetation type of the needle-leaf and broad-leaf mixed forest class has the highest mean GSN (0.78), followed by broad-leaf forest (0.77), needle-leaf forest (0.74), grassland (0.73), scrub (0.72), and cultivated vegetation (0.68). The average GSN in the non-project region (0.74) is slightly higher than in the project regions (0.71). The regions with a GSN higher than 0.8 are mainly located in Guangxi Province, such as Ninming, Baise, Tianlin, Zhaoping, Cangwu, Tiane, Yongfu County, dominated by broad-leaf forest and needle-leaf forest. The areas with a GSN less than 0.5 are associated with scrubs, cultivated vegetation and built-up areas. These regions are close to roads and railways and correspond to county seats, especially in large cities, such as Nanning, Liuzhou and Guiyang. At the county scale, we found the average GSN of the past 13 years had a significant positive relationship with the distance to roads (p < 0.01) and a significant negative relationship with the population density change from 2005 to 2010 (p = 0.016). Moreover, areas of low GSN are found in mountainous regions dominated by ethnic minorities, such as Weining County populated by the Yi nationality and the Fuchuan County inhabited by the Yao nationality. At the county scale, we found significant positive correlations (p < 0.01) between GSN and growing season climatic factors (temperature and precipitation). GSN was more strongly and positively correlated with temperature (R 2 = 0.25) than with precipitation (R 2 = 0. 11

Vegetation Patterns and Dynamics during 2001-2013
The growing season vegetation in the study area has an average GSN of 0.72 for the period 2001-2013 with a high spatial heterogeneity between the land-cover classes (Figure 3a). For the Guangxi Province, the average GSN is 0.74, which is slightly higher than in the Guizhou Province (0.69). The vegetation type of the needle-leaf and broad-leaf mixed forest class has the highest mean GSN (0.   Figure 3b shows the trends in annual GSN with distinct spatial differences. An uptrend (slope > 0) is found in 67% of the study area, with 23% of these trends being significant (p < 0.05). These pixels are mainly in Qinzhou and Beihai of Guangxi and Bijie of Guizhou, where the dominant vegetation types are scrubs and cultivated vegetation (i.e., agro-forestry, such as eucalyptus plantation and oil-tea camellia plantation). Furthermore, 33% of the pixels are detected to have a downtrend (slope < 0), with 10% being significant. These pixels are mainly located in Baise and Liuzhou, Guangxi. In addition, the dominating vegetation types of these areas are cultivated vegetation, scrubs and  Figure 3b shows the trends in annual GSN with distinct spatial differences. An uptrend (slope > 0) is found in 67% of the study area, with 23% of these trends being significant (p < 0.05). These pixels are mainly in Qinzhou and Beihai of Guangxi and Bijie of Guizhou, where the dominant vegetation types are scrubs and cultivated vegetation (i.e., agro-forestry, such as eucalyptus plantation and oil-tea camellia plantation). Furthermore, 33% of the pixels are detected to have a downtrend (slope < 0), with 10% being significant. These pixels are mainly located in Baise and Liuzhou, Guangxi. In addition, the dominating vegetation types of these areas are cultivated vegetation, scrubs and grasslands. Consequently, three vegetation trend types, increasing, decreasing and stable (not significant; p ě 0.05) accounted for 15%, 3% and 82%, respectively.

Future Vegetation Trends after 2013
Values of the H range from 0.02 to 0.86 (Figure 4a). The mean H value of the whole study area is 0.55, indicating an overall moderately persistent vegetation trend in the two provinces in Southwest China, suggesting that future vegetation trends are likely to be similar to those of the past 13 years. In particular, the H values between most of the vegetation cover types are relatively similar: broad-leaf forest has a mean H value (0.552), followed by needle-leaf forest (0.551), grassland (0.548), scrub (0.547), and cultivated vegetation (0.546). Only the class needle-leaf and broad-leaf mixed forest is below 0.5 (0.473). In 69% of the study area the H value is above 0.5. Areas with a low persistence were located in the west and the south of the study area. In the west, the regions with low H values were mostly corresponding to the road network, while in the south, the regions with low persistence are close to county seats (Figure 4a).
The vegetation dynamics and the H values were combined into different classes (Table 1) to indicate future vegetation trends (Figure 4b). The future vegetation trends in most parts of the study area (57%) are stable and steady, and particularly vegetation belonging to scrub and cultivated vegetation (Figure 1c) is detected to have no significant trend in the past 13 years and this state will be persistent (H > 0.5). However, 10% of the study area is expected to have a positive development in the future (vegetation had a persistent positive trend). These areas consist primarily of scrubs and are mainly distributed in the north and northwest Guizhou and in southern Guangxi. Areas with a negative development (where vegetation was detected to have a persistent decreasing trend), cover 3% of the study area and the main land cover of these areas is cultivated vegetation. Vegetation showing no trend and no persistence occupy 24% of the study area, where the main land cover types are scrubs and cultivated vegetation. Regions with an anti-persistent negative development cover 1% of the study area. The main land cover of these regions is cultivated vegetation. Areas with an anti-persistent positive development cover 5% of the study area, with scrub being the main land cover class. grasslands. Consequently, three vegetation trend types, increasing, decreasing and stable (not significant; p ≥ 0.05) accounted for 15%, 3% and 82%, respectively.

Future Vegetation Trends after 2013
Values of the H range from 0.02 to 0.86 (Figure 4a). The mean H value of the whole study area is 0.55, indicating an overall moderately persistent vegetation trend in the two provinces in Southwest China, suggesting that future vegetation trends are likely to be similar to those of the past 13 years. In particular, the H values between most of the vegetation cover types are relatively similar: broad-leaf forest has a mean H value (0.552), followed by needle-leaf forest (0.551), grassland (0.548), scrub (0.547), and cultivated vegetation (0.546). Only the class needle-leaf and broad-leaf mixed forest is below 0.5 (0.473). In 69% of the study area the H value is above 0.5. Areas with a low persistence were located in the west and the south of the study area. In the west, the regions with low H values were mostly corresponding to the road network, while in the south, the regions with low persistence are close to county seats (Figure 4a).
The vegetation dynamics and the H values were combined into different classes (Table 1) to indicate future vegetation trends (Figure 4b). The future vegetation trends in most parts of the study area (57%) are stable and steady, and particularly vegetation belonging to scrub and cultivated vegetation (Figure 1c) is detected to have no significant trend in the past 13 years and this state will be persistent (H > 0.5). However, 10% of the study area is expected to have a positive development in the future (vegetation had a persistent positive trend). These areas consist primarily of scrubs and are mainly distributed in the north and northwest Guizhou and in southern Guangxi. Areas with a negative development (where vegetation was detected to have a persistent decreasing trend), cover 3% of the study area and the main land cover of these areas is cultivated vegetation. Vegetation showing no trend and no persistence occupy 24% of the study area, where the main land cover types are scrubs and cultivated vegetation. Regions with an anti-persistent negative development cover 1% of the study area. The main land cover of these regions is cultivated vegetation. Areas with an anti-persistent positive development cover 5% of the study area, with scrub being the main land cover class.

Vegetation Trends and Future Vegetation Trends in Different Project Regions
Some differences in vegetation trends and future vegetation trends can be observed ( Figure 5) in different regions (Figure 2b). Apart from the karst gorge region (IV) and non-project region (VII), vegetation in 80% of other project regions is rather stable. Vegetation in 32% of the karst gorge region

Vegetation Trends and Future Vegetation Trends in Different Project Regions
Some differences in vegetation trends and future vegetation trends can be observed ( Figure 5) in different regions (Figure 2b). Apart from the karst gorge region (IV) and non-project region (VII), vegetation in 80% of other project regions is rather stable. Vegetation in 32% of the karst gorge region (IV) (having the lowest average GSN with 0.62) shows a significant increasing trend, indicating an improvement of the vegetation in this region. Overall, only a small percentage of the pixels in the project regions have a decreasing trend. With 5% in the peak forest plain region (II), 4% in the peak cluster depression region (I) and 3% in the non-project region (VII), these regions are the only ones with a noticeable share of pixels characterized by significant decreasing trends. improvement of the vegetation in this region. Overall, only a small percentage of the pixels in the project regions have a decreasing trend. With 5% in the peak forest plain region (II), 4% in the peak cluster depression region (I) and 3% in the non-project region (VII), these regions are the only ones with a noticeable share of pixels characterized by significant decreasing trends. The future vegetation trends with a stable and steady development after 2013 cover 58% and 54% of all project region and non-project region pixels, respectively. In particular, the karst trough valley region (V) and the peak forest plain region (II) account for 66% and 63% respectively, having the highest share of steady trends among the project regions. Overall, 19% of the karst gorge region (IV) and 11% of the karst plateau region (III) are predicted to have a positive development after 2013. However, 37% of the overall positive trends are anti-persistent. In particular, the karst gorge region (IV) needs attention, with 13% of the vegetation showing a positive but anti-persistent development. A total of 25% of the pixels in the project regions have an undetermined trend. More specifically, in 34% of the karst basin region (VI) and 28% of the karst gorge region (IV), the future vegetation trends cannot be determined. The highest share of negative development was found in the peak forest plain region (II) (4%). In general, the share of anti-persistent negative development is small for all project regions (0.5%), indicating some efforts in preventing degradation.

Distribution of Future Vegetation Trends for Different Terrain Conditions
The TNI values in the study area (ranging from −0.0006 to 1.4) show a pattern of decrease from the northwest to the southeast (Figure 6a), which in general is consistent with the distribution of elevation and slope. The karst gorge region (IV) has the highest TNI (0.85), followed by the karst basin region (VI) (0.84), the karst plateau region (III) (0.70), the karst trough valley region (V) (0.67), the peak cluster depression region (I) (0.54), and the peak forest plain region (II) (0.40). The distribution of future vegetation trends varies between TNI groups (i.e., with the terrain) (Figure 6b). The largest shares of negative and positive development are found in regions with low elevation and slope, which are also the areas characterized by the highest anthropogenic influence. The negative development decreases and the undetermined development increases as the TNI increases. This indicates a decline of anthropogenic impact in rough terrain and a more difficult prediction of future trends.
The distribution of the FRDI and FUDI along the seven terrain intervals (Figure 6c) shows that FRDI (representing the areas with best future prospects) is low in areas with a low TNI but is larger than 1 when the TNI is in the range between 0.4 and 1 (areas with a moderate slope and elevation). The uncertainty for future predictions (FUDI) is high in areas of low TNI, but low in areas with a TNI between 0.4 and 1, being in line with the FRDI. Thus, these terrain niche intervals are identified as regions with the best prospects for vegetation restoration in the future. The FUDI increases again The future vegetation trends with a stable and steady development after 2013 cover 58% and 54% of all project region and non-project region pixels, respectively. In particular, the karst trough valley region (V) and the peak forest plain region (II) account for 66% and 63% respectively, having the highest share of steady trends among the project regions. Overall, 19% of the karst gorge region (IV) and 11% of the karst plateau region (III) are predicted to have a positive development after 2013. However, 37% of the overall positive trends are anti-persistent. In particular, the karst gorge region (IV) needs attention, with 13% of the vegetation showing a positive but anti-persistent development. A total of 25% of the pixels in the project regions have an undetermined trend. More specifically, in 34% of the karst basin region (VI) and 28% of the karst gorge region (IV), the future vegetation trends cannot be determined. The highest share of negative development was found in the peak forest plain region (II) (4%). In general, the share of anti-persistent negative development is small for all project regions (0.5%), indicating some efforts in preventing degradation.

Distribution of Future Vegetation Trends for Different Terrain Conditions
The TNI values in the study area (ranging from´0.0006 to 1.4) show a pattern of decrease from the northwest to the southeast (Figure 6a), which in general is consistent with the distribution of elevation and slope. The karst gorge region (IV) has the highest TNI (0.85), followed by the karst basin region (VI) (0.84), the karst plateau region (III) (0.70), the karst trough valley region (V) (0.67), the peak cluster depression region (I) (0.54), and the peak forest plain region (II) (0.40). The distribution of future vegetation trends varies between TNI groups (i.e., with the terrain) (Figure 6b). The largest shares of negative and positive development are found in regions with low elevation and slope, which are also the areas characterized by the highest anthropogenic influence. The negative development decreases and the undetermined development increases as the TNI increases. This indicates a decline of anthropogenic impact in rough terrain and a more difficult prediction of future trends.
The distribution of the FRDI and FUDI along the seven terrain intervals (Figure 6c) shows that FRDI (representing the areas with best future prospects) is low in areas with a low TNI but is larger than 1 when the TNI is in the range between 0.4 and 1 (areas with a moderate slope and elevation). The uncertainty for future predictions (FUDI) is high in areas of low TNI, but low in areas with a TNI between 0.4 and 1, being in line with the FRDI. Thus, these terrain niche intervals are identified as regions with the best prospects for vegetation restoration in the future. The FUDI increases again with higher TNI, indicating that regions with a small slope and low elevation or a large slope and high elevation contain the highest uncertainty in future prospects.
Most pixels with TNI between 0.4 and 1 are found in the peak cluster depression region (I) (36%), the karst plateau region (III) (19%) and the karst trough valley region (V) (16%) (Figure 7a). Therefore, the future vegetation restoration has the best prospects in these project regions. The areas with a TNI greater than 0.8 are mainly distributed in the peak cluster depression region (I) (25%), the karst plateau region (III) (19%) and the karst gorge region (IV) (17%). The future vegetation trends in these project regions have a high uncertainty. However, based on the ratio of each terrain interval in every project region (Figure 7b), we found the largest share of TNI between 0.4 and 1 in the karst plateau region (III) (95%) and the karst trough valley region (V) (89%). Thus, the share of pixels prospecting from future vegetation restoration in the karst plateau region (III) and the karst trough valley region (V) will be greater than in other regions. The largest shares of TNI greater than 0.8 are found in the karst gorge region (IV) (62%) and the karst basin region (VI) (61%), while the largest shares of TNI smaller than 0.2 are found in the non-project region (VII) (35%) and the peak forest plain region (II) (33%). Both high and low TNI values indicate a high future vegetation uncertainty in these regions. Most pixels with TNI between 0.4 and 1 are found in the peak cluster depression region (I) (36%), the karst plateau region (III) (19%) and the karst trough valley region (V) (16%) (Figure 7a). Therefore, the future vegetation restoration has the best prospects in these project regions. The areas with a TNI greater than 0.8 are mainly distributed in the peak cluster depression region (I) (25%), the karst plateau region (III) (19%) and the karst gorge region (IV) (17%). The future vegetation trends in these project regions have a high uncertainty. However, based on the ratio of each terrain interval in every project region (Figure 7b), we found the largest share of TNI between 0.4 and 1 in the karst plateau region (III) (95%) and the karst trough valley region (V) (89%). Thus, the share of pixels prospecting from future vegetation restoration in the karst plateau region (III) and the karst trough valley region (V) will be greater than in other regions. The largest shares of TNI greater than 0.8 are found in the karst gorge region (IV) (62%) and the karst basin region (VI) (61%), while the largest shares of TNI smaller than 0.2 are found in the non-project region (VII) (35%) and the peak forest plain region (II) (33%). Both high and low TNI values indicate a high future vegetation uncertainty in these regions.

Validation of the Method
The time series was split into two parts (2001-2007 and 2008-2013) and pixels with a significant (p < 0.05) trend for both periods were chosen for validating the future trend prediction methodology ( Table 2). In spite of the short time periods, the prediction performs well for most classes-only the prediction for negative persistent trends does not apply. Table 2. For validation, the time series was split into two periods and pixels with significant (p < 0.05) trends in both periods were analyzed. Due to the short periods, the numbers of pixels meeting these assumptions are limited (shown in brackets in the first column). The mean slopes (GSN year´1) for these pixels (both periods) are shown. The predictions are based on the first period (2001 to 2007).

Drivers of Vegetation Conditions and Dynamics
Many studies have suggested that climate change is one of the main drivers of trends in vegetation productivity in China [57,58]. However, in our study the role of average growing season temperature and precipitation is limited since the coefficients of determination of both relationships are rather low (yet statistically significant), confirming the results of Hou et al. [59]. Furthermore, the relationship between GSN and precipitation is weaker than for temperature, which is consistent with the study by Wang et al. [60]. This might be due to Southwest China being located in the subtropical monsoon climate zone, with rich precipitation and moderate temperature, and thereby slight changes in precipitation will not heavily affect vegetation growth. On the other hand, the above and below-ground hydrological structures of the karst region leads to a significant drain and runoff [61], which suggests that precipitation cannot be effectively utilized for vegetation growth.
The accumulative afforested area by aerial seeding and forest plantation in Southwest China had increased quickly from 65,744 ha in 2000 to 1,709,366 ha in 2011. Furthermore, the overall deteriorated trends of rocky desertification had been transformed from net increase to net decrease since 2011 in karst regions, Southwest China [70]. With the implementation of the Karst Rocky Desertification Comprehensive Control and Restoration Project, the soil fertility and surface vegetation cover rate in karst regions have been improved [71][72][73]. Partly due to the ecological projects in the past 15 years, most vegetation communities in the present study area have achieved a stable growing stage with a relatively high GSN and show thus no significant trend (Figure 3b). However, some regions with a significant decrease in vegetation have been detected. These areas are usually around county seats or in proximity to roads, identifying rapid urbanization and road construction to lead to the degradation of vegetation. Generally, the average cover and successional rate of vegetation are lower in karst regions as compared to non-karst regions [74]. In our study, the average GSN in project regions (karst region) is still lower than in non-project region (non-karst region). However, the annual increase of GSN from 2001 to 2013 in project regions (0.0024 GSN year´1) is higher than in non-project region (0.0023 GSN year´1), indicating that vegetation cover grows more quickly under the implementation of ERPs.

Influence of Terrain Conditions on Vegetation Distribution
A framework for governmental decision makers is proposed to highlight regions that need more attention when implementing ecological projects. The FRDI and FUDI indexes, based on the future vegetation trends and the TNI, were utilized to detect specific areas where vegetation had the best prospects to be restored or the highest uncertainty in the future. The results showed that regions with a TNI between 0.4 and 1 were the most advantageous intervals for vegetation restoration in the future. In the past 13 years, the average growing season temperature and precipitation dropped as the TNI increased. In regions with TNI less than 0.4, the average growing season temperature and precipitation are 28˝C and 1721 mm respectively, whereas in regions with TNI between 0.4 and 1, averages are 25˝C and 1355 mm, and in regions with TNI greater than 1, 20˝C and 1218 mm. In the region studied, areas with a low TNI (less than 0.4) have the best combined conditions of air temperature and water favoring vegetation growth. In addition, the Returning Farmland to Forest Program has been implemented in regions with slope greater than 25˝, which supports vegetation restoration. Therefore, vegetation in these regions is likely to have reached a stable condition, and the room for further improvement is less than in other regions. The conditions of air temperature and water in regions with a TNI between 0.4 and 1 are less ideal compared to areas with a lower TNI but better than for areas with a higher TNI. Therefore, vegetation in these regions has a higher chance of improvement under the implementation of projects than in regions with a higher TNI. According to the distribution of FUDI in areas with different TNIs, vegetation development in regions with a high TNI (greater than 0.8) has a higher uncertainty. These regions with high elevation and large slope angles can induce soil erosion under heavy rainfall. Besides, vegetation in these areas suffers more easily under droughts since the steep terrain condition is not suitable for water conservation [75]. Thus, vegetation in these regions is likely to be more strongly influenced by climatic conditions than other regions and fluctuate along with climate variations. Areas with TNI smaller 0.2 are found in the non-project region (VII) and the peak forest plain region (II). These areas are well suited for vegetation growth and are characterized by high anthropogenic influence. The anthropogenic disturbances lead to a generally uncertain (reflected in a high FUDI) and negative vegetation development.

Limitations and Future Research Directions
For validation, the time series was divided into two periods. Results showed that the prediction performs well for most classes except the prediction for negative persistent trends. This could be interpreted as an indication of the successful implementation of anthropogenic projects preventing further degradation. However, both periods are rather short due to the length of the full time series based on the use of MODIS data (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). Thus, the validation is based on trends and predictions characterized by high uncertainty and a longer time series is needed for a more robust validation. Moreover, the TNI was utilized to characterize the terrain conditions and the two indexes (FRDI and FUDI) were an effective approach for identifying the most advantageous intervals for vegetation restoration in the future and the corresponding prediction uncertainty. However, we were not able to quantify the specific elevation and slope of these regions. Thus, more work is needed to analyze the relationship between the vegetation variations and topographical conditions and to explore the exact ranges of elevation and slope suited for specific purposes. The variations of the vegetation trends were affected by climate and human activities, especially national/local policies. We explored the driving forces associated with vegetation variations but did not statistically distinguish the effects caused by climate and human activities on vegetation dynamics, which facilitates a quantification of the efficiency of human-induced effort on vegetation restoration. with an increasing trend is greater than that of decreasing trends. The karst gorge region has the largest share of increasing trends amongst the 6 project regions. The peak forest plain region and the peak cluster depression region are found to have a noticeable share of pixel showing a significant decreasing trend. (3) The vegetation trends in Southwest China are rather persistent. For the whole study area, the future vegetation with a stable and steady development is dominant, especially in the karst trough valley and the peak forest plain. The percentage of vegetation with positive development is greater than that of negative development. The karst gorge has the highest share of pixels showing positive development and anti-persistent development. The karst peak forest plain region has the highest share of pixels with negative development. The karst basin has the highest share of pixels characterized by an undetermined trend. (4) Most negative future development is found in areas of high anthropogenic influence, as compared to areas of rough terrain. (5) The regions with terrain niche index (TNI) between 0.4 and 1 have the best prospects for vegetation restoration in the future (especially in the karst plateau region and the karst trough valley region). Moreover, regions with TNI greater than 1 show the highest uncertainty for future prospects (especially in the karst gorge region and the karst basin region). It is here suggested that governmental restoration projects should duly pay attention to these regions to prevent vegetation degradation. (6) We have shown that the framework of the present analysis works well for assessing restoration prospects in the study area. Due to its generic design, the method used is expected to be applicable for studies located in other areas of complex landscapes in the world, exploring future vegetation trends.