Next Article in Journal
An Investigation of the Lower Stratospheric Gravity Wave Activity in Tibetan Plateau Based on Multi-GNSS RO Dry Temperature Observations
Next Article in Special Issue
Early Warning of the Carbon-Neutral Pressure Caused by Urban Agglomeration Growth: Evidence from an Urban Network-Based Cellular Automata Model in the Greater Bay Area
Previous Article in Journal
Correction: Martinez et al. Neural Network Approaches to Reconstruct Phytoplankton Time-Series in the Global Ocean. Remote Sens. 2020, 12, 4156
Previous Article in Special Issue
Estimating Large-Scale Interannual Dynamic Impervious Surface Percentages Based on Regional Divisions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 14
Issue 22
10.3390/rs14225670
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

The Supply–Demand Budgets of Ecosystem Service Response to Urbanization: Insights from Urban–Rural Gradient and Major Function-Oriented Areas

by
Zuzheng Li
1,
Baoan Hu
1 and
Yufei Ren
2,3,*
1
School of Ecology and Nature Conservation, Beijing Forestry University, Beijing 100083, China
2
School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China
3
State Key Joint Laboratory of Environmental Simulation and Pollution Control, School of Environment, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(22), 5670; https://doi.org/10.3390/rs14225670
Submission received: 23 September 2022 / Revised: 8 November 2022 / Accepted: 8 November 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Remote Sensing Applications for Interaction between Urbanization and Eco-Environment)
Browse Figures
Review Reports Versions Notes

Abstract

:
The differentiation in the urbanization level’s impact on the supply–demand budgets of ecosystem services (ESs) from the perspective of the major function-oriented areas is of great significance for formulating sustainable development strategies at the regional level. This study first constructed the research framework of the supply, demand, and supply–demand ratios (ESDRs) of ESs responding to urbanization from the perspective of major function-oriented zoning, and then took the rapidly urbanized Beijing–Tianjin–Hebei Urban Agglomeration (BTHUA) of China as a case from 2000 to 2020. The relationships between three urbanization indicators, gross domestic production (GDP), population density (PD), and artificial land proportion (ALP), as well as ESDRs of ESs were investigated using Pearson Correlation analysis across three major functional areas. The sensitivity of ESDRs to urbanization was further evaluated using the Random Forest model. The results showed that the supply of carbon fixation, water provision, and food provision increased, whereas their demands far exceeded their supplies, resulting in an increased imbalance between ES supply and demand. With the exception of soil conservation, significantly negative relationships were observed between urbanization indicators and the other three ES supply–demand budgets. The ESDRs of water provision, carbon fixation, and food provision were the most sensitive variables that depended on the population density (PD) in almost all functional areas, whereas the ESDR of carbon fixation exhibited the highest sensitivity to GDP in developed urban areas and rural areas within the preferred development area (PDA) and key development area (KDA). This study could provide comprehensive information for decision making and ES management in different functional areas.

1. Introduction

Over the past few decades, it has become increasingly evident that human activities are disturbing natural ecosystems [1,2,3]. For human societies to achieve sustainability, it is essential to successfully understand and manage the relationship between human beings and natural ecosystems. Through mutual connections and interactions, humans and natural ecosystems form an integrated social–ecological system [4]. Ecosystem services (ES), as the link between humans and natural ecosystems, have gained increased recognition and are being used extensively [5,6,7]. ES supply is defined as the products and services provided by ecosystems to human society, whereas ES demand can be described as the consumption of the products and services by stakeholders in a particular area within a given time [8,9]. Together, ES supply and demand constitute the dynamic process by which ESs flow between natural ecosystems and social systems [9,10]. Therefore, assessments of ESs that incorporate both supply and demand can accurately reflect natural ecosystem carrying capacities and human disturbance impacts, which are used as a scientific basis for ecosystem management and resource allocation.
Urbanization, one of this century’s most important megatrends, affects the social system and natural ecosystem, resulting in a wide range of ecological issues contributing to global environmental change [11,12]. The ecological issues are frequently aggravated by growing population, economy, and artificial land in urban agglomeration, which place pressure on natural ecosystems and lead to high demand for ESs [12,13]. Due to the contradiction between the increase in urban land and population and the degradation of ESs caused by the encroachment of natural ecosystems, the relationship between supply and demand of ESs is more unbalanced [14,15]. Specifically, on the one hand, urbanization accelerates the transformation of natural and semi-natural lands into artificial lands, thus the ESs are limited through the alteration of ecological processes such as the flow of materials, energy, and biogeochemical cycle [16,17]. For instance, [18] analyzed the relationship between ES supply and urbanization in the Pearl River Delta (PRD), and found that regulating and supporting services decreased, while provisioning services increased due to urban expansion. On the other hand, the growing population and economy in urban areas tend to increase resource consumption and pollutant discharge, thus increasing the demand for ESs [19]. For instance, [20] found an increased demand for carbon storage services due to the increase in energy consumption in the process of urbanization. Therefore, sustainable regional development requires a better understanding of how urbanization impacts the supply and demand of ESs.
Most previous studies examined the relationship between urbanization and the supply of ESs [12,19]. Essentially, these studies described dynamics of ES supply caused by land use changes in the process of urbanization, and the multiple dimensions of urbanization such as economy, populations, and built-up land expansion were less comprehensively considered. In recent years, a few studies explored the changes in supply and demand of ESs under the process of urbanization and the spatial matching relationship [18,20]. Moreover, several studies considered urban agglomeration as a multilevel network system that has substantial gradient differences and heterogeneity [13,19]. For instance, [13] explored the differences in the impact of urbanization indicators on ES supply–demand budgets across regions characterized by different urbanization levels. Currently, urbanization partition only considers socioeconomic factors [19], whereas the major function-oriented areas are based on national spatial governance, considering the natural ecological system and socioeconomic system. Moreover, it has not yet appeared to bring spatial governance into the framework of ecosystem services research, especially from the perspective of major function-oriented areas, to explore the impact of urbanization on the supply and demand of ESs in different major function-oriented areas. Therefore, it is necessary to measure the urbanization level and comprehensively consider the impact of urbanization on ES supply and demand in the major function-oriented areas characterized by different urbanization levels.
The Beijing–Tianjin–Hebei urban agglomeration (BTHUA) experienced rapid urbanization, which has brought several environmental problems and further led to ES supply and demand imbalances [12,21]. Since the State Council promulgated and implemented the “National Main Functional Area Planning” in 2011, the cities of Beijing and Tianjin, and Hebei Province, have implemented the major function-oriented areas planning, respectively, aiming at realizing the sustainable development of the BTHUA. However, there are obvious differences in land space development intensity, functional orientation, and development direction among different main functional areas, which may aggravate the imbalance between social and economic development and ecological environment among different major function-oriented areas. Therefore, the BTHUA was chosen as a case in this study for evaluating the impact of three urbanization indicators, i.e., gross domestic production (GDP), population density (PD), and artificial land proportion (ALP), on ES supply–demand budgets. In detail, this study tries to discuss the following issues: (1) How did the ES supply–demand budgets of the BTHUA change during 2000–2020? (2) Across different functional areas with differing levels of urbanization, how do different urbanization indicators impact ES supply–demand budgets? Analyzing the urbanization affecting the ES supply–demand budgets from the perspective of major function-oriented zoning is of great practical significance for realizing regional sustainable development.

