Dynamic Analysis of Ecological Environment Quality Combined with Water Conservation Changes in National Key Ecological Function Areas in China

The shortage of water resources is a key factor limiting the sustainability of the economy and society. Most of the 25 National Key Ecological Function Areas (NKEFAs) in China serve as a source and supplementation for numerous rivers and playing an important role in water resource conservation. Based on the analysis of eco-environmental quality changes in NKEFAs, this study analyzed the spatial pattern of water conservation services in 2000 and 2010 by using a water balance equation. The results indicate that the land cover type of NKEFAs was dominated by grassland, and the proportion of ecological land conversion to non-ecological land (0.3%) was higher than that of non-ecological land conversion to ecological land (0.21%). The fractional vegetation coverage (FVC) and biomass density of NKEFAs gradually decreased from southeast to northwest. The FVC of the Changbai Mountain Forest Function Area (CBS) was the highest, while the biomass density and total biomass were highest in mountain areas in the Middle of Hai’nan Island (HND) and in the Great Khingan and Lesser Khingan Mountains (XAL) respectively. The FVC and biomass of NKEFAs mostly increased in 2000–2010. Water conservation amounts of NKEFAs decreased from southeast to northwest. The average water conservation and total water conservation amount of Nanling Mountain (NL), Guangxi-Guizhou-Yunnan (GQD), and the Wuling Mountain Function Area (WLS) were the highest, while the Yinshan Mountain (YS), Alkin Grassland (AEJ), and the Qilian Mountain Function Area (QLS) had the lowest values. In 2000–2010, the water conservation service of 60% of NKEFAs decreased. Spatial and temporal differences in water conservation services are the result of a combination of ecological environment quality and meteorological conditions. Protection of the ecological environment and vegetation coverage improvement should be strengthened to enhance the function of water conservation.


Introduction
National Key Ecological Function Areas (NKEFAs) refer to the areas where ecosystems play an important role in ensuring the ecological security of an entire country or a large region and where the intensity of development should be restricted because ecosystems are degrading and the capacity to supply ecological products should be maintained and improved. Since 2008, the central government has set up transfer payments for NKEFAs. From 2008 to 2014, the central government allocated 200.4 billion CNY for the transfer payments. The number of districts and counties of NKEFAs has increased to 819 and transfer payments reached 62.7 billion CNY in 2017. An ecological environment index (EI) was adopted to evaluate the annual change of the ecological environment in each city and county, but has low applicability to water conservation. Chinese researchers have conducted research on NKEFAs [1] and have tended to focus on ecological compensation mechanisms [2,3] and transfer payment issues [4][5][6][7].
With the increasing demand for global water resources and the rapid deterioration of water environments, water shortage has become a global issue as well as a key factor limiting economic and social development. China faces a large number of water problems, including river pollution, over-exploitation of groundwater and water shortages. Most of the 25 key ecological function areas in China are located in Level I and II stepped landforms at higher altitudes, serving as a source and supplementation for numerous rivers and playing an important role in water resource conservation. However, the studies primarily focused on the water conservation of forest typed key ecological function areas or single ecological function area [8][9][10] and were mostly qualitative [1]. Besides, ecosystem services (ESs) are benefits people derive from ecosystems [11] and the ecological environment of ecosystem is the basis for ESs, including the water conservation service. Ecological environment quality can be enhanced by increasing vegetation cover, which can enhance water conservation. But most studies just evaluate the spatial-temporal pattern of water conservation service and can't establish a dynamic relationship between ecological environment quality and water conservation changes. On the whole, there is a lack of systematic, targeted, comprehensive analyses and evaluations of ecological environment quality and water conservation services across the NKEFAs in China. Therefore, the scientific evaluation of ecological environment quality and water conservation in NKEFAs will provide a reliable theoretical basis for ecological compensation and targeted protection measures, and will have great practical implications.
The ecological environment quality is a combined result of multiple factors, including natural and artificial, which can be reflected by climate, terrain, biomass, land cover and vegetation data based on remote sensing interpretation. The understanding of water conservation services varies with scholars. Currently, the most comprehensive research concerns the forest ecosystem water conservation function, which involves the function of effective water conservation and runoff adjustment for the forest ecosystem through the interception and redistribution of precipitation by forest canopies, litter layers and soil layers [12]. The water conservation quantity can be calculated by numerous methods, such as the soil water storage capacity method, water balance method, underground runoff increase method, precipitation storage method and forest canopy interception residual method [13]. A large number of studies have proven that the regional water balance method is the basis for studying water conservation mechanisms. The method can accurately calculate the water conservation capacity, is easy to perform and is applicable for all temporal and spatial scales [14]. It is currently the most effective and widely used method for calculating the water conservation function [15]. Therefore, this study uses the water balance equation to calculate the water conservation quantity, including water content in soil, litter layer water holding capacity and canopy interception.
To eliminate the influence of transfer payment policy and get more objective evaluation results of ecological environment quality and water conservation changes in NKEFAs in China, we need to clarify the ecological background of NKEFAs before the implementation of transfer payments first, namely, the average state of a regional ecosystem and its changes before the implementation of ecological protection and restoration. This average state is more referential to the formulation of ecological compensation policy. Therefore, we analyzed the change and relationship of ecological environment quality and water conservation in NKEFAs in 2000 and 2010 with the aim of providing a theoretical basis for ecological protection and ecological compensation based on climate, terrain, biomass, land cover and vegetation data.

