Next Article in Journal
Theoretical Analysis on the Short-Circuit Current of Inverter-Interfaced Renewable Energy Generators with Fault-Ride-Through Capability
Previous Article in Journal
Innovation Ecosystem of CNG Vehicles: A Case Study of Its Cultivation and Characteristics in Sichuan, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Heterogeneity of Typical Ecosystem Services and Their Relationships in Different Ecological–Functional Zones in Beijing–Tianjin–Hebei Region, China

1
College of Land Science and Technology, China Agricultural University, Beijing 100193, China
2
Key Laboratory for Agricultural Land Quality, Monitoring and Control, The Ministry of Land and Resources, Beijing 100193, China
3
Center for Land Resources Research in Northwest China, Shaanxi Normal University, Xi’an 710119, China
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(1), 6; https://doi.org/10.3390/su10010006
Submission received: 4 December 2017 / Revised: 17 December 2017 / Accepted: 18 December 2017 / Published: 25 December 2017
(This article belongs to the Section Environmental Sustainability and Applications)

Abstract

:
Recognizing changes in ecosystem services (ES) and their relationships is the basis of achieving sustainable regional development. Regional collaborative development has become the core strategy of the development of the Beijing–Tianjin–Hebei (BTH) region. However, sub regions have different ecological changes and relationships. Here, we quantify and map ES, including water yield, sediment retention, carbon sequestration and grain productive capacity in 2000, 2005, 2010 and 2015, using several biophysical models and explore the relationships of spatial correction, trade-offs and synergies among multiple ES in different spatial scales. Results across the four years show that the quality and variation tendency of ES from each region are spatially heterogeneous. The relationship between ES that are not significant in the entire region shows different correlations in individual ecological–functional zones. From the perspective of regional disparity, the effect of land use factor and correlative mechanisms among ES are analyzed. To observe the spatiotemporal variations and relationships of ES in individual regions, land use management policies are proposed on the basis of the results of the relationships among ES.

1. Introduction

Ecosystems provide human beings with a wide range of material goods and services to survive and develop including a great series of material services (such as agricultural products, timber, water and other raw materials) and immaterial services (such as water and air purification, climate regulation and aesthetic benefits) [1,2]. All these material and immaterial benefits people obtain from ecosystems are generally called ecosystem services (ES) [3]. ES are generally classified into provision, regulation, support and culture services [2]. In the past 50 years, ecosystems have been destroyed in rapid urbanization. A total of 60% of ES have degraded worldwide, especially in areas with considerable human influence and dramatic land use changes [2]. As the second largest economy in the world, China has demonstrated remarkable achievements in economic growth since 1979. However, this rapid growth has come at a price and the environment is quickly deteriorating because of pollution and irrational utilization [4,5]. The Organization for Economic Cooperation and Development (OECD) pointed out that the economy of China is close to that of the developed world but its environment is similar to that of the poorest countries [6].
The Beijing–Tianjin–Hebei (BTH) region, which is the largest plain in Northern China, has also experienced rapid population growth and economic development. However, human activities, such as deforestation, mountain reclamation and excessive resource exploitation, have exacerbated the environmental problems in the region [7,8], thereby causing the decay of ES. The BTH region has a long history of agricultural reclamation. Considerable high-quality farmland has been lost due to the urban sprawl in the region and food production is projected to decrease from 109.61 × 106 t in 2013 to 106.14 × 106–108.14 × 106 t in 2040 [9]. Furthermore, the entire BTH region is classified to be under absolute scarcity level by UN standards (2013) [10]. Groundwater level has declined continuously at a rate of 0.36 m/a in the past 50 years [11]. Meanwhile, the degradation and destruction of forests and grasses has led to massive soil erosion and desertification [12]. In the past few decades, the mountains Taihang and Yan-shan have been subject to serious soil erosion because of vegetation deterioration caused by overgrazing and overfarming [13]. Photosynthesis of vegetation can improve air quality and regulate climate. However, the urbanization and industrialization in the BTH region does not only reduce urban green space, thereby lessening the capacity of carbon sequestration but also produces considerable carbon emissions, thus leading to the heat island effect and air pollution [14,15].
With respect to the dilemma between economic development and environmental deterioration, many studies have concentrated on revealing the relationship among various ES and seeking an improved solution to achieve a mutually beneficial objective. Nelson et al. [3] modeled the changes in ecosystem services, biodiversity conservation, commodity production in a spatially explicit manner, found tradeoffs between them in the Willamette Basin, Oregon and made a deal with more effective, efficient and defensible natural resource decisions. Goldstein et al. [16] evaluated the environmental and financial implications of seven planning scenarios in Kamehameha Schools and found that tradeoffs existed between carbon storage and water quality as well as between environmental improvement and financial return, which helps to guide the local land-use decisions by involving tradeoffs between private and public interests. Relevant researches were also conducted in the BTH region. By developing the relationship between net primary productivity (NPP) and grain yield, Peng et al. [17] quantified and mapped the grain production of farmland and found that the main high-production areas were mainly located in the central plain of the BTH region and the eastern piedmont plain. The plain had low soil retention because of the low potential soil loss in farmland. Li et al. [11] evaluated regional water security in the BTH region via a freshwater ES flow model and concluded that freshwater ecosystem service provision was high in mountainous areas, especially in the Yan-shan and Taihang mountains and low in the Northern China Plain and urban areas. Nian et al. [18] analyzed the net capacity of ecosystems that provide carbon fixation and oxygen release (CFOR) in the BTH region using the “provision–reception” framework and the net CFOR exhibited a declining trend from the northwest to the southeast. To fulfill people’s demand for ES in the BTH region, Zhang et al. [19] proposed an evaluation framework of landscape ecological security pattern integrating ES assessment and landscape connectivity analysis with human ecological demand to identify the local ecological characteristics. Although these studies have provided good approaches and frameworks to quantify ES in the BTH region, no unified suggestions or environmental policies can be applied in the whole region because of the significant spatial heterogeneity within BTH.
In 2015, the outline of collaborative development of Beijing, Tianjin and Hebei Province was released by the Chinese government to strengthen the ecosystem management in the region. Formulating effective policies that target regional characteristics is a major challenge in improving environmental performance, especially at the sub region level of the BTH area. Therefore, a clear exploration is needed to deepen the correlations among ES and the spatiotemporal changes over time in different physical–geographical regions [7,20].
The scientific interests of this study focused on the following: (1) determining the correlation among different ES in the BTH region and their change over time, (2) identifying variance in these correlations in the different ecological–functional zones of the BTH region and (3) learning the suggestion of the sub regional variance in these correlations to local management. Four types of ES, namely, water yield (WY), soil retention (SR), carbon sequestration (CS) and grain production (GP), in 2000, 2005, 2010 and 2015, were quantified and mapped. The correlation among the different ES in various ecological–functional zones was analyzed by Spearman nonparametric tests to reveal the spatial relationship. Finally, on the basis of the results of the correlations among ES, we proposed sustainable development advice for the different zones.

2. Background of the Study

Changes in one ES can cause changes in others through three main relationships as follows. (1) The synergy describes a positive reaction of ES to a change in the driver [21]; (2) trade-off is an antagonistic situation that involves losing one quality of an ES to gain another [22]; and (3) compatibility means that a change in ES does not significantly influence others. ES relationships have been discussed previously. Several scholars found a synergy between supply service and biodiversity [3,23] and others found a trade-off correlation between provision and regulation services [24].
At present, common research methods of analyzing trade-off or synergy correlation mainly include map comparison, scenario analysis and trade-offs using optimized landscapes [25]. The first method maps each ES and analyzes their correlation cell by cell via spatial overlay [26]. As the most commonly used method, scenario analysis sets a series of situations and analyzes trade-off or synergy correlation among ES under various conditions [27]. The third method quantifies ES with biophysical or statistical models and then compares trade-off or synergy correlation on the basis of optimized land use patterns [28].
The trade-off or synergy correlation among ES is not uniform because of the varying geographical characteristics. This correlation depends on the ecosystem compositions, structures, conditions [3,7], complex interactions between physical (e.g., topography, geology) and biological (e.g., vegetation and micro-organism) factors and land use and management [8,29]. Jiang et al. [4] and Qin et al. [30] analyzed the trade-off and synergy between ES in the Three-river Headwater area and the Guanzhong–Tianshui economic region, respectively. The studies found a trade-off between water yield provision and carbon sequestration in the former region, which was in the hinterland of the Qinghai–Tibet Plateau and had a plateau continental climate and a strong synergy in the latter region, which was in the monsoon area of transition from semi-humid to semi-arid. REDD+ policies have attempted to achieve mutually beneficial solutions to enhance the carbon sequestration service and the conservation of tropical biodiversity; however, the results demonstrated an apparent trade-off between them [31]. Therefore, economic and ecological criteria should be regionally balanced in the application of land use policies because various regions respond differently to ES.

