Next Article in Journal
Recent Accelerating Glacier Mass Loss of the Geladandong Mountain, Inner Tibetan Plateau, Estimated from ZiYuan-3 and TanDEM-X Measurements
Next Article in Special Issue
Comparison of Sentinel-2 and High-Resolution Imagery for Mapping Land Abandonment in Fragmented Areas
Previous Article in Journal
Rapid Identification and Prediction of Cadmium-Lead Cross-Stress of Different Stress Levels in Rice Canopy Based on Visible and Near-Infrared Spectroscopy
Previous Article in Special Issue
The Evolution Response of Ecosystem Cultural Services under Different Scenarios Based on System Dynamics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 12
Issue 3
10.3390/rs12030471
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Do Protected Areas Improve Ecosystem Services? A Case Study of Hoh Xil Nature Reserve in Qinghai-Tibetan Plateau

by
Jie Zeng
1,2,*,
Tianyang Chen
3,
Xiaowei Yao
4 and
Wanxu Chen
4
1
Department of Geography, School of Geography and Information Engineering, China University of Geosciences, Wuhan 430078, China
2
Research Center of Spatial Planning and Human-Environmental System Simulation, China University of Geosciences, Wuhan 430078, China
3
Department of Geography and Earth Science, The University of North Carolina at Charlotte, Charlotte, NC 28223, USA
4
Department of Land Resource Management, School of Public Administration, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(3), 471; https://doi.org/10.3390/rs12030471
Submission received: 26 December 2019 / Revised: 28 January 2020 / Accepted: 30 January 2020 / Published: 2 February 2020
(This article belongs to the Special Issue Remote Sensing Based Monitoring of Terrestrial Ecosystem Service Bundles, Trade-Offs and Synergies)
Browse Figures
Versions Notes

Abstract

:
Although there is a consensus that protected areas (PAs) can provide various ecosystem services, it is unclear whether protected areas effectively contribute to the preservation and enhancement of ecosystem services. We conducted a case study of the Hoh Xil Nature Reserve (HXNR) in Qinghai-Tibetan Plateau, China, in order to examine the effectiveness of PA in the conservation of ecosystem services. First, the dynamics of land use/land cover (LULC) were analyzed based on remotely sensed data sets. Then, the ecosystem service value (ESV) in the PA and non-PA were evaluated using a modified benefit transfer method that had been adjusted using normalized difference vegetation index (NDVI). Finally, hotspot analysis was implemented to reveal the ESV changes for the different districts of the PA. The results of the comparison experiment indicate that: (1) The ESV of the HXNR has considerably increased after it was designated as protected, which had been in decline in the previous stage. The ESVs in a near-by non-PA showed opposite results where the values initially increased but then dropped due to urban expansion and desertification. (2) The areas in HXNR with increased ESV significantly outnumbered the areas that had declining values from 1980 to 2018. For the non-PA, the areas that had increased ESV in 1980–1995 saw a decline in value in 1995–2008; moreover, new areas with decreasing ESV emerged in 2008–2018. (3) The HXNR was found to be more effective than non-PA in improving ecosystem services. (4) The core zone of the nature reserve demonstrated better effectiveness in ecosystem service preservation.

Graphical Abstract

1. Introduction

Protected areas (PAs; e.g., national parks, nature reserves, wilderness areas, landscape protected areas) are geographical spaces with unique natural and cultural resources, where specific regulations are implemented to reduce anthropogenic pressure and to achieve long-term environmental conservation alongside their associated ecosystem services and cultural values [1,2]. Since the tenth meeting of the Conference of the Parties to the Convention on Biological Diversity (CBD) in 2010, the goals of PAs have been extended from mainly biodiversity purposes into preserving ecosystem services (including maintaining biodiversity) [3,4,5]. This development incorporates ecosystem services into veritable policies of global and regional protected areas and biodiversity conservation [2,4,6,7], emphasizing PAs’ functionality in mitigating the threats of increasing human activities [8,9] and enhancing ecosystem services that contribute to human well-being [10,11,12]. Under the Convention on Biological Diversity’s Aichi Biodiversity Targets (Target 11), a new goal was established—that by 2020, 17% of terrestrial and inland waters and 10% of the coastal and marine areas would be integrated into the global protected area network [13]. As of 2018, the number of PAs around the world reached 238,563 with the total area over 2.6 × 107 km2 in 2018, covering 14.9% of the earth’s land surface and 7.3% of the world’s oceans area [14]. However, whether PAs are successful in effectively preserving ecosystem services remains unclear.
There is considerable evidence that PAs provide vital ecosystem services to humans [4,10,11,12], such that assessing the status and dynamics of ecosystem services provided by PAs could be used to indicate their effectiveness. Ecosystem services are defined as the ecological characteristics, functions, or processes derived from functioning ecosystems that directly or indirectly contribute to human wellbeing, which could be categorized into supporting, provisioning, regulating, and cultural services [10,15,16]. Ecosystem services assessment has increased substantially in recent decades [17,18], being regarded as a better approach to understand ecological processes and patterns in PAs [19,20], either in terms of biophysical metrics (e.g., tons of food produced, tons of water yield, or tons of CO2 absorbed) [21,22,23,24] or monetary units [15,25,26,27,28]. Some research and reviews have demonstrated the use of ecosystem services assessment in the confirmation and design of priority PAs [9,29], trade-offs and conflicts identification [30,31,32], and in conservation policy and decision-making [33,34,35]. Different types of ecosystem services provided by PAs have been widely evaluated, not only at national and regional scales [36,37,38] but also at the local scale [5,39,40,41]. However, such studies were mostly constrained either in the number and range of ecosystem services or in the static nature of evaluation, due to infrequent monitoring and limited data sets [5,38].
The most widely established PAs in China are nature reserves. By the end of 2016, 2750 natural reserves have been established in China, covering a total area of about 1.47 × 106 km2, spanning more than 15% of China’s land surface. Among them, 448 are national nature reserves encompassing approximately 10% of China’s land surface area (~9.7 × 105 km2), while the remaining are local reserves [42]. While a significant number of nature reserves have been established and a series of observations and assessments have been conducted on the effectiveness of biodiversity conservation [43,44], few studies have focused on the effectiveness of protecting ecosystem services, especially for Western China. One of the largest, most abundant, and best-preserved nature reserves in China is the Hoh Xil (also called Kekexili) Nature Reserve (HXNR). Located in the hinterland of the Tibetan Plateau and blessed with strategic geographical conditions and unique biodiversity [45,46], the HXNR was established in 1995 and became a national nature reserve in 1997. In 2017, it was officially recognized as a World Natural Heritage Site. Related works have examined the reserve’s habitat suitability [47,48], biodiversity [46,49], and environmental transformation [50,51], which have implicitly demonstrated the effectiveness of HXNR’s conservation policies. However, the dynamics of ecosystem services have not yet been studied in the HXNR.
In order to explore the effectiveness of ecosystem services protection in the HXNR, a systematic evaluation of the dynamic process and spatial distribution of ecosystem services would be vital to local preservation endeavors and would contribute to the global efforts towards biodiversity conservation. Given the limited monitoring datasets in the Hoh Xil area, the use of benefit transfer and expert-driven approaches has been suggested [27,28,52] to assess ecosystem services. The benefit transfer method is used to determine the regional ecosystem service value (ESV) by estimating the economic benefits (or equivalent factors) passed on from one place to another, which is commonly based on the land use/land cover (LULC) data sets. Benefiting from the use of remotely sensed imagery, this approach is appropriate for quick ESV assessment and essential for supporting national and regional decision-making [15,25,26], especially in regions with limited data sets [53,54]. Previous studies have localized equivalent factors, as proposed by Costanza et al. [15] and Xie et al. [28], based mainly on biomass and vegetation cover index [52] (e.g., the net primary productivity (NPP) of natural vegetation [55,56]) while overlooking spatial heterogeneity. However, NPP is not always positively correlated with ecosystem services [57]. In consideration of spatial heterogeneity of ecosystem services, we propose the use of the normalized difference vegetation index (NDVI) to adjust the equivalent factors, which would be used to estimate the ESV in the HXNR. In order to test the effectiveness of PAs spatially, hotspot analysis on ESV changes is employed, which is vital in integrating ESV into the implementation of conservation policies [29,58].
This study assesses the dynamics of LULC in the interior and exterior of the HXNR, before and after its formal establishment (periods of 1980–1995, 1995–2008, and 2008–2018) using land-use change indicators. The changes in the different types of ESV have been evaluated using an NDVI-based modified benefit transfer method. Hotspot analysis of ESV change was employed to examine the spatio-temporal characteristics of ESV. Finally, the ESV dynamics inside and outside of the HXNR were compared to determine the effectiveness of PA towards ecosystem service conservation. This study aims to address the following questions: (1) What were the spatio-temporal features of ESV in the HXNR due to LULC change from 1980 to 2018?; (2) Did the ESV in PA improve significantly compared to non-PA?; and, (3) Which district of PA achieved more effectiveness of ecosystem services preservation?