Study Areas
There are 25 key ecological function areas in China. The total area of NKEFAs in China is 3.81 × 10 6 km 2 , which are located in most provinces except for the eastern coastal ones (Figure 1 and Table 1), and present all climate types in China, including tropical monsoon, subtropical monsoon, temperate monsoon, temperate continental and plateau mountain climates. The annual average temperature in these areas varies from −2.8 to 27 • C and declines gradually from south to north. Annual precipitation range from 15 mm to 1600 mm and declines gradually from southeast to northwest. Vegetation in these areas includes primarily cold and temperate coniferous forest, temperate broadleaf-conifer forest, broadleaved deciduous forest, evergreen broadleaved forest, evergreen broadleaved deciduous forest, monsoon evergreen broadleaved forest, monsoon forest, rainforest, temperate steppe, temperate desert, high-cold scrub and meadow, high-cold steppe and high-cold desert. Types of soil include podzolic coniferous forest soil, dark brown forest soil, brown forest soil, yellow brown soil, red soil, chernozem, chestnut soil, grey brunisolic soil, plateau meadow soil, plateau steppe soil and plateau desert soil.

Data Sources
(1) Data on key ecological function areas were produced using 1:1,000,000 national county-level data and the directory of NKEFAs in National Main Functional Areas Planning, 2010. Network (CERN) from the study sites, and consisted of 1146 records, including forest, shrub, grass, and farmland samples. Moreover, all species of shrub and parts of herb layers were from forest communities. The survey data were collected in 2004, when the vegetation was at the accumulative process of the growth stage. The information including vegetation type, actual fresh and dry weight (including branches, limbs, leaves, and roots), selected site, date, latitude, longitude, elevation, annual average temperature, and annual precipitation were recorded.

Transition Matrix for Land Cover Change Detection
The raster land cover maps were used to detect the internal conversions of land cover in NKEFAs from 2000 to 2010 and an extended transition matrix was constructed. The maps from the initial and subsequent time periods were overlaid to produce a matrix that provided the LUCC areas by categorical transition between 2000 and 2010. The off-diagonal entries comprise the proportions of the landscape that experienced a transition from one category to another while the on-diagonal entries indicated no change in categories. The row totals at the right denote the proportion of landscape by land cover category in 2000 and the column totals at the bottom denote the proportion of landscape by category in 2010.