3. Study Area and Materials

3.1. Study Area

The BTH region (36°01′–42°37′ N, 113°04′–119°53′ E) is in Northern China and includes Beijing City, Tianjin City and Hebei Province. The region has the typical temperate semi-humid to semi-arid monsoon climate. Rainfall mainly occurs in summer and presents a unimodal distribution. Average annual rainfall is roughly 538 mm and presents a decreasing spatial pattern from eastern coastal to western inland. The northwest region has an undulating and mountainous topography with steep slopes and the southeast plain is characterized as an agricultural production base. The region has an elevation range of 0–2800 m and includes highland, mountain land, piedmont plain, low plain and littoral plain. On the basis of the physiographic conditions, the BTH region can be divided into four ecological–functional zones (Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, 2000), namely, deciduous broad-leaved forest ecoregion of the Yanshan–Taihang mountain (Zone I), agricultural ecoregion of Northern China Plain (Zone II), urban and suburban agricultural ecoregion of the Beijing–Tianjin–Tangshan region (Zone III) and typical steppe ecoregion of Central–Eastern Inner Mongolia Plateau (Zone IV) (Figure 1).

3.2. Materials

Land use datasets of the BTH region were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (http://www.resdc.cn), including farmland, forestland, grassland, water and wetland, settlement, desert and others. The boundaries of the ecological–functional zone were provided by the Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences (http://www.ecosystem.csdb.cn/). Raster-formatted soil data, that is, soil texture, soil types, organic matter content, reference soil depth, root-restricting layer depth and plant available water content (PAWC), were acquired from the Harmonized World Soil Database. Global digital elevation model data were downloaded from the United States National Aeronautics and Space Administration. Climate data were obtained from China National Meteorological Information Center. The boundaries of the BTH region were obtained from China Data Sharing Infrastructure of Earth System Science (2013). Vector data of watershed were required in water yield assessment; the entire study area was delineated into 122 watersheds using the Integrated Valuation of Ecosystem Services and Tradeoff model (InVEST model) [32]. All the data were uniformly resampled at a resolution of 1000 m. Average annual reference evapotranspiration was provided by the MOD16A3 product from the Numerical Terradynamic Simulation Group (http://www.ntsg.umt.edu/) but 2015 data was unavailable and thus replaced by 2014 data. Annual NPP was provided by the MYD17A3H Version 6 product from the Land Processes Distribution Active Archive Center (https://lpdaac.usgs.gov/).

4. Methodology

In this study, water yield and soil retention were evaluated by InVEST model, which was a model designed to map and quantify the distribution of ES through a set of standalone but linkable models and a tool for exploring how changes in ecosystems [32]. According to photosynthesis equation, NPP can be convert to the amount of carbon sequestration, therefore, NPP was used to represent the carbon sequestration, the capability of vegetation to absorb carbon dioxide and was calculated by the Carnegie–Ames–Stanford approach [33]. Grain production service was calculated by the gradation on agricultural land quality (GALQ) model, which is derived from the Agricultural Land Classification and Gradation Project by the Ministry of Land and Resources of China. Soil retention, carbon sequestration and grain productive capacity are presented at pixel scale, whereas water yield is shown at sub-watershed scale. To facilitate comparison, we indicated the values of soil retention, carbon sequestration and grain productive capacity at sub-watershed scales. Spearman nonparametric tests were performed to check for the correlations among the ES above to illustrate the regional heterogeneity in the four ecological–functional zones. The main workflow fluxes are presented in Figure 2.

4.1. Water Yield Assessment

In the InVEST model, water yield was calculated on the basis of the Budyko curve [34] and annual average precipitation as follows:
Y m = x Y x j     x m ,
Y x j = ( 1 A E T x j P x ) × P x ,
where Ym denotes the water yield in watershed m; Yxj and AETxj are the water yield and annual actual evapotranspiration on pixel x with land use type j, respectively; and Px denotes the annual precipitation on pixel x. In Equation (2), AETxj/Px approximates the Budyko curve estimated by Zhang et al. [35] as follows:
A E T x j P x = 1 + ω x R x j 1 + ω x R x j + ( 1 / R x j ) ,
where Rxj denotes the dimensionless Budyko dryness, which is the ratio of potential evapotranspiration to precipitation and ωx is the plant available water coefficient on pixel x.
The input data included annual precipitation, annual potential evapotranspiration (ET0), soil depth, PAWC, land use and land cover and root depth. ET0 was calculated by the FAO Penman–Monteith equation [36]. The root depth was referenced from Canadell et al. [37]. This empirical model used for Budyko theory has been tested at larger than pixel scales [38]; thus, the output summed or averaged at the sub-watershed scale has a higher credibility and is more interpretable than that at the pixel scale.

4.2. Soil Retention Assessment

Soil retention was calculated by the universal soil loss equation (USLE) [39], which computes the average annual amount of potential soil loss in a particular area on the basis of its land cover characteristics and topography. Soil retention was calculated as soil loss without vegetation cover and soil erosion control minus soil loss for the current land use and soil erosion management. Soil retention can be expressed mathematically as follows:
Δ A = A 0 A v = R × K × L × S × ( 1 C × P ) ,
where ΔA is the soil retention; A0 and Av represent the potential soil erosion without vegetation cover and the soil erosion under the current land cover and management conditions, respectively; R, K, L and S are the rainfall erosivity, soil erodibility, slope length and slope angle factors, respectively; and C and P are the current vegetation cover and erosion control factors, respectively.
R factor is calculated by a modified formula based on the data of annual and monthly precipitation [39,40]. Soil texture is a principal factor that affects K; the structure, organic matter and permeability contribute as well. We calculated factor K by the erosion/productivity impact calculator model [41]. The C and P factors were referenced by He et al. (2007), Bi et al. (2004) and Jiang et al. (2011) [42,43,44] according to the actual situation in the BTH region.

4.3. Grain Productive Capacity Assessment

On the basis of the relationship between soil and crop production, the GALQ model quantifies the grain productive capacity of agricultural land in China. Different from the real productivities which are not only affected by the quality factors of farmland but also restricted by the regional agricultural input (including fertilizer utilization, pest control and agricultural machinery) and agricultural technology management (including horticultural techniques and labor input), the grain productive capacity focuses on natural endowments of cropping system rather than agricultural mechanization levels. In this model, grain productive capacity is affected by natural endowments, including local climatic conditions (such as light, temperature and water) and soil quality [45]. To emphasize the natural attributes of agricultural land, the grain productive capacity calculated by the GALQ model refers to the highest yield of fine crop variety per unit area that might be obtained assuming all or part of the production factors are in the optimum state and take no account of agricultural input and agricultural technology management; therefore, the value is higher than the real grain yield. The cardinal principle of this model is building a regression equation between the physical quality of the agricultural land and the recorded highest grain yield in local history to estimate maximum grain productivity. Using a database of the experimental samples, this work assessed the grain productive capacity in the BTH region (Table 1).
Here, y denotes the grain productive capacity (per unit yield) and x denotes the physical quality of the agricultural land, including cultivated land, yard and facility farmland.

5. Results

5.1. Water Yield

The spatial heterogeneity of the water yield is visualized in Figure 3. High values of water yield were distributed in the eastern region, including Zone II, northeast of Zone III and east of Zone I. The low-value areas were in the northwest, including Zone IV, southwest of Zone III and northwest of Zone I. In general, local climate and land use affected the water yield. In comparison with the region with small forestland and settlement land, the western area of Zone I and the southeast of Zone III, which had forestland and settlement land, evaporated more water. Zone IV, which was part of the Inner Mongolia Plateau region, had a relatively high altitude and the low precipitation led to a low water yield.
The annual average water yield had a remarkable increase from 193.41 mm/km2 in 2000 to 240.35 mm/km2 in 2015 (Figure 4). The water yield in Zone II was higher than that in other regions in all years except 2005; that in Zone IV was lower than other regions in all years except 2010.

5.2. Sediment Retention

Sediment retention had high spatial distribution patterns in Zones I and IV, where the topography was mainly mountain and land coverage was mainly forest and grass. The low value was in Zones II and III, where plain was the main landform (Figure 5). The sediment erosion of the entire region increased from 990.07 ton/hm2 in 2000 to 1119.53 ton/hm2 in 2015 with a slight decline in 2015 (Figure 6). Zoning in the ecological–functional Zones I and IV had a similar increasing tendency as the entire region and peaked in 2010 (1714.97 ton/hm2 in Zone I and 1052.17 ton/hm2 in Zone IV). The soil retention in Zone II decreased and then increased and it reached a historical low in 2005. The only zone whose soil retention steadily increase from beginning to end was Zone III (from 534.19 ton/hm2 in 2000 to 698.93 ton/hm2 in 2015).

5.3. Carbon Sequestration

The carbon sequestration represented by NPP showed a similar spatial pattern as the sediment retention from 2000 to 2015. Zones I and IV, which had more forest and grass than the other zones, had high carbon sequestration values; meanwhile, low values were found in the other regions, which had less forest and grass than Zones I and IV (Figure 7). The highest annual average NPP in the entire BTH region was 324.54 g·C/m2 in 2000 and then decreased to 278.63 g·C/m2 in 2015 gradually (Figure 8). At the regional scale, the NPP in Zones I, II and III decreased by 11.64%, 18.67% and 30.88% since 2000, respectively but the tendency in Zone IV decreased initially and then increased.

5.4. Grain Productive Capacity

Significant differences in the spatial distribution of grain productive capacity are shown in Figure 9. Farmland was mainly distributed in Zones II and III; these farmland regions were of high quality because of the integration of fertile soil, sufficient sunlight and abundant rainfall. By contrast, a small amount of farmland was scattered in Zone I and the grain productive capacity was low. Most of the farmland distributed in Zone IV had medium–low yield.
The average grain productive capacity of the entire region had a slight decrease from 5.09 × 105 kg/km2 in 2000 to 4.98 × 105 kg/km2 in 2015 (Figure 10), which indicated that the quality of farmland in the whole region had degraded. Meanwhile, the grain productive capacity of Zones I, II and III all experienced a decreasing trend, especially that of Zone III, which dropped by 6.58%. By contrast, the production capability of Zone IV had a slight increase of roughly 0.67%.

5.5. Correlation among ES

Through Spearman correlation analysis, a positive correlation was found between the status quo of water yield and soil retention/grain production capability/carbon sequestration, soil retention and carbon sequestration in each year (Table 2) and the consistent correlation for the four point-in-times separated by five years can guarantee a stable relationship between each two ES; changes in water yield and soil retention in 2000–2015 (r2 = 0.469 **) were positively correlated as well. The changes in water yield in 2000–2015 were negatively correlated with those of grain production capability and carbon sequestration (r2 = −0.310 ** and −0.514 **, respectively) but an adverse relationship existed between soil retention and carbon sequestration (r2 = 0.426 **). It means that there is no uniformity between the correlation of status quo and that of changes for the same pair of ES.
In the entire study area, neither status quo nor the changes in grain production capability and soil retention/carbon sequestration were apparent in 2000–2015. However, this kind of nonobvious and unstable correlation changes at sub regional scale. Grain production capability had a positive correction with carbon sequestration in Zones II and III and a negative correction in Zones I and IV every four years (Table 3). In addition, the changes in grain production capability and carbon sequestration in 2000–2015 indicated a negative correction in Zones II and III and a positive one in Zones I and IV. A negative correlation was found between grain production capability and soil retention in Zones I and IV; no correlation existed in Zones II and III. Meanwhile, the changes in those within 2000–2015 were uncorrelated in Zones I, II and III and had a poor negative correction in Zone IV. The correlation of ES behaved differently according to the various terrains and vegetation. For instance, the grain production capability and soil retention were correlated negatively in the mountainous regions, such as Zone I and IV but not correlation in plain areas such as Zones II and III. As for the production capability and carbon sequestration, it also showed different relationships between them in mountainous and plain areas.

6. Discussion

6.1. Impact of Land Use and Land Cover Change on ES

Land use and land cover change is the main factor that affects ecosystems. Land use can significantly affect ES [46]. Figure 11, Figure 12 and Figure 13 reflect the trends of these four ES across different land types during the study period, which were as follows. (1) In farmland, water yield dropped from 2000 to 2005, reached the minimum (195.66 mm/km2) and then increased gradually. Soil retention slightly fluctuated from 422.28 ton/hm2 in 2000 to 496.02 ton/hm2 in 2010, whereas carbon sequestration decreased each year. (2) In forestland, water yield and soil retention peaked in 2010 (204.88 mm/km2 and 1549.17 ton/hm2, respectively) and then dropped slightly. Carbon sequestration decreased initially and tended to stabilize. (3) The trends of water yield and carbon sequestration for grassland were the same as those of forestland, whereas soil retention fluctuated up or down at roughly 1200 ton/hm2.
Coefficient of variation is used to reflect the trend of water yield in Figure 11. The trend decreased and then tended to stabilize within the range of 0.2–0.6. For the soil retention service (Figure 12), the coefficient of variation of farmland was considerably higher than that of forestland and grassland; therefore, in comparison with the soil retention of forestland and grassland, that of farmland significantly varied by region, being susceptible to natural regional factors. The coefficient of variation for carbon sequestration (Figure 13) remained within 0.2–0.35 and remained stable and this pattern illustrated the carbon sequestration of different land uses was not easily influenced by environmental deviation.
An index Ci was obtained to analyze the relationships among the different ES in the four ecological–functional zones and determine the intensity of land use change using land use transfer matrices as follows:
C i = A i A × 100 % ,
A = i = 1 n A i ,
where Ci is the percentage of land use change type i among all change types, Ai is the area of change type I and A is the total area of all land use change in 2000–2015. A high value of Ci demonstrates that this type of land use change plays a leading role in regional land use changes.
The results (Table 4) indicated the following. (1) The land use changes from farmland, grassland and forestland to settlement (SEL) were the main change types in Zone I and were accompanied with decreases in soil retention, water yield and carbon sequestration; (2) the conversion from farmland to settlement dominated the land use changes in Zone II, thereby leading to a decline in grain production, water yield and carbon sequestration; (3) all four types of ES in Zone III decreased because of the conversion from farmland, water and wetland (WAW) to settlement; and (4) the dominant conversion in Zone IV was slightly different from that of the other zones, which was from grass to farmland, thereby resulting in an increase in grain production and water yield but a decrease in soil retention and carbon sequestration.
In the BTH region, urban land expanded by 71% from 1990 to 2000; roughly 74% of the new urban land was converted from arable land [47]. This ratio increased to 81.72% from 2000 to 2015. The development of Zone III, which was the concentration district of the metropolis, mainly depended on secondary and tertiary industries. The urbanization progress directly led to a reduction in the average grain productive capacity in Zone III because of the encroachment of high-quality farmland (Table 3). Meanwhile, the development of Zone II relied mostly on primary industry; thus, this zone had more farmland than Zone III and had inconspicuous decreases in grain productive capacity.

6.2. Correlative Mechanisms among ES

ES are closely interlinked in a complex manner and characterized by a nonlinear relationship [48]. Bennett et al. [21] argued that ES are interrelated because of the (1) effects of common drivers or (2) interactions among ES per se. In the present research, water yield and sediment retention were correlated possibly because they shared common hydrological processes and drivers of precipitation and land use/cover change. For instance, the water retention capability of soil was correlated with water yield [32] because the soil loss might be caused by water flows. Specifically, the high water yield was consistently linked with high water retention capability in the soil. The lower the water flow, the lower the soil export. In addition, the growth of vegetation could improve the physical structure and then enhance the water retention capability of the soil. The high vegetation coverage limited soil moisture from achieving saturation and forming water flows [49,50] and consequently reduced soil erosion. This phenomenon explained the positive correlation between carbon sequestration and water yield/soil retention. In the current study, NPP was used to represent carbon sequestration; the higher the NPP, the larger the vegetation coverage [51]. Furthermore, these results also indicated that the areas with abundant vegetation except farmland consistently showed high values of water yield and soil retention (Table 2 and Table 3). The positive correlations between the differences in soil retention and water yield/carbon sequestration from 2000 to 2012 indicated the synergies between them on the Loess Plateau [4]. This study also showed negative relationships between water yield and carbon sequestration. In other words, a trade-off existed between water yield and carbon sequestration and this finding agreed with the studies in the Loess Plateau [48,52]. The Loess Plateau and the BTH region are semi-arid and semi-humid regions. With an increase in vegetation coverage, more rainwater retained by leaves evaporate and more rainwater are absorbed by vegetation on the surface and then evaporate. Water yield also decreases.
The study conducted by Li et al. [53] in 2010 on the BTH region showed a weak negative correlation between gain production and carbon sequestration/soil retention. The correlation analysis conducted by the present work in the four ecological–functional zones yielded the following results. (1) Production capability and soil retention were negatively correlated in Zones I and IV but not correlation in Zones II and III. The grain production in the mountain area easily caused soil erosion; thus, the grain production in Zones I and IV, where the terrain was predominantly mountain, was not conducive to soil retention. However, the grain production caused slight soil erosion in Zones II and III because of the flat terrain. (2) Grain production and soil retention had a weak trade-off in Zone IV and no correlation in Zones I, II and III. (3) In Zones I and IV, the forest and grass had a strong capability for carbon sequestration and a low capability for grain production. In Zones II and III, farmland had a high capability for carbon sequestration and grain production. Grain production and carbon sequestration had a negative correlation in Zones I and IV and a positive correlation in Zones II and III. (4) Grain production and carbon sequestration exhibited synergy in Zones I and IV but a trade-off in Zones II and III. The synergy in Zones I and IV agreed with those found by the study conducted on the Loess Plateau in northern part of middle China [48,52]. The trade-offs in Zones II and III were consistent with those shown in the study conducted on the Sanjiang Plain in northeastern China for the period of 1992–2012 [54].

6.3. Implications for Regional Sustainable Development

In 2014, the Chinese President Xi proposed a coordinated development of the BTH region. The generation of an eco-environmental protection strategy should eliminate the restrictions imposed by the administrative boundaries and take the ecological function in its entirety for the further regional sustainable development. Recently, the Chinese national and local governments implemented a series of measures on territory development plans and ecological environmental protection. Although certain evidence indicates that these measures have certain effects, considerable improvement is still needed for strengthening the implementation of environmental policies [55]. These policies have been applied in different ecological–functional zones and the results from other studies are listed in Table 5.
On the basis of the characteristics of land use conversion and the ecological environment problems in the different regions from 2000 to 2015, we combined the spatial correlation, trade-offs and synergies among ES and then proposed suggestions for sustainable regional development.
Zone I comprised the deciduous broad-leaved forest ecoregion of the Yanshan-Taihang Mountain, which is the eco-security shield in the northwest of the BTH region. With high mountains, steep slopes and mountainous forestland, the ecosystem in Zone I has degraded seriously over the last few decades, exhibiting ecosystem homogenization, water conservation capacity decrease and soil erosion aggravation. The area with low vegetation coverage shows land desertification and is a high-incidence area for debris flow and geological collapse. Our results demonstrated a negative spatial relationship between grain production and soil retention in Zone I. The grain production in this zone is not effective for soil retention. By contrast, the forestland and grassland can provide high-quality soil retention service. Therefore, the government must continue to implement the Sloping Land Conversion Program (SLCP) and strictly prohibit agricultural reclamation or deforestation on the high-slope areas to compensate for the encroachment in plains via the Policy of Dynamic Equilibrium of the total Cultivated Land (PDDCL). Moreover, certain engineering measures, such as the slope protection project and torrent control works, must be further improved to alleviate water and soil loss. For places with low vegetation and high risk of geologic disasters, livestock grazing and fuel gathering should be stopped to facilitate afforestation.
Zone II covered the agricultural ecoregion of the Northern China Plain and was in the southeastern plain of the BTH region. This area has a flat topography with abundant concentrated and continuous farmland and fertile soil, thereby contributing to its capability to offer grain production and carbon sequestration service. To prevent the encroachment of high-quality farmland, the Protection of Basic Agricultural Land (PBAL) should be implemented strictly and regional development must utilize land use intensively. Newly added settlements should occupy the idle land or the low-quality and scattered farmland. This region should focus on developing modern agriculture for which the High-standard Primary Farmland Construction Project (HSPFCP) should be generalized to proceed with rural land reclamation and readjustment. Our result indicated that farmland with high grain production decreases water yield. The farmland in Zone II should fallow moderately to recharge the groundwater.
Zone III comprised an urban and suburban agricultural ecoregion that included two sub-ecological regions of central urban zones and suburban agricultural ecological zones. This zone covered super urban agglomeration areas, including Beijing and Tianjin. The urban environment is deteriorating and appears as an urban heat island with degraded air quality [70]. Our results demonstrated that the farmland in this region had obviously decreased. The farmland had strong capacities of water yield and carbon sequestration and grain production in Zone III had a positive spatial correlation with carbon sequestration and water yield. As the main land cover type in the BTH region, the farmland not only improves the services of grain production, water yield and carbon sequestration but can also serve as an urban green landscape. Therefore, the laws of farmland protection and regulations should be made consistent, transparent and nondiscriminatory to prevent settlement encroachment. This work recommends that the cultivation of farmland be divided into three circle-layers, namely, circles of urban landscape agriculture, eco-leisure agriculture and high-efficiency facility agriculture.
Zone IV was a typical steppe ecoregion of the Central–Eastern Inner Mongolia Plateau that was situated mainly in the northwest of the BTH region. The grassland had been seriously degraded because of unreasonable human exploitation, exhibiting a decrease in grass cover, land desertification and depletion of soil nutrients. The results of this study showed a negative spatial relationship and weak trade-off between grain production and soil retention in Zone IV. To prevent the further deterioration of grassland, the government should control the extent of agricultural reclamation by prohibiting excessive cultivation and implement conservation tillage to prevent soil wind erosion and desertification. In addition, pasture management needs to ban grazing and implement the strategy of rotation grazing.

6.4. Limitations

Uncertainties remain in the ES assessment that is based on models and data in this study. First, various ES are generated in different spatial scales. In our research, the spatial resolution of 1 km was helpful in learning the dynamic change in the research area at the regional scale; however, this system inevitably concealed the land use changes in areas under 1 km2. Second, soil retention had not been calibrated by observation data due to limited available data, although the results were consistent with those obtained by Li et al. (2014) and Xu et al. (2011) [53,71]. Third, only the grain productive capacity of farmland was considered in the assessment of food provision service. Grassland provides fibrous feed for livestock and orchards offer fruits for people, thereby eventually providing food for human survival. In further research, the accounting scope should be further expanded to cover additional essential services. Fourth, our study analyzed the interaction of ES in four ecological function zones but spatial heterogeneity remains in each ecological–functional zone. Detailed decisions will be implemented on the basis of the natural environment and socioeconomic status within each region in the future. Finally, the spatial distribution of human habitation was mismatched with that of the ecosystem, thereby resulting in a mismatch between the demand for and supply of ES. The current decisions provided remain based on the ES supply. If possible, a policy of regulating ES for human needs that is based on human settlement distribution will be formulated by the researchers.

7. Conclusions

In this study, four ES were quantified and mapped explicitly spatially with biophysical models. The spatial correction, trade-offs and synergies among multiple ES were also analyzed in four ecological–functional regions to help determine mechanisms of ES as bases for effective decisions for sustainable regional development. The main conclusions are as follows.
  • The spatial distribution of water yield was relatively homogeneous but slightly high in Zones I, II and III and all increased generally each year. Sediment retention was high in Zones I and IV and low in Zones II and III. These values all gradually increased every year. The distribution of carbon sequestration was similar to that of soil retention but a downtrend occurred in the entire region rather than in Zone IV. High-quality farmland with high grain production was mainly distributed in Zones II and III and those in Zone III indicated a significant decline.
  • Spearman correlation analysis indicated a positive spatial relationship between the status quo of water yield and soil retention/grain production capability/carbon sequestration, soil retention and carbon sequestration in the whole region. Furthermore, a positive relationship was observed between grain production capability and carbon sequestration in Zones II and III but a negative one was observed in Zones I and IV. A negative correlation was observed between grain production capability and soil retention in Zones I and IV. Regarding the relationship between trade-off and synergy, our results indicated synergies between soil retention and water yield/carbon and trade-offs between water yield and carbon sequestration/grain production in the whole region. In addition, grain production and carbon sequestration had synergy in Zones I and IV but trade-offs in Zone II and III. A weak trade-off was observed between grain production and soil retention only in Zone IV.
  • The results of the transfer matrix demonstrated that encroachment on farmland occurred in the BTH region in 2000–2015, except for the grassland occupation by farmland in Zone IV. The government is suggested to implement sequentially the SLCP and prohibit reclamation and deforestation in Zone I. The authors emphasize the protection of high-quality farmland in Zone II to render it the “rice bag” of the BTH region. Furthermore, the farmland should be made multi-functional and its encroachment in Zone III should be prevented. Finally, the farmland should be returned for grass and grazing should be controlled in Zone IV to restore grassland ecology.

Acknowledgments

The authors gratefully acknowledge the financial support by the National Natural Science Foundation of China (41501087, 41301614).

Author Contributions

Yang Gao and Zhen Xie designed the paper. Zhen Xie collected the data and wrote the paper. Chao Li revised the paper. Jian Zhou and Tianzhu Zhang contributed to the modification. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, 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]
  2. Millennium Ecosystem Assessment. Ecosystems and Human Well-Being: Current State and Trends: Synthesis; Island Press: Washington, DC, USA, 2005. [Google Scholar]
  3. Nelson, E.; Mendoza, G.; Regetz, J.; Polasky, S.; Tallis, H.; Cameron, D.R.; Chan, K.M.; Daily, G.C.; Goldstein, J.; Kareiva, P.; et al. Modeling multiple ecosystem services, biodiversity conservation, commodity production and tradeoffs at landscape scales. Front. Ecol. Environ. 2009, 7, 4–11. [Google Scholar] [CrossRef]
  4. Jiang, C.; Wang, F.; Zhang, H.; Dong, X. Quantifying changes in multiple ecosystem services during 2000–2012 on the Loess Plateau, China, as a result of climate variability and ecological restoration. Ecol. Eng. 2016, 97, 258–271. [Google Scholar] [CrossRef]
  5. Song, W.; Liu, M. Assessment of decoupling between rural settlement area and rural population in China. Land Use Policy 2014, 39, 331–341. [Google Scholar] [CrossRef]
  6. Organization for Economic Co-Operation and Development (OECD). Environmental Performance Reviews: China; OECD Press: Paris, France, 2007. [Google Scholar]
  7. Peng, J.; Chen, X.; Liu, Y.; Lü, H.; Hu, X. Spatial identification of multifunctional landscapes and associated influencing factors in the Beijing-Tianjin-Hebei region, China. Appl. Geogr. 2016, 74, 170–181. [Google Scholar] [CrossRef]
  8. Xie, H.; He, Y.; Xie, X. Exploring the factors influencing ecological land change for China’s Beijing–Tianjin–Hebei region using big data. J. Clean. Prod. 2017, 142, 677–687. [Google Scholar] [CrossRef]
  9. Zhang, D.; Huang, Q.; He, C.; Wu, J. Impacts of urban expansion on ecosystem services in the Beijing-Tianjin-Hebei urban agglomeration, China A scenario analysis based on the Shared Socioeconomic Pathways. Resour. Conserv. Recycl. 2017, 125, 115–130. [Google Scholar] [CrossRef]
  10. Gao, Y.; Feng, Z.; Li, Y.; Li, S. Freshwater ecosystem service footprint model: A model to evaluate regional freshwater sustainable development—A case study in Beijing–Tianjin–Hebei, China. Ecol. Indic. 2014, 39, 1–9. [Google Scholar] [CrossRef]
  11. Li, X.; Zhang, Y.; Zhou, P.; Ye, S. Dynamics and evolution of the groundwater in Beijing-Tianjin-Hebei Plain in a long time scale. J. Arid Land Resour. Environ. 2017, 31, 164–170. (In Chinese) [Google Scholar]
  12. Xu, D.; Li, C.; Song, X.; Ren, H. The dynamics of desertification in the farming-pastoral region of North China over the past 10 years and their relationship to climate change and human activity. Catena 2014, 123, 11–22. [Google Scholar] [CrossRef]
  13. Zhu, Z. Concept, cause and control of desertification in China. Quat. Sci. 1998, 18, 145–155. [Google Scholar]
  14. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, X.; Du, J.; Huang, T.; Zhang, L.; Gao, H.; Zhao, Y.; Ma, J. Atmospheric removal of PM2.5 by man-made Three Northern Regions Shelter Forest in Northern China estimated using satellite retrieved PM2.5 concentration. Sci. Total Environ. 2017, 593, 713–721. [Google Scholar] [CrossRef] [PubMed]
  16. Goldstein, J.H.; Caldarone, G.; Duarte, T.K.; Ennaanay, D.; Hannahs, N.; Mendoza, G.; Polasky, S.; Wolny, S.; Daily, G.C. Integrating ecosystem-service tradeoffs into land-use decisions. Proc. Natl. Acad. Sci. USA 2012, 109, 7565–7570. [Google Scholar] [CrossRef] [PubMed]
  17. Peng, J.; Shen, H.; Wu, W.; Wang, Y. Net primary productivity (NPP) dynamics and associated urbanization driving forces in metropolitan areas: A case study in Beijing City, China. Landsc. Ecol. 2016, 31, 1077–1092. [Google Scholar] [CrossRef]
  18. Nian, W.; Chen, Y.; Gao, J.; Ma, X.; Yao, X. Relationship Between Provision and Reception of Ecological Service of Carbon Fixation and Oxygen Release in Beijing–Tianjin–Hebei Region. J. Ecol. Rural Environ. 2017, 33, 783–791. (In Chinese) [Google Scholar]
  19. Zhang, L.; Peng, J.; Liu, Y.; Wu, J. Coupling ecosystem services supply and human ecological demand to identify landscape ecological security pattern: A case study in Beijing–Tianjin–Hebei region, China. Urban Ecosyst. 2017, 20, 701–714. [Google Scholar] [CrossRef]
  20. Portman, M.E. Ecosystem services in practice: Challenges to real world implementation of ecosystem services across multiple landscapes—A critical review. Appl. Geogr. 2013, 45, 185–192. [Google Scholar] [CrossRef]
  21. Bennett, E.M.; Peterson, G.D.; Gordon, L.J. Understanding relationships among multiple ecosystem services. Ecol. Lett. 2009, 12, 1394–1404. [Google Scholar] [CrossRef] [PubMed]
  22. Cord, A.F.; Bartkowski, B.; Beckmann, M.; Dittrich, A.; Hermans-Neumann, K.; Kaim, A.; Lienhoop, N.; Locher-Krause, K.; Priess, J.; Schröter-Schlaack, C.; et al. Towards systematic analyses of ecosystem service trade-offs and synergies: Main concepts, methods and the road ahead. Ecosyst. Serv. 2017, 28, 264–272. [Google Scholar] [CrossRef]
  23. Perring, M.P.; Standish, R.J.; Hulvey, K.B.; Lach, L.; Morald, T.K.; Parsons, R.; Didham, R.K.; Hobbs, R.J. The Ridgefield Multiple Ecosystem Services Experiment: Can restoration of former agricultural land achieve multiple outcomes? Agric. Ecosyst. Environ. 2012, 163, 14–27. [Google Scholar] [CrossRef]
  24. Su, C.; Fu, B.; He, C.; Lü, Y. Variation of ecosystem services and human activities: A case study in the Yanhe Watershed of China. Acta Oecol. 2012, 44, 46–57. [Google Scholar] [CrossRef]
  25. Lautenbach, S.; Volk, M.; Gruber, B.; Dormann, C.; Strauch, M.; Seppelt, R. Quantifying ecosystem service trade-offs. In Proceedings of the International Environmental Modelling and Software Society (iEMSs) 2010 International Congress on Environmental Modelling and Software Modelling for Environment’s Sake, Ottawa, ON, Canada, 5–8 July 2010. [Google Scholar]
  26. Egoh, B.; Reyers, B.; Rouget, M.; Bode, M.; Richardson, D.M. Spatial congruence between biodiversity and ecosystem services in South Africa. Biol. Conserv. 2009, 142, 553–562. [Google Scholar] [CrossRef]
  27. Butler, J.R.A.; Wong, G.Y.; Metcalfe, D.J.; Honzák, M.; Pert, P.L.; Rao, N.; van Grieken, M.E.; Lawson, T.; Bruce, C.; Kroon, F.J.; et al. An analysis of trade-offs between multiple ecosystem services and stakeholders linked to land use and water quality management in the Great Barrier Reef, Australia. Agric. Ecosyst. Environ. 2013, 180, 176–191. [Google Scholar] [CrossRef]
  28. White, C.; Halpern, B.S.; Kappel, C.V. Ecosystem service tradeoff analysis reveals the value of marine spatial planning for multiple ocean uses. Proc. Natl. Acad. Sci. USA 2012, 109, 4696–4701. [Google Scholar] [CrossRef] [PubMed]
  29. Guswa, A.J.; Brauman, K.A.; Brown, C.; Pert, P.L.; Rao, N.; van Grieken, M.E.; Lawson, T.; Bruce, C.; Kroon, F.J.; Brodie, J.E. Ecosystem services: Challenges and opportunities for hydrologic modeling to support decision making. Water Resour. Res. 2014, 50, 4535–4544. [Google Scholar] [CrossRef]
  30. Qin, K.; Li, J.; Yang, X. Trade-Off and Synergy among Ecosystem Services in the Guanzhong–Tianshui Economic Region of China. Int. J. Environ. Res. Public Health 2015, 12, 14094–14113. [Google Scholar] [CrossRef] [PubMed]
  31. Phelps, J.; Friess, D.A.; Webb, E.L. Win–win REDD+ approaches belie carbon–biodiversity trade-offs. Biol. Conserv. 2012, 154, 53–60. [Google Scholar] [CrossRef]
  32. Tallis, H.T.; Ricketts, T.; Guerry, A.D.; Wood, S.A.; Sharp, R.; Nelson, E.; Ennaanay, D.; Wolny, S.; Olwero, N.; Vigerstol, K.; et al. InVEST 2.5.3 User’s Guide; The Natural Capital Project: Stanford, CA, USA, 2013. [Google Scholar]
  33. Potter, C.S.; Randerson, J.T.; Field, C.B.; Matson, P.A.; Vitousek, P.M.; Mooney, H.A.; Klooster, S.A. Terrestrial ecosystem production: A process model based on global satellite and surface data. Glob. Biogeochem. Cycles 1993, 7, 811–841. [Google Scholar] [CrossRef]
  34. Budyko, M.I. Climate and Life; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  35. Zhang, L.; Dawes, W.R.; Walker, G.R. Response of mean annual evapotranspiration to vegetation changes at catchment scale. Water Resour. Res. 2001, 37, 701–708. [Google Scholar] [CrossRef]
  36. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 1998. [Google Scholar]
  37. Canadell, J.; Jackson, R.B.; Ehleringer, J.R.; Mooney, H.A.; Sala, O.E.; Schulze, E. Maximum Rooting Depth of Vegetation Types at the Global Scale. Oecologia 1996, 108, 583–595. [Google Scholar] [CrossRef] [PubMed]
  38. Redhead, J.W.; Stratford, C.; Sharps, K.; Jones, L.; Ziv, G.; Clarke, D.; Oliver, T.H.; Bullock, J.M. Empirical validation of the InVEST water yield ecosystem service model at a national scale. Sci. Total Environ. 2016, 569, 1418–1426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses: A Guide to Conservation Planning; Agriculture Handbook, No. 537; Government Printing Office: Washington, DC, USA, 1978.
  40. Arnoldus, H.M.J. An approximation of the rainfall factor in the Universal Soil Loss Equation. In Assessment of Erosion; DeBoodt, M., Gabriels, D., Eds.; John Wiley: Chichester, UK, 1980; pp. 127–132. [Google Scholar]
  41. Williams, J.R.; Arnold, J.G. A system of erosion—Sediment yield models. Soil Technol. 1997, 11, 43–55. [Google Scholar] [CrossRef]
  42. He, R.; Wang, B.; Zhu, G.; Qi, S. Prediction of Soil Erosionina Water Catchment of Miyun Reservoir. J. Agro Environ. Sci. 2007, 26, 579–582. (In Chinese) [Google Scholar]
  43. Bi, X.; Duan, S.; Li, Y.; Liu, B.; Fu, S.; Ye, Z.; Yuan, A.; Lu, B. Study on soil loss equation in Beijing. Sci. Soil Water Conserv. 2004, 4, 6–13. (In Chinese) [Google Scholar]
  44. Jiang, Q.; Xie, Y.; Zhang, Y.; Zhang, H.; Hao, X. Soil erosion in a small watershed in water source areas of Beijing and Tianjin: Spatial simulation. Chin. J. Ecol. 2011, 30, 1703–1711. (In Chinese) [Google Scholar]
  45. Yun, W.; Wang, H.; Wang, G.; Zhang, L. Research of Throughput Calculation Based on Agricultural Land Classification and Agriculture Statistics. China Land Sci. 2007, 21, 32–37. (In Chinese) [Google Scholar]
  46. Costanza, R.; Mitsch, W.J.; Day, J.W. A new vision for New Orleans and the Mississippi delta: Applying ecological economics and ecological engineering. Front. Ecol. Environ. 2006, 4, 465–472. [Google Scholar] [CrossRef]
  47. Tan, M.; Li, X.; Xie, H.; Lu, C. Urban land expansion and arable land loss in China—A case study of Beijing–Tianjin–Hebei region. Land Use Policy 2005, 22, 187–196. [Google Scholar] [CrossRef]
  48. Su, C.; Fu, B. Evolution of ecosystem services in the Chinese Loess Plateau under climatic and land use changes. Glob. Planet. Chang. 2013, 101, 119–128. [Google Scholar] [CrossRef]
  49. Castillo, V.M.; Martinez-Mena, M.; Albaladejo, J. Runoff and soil loss response to vegetation removal in a semiarid environment. Soil Sci. Soc. Am. J. 1997, 61, 1116–1121. [Google Scholar] [CrossRef]
  50. Morgan, R.P.C. Soil Erosion and Conservation; Longman: London, UK, 1995. [Google Scholar]
  51. Running, S.W.; Nemani, R.R.; Heinsch, F.A.; Zhao, M.; Reeves, M.; Hashimoto, H. A continuous satellite-derived measure of global terrestrial primary production. Bioscience 2004, 54, 547–560. [Google Scholar] [CrossRef]
  52. Lv, Y.; Fu, B.; Feng, X.; Zeng, Y.; Liu, Y.; Chang, R.; Sun, G.; Wu, B. A Policy-Driven Large Scale Ecological Restoration: Quantifying Ecosystem Services Changes in the Loess Plateau of China. PLoS ONE 2012, 7, e31782. [Google Scholar] [CrossRef]
  53. Li, S. The Geography of Ecosystem Service; Science Press: Beijing, China, 2014. [Google Scholar]
  54. Wang, Z.; Mao, D.; Li, L.; Dong, Z.; Miao, Z.; Ren, C.; Song, C. Quantifying changes in multiple ecosystem services during 1992–2012 in the Sanjiang Plain of China. Sci. Total Environ. 2015, 514, 119–130. [Google Scholar] [CrossRef] [PubMed]
  55. Li, W.; Michalak, K. Environmental Compliance and Enforcement in China; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 379–391. [Google Scholar]
  56. Sun, P.; Xu, Y.; Wang, S. Terrain gradient effect analysis of land use change in poverty area around Beijing and Tianjin. Trans. Chin. Soc. Agric. Eng. 2014, 30, 277–288. (In Chinese) [Google Scholar]
  57. Zhao, F.; Yang, G.; Han, X.; Feng, Y.; Ren, G. Stratification of carbon fractions and carbon management index in deep soil affected by the Grain-to-Green program in China. PLoS ONE 2014, 9, e99657. [Google Scholar] [CrossRef] [PubMed]
  58. Yang, X.; Xu, B.; Jin, Y.; Qin, Z.; Ma, H.; Li, J.; Zhao, F.; Chen, S.; Zhu, X. Remote sensing monitoring of grassland vegetation growth in the Beijing–Tianjin sandstorm source project area from 2000 to 2010. Ecol. Indic. 2015, 51, 244–251. [Google Scholar] [CrossRef]
  59. Peng, D.; Wu, C.; Zhang, B.; Huete, A.; Zhang, X.; Sun, R.; Lei, L.; Huang, W.; Liu, L.; Liu, X.; et al. The Influences of Drought and Land-Cover Conversion on Inter-Annual Variation of NPP in the Three-North Shelterbelt Program Zone of China Based on MODIS Data. PLoS ONE 2016, 11, e0158173. [Google Scholar] [CrossRef] [PubMed]
  60. Cao, S. Why large-scale afforestation efforts in China have failed to solve the desertification problem. Environ. Sci. Technol. 2008, 42, 1828–1831. [Google Scholar] [CrossRef]
  61. Wang, X.; Zhang, X.; Hasi, E.; Dong, C. Has the Three Norths Forest Shelterbelt Program solved the desertification and dust storm problems in arid and semiarid China? J. Arid Environ. 2010, 74, 13–22. [Google Scholar] [CrossRef]
  62. Zhang, L.; Fu, B.; Lv, Y.; Dong, Z.; Li, Y.; Zeng, Y.; Wu, B. The using of composite indicators to assess the conservational effectiveness of ecosystem services in China. Acta Geogr. Sin. 2016, 71, 768–780. (In Chinese) [Google Scholar]
  63. Song, W.; Deng, X. Land-use/land-cover change and ecosystem service provision in China. Sci. Total Environ. 2017, 576, 705. [Google Scholar] [CrossRef] [PubMed]
  64. Jiang, B.; Wong, C.; Chen, Y.; Cui, Y.; Ouyang, Z. Advancing Wetland Policies Using Ecosystem Services–China’s Way Out. Wetlands 2015, 35, 983–995. [Google Scholar] [CrossRef]
  65. Jiang, M.; Tian, S.; Zheng, Z.; Zhan, Q.; He, Y. Human Activity Influences on Vegetation Cover Changes in Beijing, China, from 2000 to 2015. Remote Sens. 2017, 9, 271. [Google Scholar] [CrossRef]
  66. Feng, T.; Zhang, F.; Nie, X.; Wang, H. Spatial-temporal change characteristics of cropland, vegetable land and forest land in Beijing plain region. Trans. Chin. Soc. Agric. Eng. 2017, 33, 257–264. (In Chinese) [Google Scholar]
  67. Sun, X.; Wang, T.; Wu, J.; Ge, J. Change trend of vegetation cover in Beijing metropolitan region before and after the 2008 Olympics. Chin. J. Appl. Ecol. 2012, 23, 3133–3140. (In Chinese) [Google Scholar]
  68. Zhao, H.; Zhang, F.; Jiang, G.; Xu, Y.; Xie, Z. Analysis of the Trouble in Cultivated Land and Prime Farmland Protection in Beijing Suburbs Based on Peasents’ Household Surveys. China Land Sci. 2008, 22, 28–33. (In Chinese) [Google Scholar]
  69. Xiao, L.; Tian, G.; Qiao, Z. An Agent-based Approach for Urban Encroachment on Cropland Dynamic Model and Simulation. J. Nat. Resour. 2014, 29, 516–527. (In Chinese) [Google Scholar]
  70. Li, M.; Wang, T.; Xie, M.; Zhuang, B.; Li, S.; Han, Y.; Cheng, N. Modeling of urban heat island and its impacts on thermal circulations in the Beijing–Tianjin–Hebei region, China. Theor. Appl. Climatol. 2017, 128, 999–1013. [Google Scholar] [CrossRef]
  71. Xu, X.; Li, X.; Han, N. The Evaluation of the Ecological Service of Soil Conservation Based on Multisource Remote Sensing Data: A Case Study of Four Areas in Northern Hebei. Remote Sens. Land Resour. 2011, 19, 123–129. (In Chinese) [Google Scholar]
