Eco-Environment Evaluation of Grassland Based on Remote Sensing Ecological Index: A Case in Hulunbuir Area, China

This research is based on Landsat5 TM, Landsat8 OLI/TIRS remote sensing data using RSEI model to analyze and monitor the ecological environment and its temporal and spatial changes in the forest-grass transition zone in Northeast China from 2004 to 2019. The change characteristics of the ecological environment of different types of land cover types are monitored by RSEI method, and the response of different land cover types to natural factors such as precipitation and temperature is analyzed at the same time. The distribution of RSEI in the study area presents the characteristics of high in the east and low in the west. The eastern mountainous area is densely covered with woodland, which is the area with the best ecological environment quality in the study area. The grassland in the western plain and the saline-alkali land around the river are the areas with poor ecological environment in the study area. Climate, precipitation, topography and other natural elements work to-gether to form the quality of the ecological environment in the study area roughly bounded by 120˚E. In years with poor natural conditions, this divid-ing line will have a clear eastward shifting trend, especially in the northern part of the study area. The spatial distribution of RSEI in the study area has a high degree of spatial autocorrelation, and Global Moran’s I has been above 0.8 over the years. In terms of temporal changes in ecological conditions, the ecological environment in the study area was basically stable from 2004 to


Introduction
Grassland resources are the largest ecological system on the land of China, and play an important role and value in developing animal husbandry, protecting biodiversity, maintaining water and soil, and maintaining ecological balance. At present, the ecological quality of China's grassland is facing serious risks of degradation and desertification [1] [2] [3]. Many studies have shown that overgrazing has become one of the main causes of grassland degradation [4] [5] [6] [7] [8]. From 2000 to 2012, the total amount of livestock in Hulunbuir grassland and the overload rate of livestock continued to rise, and the overload rate of livestock in 2012 exceeded 30% [9]. In the context of the increasingly serious grazing problem in the Hulunbuir grassland area, it is necessary to evaluate and monitor the ecological conditions of the Hulunbuir grassland area. At present, the commonly used methods for grassland ecological evaluation and monitoring include NDVI change monitoring based on remote sensing data [10], vegetation change monitoring based on drone data [11], and comprehensive evaluation based on multiple data sources [12] and so on.
Remote sensing data has the characteristics of rich information, fast acquisition speed, periodicity, etc., and is widely used in the field of ecological environment monitoring (LI, 2001). Many remote sensing ecological indexes have been proposed to reflect the surface coverage and ecological quality, such as: NDVI, IBI [13], NDWI [14], SWDI [15], etc. Most of the above remote sensing indexes are extracted and studied for a single object and cannot comprehensively reflect the regional ecological environment. Some scholars have constructed a comprehensive remote sensing ecological index for a variety of objects and used multiple remote sensing ecological indicators to measure the quality of the regional ecological environment. Xu Hanqiu and others proposed a comprehensive remote sensing ecological index RSEI [16] that combines surface humidity, greenness, dryness, and temperature information. Shen Ge et al. proposed WDCI for wetland destruction [17]; Yang et al. added a parameter to measure the health of vegetation on the basis of RSEI, and proposed the CEEI index [18]. Among them, the RSEI model is widely used in regional ecological environment assessment and monitoring due to its strong comprehensiveness, simple data acquisition, and high degree of data visualization [19] [20] [21] [22].
In order to analyze the problem of the low contribution rate of the first principal component in the process, this paper adopts the RSEI model improved by Pan et al. [23] to quantitatively evaluate and monitor the characteristics of the temporal and spatial changes of the ecological environment quality in the Hulunbuir grassland area, and analyze the Hulunbuir grassland. The clustering cha-racteristics of the ecological environment of Hulunbuir grassland area and the change characteristics of the ecological environment of different land types in the grassland area are also analyzed. At the same time, the main reasons for the changes in ecological quality are analyzed in combination with precipitation and climatic factors.