Estimation of Biomass
In this study, the above-and below-ground biomass (ABGB) across NKEFAs was estimated with MLR to establish the correlation between field measurements and MODIS observations, meteorological data, coordinates, terrain data, and statistical data [16]. MLR is a supervised method that aim to establish a mathematical relationship between a property of a given system and a set of molecular characteristics or descriptors that encode information and is expressed in Equation (1): where ε is an n × 1 residuals vector; X is a known n × k matrix of description; A is a K × 1 vector of adjusted parameters; and Y is a n × 1 vector of the response variable related to either the activity or other system property [17]. The model performance for MLR was assessed based on the agreements between the predicted value and the observed value. The agreements were quantified using relative estimation error (REE), which is calculated as Equation (2) [18]: where Y i is the observed data, Y i is the predicted value, and N is the number of validation points. Ninety percent of the biomass files were selected through a hierarchical sample method according to different classes for establishing regression models and the remainder 10% of samples were used to calculate REE. MLR was calculated for the terrestrial ecosystem biomass of China in 2010. The simulation results were compared with the survey data, and the simulation precision was calculated using Equation (3) [19]: where V m is the value of simulated results by the regression model, and V s is the value of survey data. The selected variables were considered as influence factors, such as geographic elements (including longitude, latitude, and elevation) and environmental factors (including annual precipitation and annual average temperature) influencing the growth of vegetation, statistical data, and remote-sensing vegetation index (including NDVI and LAI) that directly showed the status of vegetation on a large scale. The construction of MLR for biomass assessment was executed through a step-wise regression analysis using the eight predictive variables and the biomass variable. Different types have variable correlations. Finally, before variables can be considered, the correlation coefficients between biomass data and other factors must be tested. We divided them into eight groups according to different vegetation types, including evergreen needle-leaved forest (EN), evergreen broadleaved forest (EB), deciduous coniferous forest (DC), deciduous broad-leaved forest (DB), mixed forest (MF), shrub, grass, and farmland. We used the data analysis add-in, available in the software of SPSS, and different types have different correlations with different factors (Table 2). Correlation analysis was conducted between biomass and eight factors. The linear relationships for the four primary types are given in Table 3. The R 2 value of DC forest was higher than other types, and grass had a lower R 2 value. More than 95% of the variance in ABGB density was explained, with a relative estimation error of 67 g·m −2 for a range of biomass density in 2010 (Figure 2), when the residual data were used for training and cross-validation ( Figure 2).  Figure 2. Relationship between observed biomass (g/m 2 ) and multiple regression model predictions (g/m 2 ).

Fractional Vegetation Cover
The annual fractional vegetation coverage (FVC) data were obtained according to the dimidiate pixel model based on NDVI. Monthly NDVI values were calculated using the maximum value composite (MVC) method, which minimizes the following factors: atmospheric effects, scan angle effects, cloud contamination, and solar zenith angle effects [20]. Annual NDVI values comprised the maximum value of the monthly NDVI datasets. According to the dimidiate pixel model, the NDVI value of a pixel is a combination of information contributed by vegetation and non-vegetation. The calculation formula of FVC is as follows: where NDVI veg is the NDVI value of a pure vegetation pixel, and NDVI soil is the NDVI value of a non-vegetation pixel. In this study, the NDVI value of 95% and 5% cumulative frequency was NDVI veg and NDVI soil , respectively. In this study, FVC was divided into five grades: 0-0.2, 0.2-0.4, 0.4-0.6, 0.6-0.8, and 0.8-1.0, representing extremely low coverage, low coverage, medium coverage, high coverage and extremely high coverage respectively.

