Proﬁling Human-Induced Vegetation Change in the Horqin Sandy Land of China Using Time Series Datasets

: Discriminating the signiﬁcant human-induced vegetation changes over the past 15 years could help local governments review the effects of eco-programs and develop sustainable land use policies in arid/semi-arid ecosystems. We used the residual trends method (RESTREND) to estimate the human-induced and climate-induced vegetation changes. Two typical regions in the Horqin Sandy Land of China were selected as study areas. We ﬁrst detected vegetation dynamics between 2000–2014 using Sen’s slope estimation and the Mann–Kendall test detection method (SMK) based on the Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) time series, then used RESTREND to proﬁle human modiﬁcations in areas of signiﬁcant vegetation change. RESTREND was optimized using statistical and trajectory analysis to automatically identify ﬂexible spatially homogeneous neighborhoods, which were essential for determining the reference areas. The results indicated the following. (1) Obvious vegetation increases happened in both regions, but Naiman (64.1%) increased more than Ar Horqin (16.8%). (2) Climate and human drivers both contributed to signiﬁcant changes. The two factors contributed equally to vegetation change in Ar Horqin, while human drivers contributed more in Naiman. (3) Human factors had a stronger inﬂuence on ecosystems, and were more responsible for vegetation decreases in both regions. Further evidences showed that the primary human drivers varied in regions. Grassland eco-management was the key driver in Ar Horqin, while farming was the key factor for vegetation change in Naiman.


Introduction
Discriminating between the driving forces of vegetation change is necessary in order to understand the interacting mechanisms between ecosystems and external drivers, and help environmental managers make effective decisions that maintain the sustainable development of ecosystems [1,2]. Complex external driving forces influence the processes of vegetation change, of which climate change and human activities are two main factors within decadal time scales [3][4][5]. Climate variations, such as increases or decrease in precipitation and temperature would result in a corresponding trend of greening and browning in vegetation cover [3,5,6]. Various human activities may also influence the vegetation growth by controlling or modifying ecosystem composition and distribution [3]. Worldwide, urbanization usually causes local vegetation to decrease [7][8][9], while eco-managements such as encouraging afforestation, grazing prohibition, and rest-rotation grazing usually cause regional-dependent task. Furthermore, heterogeneous vegetation patterns usually exist within a regular and limited distance, which weakens the applicability of this strategy. An irregular homogeneous neighborhood is a better choice. Prince [32] used a land capability class map, which was defined by a k-prototypes clustering technique using maps of rainfall, soils, and land use, to obtain a homogeneous spatial neighborhood. However, this method still faces two obvious problems. (1) It is hard to control the number of classes in the k-prototypes clustering algorithm, which results in uncertainty in land capability classes. (2) The type of land use within a time series might change from year to year, so it is not reasonable to use a stationary land use map to define a homogeneous area. However, it is difficult to obtain a land use map each year to update the land capability class map.
To overcome the problems of these two methods, this paper optimized the strategy of determining homogeneous spatial neighborhoods in RESTREND by statistical analysis and trajectory analysis to automatically define flexible spatially homogeneous neighborhoods. Based on this, we first used a trend analysis method to find significant vegetation change from 2000 to 2014 in the study areas, and then used the optimized RESTREND in areas with significant changes to identify significant human-induced changes during the past 15 years. By choosing two typical regions in Horqin Sandy Land as study areas, this paper aimed to discover different human-climate-vegetation change mechanisms in arid and semi-arid ecosystems with different levels of human influence.