Study Area
The study areais located in Hulunbuir City, Inner Mongolia Autonomous Region, Northeast China. The terrain is high in the west and low in the east, with an altitude of 600 -650 m. Taking the meridian of 120˚E as the boundary, the climate of the study area is a mid-temperate monsoon mixed forest climate to the east, and a mid-temperate monsoon forest grassland climate to the west. The climate here is characterized by large temperature difference between day and night, lack of heat, frost-free period, abundant sunshine, large difference in monthly precipitation, and precipitation is concentrated in July-August. Woodland and grassland are widely distributed in Hulunbuir City, with typical characteristics of forest-grass transition zone. Figure 1 is a schematic diagram of the ecological location and elevation of the study area.

Data and Processing
We selected 5 periods of remote sensing data covering the study area, including Landsat5 TM data in 2004, 2007, and 2011 and Landsat8 OLI/TIRS data in 2015 and 2019. Table 1 is the detailed information of remote sensing data. All data were download from the https://earthexplorer.usgs.gov/.
In order to use satellite radiation data to reflect surface information, remote  sensing data needs to be properly processed. Using the calibration parameters provided by metadata, the image gray values of reflection and thermal bands are converted into absolute radiance. Then, the Fast Line-of sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) model was used to perform atmospheric correction on the reflection and thermal bands to eliminate radiation errors caused by atmospheric influences.

Remote Sensing-Based Ecological Index (RSEI)
The Remote Sensing Ecological Index (RSEI) is a comprehensive ecological evaluation index proposed by Xu Hanqiu in 2013 that is completely based on remote sensing technology, indicators to obtain, without artificially set weights, and a high degree of visualization of the results. Xu et al. selected four commonly used ecological elements that can reflect the quality of the ecological environment and can be retrieved by remote sensing bands, including greenness (G), wetness (W), dryness (D) and temperature (T) [13].
Then, we selected four commonly used remote sensing ecological indexes to reflect the above ecological elements, namely The Normalized Difference Vegetation Index (NDVI), Wetness (WET), Normalized Difference Soil Index (NDSI), and Land Surface Temperature (LST). The four of them respectively represent G, W, D, T in Equation (1) Principal component analysis (PCA) is used to compress the obtained index data and extract their main information. The dimension of each index is different, so it is necessary to use min-max normalization to process the obtained data before PCA.
Combined with the situation in the study area, PC 1 cannot concentrate most of the information as a result of principal component analysis, and PC 2 needs to be supplemented. Using the eigenvalue contribution rate obtained by principal component analysis as the weight, RSEI 0 can be calculated using PC 1 and PC 2 for weighted summation [23]. Then, the RSEI 0 is normalized to obtain the Remote Sensing Ecological Index (RSEI), so as to classify the ecological environment. The calculation methods of NDVI, WET, NDSI, LST can be found in the references.

Spatial Autocorrelation Analysis
Spatial autocorrelation analysis is a global spatial autocorrelation analysis of elements in space, the purpose is to judge whether the research object is cluster mode, discrete mode or random mode. Spatial autocorrelation is usually expressed by Global Moran's I. The closer its value is to 1, the stronger the spatial autocorrelation [24]: where i x is the observation value; ij w is the spatial weight matrix, which is an n*n symmetric matrix. When the spatial elements are adjacent, ij w takes 1, and when they are not adjacent, ij w is 0.
By calculating the statistically significant Z score (multiple of standard deviation) and p value (probability), the two are judged one by one the null hypothesis, so as to decide whether to reject the null hypothesis. When p is very small, it means that the observation value has a small probability of a null hypothesis, that is, it is unlikely to produce randomness.
The Local Moran Index (LISA) can effectively reflect the correlation between the quality of each ecological unit in the study area. When the global correlation is weak or absent, this tool can be used to reflect the location of the study area; when the global correlation is strong, this tool can be used to reflect the presence or absence of spatial heterogeneity [25]: The meaning of each symbol is the same as equation (3).