Water Conservation Amount
This study uses the following water balance equation to calculate the water conservation capacity: where WCON denotes the total quantity of water conservation (m 3 ), P i denotes precipitation (mm), R i denotes surface runoff (mm), ET i denotes evapotranspiration (mm), A i denotes the area of the ecosystem, i denotes the ith ecosystem type and j denotes the ecosystem type number. For areas where surface runoff (R i ) is the product of precipitation and the surface runoff coefficient, the formula is as follows: where R is surface runoff (mm), P is precipitation (mm) and α is the average surface runoff coefficient. Based on the regional conditions of NKEFAs in combination with data from the literature, we obtain the average runoff coefficient for the various ecosystems in NKEFAs (Appendix A) [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39]. Actual evapotranspiration is calculated using the Technical Guidelines for the Redline Delimitation of The Ecological Protection (Formula (7)): where P is multi-annual mean precipitation (mm), ET is actual evapotranspiration (mm), ET 0 is multi-annual mean latent evapotranspiration, and ω is the underlying surface (land cover) impact coefficient from the Technical Guidelines for the Redline Delimitation of The Ecological Protection (Table 4). The latent evapotranspiration is calibrated with the FAO56Penman-Monteith formula [40] based on daily meteorological data: where R n is surface net radiation (MJ·m −2 ·d −1 ) from crops, G is the soil heat flux density (MJ·m −2 ·d −1 ), T is monthly mean temperature ( • C), U 2 is 2 m wind velocity (m·s −1 ), e s is saturated water vapour pressure (kPa), e a is the actual water vapour pressure (kPa), ∆ is the slope of the saturated water vapour-temperature curve (kPa· • C −1 ), and γ is the psychrometric constant (kPa· • C −1 ).

Land Cover Change Analysis
The land cover type of NKEFAs is dominated by grassland, covering 40.83% and 40.7% of NKEFAs in 2000 and 2010, respectively, followed by sand and bare land (22.79%, 22.77%), forest land (18.53%, 18.58%), farmland (8.31%, 8.38%), bush land (4.74%, 4.75%), wetland and water areas (4.47%, 4.46%), rural settlements (0.26%, 0.27%), and urban settlements (0.06%, 0.09%). Among them, forest land, bush land, and grassland are categorized as ecological land, and their ecological levels decrease successively. Farmland, rural settlements, urban settlements, sand and bare land are categorized as non-ecological land, and the ecological level decreases successively [41]. 99.17% of land cover types did not change between 2000 and 2010 (Table 5). Although the land cover change is relatively small, the conversion of ecological land to non-ecological land accounts for 0.3% of the total area, which is higher than the proportion of non-ecological land conversion to ecological land (0.21%). For ecological land, urban settlements, rural settlements and farmland increased by 41.73%, 3.24% and 0.84% respectively. The proportion of urban settlements increased significantly, which is mainly from the conversion of farmland, forest land and grassland. In non-ecological land, wetland and water areas and grassland were both reduced by 0.32%. The wetland and water areas were mainly converted into grassland and farmland, and the grassland was mainly converted into farmland, sand areas and bare land. Forest land and bush land increased by 0.26% and 0.19% respectively, which mainly came from the conversion of grassland, farmland, and bush land and grassland, and forest respectively.
The land cover of 21 NKEFAs were dominated by ecological land (Table 6), among which the land cover types did not change between 2000 and 2010 in ZDN and ZXB. The proportion of non-ecological land conversion to ecological land was higher than that of ecological land conversion to non-ecological land in GN, GQD, HTGY, SX, YS, CD, QB, and WLS. However, the urban settlements in these function areas all increased, and the growth rate was highest for all land cover change rates, at over 40%. The growth rates of forest land in SX, GQD, WLS was 81.62%, 33.55%, and 22.46% respectively, which were higher than in other function areas. The proportion of ecological land conversion to non-ecological land was higher than that of non-ecological land conversion to ecological land in XAL, REG, CBS, SJY, DBS, HLBR, HSDK, KEQ, HND, NL, and SJPY and these area had the highest growth rate of urban settlements, at over 20%. Among them, the increase of urban settlements in SJPY and HLBR were lower. The sand areas and bare land in SJPY increased by 12.5%, and grassland decreased by 11.56%, while bush land decreased by 10.51% in HLBR. Land cover type was mainly non-ecological land in QLS, AEJ, AET, and TLM, among which the proportion of non-ecological land conversion to ecological land was higher than that of ecological land conversion to non-ecological land in QLS and AEJ. The increase of urban settlements was the highest, with the increase of the urban settlements of AEJ the most substantial, up 38.75 times. The proportion of ecological land conversion to non-ecological land was higher than that of non-ecological land conversion to ecological land in AET and TLM. Farmland increased up to 28.27% in AET, and the increase of urban settlements dominated land cover change in TLM.  AET TLM Note: 1-Land cover was dominated by ecological land and the proportion of non-ecological land conversion to ecological land was higher than that of ecological land conversion to non-ecological land between 2000 and 2010. 2-Land cover was dominated by ecological land and had no change between 2000 and 2010. 3-Land cover was dominated by ecological land, and the proportion of ecological land conversion to non-ecological land was higher than that of non-ecological land conversion to ecological land between 2000 and 2010. 4-Land cover was dominated by non-ecological land, and the proportion of non-ecological land conversion to ecological land was higher than that of ecological land conversion to non-ecological land between 2000 and 2010. 5-Land cover was dominated by non-ecological land, and the proportion of ecological land conversion to non-ecological land was higher than that of non-ecological land conversion to ecological land between 2000 and 2010.