Figure 1. Ecological–functional zones of the BTH region.
Figure 1. Ecological–functional zones of the BTH region.
Sustainability 10 00006 g001
Figure 2. Workflow diagram of this research.
Figure 2. Workflow diagram of this research.
Sustainability 10 00006 g002
Figure 3. Water yield in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Figure 3. Water yield in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Sustainability 10 00006 g003
Figure 4. Average water yield in ecological–functional zones in 2000–2015.
Figure 4. Average water yield in ecological–functional zones in 2000–2015.
Sustainability 10 00006 g004
Figure 5. Soil retention service in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Figure 5. Soil retention service in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Sustainability 10 00006 g005
Figure 6. Average soil retention in ecological–functional zones in 2000–2015.
Figure 6. Average soil retention in ecological–functional zones in 2000–2015.
Sustainability 10 00006 g006
Figure 7. NPP in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Figure 7. NPP in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Sustainability 10 00006 g007
Figure 8. Annual average NPP in ecological–functional zones in 2000–2015.
Figure 8. Annual average NPP in ecological–functional zones in 2000–2015.
Sustainability 10 00006 g008
Figure 9. Grain productive capacity in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Figure 9. Grain productive capacity in (a) 2000, (b) 2005, (c) 2010 and (d) 2015.
Sustainability 10 00006 g009
Figure 10. Average grain productive capacity in ecological–functional zones in 2000–2015.
Figure 10. Average grain productive capacity in ecological–functional zones in 2000–2015.
Sustainability 10 00006 g010
Figure 11. Mean value and variability coefficient of water yield in farmland (FAL), forestland (FOL) and grassland (GRL) in 2000–2015.
Figure 11. Mean value and variability coefficient of water yield in farmland (FAL), forestland (FOL) and grassland (GRL) in 2000–2015.
Sustainability 10 00006 g011
Figure 12. Mean value and variability coefficient of soil retention in FAL, FOL and GRL in 2000–2015.
Figure 12. Mean value and variability coefficient of soil retention in FAL, FOL and GRL in 2000–2015.
Sustainability 10 00006 g012
Figure 13. Mean value and variability coefficient of carbon sequestration in FAL, FOL and GRL in 2000–2015.
Figure 13. Mean value and variability coefficient of carbon sequestration in FAL, FOL and GRL in 2000–2015.
Sustainability 10 00006 g013
Table 1. Regression Equations for Grain Productive Capacity with Physical Quality of Agricultural Land in Beijing, Tianjin and Hebei.
Table 1. Regression Equations for Grain Productive Capacity with Physical Quality of Agricultural Land in Beijing, Tianjin and Hebei.
RegionRegression EquationsR-Squared
Beijingy = 0.1576x + 332.440.5142
Tianjiny = 0.3076x − 101.230.7803
Hebeiy = 0.302x − 34.460.873
Table 2. Correlations of Four ESs in the BTH region.
Table 2. Correlations of Four ESs in the BTH region.
Water Yield Grain Production Capability
20002005201020152000–201520002005201020152000–2015
Soil retention20000.312 **----Water yield0.811 **----
2005-0.244 **----0.701 **---
2010--0.297 **----0.634 **--
2015---0.274 **----0.800 **-
Change during 2000–2015----0.469 **----−0.310 **
Grain Production capability Carbon Sequestration
20002005201020152000–201520002005201020152000–2015
Soil retention2000−0.144----Water yield0.594 **----
2005-−0.183 *----0.790 **---
2010--−0.178 *----0.830 **--
2015---−0.16----0.678 **-
2000–2015----−0.018----−0.514 **
Carbon Sequestration Grain Production Capability
20002005201020152000–201520002005201020152000–2015
Soil retention20000.606 **----Carbon sequestration0.021----
2005-0.674 **----0.180 *---
2010--0.677 **----0.036--
2015---0.645 **----0.187 *-
2000–2015----0.426 **----0.065
* Correlation is significant at the 0.05 level (two-tailed). ** Correlation is significant at the 0.01 level (two-tailed).
Table 3. Correlations of Three ESs in Ecological–Functional Zones.
Table 3. Correlations of Three ESs in Ecological–Functional Zones.
Ecosystem ServicesIII
20002005201020152000–201520002005201020152000–2015
GP–CS−0.206 **−0.266 **−0.249 **−0.246 **0.128 **0.191 **0.188 **0.196 **0.189 **−0.023 **
GP–SR−0.279 **−0.288 **−0.290 **−0.286 **0.004−0.008−0.008−0.009−0.0080.006
Ecosystem ServicesIIIIV
20002005201020152000–201520002005201020152000–2015
GP–CS0.347 **0.370 **0.365 **0.331 **−0.017 **−0.213 **−0.288 **−0.235 **−0.328 **0.154 **
GP–SR−0.007−0.005−0.004−0.0010.007−0.186 **−0.192 **−0.190 **−0.194 **−0.016 *
* Correlation is significant at the 0.05 level (two-tailed). ** Correlation is significant at the 0.01 level (two-tailed).
Table 4. Main Land Type and Ecosystem Service Changes in Ecological–Functional Zones.
Table 4. Main Land Type and Ecosystem Service Changes in Ecological–Functional Zones.
Ecological–Functional ZonesIDMain Land Type ChangesArea/km2Ci/%Ecosystem Service Changes
+
I1FAL–SEL30943.34 GP, WY, CS
2GRL–SEL11115.57 WY, SR, CS
3FOL–SEL7710.80 WY, SR, CS
Accumulation above42069.71
Other changes29330.29
II1FAL–SEL70383.29 GP, WY, CS
Accumulation above70383.29
Other changes14116.71
III1FAL–SEL125362.28 GP, WY, CS
2WAW–SEL24111.98 WY, SR, CS
Accumulation above125374.25
Other changes75925.75
IV1GRL–FAL6849.28GP, WYSR, CS
2FAL–SEL1913.77 GP, WY, CS
Accumulation above8763.04
Other changes5136.96
Table 5. Territory Development Plans and Ecological Environmental Protection Policies in the BTH region.
Table 5. Territory Development Plans and Ecological Environmental Protection Policies in the BTH region.
Ecological Function ZoneRepresentative CitiesDominant PolicyThe Content of PolicyImpacts for Land Use TypesResults and Related Research
I/IVBeijing, Chengde, Zhangjiakou and Baoding (mountain part)The Sloping Land Conversion Program (2002–2022) (SLCP)To protect and improve the ecological environment, SLCP plans to stop the farming of sloping farmlands, which is likely to cause soil and water loss and plant trees to restore forest vegetation.Increase forestland and grassland but reduce farmlandThe forestland was restored well in the high–terrain-gradient area and the predominant distribution area of unused land and grassland was gradually diminished [56]. The project also had strong influences on the soil quantity improvement and vegetation coverage increase [57].
Beijing–Tianjin Sandstorm Source Control Project (2000–2010) (BTSSCP)BTSSCP is a sandstorm control policy implemented in Beijing, Tianjin, Hebei, Shanxi and Inner Mongolia to curb the expansion of desertification land and improve the ecological environment.Increase forestland and grasslandThe implementation of BTSSCP has significantly improved the growth conditions of grassland and forestland vegetation [58].
The Three-north Shelterbelt Program (1978–2050) (TNSP)TNSP zone covers the northwest, north and northeast regions of China and consists of three stages (1978–2000, 2001–2020 and 2021–2050); the key objective is to increase regional forest coverage in arid and semi-arid China from 5% to 15%.Increase forestlandThe Chinese government claims that TNSP increased forest cover from 5% in 1978 to more than 13% in 2017. TNSP improved the forest cover of the area effectively and positively affected carbon sequestration [59]. However, Cao et al. and Wang et al. [60,61] strongly argued that the large-scale afforestation failed to address the desertification in some arid and semi-arid regions.
The Policy of Dynamic Equilibrium of the total Cultivated Land (1997–) (PDDCL)To maintain a dynamic balance of cultivated land, provincial governments are required to reclaim new farmland, rehabilitate damaged land, or reuse deserted land to compensate the losses of farmland.Protect farmlandThe farmland expanded to high-terrain-gradient area to compensate the region occupied by the construction of infrastructure facilities in the plain of BTH [56].
IIHengshui, Cangzhou, Xingtai and Handan(1) The Protection of Basic Agricultural Land (1994–) (PBAL);
(2) High-standard Primary Farmland Construction
Project (2011–) (HSPFCP)
(1) According to the demand of population growth for agricultural products, the government set a certain amount of farmland that cannot be occupied.
(2) HSPFCP plans to construct farmlands with concentration, fully supported facilities, stable and high yield, friendly environment and high anti-disaster capability and can adopt to the modern agricultural production and operation modes by the rural land reclamation and readjustment.
Protect farmland and improve the quality of farmlandAlthough PBAL and HSPFCP have protected a certain amount of high-quality farmland and improved their quality to some extent, the fast urbanization in this area has led to significant farmland loss [62,63]. According to Zhang et al. [62], the main reason of the ecosystem service decrease in the BTH region is the increase in artificial land and loss of cropland.
National Wetland Protection (2004–2010) (NWP)NWP aimed to strengthen the protection of wetland by improving their ecological environment.Protect wetlandThe wetlands remained threatened because of unsustainable usage. The enforcement of the wetland protection law was weak [64].
IIIBeijing and Tianjin (plain part)Beijing Plain Afforestation Project (2012–2015)The municipal government afforested to make the capital an “ecological livable city” by improving the urban ecological environment in the Beijing plain areas.Increase forestland and grassland but decrease farmlandIn Beijing, the afforestation efforts covered 1.05 million acres and planted more than 5400 million trees but the farmland decreased at the same time [65,66].
Ecological engineering projects during preparation for 2008 Olympics (2000–2008)To improve the ecological environment during the Beijing Olympic Games period, the capital conducted a series of urban green space construction efforts.Increase forestland and grassland but decrease farmlandThe vegetation increased significantly by occupying the farmland [66,67].
(1) The Protection of Basic Agricultural Land (1994–) (PBAL);
(2) High-standard Primary Farmland Construction
Project (2011–) (HSPFCP)
(1) According to the demand of population growth for agricultural products, the government set a certain amount of farmland that cannot be occupied.
(2) HSPFCP plans to construct farmlands with concentration, fully supported facilities, stable and high yield, friendly environment and high anti-disaster capability and can adopt to the modern agricultural production and operation modes by the rural land reclamation and readjustment.
Protect and improve quality of farmlandFarmland protection in Beijing and Tianjin is difficult to implement thoroughly, hence, the farmlands have been nearly depleted [68,69].