2. Study Area

The Hoh Xil Nature Reserve (34°19’–36°16’N and 89°25’–94°05’E) is located in the hinterland of the Qinghai-Tibetan Plateau, with an official size of 4.5 × 104 km2 and includes part of the Zhiduo country in Qinghai Province, China (Figure 1). It extends to the Kunlun Mountains at the south and the Ulan Ula Mountain at the north and runs from the Qinghai-Tibet Railway in the east to Tibet Autonomous Region in the west. The average altitude of the HXNR is above 4600 m and is geographically composed of gentle and low undulating lake planes, plains, and terraces [59]. The climate in the region is a typical alpine climate characterized by low temperatures and low precipitation [60]. The annual average temperature of the Hoh Xil Reserve is −10.0 to 4.1°C, and the lowest recorded temperature is −46.4°C; the average annual precipitation is 173.0 to 494.9 mm, which is mainly concentrated in summer. These conditions make living in Hoh Xil region extremely difficult, and given the region’s inhospitable climate, the HXNR has preserved its original landscape and rare special species [45,46], providing unique and vital ecosystem services and contributing to ecosystem conservation and human well-being at the regional and global scales. However, while the region has been designated as a protected area since 1995, maintaining and improving the ecological environment remains to be a major challenge due to the changing climate [60] and pressures from human activities (e.g., grazing and infrastructure construction) [61].
For comparison, we selected Zhiduo county as control, which is a region not designated as a protected area and is located in the southeast of the HXNR with an area of approximately 1.08 × 104 km2. The area between the control area and HXNR is part of the Sanjiangyuan Nature Reserve, and will not be discussed in this study due to the limitation of data.

3. Materials and Methods

3.1. Data Sources and Processing