2. Materials and Methods

This study proposed a framework to assess the effects of urbanization on ESDRs and support ES management for each functional area in the Beijing–Tianjin–Hebei urban agglomeration (BTHUA) (Figure 1). Firstly, natural and socioeconomic data were collected, and urbanization of population, economy, and land were measured. Secondly, the spatial variations of urbanization indicators (i.e., GDP, PD, and ALP) and four ES supply–demand ratios (ESDRs) were evaluated from 2000 to 2020. Thirdly, relationships between ESDRs and urbanization indicators were analyzed by using Pearson Correlation and the Random Forest model. Finally, several policy recommendations and measures of functional areas were put forward to improve the coordination development between ESs and urbanization.

2.1. Study Area

The Beijing–Tianjin–Hebei urban agglomeration (BTHUA) consists of Beijing, Tianjin, and Hebei Province, and is located in Northern China (35°03′–42°40′N, 113°27′–119°50′E, Figure 2). It covers a total land area of 212,962 km2, accounting for 2.2% of the total land area of China. The northwestern part of this region has a higher elevation and is mostly hilly, while the southeast is relatively flat with an elevation of less than 100 m. This region has a temperate continental climate, with annual mean precipitation of 420–550 mm and an annual mean temperature of −3.5 °C to 24.5 °C [22]. The southeastern plain is an important grain production area in China, growing corn, wheat, and peanuts.
The BTHUA is a political, cultural, and economic center in China, with 13 cities and 173 counties. In 2019, the BTHUA created a GDP of CNY 8.46 trillion, accounting for 8.5% of the country’s GDP [23,24,25]. From 2000 to 2019, the urbanization of this region developed rapidly. The urbanization rate increased from 38.50% in 2000 to 66.70% in 2019. The urban population of the BTHUA increased from 70.91 million to 113.07 million. With its high-density population and high rate of urbanization, the region has experienced long-term water resource shortages and unbalanced development between the sown area and food yield [26]. Moreover, the imbalance between resource supply and demand has triggered a series of regional problems, including surface runoff decrease, land desertification, groundwater over-extraction, air quality deterioration, biodiversity reduction, and ecosystem degradation [12,17]. Therefore, a comprehensive diagnosis of ESs from the supply–demand perspective is highly significant for the BTHUA’s sustainable ES management.

2.2. Data Sources and Processing

This study used eight sets of geospatial and statistical data, including land use, socioeconomic, Digital Elevation Model (DEM), meteorological, energy consumption, soil properties, watersheds, and Net Primary Productivity (NPP). The sources and descriptions of each data used in this study are presented in Table 1. All vector and raster data were converted to the same projection coordinate system (Bei-jing_1954_3_Degree_GK_CM_ 114E), and the spatial accuracy of all raster data was modified to 1000 m by resampling in ArcGIS 10.3.

2.3. Quantification of Ecosystem Services Supply and Demand

