Spatio-Temporal Grassland Development in Inner Mongolia after Implementation of the First Comprehensive Nation-Wide Grassland Conservation Program

: Protection of the grassland’s ecological environment and improvement of people’s livelihoods are major tasks for the management of pastoral areas in Inner Mongolia. The comprehensive program for grassland conservation in China, the Subsidy and Incentive System for Grassland Conservation (SISGC), was launched in 2011. To comprehend the effects of this major step towards sustainable grassland development, this study focuses on the spatio-temporal development of grasslands in Inner Mongolia since 2011. Through the combination of MODIS (Moderate-resolution Imaging Spectroradiometer) satellite data with up to date meteorological data, we used the indicators of Fractional Vegetation Cover (FVC) and Net Primary Productivity (NPP) to analyze qualitative and quantitative grassland changes. A classification system on the pixel level, reflecting change trends and fluctuations of both FVC and NPP, was applied to monitor and analyze the grassland development from 2011 to 2019. In particular, the spatial transfer matrix of the recent two years (2018 to 2019) was analyzed to reveal the latest potential issues and random impact factors. The results show a positive overall but spatially unbalanced effect of SISGC, with a prominent positive impact in the semi-desert grassland area. The potential threats from both social and natural aspects as well as the importance of a forecast system for local stakeholders in the pastoral area are discussed.