The LULC data for the period 1980–2018 were retrieved from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn) [62,63]. The LULC data initially had seven categories and 25 sub-categories (see Table 1). The LULC data for 1980, 1995, and 2008 were generated using Landsat Thematic Mapper™/Enhanced Thematic Mapper (ETM) remotely sensed imagery, while the LULC data for 2018 were updated based on the 2015 LULC data from Landsat 8 remotely sensed imagery and generated by manual interpretation. Using field surveys and records for comparison, the overall identified accuracy of the LULC data sets for the seven land use/cover types reached 94.3%, while for the accuracy for the 25 subclasses types reached 91.2% [62,63]. The spatial resolution was 30 × 30 m, and the land use/land cover were reclustered into six primary types: forest, grassland, water bodies, wetland, urban, and desert.
The NDVI data set, which was used to revise the equivalent factors in this study, was downloaded from the Geospatial Data Cloud site of the Computer Network Information Center at the Chinese Academy of Sciences (http://www.gscloud.cn) [64].

3.2. LULC Dynamic Analysis

The spatiotemporal characteristics of the LULC in HXNR and the control area were analyzed using land-use dynamic degree index (LUDD), comprising the single LUDD (SLUDD) and comprehensive LUDD (CLUDD) [65]. SLUDD calculates the change rates of certain single land-use types, while CLUDD determines the overall change rate of LULC. SLUDD and CLUDD can be calculated using Equations (1) and (2), respectively:
S L U D D = A i , t 2 − A i , t 1 A i , t 1 × 1 T × 100 %
C L U D D = [ ∑ i = 1 n Δ A i − j ∑ i = 1 n A i , t 1 ] × 1 T × 100 %
where Ai,t1, and Ai,t2 represent the areas of land use type i at the start (t1) and end (t2) of the study period, respectively; ΔAi−j is the area of land use type i converted to land use type j (j = 1,2,…n, i ≠ j) during the study period; n is the number of land-use types in the study area; and, T is the study period.

3.3. Calculation of ESVs

In this study, we adopted the value transfer method developed by Xie et al. [28] to estimate the ESVs inside and outside of the HXNR. The benefit transfer model proposed by Costanza et al. [15] and Xie et al. [28] has been commonly used to evaluate ESVs using monetary values designated to the different LULC types. Considering the inapplicability of the ecosystem services classification and equivalent coefficients proposed by Costanza et al. [15] for China, Xie et al. [28] proposed modifications to the ecosystem services classification and the equivalent factors table based on a survey of more than 700 ecological experts. In the modified version, the ESV for food production of croplands was set to 1, while the equivalent ESV coefficients for the other ecosystems were based on their relative importance to croplands. The value of each equivalent ESV factor was identified as 1/7 of the market value of the annual regional average grain yield, equal to 58.5 USD/ha in China. The equation to calculate ESV is given as follows:
E S V = ∑ k = 1 m ∑ i = 1 n ( V C i , k × A i )
where VCi,k is the kth type of ESV coefficient for land use type i; Ai represents the ith type of LULC; m and n indicate the number of sub-categories of ecosystem services and the LULC types, respectively.
However, the value transfer method proposed by Xie et al. [28] overlooked the spatial heterogeneity [52,57]. Previous studies have enhanced Xie’s model by considering biomass and vegetation cover index; however, they still have deviations on ESV evaluation [55,56]. In contrast, the NDVI index is shown to have a strong correlation with meteorological condition and local resource background, especially in areas with a large number of different vegetation covers [66], which is believed to reflect the spatial patterns of ecosystem background and associated with ecosystem services. Thus, we use NDVI to adjust the ESV coefficients to localize the equivalent factors in the Hoh Xil region compensating for the limitations of the previous model from neglecting spatial heterogeneity. The corrected ESVs were evaluated by Equations (4) and (5).
V f = ( N D V I f − N D V I m i n N D V I m a x − N D V I m i n ) × 100 %
E S V c o r r e c t e d = E S V × V f V
where NDVIf, is the average annual normalized difference vegetation index of grid f; NDVImax and NDVImin represent the maximum and minimum of the average annual NDVI in China; Vf is the normalized NDVIf; and, V is the average of normalized NDVI in China.

3.4. Hotspot Analysis of ESV Change

Hotspot analysis is a widely used to reveal the spatial locations of statistically significant high-value clusters (hotspots) and low-value clusters (cold spots). Theoretically, a hotspot indicates that its surrounding areas are statistically in high values, while a cold spot is surrounded by low value spots. Therefore, the effects of abnormally high and low values can be excluded [57,58]. In this study, we apply hotspot analysis to identify the distribution of hotspots and cold spots showing ESV changes in the HXNR. As the hotspots and the cold spots of ESV change illustrate the significant points of ESV increase and decrease, it would reveal the effectiveness of PA spatially. Based on the calculation of adjusted ESV for each LULC patch in different years, the change in ESV was assigned to each patch. To enhance feasibility and flexibly, the patches with modified ESV were selected and vectorized and were then applied to the Hotspot Analysis (using the Getis-Ord Gi* statistic) to reveal the statistically significant ESV change areas. The z-scores and p-values within the analysis indicate standard deviations and probabilities. The typical probabilities are 0.01, 0.05, and 0.1, and the z-scores at < −1.65 or > 1.65, < −1.96 or > 1.96, and < −2.58 or > 2.58 are critical for 90%, 95%, and 99% confidence levels. A positive and significant z-score indicates that the change in ESV in the vicinity is relatively in high value (hotspots of the ESV changed points). In contrast, a z-score that is negative and significant suggests that the change in ESV in the surrounding area is relatively in low value (cold spots of the ESV changed points).

4. Results

4.1. Land Use/Land Cover Dynamics

4.1.1. LULC Change from 1980 to 2018

For the HXNR, due to its high altitude and harsh climate, the region is considered as not suitable for human habitation. As a result, there were no croplands or construction lands, as shown in Table 2. The main LULC types from 1980 to 2018 were desert and grassland (Figure 2). Desert accounted for 55.98%, 55.58%, 56.03%, and 55.08% in 1980, 1995, 2008, and 2018, respectively; grassland constituted 35.31%, 36.23%, 35.27%, and 35.91% for the same periods. The average proportions of waters and wetlands in the HXNR were about 7.1% and 1.4%, while dispersed patches of forest had a total area of approximately 7.20 km2.
For the control area, there were no cultivated lands similar to the HXNR, as shown in Figure 2. The area of grassland was the largest among the six LULC types, comprising around 90% of the total extent from 1980 to 2018 (Table 2). Desert in the control area was significantly smaller than in the HXNR but still was the second-largest LULC type, covering approximately 10% of the control area. And despite the control area being the capital of Zhiduo county where the population and industries are located, the urban area still accounts for less than 0.1%.
During the period 1980–2018, the areas of grassland, wetland, and waterbodies in the HXNR increased significantly by 293.36 km2, 110.43 km2, and 37.42 km2, respectively, with the SLUDD of 0.05%, 0.03% and 0.44% (see Table 3). Meanwhile, the SLUDD of the desert in the HXNR is −0.05%, with a total reduction of 441.22 km2. In contrast, the grassland in the control area had been reduced considerably, which diminished by 99.35 km2 (the SLUDD was −0.03%) from 1980 to 2018. The urban and desert areas showed relatively significant expansion, increasing by 5.36 km2 and 94.39 km2 and having the SLUDD of 0.24% and 19.08%, respectively.

4.1.2. LULC Change before and after 1995

For the period 1980 to 1995 (Figure 3a), the CLUDD of the HXNR was estimated to be 6.05%, compared to 8.43% in the control area. During this period, the conversion of desert to green spaces (i.e., forests, grasslands, and wetlands) was substantial, with a total area of 1508.27 km2 (Table 4). However, newly desertified lands had an aggregated area of 1312.6 km2, mainly due to grassland degradation. Similarly, the area of desert greening in the control area was also pronounced. The desert area converted into grasslands, waterbodies, and wetlands reached 778.86 km2, while the area of newly desertified lands was relatively small, with a total area of 70.89 km2. Another noticeable change in the control area was the conversion of waterbodies and wetlands into grasslands with converted areas of 79.12 km2 and 44.36 km2, respectively. At this stage, urban expansion was slow, expanding only by 0.96 km2.
From 1995 to 2008, large areas of land in the HXNR were converted from and into deserts (Figure 3b). While 1326.43 km2 of desert lands were reclaimed through desert greening, this area was exceeded by desertified land, which reached 1548.07 km2 (Table 4). Previous conversions of waterbodies and wetlands into grasslands had been reversed to a certain extent at this period. The area of HXNR grasslands converted into waterbodies and wetlands reached 84.83 km2. As for the control area, desertified lands (781.54 km2) were substantially more extensive than those that have been reclaimed through desert regreening (63.32 km2). Similar to the HXNR, the grasslands converted into waterbodies and wetlands (44.98 km2) were also significantly larger than the waterbodies and wetlands being converted into grasslands (5.11 km2). The CLUDD of the HXNR and the control area were 6.14% and 8.39%, respectively.
From 2008 to 2018 (Figure 3c), protection policies and actions had been greatly improved, and human activities have been effectively monitored and regulated in the HXNR (e.g., poaching activities in the HXNR have been eliminated since 2008). The CLUDD in the HXNR was 3.55%, exceeding that in the control area (2.34%). The reclaimed area from desert greening (i.e., transformed into forests, grasslands, wetlands, and waterbodies) totaled 839.86 km2, which is much larger than the desertified regions (368.16 km2) (Table 4). In the same period, the area of wetlands also increased substantially. The area converted between wetlands and waterbodies reached 400.55 km2, which includes water-to-wetland (271.81 km2) and wetland-to-water (128.72 km2). In contrast, the control area still reflected a distinct trend of desertification. The desert expanded by 159.84 km2, while desert greening was only able to reclaim 73.21 km2. Urban expansion was significant at this stage, with the SLUDD of 66.51% (Table 1) and an expansion area of 5.51 km2, which mainly encroached on grasslands, wetlands, and deserts.

4.2. Ecosystem Services Value Change

4.2.1. ESV Change in Different Periods

From 1980 to 2018, the total ESV of the HXNR showed an overall growth trend, rising from 764.11 million USD in 1980 to 777.92 million USD in 2018, which is an increase of 13.81 million USD (Table 5). In contrast, the total ESV of the control area declined, with a total reduction of 5.33 million USD and an average annual decrease of 0.14% during the study period. The values of HXNR’s provisioning, regulating, supporting, and cultural services increased by 0.59, 9.54, 2.48, and 1.19 million USD from 1980 to 2018, while the values for the control area decreased by 0.38, 2.77, 1.84, and 0.33 million USD respectively.
Even though the ESVs in the HXNR have shown upward trends, their changes varied in the different stages (Table 5). From 1980 to 1995, the total ESV decreased annually by 0.01%, while the values of provisioning and supporting services increased by 0.15% and 0.13%, respectively. During 1995–2008, the total ESV also declined, but the yearly rate of decline was reduced to 0.004. The values of the provisioning and supporting services decreased by 0.18% and 0.16%, while the values of regulating and cultural services increased by 0.07% and 0.05%. In the period of 2008–2018, the ESVs of the HXNR generally expanded, with the total ESV increasing annually by 0.20%. The values of provisioning, regulating, supporting, and cultural services grew by 0.17%, 0.24%, 0.13%, and 0.19%.In contrast to the HXNR, the ESVs in the four ecosystem categories simultaneously increased or decreased for each period. From 1980 to 1995, the total ESV rose by 0.19%, with the provisioning, regulating, supporting, and cultural ecosystem services increasing in value by 0.34%, 0.11%, 0.30%, and 0.10%, respectively (Table 5). For 1995–2008, the total ESV showed an annual decline of 0.22%, with all the ecosystem service types decreasing in value. For 2008–2108, the downtrend trend in the total ESV continued having a yearly decline of 0.07%, and all ecosystem service types registering negative rates of change.

4.2.2. Spatial Change of ESV Per Unit

In both the HXNR and in the control area, the total ESV per unit have been spatially high in the south while low in the north (see Figure 4), which is similar to the spatial configuration of the average NDVI. The highest value per unit in the HXNR reached 2785 USD/ha between 1980 and 2018, while the control area had the highest value at 4797 USD/ha. High-value zones of the HXNR were mainly distributed in the southeast of the reserve, which were mostly situated around Trache Lake, Rigachi Mountain, and Dongduoqu River, while the rest were located along the route of Lake Trache-Dorge Co-Lake Hoh Xil. Due to high NDVI values, the adjusted ESV per unit in the control area was approximately four times higher, on average, than in the HXNR (Table 6). There were still low-value zones in the control area, mainly located along the 313 and 408 Provincial Highway. Table 6 shows the change in average ESV per unit from 1980 to 2018. Although the average value per unit slightly dropped for two periods (1980–1995 and 1995–2008), the average total ESV per unit had increased from 155.74 USD/ha in 1980 to 158.59 USD/ha in 2018 for the HXNR. In the control area, the average total ESV per unit decreasing from 653.97 USD/ha in 1980 to 648.98 USD/ha in 2018.
Spatially, the changes in the total ESV per unit at the HXNR for 1980–1995 (Figure 5a) and 1995–2008 (Figure 5b) were mainly situated around the Kunlun Mountain Pass, south of the Rigachi Mountain, and the line from Dorge Co Lake to Chuma’er River, primarily due to the conversion of deserts and grasslands. In some of these places, value losses occurred in 1980–1995 but were regained in 1995–2008 (e.g., the areas around the Kunlun Mountain Pass and Dorge Co Lake), while other areas had the reverse trend. From 2008 to 2018 (Figure 5c), the areas that increased of the total ESV per unit in the HXNR were significantly more than those that declined, and are mainly concentrated around Xijir Ulan Lake, north of Lake Hoh Xil and north of Daxue Peak. The ESV increases in these areas can be attributed to the progress of desert greening and the transformation of lakes and wetlands. In contrast, areas with increased ESV could scarcely be found in the control area, while those with decreased ESV are situated around the towns of Jiajiboluoge and Lixin due primarily to urban growth and desert expansion.

4.3. Hotspots of ESV Change

Using the Hot Spot Analysis (Getis-Ord Gi*) tool in ArcGIS Pro 2.1, the statistically significant high values (hotspots) and low values (cold spots) of ESV change in the HXNR were identified. Table 7 summarizes the statistical results of hotspot analysis at 95% confidence level. From 1980 to 1995, the hotspots (49% of the total ESV changed patches) outnumbered the cold spots (44%), but there were more cold spots (31%) in the core zone than hotspots (9%), which indicates that at this stage, the decline of ESV in the core zone had played a crucial role in reducing the total ESV in HXNR. From 1995 to 2008, while cold spots (50%) exceeded the number of hotspots (44%) in total, there were more hotspots found in the core area (31%) than cold spots (10%). Since the total ESV was still in a downward trend at this period, this means that the increase in ESV in the core zone was not enough to offset losses in the buffer zone. From 2008 to 2018, the proportion of hotspots was significantly higher than that of cold spots, accounting for 52% and 28% of the total changed patches, respectively, prompting a notable increase in HXNR’s aggregated ESV. For the entire period of 1980–2018, the hotspots in the HXNR (51%) were far greater than cold spots(30%). Hotspots (28%) in the core zone outnumbered those in the buffer zone (23%), while cold spots in the core zone (11%) were less than those found in the buffer zone (19%). This suggests that the HXNR‘s core zone had been more effective in improving ecosystem services.
In terms of spatial configuration (Figure 6), for 1980–1995, the hotspots of ESV change were mainly concentrated in the south of Fengxue Beach (located in the southwest corner of the core zone of HXNR), the south of Trache Lake-Rigachi Mountain, and the north of Chuma’er River in the buffer zone. Around Dorge Co Lake, there were also some scattered hot spots. Meanwhile, cold spots were mainly located around Cuodarima Lake-Dorge Co Lake, Xijir Ulan Lake, and Lianhu Lake in the core zone, as well as in the Chuma’er River section and Lubei Beach. From 1995 to 2008, the hotspot and cold spot regions were noticeably in reverse of the preceding period, and both distributions expanded significantly, either in the core area or in the buffer zone. For the succeeding period of 2008–2018, the distribution of hotspots and cold spots was more widespread and uniform. The regions with relatively concentrated hotspots were the south of Trache Lake-Rigachi Mountain and the areas surrounding Hoh Sai Lake, Xijir Ulan Lake, Lexie Wudan Lake, Lianhu Lake, and Moon Lake. The cold spots were more concentrated in the areas around Haidingnuo’er Lake and Yanhu lake and the line of Daxue Peak-Buka Daban Mountain Peak.

5. Discussion

5.1. Characteristics of ESV Change Due to LULC Dynamics

Based on the LULC data sets, we used NDVI to enhance the equivalent ESV coefficients proposed by Xie et al. [28] in order to evaluate the ESVs from 1980 to 2018 in the study area. In the HXNR, although declining slightly in the first two periods (1980–1995 and 1995–2008), the total ESV exhibited a general upward trend, increasing from 764.11 million USD in 1980 to 777.92 million USD in 2018. Previous studies have explored the ESV change of the Sanjiangyuan area (including parts of the HXNR) and found that the ESV decreased initially and then increased considerably [67]. Our findings support this conclusion that there is a downward trend from the 1980s to the mid-1990s.
Spatially, the ESV change in HXNR can be found in the Kunlun Mountain Pass, south of the Rigachi Mountain, and the line from Dorge Co Lake to Chuma’er River. The areas with increased ESV can be found in the north of Hoh Xil Lake, north of Daxue Peak, and the surrounding regions of Xijir Ulan Lake, which could be attributed with the conversion of desert to grasslands in north of Hoh Xil Lake and north of Daxue Peak and in the expansion of lakes and wetlands around the Hoh Xil Lake and Xijir Ulan Lake. Yao et al. [50] have reported the continuous expansion of Hoh Xil Lake and Xijir Ulan Lake from 1970 to 2010 due to the increase in precipitation and temperature, which is consistent with the findings of this study. Moreover, given the almost opposing change trends in ESVs before and after the establishment of PA, the restoration and maintenance of ecosystem services in the HXNR have clearly benefited with the establishment of the PA after 1995.

5.2. Effectiveness of the HXNR on Ecosystem Services Protection

In this study, we chose a non-protected region as a control area, situated in the southeast of the HXNR. By comparing the ESV changes between PA and non-PA, We determined the effectiveness of the HXNR in promoting ecosystem services. Compared with the HXNR, the ESV in non-PA demonstrated a contrary trend change, which increased by 2.85 million USD in 1980–1995, decreased sharply in 1995–2008 (decreased by 2.87 million USD), and continued dropping by 0.66 million USD in 2008–2018, mainly due to the impact of urban expansion and desertification. Without protection, the non-PA suffered more ESV losses, even with an initial upward trend. Additionally, the results of hotspot analysis showed that most of the sites with increased ESV belong to the HXNR’s core zone, significantly contributing to the growth of the total ESV. Given the characteristics of ESV change through the different phases in establishing the PA and the contrasting developments between PA and non-PA, the establishment of the HXNR has clearly mitigated the degradation of the ecosystem and has been effective in enhancing ecosystem services. While there has been limited research evaluating the conservation effectiveness in the HXNR, several studies have observed improvements in habitat suitability and biological population [46,47,48], and the effective regulation of human activities in the HXNR after the 1990s [61], which are consistent with the conclusions of this study.
Even though our findings suggest that the HXNR was effective in improving ecosystem services, it is still less clear whether all PAs can successfully contribute to the enhancement of ecosystem services. Palomo et al. [68] demonstrated that high levels of PAs, such as the National Parks and Natural Reserves, did not ensure the provision of multiple ecosystem services due to the disconnect between protected areas and society. Castro et al. [36] suggested that the current PA networks supply only slightly higher levels of regulating services than non-PAs at the landscape level. On the other hand, based on nine individual cases, Eastwood et al. [5] assert that PAs deliver higher levels of ecosystem services than non-PAs. Given the variations in research scales, methods, and settings, there is still considerable uncertainty whether PAs are effective in promoting ecosystem services. The differences in outcomes would be largely dependent on geographical conditions, targets, and policies of PAs [68]. Thus, to ensure the effectiveness of PAs in promoting ecosystem services, ecosystem assessments and environmental evaluations should integrate more the use ecosystem services, particularly at the local scale.

5.3. Application of ESV on the Effectiveness Assessment of PAs

Based on LULC data from remotely sensed imagery, this study evaluated the ESV using an NDVI-based modified version of the benefit transfer method. The ESV evaluation was integrated into the effectiveness assessment of PAs by comparing the ESV changes at different phases between PA and non-PA. This approach exhibited feasibility and practicality in various aspects that could be implemented in other PAs. First, with the use of remotely sensed imagery, which has wide observation range, fast acquisition speed, and short update period, a quick evaluation of ESV can be implemented on a large-scale, long-term basis, and at low costs [69], particularly in remote and data-scarce areas.
Second, this approach assesses the various types of ecosystem services in monetary units using benefit transfer method, which demonstrates flexibility in analyzing the levels of different ecosystem services and practicability in comparing economic benefits of PAs. Such an approach that analyzes benefits in pecuniary terms can be integrated into evaluation studies of PAs to provide a more effective platform that would be easily understood by decision-makers and other stakeholders [9,70]. For example, based on the ESV change in monetary units, it is workable for the government and environmental protection organizations to evaluate the income on the investment of ecosystem protection projects, and also provides a basis for the improvement of ecological compensation standards for local residents. Since the benefit transfer method proposed by Costanza et al. [15] and Xie et al. [28] overlooked spatial heterogeneity [52,57], our revised approach uses NDVI to modify equivalent coefficients in order to overcome this limitation and lessen deviations inherent in other techniques (e.g., modification based on NPP) [52,55,56]. In vegetation monitoring, the NDVI varied in space with the meteorological condition and local resource background [66], especially in the regions with numerous different vegetation covers [71]. Thus, NDVI is adequate to represent the spatial heterogeneity of ecosystems and to adjust the benefit transfer method.
Additionally, the ESV evaluation based on grid cells and the hotspot analysis of ESV change provide more spatial details of ESV and allows the spatial evaluation of ecosystem service improvements in PAs. This offers additional value in the decision-making of PAs’ protection and governance [53] since strategies and policies promoting ecosystem services need to be adjusted in accordance with the varying spatial conditions.

5.4. Limitations and Future Directions

This study reclassified land use/land cover into six primary categories based on remotely sensed LULC data and evaluated the ESV of the HXNR by using an NDVI-revised benefit transfer method. This approach is preferred when quick ESV evaluation is required but the available data are limited. However, one major limitation in this approach is that the differences in areas having the same LULC type were largely ignored, which may influence the accuracy of the evaluation results. In fact, the LULC primary type can be further subdivided into multiple subclasses (e.g., dry lands and paddy fields in croplands, broad-leaved woodlands, and shrublands in the forest), such that the ecosystem services provided by each LULC subclass could significantly vary [16,55]. Considering the diversity of bio-species and landscape in PAs, more detailed research on equivalent coefficients in benefit transfer models is essential, particularly those focusing on variations in LULC subclasses. Moreover, this method is also limited by ground truthing as the study focused on a long-term dynamic process covering a large area. Hence, incorporating reliable field verification into the ESV assessment needs to be fully considered in future studies.
Considering the strong correlation with meteorological condition and local resource background [66], the NDVI index was used to modify the benefit transfer method, which make up for the shortage of ignoring the spatial heterogeneity of ESV. However, the uncertainty of the NDVI index was overlooked. For instance, the changeable relationship between NDVI and temperature in different regions [71] may affect the accuracy of ESV evaluation. Accordingly, the NDVI-based localized benefit transfer method needs to integrate multiple-indicates in future ESV evaluation studies. Previous studies developed various satellite-based vegetation indices applied in monitoring the vegetation [69,71]. Such indices (e.g., the Vegetation Condition Index (VCI), the Vegetation Health Index (VHI), the Perpendicular Vegetation Index (PVI), and the two-band version of the enhanced Vegetation Index (EVI2) et al.) representing a range of spectral responses of vegetation conditions [69], would be appropriate to adjust the value transfer method. However, each vegetation index varies in performance due to environmental changes (e.g., weather) [71]. Therefore, there is still a call for studies on the application scenarios of different vegetation indices in modifying value transfer method.
In addition, by comparing the ESV change in PA and non-PA at different periods, our findings demonstrated the effectiveness of PA in improving the different types of ecosystem services. Although the spatial effectiveness was evaluated using hotspots analysis, it remains incomplete as spatial conflicts and trade-offs within different stakeholders on ecosystem services were not considered. Trade-offs of ecosystem services, arising from human preferences for specific sets of ecosystem services (e.g., human preferences often prioritize provisioning services over regulating services), have shown significant influence on policies and strategies in environmental management [70,72]. Thus, future research focusing the effectiveness of PAs would need to consider ecosystem service trade-offs, in order to provide a comprehensive understanding of ecosystem services and to ensure that the target of ecosystem conversation is not dominated by short-term needs [70], which could be more valuable in policy-making and governance.

6. Conclusions

We analyzed the dynamics of LULC in the HXNR and the control area using remote sensing techniques and evaluated the ESV change due to LULC using the NDVI-based modified value transfer method. Comparing ESV changes between PA and non-PA at different periods, we confirmed the effectiveness of the HXNR in improving ecosystem services and identified the spatial effectiveness at the different districts. The main findings of this study are as follows:
(1) The ESV of the HXNR initially decreased in value and then increased considerably over time, with the total ESV of 764.11 million USD in 1980 and increasing to 777.92 million USD in 2018. The four categories of ecosystem services (provisioning, regulating, supporting, and cultural) all increased in value, largely due to the conversions of desert-to-green and waterbodies-to-wetland. In contrast, the ESVs of the control area improved in 1980–1995 and declined in the succeeding years, with the total ESV decreasing by 5.33 million USD mainly due to urban expansion and desertification.
(2) The change of ESV per unit in the HXNR were concentrated in regions around the Kunlun Mountain Pass, Xijir Ulan Lake, north of Hoh Xil Lake, south of the Rigachi Mountain, north of Daxue Peak, and the line from Dorge Co Lake to Chuma’er River. Generally, the areas that achieved ESV growth significantly outnumbered areas with declining values. As for the control area, the regions that saw an increase in ESV in 1980–1995 experienced a decline in values in 1995–2008; also, the areas around the towns of Jiajiboluoge and Lixin experienced a value decline in 2008–2018.
(3) Comparing the ESV changes between the HXNR and the control area at different periods, our results show the HXNR has been more effective in improving the value of ecosystem services than the control area.
(4) The core zone of the nature reserve achieved a higher level of effectiveness in improving ecosystem services, given that the hotspots of ESV change were more prominent in the core zone, while the cold spots were much concentrated in the buffer zone.
The methodology used in this study can be widely applied for prompt and effective evaluation of PA in enhancing ecosystem services, especially for long-term conservation assessment in areas with limited data. In order to comprehensively understand the effectiveness of PAs on ecosystem services conservation, more research and case studies from around the world have to be undertaken, particularly at the local scale.

Author Contributions

Conceptualization, J.Z. and X.Y.; methodology, T.C. and W.C.; formal analysis, J.Z., T.C., and X.Y.; resources, J.Z. and W.C.; writing—original draft preparation, J.Z., and T.C.; writing—review and editing, J.Z. and X.Y.; visualization, J.Z. and T.C.; supervision, J.Z.; project administration, J.Z.; funding acquisition, J.Z. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 41701629), the MOE(Ministry of Education of China) Project of Humanities and Social Sciences (No. 19YJCZH224), and Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan), grant number CUGW170206.

