An Integrated Approach for Detection and Prediction of Greening Situation in a Typical Desert Area in China and Its Human and Climatic Factors Analysis

The combined study of vegetation coverage (VC) and land use change provides important scientific guidance for the restoration and protection of arid regions. Taking Hongjian Nur (HJN) Lake in the desert region as a case study, the VC of this area was calculated using a normalized difference vegetation index (NDVI), which is based on a mixed pixel decomposition method. A grey forecasting model (GM) (1, 1) was used to predict future VC. The driving factors of VC and land use change were analyzed. The results indicate that the average VC of the whole watershed showed a gradual increase from 0.29 to 0.49 during 2000–2017. The prediction results of the GM VC showed that the greening trend is projected to continue until 2027. The area of farmland in the watershed increased significantly and its area was mainly converted from unused land, grassland, and forest. The reason for increased VC may be that the combination of the exploitation of unused land and climate change, which is contrary to the country’s sustainable development goals (SDG; goal 15). Therefore, the particularities of the local ecological environment in China’s desert area needs to be considered in the development of ecological engineering projects.


Introduction
Since the 20th century, the global average temperature has significantly increased. Global warming has brought about serious threats to human society, the economy, and the ecosystem [1]. Natural lakes in desert areas often provide important ecosystem functions such as conservation of biological diversity and carbon storage [2]. The ecosystems surrounding lakes in arid regions are important because they produce numerous beneficial services [3,4]. Vegetation is the main support system (terrestrial ecosystem) on which human beings depend for survival and sustainable development. The change of vegetation cover can reflect the situation of a regional ecological environment [5]. The changes in vegetation coverage and land cover type surrounding lakes in arid regions affect whole wetland ecosystems, which may be serious impacted by surrounding human life, biodiversity, and ecological environment [6][7][8]. Seasonality in arid Africa requires the assessment and prediction of vegetation response due to climate change. Scholars have used remote sensing to indirectly measure vegetation growth by calculating the vegetation index [9]. Research in the Loess Plateau demonstrated that, in addition to improving vegetation coverage, vegetation restoration projects may also bring negative impacts on local water resources [10]. In desert or arid regions, the greening situation followed by The studies above have partially revealed the characteristics and trends of VC and land use change in the HJN region. The possible factors driving vegetation change in the HJN region were also discussed by previous studies [11]. However, the changes of the overall ecological environment of HJN is inadequate. This may be because the relevant research data sources are varied and the time series is short. Most of the existing studies separately analyzed VC and land use in the HJN region. Remote sensing data sources for vegetation cover and land use change studies include hyperspectral data, multispectral data, microwave data, and LiDAR data [29]. Moderate-resolution imaging spectroradiometer (MODIS) products have been widely used in regional VC and land use change investigations due to their wide coverage and high temporal resolution. Many models have been developed to detect VC and land use based on MODIS data [30]. The commonly used models include regression models, mixed pixel decomposition methods, and machine learning methods [18,31,32]. The application of the regression method as a regional empirical model to estimate VC at a large scale may cause problems because it is unable to accurately describe the complex land surface conditions [33][34][35]. However, machine learning methods are somewhat limited in practical applications due to their computational complexity [36]. The mixed pixel decomposition model is a commonly used method for calculating regional VC based on hyperspectral data. The binary pixel model is the most common model used in the linear mixed pixel decomposition method [37]. It is widely used to estimate VC due to its simplicity of application and the way it represents physical mechanisms [29]. Zribi et al. used the binary pixel model to decompose the radar ERS2/SAR (European Sensing Satellite-2/Synthetic Aperture Radar) signal and retrieved the vegetation coverage in the semiarid area [38]. Qi j et al. combined NDVI with a binary pixel model and analyzed the spatiotemporal dynamic change characteristics of vegetation coverage in the San Pedro basin in the southwestern United States [39]. Many other authors have used the binary pixel model to retrieve vegetation coverage and to explore its change trends over the years as well as the driving forces that affect these changes [40]. In addition, we aimed to describe the rapid reduction in lake area in the HJN region in recent years and identify any relevant relationships. We used GIS (geographic information science) technology to analyze the land use map overlay to obtain the land use transfer matrix and studied the mutual conversion relationships of land use. This method has been applied to the study of land use change in various regions [41,42].
In summary, it is still unclear why there is a long-time trend of comprehensive vegetation environment around the desert wetland. Thus, the main objectives of this study are as follows.
(1) Herein, we try to reveal the dynamic change characteristics of the greening situation followed by vegetation restoration projects and their human and climatic factors based on both VC and land use monitoring. (2) We assume that the influence factors remain unchanged, simulating the long-time trend of the VC that will help provide reference for future desert wetland ecological environment protection.
(3) We investigate the influencing factors for the greening situation around the typical desert wetland from the perspective of human and climate change. This study will help provide a scientific basis for the restoration of the ecological and economic functions of desert wetland area under vegetation restoration projects and climate change backgrounds.

