Assessment of Changes in Key Ecosystem Factors and Water Conservation with Remote Sensing in the Zoige

: As the largest alpine peat swamp wetland distribution area in the world, the Zoige has important ecological functions, including water conservation and biodiversity maintenance. In the past 20 years, the regional ecological protection and restoration measures continuously strengthened under the leadership of the local government have led to gradual improvements in the ecological environment of the region. In this study, multisource satellite remote-sensing image data were used to carry out quantitative monitoring and assessment of the main ecological elements (vegetation and water), as well as the regional leading ecosystem service function in the Zoige. Combined with local ecological protection management policies and measures, we analyzed the characteristics and effectiveness of ecological protection. We compared the ecosystem change trends of the Zoige reserve and the county, from 2001 to 2020, and found that the fractional vegetation cover (FVC) of Zoige county has increased at a rate of 0.25%/year. The growth rate was highest between 2015 and 2020, and the growth rate of FVC in the Zoige Wetland National Nature Reserve is approximately 1.89-fold that of the whole county. The water area also shows similar variation characteristics. On the whole, the water conservation capacity of the Zoige showed a signiﬁcant increase from 2001 to 2020. We used high-resolution satellite remote-sensing images to capture the details of land use changes brought about by local ecological protection policies and measures, and together with macroecological indicators, we reﬂected on the effectiveness of regional ecological protection measures. We observed that the ecological effects of nature reserves are more direct and rapid, and the amount of water conservation within the nature reserve is about 1 × 10 4 m 3 /km 2 higher than that of the surrounding grasslands. Satellite remote-sensing images can not only capture the multiscale change information of ecological indicators, such as vegetation and water, in a timely manner, but can also help us to identify the effectiveness of conservation measures by distinguishing and analyzing the causes of these changes.


Introduction
Wetlands represent the transitional zone between terrestrial and aquatic ecosystems, and possess a unique ecosystem and functional characteristics [1,2]. At present, changes in vegetation and water bodies are basic factors that reflect the status of and changes in wetland ecosystems, and are key indicators of wetland shrinkage and restoration monitoring based on satellite remote sensing [3,4].
The Zoige wetland is located on the eastern edge of the Qinghai Tibet Plateau and is listed as an "International Important Wetland" by UNESCO [5,6]. It is an important water conservation area in the upper reaches of the Yellow and Yangtze Rivers, with the largest plateau peat swamp in the world [7]. It is a key area for biodiversity protection both in China and worldwide [8]; it is also an important habitat for wildlife species, such as the black-necked cranes and white-tailed sea eagles, and is highly sensitive to climate change [9,10]. The Zoige Wetland National Nature Reserve was established in 1998 and covers an area of 1665.70 km 2 . The alpine swamp wetland ecosystems, black-necked cranes, and other rare animals are especially targeted for protective measures.
Owing to the joint impact of climate change and human activities, the Zoige wetland has been shrinking since the 1950s, resulting in the decline in the groundwater level, land degradation, reduction in biodiversity, and increased carbon emissions. Relevant research has mainly focused on regional climate change analysis [11][12][13], wetland area change monitoring [14], grassland degradation, land desertification [15][16][17], soil organic carbon change [18,19], and ecosystem services and health [20,21]. Most of these studies were published before 2016, and the conclusions were mainly based on the continuous degradation of the ecological environment. There is still a lack of analysis on the effectiveness of conservation measures.
Since the 1970s, China has successively launched a series of ecological protection and restoration projects, including the Natural Forest Protect Projection, the Grain for Green program, and the Returning Rangeland to Grassland, which have increased regional ecological functions [22,23]. Since 2003, ecological protection and restoration projects, such as grassland restoration and desertification land management, have significantly increased in the Zoige area [24][25][26], and a series of restoration and rectification works of nature reserves have been carried out in succession. In the past 30 years, the growth rate of regional desertification land has significantly slowed [17,25]. The project of returning grazing land to grassland in the county has had significant benefits: the vegetation height in the project area has elevated by 69.1% on average, compared with that in the nonproject area [26].
From the perspective of satellite remote sensing, we used the main elements of the ecosystem (vegetation, water bodies, etc.) and the service function of water conservation, as indicators to analyze the changing characteristics of the ecological environment in Zoige county and nature reserves [27]. Simultaneously, combined with high-resolution remotesensing images and field investigations, we attempted to link the direct image evidence of ecological restoration with the above-mentioned remote-sensing indicators of ecological change, thus further elaborating the obvious ecological benefits brought about by ecological protection and supervision.