Acknowledgments

Much thanks to Wei-Ning Xiang and Wenwu Tang (Department of Geography and Earth Science, the University of North Carolina at Charlotte) for their valuable assistance and suggestions to improve this study. We are grateful to the comments provided by reviewers and the editor which were of great value for improving the final version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chape, S.; Spalding, M.; Jenkins, M. The World’s Protected Areas: Status, Values and Prospects in the 21st Century; University of California Press: Berkeley, CA, USA, 2008. [Google Scholar]
  2. Dudley, N. Guidelines for Applying Protected Area Management Categories; IUCN: Gland, Switzerland, 2008. [Google Scholar]
  3. Rands, M.R.; Adams, W.M.; Bennun, L.; Butchart, S.H.; Clements, A.; Coomes, D.; Entwistle, A.; Hodge, I.; Kapos, V.; Scharlemann, J.P. Biodiversity conservation: Challenges beyond 2010. Science 2010, 329, 1298–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wittmer, H.; Gundimeda, H. The Economics of Ecosystems and Biodiversity in Local and Regional Policy and Management; Earthscan: London, UK; Washington, DC, USA, 2012. [Google Scholar]
  5. Eastwood, A.; Brooker, R.; Irvine, R.J.; Artz, R.R.E.; Norton, L.R.; Bullock, J.M.; Ross, L.; Fielding, D.; Ramsay, S.; Roberts, J.; et al. Does nature conservation enhance ecosystem services delivery? Ecosyst. Serv. 2016, 17, 152–162. [Google Scholar] [CrossRef] [Green Version]
  6. European Commission. Our Life Insurance, Our Natural Capital: An EU Biodiversity Strategy to 2020; European Commission: Brussels, Belgium, 2011. [Google Scholar]
  7. Claret, C.; Metzger, M.J.; Kettunen, M.; ten Brink, P. Understanding the integration of ecosystem services and natural capital in Scottish policy. Environ. Sci. Policy 2018, 88, 32–38. [Google Scholar] [CrossRef] [Green Version]
  8. Jones, K.R.; Venter, O.; Fuller, R.A.; Allan, J.R.; Maxwell, S.L.; Negret, P.J.; Watson, J.E. One-third of global protected land is under intense human pressure. Science 2018, 360, 788–791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Xu, W.; Xiao, Y.; Zhang, J.; Yang, W.; Zhang, L.; Hull, V.; Wang, Z.; Zheng, H.; Liu, J.; Polasky, S.; et al. Strengthening protected areas for biodiversity and ecosystem services in China. Proc. Natl. Acad. Sci. USA 2017, 114, 1601–1606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Millennium Ecosystem Assessment. Ecosystems and Human Well-Being; Island Press: Washington, DC, USA, 2005. [Google Scholar]
  11. Anton, C.; Young, J.; Harrison, P.A.; Musche, M.; Bela, G.; Feld, C.K.; Harrington, R.; Haslett, J.R.; Pataki, G.; Rounsevell, M.D.A.; et al. Research needs for incorporating the ecosystem service approach into EU biodiversity conservation policy. Biodivers. Conserv. 2010, 19, 2979–2994. [Google Scholar] [CrossRef]
  12. Jiricka-Purrer, A.; Tadini, V.; Salak, B.; Taczanowska, K.; Tucki, A.; Senes, G. Do Protected Areas Contribute to Health and Well-Being? A Cross-Cultural Comparison. Int. J Environ. Res. Public Health 2019, 16, 1172. [Google Scholar] [CrossRef] [Green Version]
  13. Convention on Biological Diversity. Aichi Biodiversity Targets. Available online: https://www.cbd.int/sp/targets (accessed on 1 November 2019).
  14. UNEP-WCMC; IUCN; NGS. Protected Planet Report 2018; World Conservation Monitoring Centre (UN Environment Conservation Monitoring Centre): Cambridge, UK, 2018. [Google Scholar]
  15. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; Oneill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  16. Costanza, R.; de Groot, R.; Braat, L.; Kubiszewski, I.; Fioramonti, L.; Sutton, P.; Farber, S.; Grasso, M. Twenty years of ecosystem services: How far have we come and how far do we still need to go? Ecosyst. Serv. 2017, 28, 1–16. [Google Scholar] [CrossRef]
  17. Zhao, J.; Xiao, H.; Wu, G. Comparison analysis on physical and value assessment methods for ecosystems services. Chin. J. Appl. Ecol. 2000, 11, 290–292. [Google Scholar]
  18. Häyhä, T.; Franzese, P.P. Ecosystem services assessment: A review under an ecological-economic and systems perspective. Ecol. Modell. 2014, 289, 124–132. [Google Scholar] [CrossRef]
  19. Armsworth, P.R.; Chan, K.M.A.; Daily, G.C.; Ehrlich, P.R.; Kremen, C.; Ricketts, T.H.; Sanjayan, M.A. Ecosystem-service science and the way forward for conservation. Conserv. Biol. 2007, 21, 1383–1384. [Google Scholar] [CrossRef] [PubMed]
  20. Cumming, G.S.; Allen, C.R.; Ban, N.C.; Biggs, D.; Biggs, H.C.; Cumming, D.H.M.; Vos, A.D.; Epstein, G.; Etienne, M.; Maciejewski, K.; et al. Understanding protected area resilience: A multi-scale, social-ecological approach. Ecol. Appl. 2015, 25, 299–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Burkhard, B.; Kroll, F.; Nedkov, S.; Müller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic. 2012, 21, 17–29. [Google Scholar] [CrossRef]
  22. Kandziora, M.; Burkhard, B.; Müller, F. Interactions of ecosystem proper-ties, ecosystem integrity and ecosystem service indicators—A theoretical matrix exercise. Ecol. Indic. 2013, 28, 54–78. [Google Scholar] [CrossRef]
  23. Stürck, J.; Poortinga, A.; Verburg, P.H. Mapping ecosystem services: The supply and demand of flood regulation services in Europe. Ecol. Indic. 2014, 38, 198–211. [Google Scholar] [CrossRef]
  24. Bao, Y.; Li, T.; Liu, H. Spatial and temporal changes of water conservation of Loess Plateau in northern Shaanxi province by InVEST model. Geogr. Res. 2016, 35, 664–676. [Google Scholar]
  25. Costanza, R.; de Groot, R.; Sutton, P.; van der Ploeg, S.; Anderson, S.J.; Kubiszewski, I.; Farber, S.; Turner, R.K. Changes in the global value of ecosystem services. Glob. Environ. Change 2014, 26, 152–158. [Google Scholar] [CrossRef]
  26. de Groot, R.; Brander, L.; van der Ploeg, S.; Costanza, R.; Bernard, F.; Braat, L.; Christie, M.; Crossman, N.; Ghermandi, A.; Hein, L.; et al. Global estimates of the value of ecosystems and their services in monetary units. Ecosyst. Serv. 2012, 1, 50–61. [Google Scholar] [CrossRef]
  27. Dasgupta, P. Nature in economics. Environ. Resour. Econ. 2008, 39, 1–54. [Google Scholar] [CrossRef] [Green Version]
  28. Xie, G.; Zhen, L.; Lu, C.; Yu, X.; Cao, C. Expert Knowledge Based Valuation Method of Ecosystem Services in China. J. Nat. Resour. 2008, 23, 911–919. [Google Scholar]
  29. Naidoo, R.; Balmford, A.; Costanza, R.; Fisher, B.; Green, R.E.; Lehner, B.; Ricketts, T.H. Global mapping of ecosystem services and conservation priorities. Proc. Natl. Acad. Sci. USA 2008, 105, 9495–9500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Kovács, E.; Kelemen, E.; Kalóczkai, Á.; Margóczi, K.; Pataki, G.; Gébert, J.; Málovics, G.; Balázs, B.; Roboz, Á.; Krasznai Kovács, E.; et al. Understanding the links between ecosystem service trade-offs and conflicts in protected areas. Ecosyst. Serv. 2015, 12, 117–127. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, Y.; Tang, M.; Zhang, Z.; Wu, B. Whether Urban Development and Ecological Protection Can Achieve a Win-Win Situation—The Nonlinear Relationship between Urbanization and Ecosystem Service Value in China. Sustainability 2019, 11, 3277. [Google Scholar] [CrossRef] [Green Version]
  32. Karp, D.S.; Mendenhall, C.D.; Callaway, E.; Frishkoff, L.O.; Kareiva, P.M.; Ehrlich, P.R.; Daily, G.C. Confronting and resolving competing values behind conservation objectives. Proc. Natl. Acad. Sci. USA 2015, 112, 11132–11137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Heather Tallis, P.K.; Michelle, M.; Chang, A. An ecosystem services framework to support both practical conservation and economic development. Proc. Natl. Acad. Sci. USA 2008, 105, 9457–9464. [Google Scholar] [CrossRef] [Green Version]
  34. Scolozzi, R.; Schirpke, U.; Morri, E.; D’Amato, D.; Santolini, R. Ecosystem services-based SWOT analysis of protected areas for conservation strategies. J. Environ. Manage. 2014, 146, 543–551. [Google Scholar] [CrossRef]
  35. García-Llorente, M.; Harrison, P.A.; Berry, P.; Palomo, I.; Gómez-Baggethun, E.; Iniesta-Arandia, I.; Montes, C.; García del Amo, D.; Martín-López, B. What can conservation strategies learn from the ecosystem services approach? Insights from ecosystem assessments in two Spanish protected areas. Biodivers. Conserv. 2016, 27, 1575–1597. [Google Scholar] [CrossRef] [Green Version]
  36. Castro, A.J.; Martín-López, B.; López, E.; Plieninger, T.; Alcaraz-Segura, D.; Vaughn, C.C.; Cabello, J. Do protected areas networks ensure the supply of ecosystem services? Spatial patterns of two nature reserve systems in semi-arid Spain. Appl. Geogr. 2015, 60, 1–9. [Google Scholar] [CrossRef]
  37. Carvalho-Santos, C.; Monteiro, A.; Arenas-Castro, S.; Greifeneder, F.; Marcos, B.; Portela, A.; Honrado, J. Ecosystem Services in a Protected Mountain Range of Portugal: Satellite-Based Products for State and Trend Analysis. Remote Sens. 2018, 10, 1573. [Google Scholar] [CrossRef] [Green Version]
  38. Yu, B.; Rao, E.; Chao, X.; Shi, J.; Zhang, C.; Xu, W.; Xiao, Y.; Ouyang, Z. Evaluating the effectiveness of nature reserves in soil conservation on Hainan Island. Acta Ecol. Sin. 2016, 36, 3694–3702. [Google Scholar]
  39. Sharma, B.; Rasul, G.; Chettri, N. The economic value of wetland ecosystem services: Evidence from the Koshi Tappu Wildlife Reserve, Nepal. Ecosyst. Serv. 2015, 12, 84–93. [Google Scholar] [CrossRef]
  40. Yu, D.; Han, S. Ecosystem service status and changes of degraded natural reserves – A study from the Changbai Mountain Natural Reserve, China. Ecosyst. Serv. 2016, 20, 56–65. [Google Scholar] [CrossRef]
  41. Kibria, A.S.M.G.; Behie, A.; Costanza, R.; Groves, C.; Farrell, T. The value of ecosystem services obtained from the protected forest of Cambodia: The case of Veun Sai-Siem Pang National Park. Ecosyst. Serv. 2017, 26, 27–36. [Google Scholar] [CrossRef]
  42. Jia, J. The characteristics of formation, development and evolution of National Protected Areas in China. Int. J. Geoherit. Parks 2019, 7, 24–32. [Google Scholar] [CrossRef]
  43. Wang, W.; Xin, L.; Du, J.; Chen, B.; Liu, F.; Zhang, L.; Li, J. Evaluating conservation effectiveness of protected areas: Advances and new perspectives. Biodiv. Sci. 2016, 24, 1177–1188. [Google Scholar] [CrossRef] [Green Version]
  44. Zhang, Y.; Hu, Z.; Qi, W.; Wu, X.; Bai, W.; Li, L.; Ding, M.; Liu, L.; Wang, Z.; Zhang, D. Assessment of protection effectiveness of nature reserves on the Tibetan Plateau based on net primary production and the large-sample-comparison method. Acta Geogr. Sin. 2015, 70, 1027–1040. [Google Scholar]
  45. Xu, W.; Ouyang, Z.; Huang, H.; Wang, X.; Miao, H.; Hu, Z. Priority analysis on conserving China’s terrestrial ecosystems China. Acta Ecol. Sin. 2006, 26, 271–280. [Google Scholar]
  46. George, B.S.; Kang, A.; Hashi, T.; Cai, P. A winter wildlife survey in the northern Qiangtang of Tibet Autonomous Region and Qinghai Province, China. Acta Theriol. Sin. 2007, 27, 309–316. [Google Scholar]
  47. Song, X. Dynamic monitoring of Tibetan antelope habitat suitability in the Hoh Xil Nature Reserve using remote sensing images. Resour. Sci. 2016, 38, 17–23. [Google Scholar] [CrossRef]
  48. Hu, Y.; Li, W.; Jiang, Z.; Liu, W.; Liang, J.; Lin, Y.; Huang, Z.; Qin, H.; Jin, K.; Hu, H. A wild yak survey in Chang Tang of Tibet Autonomous Region and Hoh Xil of Qinghai Province. Biodiv. Sci. 2018, 26, 185–190. [Google Scholar] [CrossRef] [Green Version]
  49. Ma, S. Current situation investigation of grassland resource in Tanggula Mountain and Kekexili in Qinghai-Tibet Plateau. Pratacul. Sci. 2007, 9, 15–19. [Google Scholar]
  50. Yao, X.; Liu, S.; Li, L.; Sun, M.; Luo, J. Spatial-temporal characteristics of lake area variations in Hoh Xil region from 1970 to 2011. J. Geogr. Sci. 2014, 24, 689–702. [Google Scholar] [CrossRef] [Green Version]
  51. Liu, W.; Xie, C.; Wang, W. Analysis on expansion trend and outburst risk of the Yanhu Lake in Hoh Xil region, Qinghai-Tibet Plateau. J. Geogr. Sci. 2019, 41, 1–12. [Google Scholar]
  52. Xie, G.; Zhang, C.; Zhen, L.; Zhang, L. Dynamic changes in the value of China’s ecosystem services. Ecosyst. Serv. 2017, 26, 146–154. [Google Scholar] [CrossRef]
  53. Li, S.; Bing, Z.; Jin, G. Spatially Explicit Mapping of Soil Conservation Service in Monetary Units Due to Land Use/Cover Change for the Three Gorges Reservoir Area. China Remote Sens. 2019, 11, 468. [Google Scholar] [CrossRef] [Green Version]
  54. Sharp, R.; Tallis, H.T.; Ricketts, T.H.; Guerry, A.D. InVEST 3.2.0 User’s Guide; The Natural Capital Project, Stanford University, University of Minnesota, The Nature Conservancy, and World Wildlife Fund: Gland, Switzerland, 2015. [Google Scholar]
  55. Xie, G.; Zhang, C.; Zhang, L.; Chen, W.; Li, S. Improvement of the Evaluation Method for Ecosystem Service Value Based on Per Unit Area. J. Nat. Resour. 2015, 30, 1243–1254. [Google Scholar]
  56. Li, F.; Zhang, S.; Yang, J.; Chang, L.; Yang, H.; Bu, K. Effects of land use change on ecosystem services value in West Jilin since the reform and opening of China. Ecosyst. Serv. 2018, 31, 12–20. [Google Scholar]
  57. Chen, W.; Zhao, H.; Li, J.; Zhu, L.; Wang, Z.; Zeng, J. Land use transitions and the associated impacts on ecosystem services in the Middle Reaches of the Yangtze River Economic Belt in China based on the geo-informatic Tupu method. Sci. Total. Environ. 2019, 701, 134690. [Google Scholar] [CrossRef]
  58. Li, G.; Fang, C.; Wang, S. Exploring spatiotemporal changes in ecosystem-service values and hotspots in China. Sci. Total. Environ. 2016, 545, 609–620. [Google Scholar] [CrossRef]
  59. Xu, L.; Li, J.; Liu, C.; Cui, X. Research on Geomorphological Morphology and Regionalization of Hoh Xil Based on Digital Elevation Model (DEM). Acta Sci. Nat. Univ. Pekin. 2017, 53, 833–842. [Google Scholar]
  60. Li, R.; Duo, H.; Shi, L.; Zhou, Y.; Jiao, S.; Qiupei, Z.; Lei, G. Characteristics of climate change in Hoh Xil and its impact on landscape pattern during 1970—2013. J. Beijing Fores. Univ. 2015, 37, 59–68. [Google Scholar]
  61. Wu, D.; Li, Y.; Jiang, D.; Huang, Y.; Yan, H. Research of Remote Sensing Monitoring on Human Activity Impacts on Hoh Xil Nature Reserve. J. Gansu Sci. 2015, 27, 73–84. [Google Scholar]
  62. Liu, J.; Kuang, W.; Zhang, Z.; Xu, X.; Qin, Y.; Ning, J.; Zhou, W.; Zhang, S.; Li, R.; Yan, C.; et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
  63. Ning, J.; Liu, J.; Kuang, W.; Xu, X.; Zhang, S.; Yan, C.; Li, R.; Wu, S.; Hu, Y.; Du, G.; et al. Spatiotemporal patterns and characteristics of land-use change in China during 2010–2015. J. Geogr. Sci. 2018, 28, 547–562. [Google Scholar] [CrossRef] [Green Version]
  64. Xu, X. The Annually Normalized Difference Vegetation Index (NDVI) Spatial Distribution Datasets for China; Data Registration and Publishing System of the Resource and Environmental Science Data Center of the Chinese Academy of Sciences: Beijing, China, 2018. [Google Scholar]
  65. Wang, X.; Bao, Y. Study on the methods of land use dynamic change research. Prog. Geogr. 1999, 18, 83–89. [Google Scholar]
  66. Kogan, F.N. Remote sensing of weather impacts on vegetation in non-homogeneous areas. Int. J. Remote Sens. 1990, 11, 1405–1419. [Google Scholar] [CrossRef]
  67. Chen, C.; Dai, J.; Wang, H.; Liu, Y. Changes of the Value of Ecosystem Services in the Sanjiangyuan Region Based on Land Use Data. Prog. Geogr. 2012, 31, 970–977. [Google Scholar]
  68. Palomo, I.; Montes, C.; Martín-López, B.; González, J.A.; García-Llorente, M.; Alcorlo, P.; Mora, M.R.G. Incorporating the Social–Ecological Approach in Protected Areas in the Anthropocene. BioScience 2014, 64, 181–191. [Google Scholar] [CrossRef]
  69. Zaimes, G.N.; Gounaridis, D.; Symeonakis, E. Assessing the impact of dams on riparian and deltaic vegetation using remotely-sensed vegetation indices and Random Forests modelling. Ecol. Indic. 2019, 103, 630–641. [Google Scholar] [CrossRef]
  70. Peng, J.; Hu, X.; Zhao, M.; Liu, Y.; Lu, T. Research progress on ecosystem service trade-offs: From cognition to decision-making. Acta Geogr. Sin. 2017, 72, 960–973. [Google Scholar]
  71. Pei, F.; Wu, C.; Liu, X.; Li, X.; Yang, K.; Zhou, Y.; Wang, K.; Xu, L.; Xia, G. Monitoring the vegetation activity in China using vegetation health indices. Agr. Forest Meteorol. 2018, 248, 215–227. [Google Scholar] [CrossRef]
  72. Biggs, R.; Schluter, M.; Schoon, M.L. Principles for Building Resilience: Sustaining Ecosystem Services in Social-Ecological Systems; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar]