Study Area
HJN Lake is located at 39 • 04 ~39 • 08 N and 109 • 50 ~109 • 56 E (Figure 1). In recent years, due to the dryness of the Taolimiao Alashan Nur in Ordos, Inner Mongolia, China, relict gulls have gradually chosen to move to and roost in HJN Lake. HJN Lake has become a crucial breeding and habitat for relict gulls. The HJN nature reserve is in northern China, at the southern edge of the Mu Us Sandy Land and between the Yike Zhao League of Inner Mongolia and the North Shaanxi Loess Plateau. It has the characteristics of a plateau erosive hilly landform. The HJN watershed is mainly controlled by the cold continental polar air mass and the influence of the marine hot air mass is low. The HJN watershed has continental monsoon climate characteristics because it is in the desert and in a loess boundary zone on the northwest edge of the East Asian monsoon region. This region has four typical seasons. Wind evaporation is strong in this region, which has low precipitation and rapid warming in spring. The summer is short with large temperature differences and more thunderstorms than other seasons. It is relatively cool in autumn but cold in winter, with less snow and rain than in other seasons. The HJN region has a grassland-desert transition landscape as a consequence of the low precipitation and high evaporation with sparse vegetation. Because of the specific geographical environment and climatic conditions of the HJN region, the vegetation and land use conditions are unique. the dryness of the Taolimiao Alashan Nur in Ordos, Inner Mongolia, China, relict gulls have gradually chosen to move to and roost in HJN Lake. HJN Lake has become a crucial breeding and habitat for relict gulls. The HJN nature reserve is in northern China, at the southern edge of the Mu Us Sandy Land and between the Yike Zhao League of Inner Mongolia and the North Shaanxi Loess Plateau. It has the characteristics of a plateau erosive hilly landform. The HJN watershed is mainly controlled by the cold continental polar air mass and the influence of the marine hot air mass is low. The HJN watershed has continental monsoon climate characteristics because it is in the desert and in a loess boundary zone on the northwest edge of the East Asian monsoon region. This region has four typical seasons. Wind evaporation is strong in this region, which has low precipitation and rapid warming in spring. The summer is short with large temperature differences and more thunderstorms than other seasons. It is relatively cool in autumn but cold in winter, with less snow and rain than in other seasons. The HJN region has a grassland-desert transition landscape as a consequence of the low precipitation and high evaporation with sparse vegetation. Because of the specific geographical environment and climatic conditions of the HJN region, the vegetation and land use conditions are unique. Figure 1. Study area. In the figure, the elevation increases gradually from blue to brown. Blue is low and brown is high.