Study Area
The Horqin Sandy Land spans from 41 • 41 N to 46 • 12 N latitude and 117 • 49 E to 123 • 42 E longitude, and is located in the southeast of the Inner Mongolian Autonomous Region of China ( Figure 1). It is a typical transitional zone between the Inner Mongolian Plateau and the Northeast China Plain, with a temperate sub-humid and semi-arid monsoon climate [11]. The annual precipitation varies between 340-450 mm. The annual average temperature is about 5.8-6.4 • C, and the annual average wind velocity is up to 3.5-4.5 m/s. The landscape of this area is mainly consisted of sand dunes alternating with gently undulating meadows. The soil of this area is sandy, light, loose, and especially susceptible to wind erosion. The Horqin Sandy Land is a typical transitional zone between semi-pastoral areas located in the south of the study area, and pastoral areas located in the north of the study area.
We used two regions, i.e., Ar Horqin and Naiman, with different levels of human influence in the Horqin Sandy Land as our study areas to discover different human-climate-vegetation change mechanisms. The Ar Horqin region is a typical pastoral area that is in the northwest of the Horqin Sandy Land, with 71.0% of the land covered by grassland. The Naiman region is in the southeast of the Horqin Sandy Land. It is a typical agricultural area with the majority (46.5%) covered by farmland ( Figure 1). Traditionally, people in the pastoral area prefer to feed sheep and cows for a living, while people in the semi-pastoral area live on farming. According to the Inner Mongolia Statistics Bureau in 2013, the population density in Naiman is 53 people/km 2 , which is much more than in Ar Horqin (21 people/km 2 ). It is believed that human activities in Naiman are much more intense than in Ar Horqin. Over the past two decades, the ecosystems in these two regions have changed markedly because of frequent human activities [34][35][36][37]. The Chinese government has implemented several region-wide programs in the two locations; e.g., the West Development Strategies (introduced in 1999) and the Grain for Green Program (introduced in 1999). The former strategy promotes economic development in western China, and the latter encourages farmers to tree plant in order to restore local ecological environment. Additionally, three grassland conservation programs have also been implemented in the study areas. The Beijing-Tianjin Sandstorm Source Controlling Program (introduced in 2001) and the Ecological Subsidy and Award System (introduced in 2011) were also implemented in the Ar Horqin region. The Grazing Withdrawal Program (introduced in 2003) and the Ecological Subsidy and Award System (introduced in 2011) were implemented in Naiman [38].
Discriminating the significant human-induced vegetation changes over the last 15 years could help local governments review the effects of the eco-programs and develop sustainable land use policies in the future.