In this work, four key ESs were selected according to the latest version of the Common International Classification of ES (CICES) (https://cices.eu/resources/ (accessed on 6 March 2021)). The selected ESs are very important for the sustainable development of the BTHUA and are sensitive to global climate changes, land use changes, and human activities [12,22,28]. Table 2 provides an overview of the ESs evaluated in the study area and the reasons for their selections by literature reviews.

2.3.1. Carbon Fixation

(1)
Supply
Carbon fixation (CF) refers to the capacity of vegetation to store carbon dioxide in the atmosphere through photosynthesis, which is essential for climate change mitigation [29]. Here, we measured CF’s supply capacity by using the amount of carbon dioxide consumed in photosynthesis [30]. According to the photosynthesis formula, producing 1 g dry matter needs 1.63 g carbon dioxide, then the supply of CF was calculated as follows:
SCS = 1.63 × NPPi
where NPPi is the value of NPP for pixel i.
(2)
Demand
The total amount of carbon emissions was used as the demand for carbon fixation in this study, as the changes in per capita emissions will lead to an increase or decrease in the demand for carbon fixation [31]. The total amount of carbon emission is calculated by multiplying the total energy consumption by the standard carbon emission coefficient according to the data provided by CO2 emission inventory of the BTHUA (e.g., residential, industry, and agriculture) [23,24,25]. The amount of carbon contained in CO2 was estimated to be about 27% of total carbon. We multiplied the average value of emissions per capita by the population density to obtain the total carbon emitted per pixel [32]. The calculation of the total amount of carbon emissions is as follows:
DCP = Dpcfc × Ppop
where DCP is the demand for carbon fixation (t), Dpcfc is the per capita carbon emissions (t), and Ppop is the population density (person·km2).

2.3.2. Water Provision

(1)
Supply
Water provision (WP) refers to the annual quantity of water yield available to humans within a given region [33]. The “Hydropower Water Yield module” of InVEST was used to quantify water provision based on the Budyko curve, with the data including average annual precipitation, root restricting layer depth (mm), plant available water content, annual reference evapotranspiration (mm), and land use maps [34]. The calculations of annual water provision Y for each pixel are as follows:
Y = (1 − AET/P) × P
AET/P = (1 + PET/P) − [1 + (PET/P)ω]1/ω
PET = K × ET0/P
ω = Z × AWC/P + 1.25
AWC = Min (Rest. layer. Soil Depth, Root.Depth) × PAWC
where AET is the annual actual evapotranspiration (mm); P is the annual precipitation (mm); AET/P is based on an expression of the Budyko curve proposed by [35,36]; PET is the potential evapotranspiration; ω is an empirical parameter that characterizes the natural climatic–soil properties; ET0 is the reference evapotranspiration; K is the vegetation evapotranspiration coefficient associated with the land use [34]; AWC is the volumetric plant available water content; Z is the empirical constant, and sometimes referred to as “seasonality factor” (1–30); and PAWC is the plant available water content fraction (0–1) [34].
(2)
Demand
Water demand (DWY) refers to the total amount of water consumption for agricultural and industrial production, inhabitants, and ecological purposes [37,38]. For this estimate, we collected the spatially explicit population density data and the water resource bulletins from each county which provide the water consumption per inhabitant per studied year [23,24,25]. The calculation of water demand DWY is as follows:
DWY = Dpcwc × Ppop
where Dpcwc is the per capita water consumption and Ppop is the population density (person·km2). For details, please refer to Table S1.

2.3.3. Food Provision

Both the supply and demand for food provision (FP) were estimated through statistical data. Here, for FP supply, we first added up each county’s production of grain, vegetables, and fruit products, which were the three main types of foods produced in cropland (https://data.cnki.net/Yearbook/Navi?type=type&code=A, accessed on 8 April 2022). Then, we estimated the demand for FP by multiplying the per capita food consumption by the population density. According to [39], the per capita food demand was 322.07 kg·a−1, which can basically meet China’s per capita food security. For counties without per capita food consumption data, city or provincial-scale per capita food consumption was used as an alternative. Food provision supply and demand can be calculated using the following equations:
SFPi = P(i,j)
DFPi = Dpcfp × Ppop
where SFPi is the FP supply for county i; P(i, j) is the annual provision of j type food for each county, including grains, vegetables, and fruits; DFPi is the FP demand for county i; Dpcfp is the per capita food demand; and Ppop is the population density (person·km2).

2.3.4. Soil Conservation

(1)
Supply
The “Sediment Delivery Ratio module” of InVEST was used to map and calculate the total amount of soil conservation per pixel based on the universal soil loss equation (RUSLE) [34]. Soil conservation equals the difference between the actual amount of soil erosion (USLE) and the maximum potential amount of soil erosion (RKLS), assuming the original land cover without the C or P factors [40,41], and is calculated as:
SC = RKLS − usle
usle = R × K × LS × C × P
RKLS = R × K × LS
where SC is the amount of soil conservation (tons·ha−1·yr−1); usle is the amount of actual soil loss (tons·ha−1·yr−1); RKLS is the amount of potential soil loss (tons·ha−1·yr−1); R is rainfall erosivity (MJ·mm (ha × hr)−1); K is the soil erodibility factor for each pixel (MJ·mm (ha × hr)−1); LS is the slope length-gradient factor (dimensionless); and C and P are the crop management and support practice factors for each pixel (dimensionless), respectively.
(2)
Demand
The amount of actual soil loss was used to define the demand for soil conservation, which is based on the quantity of actual soil erosion that human beings are expected to deal with [42]. The calculation of the amount of actual soil loss is below:
usle = R × K × LS × C × P
where usle is the amount of actual soil loss (tons·ha−1·yr−1). For details, please refer to Table S2.

2.3.5. Ecosystem Service Supply–Demand Ratio (ESDR)

The ecosystem service supply–demand ratio (ESDR) index was used to quantify the relationship between the actual ES supply and human demand, identifying ES supply–demand shortages and mismatches [37]. The ESDR index is calculated as follows:
ESDRi = (Si − Di) / (Smax + Dmax) / 2
where Si and Di refer to the actual ES supply and demand for pixel i, respectively, and Smax and Dmax indicate the maximum value of actual ES supply and human demand in the county, respectively. A value greater than 0 indicates an ES surplus, meaning that supply can meet demand, a value of 0 indicates ES supply–demand balance, and a value lower than 0 indicates a deficit—supply cannot meet demand.

2.4. Urbanization Classification

As a long-term and complicated process of social and economic development, urbanization has often been measured from three perspectives, i.e., economic urbanization, population urbanization, and land urbanization [13]. To be specific, gross domestic production (GDP, ten thousand yuan/km2), population density (PD, persons/km2), and artificial land proportion (ALP, %) were used to represent economic urbanization, population urbanization, and land urbanization, respectively [13,19]. Moreover, three sub-regions with varying levels of urbanization were identified based on nighttime light data in 2000 and 2020 (Figure 3a). Specifically, the nighttime light data in 2020 was first equally divided into 8 categories according to the quantile method, and then the same method was applied to the nighttime light data in 2000. The second to eighth categories in 2000 were identified as developed urban areas (E). Areas other than urban areas of the nighttime light data in 2020 were identified as developing urban areas (I), and the last remaining areas were rural areas (R). There was a significant stratification effect among the three regions of urbanization in 2020 (Figure 3b).

2.5. Function Division

According to the major function-oriented area planning issued by Beijing, Tianjin, and Hebei provincial governments, this study divided the BTHUA into key development area (KDA), preferred development area (PDA), and ecological function area (EFA) (Figure 2b) [43]. KDA is defined as a fully urbanized area with the highest development intensity, and its main function is to optimize development. PDA is defined as the area with the greatest development potential and the level of urbanization that needs to be improved, in which the main function is to focus on the development of the construction land. EFA is an important area to ensure regional ecological security and water resource conservation, in which the main function is to restrict development and to restrict large-scale and high-intensity industrialized urbanization development.

2.6. Statistical Analysis

To explore how urbanization impacts different ESDRs, Pearson Correlation analysis was applied to assess the relationships between three urbanization indicators and ESDRs across three functional areas. Before the analysis, in order to ensure the unity of all variable dimensions, the variables of PD and GDP were log10 transformed. Random Forest (RF) model is a powerful machine learning tool with high prediction accuracy by using an ensemble of decision trees based on bootstrapped samples from a dataset [44,45]. With the help of this model, multiple collinear problems can be avoided, and independent variables can be assessed separately [46]. Then, the RF model was performed to identify the dominant factors within PD, GDP, and ALP driving ESDRs. Accordingly, the proportion of increased MSE (IncMSE) was used as a measure of the relative importance of each independent variable, in which the IncMSE was calculated as the ratio of the increased mean square error (MSE) to its initial MSE values [47]. A higher IncMSE indicates greater importance of the independent variable. The significance of each predictor was determined using the “rfPermute” package, and the RF model was performed using the “randomForest” package in the R version 4.1.2 (https://cran.r-project.org/, accessed on 3 June 2022).

3. Results

3.1. Spatiotemporal Dynamics of Mismatch of ES Supply and Demand

The supply, demand, and ESDRs of the carbon fixation, water provision, food provision, and soil conservation were distributed unevenly and had significant changes in the BTHUA from 2000 to 2020 (Figure 4, Figures S1 and S2). The high supply areas of soil conservation, water provision, and carbon fixation service were mainly located in the northwestern parts, whereas the food provision supply was larger in the southeastern parts (Figure S1). The high demand for carbon fixation, water provision, and food provision was mainly distributed in the central area of these cities and increased rapidly (Figure S2). The spatial distribution of supply and demand of water yield, food provision, and carbon fixation showed a serious mismatch, especially in the southeastern parts and central areas of these cities. The soil conservation service has always been in balance in total and in surplus in the mountains, with a slight increase in the ESDR in the southeast parts and a decrease in the other areas. From 2000 to 2020, the ESDRs of carbon fixation, water provision, and food provision were increasing in the city centers of Hebei province and were decreasing in Beijing and Tianjin city centers.
The total supply of carbon fixation, water provision, and food provision increased from 53.32 billion tons, 58.43 billion m3, and CNY 113.66 billion in 2000 to 72.29 billion tons, 59.25 billion m3, and CNY 408.10 billion in 2020, respectively (Table S3). The food provision supply has increased substantially, while water production and soil conservation have changed very little. In contrast, the demand for water provision decreased from 35,348.99 billion m3 to 31,065.73 billion m3. The demand for carbon fixation and food provision increased from 242.49 billion tons and CNY 286,845.22 billion in 2000 to 1047.48 billion tons and CNY 1,148,838.66 billion in 2020, with a dramatic increase of 331.96% and 300.51, respectively (Table S4). Hence, although carbon fixation, water provision, and food provision have been changed to various extents from 2000 to 2020, their demand far exceeded their supply in total (Table S5).

3.2. Relationships between Urbanization and ESDRs

In general, urbanization had a significant negative influence on ESDRs of ESs except soil conservation (SC) in all functional areas (Figure 5). The relationships between urbanization indicators and ESDRs were different under different urbanization levels within different functional areas. Population density (PD) has the strongest correlation with ESDRs of water provision (WP), food provision (FP), and soil conservation (SC) in all functional areas. GDP and ALP had the strongest correlation with ESDRs of water provision (WP), food provision (FP), and carbon fixation (CF) in developed urban areas within preferred development (PDA) and key development areas (KDA). In rural areas within ecological function areas, GDP had the strongest correlation with ESDRs of water provision (WP), food provision (FP), and carbon fixation (CF). Moreover, all urbanization indicators had a low correlation with ESDRs of soil conservation (SC).

3.3. Sensitivity of ESDR to Urbanization

The ESDRs of water provision and food provision services were more sensitive to PD in almost all functional areas, with an increased MSE (IncMSE from 47.08 to 96.72%), followed by GDP and ALP (Figure 6). The results also showed that the ESDR of carbon fixation was most sensitive to PD in developed urban areas and rural areas within the PDA and KDA, with the IncMSE ranging from 33.62% to 89.84%. The ESDR of carbon fixation was most sensitive to GDP in developing urban areas within the PDA and KDA, with IncMSE of 72.18% and 80.35%, respectively. Urbanization indicators had the highest explanation for ESDR of water provision, food provision, and soil conservation services of developed urban areas and developing urban areas within the preferred development areas (variance explained from 6.69% to 98.85%). Urbanization indicators had the highest explanation for the ESDR of water provision and food provision services of rural areas in ecological function areas (variance explained from 36.81% to 99.21%).

4. Discussion

4.1. The Impact of Urbanization on ESDRs in Different Functional Areas

As a result of sensitivity analysis, we identified the most sensitive urbanization indicators of ESDRs. ESDRs of water provision, food provision, and carbon fixation were most sensitive to PD, with the exception of developed urban and rural areas within the PDA and KDA. This is attributed to the development of urbanization led to the aggregation of the population in most areas, thus leading to the insufficient ES supply [13,19]. The most important problem appears to be related to the dense population, which has led to a high rise in demand and a steep decline in supply for ESs, creating an imbalanced state between ES supply and demand. Particularly, it was clear that the ESDR of carbon fixation was more sensitive to GDP and ALP in developing urban areas within the PDA and KDA. Since the developing urban areas of PDA and KDA were experiencing rapid urbanization, economic growth and land use changes resulted in a large number of carbon emissions and affected the ESDR of carbon fixation, which is consistent with the research by [19]. The ESDR of soil conservation in rural areas within the KDA and PDA was more sensitive to ALP than in other areas. This is because rural areas have been in the primary stage of urbanization in the past few decades, which is dominated by the expansion of artificial land; with a growing number of natural and semi-natural lands that have been converted into artificial lands, erosion occurs easily [48,49].
In general, ESDRs were more sensitive to urbanization indicators in developed urban and rural areas within the PDA than in KDA and EFA, which is consistent with the strategic positioning of the main functional areas. On the one hand, the KDA and PDA have a high degree of land space development, which undertake the tasks of absorbing population, economic development, and industrial agglomeration, and are the main pressure areas of water resources, environmental pollution, and social environmental problems [50]. On the other hand, the EFA is responsible for regulating regional climate, conserving water sources, and preventing soil desertification, and it is the type of area with the largest carbon pool [51]. Moreover, the EFA has a wide planting area of grain and vegetables, a large amount of forest land in the field, a high forest coverage rate, and therefore the second largest carbon pool. The ESDR of carbon fixation in developing urban areas within the KDA was more sensitive to urbanization indicators than that of the EFA and PDA. This is mainly because the KDA’s main development goal is to promote industrialization and urbanization, while the PDA’s goal is to promote energy conservation and emission reduction, strictly control the production of pollutants, and strengthen ecological construction [43].

4.2. Suggestions for ES Management and Urban Planning

Several suggestions were proposed to mitigate the imbalance between ES supply and demand, which would contribute to the sustainable development of the BTHUA. The PDA, which is characterized by extremely high urban population density, should first advocate intensive land use and restrict the development of land. For instance, some of the original industrial storage land in developing urban areas has been left idle for the purpose of industrial structure adjustment, and urban planning should be recovered and reallocated [52]. Secondly, the industrial workshop should be three-dimensional, and the plane layout and vertical layout of the workshop should be overlapped to form a complex industrial building. Thirdly, high-rise and small high-rise residential buildings by developers are encouraged [53].
The KDA, as a key development and construction area of urbanization, should first coordinate space development, for instance, strengthening the integration of land and space resources and increasing the space for urban construction in developing urban areas [54]. Secondly, accelerating population agglomeration is necessary, such as expanding the population gathering capacity of developed urban areas and promoting the concentration of rural populations in developing urban areas [55]. Thirdly, protecting the ecological environment is needed, such as coordinating urban and rural environmental protection, strengthening the protection of shelterbelts, increasing urban green space, and protecting and restoring ecological functions [56].
The EFA, as the main supply area of multiple ESs, should first carry out ecological restoration and construction, for instance, strengthening the construction of ecological function zones, effectively restoring and upgrading ecological functions, and improving the production capacity of ecological products [57]. Secondly, it is necessary to maintain the integrity of the ecosystem and improve ecological functions such as water conservation, soil and water conservation, and wind and sand prevention. Thirdly, optimizing the industrial structure is needed, such as strictly controlling the development intensity in developing urban areas, developing the ecological economy, and developing other suitable industries that do not affect the positioning of main functions [58]. Fourthly, orderly guiding of the population transfer to urban areas is needed, such as accelerating the ecological migration and rationally guiding the population in plateau mountainous areas to move to developed and developing urban areas [59,60].
However, the major function-oriented area can play a variety of functions, because the main functions of a certain area and the main contents and tasks of its development do not exclude the area from playing other functions. For instance, as urbanized areas, the major functions of the PDA and KDA are to provide industrial products and service products and to gather the population and economy. Meanwhile, it is also necessary to protect agricultural space such as basic farmland and ecological space such as forests, grasslands, water surfaces, and wetlands, and to provide a certain amount of agricultural products and ecological products [61]. In addition, the major function of the EFA is to provide agricultural products and ecological products, to ensure the stability of the ecological system, to allow the moderate development of industries that do not affect the major functional orientation, and to allow the necessary urban construction. Therefore, we should carry out orderly development based on the carrying capacity of natural resources and environmental capacity, vigorously develop a green economy and low-carbon economy, and promote the coordination of population distribution and town and economic layout with resources and the environment.

5. Conclusions

In this study, the differentiation in urbanization level’s impact on the ESDRs was explored across three functional areas in the BTHUA. The results showed that the demands for carbon fixation, water provision, and food provision were growing much faster in urban areas than their supplies from 2000 to 2020, indicating a mismatch of ES supply and demand both in quantity and space. With the exception of the ESDR of soil conservation, urbanization indicators were negatively related to the other three ESDRs. In addition, sensitivity analysis showed that the ESDRs of carbon fixation, water provision, and food provision were most sensitive to population density (PD) in almost all functional areas, whereas carbon fixation also exhibited the highest sensitivity to GDP in developed urban areas and rural areas within the preferred development area (PDA) and key development area (KDA). Five suggestions were put forward for decision making and ES management in different functional areas.
Despite its positive findings, this study does have some limitations that should be highlighted and explored further. First, this study did not comprehensively analyze all types of ecosystem supply and demand, since the coupled system of human and natural factors is incredibly complicated. Second, the supply of ESs, such as water provision and soil conservation, is influenced by land use and biophysical parameters, which are modified according to the actual situation in the BTHUA. Thus, the calculation of ES supply may not be accurate, and the regional differences in ES supply patterns cannot be precisely described. Third, local stakeholders with short-term or long-term interests in ESs were not considered in this study. The future could therefore be devoted to quantifying ES supply and demand from multiple stakeholders. Furthermore, it is imperative to identify who provides these ESs, which areas are entitled to consume them, and whether supplies and consumption are limited to a particular area or extend to other regions. Therefore, ES flow identifications should be incorporated into ES assessments to reveal more details about the interactions between natural and social systems and provide more information on ES management.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs14225670/s1, Table S1: Input data for each LULC class in the InVEST 3.8.0 water yield model. Table S2: Input data for each LULC class in the InVEST 3.8.0 sediment delivery ratio model. Table S3: Total supply of ESs in the BTHUA. Table S4: Total demand of ESs in the BTHUA. Table S5: Total ESDRs of ESs in the BTHUA. Figure S1: Spatial patterns of the supply of carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) from 2000 to 2020. Figure S2: Spatial patterns of the demand of carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) from 2000 to 2020.

Author Contributions

Conceptualization, Methodology, and Writing—original draft, Z.L.; Formal analysis, B.H.; Validation and Resources, B.H.; Software and Data curation, B.H.; Project administration, Y.R.; Conceptualization, Y.R.; Funding acquisition, Y.R.; Supervision, Y.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 42101190) and the China Postdoctoral Science Foundation (Grant No. 2020M680018).