Ecological Environment Quality Change Trends
Vegetation is the main body of terrestrial ecosystem and has an important role in soil and water conservation, atmospheric regulation, mitigation of greenhouse gas emissions, and in maintaining the climate and ecosystem stability. Change of vegetation cover is a direct result of change in the ecological environment, and is significant in evaluating the area's ecological environment quality. The FVC of NKEFAs gradually decreased from southeast to northwest due to the influence of climate and other factors ( Figure 3) and was highest in CBS, reaching 0.962 in 2000 and 0.959 in 2010. The FVC of XAL, NL, DBS, WLS, SX, QB, and HND were higher than 0.9 (listed here in descending order), areas that were primarily water conservation and biodiversity maintenance ecological function areas in the northeast or southeast. The FVC was 0.039 in 2000 and 0.059 in 2010 in AEJ, which was the lowest, and the FVC was also lower than 0.2 in ZXB and TLM, a value that corresponds to the extremely low coverage area. In 2000-2010, the FVC of 80% of NKEFAs increased, and KEQ and HSDK had the highest growth rate in vegetation coverage with increases of 48.7% and 48.2%, respectively. The FVC of CBS, YS, ZDN, CD and HND were slightly reduced by less than 5.1%. Biomass is a key variable that represents the vegetation in a biological community, and directly reflects the supply capacity of the ecosystem in the natural environment. The spatial distribution of biomass in NKEFAs is not uniform, and the biomass density was higher in the southeast than the northwest (Figure 4). The total biomass of terrestrial vegetation in NKEFAs was 7.61 Pg (1 Pg = 1015 g) in 2000, with an average biomass density of 2.22 kg/m 2 . The biomass of XALwas the largest in 2000, reaching 2.22 Pg and accounting for 29.2% of the total biomass of the whole NKEFAs, followed by CD (1.03 Pg). The biomass of YS was the lowest at 0.02 Pg. The function areas with large biomass density were mainly located in the south. The biomass density of HND was the highest in 2000, reaching 11.95 kg/m 2 , followed by ZDN (8.68 kg/m 2 ). The biomass density of AEJ was 0.10 kg/m 2 , which was the lowest. In 2010, the total biomass of terrestrial vegetation in NKEFAs was 8.96 Pg, with an average biomass density of 2.42 kg/m 2 , which was higher than 2000 indicating that the total biomass quality of NKEFAs improved from 2000 to 2010. The total biomass levels of XAL and CD were the highest, reaching 2.65 Pg and 1.19 Pg respectively. The biomass of YS was the lowest at 0.03 Pg, which was higher than in 2000. Biomass density was high in HND and ZDN, which were 12.67 kg/m 2 and 8.68 kg/m 2 , respectively in 2010, and the lowest biomass density in 2010 was 0.14 kg/m 2 in AEJ in 2010. According to the types of ecological function areas, the vegetation coverage of soil and water conservation function areas was the highest, the average biomass density of the biodiversity maintenance function areas was the highest, and the total biomass of water conservation function areas was the highest. Windbreak and sand-fixation function areas are mainly located in the northwest inlands, where desert is the main land cover type, leading to the lowest vegetation coverage, biomass density and total biomass ( Figure 5). On the whole, the ecosystem quality of NKEFAs improved from 2000 to 2010.