Data and Data Preprocessing
This study selected NDVI as the original data, because the index is very sensitive to different vegetation covers in arid/semi-arid areas. It has a strong capability for monitoring vegetation changes, land cover classifications, and the change analysis of land cover in arid and semi-arid areas [21,22,39]. The NDVI, which is a proxy for NPP, is the most commonly used indicator in RESTREND [1,[40][41][42]. Time series Moderate Resolution Imaging Spectroradiometer (MODIS) 13Q1 NDVI datasets were downloaded from the MODIS website (http://modis.gsfc.nasa.gov/). The datasets covers 342 images from 2000 to 2014 with the same temporal resolution (16 days) and spatial resolution (250 m).
The original dataset had a hierarchical data format, within which the layer 1 and the layer 12 reflected the NDVI value and the reliability and quality of the data for each pixel, respectively.
We first removed noise from clipped and re-projected NDVI datasets using a Savitzky-Golay filter [43]. Then, we replaced high ranking pixels (i.e., ranking 0 or 1) with the original pixel value [44] and masked non-vegetated areas (e.g., deserts, water bodies, and other areas with mean annual NDVI <0.1) [45]. Here, we changed the NDVI values of high-ranking pixels, because these two types of pixels were confirmed as high reliability pixels. In layer 12 of dataset, a high ranking of '0 for pixels indicates high quality. The ranking of '1 indicates that their value was useful and could be used with reference to other QA information. Finally, the time series ∑ NDVI trajectory, i.e., the annual accumulated NDVIs (AA-NDVIs), was aggregated from Julian day 145 to Julian day 273 in the growing season of each year. The time series ∑ NDV I dataset was then produced so that the impact of frequent and strong interannual variations could be reduced.
Although the standard RESTREND only used rainfall as the climate factor, many studies suggest that other climate factors such as temperature might be key controls on vegetation growth in some regions [19]. Therefore, it was necessary to determine the principal climate factors and the suitable climate variables to build the climate-vegetation model [33,46,47]. In this paper, precipitation, temperature, and wind speed were selected as three representative climate factors to pre-analyze the best climate factors to use in the study area. Datasets for monthly average precipitation (P), monthly average temperature (T), and monthly maximum wind speed (W) were produced based on time series data (2000-2014) from 46 stations acquired from the website of the China Meteorological Data Service Center (http://data.cma.cn/). Due to the lag effect of vegetation growth, 45 indices (Table 1) with different temporal combinations were calculated for each meteorological dataset (P/T/W); then, the indices in each station were interpolated by the inverse distance weighting method to obtain 135 climate raster grids with the same resolution as the NDVI dataset. The 1:1,000,000 ecological map for Inner Mongolia, produced by the Chinese Academy of Sciences (2001), was used to identify ecological zones. An ecological zone is determined by the climate, soil, landform, and vegetation. The combination of land cover types in an ecological zone is unique, and the same land cover type in an ecological zone indicates a similar combination of plants and their growing processes.
A Level 1 Terrain-corrected (L1T) Landsat Thematic Mapper and Enhanced Thematic Mapper Plus (TM/ETM+) time series image dataset covering the study areas in summer of 2000,2005,2009, and 2014 was downloaded from the United States Geological Survey website (http://earthexplorer.usgs.gov/). The datasets had been processed for systematic radiometric calibration and geometric correction by incorporating ground control points while employing a digital elevation model for topographic correction. Therefore, it could be used to visually validate the results of human-induced vegetation change. The 2010 land cover data for the study area was obtained from the China Environmental Protection Agency and the Chinese Academy of Sciences. It was produced by object-based classification using charge coupled device camera data from the Chinese Huanjing-1/A/B/C satellite. The classification of land cover in the data was based on first-level categories, and the overall accuracy was up to 85% [48]. It could indicate farmland, forest, grassland, wetland, sand land/unused land, and urban construction areas, so that can be used to further target specific driven factors. In addition, the amount of chemical fertilizer and number of livestock at the end of each year in the study areas were obtained from local economic statistical datasets, i.e., the Inner Mongolia Statistics Bureau, over the study period (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015). The price of corn was obtained from China Animal Industry Yearbooks (2001-2015).

Methodology
To discover different human-climate-vegetation change mechanisms in the study areas, two general steps were taken to profile human-induced vegetation change ( Figure 2). (1) We detected vegetation dynamics from 2000 to 2014 to find out areas of significant change during the past 15 years, and (2) we used optimized RESTREND to discriminate human-induced changes in areas with significant change. Areas with no significant vegetation changes determined in the first step were not analyzed in the second step.

Detecting Vegetation Dynamics Using the SMK Method
Statistical methods are used frequently to detect trend dynamics of vegetation [49,50], including parametric methods such as the least squares linear regression method [1,51] and non-parametric methods such as Sen's slope estimation coordinated with the Mann-Kendall test detection (SMK) method [52,53]. In this study, the SMK method was used to determine whether a time series trajectory of ∑ NDV I indicated a statistically significant trend (p < 0.05). It is a non-parametric method that is robust against non-normality in a dataset and has good performance in error resistance [54,55], which is suitable for use in arid and semi-arid ecosystems with strong and frequent interannual vegetation variability. The trend slope can be defined by Sen's slope estimation, and the significance of the slope can be determined using the Mann-Kendall test [52,53].
Besides, the ∑ NDV I change rate over the study period was determined using Equations (1) and (2), and can directly express the magnitude of trend change: where S ∆NDV I is the ∑ NDV I change rate (percentage of ∑ NDV I increase over the study period), X is the mean of the time series of ∑ NDV Is, t is the mean value of the dates, and t 1 and t n are the start and end years of the time series, respectively. The significant change areas can be acquired using Equation (3): where p is the change rate threshold, which was set as 10% in this study.

Discriminating Human-Induced Changes Using Optimized RESTREND
In the significant change areas, we used optimized RESTREND to profile human-induced changes. The first step for optimized RESTREND was pre-analyzing the relationship between climate factors and AA-NDVIs to determine suitable climate factors and their best related period in the study areas. Second, a multiple regression model was built at the pixel level to simulate ideal climate-vegetation relationships during the study period. Finally, the residual analysis was conducted. If the changes in NDVI response were only due to climatic effects, the residuals of the predicted NDVIs and actual NDVIs in the regression model should be close to zero, and only fluctuate randomly without significant trend signals. Any significant trend in the residuals will indicate changes that are not due to climatic effects [1]. A decreasing trend in residuals represented human-induced vegetation degradation, while an increasing trend in residuals indicated that vegetation improvements were happening in the study area because of human activity. If there was no trend in the residual series, all of the vegetation changes that happened during study period were driven by climate change [28].
(1) Find the most suitable climate variables and their best-related period.
where r is the correlation coefficient; x i is the value of AA-NDVIs in the ith year; x is the average value of x i from 2000 to 2014; y i is the value of the specific climate index (Table 1) in the ith year; and y is the average value of y i from 2000 to 2014. In each pixel, there were a total of 135 (3 × 45) r results for climate factors (P/T/W). The best-related period of T, P, and W (hereafter T best , P best , and W best , respectively) were the most suitable climate variables to build the climate-vegetation model.
Treating T best , P best , and W best as the representatives of the climate factors, we first normalized each variable using the min-max normalization method as in Equation (5), and then used multiple regression analysis to build the pixel-level climate-vegetation model as in Equation (6).
where v is a specific time series variable (i.e., T best , P best , and W best ); v i is the value of variable v at the ith time; Norv i is the normalized value of variable v at the ith time; and the max(v) and min(v) is the maximum and minimum values of time series variable v, respectively.
where y is the time series ∑ NDV Is of the 'reference pixel', which refers to the maximum value of ∑ NDV Is in the homogeneous spatial neighborhood of the target pixel; P, T, and W are the normalized values of the time series T best , P best , and W best respectively; and β 0 , β 1 , β 2 , and β 3 , are the parameters of the model. In this paper, the location of a corresponding reference pixel for a specific pixel x was determined by three steps (Figure 3).
Step 1 found all of the neighborhood pixels that were located within the same ecological zone of x, and marked them as candidates. In step 2, candidates that could not simultaneously satisfy the following two conditions lost their qualifications: (1) there must be a significant correlation of ∑ NDV I trajectories between x and the candidates; (2) the differential of mean value of the ∑ NDV I trajectory between x and candidates must be below 3δ. The 3δ method has been used to detect land cover change in some studies [14,56], and helped determine homogeneous spatial neighborhoods in this study. All of the candidates selected by this step constituted the homogeneous spatial neighborhood of x. In step 3, the candidate with the maximum value of ∑ NDV Is was defined as the reference pixel. (3) Residual analysis to discriminate human-induced changes Residual analysis was conducted using the SMK method though Equations (1)- (6). Residuals refer to the differential value between the potential ∑ NDV Is and the actual ∑ NDV Is. A significant decreasing or increasing trend in residuals indicated that vegetation change is driven by human activities.

Vegetation Dynamics during the Past 15 Years
Over the past 15 years, obvious vegetation increases happened in both regions, and Naiman increased more dramatically than Ar Horqin. In the Ar Horqin region, about 18.6% of the total area (2371.1 km 2 ) had a significant trend of change (p < 0.05), with 16.8% of the area increasing by more than 10%, and only 1.3% of the area decreasing more than −10%. About 66.3% of the Naiman region (5301.9 km 2 ) showed a significant trend change during this period, with much of the change area (64.1%) increasing by more than 10%, and only 1.1% of the area decreasing by more than −10%. About half of the change area increased between 20-40% during the study period ( Table 2). In both regions, most positive changes happened in the southern study area showing a zonal or block distribution. The negative change areas showed a sporadic distribution in the northern study area (Figures 4a and 5a). Specifically, positive changes always happened along rivers and roads and negative changes mainly happened in the south of the Greater Khingan Range Mountain in the Ar Horqin region. In the Naiman region, negative changes were mainly located around the county center.

Significant Human-Induced Change in the Ar Horqin Region
In the Ar Horqin region, both climate change and human activities contributed equally to the vegetation changes over the past 15 years, according to the trend of residuals shown in Figure 4b. About 50.1% of the significant changes resulted from human activities, while 49.9% of the significant changes were caused by climate change.
However, human activities caused more dramatic changes than climate in Ar Horqin. The results showed that vegetation in most of the human-induced change area (about 57.2%) increased by more than 40%, while climate change resulted in less rapid changes, with most of the climate-induced change area (about 56.7%) increasing between 20-40%.
The results also showed that human activities were more responsible for a vegetation decrease than climate change in the Ar Horqin region. About 85.3% of the decreases were caused by human activities. Vegetation increases over the past 15 years may be equally induced by human activities (47.3%) and climate changes (52.7%).

Significant Human-Induced Change in the Naiman Region
In the Naiman region, results demonstrated that human activities were the main drivers of vegetation change between 2000-2014, contributing up to 60.3% of the total significant changes ( Figure 5b). Furthermore, human activities were more responsible for both negative and positive vegetation change than climate. About 84.7% of the decreases were caused by human activities, and about 59.8% of the vegetation increases were caused by human activities. Climate change resulted in more vegetation increases (40.2%) than decreases (14.3%). The results indicated that human activities had deeper and broader influences in the Naiman region compared with the Ar Horqin region.

Validation of Human-Induced Vegetation Change
This study aimed to discriminate significant human-induced changes in two different regions. However, the validation of vegetation change at regional scales is notoriously hard because there is seldom ground-truthed data of historical land cover changes, as these data are difficult to collect over long periods [57]. Here, we used time series Landsat images with a high spatial resolution (30 m) to validate vegetation changes and try to find further evidence of human-induced changes through visual interpretation of the images [58]. Visual interpretation of false color composite images showed that our methods effectively discriminated significant human-induced vegetation change in both regions. Figures 6 and 7 show evidence that indicates land cover changes and driving forces during the study period in a typical change area.   An in situ investigation of vegetation increases in grassland showed that in Figures 6f-j and 7f-j, all of the increases were caused by human activity, e.g., grassland management. The conversion of natural grassland into artificial pastures was verified by the appearance of several circular grasslands, as found in Figures 6j and 7i,j. Two photographs of the circular patches are shown in Figure 8, demonstrating clear signs of the human management that resulted in a vegetation increase in these areas.

Different Driving Factors in the Study Areas
Although human factors contributed more to these changes than climate change in both regions, the specific human activities that caused the changes in the two regions were different.
In the Ar Horqin region, about 71.0% of the total area was grassland. Human activities caused 47.3% of the vegetation increase and 85.3% of the vegetation decrease. Human eco-managements (i.e., grass land restoration in Figure 8) and overgrazing were the two main human activities with positive and negative effects, respectively. Figure 9a further supported this point of view. As we mosaicked the land cover map in 2010 by human-induced vegetation change area, it was clear that 57.9% of the human-induced changes were happening on grassland. Figure 4a also showed that most of the human-induced grassland changes were positive changes, which can be attributed to human management such as the Beijing-Tianjin Sandstorm Source Controlling Program (introduced in 2001) and the Ecological Subsidy and Award System (introduced in 2011) that were carried out in this region. Figures 5f-j and 4b indicated the effects of human management on the grassland located in the southeast corner of the Ar Horqin region. However, most of the human-induced grassland decrease that happened in the northern area shown in Figure 4b was also caused by human activities according to Figure 9a, and these decreases can be attributed to overgrazing. We believe the effect of positive human management of grassland is greater than the negative outcomes of overgrazing, although the grazing activities increased in the Ar Horqin region during the study period (Figure 10b).
Farming is not a key human factor of change in this area, because only 33.9% of human-induced changes happened on farmland, as shown in Figure 9a. The majority of the significantly changed areas of farmland are located around the county center and spread in a north-south direction along the valley (Figure 4a). However, Figures 4b and 9a showed that only part of the significant farmland changes in this area was attributed to human activities. The amount of chemical fertilizer used only had a small increase (Figure 10a), although 13.6% of the total area is farmland.
In the Naiman region, farming was the key human driver of significant vegetation change. About 46.5% of total area is farmland, and 66.3% of the human-induced significant changes happened in farmland, as shown in Figure 9b. The majority of farmland was in the southern area, and some was located along the main road in a northeast-southwest direction. Human-induced farmland changes were almost all positive, as shown in Figures 5b and 9b. The Naiman region is a typical producing area of corn, and 90% of the farmland in this region is corn. The vegetation increase that occurred in farmlands in Naiman was mainly because of market incentives and the sharp increase in using fertilizer. According to the China Animal Industry Yearbook (2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015), the corn price in Inner Mongolia kept rising from 0.81 yuan/kg to 2.32 yuan/kg, with an increasing rate of up to 186.4% from 2000 to 2014. We believe that it greatly increased local people's willingness to plant corn and improve corn output. In addition, statistical data also showed that the amount of chemical fertilizer used in 2013 was up to 6.7 times greater than in 2000 (Figure 10a), while the area of farmland in 2011 was up to 1.5 times more than in 2000. This could be further evidence of the human activities that positively influenced farmland changes.   (Figure 10b). Meanwhile, the area of decreased grassland caused by overgrazing was still small. Therefore, we believe that the human management of grassland was the other important factor of vegetation change in the Naiman region, which means the Grazing Withdrawal Program (introduced in 2003) and the Ecological Subsidy and Award System (introduced in 2011) implemented in this region had positive eco-effects.
The increased number of livestock in the study area is also noteworthy (Figure 10b). Previous study had proved that the increased number of livestock at the end of each year was mainly because of market incentives [59]. It is reported that the market prices of cattle, sheep, and pork had grown rapidly in Inner Mongolia. The price of sheep raised from 14.77 yuan/kg to 65.41 yuan/kg. The price of pigs raised from 10.1 yuan/kg to 22.48 yuan/kg. Especially, the price of cattle raised from 12.88 yuan/kg to 63.29 yuan/kg, with an increasing rate up to 26.09% per year. These price increases were a great incentive for people in the Naiman region, because the Naiman region is a typical producing area of cattle. According to our field investigation, the increasing supply of forage, which was caused by the increasing number of livestock, was no longer from overgrazing. People used the subsidy given by Grazing Withdrawal Program from the government to buy forage or planted silage corn and alfalfa to give fodder to livestock.

Limitations and Further Research
This study profiled significant human-induced vegetation change in the study area. We first used the SMK method to find areas of significant change, and then used RESTREND to discriminate human-induced from climate-induced change. We used SMK to mask slightly changed and no change areas, i.e., change rates between −10% and 10%, because the variability of vegetation growth in arid/semi-arid areas is commonly high, and this extent of change was regarded as a relatively stable growth rate in the study area. However, this kind of area can also be influenced by human activities, and even further by the complicated interaction and offset mechanisms between climate change and human activities, such as increased irrigation in less rainy years, which need to be analyzed in further studies.
The strategy in optimized RESTREND had two advantages compared with previous studies [31][32][33]. First, the algorithm flow was automatic, and did not require human intervention, which could strengthen its applicability and ease-of-use in large areas. Second, we used a time series vegetation index to finish defining the flexible homogeneous areas, which avoided uncertainties caused by adding external time-varying data such as land use maps. Although we have detailed analysis of the human-induced vegetation change in the two regions, it is hard to validate the result directly and compare it with other methods. Here, we took a time series Landsat image dataset as the main source of validation data, and the attribution analysis mainly relied on visual interpretation. Thus, the results may contain some uncertainty. Accumulating a site investigation dataset over a long period and using a questionnaire survey about the effects of specific eco-programs are both very important in future studies.
The uncertainty of the results could be also caused by the coarse temporal resolution of meteorology datasets (i.e., one month), which was different form the time series NDVI datasets (i.e., 16 days). Using the highest correlation coefficient to choose the best-related period between NDVI and each climate factor might not be the best strategy. The partial correlation coefficient could be used in the future study to improve the process. The suddenly appearances of wildfire, diseases, and insects invasion [4,60,61] in ecosystems are neither direct human-related nor climate-related factors. Meanwhile, the drought and flood that happened during the study period may influence the linear regression relationship between climate and vegetation, which is the basic of RESTREND [5]. Results would be unreliable when these interferences happened. Breakpoints detection methods [56][57][58] can be used in combination with RESTREND to accurately define specific human-induced changes in the future study.
In addition, this study did not consider the possibility of changes in dominant driving factors over time as we used linear analysis methods. However, vegetation increases during past 15 years might reflect climate-induced negative changes in the first five years, and then human-induced positive changes in the following 10 years, and vice versa. A finer-scale temporal analysis would help define human-induced changes. Piecewise analysis methods [11] or other time series detection methods [62,63] can be used in further research to resolve this problem.

Conclusions
This paper first detected vegetation dynamics using SMK in the Ar Horqin and Naiman regions of China from 2000 to 2014, and then applied optimized RESTREND to profile significant human-induced changes in vegetation. We optimized strategies to define homogeneous spatial neighborhoods in RESTREND. An automatic algorithm was proposed to define flexible homogeneous areas based on the statistical analysis and trajectory analysis of a time series NDVI dataset.
Results showed that during past 15 years, significant vegetation increases happened in both regions, and climate change and human activities both contributed to the changes. However, more significant vegetation increases happened in the Naiman region (64.1%) than in the Ar Horqin region (16.8%). In the Ar Horqin region, climate and human activities equally contributed to vegetation change. Meanwhile, in the Naimen region, it was proved that human activities contributed more than climate factors to vegetation change. Human factors had a stronger influence on ecosystems and were responsible for the vegetation decrease in both regions. In different study areas, dominant human driving forces varied. In the Ar Horqin region, human management of grassland was the key factor, and the secondary human factor was overgrazing. Farming was not a key human factor for change in this area. In the Naiman region, farming was the key human factor for vegetation change. Positive human management in grassland was the other key driver of change.
In future studies, the improvement of RESTREND and finer-scale temporal analysis would help to accurately define specific human-induced changes. The study of a very long series of palaeoenvironmental records can also improve knowledge on the drivers of changes over long time periods in these regions [64]. The complex interactions and offset mechanisms between climate change and human activities need to be clarified so that the underestimates or overestimates of human effects can be avoided. In addition, we did not consider the trend-reversal phenomenon, which is very likely to happen over time, when analyzing the dominant driving factors. Some abrupt change detection methods can be used in further research.