Remotesensing 12 00471 g001 550
Figure 1. Location of the study area in Zhiduo county, China. The Hoh Xil Nature Reserve (HXNR) is categorized into the core zone (A), the buffer zone (B), and the experimental zone by policies. The normalized difference vegetation index (NDVI) data set was downloaded from the Geospatial Data Cloud.
Figure 1. Location of the study area in Zhiduo county, China. The Hoh Xil Nature Reserve (HXNR) is categorized into the core zone (A), the buffer zone (B), and the experimental zone by policies. The normalized difference vegetation index (NDVI) data set was downloaded from the Geospatial Data Cloud.
Remotesensing 12 00471 g001
Remotesensing 12 00471 g002 550
Figure 2. Comparison of land use/land cover pattern in the Hoh Xil Nature Reserve and the control area for (a) 1980, (b) 1995, (c) 2008, and (d) 2018.
Figure 2. Comparison of land use/land cover pattern in the Hoh Xil Nature Reserve and the control area for (a) 1980, (b) 1995, (c) 2008, and (d) 2018.
Remotesensing 12 00471 g002
Remotesensing 12 00471 g003 550
Figure 3. Comparison of land use/land cover change in HXNR and the control area for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Figure 3. Comparison of land use/land cover change in HXNR and the control area for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Remotesensing 12 00471 g003
Remotesensing 12 00471 g004 550
Figure 4. Distribution of the total ecosystem service value (ESV) per unit in the study area for (a) 1980, (b) 1995, (c) 2008, and (d) 2018.
Figure 4. Distribution of the total ecosystem service value (ESV) per unit in the study area for (a) 1980, (b) 1995, (c) 2008, and (d) 2018.
Remotesensing 12 00471 g004
Remotesensing 12 00471 g005 550
Figure 5. Comparison of the changes of ecosystem service value per unit in HXNR and the control area for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Figure 5. Comparison of the changes of ecosystem service value per unit in HXNR and the control area for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Remotesensing 12 00471 g005
Remotesensing 12 00471 g006 550
Figure 6. Hotspot of the total ESV change for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Figure 6. Hotspot of the total ESV change for (a) 1980–1995, (b) 1995–2008, and (c) 2008–2018.
Remotesensing 12 00471 g006
Table
Table 1. The initial classification of Land use/land cover data.
Table 1. The initial classification of Land use/land cover data.
Primary Categories (Code)Subcategories Types (Code)
Cropland (1)Paddy field (11), dry land (12)
Woodland (2)Forested land (21), Shrubland (22), Sparse land (23), Other woodlands (24)
Grassland (3)High dense grassland (31), Moderate dense grassland (32), Sparse grassland (33)
Water bodies (4)River and canals (41), Lakes (42), Reservoir and pond (43), Permanent glaciers and snow (44), Intertidal zone (45), Floodplain (46)
Built-up areas (5)Urban areas (51), Rural housing areas (52), Other construction land (53)
Unused land (6)Sandy desert (61), Harsh/Gobi desert (62), Saline and alkaline land (63), Marshy land (64), Bare land (65), Bare rock (66), Other unused land (67)
Ocean (99)-
Data source: The Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn) [62,63].
Table
Table 2. The area of different land use/land cover types in the study area during 1980–2018.
Table 2. The area of different land use/land cover types in the study area during 1980–2018.
AreasLULC Area (km2)1980199520082018
The HXNRForest 7.20 (0.01) 10.19 (0)7.21 (0.01)7.22 (0.01)
Grassland 17,322.98 (35.31)17,773.43 (36.23)17,304.22 (35.27)17,616.34 (35.91)
Waterbodies 3543.29 (7.22)3412.15 (6.95)3584.37 (7.31)3580.71 (7.30)
Wetland723.93 (1.48)607.31 (1.24)675.65 (1.38)834.36 (1.70)
Urban0000
Desert27,466.26 (55.98)27,270.60 (55.58)27,492.25 (56.03)27,025.04 (55.08)
The Control AreaForest 4.03 (0.04)4.04 (0.04)4.06 (0.04)3.28 (0.03)
Grassland 9389.02 (87.95)10,133.31 (94.93)9379.32 (87.86)9289.67 (87.02)
Waterbodies 126.70 (1.19)121.27 (1.14)126.88 (1.19)128.44 (1.20)
Wetland50.75 (0.48)19.51 (0.18)50.60 (0.47)49.37 (0.46)
Urban0.80 (0.01)1.74 (0.02)0.81 (0.01)6.16 (0.06)
Desert1103.58 (10.34)395.01 (3.70)1113.22 (10.43)1197.97 (11.22)
1 (0.01): Percentage of each land use/land cover type (%), others as same. The total area of the HXNR and the control area is 49,063.68 km2 and 10,674.88 km2, respectively.
Table
Table 3. Land use dynamics in the study area during 1980–2018 and for different periods before and after the establishment of the HXNR.
Table 3. Land use dynamics in the study area during 1980–2018 and for different periods before and after the establishment of the HXNR.
LULC Types (%)1980–19951995–20082008–20181980–2018
SoH 1SoC 2SoHSoCSoHSoCSoHSoC
Forest−6.490284.510.040.01−1.920.01−0.54
Grassland0.170.53−0.20−0.570.18−0.100.05−0.03
Waterbodies−0.25−0.290.390.36−0.010.120.030.04
Wetland−1.07−4.100.8712.252.35−0.240.44−0.08
Urban07.800−4.14066.51019.08
Desert−0.05−4.280.0613.99−0.170.76−0.050.24
1 SoH: Single land use dynamic degree of the HXNR; 2 SoC: Single land use dynamic degree in the control area.
Table
Table 4. Land use/land cover conversions in the study area for different periods before and after the establishment of the HXNR.
Table 4. Land use/land cover conversions in the study area for different periods before and after the establishment of the HXNR.
Conversion Types (km2)The HXNRThe Control Area
1980–19951995–20082008–20181980–19951995–20082008–2018
Forest Translations7.087.080.510.010.120.61
Grassland to Waterbodies/Wetland10.2684.8354.804.5944.985.67
Urban Expansion---0.960.035.51
Desert Expansion1312.601548.07368.1670.89781.54159.84
Waterbodies/Wetland to Grassland 79.1214.2677.2644.355.114.54
Transitions between Waterbodiesaters & Wetland51.8933.05400.550.190.190.35
Desert to Green1508.271326.43839.86778.8663.3273.21
Table
Table 5. Ecosystem service value in the study area during 1980–2018 and the annual rates of change for three periods before and after the establishment of the HXNR.
Table 5. Ecosystem service value in the study area during 1980–2018 and the annual rates of change for three periods before and after the establishment of the HXNR.
AreasESV
(million USD)		19801995200820181980–1995
(%)		1995–2008
(%)		2008–2018
(%)
The HXNRProvisioning38.0038.8637.9438.600.15−0.180.17
Regulating447.07441.41445.77456.61−0.080.070.24
Supporting210.67214.75210.39213.150.13−0.160.13
Cultural68.3767.8268.2569.56−0.050.050.19
Total764.11762.85762.36777.92−0.01−0.0040.20
The Control AreaProvisioning45.1747.5045.1344.790.34−0.40−0.08
Regulating362.53368.49362.14359.760.11−0.13−0.06
Supporting236.97247.73236.73235.130.30−0.35−0.07
Cultural53.4454.2753.3953.110.10−0.13−0.06
Total698.11718.00697.38692.780.19−0.22−0.07
Table
Table 6. The average ESV per unit in the study area during 1980–2018.
Table 6. The average ESV per unit in the study area during 1980–2018.
ESV per Unit
(USD per ha)		The HXNRThe Control Area
19801995200820181980199520082018
Provisioning7.757.927.737.8742.3144.5042.2741.96
Regulating91.1289.9790.8693.09339.61345.20339.25337.01
Supporting42.9443.7742.8843.45221.99232.06221.76220.26
Cultural13.9313.8213.9114.1850.0650.8450.0149.75
Total155.74155.48155.38158.59653.97672.60653.29648.98
Table
Table 7. Statistics for hotspots and cold spots of ESV increases at a confidence level of 95% in the HXNR during 1980–2018.
Table 7. Statistics for hotspots and cold spots of ESV increases at a confidence level of 95% in the HXNR during 1980–2018.
Periods ESV Changed PatchesNumber of Hot SpotsNumber of Cold Spots
A 1B 2AB
1980–19953,281,059307,324 (9) 31,311,912 (40)1,016,464 (31)437,961 (13)
1995–20083,334,3181,032,753 (31)435,893 (13)330,902 (10)1,338,828 (40)
2008–20181,851,733527,571 (28)450,687 (24)189,580 (10)337,833 (18)
1980–20181,939,336539,118 (28)455,401 (23)211,097 (11)373,246 (19)
1 A: Core zone, 2 B: Buffer zone, 3 (9): Percentage of the ESV changed patches (%), others as same.