MODIS Data
The MODIS data used in this study are derived from the Land Processes Distributed Active Archive Center (LP DAAC) managed by the NASA Earth Science Data and Information System (ESDIS) project [43]. The MOD13A2 vegetation index with 16-day 1 km SIN Grid V006 (standard Figure 1. Study area. In the figure, the elevation increases gradually from blue to brown. Blue is low and brown is high.

MODIS Data
The MODIS data used in this study are derived from the Land Processes Distributed Active Archive Center (LP DAAC) managed by the NASA Earth Science Data and Information System (ESDIS) project [43]. The MOD13A2 vegetation index with 16-day 1 km SIN Grid V006 (standard resolution/globally produced) was used in this study to calculate the VC in the study area. There are two vegetation index (VI) values in these products: NDVI and the enhanced vegetation index (EVI). The quality assurance (QA) layers associated with NDVI and EVI were also provided with the products. As the HJN watershed has a relatively low VC, only the NDVI layer was used in our study. All cloud-contaminated pixels were removed based on their corresponding QA value.

Landsat Data
Landsat satellite data with a multispectral resolution of 30 m and a revisit period of 16 days were mainly used to detect land use/coverage change in the HJN watershed. In this investigation, Landsat TM/ETM+ images mainly from July to August during 2006-2016 were used to detect the land use types and their changes in the HJN region. The data used were primarily downloaded from the GloVis database of the USGS (United States geological survey) [44]. The downloaded data are the L1T product (Level 1T). According to the evaluation results of geometric accuracy stability, the L1T product has high geometric accuracy [45]. We used the Landsat Ecosystem Distribution Adaptive Processing System (LEDAPS) software to convert the DN value of each pixel in the image into the surface reflectance. The atmosphere correction mainly used the FLAASH Atmospheric correction of atmospheric correction module in envi5.3. We used a newly developed algorithm, F-mask, to remove clouds and shadows in the image [46]. All pixels affected by clouds and shadows were removed and only clear-sky observations were used for our study.

Meteorological Data
There are two meteorological stations Ejin Horo and Shenmu located around the study area. The daily climate data from the two stations was derived from National Meteorological Information Center [47]. Both stations had meteorological data from 1961 to 2016. However, data from many days were missed in November 1968. The missing data were replaced by an average value from two adjacent years to maintain the integrity of the data sequence.
All input data related information is shown in Table 1.

VC Calculation
The mixed pixel decomposition model is a method for calculating VC [48]. It assumes that each pixel of a remote sensing image is composed of vegetation and soil. The value of each pixel contains information about its soil and vegetation composition, while the NDVI value of the mixed pixel is the weighted average sum of the VI (vegetation index) values based on the soil and vegetation composition of the pixel. The expression for calculating NDVI is as follows [49]: where is the value of the mixed pixel, NDVI veg is the VI value of the pure vegetation pixel, NDVI soil is the vegetation index value of a pure soil pixel, and f v is the VC. Therefore, the VC can be obtained as follows [49]: where NDVI soil is approximately equal to the minimum value of a pure soil pixel (close to 0, in theory) and NDVI veg is equal to the maximum value of a pure vegetation pixels (close to 1, in theory). NDVI soil and NDVI veg will change in different times and places due to the influence of meteorological conditions, vegetation types, and their distributions. Seasonal changes and other factors may also influence NDVI soil and NDVI veg . The preprocessed NDVI data were treated as mixed pixel vegetation index values. First, the yearly NDVI values were obtained from NDVI images of the vegetation-growing season from May to October in 2000-2017 through the maximum synthesis method. The overall NDVI values in the HJN region were low. Then, to avoid overestimating the VC result, the maximum and minimum NDVI values over the whole tile covering the study area were calculated as NDVI veg and NDVI soil , respectively. Finally, the yearly VC of the HJN watershed and the surrounding counties was calculated pixel by pixel with formula (2).

Land Use Analysis
Land use and cover change (LUCC) data were mainly obtained through the interpretation of remote sensing images per biennium during 2006-2016. All Landsat data were processed by radiometric correction and atmospheric correction. Interpretation keys were established by artificial methods. The land cover system of the International Geosphere-Biosphere Program (IGBP) was used in our study. Based on the actual land use/cover situation in the HJN region and the research results from landscape classification, the land use types were divided into woodland, grassland, water body, urban land, rural land, other construction land, unutilized land, and farmland. The satellite data were interpreted by a combination of manual interpretation and supervised classification methods. Too-small or surplus spots were eliminated or merged. The interpretation accuracy of Kappa index was calculated based on field observation data and interpretation results. Kappa coefficient is a multivariate statistical method to evaluate the accuracy of classification, which represents the proportion of error reduction in the evaluated classification compared with the completely random classification. The calculation formula is as follows:ˆK whereˆK is kappa coefficient, r is the number of rows of error matrix, x ii is the value on the value (main diagonal) of row i, x i+ , and x +i are the sum of row i and column i, respectively. n is the total number of samples. The land use interpretation results were consistent with the requirements for the analysis of land use dynamics in this study.

VC Prediction
VC was forecasted using data from the HJN watershed, Ejin Horo and Shenmu City from 2000 to 2017. VC in 2000-2027 was predicted by a grey model (GM) and a chart of the relationships between real values and predicted values were obtained. The grey model is used to predict the grey process that changes in a certain range and is related to time. Grey models are based on grey systems, which have the fuzziness of hierarchical and structural relations, the randomness of dynamic change, and the incompleteness or uncertainty of index data. This model does not require a large amount of regular sample data is easy to use and has high prediction accuracy.
GM (1, 1) model represents the first order differential equation model with one variable. The prediction principle of the GM (1, 1) model is the generation of a set of new data sequences with obvious trends accumulated from certain data sequences. The growth trend of this new data sequence is used to build the prediction model. Then, the accumulative method is used for reverse calculation to restore the original data sequence. Finally, the VC series can be forecasted.
The modeling principle of the grey model is as follows: (1) Set the known reference data as: x (0) , then sum once to get, Among them, (2) Take the average sequence Then, (3) Modelling The albinism equation of sequence x (1) is, The above formula is discretized to obtain the grey differential equation of GM (1,1), (4) Solve for the equation parameters a and b Using the least square method, is obtained, Achieve the minimum value,û (5) Establishment of prediction formula The prediction formula of x (1) , The prediction sequencex (1) of sequence x (1) is obtained and the prediction sequencex (0) of x (0) can be obtained by reducingx (1) , For the grey prediction model, the mean square deviation ratio C and the small error probability p are generally needed. The smaller the C value is the closer the p value is to 1 and the better the prediction effect of the model is. Otherwise, the worse the prediction effect is.

Zonal Statistics Method
Zonal statistics were used to analyze the general VC and land use conditions. With the zonal statistics method, a statistic is calculated for each zone defined by a zone dataset, based on values from another dataset (a value raster). A single output value is computed for every zone in the input zone dataset. Due to the regional difference of NDVI value, the mean values of VC were calculated for the study area and the surrounding districts and counties. These values are used to establish a discount chart for further analysis. The total area of each land use type in the HJN watershed was calculated to analyze the temporal and spatial land use change by the thematic map. In order to prove the accuracy of the average value of VC for more in-depth analysis, we analyzed the normality of the values in different regions for many years. A conventional descriptive statistical method was used to calculate the meteorological data for climate change analysis and to study the drivers of VC and land use variation.

The State Transition Matrix (STM)
STM has been widely used to analyze changes in land use types. This method was developed by Markov and was originally used to describe the probability of state transition in mathematical analysis [50]. STM not only relates the state of each land use type from the initial time to the current time but also propagates the covariance of the land use state. The land use change transition matrix can be defined as follows: where → L is the transfer matrix of land use change, with superscript T transposing each row vector to a column vector in a given time, and → x 1 and → x 2 are the different land use types. In order to calculate the transfer matrix, we need to dissolve the two images rely on GIS technology. Then analyze the overlay of the two images and calculate the area of the intersection area. Finally, the transfer matrix is calculated by using the statistical tools.

Linear Regression
Linear regression was used to analyze the variation trends in land coverage and climate data. Linear regression is a robust model for obtaining the relationship between independent variables and dependent variables based on a least squares function. In this study, the precision of the regression or fitting was determined by the correlation coefficient and a T-test. The correlation coefficient we used was the Spearman's rank correlation coefficient. In the correlation analysis, we used the Spearman's rank correlation coefficient. Before calculating the coefficients, we first used SPSS 25 and Origin 2018 to verify the normality of the data, and then chose a more appropriate correlation coefficient, the Spearman's rank correlation coefficient does not need to assume that the relationship between variables is linear or measure variables on an interval scale. It can be used for variables measured at the ordinal level. It has no strict requirements on the data and has a wider range of application [51].

Results
In our study, we calculated the change in vegetation coverage from 2000-2017 in HJN Lake and its surrounding areas and the change in land use from 2006 to 2016 and predicted the vegetation coverage from 2017-2027. The results showed that the vegetation coverage in HJN Lake and the surrounding areas increased significantly, and the water area of HJN Lake showed an obvious decrease.

VC Change in the Region Surrounding HJN Lake
To understand the overall vegetation condition, the VCs of six surrounding counties were calculated and extracted. The overall VC (0.2-0.6) of the HJN watershed is lower than that of the whole calculation region (Figure 2). its surrounding areas and the change in land use from 2006 to 2016 and predicted the vegetation coverage from 2017-2027. The results showed that the vegetation coverage in HJN Lake and the surrounding areas increased significantly, and the water area of HJN Lake showed an obvious decrease.

VC Change in the Region Surrounding HJN Lake
To understand the overall vegetation condition, the VCs of six surrounding counties were calculated and extracted. The overall VC (0.2-0.6) of the HJN watershed is lower than that of the whole calculation region (Figure 2). The VC in the six counties around HJN Lake shows obvious spatial heterogeneity. The overall VC of Ejin Horo Banner, Uxin Banner, and Dongsheng District in Ordos city is lower than that of Shenmu County, Fugu County, and Yuyang District in Yulin city. Although they are located in the same city, Shenmu, Fugu, and Yuyang have different VC. Generally, the VC in the east is higher than that in the west in those counties. Shenmu and Fugu have higher VC than Yuyang. The eastern and northern counties as well as most parts of Fugu County have significantly higher VC than the other counties. In Yuyang County, the VC is higher only in the southeastern hilly area. The VC of Ejin Horo Banner, Dongsheng District and Wushen Banner in Ordos City is also uneven. The VC of Ejin Horo The VC in the six counties around HJN Lake shows obvious spatial heterogeneity. The overall VC of Ejin Horo Banner, Uxin Banner, and Dongsheng District in Ordos city is lower than that of Shenmu County, Fugu County, and Yuyang District in Yulin city. Although they are located in the same city, Shenmu, Fugu, and Yuyang have different VC. Generally, the VC in the east is higher than that in the west in those counties. Shenmu and Fugu have higher VC than Yuyang. The eastern and northern counties as well as most parts of Fugu County have significantly higher VC than the other counties. In Yuyang County, the VC is higher only in the southeastern hilly area. The VC of Ejin Horo Banner, Dongsheng District and Wushen Banner in Ordos City is also uneven. The VC of Ejin Horo Banner is the highest, followed by that of Dongsheng district, and Wushen Banner has the lowest VC of the three.
According to the statistical analysis of the average vegetation coverage of Shenmu and Ejin Horo, the two groups of data conform to a normal distribution; the distributions are shown in the Figure 3. Banner is the highest, followed by that of Dongsheng district, and Wushen Banner has the lowest VC of the three.
According to the statistical analysis of the average vegetation coverage of Shenmu and Ejin Horo, the two groups of data conform to a normal distribution; the distributions are shown in the Figure 3. From 2000 to 2017, the VC in the six county-level administrative divisions significantly changed. The overall VC in the whole region increased and the growth trend was obvious. The VC significantly increased in Shenmu and Ejin Horo Banner, which are surrounded by the HJN watershed (Figure 4). The VC in Dongsheng, Wushen and Yuyang also showed an increasing trend, but the change was not as obvious as that in the first two areas. From 2000 to 2017, the VC in the six county-level administrative divisions significantly changed. The overall VC in the whole region increased and the growth trend was obvious. The VC significantly increased in Shenmu and Ejin Horo Banner, which are surrounded by the HJN watershed (Figure 4). The VC in Dongsheng, Wushen and Yuyang also showed an increasing trend, but the change was not as obvious as that in the first two areas.

VC Change in HJN Watershed
The VC values in the HJN watershed from 2000 to 2017 were extracted as shown below ( Figure  5).

VC Prediction in the HJN Watershed and the Surrounding Regions
According to the calculated VC of HJN Lake and its surrounding areas during 2000-2017, a grey model was used to predict the VC from 2000 to 2027. The grey model can be used for short-, medium-, and long-term predictions. It also has good prediction accuracy for increase data. Figure 7 clearly shows that the predicted VC is consistent with the known VC in these areas. Therefore, the predicted value is highly credible and basically conforms to the development trend for vegetation coverage in the area. ISPRS Int. J. Geo-Inf. 2020, 9, x FOR PEER REVIEW 17 of 25 According to the calculated VC of HJN Lake and its surrounding areas during 2000-2017, a grey model was used to predict the VC from 2000 to 2027. The grey model can be used for short-, medium-, and long-term predictions. It also has good prediction accuracy for increase data. Figure 7 clearly shows that the predicted VC is consistent with the known VC in these areas. Therefore, the predicted value is highly credible and basically conforms to the development trend for vegetation coverage in the area. To verify the accuracy of the predictions based on the grey model, the model was tested. A smallerror probability P test and a variance ratio C test were performed on the model. As shown in Table  2, according to the existing rating scale that corresponds to the inspection results, the comprehensive inspection results are rated as good. Table 2. Prediction accuracy and level and ratings of the grey model. To verify the accuracy of the predictions based on the grey model, the model was tested. A small-error probability P test and a variance ratio C test were performed on the model. As shown in Table 2, according to the existing rating scale that corresponds to the inspection results, the comprehensive inspection results are rated as good.

Land Use in the HJN Watershed
In 2006-2016, the area of farmland in the HJN watershed increased the most of all land use types, with the area increasing significantly from 227.31 km 2 to 441.34 km 2 in 2010 to 2012, respectively, and remaining basically unchanged in the following years (Table 3). At the same time, some of the area around HJN Lake was used as urban land and other construction land. The areas of woodland, grassland, water body, and unutilized land decreased gradually in this decade. The water body area decreased obviously from 86.82 km 2 in 2006 to 45.35 km 2 in 2016 and the area of unutilized land decreased by up to 111.26 km 2 . Rural land changed the least during this decade. As shown in the Figure 8, there has been a remarkable change in the land use in the HJN watershed during 2010-2012. The area of farmland increased considerably, while the area of unutilized land decreased simultaneously. Part of the unutilized land was converted into grassland, and part of the grassland was converted into farmland. The water body also changed significantly from 2006 to 2016.  Table 4 shows that the total area of cropland and other construction land increased, while the total area of grassland, woodland, and water body decreased. The other land use types showed little change from 2006 to 2010. Among them, the total area of farmland increased by 25.52 km 2 , with a growth rate of 12.65%, which was mainly transformed from parts of the grassland and water body  Table 4 shows that the total area of cropland and other construction land increased, while the total area of grassland, woodland, and water body decreased. The other land use types showed little change from 2006 to 2010. Among them, the total area of farmland increased by 25.52 km 2 , with a growth rate of 12.65%, which was mainly transformed from parts of the grassland and water body areas. The total area of other construction land increased to 14.04 km 2 . The total area of grassland decreased from 788.18 km 2 to 771.37 km 2 . In addition, the total area of woodland and water bodies decreased by 1.61 km 2 and 16.58 km 2 , respectively, and the reduction rate of water bodies was 19.10%. It is also shown clearly in Table 4 that the decrease in the total grassland and water body area is due to their transformation into farmland. As shown in Table 5, the increase in total farmland and the decrease in unutilized land are the primary changes in land use types. The total area of farmland increased by 193.19 km 2 and the increase rate was 84.99%. The added farmland was mainly transformed from grassland (162.85 km 2 ). The total area of unutilized land decreased by 116.31 km 2 , with a decrease rate of 49.85%. The change in unutilized land and farmland was related to the transformation between these land use types and the grassland. The total grassland, woodland and water body area decreased by 38.7 km 2 , 13.38 km 2 , and 24.89 km 2 , respectively. The decrease rate of the water body area was 35.44% and it was mainly transformed into grassland, farmland, and woodland.

Impact of Climate Change on VC and Land Use
Related studies have shown that the HJN watershed is a typical inland watershed [23]. VC is significantly correlated with precipitation and temperature [52]. VC is greatly affected by precipitation. However, VC has different responses to temperature and precipitation in different seasons [11]. In spring and winter, the VC is mainly affected by temperature and elevated temperature can promote vegetation growth. In contrast, the VC decreases with increasing precipitation in winter. The vegetation is affected mainly by precipitation in summer and autumn. An increase in temperature can also promote the growth of vegetation [53].
As shown in Figure 9a, the temperature showed an obvious increasing trend from 1961 to 2016 at both stations. The cumulative temperature anomaly showed that the temperature in the study area has been rising continuously since 1990. The precipitation trend is not obvious compared with that of temperature. Precipitation fluctuated before 2010, but the cumulative precipitation anomaly showed a slow downward trend. After 2010, although the HJN watershed also experienced a brief drought, the overall cumulative precipitation anomaly increased significantly. The increased precipitation and temperature promoted the increase in VC in the study area. VC increased significantly with increasing precipitation and temperature.
with 218.71 km 2 growth. The transformation of woodland and grassland to farmland was the main reason for the change in the total area of woodland and grassland.

Impact of Climate Change on VC and Land Use
Related studies have shown that the HJN watershed is a typical inland watershed [23]. VC is significantly correlated with precipitation and temperature [52]. VC is greatly affected by precipitation. However, VC has different responses to temperature and precipitation in different seasons [11]. In spring and winter, the VC is mainly affected by temperature and elevated temperature can promote vegetation growth. In contrast, the VC decreases with increasing precipitation in winter. The vegetation is affected mainly by precipitation in summer and autumn. An increase in temperature can also promote the growth of vegetation [53].
As shown in Figure 9 a (1), the temperature showed an obvious increasing trend from 1961 to 2016 at both stations. The cumulative temperature anomaly showed that the temperature in the study area has been rising continuously since 1990. The precipitation trend is not obvious compared with that of temperature. Precipitation fluctuated before 2010, but the cumulative precipitation anomaly showed a slow downward trend. After 2010, although the HJN watershed also experienced a brief drought, the overall cumulative precipitation anomaly increased significantly. The increased precipitation and temperature promoted the increase in VC in the study area. VC increased significantly with increasing precipitation and temperature. Before calculating the correlation between temperature and humidity and VC value, since the sample size is less than 1000, Shapiro Wilk is selected to test the normality first. Through the data Before calculating the correlation between temperature and humidity and VC value, since the sample size is less than 1000, Shapiro Wilk is selected to test the normality first. Through the data test, we found that part of the data P value is low, the degree of conforming to normal distribution is low. Therefore, we used spearman coefficient to discuss its correlation. The normality test results of the data are as follows (Table 6). If P is greater than 0.05, the data is considered to be normal distribution.
The correlation coefficients between VC and the cumulative temperature anomaly and cumulative precipitation anomaly for the Shenmu and Ejin Horo stations are provided below.
As shown in Table 7, VC has a strong correlation with the cumulative temperature anomaly at both stations. The correlation coefficients between VC and the cumulative temperature anomaly reached 0.92 and 0.91 (significant at the 0.01 level) at the Ejin Horo site and Shenmu site, respectively. This indicates that the increase temperature plays an important role in the increase in vegetation in the study area. The precipitation increases also promote the increase in vegetation, with correlation coefficients of 0.73 (significant at the 0.05 level) and 0.43 between VC and the cumulative precipitation anomalies at both stations, respectively. As the increasing trend for precipitation is not obvious compared to the trend for temperature, the correlation coefficients between VC and the cumulative precipitation anomaly are lower than those between VC and temperature.

Impact of Climate Change on VC and Land Use
The data from the Shenmu and Ejin Horo meteorological stations near the HJN watershed from 1961 to 2016 show that the temperature of the two stations had an obvious upward trend during this period. Precipitation fluctuated greatly before 2010. After 2010, the total precipitation accumulation increased significantly. The correlation coefficient of VC with temperature and precipitation at the Shenmu and Ejin Horo stations shows that the increase in temperature plays an important role in increasing vegetation coverage in the study area. The increase in precipitation also promotes the spread of vegetation. This result is consistent with the existing research [11]. Xiu et al. (2019) found that there was a positive correlation between average annual temperature and NDVI during 1982-2015 [11]. The investigation in other regions across the world also demonstrated the obvious relationship between climatic factors. The study over Heihe River Baisin in Northwest China found that the relationships between the changes of NDVI and temperature were significantly positive. In addition, the relationships between NDVI and temperature are higher than the relationships between NDVI and precipitation [54]. Vali et al. (2020) found that the temperature and precipitation are the main climatic factors along with topographic factors which significantly affect the vegetation distribution and change trend [55]. The spatiotemporal vegetation cover variations were also found associated with climate change in the Loess Plateau [56,57]. All these regions are located in arid areas and have many similarities with HJN in climatic or geographical divisions. Therefore, the consistent research results in various area over the world proved the credibility of the conclusion of association between climatic factors and VC.

Effect of Afforestation on VC in HJN Watershed
In the '12th Five-Year Plan of China (2011-2015)', it was reported that the area of afforestation was the largest in the history of the HJN watershed and that afforestation projects had achieved remarkable results. The area preserved for various types of forest trees reached 3066 km 2 , with an increase of 266 km 2 compared to that at the end of the '11th Five-Year Plan of China (2006-2010)'. A total of 863 km 2 of afforestation and 441 km 2 of greening mileage were completed [58]. Based on GIMMS AVHRR normalized difference vegetation index (NDVI) datasets, Wang et al. (2011) found that the vegetation coverage area of the four construction areas increased significantly after the start of the Shelter Forest System Program in the three north regions of China [59]. The greening projects mainly include greening in and around urban areas and industrial parks, the construction of green promenades and camphor pine stands, desertification control, and greening and beautification in key villages and towns. In the monitoring of vegetation cover change based on NDVI in Patuakhali District, it is also found that afforestation has an impact on vegetation cover area [60]. In our study, we also found that after afforestation, the increase in VC in the HJN watershed occurred in the context of increased greening in China from 2000 to 2017 [17,34]. This is consistent with a previous study showing that forest and farmland are the most important contributors to greening in China [17].

Relationship among VC, Land Use, and Lake Area
Human activities have a high impact on land use and vegetation cover dynamics and the soil erosion around the watershed will directly lead to the shrinkage of lake area [61]. Different research areas have shown that there is a very close relationship between vegetation coverage, land use, and lake area, which is reflected in the studies of Tso MorIri Lake, the Hulun Lake (HLL) in China, and Poyang Lake in China [62][63][64].
From 2000 to 2017, the VC of the HJN watershed was generally on the rise. In 2017, the VC of the HJN watershed reached its highest value. In addition, the vegetation coverage in the HJN watershed and its surrounding areas in 2017-2027 predicted by the grey prediction model shows that it will continue to increase. According to the land use data from 2006-2016, the increase in VC is due to the increasing areas of grassland, woodland, and cropland. In 2010, the lake area decreased significantly, which is the same finding as in a study of water areas from 1957 to 2015 by Wangying et al. However, VC continued to increase. Figure 8 and the two transfer matrices show that the reduction in the water body area of Shenmu County and Ejin Horo Banner has a close relationship with the increase in farmland. In the vicinity of the HJN watershed in 2010, the lake area clearly decreased, and the farmland area increased. This trend is consistent with that found in Liuying et al.'s research on the change in land use types in the HJN region [20]. After 2010, the VC continued to increase, while the area of the lake in the HJN watershed continued to decrease. On the whole, the increase in cultivated land area leads to an increase in VC, but it may accelerate the shrinking of lake areas [11]. As shown in Figure 8, the water area increased obviously, while the lake area decreased year by year. This means that upstream water storage aggravates the reduction of the HJN Lake area.

Conclusions
Our study has shown the overall VC has increased in both the HJN watershed and its surrounding region, and the upward trend in VC was obvious. This trend in VC is also revealed in the results predicted by the grey model. The results show that VC will continue to increase in the next decade. According to the statistical results for land use in the HJN watershed in 2006-2016, the areas of farmland and construction land increased, while the areas of unused land and water bodies significantly decreased. The total area of unutilized land decreased by 111.25 km 2 and the reduction area was mainly transformed to farmland, woodland, or grassland. This indicates that afforestation and reclamation of farmland contribute more than other factors to the increase in VC in the HJN watershed. The increased precipitation and temperature are also beneficial to the increase in VC in the study area.
From 2000 to 2017, the VC of the HJN watershed showed an obvious upward trend, but the area of the lake decreased year by year. The increase in farmland near the HJN watershed may be the reason for the decrease in the lake size. The increase in farmland leads to an increase in VC, but it may also accelerate the decrease in lake area. The storage of water in the upstream dam further aggravates the reduction in lake size. Although human activity is an important contributor to greening in the HJN watershed, it has also led to the shrinking of HJN Lake. The ecological balance is broken in this watershed, which is located in an arid desert area. Therefore, the particularities of the local ecological environment in this desert area of China need to be considered in ecological engineering projects. For the HJN watershed, it is also important to maintain the lake area to protect relict gull habitat.