Acknowledgments

We would like to thank Yangyi Qin for his assistance with technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paruelo, J.M.; Burke, I.C.; Lauenroth, W.K. Land-use impact on ecosystem functioning in eastern Colorado, USA. Glob. Chang. Biol. 2001, 7, 631–639. [Google Scholar]
  2. Tollefson, J. Humans are driving one million species to extinction. Nature 2019, 569, 171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kuang, W.; Zhang, S.; Du, G.; Yan, C.; Wu, S.; Li, R.; Lu, D.; Pan, T.; Ning, J.; Guo, C.; et al. Monitoring periodically national land use changes and analyzing their spatiotemporal patterns in China during 2015–2020. J. Geogr. Sci. 2022, 32, 1705–1723. [Google Scholar] [CrossRef]
  4. Bruckmeier, K. Social-Ecological Systems and Ecosystem Services. In Social Ecological Transformation; Palgrave Macmillan: London, UK, 2016. [Google Scholar]
  5. Costanza, R.; Arge, R.; De Groot, R.; Farber, S.; Hannon, B.; Limburg, K.; Suttonkk, P. The Value of the World’s Ecosystem Services and Natural Capital by. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  6. Daily, G.C.; Polasky, S.; Goldstein, J.; Kareiva, P.M.; Mooney, H.A.; Pejchar, L.; Ricketts, T.H.; Salzman, J.; Shallenberger, R. Ecosystem services in decision making: Time to deliver. Front. Ecol. Environ. 2009, 7, 21–28. [Google Scholar] [CrossRef] [Green Version]
  7. Rao, Y.; Zhou, M.; Ou, G.; Dai, D.; Zhang, L.; Zhang, Z.; Yang, C. Integrating ecosystem services value for sustainable land-use management in semi-arid region. J. Clean. Prod. 2018, 186, 662–672. [Google Scholar] [CrossRef]
  8. Burkhard, B.; Kroll, F.; Nedkov, S.; Müller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic. 2012, 21, 17–29. [Google Scholar] [CrossRef]
  9. Schirpke, U.; Candiago, S.; Egarter Vigl, L.; Jäger, H.; Labadini, A.; Marsoner, T.; Meisch, C.; Tasser, E.; Tappeiner, U. Integrating supply, flow and demand to enhance the understanding of interactions among multiple ecosystem services. Sci. Total Environ. 2019, 651, 928–941. [Google Scholar] [CrossRef]
  10. Goldenberg, R.; Kalantari, Z.; Cvetkovic, V.; Mörtberg, U.; Deal, B.; Destouni, G. Distinction, quantification and mapping of potential and realized supply-demand of flow-dependent ecosystem services. Sci. Total Environ. 2017, 593–594, 599–609. [Google Scholar] [CrossRef]
  11. Bai, X.M.; Chen, J.; Shi, P.J. Landscape urbanization and economic growth in China: Positive feedbacks and sustainability dilemmas. Environ. Sci. Technol. 2012, 46, 132–139. [Google Scholar] [CrossRef]
  12. 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]
  13. Zhang, Z.; Peng, J.; Xu, Z.H.; Wang, X.Y.; Meersmans, J. Ecosystem services supply and demand response to urbanization: A case study of the Pearl River Delta, China. Ecosyst. Serv. 2021, 49, 101274. [Google Scholar] [CrossRef]
  14. Chen, G.; Li, X.; Liu, X.; Chen, Y.; Liang, X.; Leng, J.; Xu, X.; Liao, W.; Qiu, Y.; Wu, Q.; et al. Global projections of future urban land expansion under shared socioeconomic pathways. Nat. Commun. 2020, 11, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yu, L.; Lyu, Y.; Chen, C.; Choguill, C.L. Environmental deterioration in rapid urbanisation: Evidence from assessment of ecosystem service value in Wujiang, Suzhou. Environ. Dev. Sustain. 2021, 23, 331–349. [Google Scholar] [CrossRef] [Green Version]
  16. Kuang, W. 70 years of urban expansion across China: Trajectory, pattern, and national policies. Sci. Bull. 2020, 65, 1970–1974. [Google Scholar] [CrossRef]
  17. Yang, Y.; Zheng, H.; Kong, L.; Huang, B.; Xu, W.; Ouyang, Z. Mapping ecosystem services bundles to detect high- and low-value ecosystem services areas for land use management. J. Clean. Prod. 2019, 225, 11–17. [Google Scholar] [CrossRef]
  18. Wang, S.; Hu, M.; Wang, Y.; Xia, B. Dynamics of ecosystem services in response to urbanization across temporal and spatial scales in a mega metropolitan area. Sustain. Cities Soc. 2021, 77, 103561. [Google Scholar] [CrossRef]
  19. Kuang, W.; Du, G.; Lu, D.; Dou, Y.; Li, X.; Zhang, S.; Chi, W.; Dong, J.; Chen, G.; Yin, Z.; et al. Global observation of urban expansion and land-cover dynamics using satellite big-data. Sci. Bull. 2021, 66, 297–300. [Google Scholar] [CrossRef]
  20. Deng, C.; Liu, J.; Liu, Y.; Li, Z.; Nie, X.; Hu, X.; Wang, L.; Zhang, Y.; Zhang, G.; Zhu, D.; et al. Spatiotemporal dislocation of urbanization and ecological construction increased the ecosystem service supply and demand imbalance. J. Environ. Manag. 2021, 288, 112478. [Google Scholar] [CrossRef]
  21. Ren, Y.F.; Li, Z.Z. Unraveling the dynamics, heterogeneity, determinants of eco-efficiency in Beijing-Tianjin-Hebei urban agglomeration, China. J. Environ. Manag. 2022, 317, 115407. [Google Scholar] [CrossRef]
  22. Shen, J.; Li, S.; Liang, Z.; Liu, L.; Li, D.; Wu, S. Exploring the heterogeneity and nonlinearity of trade-offs and synergies among ecosystem services bundles in the Beijing-Tianjin-Hebei urban agglomeration. Ecosyst. Serv. 2020, 43, 101103. [Google Scholar] [CrossRef]
  23. National Bureau of Statistics. China Statistical Yearbook; China Statistics Press: Beijing, China, 2020.
  24. National Bureau of Statistics. China Statistical Yearbook; China Statistics Press: Tianjin, China, 2020.
  25. National Bureau of Statistics. China Statistical Yearbook; China Statistics Press: Hebei, China, 2020.
  26. Li, D.; Wu, S.; Liu, L.; Liang, Z.; Li, S. Evaluating regional water security through a freshwater ecosystem service flow model: A case study in Beijing-Tianjian-Hebei region, China. Ecol. Indic. 2017, 81, 159–170. [Google Scholar] [CrossRef]
  27. Ministry of Natural and Resources of R.P.C (MNRC). GB/T21010-2017, Current Land Use Classification. Available online: http://std.samr.gov.cn/ (accessed on 8 June 2021).
  28. Kang, P.; Chen, W.; Hou, Y.; Li, Y. Linking ecosystem services and ecosystem health to ecological risk assessment: A case study of the Beijing-Tianjin-Hebei urban agglomeration. Sci. Total Environ. 2018, 636, 1442–1454. [Google Scholar] [CrossRef] [PubMed]
  29. Onaindia, M.; de Manuel, B.F.; Madariaga, I.; Rodríguez-Loinaz, G. Co-benefits and trade-offs between biodiversity, carbon storage and water flow regulation. For. Ecol. Manag. 2013, 289, 1–9. [Google Scholar] [CrossRef]
  30. Tallis, H.; Polasky, S. Mapping and valuing ecosystem services as an approach for conservation and natural-resource management. Ann. N. Y. Acad. Sci. 2009, 1162, 265–283. [Google Scholar] [CrossRef] [PubMed]
  31. Bateman, I.J.; Harwood, A.R.; Mace, G.M.; Watson, R.T.; Abson, D.J.; Andrews, B.; Binner, A.; Crowe, A.; Day, B.H.; Dugdale, S.; et al. Bringing ecosystem services into economic decision-making: Land use in the United Kingdom. Science 2013, 341, 45–50. [Google Scholar] [CrossRef] [PubMed]
  32. González-García, A.; Palomo, I.; González, J.A.; López, C.A.; Montes, C. Quantifying spatial supply-demand mismatches in ecosystem services provides insights for land-use planning. Land Use Policy 2020, 94, 104493. [Google Scholar] [CrossRef]
  33. Boithias, L.; Acuña, V.; Vergoñós, L.; Ziv, G.; Marcé, R.; Sabater, S. Assessment of the water supply: Demand ratios in a Mediterranean basin under different global change scenarios and mitigation alternatives. Sci. Total Environ. 2014, 470, 567–577. [Google Scholar] [CrossRef]
  34. Sharp, R.; Tallis, H.T.; Ricketts, T.; Guerry, A.D.; Wood, S.A.; Chaplin-Kramer, R.; Nelson, E.; Ennaanay, D.; Wolny, S.; Olwero, N.; et al. InVEST 3.8.0 User’s Guide; The Natural Capital Project: Stanford University, University of Minnesota, The Nature Conservancy, and World Wildlife Fund. 2020; Available online: http://releases.naturalcapitalproject.org/invest-userguide/latest/#supporting-tools (accessed on 6 May 2022).
  35. Fu, B.P. On the calculation of the evaporation from land surface (in Chinese). Chin. J. Atmos. Sci. 1981, 5, 23 31. [Google Scholar]
  36. Zhang, L.; Hickel, K.; Dawes, W.R.; Chiew, F.H.S.; Western, A.W.; Briggs, P.R. A rational function approach for estimating mean annual evapotranspiration. Water Resour. Res. 2004, 40, 2. [Google Scholar] [CrossRef]
  37. Chen, J.Y.; Jiang, B.; Bai, Y.; Xu, X.B.; Alatalo, J.M. Quantifying ecosystem services supply and demand shortfalls and mismatches for management optimisation. Sci. Total Environ. 2019, 650, 1426–1439. [Google Scholar] [CrossRef] [PubMed]
  38. Burkhard, B.; Crossman, N.; Nedkov, S.; Petz, K.; Alkemade, R. Mapping and modelling ecosystem services for science, policy and practice. Ecosyst. Serv. 2013, 4, 1–3. [Google Scholar] [CrossRef]
  39. Tang, H.J.; Li, Z.M. Study on Per Capita Grain Demand Based on Chinese Reasonable Dietary Pattern. Sci. Agric. Sin. 2012, 45, 2315–2327. [Google Scholar]
  40. Kareiva, P.; Tallis, H.; Ricketts, T.H.; Daily, G.C.; Polasky, S. Natural Capital: Theory and Practice of Mapping Ecosystem Services, 1st ed.; Oxford University: Oxford, UK, 2011; pp. 3–128. [Google Scholar]
  41. Ouyang, Z.; Fan, P.; Chen, J. Urban built-up areas in transitional economies of southeast Asia: Spatial extent and dynamics. Remote Sens. 2016, 8, 819. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, L.; Liu, C.; Wang, C.; Li, P. Supply and demand matching of ecosystem services in loess hilly region: A case study of Lanzhou. Acta Geogr. Sin. 2019, 74, 1921–1937. (In Chinese) [Google Scholar]
  43. Xia, S.; Yang, Y. Spatio-temporal differentiation of carbon budget and carbon compensation zoning in Beijing-Tianjin-Hebei Urban Agglomeration based on the Plan for Major Function-oriented Zones. Acta Geogr. Sin. 2022, 77, 679–696. (In Chinese) [Google Scholar]
  44. Biau, G.; Scornet, E. A random forest guided tour. TEST 2016, 25, 197–227. [Google Scholar] [CrossRef] [Green Version]
  45. Carvalho, T.P.; Soares, F.; Vita, R.; Francisco, R.D.; Basto, J.P.; Alcalá, S.G. A systematic literature review of machine learning methods applied to predictive maintenance. Comput. Ind. Eng. 2019, 137, 106024. [Google Scholar] [CrossRef]
  46. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  47. Liaw, A.; Wiener, M. Classification and regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  48. Tian, G.; Qiao, Z. Assessing the impact of the urbanization process on net primary productivity in China in 1989–2000. Environ. Pollut. 2014, 184, 320–326. [Google Scholar] [CrossRef] [PubMed]
  49. Yun, G.; Zhao, S. The imprint of urbanization on PM2.5 concentrations in China: The urban-rural gradient study. Sustain. Cities Soc. 2022, 86, 104103. [Google Scholar] [CrossRef]
  50. Yuan, Z.; Chen, M.; Tang, Z.; Mei, Z. Urbanization, land use change, and carbon emissions: Quantitative assessments for city-level carbon emissions in Beijing-Tianjin-Hebei region. Sustain. Cities Soc. 2021, 66, 102701. [Google Scholar]
  51. Xu, K.; Wang, J.; Wang, J.; Wang, X.; Chi, Y.; Zhang, X. Environmental function zoning for spatially differentiated environmental policies in China. J. Environ. Manag. 2020, 1, 255, 109485. [Google Scholar] [CrossRef]
  52. Qu, Y.; Jiang, G.; Li, Z.; Shang, R.; Zhou, D. Understanding the multidimensional morphological characteristics of urban idle land: Stage, subject, and spatial heterogeneity. Cities 2020, 97, 102492. [Google Scholar] [CrossRef]
  53. Yuan, M.; Li, Z.; Li, X.; Li, L.; Zhang, S.; Luo, X. How to promote the sustainable development of prefabricated residential buildings in China: A tripartite evolutionary game analysis. J. Clean. Prod. 2022, 349, 131423. [Google Scholar] [CrossRef]
  54. Liu, X.; Li, X.; Yang, J.; Fan, H.; Zhang, J.; Zhang, Y. How to resolve the conflicts of urban functional space in planning: A perspective of urban moderate boundary. Ecol. Indic. 2022, 144, 109495. [Google Scholar] [CrossRef]
  55. Wang, Y.; Niu, Y.; Li, M.; Yu, Q.; Chen, W. Spatial structure and carbon emission of urban agglomerations: Spatiotemporal characteristics and driving forces. Sustain. Cities Soc. 2022, 78, 103600. [Google Scholar] [CrossRef]
  56. Xing, L.; Hu, M.; Xue, M. The effect of urban–rural construction land transition on ecosystem services: A theoretical framework and empirical evidence for China. Habitat Int. 2022, 124, 102576. [Google Scholar] [CrossRef]
  57. Yang, Y.; Zhang, Y.; Yang, H.; Yang, F. Horizontal ecological compensation as a tool for sustainable development of urban agglomerations: Exploration of the realization mechanism of Guanzhong Plain urban agglomeration in China. Environ. Sci. Policy 2022, 137, 301–313. [Google Scholar] [CrossRef]
  58. Wang, Y.; Yao, Y.; Chen, S.; Ni, Z.; Xia, B. Spatiotemporal evolution of urban development and surface urban heat island in Guangdong-Hong Kong-Macau greater bay area of China from 2013 to 2019. Resour. Conserv. Recycl. 2022, 179, 106063. [Google Scholar] [CrossRef]
  59. Wang, W.; Jiao, L.; Jia, Q.; Liu, J.; Mao, W.; Xu, Z.; Li, W. Land use optimization modelling with ecological priority perspective for large-scale spatial planning. Sustain. Cities Soc. 2021, 65, 102575. [Google Scholar] [CrossRef]
  60. Li, Z.; Cheng, X.; Han, H. Future impacts of land use change on ecosystem services under different scenarios in the ecological conservation area, Beijing, China. Forests 2020, 11, 584. [Google Scholar] [CrossRef]
  61. Kuang, W.; Liu, J.; Tian, H.; Shi, H.; Dong, J.; Song, C.; Li, X.; Du, G.; Hou, Y.; Lu, D.; et al. Cropland redistribution to marginal lands undermines environmental sustainability. Natl. Sci. Rev. 2022, 9, 66–78. [Google Scholar] [CrossRef] [PubMed]