Water Conservation Services Change Trends
The water conservation quantity of NKEFAs showed a trend of significant reduction from southeast to northwest (Figures 6 and 7). The average water conservation amount of each function area was different, which was mainly caused by various factors such as ecosystem type, precipitation, vegetation coverage, vegetation type and so on. Precipitation decreased significantly from the southeast to the northwest, and, correspondingly, vegetation coverage gradually decreased as the vegetation type changed from forest and bush to desert grassland. In 2000, the average water conservation amount of each NKEFAs was 124.28 mm, and the total amount of water conservation was 4674.98 × 10 8 m 3 . The average water conservation of HND was the highest, reaching 1753.07 mm, followed by NL, GQD, WLS, DBS and SX, which had average water conservation levels between 500 mm and 1000 mm. The average water conservation of other function areas was lower, especially for YS, AEJ and QLS, which had negative average water conservation levels. The total water conservation amount of CD, 1036.28 × 10 8 m 3 , was the highest, which was almost twice as much as NL at 652.47 × 10 8 m 3 , followed by GQD, WLS, QB, XAL, ZXB, and SJY, whose total water conservation amounts were between 200 × 10 8 m 3 and 600 × 10 8 m 3 . The total water conservation amount of other function areas was small, especially YS, AEJ, QLS. In 2010, the average water conservation in each function area was reduced to 107.39 mm, and the total water conservation amount was 4048.52 × 10 8 m 3 . Similar to 2000, the average water conservation of HND was the highest at 1095.99 mm, followed by NL and WLS and the negative average water conservation of YS, AEJ, TLM, and AET were the lowest of the function areas. The total water conservation amount of CD, 888.31 × 10 8 m 3 , was the highest. The total water conservation amount in QB, XAL, NL, SJY, WLS, and GQD decreased successively, and the total water conservation amounts of AEJ, YS, TLM, AET were the lowest. The difference of water conservation in each function area is mainly due to the different water conservation capacity of various land cover types and the land cover type composition.
On the whole for 2010, the average water conservation and total water conservation of NL, GQD and WLS was higher than in 2000, while for YS, AEJ, and QLS both were lower and the rank of both in HND and CD differed significantly. The ranking order of average water conservation and total water conservation in other NKEFAs also shows a degree of inconsistency, mainly because of the different size of the total water conservation and the difference of average water conservation amount. In 2000 and 2010, the function of water conservation service in 60% of NKEFAs decreased.