Study Area
Zoige county (102 • 08-103 • 39 E, 32 • 56-34 • N) is located in the northern part of the Northwest Sichuan Plateau, which is the core area of the Zoige wetland. The topography is low in the middle and high in the surrounding areas. It belongs to the Aba Tibetan and Qiang Autonomous Prefecture of Sichuan province, and has a total area of 1.04 × 10 4 km 2 and an average altitude of 3471 m ( Figure 1). The regional altitude difference is largely attributed to the depth of the river. The study area has a plateau cold temperate humid monsoon climate, with an annual average temperature of 1.1 • C, and annual precipitation of 650 mm [28]. However, in recent years, there has been an obvious trend of warming and drying [29].

Vegetation Data Acquisition and Calculation Method
Fractional vegetation cover (FVC) reflects the vegetation coverage of a certain area, and is an important indicator for measuring the quality of vegetation and ecosystem change [30]. Based on the remote-sensing cloud platform of Google Earth Engine and the Moderate Resolution Imaging Spectroradiometer (MODIS) 13Q1 normalized difference vegetation index (NDVI) product (https://ladsweb.modaps.eosdis.nasa.gov, (accessed on 5 June 2021)), we used the pixel dichotomy model to calculate the vegetation coverage [31,32]. The MODIS Diversity 2022, 14, 552 3 of 13 13Q1 NDVI product has a temporal resolution of 16 days with 23 cycles of data per year. In this study, FVC data for 23 periods were generated annually, and then aggregated into the average annual FVC. The specific formula of FVC is as follows: where NDV I v and NDV I s are the NDVI pixel values of pure vegetation and bare land, respectively. With the NDVI value statistics, we determined the cumulative frequency of 1% NDV I s and cumulative frequency of 99% NDV I v [33].

Vegetation Data Acquisition and Calculation Method
Fractional vegetation cover (FVC) reflects the vegetation coverage of a certain area, and is an important indicator for measuring the quality of vegetation and ecosystem change [30]. Based on the remote-sensing cloud platform of Google Earth Engine and the Moderate Resolution Imaging Spectroradiometer (MODIS) 13Q1 normalized difference vegetation index (NDVI) product (https://ladsweb.modaps.eosdis.nasa.gov, (accessed on 5 June 2021)), we used the pixel dichotomy model to calculate the vegetation coverage [31,32]. The MODIS 13Q1 NDVI product has a temporal resolution of 16 days with 23 cycles of data per year. In this study, FVC data for 23 periods were generated annually, and then aggregated into the average annual FVC. The specific formula of FVC is as follows: where and are the NDVI pixel values of pure vegetation and bare land, respectively. With the NDVI value statistics, we determined the cumulative frequency of 1% and cumulative frequency of 99% [33].

Water Data and Acquisition Methods
Water information in this study mainly included open water surfaces, such as rivers and lakes. Using information based on the Google Earth Engine platform, long time-series Landsat remote-sensing image data (https://ladsweb.modaps.eosdis.nasa.gov, (accessed on 24 May 2021)) corresponding to the study area and study period were directly called through the database, from which the vegetation, water, and dry naked index were calcu-

Water Data and Acquisition Methods
Water information in this study mainly included open water surfaces, such as rivers and lakes. Using information based on the Google Earth Engine platform, long time-series Landsat remote-sensing image data (https://ladsweb.modaps.eosdis.nasa.gov, (accessed on 24 May 2021)) corresponding to the study area and study period were directly called through the database, from which the vegetation, water, and dry naked index were calculated [34]. Each index is calculated as follows: where NDVI is the normalized vegetation index, MNDWI is the improved normalized difference water body index, AWEI is the automatic water body extraction index, DBSI is the dry bare index, ρ Blue is the blue band, ρ Green is the green band, ρ N IR is the near-infrared band, ρ red is the red band, ρ SW IR1 is the short-wave infrared band 1, and ρ SW IR2 is the short-wave infrared band 2. Secondly, the regional water body area was extracted by combining the remote-sensing water body index and the basin relative elevation setting, and the non-water-body information was filtered by combining the vegetation index and drought exposure index [34]. Among them, 30 m spatial resolution digital elevation model (DEM) products came from the geospatial data cloud platform (http://www.gsclod.cn, (accessed on 24 May 2021)).
Finally, the annual data product of flood frequency (the ratio of the number of times a place is submerged to the number of effective observation) was constructed using annual multiperiod water body data to determine the flood frequency [34,35]. Then, the human-computer interaction method was used to determine the data threshold of flooding frequency, which could maximize the avoidance of seasonal water and noise (the threshold in the study is 70%), and obtain the stable water distribution of the study area in the current year.

Ecosystem Function Assessment
Water conservation means that ecosystems (such as forests, grassland, etc.) intercept infiltrate, accumulate precipitation through their unique structures and water interactions, and control water flow and the water cycle through evapotranspiration. We used the water balance equation to calculate water conservation capacity [36]: where TQ is water conservation (m 3 ), P i is rainfall (mm), R i is surface runoff (mm), ET i is evapotranspiration (mm), A i is the type i ecosystem area (km 2 ), i is the type i ecosystem in the study area, and j is the number of ecosystem types in the study area. Rainfall data were obtained from a TerraClimate dataset with a spatial resolution of 1 km (https://climate.northwestknowledge.net/TERRACLIMATE/index, (accessed on 10 June 2021)). Evapotranspiration data were the result of the MODIS MOD16A3 product, with a spatial resolution of 1 km (https://ladsweb.modaps.eosdis.nasa.gov, (accessed on 24 June 2021))). The surface runoff was obtained by multiplying the rainfall by the surface runoff coefficient, and the surface runoff coefficient was obtained from the literature [36]. The ecosystem area was derived from data comprising remote-sensing surveys and assessments of the national ecological status, which mainly included forest, shrub, grassland, garden, and wetland types [37].
To avoid the influence of the interannual precipitation fluctuations, based on the calculation of annual water conservation, we took the average every five years to represent the regional water conservation capacity in different periods. Then, we comprehensively analyzed the change characteristics of water conservation function over the past 20 years.

Change Trend Analysis Method
For trend analysis, we used the Theil-Sen median method, a trend analysis method for nonparametric statistics. The calculation method of the Sen trend degree is as follows [38]: where x j and x i are time series data, a β greater than 0 indicates an increasing trend in the time series, and a β of less than 0 indicates a decreasing trend in the time series. The Mann-Kendall trend test is a nonparametric statistical test method for which measured values do not need to follow a normal distribution. The trend is not required to be linear, and is not affected by missing values and outliers [39,40]. The statistical test method used is as follows: where sgn is a symbolic function, x j and x i are sequential data sets, and n is the length of the data samples. Z c followed a standard normal distribution. If Z c > Z (1−α)/2 , there is a significant trend change. Z (1−α)/2 is the corresponding value of the standard normal function distribution table at confidence level α. The confidence level α was set at 0.05.

Vegetation Change Characteristics
The spatial distribution of fractional vegetation cover in Zoige county is shown in Figure 2. The FVC of the county in 2020 was 57.56%. Areas with large FVC (>70%) were mainly distributed in the eastern Zoige forest ecosystem, while areas of low FVC (<20%) were distributed in the northwest of the county where sandy, urban, and rural land were concentrated.
be linear, and is not affected by missing values and outliers [39,40]. The statistical test method used is as follows: where sgn is a symbolic function, and are sequential data sets, and n is the length of the data samples.
followed a standard normal distribution. If ⁄ is the corresponding value of the standard normal function distribution table at confidence level α. The confidence level α was set at 0.05.

Vegetation Change Characteristics
The spatial distribution of fractional vegetation cover in Zoige county is shown in Figure 2. The FVC of the county in 2020 was 57.56%. Areas with large FVC (>70%) were mainly distributed in the eastern Zoige forest ecosystem, while areas of low FVC (<20%) were distributed in the northwest of the county where sandy, urban, and rural land were concentrated.
In 2020, the average FVC of the nature reserve was 53.22%. Because the terrain was relatively flat, the grassland in the area was evenly distributed. Although the average vegetation coverage in the nature reserve was 4% lower than that of the whole county, the vegetation coverage of grassland ecosystems inside the reserve was 1% higher than that outside of the reserve. In 2020, the average FVC of the nature reserve was 53.22%. Because the terrain was relatively flat, the grassland in the area was evenly distributed. Although the average vegetation coverage in the nature reserve was 4% lower than that of the whole county, the vegetation coverage of grassland ecosystems inside the reserve was 1% higher than that outside of the reserve.
Between 2001 and 2020, the FVC in Zoige demonstrated an overall upward trend, with an increase of 5.38% in 2020 compared with 2001. In particular, the rate of increase from 2015 to 2020 was significant compared with previous years (an increase of 1.35%). The average annual increase in FVC in this period was approximately 1.64-fold higher than that of the previous 15 years (2001-2015; Figure 3). demonstrated an increasing trend, consistent with the changing trend of the whole county, and reflecting the overall improvement in vegetation growth. The mean FVC in the protected area was smaller than the average FVC of the whole county, mainly due to the existence of a forest ecosystem in the east of the county. However, from 2015 to 2020, the FVC in the reserve significantly improved, and was approximately 1.35-fold larger than the increase of the county value. This suggests stronger ecological resilience in the reserve. From 2001 to 2020, the FVC in Zoige county increased at an annual rate of 0.25%; the FVC of 95.23% of the county increased, of which 63.6% of the area experienced significant increases (Figure 4). The rapidly increasing area was concentrated in the east, where forest and grassland met, and the nature reserve is located. Across the county, the area of FVC decrease accounted for 4.5%, mainly distributed in the northwest, where urban and rural land was concentrated. By comparing high-resolution remote-sensing images, some illegal facilities in the reserve were demolished ( Figure 5).
From 2015 to 2020, the rate of FVC increase escalated in both Zoige county and the reserve. The growth rate of the reserve (0.68% per year) was greater than that of the county as a whole (0.36%/a), reflecting the positive effects of targeted vegetation protection.  It can be seen that the FVC of the Zoige Wetland National Nature Reserve also demonstrated an increasing trend, consistent with the changing trend of the whole county, and reflecting the overall improvement in vegetation growth. The mean FVC in the protected area was smaller than the average FVC of the whole county, mainly due to the existence of a forest ecosystem in the east of the county. However, from 2015 to 2020, the FVC in the reserve significantly improved, and was approximately 1.35-fold larger than the increase of the county value. This suggests stronger ecological resilience in the reserve.
From 2001 to 2020, the FVC in Zoige county increased at an annual rate of 0.25%; the FVC of 95.23% of the county increased, of which 63.6% of the area experienced significant increases ( Figure 4). The rapidly increasing area was concentrated in the east, where forest and grassland met, and the nature reserve is located. Across the county, the area of FVC decrease accounted for 4.5%, mainly distributed in the northwest, where urban and rural land was concentrated. By comparing high-resolution remote-sensing images, some illegal facilities in the reserve were demolished ( Figure 5).

Water Body Change Characteristics
Remote-sensing data revealed that the water area of Zoige county in 2020 was 71.8 km 2 , and mainly comprising the Yellow, Baihe, and Heihe Rivers in the west, along with the Huahu and Cuorewajian Lakes, and surrounding water bodies in the northwest.
The average water area of Zoige county from 2001 to 2020 was approximately 42.77 km 2 . It exhibited an increasing trend. The water body area was lowest in 2002 (29.7 km 2 ). In 2020, the water body area was approximately 32.66 km 2 larger than that in 2001, an increase of 83.44% (Figure 6). By comparing high-resolution remote-sensing images, we noted that the area of some water bodies in Zoige county has surged (Figure 7).
Viewed by stages, the water body area in Zoige county was basically stable from 2001 to 2015, and continued to increase after 2015. From 2015 to 2020, the county's water body area rapidly increased, reaching nearly 1.83 times the level of early 2001. From 2015 to 2020, the rate of FVC increase escalated in both Zoige county and the reserve. The growth rate of the reserve (0.68% per year) was greater than that of the county as a whole (0.36%/a), reflecting the positive effects of targeted vegetation protection.

Water Body Change Characteristics
Remote-sensing data revealed that the water area of Zoige county in 2020 was 71.8 km 2 , and mainly comprising the Yellow, Baihe, and Heihe Rivers in the west, along with the Huahu and Cuorewajian Lakes, and surrounding water bodies in the northwest.
The average water area of Zoige county from 2001 to 2020 was approximately 42.77 km 2 . It exhibited an increasing trend. The water body area was lowest in 2002 (29.7 km 2 ). In 2020, the water body area was approximately 32.66 km 2 larger than that in 2001, an increase of 83.44% (Figure 6). By comparing high-resolution remote-sensing images, we noted that the area of some water bodies in Zoige county has surged (Figure 7).

Figure 5.
Comparison of remote-sensing images of demolition points of typical illegal facilities (high-resolution remote-sensing images from Chinese GF-1 and GF-2 satellites).

Water Body Change Characteristics
Remote-sensing data revealed that the water area of Zoige county in 2020 was 71.8 km 2 , and mainly comprising the Yellow, Baihe, and Heihe Rivers in the west, along with the Huahu and Cuorewajian Lakes, and surrounding water bodies in the northwest.
The average water area of Zoige county from 2001 to 2020 was approximately 42.77 km 2 . It exhibited an increasing trend. The water body area was lowest in 2002 (29.7 km 2 ). In 2020, the water body area was approximately 32.66 km 2 larger than that in 2001, an increase of 83.44% (Figure 6). By comparing high-resolution remote-sensing images, we noted that the area of some water bodies in Zoige county has surged (Figure 7).
Viewed by stages, the water body area in Zoige county was basically stable from 2001 to 2015, and continued to increase after 2015. From 2015 to 2020, the county's water body area rapidly increased, reaching nearly 1.83 times the level of early 2001.

Characteristics of Water Conservation Function
From 2016 to 2020, the mean annual water conservation capacity of Zoige county was approximately 21.2 × 10 4 m 3 /km 2 , and high-value areas were mainly distributed in the southern herbaceous swamp and lake area and eastern mountainous forest area ( Figure  8). The mean annual water conservation capacity in the Nature Reserve from 2016 to 2020 was 19 × 10 4 m 3 /km 2 . Viewed by stages, the water body area in Zoige county was basically stable from 2001 to 2015, and continued to increase after 2015. From 2015 to 2020, the county's water body area rapidly increased, reaching nearly 1.83 times the level of early 2001.

Characteristics of Water Conservation Function
From 2016 to 2020, the mean annual water conservation capacity of Zoige county was approximately 21.2 × 10 4 m 3 /km 2 , and high-value areas were mainly distributed in the southern herbaceous swamp and lake area and eastern mountainous forest area (Figure 8). The mean annual water conservation capacity in the Nature Reserve from 2016 to 2020 was 19 × 10 4 m 3 /km 2 .

Characteristics of Water Conservation Function
From 2016 to 2020, the mean annual water conservation capacity of Zoige county was approximately 21.2 × 10 4 m 3 /km 2 , and high-value areas were mainly distributed in the southern herbaceous swamp and lake area and eastern mountainous forest area ( Figure  8). The mean annual water conservation capacity in the Nature Reserve from 2016 to 2020 was 19 × 10 4 m 3 /km 2 .  fastest growth rates were mainly located in the northern part of the reserve and the forested areas in the eastern part of Zoige county. Overall, water conservation in 15.02% of the county decreased, and the areas with the fastest decrease were mainly distributed in the west and south of Zoige county (Figure 9). Between 2001 and 2020, the annual average water conservation of the nature reserve was 16.89 × 10 4 m 3 /km 2 ; water conservation of the reserve increased at a rate of approximately 0.2 × 10 4 m 3 /km 2 /year, with 80.9% of the area experiencing an increasing trend. The reserve consisted of mainly grassland ecosystems. Although the average water conservation in the entire reserve was slightly lower than that in the whole county, the average water conservation in the grasslands of the reserve was about 1 × 10 4 m 3 /km 2 higher than that of the surrounding grasslands ( Figure 10).  From 2001 to 2020, water conservation in Zoige county increased at a rate of 0.25 × 10 4 m 3 /km 2 /year, and 84.55% of the region exhibited an increasing trend. Areas with the fastest growth rates were mainly located in the northern part of the reserve and the forested areas in the eastern part of Zoige county. Overall, water conservation in 15.02% of the county decreased, and the areas with the fastest decrease were mainly distributed in the west and south of Zoige county (Figure 9).
Between 2001 and 2020, the annual average water conservation of the nature reserve was 16.89 × 10 4 m 3 /km 2 ; water conservation of the reserve increased at a rate of approximately 0.2 × 10 4 m 3 /km 2 /year, with 80.9% of the area experiencing an increasing trend. The reserve consisted of mainly grassland ecosystems. Although the average water conservation in the entire reserve was slightly lower than that in the whole county, the average water conservation in the grasslands of the reserve was about 1 × 10 4 m 3 /km 2 higher than that of the surrounding grasslands ( Figure 10).

Discussion
The purpose of this study was to make full use of satellite remote-sensing images with different resolutions to try to capture the changes in the ecosystem in Zoige from different scales, in order to corroborate changes in statistical analysis indicators and determine the intuitive characteristics of surface information. Furthermore, the objective of this study was also to analyze and evaluate the causes of ecosystem changes at the regional scale, and help us put forward specific and quantitative measures for regional ecological protection and restoration. Remote sensing has proven to be both feasible and effective. However, the existing indicators still require improvement.
The FVC describes key parameters of land surface vegetation, and reflects the structural characteristics of vegetation [30]. For grassland ecosystems, FVC is the basic index that reflects changes in ecological status. Here, we used FVC to monitor and evaluate the basic vegetation status of Zoige county. Furthermore, it was utilized to describe the regional ecological protection status by comparing high-resolution images of grassland and sand variations. However, FVC has obvious limitations in identifying the succession direction of the vegetation community. For example, FVC cannot directly reflect the characteristics of forage productivity, nor can it reflect the quantitative characteristics of dominant forage. Therefore, for future monitoring and evaluation of grassland ecosystem restoration, as well as for desertification control effectiveness analysis, it will be necessary to continuously carry out long-term tracking of ground communities in the Zoige area, especially in the grassland ecosystem restoration area. Additionally, it is necessary to comprehensively evaluate the effectiveness of regional ecological protection in terms of quantity, productivity, and grassland quality.
Due to the influence of climate warming and humidification, glacier melting and surface runoff have increased to a certain extent, which has contributed to the increase in the area of regional lakes and other water bodies, and consequently to the amount of water conservation [35,41]. We used Landsat satellite imagery to extract the water body area with a spatial resolution of 30 m. However, water bodies of subpixel size are widely distributed in the peat swamp area, which is not only an important indicator of regional water conservation function, but also very important for climate change research, including regional greenhouse gas change monitoring [41]. Therefore, to carry out a more refined regional ecological function assessment of the Zoige wetland, it will be necessary to perform higher-resolution identification of water information. In the future, we will explore the driving effects of human activities and climate change on ecological conditions. Water conservation is the dominant ecosystem function in the Zoige area. As an internationally important wetland, the Zoige wetland is also a habitat for rare species such as black-necked cranes, and the function of maintaining biodiversity is also worthy of attention. Therefore, in the future, we will also evaluate ecosystem protection effectiveness and develop technical methods for researching regional biodiversity based on satellite remote-sensing images [41].

Conclusions
From 2001 to 2020, the changes in vegetation and water bodies showed that the ecological status of Zoige improved, especially within the last five years (2015-2020). The vegetation coverage of the county showed an increasing trend; about 95.23% of the county experienced increasing vegetation coverage. During 2015-2020, the growth rate of vegetation coverage was 1.64 times higher than that of the previous 15 years (2001-2015). In the last five years, the water area increased by 77.35%, which is the fastest period of increase in the past two decades.
Local land desertification control, wetland restoration, restrictions on illegal construction in nature reserves, facility demolition, and other protective measures not only reflect the local government's determination and willingness to protect the ecological environment, but also to improve the quality of the environment. Based on an interannual comparison of high-resolution satellite images, some land desertification control areas in the county have been restored to green. For example, the vegetation coverage of typical patches for desertification control had increased by about 5-10% since the implementation of the project in 2016. Some dried lakes and marshes have also been significantly restored. For example, multiphase remote-sensing images showed that the water body of Huahu Lake in the reserve has increased, and the surrounding land has become more swampy.
Nature reserves generally represent areas with good ecological background and fragile ecology. Driven by regional ecological protection policies, the protection effects of nature reserves will be more prominent and direct (for example, faster vegetation restoration, priority of water body restoration, greater force of human interference reduction measures, etc.) In addition, nature reserves can play a leading and exemplary role in the protection of the overall ecological function of the region. Taking the water conservation function of grassland ecosystem as an example, the amount of water conservation within the nature reserve was about 1 × 10 4 m 3 /km 2 higher than that of the surrounding grasslands.
In the future, with the continuous enrichment in satellite remote-sensing data sources, especially with the increase in high-resolution images, it will be possible to extract highspatial resolution and long time series of vegetation, water and other information, which will be very useful to improving the accuracy of regional ecological protection and assessment, and even replace most of the ground investigation and engineering implementation assessment tasks. Information pertaining to regional ecological function changes will be further refined with this new technology.