Remotesensing 14 05670 g001 550
Figure 1. The framework and procedures of this study.
Figure 1. The framework and procedures of this study.
Remotesensing 14 05670 g001
Remotesensing 14 05670 g002 550
Figure 2. Study area. (a) LULC distribution of cover types. (b) Map of major function-oriented areas, including preferred development area (PDA), ecological function area (EFA), and key development area (KDA).
Figure 2. Study area. (a) LULC distribution of cover types. (b) Map of major function-oriented areas, including preferred development area (PDA), ecological function area (EFA), and key development area (KDA).
Remotesensing 14 05670 g002
Remotesensing 14 05670 g003 550
Figure 3. Spatial characteristics of different urbanization levels of the BTHUA. (a) Distribution of urbanization levels. (b) Urbanization indicators, i.e., log10 (PD), log10 (GDP), and log10 (ALP) of E, I, and R. E, I, and R represent developed urban areas, developing urban areas, and rural areas, respectively.
Figure 3. Spatial characteristics of different urbanization levels of the BTHUA. (a) Distribution of urbanization levels. (b) Urbanization indicators, i.e., log10 (PD), log10 (GDP), and log10 (ALP) of E, I, and R. E, I, and R represent developed urban areas, developing urban areas, and rural areas, respectively.
Remotesensing 14 05670 g003
Remotesensing 14 05670 g004 550
Figure 4. Spatial pattern and change of four ESDRs for carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) from 2000 to 2020. (a) Spatial patterns of ESDRs, (b) spatial changes of ESDRs, and (c) changes in mean values of ESDRs.
Figure 4. Spatial pattern and change of four ESDRs for carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) from 2000 to 2020. (a) Spatial patterns of ESDRs, (b) spatial changes of ESDRs, and (c) changes in mean values of ESDRs.
Remotesensing 14 05670 g004
Remotesensing 14 05670 g005 550
Figure 5. Relationships between urbanization and ES supply–demand ratios (ESDRs) in BTHUA in 2020. Three functional areas were distinguished, i.e., key development area (KDA), preferred development area (PDA), and ecological function area (EFA). Four ESs are water provision (WP), food provision (FP), carbon fixation (CF), and soil conservation (SC). Urbanization indicators include population density (PD), gross domestic production (GDP), and artificial land proportion (ALP).
Figure 5. Relationships between urbanization and ES supply–demand ratios (ESDRs) in BTHUA in 2020. Three functional areas were distinguished, i.e., key development area (KDA), preferred development area (PDA), and ecological function area (EFA). Four ESs are water provision (WP), food provision (FP), carbon fixation (CF), and soil conservation (SC). Urbanization indicators include population density (PD), gross domestic production (GDP), and artificial land proportion (ALP).
Remotesensing 14 05670 g005
Remotesensing 14 05670 g006 550
Figure 6. Sensitivity of ES supply–demand ratios (ESDRs) for carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) to urbanization in BTHUA in 2020. Three functional areas were distinguished, i.e., key development area (KDA), preferred development area (PDA), and ecological function area (EFA).
Figure 6. Sensitivity of ES supply–demand ratios (ESDRs) for carbon fixation (CF), water provision (WP), food provision (FP), and soil conservation (SC) to urbanization in BTHUA in 2020. Three functional areas were distinguished, i.e., key development area (KDA), preferred development area (PDA), and ecological function area (EFA).
Remotesensing 14 05670 g006
Table
Table 1. The data requirements, sources, and descriptions.
Table 1. The data requirements, sources, and descriptions.
Data TypesData SourcesSources and Descriptions
Land useCASEarth (https://data.casearth.cn/, accessed on 2 March 2022)Raster (30 × 30 m). Land use types in 2000, 2010, and 2020 were divided into six categories: grassland, water body, cultivated land, artificial surface, unused land, forest land, and shrub land [27].
Socioeconomic dataWordPop (https://www.worldpop.org/, accessed on 4 March 2022)Raster (1000 × 1000 m). Including population density (PD) and gross domestic product (GDP).
Digital Elevation Model data (DEM)Geospatial Data Cloud (http://www.gscloud.cn/, accessed on 6 March 2022)Raster (30 × 30 m). A value of elevation for each grid cell.
Meteorological data NCAR (https://climatedataguide.ucar.edu/, accessed on 8 March 2022)Raster (1000 × 1000 m). Including annual average precipitation and evapotranspiration.
Statistical dataEconomic Yearbook of Beijing, Hebei, and Tianjin, and the China Energy Statistical YearbookExcel format. Per capita annual energy consumption and water use data.
Soil propertiesHarmonized World Soil Database (HWSD; https://www.fao.org/soils-portal/data-hub/soil-properties/en/, accessed on 8 March 2022)Raster (1000 × 1000 m). Including root restricting layer depth, plant available water content (PAWC) range, etc.
WatershedsHydroSHEDS (https://hydrosheds.org/, accessed on 8 March 2022)Vector format. A number is assigned for each watershed.
Net Primary Productivity (NPP)USGS (https://www.usgs.gov/, accessed on 8 March 2022)Raster (500 × 500 m). The amount of organic matter accumulated by plants in unit area and time.
Table
Table 2. Ecosystem services evaluated in the Beijing–Tianjin–Hebei urban agglomeration.
Table 2. Ecosystem services evaluated in the Beijing–Tianjin–Hebei urban agglomeration.
Ecosystem ServicesSelection Reasons
Carbon storage (CS)The absorption of CO2 by vegetation is of great significance to regional climate change and directly affects human health.
Water provision (WP)Water resources are recharged by terrestrial and aquatic ecosystems, affecting the growth of vegetation.
Food provision (FP)Food production is mainly provided by cultivated land and is the basic material for human survival.
Soil conservation (SC)Reduction in soil erosion caused by storm runoff and topography is important in the BTHUA.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Hu, B.; Ren, Y. The Supply–Demand Budgets of Ecosystem Service Response to Urbanization: Insights from Urban–Rural Gradient and Major Function-Oriented Areas. Remote Sens. 2022, 14, 5670. https://doi.org/10.3390/rs14225670

AMA Style

Li Z, Hu B, Ren Y. The Supply–Demand Budgets of Ecosystem Service Response to Urbanization: Insights from Urban–Rural Gradient and Major Function-Oriented Areas. Remote Sensing. 2022; 14(22):5670. https://doi.org/10.3390/rs14225670

Chicago/Turabian Style

Li, Zuzheng, Baoan Hu, and Yufei Ren. 2022. "The Supply–Demand Budgets of Ecosystem Service Response to Urbanization: Insights from Urban–Rural Gradient and Major Function-Oriented Areas" Remote Sensing 14, no. 22: 5670. https://doi.org/10.3390/rs14225670

APA Style

Li, Z., Hu, B., & Ren, Y. (2022). The Supply–Demand Budgets of Ecosystem Service Response to Urbanization: Insights from Urban–Rural Gradient and Major Function-Oriented Areas. Remote Sensing, 14(22), 5670. https://doi.org/10.3390/rs14225670

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