Discussion
Water conservation services are determined by both ecological environment quality and meteorological conditions. The FVC of CBS and ZDN declined from 2000 to 2010, meanwhile, precipitation increased, and actual evapotranspiration decreased. However, the average and total water conservation amount both decreased, indicating that vegetation cover played a leading role in water conservation service in this time period. In contrast, the FVC, biomass, and the proportion of ecological land in GQD and SX all increased from 2000 to 2010, while precipitation decreased, and actual evapotranspiration increased for the same time period. However, the average and total water conservation amount both decreased, indicating that meteorological conditions played a dominant role in water conservation service in this time period.
In 2000 and 2010, the average water conservation amount of 25 key ecological function areas in China was 270.2 mm and 210.9 mm respectively. However, among the eight water-source conservation functional areas, only the average water conservation amount of NL and REG exceeded the average, which indicated that the water conservation function of water-source conservation functional areas was not the highest among NKEFAs and that the significance of these areas water conservation was reflected in the importance to the region and its surrounding areas. However, some water-source conservation functional areas faced serious threats. For example, the water conservation amounts of QLS and AET were negative in 2000 and 2010, respectively. It is necessary to strengthen the protection of the ecological environment and improve water conservation service especially in the eight water-source conservation functional areas.
In 2000 and 2010, the total water conservation amount of NKEFAs was 4674.98 × 10 8 m 3 and 4048.51 × 10 8 m 3 , respectively and the water conservation amount of the eight water-source conservation functional areas was 1381.95 × 10 8 m 3 and 1350.04 × 10 8 m 3 , respectively. Gong et al. [21] calculated that the total amount of water conserved by forest, shrub, grassland, garden, and wetland ecosystems in China in 2010 was 12224.33 × 10 8 m 3 . Huang et al. [42] calculated that the water conservation amount of eight water-source conservation functional areas was 1050.72 × 10 8 m 3 and 1081.56 × 10 8 m 3 in 2000 and 2010, respectively. The above results were comparable and had no distinct difference (Table 7). In addition, the spatial distribution of water conservation of NKEFAs was consistent with previous studies, showing a trend of significant reduction from southeast to northwest [21]. The water conservation amounts of XAL and SJY were much higher than those of GN, REG, QLS, and AET [42], which was consistent with this study. In this study, different land cover types were considered in the evaluation of the spatial distribution of water conservation, and the results were closer to the actual results. In addition, an analysis of land cover, FVC, biomass change combined with the water conservation calculation can provide a more scientific basis for ecological environment monitoring and natural resources management.

Conclusions
The ecological environment quality of NKEFAs improved between 2000 and 2010, and the average vegetation coverage and biomass increased, but the proportion of non-ecological land increased slightly. Ecological land serves as an important basis for the improvement of ecological environment quality and ecosystem function realization in NKEFAs, especially in AET and TLM whose non-ecological land are not only the dominating land cover types but also present an increasing trend. In terms of ecological environment quality of different function areas, the FVC of soil and water conservation ecological function areas were the highest, the biodiversity maintenance ecological function areas had the highest average biomass density, and the total biomass of water conservation ecological function areas were the highest. The wind prevention and sand fixation ecological function areas were mainly located in the northwest inland area, where the desert coverage was relatively larger than in other areas, leading to the lowest levels of vegetation coverage, biomass density and total biomass. In general, the spatial distribution of water conservation services showed a trend of significant reduction from southeast to northwest, which was similar to the spatial distribution pattern of ecological environment quality. NL, GQD and WLS had high water conservation amounts and good ecological environment quality, while in comparison, the levels of both for YS, AEJ and QLS were significantly inferior. Among the eight water-source conservation functional areas, SJY, QLS and AET had low water conservation amounts and effective environmental protection engineering measures should be taken in these areas. It was obvious that not all the water-source conservation functional areas had higher water conservation amount, and the water conservation services was the result of a combination of factors, including climate condition and ecological environment quality. However, water-source conservation functional areas are mainly located on the upstream of important river in China, the water conservation services of which are more important than other services for local and their downstream area. This can reflect the orientation and target for different types of NKEFAs. Besides, we should notice that the water conservation service in each type of functional areas are important that can't be ignored and the NKEFAs is an organic whole which is crucial to China's ecological environment quality improvement. These results provide specific information that may serve to strengthen necessary public awareness about protecting and restoring water conservation services and improving ecological environment quality of NKEFAs, issues that are also significant for natural resources management.