Introduction
The development of the grassland ecosystem, one of the most vulnerable ecosystems, has aroused a critical worldwide concern in recent decades [1][2][3], especially under the gradually consistent perception of global climate change. Grassland in the Inner Mongolia Autonomous Region (IMAR) of China belongs to the Eurasian grassland habitat [4] and shows the second highest diversity of indigenous plants in the world [5]. Towards the beginning of the 21st century, IMAR suffered great grassland degradation [6,7]. Protecting the grassland's ecological environment and improving people's livelihoods are major governmental tasks. Recently, it also became an essential requirement for the Beautiful China Construction (BCC) initiative, which was initiatively planned by the 18th Chinese National Congress in 2015. BCC is the blueprint for ecological civilization (a construction goal for Chinese government to emphasize harmony between humans and nature [8,9]), of fundamental importance for the sustainable development of the Chinese nation and a Chinese practice of the 2030 UN sustainable development agenda.
International grassland conservation and protection trials continuously contribute towards valuable experiences and knowledge of grassland ecosystems on a global scale. An example is the Conservation Reserve Program (CRP), developed by the US federal government. Through rental payment and voluntary participation, the CRP proved to be beneficial for the conservation of grassland, control of soil erosion and even wildlife habitat conservation [10]. The Chinese Government has also implemented a series of grassland conservation projects since 2000 [11,12]. These projects are moving from macro-level planning towards local condition adaptation. The Beijing-Tianjin sandstorm source control project launched in 2000, targeted areas mainly located in the sandy and severely degraded dry grassland areas. Key measures included fencing, reseeding and the establishment of grazing bans. In 2003, the Returning Grazing Land to Grassland project was implemented which also focused on severely degraded grassland areas. Seasonal, semi-annual or year-round grazing suspension is implemented regionally. Since 2011, the Subsidy and Incentive System for Grassland Conservation (SISGC) was launched. This is the first all herder inclusive program for grassland conservation in China. Previous projects mainly focused on the degraded regions and provided a subsidy for those who sacrificed their benefits for local grassland conservation. The SISGC would also reward herders who maintain the balance of grassland carrying capacity and livestock, even in a good and undegraded grassland region. Herders switch from passively accepting the implementation of grassland ecological measures to actively carrying out grass-livestock balance and grazing prohibitions followed by corresponding rewards. This positive feedback chain would further enhance the cognition and consciousness of herdsmen and reduce the threat caused by overuse of grasslands under good precipitation conditions.
The first round of SISGC was from 2011 to 2015 and consisted of four regulations. (1) For grasslands that are seriously degraded and not suitable for grazing, a grazing ban is implemented. Central government grants a subsidy to herders of 14.44 USD/ha (1 USD = 6.2284 CNY, in 2015). (2) A grass-livestock balancing reward: for other available grasslands which show a reasonable carrying capacity, and checked by local institutions, the central government will reward the herders who have not overloaded grazing with the standard 3.61 USD/ha. (3) Herders' productivity subsidy: the subsidies for the improvement of livestock variety and herders' production. This subsidy includes 80.28 USD per household for about 2 million herdsmen in eight provinces and autonomous regions. (4) A performance appraisal regulation to improve verification mechanisms. The second round of SISGC started in 2016. The reward and subsidy standard should be adjusted each round by central government, and the specific regions and local accordance subsidy need to be adjusted every year by local government according to the changing condition of grassland. Continuous observation and real-time monitoring of grassland change and pixel-specific analysis are important for an optimal implementation of policy and to balance the interests of relevant stakeholders.
To identify the spatio-temporal grassland cover change, the traditional and effective method is based on in situ surveys, which can provide creditable and precise results but, considering the wide expanses of grassland, this also goes along with costly and timeconsuming intensive fieldwork, limited sites and lagging data processing [13]. With the improvement of spatio-temporal resolution of remote sensing (RS) and data processing technologies, mapping and monitoring of grassland cover change becomes more efficient and is often combined with in situ data for validation [14]. Many researches implement vegetation indexes (VIs) extracted from spectral sensors for grassland monitoring, including normalized difference vegetation index (NDVI) [15], enhanced vegetation index (EVI) [16], net primary productivity (NPP) and surface vegetation coverage [17], etc. There are also synthetic researches based on the analysis of multiple VIs addressing the importance of multiple perspectives [18] through analysis of variations in species composition, grassland desertification and above ground biomass of grassland degradation in semiarid grasslands. Although the VIs are abundant and flexible, and the information on grassland degradation and condition classification show large differences [19], grassland monitoring activities on a regional and global level are gradually providing more specific information to support sustainable grassland management.
In recent decades, many studies have concentrated on the monitoring of vegetation dynamics in IMAR. Some results indicate that vegetation growth has been poor since the beginning of the 21st century and that vegetation cover is lower in comparison to the 1980s or 1990s [20]. Other contradictory results have shown an increase of vegetation cover or productivity since the year 2000, possibly due to governmental decisions in the field of grassland management and conservation [21,22]. Historic change is always under debate. Considering that SISGC is an annual-based adjustment project, the grassland development analysis during the implementation period seems more practical for local stakeholders. In this research, we attempt to use open access RS data to monitor changes in grassland coverage at pixel level, combining Fractional Vegetation Cover (FVC) and NPP to qualitatively and quantitatively explore the annual grassland change in IMAR. The main objectives of this research are to understand grassland development and its trend since SISGC implementation in 2011. The results can help to support sustainable grassland management strategies for local governments when adjusting to SISGC on an annual basis and provide management references for other global grassland areas.