Method of Accuracy Evaluation
We randomly sampled the RSEI data and each sub-indicator for each period 500 times, and obtained 2500 sampling points in total. The method of multiple linear regression analysis was used to evaluate the data accuracy of the 2500 sets of data to the RSEI value. The sampled value of RSEI was used as the dependent variable, and the sampled value of NDVI, WET, NDSI, and WET were used as independent variables to perform multiple linear regression analysis.

RSEI Data Accuracy Evaluation
No independent variables are excluded in the regression analysis, indicating that the four dependent variables and independent variables all have a significant correlation. According to multiple linear regression analysis, Equation (5) is obtained. The results show that: RSEI has a linear correlation with the four sub-indices. Among them, NDVI, WET and RSEI are positively correlated, NDSI, LST and RSEI are negatively correlated. RSEI has a significant positive correlation with ecological environment quality. Table 2 shows the eigenvalues and percentages of each index in each principal component obtained after principal component analysis. The sum of the contribution rates of PC1 and PC2 can reach more than 80%, which can meet the research requirements, so Pan's improved method is adopted in order to obtain more comprehensive information [23]. According to the RSEI value, the ecological environment of the study area was graded over the years.  Figure 2 shows the spatial distribution of the ecological environment based on the above-mentioned level standards.

Analysis on the Changes of Ecological Environment Quality in the Study Area
It can be seen from Figure 2 that the overall ecological environment quality of the study area is poor, and it has the basic characteristics of high in the east and low in the west. Lv5 and lv4 are mainly distributed in the woodlands and floodplains of the mountainous areas in the southeast. lv1 and lv2 are mainly distributed in the vast grasslands in the west and saline-alkali areas near the rivers in the middle. The boundary between the good and bad of the ecological environment is roughly around the 120˚E meridian. In a year of good ecological environment (2011), this boundary will move to the west, otherwise it will move to the east. Figure 3 shows the spatial distribution of the rise and fall of the ecological environment level in the study area in each period. It can be seen from Figure 3   that before 2015, the main areas where the ecological environment deteriorated were concentrated in the western grasslands of the study area and the cultivated land areas in the northeast, which were greatly affected by natural elements. As the ecological units with the best and worst ecological environment, floodplains and woodlands are areas with the most stable ecological environment. However, after 2015, large-scale deterioration of the ecological environment occurred in the woodlands in the southeast of the study area. According to high-precision remote sensing images, the main reason for the deterioration is deforestation and land reclamation.

Spatial Autocorrelation Analysis
In order to understand the spatial and temporal distribution pattern of the ecological environment in the study area, the study is divided into 1 km * 1 km grids. The RSEI average value in the grid is assigned to the grid to calculate the global Moran's I of the study area (Table 3). It can be seen from the results that all the images of the 5 phases have passed this important test. From 2004 to 2019, the P-value is very close to 0, and the Z-score value is large, indicating that the RSEI distribution in the study area has a significant spatial autocorrelation. During the study period, the highest global Moran's I was in 2015 when the ecological environment was the best. The change trend of this index is consistent with the change trend of the ecological environment quality. The above image is used to calculate the local Moran I to generate a LISA cluster map (Figure 4). There are five clustering modes in the LISA graph, which

Conclusion
We selected the remote sensing impact of Hulunbuir grassland in 2004,2008,2011,2015, and 2019, and adopted the RSEI model to quantitatively evaluate and classify its ecological environment. The results show that the ecological environment is generally high in the east and low in the west. It gradually got better before 2011, but it kept getting worse since then. The floodplain and the surrounding saline-alkali land are the areas with the most stable ecological environment quality in the study area. Then, we conducted a cluster analysis on the RSEI value of the study area, and the results showed that there is a strong spatial autocorrelation, and the deterioration of the ecological environment will have a high probability of affecting the surrounding ecological environment. RSEI can effectively show the quality of the ecological environment in the study area. In this study, the RSEI model was used to evaluate and monitor the 15-year ecological environment of the Hulunbuir grassland area, which verified a strong feasibility method for the evaluation and monitoring of the grassland ecological environment.