Share and Cite

MDPI and ACS Style

Xie, Z.; Gao, Y.; Li, C.; Zhou, J.; Zhang, T. Spatial Heterogeneity of Typical Ecosystem Services and Their Relationships in Different Ecological–Functional Zones in Beijing–Tianjin–Hebei Region, China. Sustainability 2018, 10, 6. https://doi.org/10.3390/su10010006

AMA Style

Xie Z, Gao Y, Li C, Zhou J, Zhang T. Spatial Heterogeneity of Typical Ecosystem Services and Their Relationships in Different Ecological–Functional Zones in Beijing–Tianjin–Hebei Region, China. Sustainability. 2018; 10(1):6. https://doi.org/10.3390/su10010006

Chicago/Turabian Style

Xie, Zhen, Yang Gao, Chao Li, Jian Zhou, and Tianzhu Zhang. 2018. "Spatial Heterogeneity of Typical Ecosystem Services and Their Relationships in Different Ecological–Functional Zones in Beijing–Tianjin–Hebei Region, China" Sustainability 10, no. 1: 6. https://doi.org/10.3390/su10010006

APA Style

Xie, Z., Gao, Y., Li, C., Zhou, J., & Zhang, T. (2018). Spatial Heterogeneity of Typical Ecosystem Services and Their Relationships in Different Ecological–Functional Zones in Beijing–Tianjin–Hebei Region, China. Sustainability, 10(1), 6. https://doi.org/10.3390/su10010006

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