Study Area
The IMAR, located at 37 • 24 N to 53 • 23 N and 97 • 12 E to 126 • 04 E (Figure 1), is dominated by a temperate continental monsoon climate. The annual precipitation is in the range of 100 mm to 500 mm, with a majority of precipitation occurring during the growing season between July and September. The mean annual temperature ranges from −2 • C to 6 • C. The frost-free period is 70 to 160 days [23]. Grasslands are the main land cover type in IMAR (21.7% of total China's natural grassland) [24]. Soil types in the region consist of chernozem, chestnut, and calcic brown. Meadow steppe, typical steppe, and desert steppe [25] are spread from east to west, and the production of grassland also gradually decreases accordingly. The average yield of various pastures is 2863.57 kg/ha to 344.83 kg/ha, and the animal carrying capacity is between 104.95 ha and 1589.21 ha per sheep unit [26]. The main grassland types include Stipa krylovii, Leymus chinensis and Artemisia spp. [27]. Animal husbandry is the basic industry of IMAR and occupies a key position in the country. The output of lamb, milk, cashmere and the production capacity of grassland animal husbandry of IMAR are all ranked first in China [26]. Grasslands are closely related to the production, life and cultural development of herders. Over 90% of herdsmen's income comes from grasslands [26]. In 2019, the per capita disposable income of permanent residents in rural pastoral areas was 2454 USD, about half of the per capita disposable income of all residents in IMAR (4906 USD) [29].

Data Sources
Vegetation index production, meteorological data and the grassland conservation and protection policy in IMAR are basic data sets in this research (Table 1). NDVI and NPP were obtained from moderate-resolution imaging spectroradiometer (MODIS). IMAR is covered by seven tiles of global MODIS data. We collected 1392 tiles of MOD13A1 NDVI data in a temporal resolution of 16 days and spatial resolution of 500 m, and 63 annual MOD17A3HGF NPP data products with 500 m resolution [30] covering the period from 2011 to 2019. MOD13A1 NDVI data provides VIs, including NDVI and EVI. The algorithm for this product chooses the best available pixel value from the 16-day image composites. The criteria are low cloud cover, a low view angle, and the highest NDVI value. Annual NPP (MOD17A3HGF NPP) is derived from Net Photosynthesis (PSN), which is the difference of the Gross Primary Productivity (GPP) and the Maintenance Respiration (MR). We conducted data preprocessing including data format transformation, reprojection, mosaicking and resampling with the MODIS Reprojection Tool (MRT). To further reduce the influence of clouds, atmosphere and solar altitude angle, we transferred the 16-day NDVI composites to yearly composites using Maximum Value Compositions (MVC), extracted by the Inner Mongolia administrative boundary and thus obtained yearly NDVI products of IMAR from 2011 to 2019.
The Climate Prediction Center (CPC) Global Temperature and Precipitation meteorological dataset from 2011 to 2019 (Table 1), included global daily precipitation and daily maximum and minimum temperatures in a 0.50-degree spatial resolution. We transformed the NetCDF (network Common Data Format) daily value data from 2011 to 2019 into raster sets through raster calculation, projection and resampling. The resulting annual average precipitation and temperature products in a resolution of 500 m were finally used in conjunction with the processed NPP and NDVI data.

Methods
There are two key lines for grassland development analysis in this study. One is the indicator of FVC from a qualitative aspect, another is the indicator of NPP from a quantitative aspect. To reveal more information on impact factors for grassland productivity, we developed Human-induced NPP (HNPP) and Potential NPP (PNPP), which will be explained in more detail in Section 2.3.2. Focusing on these four indicators, we analyzed their trends and fluctuations based on a time series change analysis conducted from 2011 until 2019 and between 2018 and 2019. The results are finally discussed and put in context to local sustainable grassland management strategies, SISGC and impacts of climate change and human activities.

Fractional Vegetation Cover
FVC is the projected area of vegetation per unit of ground area and an important ecological parameter to evaluate grassland change and degradation [31]. It is assumed that the NDVI value of each pixel is composed of two parts, vegetation NDVI (NDVI v ) and soil NDVI (NDVI s ) (see Formula (1)). The linear mixture model [32] is used to estimate FVC. The model has been widely used and has shown to be efficient for global [33], national [17] and regional (IMAR) applications [34]. For instance, Tong et al. used the dimidiate pixel model to analysis FVC change from 1982 to 2010 in IMAR. The results have explained well the correlation with precipitation and the complex impact of human activities [20]. Furthermore, research based on unmanned aerial vehicle imagery in desert steppe of Alxa in the western part of IMAR, proved that the linear mixture model can obtain an accuracy of 83.1% [35]. In this model, NDVI is not a fixed value due to the influence of atmospheric conditions, humidity and sunlight. The different composition, spatial distribution and seasonal changes in vegetation growth also lead to a spatio-temporal variation of NDVI v . Due to the lack of large-area measured data, a confidence interval of the NDVI statistical histogram is usually used to obtain the maximum and minimum values and determine the vegetation index of the whole vegetation coverage as well as the bare area [36]. In this study, the cumulative percentages of 5% and 95% were selected as confidence intervals, which are widely used in northern China [20,37]. The cumulative percentage of less than 5% is the approximate pure soil cover vegetation index (NDVI s ), and the grassland area with a cumulative percentage greater than 95% is the pure vegetation cover vegetation index (NDVI v ). FVC were calculated annually from 2011 to 2019.

PNPP and HNPP
PNPP symbolizes potential vegetation conditions under the hypothesis that NPP is affected by climate change only. The Miami model [38] is empirical, derived through correlated mean annual temperature and precipitation values in conjunction with NPP and widely used for analyzing global productivity patterns. The model also acts as a baseline for comparisons [17]. The employed Miami model was used to simulate PNPP. The basic principles of the Miami Model are described by Formulas (2) and (3).
With NPP T = 3000/ 1 + e 1.315−0.119t where t is the annual mean temperature ( • C), and p is annual precipitation (mm). HNPP, reflecting the impact of human activities on ecosystems, can be quantified as the difference between PNPP and actual NPP (see Formula (4)) [39]. Positive HNPP represents negative impacts of human activities; whereas negative HNPP indicates that human activities promote NPP. The relative contributions of climate variation and human activities on NPP can be identified by comparing the slope of NPP, PNPP and HNPP [39,40].
2.3.3. Changing Trend Analysis of FVC, NPP, PNPP, and HNPP Linear regression can simulate the changing trend of each pixel in space, and reflect the changing trend of the entire space (see Formula (5)). If the slope is positive, it reflects an increasing trend, otherwise, it indicates a decreasing trend. The larger the slope, the more obvious the changing trend. Linear regression analysis was used to detect the spatialtemporal variation of vegetation indexes (including NPP, FVC, PNPP, and HNPP) from 2011 to 2019 on pixel level.
where X are the indexes (NPP, FVC, PNPP, or HNPP); n is the total of monitoring years; and i is the ith year.

Variation Analysis of FVC, NPP, PNPP, and HNPP
The coefficient of variation (CV) is used to measure the degree of variation of observations (see Formula (6)). In this study, the model is used to analyze the variation of VIs (NPP, FVC, PNPP, and HNPP) in the time series on pixel level. This reflects the discreteness and volatility of the data. Higher values of CV indicate a more discrete data distribution, higher volatility and drastic changes.
where σ is standard deviation; x is the average value.

Grassland Change Classification
With a focus on the indexes that reflect grassland conditions from a qualitative aspect (FVC) and quantitative aspect (NPP), more analysis in respect to their changing trend (X slope ) and variation (CV) over a time series is possible. Based on the analysis, we classified the grassland change trend in IMAR into four categories. If both NPP and FVC show increasing trends (X slope > 0), pixels are classified as synthetic increase; if NPP or FVC show a decreasing trend (X slope < 0), pixels are classified as decreasing trend, using Slope nf to represent the synthetic slope of NPP (Slope NPP ) and FVC (Slope FVC ) for each pixel. With the study region as a reference, we get the mean value of CV of the whole region (i.e., 0.20, in this study). Relevant variation degree would be classified into fluctuating (CV > 0.20) and persistent (CV < 0.20) categories. To highlight the fluctuation factor and reveal the potential threats in advance, the maximum value of CV for NPP and FVC in each pixel was compared (CV nf = Max (CV NPP , CV FVC )). Based on the above, the four following categories can be derived (see Table 2): (1) fluctuating increasing trend (Slope NPP > 0 and Slope FVC > 0; CV nf ≥ MeanCV nf ); (2) fluctuating decreasing trend (Slope NPP < 0 or Slope FVC < 0; CV nf ≥ MeanCV nf ); (3) persistent increasing trend (Slope NPP > 0 and Slope FVC > 0; CV nf < MeanCV nf ); (4) persistent decreasing trend (Slope NPP < 0 or Slope FVC < 0; CV nf < MeanCV nf ).

Grassland Change Matrix in the Last Two Years
Due to climate change and increase of extreme events, random and unpredictable factors are one of the major threats for pastoralists in grassland areas [41]. Therefore, it is essential to consider these factors to support sustainable grassland management strategies. In addition to the time series analyses in the above sections, we explored grassland change only within two adjacent years to reflect the random factor. Considering the yearly-based adjustment feature of SISGC, the latest years of 2018 and 2019 were set as an exploration. A spatial transfer matrix based on FVC and NPP results was designed.

Spatial Characteristics of Grassland
The FVC in IMAR is increasing from west to east (see Figure 2 left). The highest grasslands coverage is in the range of 40% to 60%, accounting for 24.83% of the total area. This area is mainly concentrated in the central part of IMAR, covering the western part of Hulunbuir and the central part of Xilin Gol. The region with less than 20% of grass vegetation covers 18.62% of IMAR and is concentrated in the west, including the western part of Xilin Gol and most area of Alxa, Bayan Nur, Bautou and Ulanqab. The average NPP of grasslands in IMAR shows an increasing, then decreasing trend from west to east (see Figure 2 right). The area with the highest NPP is concentrated to the west of Hulunbuir and Hinggan with mainly distributed meadow steppe and typical steppe. About 50.00% of grassland within the whole area has an average NPP in the range of 150 gCm −2 yr −1 to 300 gCm −2 yr −1 . However, the grassland production under 150 gCm −2 yr −1 is the second-highest category covering 37.08% of grassland in IMAR.

Spatial-Temporal Change of Grassland from 2011 to 2019
Both FVC and NPP in IMAR showed an increasing trend (see Table 3). The multi-year average FVC is 0.48 and the multi-year average NPP is 124.83 gCm −2 yr −1 . The trends of FVC and NPP in IMAR are positive. More than half of the regions showed an increasing trend in FVC and NPP, with proportions of 52.71% and 50.43%, respectively. The increasing trend of FVC is not apparent, with a mean slope value of 0.001. However, NPP reflects an obviously relevant increasing trend (mean value is 10.45). Due to better rain conditions, 65.99% PNPP of grassland shows a relevant insignificant increasing trend. The multi-year average of HNPP is 166.94 gCm −2 yr −1 , the positive value of HNPP reflecting a negative impact of human activities. Furthermore, HNPP has increased in 63.32% of grassland areas. Considering the fact that the actual NPP shows an increasing trend in the whole region of IMAR, the impact of multi-factors, including the positive effects of precipitation and grassland restoration, as well as disturbing effects from human activities, should be considered when assessing local grassland use conditions. In addition, both FVC and NPP of grassland show fluctuations from 2011 to 2019. Especially for FVC, the average CV for annual FVC is 0.22. Fluctuating FVC reflects the unpredictability of grassland condition and challenges sustainable grassland management.

Grassland Cover Change Trend Classification from 2011 to 2019
Combining FVC with NPP, the northern part of IMAR reflected an obviously decreasing trend, while the western part of IMAR fluctuates more strongly (see Figure 3). Around 33.69% of grasslands in IMAR are above the average coefficient of variation for the whole area and also show a decreasing trend. These areas are mainly located in the meadow steppe area in the west of Hulunbuir, Hinggan and east of Xilin Gol. Under good climatic conditions, the persistent decreasing trend reveals a constant human pressure on these meadow and typical steppe areas. With a lower precipitation rate, the western area is more sensitive to climate factors and reflects both increasing (19.21%) and decreasing trends (25.69%). Areas with an increasing trend are concentrated in the west and southeast area, covering the main part of Ordos, Hohhot, Bayan Nur, Baotou and Ulanqab. With mainly semi-desert grassland in these areas, the fluctuating increasing trend of 19.21% and persistent increasing trend of 21.42% reflects an overall improvement of grassland condition both in vegetation coverage and production.

Grassland Change Matrix from 2018 to 2019
Compared to 2018, only 15.03% of grassland shows an increase in both NPP and FVC in 2019, while 41.60% of the grassland area shows a decrease in NPP and FVC. The decreasing areas are mainly distributed in typical steppe zones including areas north of Xilin Gol, and west of Hulunbair and Ordos (see Figure 4). According to the meteorological data, annual average precipitation in IMAR was 346 mm in 2019, which is comparable to the precipitation of 369 mm in 2018. Nevertheless, there is high spatial variability throughout the year. In Xilin Gol league, where the majority of grasslands are typical grasslands, the precipitation was 261 mm during the grassland growing season in 2019, which was 89 mm less than the previous year (percentage of precipitation anomaly, 34%). The average temperature was 14.5 • C, 1.2 • C higher than in the previous years. Less precipitation and higher temperatures in mid-season are the main factors affecting the formation and growth of forage production in Xilin Gol League in 2019 [42]. This phenomenon reveals the fact that areas with better precipitation conditions and grass production show a more obvious amplitude of vegetation change. On the other hand, the unpredictable annual meteorological conditions are a challenge for a precise livestock-forge balance policy, as well as effective implementation of SISGC. These projects need to set rules for the size of households' livestock according to the grassland condition of the exact current year.

Climate and Human Activities Impact
Annual mean NPP, PNPP and HNPP show an overall increasing trend (see Figure 5). PNPP, affected by meteorological factors, shows an obvious trend with a slope of 2.4432, reflecting a positive impact of precipitation and temperature for the grassland system. However, HNPP shows a higher increasing trend (slope = 2.658). According to the statistical data, the number of livestock in IMAR increased from 68.06 million in 2011 to 71.92 million in 2019 [43]. Although the actual NPP in the past decade shows a gradually increasing trend (slope = 0.3115), pressure from human activities also exists.

Discussion and Conclusions
Through the combination of MODIS data with up to date meteorological data, and employing the indicators of FVC and NPP in both vegetation's qualitative and quantitative aspects, we analyzed the grassland development in IMAR since 2011, when the comprehensive grassland conservation program (SISGC) was launched.

Discussion and Conclusions
Through the combination of MODIS data with up to date meteorological data, and employing the indicators of FVC and NPP in both vegetation's qualitative and quantitative aspects, we analyzed the grassland development in IMAR since 2011, when the comprehensive grassland conservation program (SISGC) was launched.

The Positive Effect of SISGC, Especially in the Western Semi-Desert Grassland Area
Many researchers have tried to reveal and explore the change of grassland development, especially the degradation condition after the implementation of environmental conservation projects, such as the improvement of height, coverage, and biomass of natural grassland. However, results are often spatially and temporally inconsistent. Some studies showed that grassland condition has recovered to some extent in the past decades [24], while others claim that overall grassland condition has continued to deteriorate and grassland degradation is still the major ecological issue [25]. However, recent researches revealed that the year 2011 can be seen as a turning point for grassland dynamic trends in IMAR. By empirically analyzing the NDVI developing trend based on data of 52 counties across a 15-year timespan, results showed that there was a deteriorating trend between 2001 and 2010, followed by an improving trend since 2011 [44]. Studies of vegetation growth dynamics based on EVI since 2000 also reflected an obvious turning point in 2011 ( Figure 3 in the reference research) showing a significant standard deviation in 2012 [45]. NDVI studies based on datasets of NASA's Earth Science Data Systems Programs from 2001 to 2014 also showed a visible increasing trend since 2011 ( Figure 2 in the reference research) [24].
Based on these studies, the SISGC implementation in 2011 seems to have positively contributed to local grassland conservation in IMAR. The continuous effects in the past ten years and the potential yearly adjustment for the SISGC seem more practically essential. We chose to use FVC and NPP as two main indicators for detection of grassland change in IMAR, using FVC to visualize grassland quantity and NPP to assess product quality. These two aspects are in line with the national standards of Peoples' Republic of China "parameters for degradation, sandification and salification of rangelands (GB19377-2003)" [46]. The indicators proved to be effective for appropriate grassland degradation monitoring in IMAR [17]. Our study results reflected the positive effects of grassland recovery since 2011. Both FVC and NPP showed an increasing trend (Section 3.1.2). Su's research based on field survey data in IMAR also indicated that the aboveground NPP of different grassland types increased significantly between 2011 to 2013, in line with the progress of the national ecological program [47].
Besides, our study shows that the positive effect since 2011 is spatially unbalanced. The grassland change trend classification ( Figure 3) shows that most green pixels (increase) are concentrated in the western part of IMAR's grassland area, while most of the red pixels (decrease) are distributed in the upper eastern part. The grassland types in IMAR are desert steppe, typical steppe, and meadow steppe [25] following a precipitation gradient from west to east. Thus, the increasing trend or the positive effect is more obvious in the western part with desert steppe than in the eastern part with typical and meadow steppe. For example, in Baotou, where semi-desert grasslands are characteristic, the combined NPP and FVC classification system indicates that 86.03% of grassland shows a fluctuating or a persistent increasing trend. However, in Hulunbuir, characterized by meadow and typical steppe, only 16.99% of grassland shows an increasing trend. The possible reason might be that semi-desert grassland with poor ecological conditions has a higher potential for improvement. However, the significantly lower increasing trend in typical and meadow grassland possibly reflects that these areas receive less attention.

Potential Threats and Possible Reasons: Social and Natural Effects on Grassland Change
From 2011 to 2019, the annual average precipitation in IMAR increased from 110.39 mm to 120.00 mm; The annual average temperature also increased from 3.78 • C to 5.02 • C (See supplement file, Table S1). These positive climatic factors for grassland production result in an increase of PNPP (See Formula (3), Section 2.3.2). However, HNPP also shows an increasing trend (Section 3.1.2), implicating the intensity of human activities. Spatially, the annual average precipitation decreases from east to west of IMAR (See supplement file, Figure S1). The analysis of spatial-temporal precipitation change, based on available data from 2010 to 2015 in a 1 km spatial resolution (See supplement file, Figure S2), shows an increasing trend of precipitation in the eastern part, where typical and meadow steppe dominated. Improving natural climatic conditions will most probably improve grassland production in these areas. However, the grassland classification results show a decreasing trend in most of these areas (Figure 3). The contrasting results are alarming for the implementation of the grassland conservation program. The additional increase of PNPP and NPP in the whole region might cover up the negative effects, especially in areas under better precipitation conditions.
Considering SISGC implementation, Hu's field survey in 2015 revealed that SISGC had limited impact on herders' livestock production and grazing management, especially for herders in good grassland regions [48]. The compensation is much lower compared to possible market benefits and seems to be less attractive, although the subsidy and reward from central government on SISGC in IMAR are approximately 3.2 billion USD for the first round (2011-2015) [48]. Statistical data shows that stocking livestock in IMAR increased from 68.06 million to 91.92 million between 2011 and 2019 [43]. A household survey research also revealed that the overall stocking rate marginally increased significantly [49]. Herders were found to be more sensitive to market signals rather than the subsidy [48]. Local government encountered increasing supervision costs when the producer price of live sheep soared, which prompted local herders to gain more economic benefits by increasing illegal grazing [24]. The difficult tasks for local government may also lead to the malfunction of SISGC implementation to some extent. Some researchers have pointed out that the determination of a maximum stocking rate and various grazing bans have a less feasible effect under the spatial heterogeneity of grasslands and climate variability [50], and the effectiveness of conservation strategies in arid and semi-arid areas need to be discussed thoroughly [51]. Integrated management considering climate conditions, geographical characteristics and socioeconomic factors is prospected to be beneficial for regional grassland conservation [52].
The grassland ecosystems are ecologically vulnerable and sensitive to climate change [53]. Vegetation phenology in temperate grasslands is highly sensitive [54]. Gong et al. observed that the length of the growing season in grassland regions of IMAR has increased by 10.86 days from 2002 to 2014 [21]. Ren et al. also revealed that during 2000 to 2015, the start of the growing season in typical steppe and desert steppe has significantly advanced by 2.2 days and 10.6 days, and the end of the growing season in desert steppe has also significantly advanced by 6.8 days [54]. Combining the phenology changing results under the background of unavoidable global climate change could help to guide grazing systems and to develop more efficient policy frameworks for grassland protection.

The Importance of a Meteorological Prediction System for Local Stakeholders
The difference between 2018 and 2019 (Section 3.2.2) reveals the impact of unpredictable meteorological conditions. Although conservation programs and human efforts can make a positive difference in the grassland region, the sensitivity of vulnerable grassland systems to environmental factors still acts as the main driver. This also sets obstacles for a successful implementation of grassland conservation programs, balancing the capacity of grassland and livestock on a yearly basis. Under these circumstances, an earlier and up-to-date grassland growing prediction system is of great importance, not only for the equitable implementation of conservation projects but also for efficient grassland use strategies of herders. Local government could develop improved planning for implementation of grazing bans or forage-livestock balancing measures in the region; herders could decide whether natural grassland needs to be reserved or could rationally be used during a predictable good production year. A similar project has already been employed in Mongolia, where natural hazards are the main issues in grassland systems [55], through the use of near real-time satellite-based climate data, computer simulation and modelling. Mongolia's livestock early warning system, accessible to all people, monitors and forecasts livestock forage conditions to help pastoralists and decision makers make timely decisions while facing natural hazards. This system can be a good reference for IMAR and other pastoral areas in the world in need of balanced grassland management.
Some issues still need to be addressed beyond this research. Firstly, in the binary pixel model for calculation of FVC, the impact of soils on vegetation signatures is an important factor which often causes errors. The parameters of NDVI s and NDVI v are always under discussion and the selection principles vary throughout the research, especially concerning the determination of NDVI s . Basically, the NDVI s vary from −0.1 to 0.2 [56]. The estimation study of FVC in Northeast China pointed out that NDVI s ranges between 0.07 and 0.22 [37]. Research on FVC estimation in arid and semi-arid environments recommended NDVI s and NDVI v in grassland are 0.045 and 0.593, respectively [57]. Estimation research of FVC in IMAR by MODIS data in 2013 obtained the result that NDVI s and NDVI v in grassland are 0.118 and 0.806, while with in situ measurement the results are 0.203 and 0.891 [34]. In this research, the mean value of NDVI s and NDVI v from 2011 to 2019 are 0.07 (CV = 0.0%) and 0.89 (CV = 0.7%). The value of NDVI s shows in the normal ranges as mentioned above, but is lower than values obtained through in situ data in 2013. The value of NDVI v is accordance with the MODIS and in situ research result in IMAR [34]. Secondly, although the study has qualitatively and quantitatively analyzed the change of grassland based on FVC and NPP indicators, grassland species succession during the process of grassland degradation may also contain the increase of mixed shrub types, which will also lead to an increase in both FVC and NPP. As such it will be important to develop a method to differentiate between shrubs and grassland. With the development of LiDAR, hyperspectral data and unmanned aerial vehicles, attempts are arising in discriminate specific grassland species [58]. Although limited to regional and short-term studies due to limited data access and high-cost data, this discrimination process should be incorporated into future research to provide more precise grassland change results and reliable support for local decision making. Besides, livestock always act as the basic conjunction between human society and environment. Individual animal health and performance change (i.e., body condition, reproduction) need more attention for a thorough understanding of the effectiveness, impaction and sustainability of grassland conservation programs.
Author Contributions: Conceptualization, M.K. and Z.X.; methodology and analysis, Z.X.; writingoriginal draft preparation, Z.X.; writing-review and editing, D.W., M.K. and Z.X.; All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.