Share and Cite

MDPI and ACS Style

Zeng, J.; Chen, T.; Yao, X.; Chen, W. Do Protected Areas Improve Ecosystem Services? A Case Study of Hoh Xil Nature Reserve in Qinghai-Tibetan Plateau. Remote Sens. 2020, 12, 471. https://doi.org/10.3390/rs12030471

AMA Style

Zeng J, Chen T, Yao X, Chen W. Do Protected Areas Improve Ecosystem Services? A Case Study of Hoh Xil Nature Reserve in Qinghai-Tibetan Plateau. Remote Sensing. 2020; 12(3):471. https://doi.org/10.3390/rs12030471

Chicago/Turabian Style

Zeng, Jie, Tianyang Chen, Xiaowei Yao, and Wanxu Chen. 2020. "Do Protected Areas Improve Ecosystem Services? A Case Study of Hoh Xil Nature Reserve in Qinghai-Tibetan Plateau" Remote Sensing 12, no. 3: 471. https://doi.org/10.3390/rs12030471

APA Style

Zeng, J., Chen, T., Yao, X., & Chen, W. (2020). Do Protected Areas Improve Ecosystem Services? A Case Study of Hoh Xil Nature Reserve in Qinghai-Tibetan Plateau. Remote Sensing, 12(3), 471. https://doi.org/10.3390/rs12030471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop