Next Article in Journal
Effect of Performance by Excessive Advanced Fuel Injection Timing on Marine Diesel Engine
Next Article in Special Issue
Using Public Landslide Inventories for Landslide Susceptibility Assessment at the Basin Scale: Application to the Torto River Basin (Central-Northern Sicily, Italy)
Previous Article in Journal
Shadow Hunter: Low-Illumination Object-Detection Algorithm
Previous Article in Special Issue
Formal Matters on the Topic of Risk Mitigation: A Mathematical Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation and Driving Determinants of the Coordination between Ecosystem Service Supply and Demand: A Case Study in Shanxi Province

1
School of Culture Tourism and Journalism Arts, Shanxi University of Finance and Economics, Taiyuan 030006, China
2
College of Mining Engineer, Taiyuan University of Technology, Taiyuan 030024, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(16), 9262; https://doi.org/10.3390/app13169262
Submission received: 22 July 2023 / Revised: 10 August 2023 / Accepted: 12 August 2023 / Published: 15 August 2023
(This article belongs to the Special Issue Advancing Complexity Research in Earth Sciences and Geography)

Abstract

:
Understanding the coordination relationship between ecosystem service (ES) supply and demand and elucidating the impact of driving factors is critical for regional land use planning and ecological sustainability. We use a large watershed area as a case to map and analyze ES supply, demand and the coordination relationship, and identify the associated socio-ecological driving variables. This study assessed the supply and demand of five ESs (crop production, water retention, soil conservation, carbon sequestration, and outdoor recreation) in 2000 and 2020, and evaluated the coordination between them employing the coupling coordination degree model (CCDM). Additionally, we utilized the geo-detector model (GDM) to identify driving determinants and their interactive effects on the spatial pattern of the coupling coordination degree (CCD) between ES supply and demand. The results showed that mountainous regions with abundant forest coverage were high-value areas for ES supply, while the ESs were predominantly required in city center areas within each basin area. From 2000 to 2020, there was a slight decline in ES supply and a significant increase in ES demand. Counties were grouped into four coordination zones in the study area: extreme incoordination, moderate incoordination, reluctant coordination, and moderate coordination. The number of counties with extreme incoordination linked to regions with a mountain ecosystem is increasing, where the ES supply is much greater than the demand. The moderate incoordination counties dominated by a cropland ecosystem exhibited slightly higher levels of ES supply than demand. The moderate and reluctant coordination were linked to counties with distinct ecological characteristics. Construction land played a major role in the characteristics of the CCD, followed by grassland. The interaction between construction land and all other factors significantly increased the influence on the CCD. These findings offered valuable insights for land managers to identify areas characterized by incoordination between ES supply and demand and understand associated factors to develop optimal ES management strategies.

1. Introduction

Ecosystem services (ESs) are defined as the benefits that humans directly or indirectly receive from ecosystems [1]. The Millennium Ecosystem Assessment (MA) [2], first conducted in 2005, established a framework for the global assessment of ecosystems. This framework divides indicators of ESs into four categories: provisioning, regulating, cultural, and supporting services. These categories are based on the connection between ESs and human well-being. Due to regional socio-economic factors and the rapid, high-intensity expansion of land for human use, ecological systems are facing continuous destruction [3]. This has caused a significant reduction in the ecosystem service supply (ESS), while the demand for a better living environment continues to grow [4]. Consequently, this has intensified the incoordination relationship between ESS and ecosystem service demand (ESD), negatively affecting sustainable development and human well-being [5]. Thus, it will be difficult to effectively manage and optimize regional ecosystems and encourage sustainable development via only the supply of ESs, while ignoring the human demand for ESs. Investigating the relationship between ESS and ESD not only addresses the challenges of sustainability arising from increasing human demand, but also establishes a solid theoretical foundation for adept and efficient ES management and utilization practices [6,7]. Therefore, understanding the relationship between ESS and ESD is essential for effective land use planning and decision-making processes that ensure the long-term resilience of ecosystems and social well-being.
As a closer relationship is established between ESs and human well-being, ES demand (ESD) gradually integrates with ESS. ESS is generally defined as the beneficial effect that ecosystems have on society [8]. The structure and function of ecosystems often result in simultaneous positive and negative changes in ESS [9]. However, there is no universally accepted definition of ESD. Currently, it is mainly defined from two perspectives. From a consumption perspective, ESD refers to the services provided by ecosystems that are useful to consumers, reflecting the actual demand for ESs [10]. From a preference perspective, ESD refers to the ecosystem services requested by social groups [11], reflecting not only the actual demand for ESs, but also the potential demand that cannot be met due to certain conditions. ESD includes the requirements, desires, or aspirations of human societies in relation to the benefits and contributions provided by ecosystems [12]. These needs arise from the dependence of human well-being and quality of life on the services and resources provided by natural ecosystems. ESD is driven by various factors, including population growth, economic activities, urbanization, and societal preferences. When the demand for ESs exceeds the natural capacity of ecosystems to provide them, it can lead to the overexploitation, degradation, or loss of those services. Conversely, if the supply of ES surpasses the demand, it can result in the underutilization or inefficient allocation of resources. Therefore, the relationship between ESS and ESD is crucial for creating a dynamic balance through which ecosystem products and services are transferred from ecosystems to social systems [13,14].
In recent years, an increasing number of studies have focused on evaluating the relationship between ESS and ESD, considering the coordination or conflict between regional social systems and ecosystems [15,16]. Thus far, these studies have primarily concentrated on quantifying and comparing the supply and demand of distinct ES indicators, such as water supply [17], flood regulation [18], air purification [19], and erosion control [20], and identifying the degree of mismatch, imbalance or incoordination between ESS and ESD. The majority of these investigations are centered around specific types of ecosystems, including forests [21], croplands [22], and urban regions [23]. These studies were mainly conducted in European contexts; however, in recent years, an increasing number of studies have been conducted in China, driven by ecological regionalization policies at different scales. Most previous studies assessed ESS, mainly focusing on quantifying and identifying the patterns and functions of ESs, as well as the effects of land use and land cover changes on ESs [24]. Compared to ESS, there are two main distinct approaches to ESD assessment: One is to evaluate the demand for each individual ES, and the other involves conducting a comprehensive assessment of the overall demand for ecosystem services. The land development index (LDI) has been widely used to comprehensively assess the ESD [25]. The rationale behind using the LDI to assess demand for ESD lies in its capacity to gauge the intensity of land development and its corresponding impact on the demand for these services [26]. Thus, ESS and ESD can be expressed from the perspective of the interaction between ecosystems and social systems; hence, the coordination relationship between ESS and ESD is a significant reflection of whether regional socio-economic structure and natural ecological backgrounds can develop in a coordinated manner [27].
Previous studies mainly analyzed temporal dynamic variations and spatial imbalances or mismatches in ESS and ESD [28,29], using methods such as modeling, mapping, participatory methods, etc. [30]. The concept of “coupling coordination” provides a framework for quantitating the coordination between ESS and ES [31,32]. The coupling coordination degree (CCD) is a measure used to assess the level of coordination and interdependence between different components or subsystems within a larger system [33]. In the context of ES, the coupling coordination degree model (CCDM) can be used to evaluate the level of coordination between the ecological system that provides services and the social system that demands and utilizes such services. It can also be applied to assess how effectively the ESS and ESD function together and whether their interactions are balanced and mutually beneficial. In recent years, many studies have applied the CCDM to analyze the relationship between ESS and ESD. For example, Guan et al. [34] provided valuable insights into the evolving characteristics of the spatial coupling between ESS and ESD using the CCDM. Li et al. [35] analyzed the dynamic characteristics of the supply and demand coupling of ESs in Lanzhou, China. Yang et al. [36] identified the coupling coordination relationship between sustainable development and ESs in Shanxi Province, China. Therefore, CCD is regarded as an efficient tool for researchers and policymakers to gain insights into the functioning of the coupled ecological and social systems and identify areas that require attention or intervention. Consequently, this helps us to understand the complex interactions between systems as well as design strategies for the sustainable management and conservation of ESs.
The coordination between ESS and ESD is influenced by the rapid development of regional socio-economic factors and the ecological factors, especially when socio-economic factors interact with natural and ecological factors [37]. Natural factors, land use/land cover, and socio-economic factors have been identified as the primary determinants of the relationship between ESS and ESD. For instance, Sun et al. [38] conducted an empirical study examining the correlations among 12 natural and socio-economic variables related to both ESS and ESD. Their study sheds light on the disparity between these two aspects within the United States. Wu et al. [16] analyzed the relationships between ES supply and demand and identified the effect of forest area on ESS, as well as the effects of per capita GDP, energy consumption per unit of GDP, and permanent population on ES demand in China. Peng et al. [39] systematically analyzed the impact of urbanization on ESs in metropolitan areas. Previous studies identified multiple influencing factors and encompassed different ESs. However, the majority of studies primarily focus on identifying the individual influencing factors of ESS or ESD. Alternatively, some studies have solely examined the influencing factors of ESS and ESD as substitutes for the actual relationship between the two [40]. As a result, there is a lack of studies investigating factors that directly impact the relationship between ESS and ESD. Research on the relative importance of socio-ecological drivers of coordination between ESS and ESD remains limited, and little attention has been given to the relationships between various drivers and the coordination between ESS and ESD, as well as associated spatial influences [41,42].
Shanxi Province is a typical resource-based area heavily reliant on coal and other resources for rapid economic growth. However, this excessive consumption of resources has created an exceedingly fragile ecosystem across the province. This increase in ecosystem degradation is having a negative effect on the economy and society, posing significant challenges to the sustainable development of Shanxi Province [43]. For example, intensive agriculture and mining activities are contributing to soil degradation in Shanxi. Erosion, the loss of topsoil, and contamination from pollutants can have serious consequences for agricultural productivity and ecosystem health. Unsustainable mining practices, particularly in coal mining, are causing land subsidence issues in Shanxi [44,45]. This phenomenon can lead to infrastructure damage, waterlogging, and the disruption of ecosystems [46]. In addition, some parts of Shanxi are vulnerable to desertification due to factors such as soil erosion, overgrazing, and unsustainable land use practices. This threatens agricultural productivity and ecosystem stability [47]. Since 2000, large-scale ecological restoration projects, such as the Natural Forest Protection Projects (NFPP) and the Grain for Green Program (GFGP), have significantly bolstered vegetation restoration in Shanxi Province. Efforts to address these ecosystem problems likely involve a combination of policy interventions, stricter environmental regulations, technological innovations, public awareness campaigns, and sustainable development practices.
Based on this context, this study aims to analyze the coupling coordination relationship between ESS and ESD and identify the associated socio-ecological driving variables in Shanxi Province. The objectives of this study are to (i) quantify and map the spatial distribution of the ES supply and demand, respectively; (ii) analyze the spatial–temporal characteristics of CCD between ESS and ESD based on CCDM; and (iii) determine the decisive influencing factors and the effects of interactions between factors using GDM. The results are anticipated to provide valuable information for achieving a harmonious balance between economic development and ecological restoration on a national scale within the provincial administrative units of China.

2. Materials and Methods

2.1. Study Area

Shanxi Province, located in the northern part of China (110°14′–114°33′ E, 34°34′–40°44′ N), covers an area of 156,700 km2, accounting for 1.6% of the country’s territory (Figure 1a). It consists of 107 counties and is characterized by a typical mountain plateau terrain. Its topography is complex and diverse, including mountains, hills, plateaus, basins, and platforms. Mountains and hills make up 80% of its area, with altitudes ranging from 208 m to 2988 m above sea level (Figure 1b). The study area falls within the temperate continental monsoon climate region, the annual average temperature is between 4 °C and 14 °C, and the average annual rainfall is 468 mm. The dominant land use/land cover types include cropland, grassland, forest, and construction land (Figure 1c). Shanxi Province straddles the Yellow River basin and the Haihe River basin, and the river system is a self-generated outflow. The total population of the study area accounts for 2.48% of the national population and 1.71% of China’s total GDP [48].
Shanxi Province is one of the most important provinces with coal and mineral resources in China. As a region highly dependent on coal and mineral resources, its economic and social development is closely linked to its ecological environment [49]. However, excessive coal and mineral resource exploitation, the overexploitation of groundwater, rapid urbanization, and accelerated water and land resource exploitation have made the ecosystem in the study area extremely fragile. In 2019, the Chinese Government put forth the objectives of “ecological protection and high-quality development” for the Yellow River basin. As part of the Yellow River basin, in recent decades, Shanxi Province has faced daunting challenges in coordinating population, resources, ecosystems, and economic development.

2.2. Data Collection

In this study, we selected 2000 and 2020 as representative years to collect and analyze data. The data encompass both spatial and statistical information. The spatial data were processed at a grid cell resolution of 1 km × 1 km, while statistical data were aggregated at the county level. Table 1 provides a comprehensive list of the primary data required to calculate the ES indicators. To ensure consistency, all spatial data were transformed to a common spatial reference system, specifically the WGS84 coordinate system and Albers equal-area conic projection. The flowchart depicted in Figure 2 illustrates the methodology that we employed to achieve our study objectives.

2.3. Quantifying ES Supply

2.3.1. Selection of ES Indicators

The selection of appropriate ES indicators is a critical step in assessing ecosystem services. In this study, we followed certain criteria for selecting ES indicators: (i) aligning with the classification of ES according to the Millennium Ecosystem Assessment [2] to ensure comparability with other studies; (ii) considering existing case studies conducted in Shanxi Province and selecting ES indicators that are closely related to the natural, ecological, social, and economic conditions of the study area; and (iii) the availability of the primary data required for evaluating ES indicators. Based on these criteria, our study focused on five key ES indicators relevant to the study area. These included one provisioning service (crop production), three regulating services (water retention, soil conservation, and carbon sequestration), and one cultural service (outdoor recreation).

2.3.2. Calculation of the ES Indicators

To quantify the selected ES indicators, we employed existing and widely used assessment models originally developed for this purpose. Specifically, the assessment of crop production was based on annual crop yield data [50]. The water balance equation served as a proxy for measuring water retention. The Universal Soil Loss Equation (USLE) model was used to calculate soil conservation [51]. Net primary production (NPP) was used as a proxy for carbon sequestration [52,53] and was assessed using the Carnegie–Ames–Stanford Approach (CASA) model, a widely adopted approach for NPP estimation [54,55]. The spatial distribution of individual ES indicators was visualized via mapping in ArcGIS.
To analyze the relationships between the five ES indicators, we employed ArcGIS 10.2, which serves as a common spatial unit, to aggregate all ES indicators at a national level. According to Raudsepp-Hearne et al. [56], administrative boundaries are suitable for identifying socio-ecological systems in a landscape, as management decisions at this level influence the provision and consumption of ES. The specific models and processes used for assessing the ES indicators are summarized in Table 2.

2.3.3. Assessment of ES Supply Index

Since each ES has its own measurement unit, we individually standardized the ES values and then summarized them within each county to mitigate the influence of magnitude and variability. We employed min–max normalization to standardize the values of the five ES indicators [63,64]. This normalization method removes the units of the input data and scales them to a common range. After standardization, the standardized values were accumulated to obtain the ES supply index (ESSI), which represents the total ecosystem service supply. The calculation equation for ESSI is as follows:
E S S i j = E S i j m i n E S j m a x E S j m i n E S j
E S S I i = j = 1 n E S S i j
where ESSIi is the ES supply index of county i; ESSij is the standardized value for ES j of county i; ESij is the initial value for ES j of county i; max ESj is the maximum value of ES over 107 counties, and min ESj denotes the minimum value of ESj; n = 5.

2.4. Quantifying ES Demand

The ES demand represents the human demand and preference for ecosystem products and services within a specific time period. In this study, we used a research method [26,65] to quantify ESD by considering land development intensity, population density, and gross domestic product (GDP) per area. This helped us to understand the coordination between the development and preservation of essential natural processes that sustain human well-being and environmental quality [29,66,67]. Specifically, land development intensity was measured as the percentage of construction land in the total land area. It reflects the intensity of the human consumption of ES. A higher percentage of construction land indicates a greater intensity of human land development in a given area, and consequently, a higher demand for ES. Population density serves as an indicator of the amount of ES demand. A higher population density corresponds to a greater ESD. GDP per area reflects the economic development of the region and indirectly indicates how much humans wish to consume or utilize ESs. Logarithmic methods were employed to remove fluctuations in the data. The ESD index is calculated using the following formula:
E S D I i = D i × l g P i × l g G i
where ESDIi is the ES demand index of each county i; Di, Pi, and Gi are the land development intensity (%), population density (person/km2), and GDP per area (yuan/km2) of county i, respectively.

2.5. Assessing Coordination between ES Supply and Demand

In this study, we employed the CCDM to investigate the interactive coordination relationship between ESS and ESD. The CCDM highlights the interdependence of ecological and social systems and aims to understand how they interact and mutually influence each other. By adopting this model, we can develop a holistic understanding of the relationship between ecosystems and human societies.
The CCDM recognizes the interconnectedness of ecological and social systems and emphasizes the importance of studying them together. This enables us to analyze the developmental pattern of these systems or indicators, progressing from disorder to order [68]. This representation provides insights into the overall effectiveness and synergistic impact between systems [69].
Mathematically, the CCDM is expressed as follows:
C C D i = C i · T i
C i = S E S S C I i · D E S S D I i / S E S S C I i + D E S S D I i / 2 2 1 / 2
T i = α · S E S S C I i + β · D E S S D I i
where CCDi represents the coupling coordination degree of county i (0 ≤ CCDi ≤ 1) between ESS and the ESSD; Ci refers to the coupling degree between ESS and the ESSD; Ti is the comprehensive development index of ESS and the ESSD; and SESSIi and DESDIi are the values of standardized ESSI and ESDI (0 ≤ SESSIi ≤ 1, 0 ≤ DESDIi ≤ 1). α and β are the weights to be determined; due to the equal importance of ESS and the ESSD in the coordination, α and β are given the same weight, that is, α = β = 0.5. Referring to previous research [32], we divide the CCD into five levels: When 0 ≤ CCDi ≤ 0.20, the ESS and ESSD are in extreme incoordination; when 0.20 < CCDi ≤ 0.35, they are in moderate incoordination; when 0.35 < CCDi ≤ 0.55, they are in reluctant coordination; when 0.55 < CCDi ≤ 0.70, they are in moderate coordination; and when 0.70 < CCDi ≤ 1, they are in superior coordination.

2.6. Driving Variables of ES Coordination

2.6.1. Critical Driving Variables

In this study, we selected socio-ecological variables to explain the spatiotemporal differences between ESS and ESD, based on relevant research [70,71]. Potential explanatory variables were chosen from three sources: (1) the variables used to quantify ESS or ESD in our study, (2) variables identified in the literature as directly or indirectly driving individual ESs and/or their associations [72], and (3) variables for which quantitative data were available. After considering these factors, we ultimately selected thirteen potential socio-ecological variables, including natural variables such as elevation (DEM), slope (SLOPE), average annual precipitation (PRE), and average annual temperature (TEM); ecological variables such as NDVI, percentage of crop land (CROP), percentage of forestland (FOREST), and percentage of grassland (GRASS); and socio-economic variables such as percentage of construction land (CON), total population (POP), GDP, proportion of urban population (URBAN) and distance to the nearest county center (COUNTY) (Table 3).

2.6.2. Effects of Driving Variables on Coordination via Geo-Detector Model

The geo-detector model (GDM) can help identify the most influential factors or variables that contribute to specific spatial patterns, and reveal how different factors interact in a spatial context. It offers insights into the potential effects of human activities on ecosystems, water resources, or air quality. In this study, we employed the GDM to assess the spatial correlation between the explanatory variables and the dependent variables through spatial variance analysis (SVA) [73,74]. The GDM is a valuable analytical tool to identify and quantify the spatial associations between driving factors and specific outcomes. This insight can guide decision making by highlighting where interventions or ecological resource allocation should be focused for maximum impact. By employing statistical and spatial analytical techniques, the GDM enables researchers to identify dominant driving factors and their interactive effects, as well as explore spatial patterns and trends in complex geographical processes [75]. The fundamental assumption of the GDM is as follows: if an independent variable X significantly affects a dependent variable Y, then the spatial distributions of X and Y should exhibit similarity. SVA is used to compare the spatial consistency between the dependent variable and independent variables. Based on this comparison, the interpretation of independent variables in relation to the dependent variable can be quantified.
In this study, we utilized the “factor detector” module of GDM to identify the driving factor(s) that determine the distribution of CCD. This module identifies the extent to which the driving variables explain the spatial differentiation of CCD. The calculation results of the factor detector include the q-statistic and p-value. The q-statistic represents the influencing coefficient of the driving variable on CCD, with larger values indicating a stronger impact of the driving variable on CCD. The p-value indicates the significance level of the explanation, and a significance level of 0.1 (p-value < 0.1) is considered statistically significant. The formula for the factor detector is as follows:
q = 1 h = 1 l N h σ h 2 N σ 2
where q signifies the influencing coefficients of the driving variables for the ES (q -statistic), the values of which range from 0 to 1, where 0 corresponds to no correlation between the two and 1 to CCD’s complete dependence on a driving variable. σ2 is the variance of the CCD, and N is the size of CCD. The superposition of the driving variables and CCD forms L layers in CCD, which are indexed by h = 1, 2… l, and Nh and σ h 2 represent the scale and variance of layer h, respectively.
The “interaction detector” module of GDM was used to examine whether two factors have a stronger or weaker effect on ESs than they do independently. The types of interactions between two variables are as follows:
Enhance: if q (D1 ∩ D2) > q (D1) or q (D2)
Enhance, bivariate: if q (D1 ∩ D2) > q (D1) and q (D2)
Enhance, nonlinear: if q (D1 ∩ D2) > q (D1) + q (D2)
Weaken: if q (D1 ∩ D2) < q (D1) + q (D2)
Weaken, univariate: q (D1 ∩ D2) < q (D1) or q (D2)
Weaken, nonlinear: if q (D1 ∩ D2) < q (D1) and q (D2)
Independent: if q (D1 ∩ D2) = q (D1) + q (D2)
where the symbol “∩” denotes the intersection between the layers D1 and D2. The attributes of layer (D1 ∩ D2) are determined by the combination of the attributes of layer D1 and D2 using a spatial overlay to form a new layer. q (D1), q (D2), and q (D1 ∩ D2) were calculated using Equation (1). By comparing the sum (q (D1) + q (D2)) of the factors’ contribution to two individual attributes (q (D1), q (D2)) with the contribution of the two attributes when combined (q (D1 ∩ D2)), the interactive effects of the two factors can be defined using the above seven types.

3. Results

3.1. Spatial–Temporal Patterns of ES Supply and Demand

Using a quantitative method, we calculated ESSI and ESDI data for the 107 counties of Shanxi Province for the reference years 2000 and 2020. We normalized the two indices to a range of 0–1 and mapped them to facilitate comparisons (Figure 3 and Figure 4).
The unique geographical location of Shanxi Province resulted in significant variations in natural conditions, leading to considerable heterogeneity in the counties’ ability to provide ESs. The ESSI exhibited substantial variation across Shanxi Province in both 2000 and 2020, which was similar to the distribution characteristics of outdoor recreation (Figure 3). Areas with high ESS were dispersed across the mountainous regions of the study area, surrounded by counties with higher ESSI. Specifically, Mount Taiyue, Mount Zhongtiao, and Mount Wangwu, located in the southern parts of the province, have a particularly abundant supply of ESs. The areas with low supply were concentrated in the northwest edge of Shanxi and the western part of Mount Luliang. From 2000 to 2020, the spatial pattern of ESSI remained largely unchanged in most counties. However, there was a slight overall decrease in the level of ESS, and the ESSI of Yuanping and Xinzhou, located in the northwestern part of Shanxi Province, experienced significant increases. There was a clear upward trend in the number of counties with the lowest supply, increasing from 20.6% for the 107 counties in 2000 to 32.7% in 2020. However, the number of counties with the highest supply remained static.
The spatial distribution hierarchy of ESDI was weaker compared to ESSI (Figure 4). In both 2000 and 2020, ESDI displayed spatial distribution characteristics with higher values in the central areas and lower values in the outer areas of the study area, which closely resemble the spatial distribution of construction land. The areas with the highest and higher ESS grades were primarily concentrated in central Shanxi, southeastern areas, and the northern plains, indicating a relatively concentrated distribution. The areas with the lowest demand were contiguous and distributed in mountainous regions such as Mount Luliang in the west and Mount Taihang in the east. From 2000 to 2020, due to population growth and economic development in the study area, there was a clear increase in ESDI. The number of counties with higher grades of ESS increased, while the counties with the lowest demand decreased. Some areas with medium demand in 2000 shifted to higher demand categories by 2020, indicating a transformation from relatively lower to higher grades. Overall, the spatial distribution of ESDI displayed noticeable differences between the outskirts and the middle regions of the study area.

3.2. Coupling Coordination Characteristics of ES Supply and Demand

Using the CCDM, we measured and mapped the CCD of ESSI and ESDI for the 107 counties of Shanxi Province in 2000 and 2020 (Figure 5). The CCD values ranged from 0 to 0.57 in 2000 and from 0 to 0.43 in 2020 (Table 4), indicating a relatively low level of coupling coordination.
There are four main types of coupling coordination relationships between ESSI and ESDI: extreme incoordination, moderate incoordination, reluctant coordination, and moderate coordination (Figure 5). Most counties belong to the extreme incoordination and moderate incoordination patterns, mainly located in the western and eastern parts of Shanxi. In 2000, the areas with reluctant coordination were primarily found in the Taiyuan Basin in central Shanxi, Linfen Basin, and the Yuncheng Basin in the southeastern parts. However, in 2020, reluctant coordination areas were sparsely distributed in only a few municipalities. Overall, there was a clear incoordination relationship between ES supply and demand in Shanxi.
Between 2000 and 2020, the number of counties with extreme incoordination significantly increased from 29.0% to 36.5%. The number of reluctant coordination areas decreased from 34.6% to 24.3%, with fourteen counties transitioning from reluctant coordination to moderate incoordination. This indicated a substantial decline in the coupling coordination between ESSI and ESDI in Shanxi over the 20-year period. The main change characteristics were the negative transitions from relatively high coordination grades to incoordination grades. The most significant changes were observed in the shifts from moderate incoordination to extreme incoordination and from reluctant coordination to moderate incoordination, accounting for 20.31% and 16.28% of the counties, respectively.
In general, the analysis of coordination between ES supply and demand revealed coexisting states of coordination and incoordination, with incoordination being predominant in most counties. The difference in the spatial polarization of CCD in 2020 was more significant than in 2000. Overall, the coordination relationship between ES supply and demand in Shanxi Province deteriorated between 2000 and 2020.

3.3. Determining Drivers for the Coupling Coordination Degree between ES Supply and Demand

The GDM was utilized to identify the most influential socio-ecological drivers for the coupling coordination between ES supply and demand. The factor detection results for thirteen socio-ecological variables yielded the influencing coefficients (q-statistic values) and significance levels (p-values) (Figure 6).
In 2000 and 2020, the variables SLOPE, PRE, and COUNTY did not pass the significance tests. Meanwhile, the variable TEM was a statistically significant driving variable (p-value < 0.001) in 2000, but its significance level was greater than 0.1 in 2020. Conversely, the variable URBAN was not significant in 2000 (p-value > 0.1) but became a significant driving variable in 2020 (p-value < 0.05). Figure 6 presents the sorting results for the q-statistic values for significant variables (p-value ≤ 0.05), revealing their influencing coefficients on CCD in 2000 and 2020. The influencing coefficients of CON, POP, and GDP were greater than 0.5 and higher than the other factors in both 2000 and 2020. This demonstrates that these variables were the main drivers of the spatial pattern of CCD during both time periods. GRASS, CROP, and DEM were considered sub-high determinate variables based on their q-statistic values. Notably, the variables TEM, NDVI, and FOREST had weak effects on CCD in 2000, whereas in 2020, the variables URBAN, FOREST, and NDVI had a weak effect. This suggests that the impact of forest coverage on the coordination relationship was relatively small in the study area.
Furthermore, the interaction detector module of GDM assessed the influencing coefficients of any two socio-ecological variables (Figure 7) and compared them with their separate influencing coefficients (q-statistic values). The results reveal two interaction modes of socio-ecological variables on CCD: nonlinear enhancement and mutual enhancement. This indicates that the explanatory power of the interaction between any two variables for CCD is greater than that of any single variable. The interactive effects on CCD between CON and the other variables were the strongest, with q-statistic values exceeding 0.85. After interacting with POP and GDP in both 2000 and 2020, the q-statistic values of all variables were above 0.65, which was higher than the values for the separate effects of POP and GDP on CCD. It is important to note that SLOPE, PRE, and COUNTY did not exhibit significant effects on CCD in the single-factor detection results (Figure 6). However, after interacting with CON, POP, and GDP, the influencing coefficients of SLOPE, PRE, and COUNTY were 0.88, 0.67, and 0.68 in 2000, respectively, and similar results were observed in 2020. This suggests that even if individual socio-economic variables do not have a significant effect on the spatial distribution of CCD, they may play a key role via interaction with variables that have high influencing coefficients. Thus, SLOPE, PRE, and COUNTY were identified as important external driving factors for CCD between ESSI and ESDI. Additionally, this highlights the importance of considering the interactions between variables in understanding the spatial distribution of CCD.

4. Discussion

4.1. Spatial–Temporal Characteristics of ESSI and ESDI

Our framework of analysis was used to make sense of relationships between ES supply (ESS) and ES demand (ESD) in a complex social–ecological system. Utilizing the ES supply index (ESSI), we have derived the spatiotemporal variation characteristics of ESS (Figure 3). These characteristics provide valuable insights into assessing the level of ES supply, as well as tracking its enhancement or decline in every county within Shanxi Province.
We found that high values of ESSI were predominantly observed in mountainous areas with dense forest coverage. In contrast, low values were mainly found in the Fen River basin, which is dominated by cropland and stretches from northeast to southwest. The Fen River basin is the most active area in terms of social and economic activities in Shanxi Province. The extensive human disturbance in this basin, characterized by high population density and economic growth, has led to the rapid expansion of construction land, causing a reduction in ecological land. As a result, the supply of ESs in the counties within this watershed is relatively low. Similarly, the northwestern parts of the study area, dominated by cropland and grassland, also exhibited relatively low ESSI values. This aligns with findings from previous studies [76]. This region represents a transition zone from scrub steppe to typical steppe, and factors such as increasing population, livestock, and desertification posed significant threats to the area for an extended period. Despite the implementation of the GFGP in these northwestern counties of Shanxi Province since 2000, the supply of local ESs continues to deteriorate [77]. Hence, the grass-planting-based GFGP in this region is required to improve efficiency and achieve ecological protection. Additionally, the western area of Mount Luliang exhibited relatively low ESSI values. This can be attributed to Mount Luliang’s location in the loess hilly and gully region, which is characterized by severe soil and water loss [78]. Moreover, these mountainous areas are contiguous to regions with concentrated mineral resources and frequent human activities. As of 2021, there are 91 coal mines in Luliang, accounting for 13.62% of the total number of coal mines in Shanxi. Extensive mining methods and inadequate management measures have imposed a burden on ecosystems, leading to a diminished supply level of regional ESs.
Concerning the ES demand index (ESDI), we carefully curated three indicators representing land use, population dynamics, and economic factors (Figure 4). These indicators collectively offer a comprehensive representation of ESD, encompassing preferences and requirements. We found that high values of ESDI were primarily distributed in the central Taiyuan Basin, northern Datong Basin, southern Linfen Basin and Yuncheng Basin, and southeastern Changzhi Basin, while low values were observed in the mountainous areas on the east and west sides. These basin regions, with Taiyuan, Datong, Linfen, Yuncheng, and Changzhi as their central cities, have experienced significant population and industrial concentration, resulting in high demand for ESs. Notably, the spatial pattern of ESDI closely resembles that of the degree of land use development (Figure 3). This finding is consistent with existing studies [79] that indicate a strong correlation between ESDI distribution and land use development degree, population density, and per capita GDP. It is a common phenomenon in many parts of China to expand the area of construction land in order to meet social and economic needs, particularly in areas with extensive human activity [16]. From 2000 to 2020, the average ESDI values significantly increased in the Datong Basin, decreased in the Taiyuan Basin and Yuncheng Basin, and the disparity in ESD values among the counties within these basins significantly decreased. Additionally, ESDI levels in Taiyuan city and Changzhi city district, which had the highest ESDI values in the whole study area, remained relatively stable between 2000 and 2020. This stability can be attributed to the slower economic and population growth in Taiyuan and Changzhi, reaching a stable socio-economic agglomeration state within Shanxi Province [16]. Overall, Shanxi Province exhibits clear spatial mismatch characteristics between ESS and ESD.

4.2. Spatial–Temporal Characteristics of CCD between ESSI and ESDI

The coordination mechanism between ESS and ESD primarily revolves around the harmonization of the ecosystem and social system [80]. It involves a feedback loop, where societal choices, such as land use and economic activities, impact the provision of ESs. In turn, the condition and health of the ecosystem shape the quality and quantity of services that can be supplied to meet social needs. Thus, it is rare to gain a clear spatial relationship between the biophysical supply of ESs and their demand; it is a rare achievement. This rarity primarily stems from the utilization of distinct measurement units for assessing supply and demand [78]. Furthermore, this situation underscores a current challenge within this field of study. In this study, we introduced the CCD model to address this challenge, and obtained the spatiotemporal variation characteristics of the coordination relationship between ESS and ESD in various counties of the Shanxi Province.
The majority of counties exhibited a state of incoordination in terms of CCD (Figure 5). Extreme incoordination was predominantly observed in Mount Taihang in the east, Mount Luliang in the west, and Mount Taiyue in the south of Shanxi Province. This can be attributed to the higher supply of the five key ESs compared to the demand in these counties (Figure 3 and Figure 4). Most of these counties are situated in mountainous areas characterized by forest and grassland with high vegetation coverage [36]. Due to the implementation of ecological protection projects, there is minimal interference from human activities in these counties, resulting in little change in the type and quantity of land use. Furthermore, NFPP and GFGP contribute to increased vegetation coverage in these mountainous areas, enhancing the types and capacities of ESs. As a result, a high level of “lock-in effect” of regional ESS has been achieved [35,76], and ESS in these mountains is significantly weaker compared to the plain and basin areas. On the other hand, moderate incoordination was primarily observed in the Datong Basin, Xinding Basin, Taiyuan Basin, Linfen Basin, Yuncheng Basin, and Changzhi Basin, spanning from the north to south of the study area. This can be attributed to the extensive human activity in these basin regions and indicates the need for additional efforts to improve the coordination between ESS and ESD in these areas.
The coordination relationship between ESSI and ESDI exhibited a decrease in incoordination between 2000 and 2020, as evidenced by an increase in the number of extreme incoordination counties and a decrease in reluctant coordination and moderate coordination counties. For example, the relationship between ESSI and ESDI in Taiyuan, Linfen, Yuncheng, and Changzhi cities shifted from moderate coordination in 2000 to reluctant coordination in 2020. No other counties changed to a moderate coordination, resulting in the absence of moderate coordination counties/districts in the study area in 2020. This shift can be attributed to the slower growth rate of ESS compared to ESD in these cities. The increase in population, expansion of construction land, and economic growth have led to higher ESD in these districts [81]. Despite the increase in green spaces within the cities between 2000 and 2020 due to China’s ecological civilization construction projects, it falls short of meeting the substantial demand for ES in urban production and daily life [82].
Overall, the distribution of CCD in 2020 exhibited noticeable spatial differences between the basin region and the mountain region compared to 2000. The coordination relationship between ESSI and ESDI demonstrates incoordination and spatially varies across Shanxi Province. The degree of incoordination intensified in nearly half of the counties from 2000 to 2020. In 2020, the CCD highlighted the spatial disparities between urban areas, agricultural areas, and forest–grassland areas from the perspective of land use, as well as the differences between valleys, basins, and mountains from the perspective of terrain.

4.3. Associations between CCD and Driving Covariates

In our study, we employed a geo-detector model (GDM) to capture the spatially response characteristics of the dependent variable in relation to the independent variable. This approach allowed us to analyze how different factors contribute to the observed spatial patterns of the coordination between ESS and ESD across the study area (Figure 6 and Figure 7). We found that socio-economic factors had a greater impact on the coordination relationship between ESS and ESD than natural and ecological factors, emphasizing the significance of socio-economic factors in shaping the spatial pattern of this relationship. This finding aligns with previous research that highlighted the role of socio-economic variables as determinants of ES supply and demand distribution [43,83]. The dominant factors influencing the spatial pattern of CCD were identified as construction land, followed by population and GDP. The expansion of construction land, driven by extensive human activity, has significantly influenced land use patterns in terms of magnitude, type, and distribution [84]. The complex nature of human activity further complicates these relationships. However, these findings contrast with a study by Yang et al. [82] in China’s Loess Plateau, which found that vegetation cover had the greatest positive effect on the relationship between ES supply and demand. This discrepancy could be attributed to the dominant influence of vegetation coverage on both ESS and ESD in the Loess Plateau.
Interestingly, the effect of grassland on CCD was second only to socio-economic factors, ranking below elevation and cropland. Grassland emerged as the vegetation cover factor with the strongest impact on CCD, indicating its significant role in shaping the spatial pattern of the coordination relationship between ESS and ESD in Shanxi Province. This finding is supported by previous studies that emphasized the influence of land use changes caused by socio-economic factors on ESs [85]. Notably, the effect of forestland and NDVI on CCD was much smaller compared to grassland, suggesting that the increase in grassland area through GFGP mainly influenced the relationship between ES supply and demand in Shanxi Province. According to data from the 2020 “Shanxi Province Third Land Survey Main Data Bulletin”, grassland covers an area of 3.11 million hectares in Shanxi Province, and is primarily distributed in Datong, Xinzhou, Luliang, Jinzhong, and Linfen cities, accounting for 73% of the province’s grassland [86]. However, the grassland ecosystem in the study area remains fragile, with 70% of the grassland experiencing varying degrees of degradation and facing challenges such as insufficient protection and restoration, low utilization and management efficiency, and a lack of effective technological support. Therefore, the coordination degree between ESS and ESD could be improved by implementing quantitative and spatial adjustments to grassland planting policies.
In this study, the slope and precipitation factors were found to have no significant direct effect on CCD. However, after interacting with construction land, population, and GDP, slope and precipitation played an important role in shaping the spatial pattern of CCD. This indicates that socio-economic factors enhance the influence of natural factors on CCD. Generally, slope and precipitation influence the supply capacity and demand preferences for ESs by controlling the spatial distribution of human activity and landscapes [76], and previous studies demonstrated that precipitation facilitates coordination between ESS and ESD in arid and semi-arid areas [79]. The changes in terrain and precipitation affect both social and economic processes, leading to changes in coordination relationships. Previous research provided guidance on identifying factors that contribute to ES supply and demand mismatch [87], ESS, and ESD, including terrain ruggedness for supply and population density for demand [88].
The spatial distribution of the coordination relationship between ESS and ESD is primarily supported by the expansion of construction land, population growth, and economic benefits. However, this puts immense pressure on ES supply. Although forest and grassland coverage significantly increased between 2000 and 2020, the expansion of construction land has outpaced these gains. The growth rates of forestland and grassland in the study area were 4.2% and 4.3%, respectively, while the growth rate of construction land was as high as 105.4%. Since 2000, Shanxi Province has been the subject of economic system reforms, and thus rapid economic development, but the lack of awareness regarding ecological protection has resulted in the degradation of ESs. Consequently, the coordination relationship between ES supply and demand deteriorated somewhat between 2000 and 2020, indicating that the expansion of construction land and the concentration of population and industry have threatened the coordination between ESS and ESD in the counties of Shanxi. To address this issue, policymakers in Shanxi Province must make significant progress in promoting the coordination between ESS and ESD, implementing measures to ensure a dynamic balance between ES supply and demand.

4.4. Limitations

There are several limitations and uncertainties in this study. Firstly, the evaluation of five types of ES was limited in terms of reflecting the overall ESS level due to data availability and quality constraints [89]. In addition, the equal-weight superposition calculation method used to calculate total ESS may have overlooked significance of various ES types in Shanxi Province. Future studies should strive to include a more comprehensive range of ES indicators and consider their relative importance. Secondly, the selection of indicators to represent ESD focused on land development intensity, population density, and GDP per area, assuming that all types of ESs have the same demand. This oversimplification may not accurately reflect spatiotemporal changes in ESD. Further research is needed to refine the measurement of ESD and capture its dynamics. Thirdly, this study compared the changes in ESS and ESD between 2000 and 2020, overlooking the temporal volatility of ES. Incorporating temporal dynamics would enhance our understanding of the coordination relationship between ES supply and demand.
The use of county-level administrative boundaries was advantageous [90], providing official statistical data; however, this limits our understanding of causality and spatial heterogeneity in coordination relationships [91]. Shanxi Province has diverse topography, vegetation and natural geographical features, and there is a significant spatial mismatch between ESS and demand. Thus, it is difficult to reveal the variation of ESS at a local scale in the spatial units of county. In addition, spatial heterogeneity and scale effects impact the relationship between ESS and ESD [92]. Focusing only on a single scale tends to miss information about the correlation between scales, and the influence mechanism is inevitably one-sided. Future research should use a multi-scale analysis to capture the correlation between scales and comprehensively investigate influence mechanisms.
Additionally, our study assumed that ESD in a county is solely provided by the local ecosystem without considering the flow of ecosystem services across county boundaries. For example, water resources can originate from upstream regions, and food shortages in a city can be mitigated via food trade [93]. Future studies should account for the cross-boundary flow of ES, as services from neighboring ecosystems can also contribute to ESD. Furthermore, while GDM helped identify the strength of influence of socio-ecological factors on the spatial pattern of CCD, it was unable to capture whether this effect was positive or negative. Exploring the positive or negative nature of these effects would provide a more comprehensive understanding of the coordination relationship between ESS and ESD.

5. Conclusions

Shanxi Province has experienced a stage of rapid expansion of construction land; there has been a rapid transformation in the intensity, type, and pattern of land use, which has created an urgent need to optimize ESs, social progress, and economic development. Under this background, this study analyzed spatiotemporal changes in ES supply, demand, and their coordination relationship, and identified the relevant socio-ecological driving factors across 107 counties in Shanxi from 2000 to 2020.
The results reveal that the spatial pattern of ESS was closely linked to forest coverage, while ESD was closely related to the degree of land use development. The changes in ESS and ESD exhibited spatial heterogeneity. Over the study period, Shanxi Province experienced a slight decrease in ESS and an increase in ESD. The evaluation using the CCDM demonstrated significant incoordination between ESS and ESD in Shanxi, which worsened between 2000 and 2020. Based on the CCD of ESS and ESD, Shanxi Province was divided into four zones: extreme incoordination, moderate incoordination, reluctant coordination, and moderate coordination. The extreme incoordination zone was mainly located in mountainous regions, where ecosystems dominated by the eastern Taihang Mountain and western Luliang Mountain provided high levels of ESs but low ESD. The moderate incoordination zone was primarily found in basins with intensive human activity, where ecosystems dominated by cropland exhibited slightly higher levels of ESs compared to ESD. The reluctant coordination and moderate coordination zones were mainly located in central cities within basins, where the CCD between ESS and ESD significantly decreased. The spatial distribution of CCD was primarily influenced by construction land, population, and GDP, with grassland playing a secondary role, largely driven by the GFGP policy.
Accordingly, in order to promote the sustainable development of the counties, we propose the following recommendations: Firstly, the government should enhance ecological compensation policies for residents in mountainous regions, with a special focus on areas like Taihang Mountain and Luliang Mountain. This entails raising compensation for ecosystem service providers. Secondly, county and municipal district administrations should adopt advanced strategies for balanced ecological, social, and economic development. This will bolster land use efficiency and curb haphazard expansion of construction projects. Thirdly, policymakers and governments need to comprehensively assess the distinct impacts of various drivers on the interplay between ESS and ESD at both the local and county scales. This is particularly pertinent when formulating strategies for regional management approaches. These policy insights are applicable not just to Shanxi Province, but also to other regions endowed with coal and mineral resources. In future research, it is imperative to prioritize effective land management decisions that account for the interrelation between socio-ecological factors and the supply and demand of ecosystem services.

Author Contributions

Y.Z. and B.L. designed the study, collected data and conducted the research. B.L. analyzed data and participated in data processing and model calculation. Y.Z. wrote the first draft, and B.L. revised and edited the first draft. R.S. collected and processed the data. The remaining authors contributed to the discussion of results and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Shanxi Province Basic Research Plan Project under grant No. 202203021212496, and the Science and Technology Innovation Project of University in Shanxi Province under grant No. 2020L0248. We thank the academic editors and anonymous reviewers for their kind suggestions and valuable comments.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Costanza, R.; D’Arge, R.; Groot, R.D.; 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. Ecol. Econ. 1997, 25, 3–15. [Google Scholar] [CrossRef]
  2. MA (Millennium Ecosystem Assessment). Overview of the Millennium Ecosystem Assessment. 2005. Available online: http://www.millenniumassessment.org/en/About.html# (accessed on 7 July 2017).
  3. Chen, W.X.; Chi, G.Q.; Li, J.F. 2019. The spatial association of ecosystem services with land use and land cover change at the County level in China, 1995–2015. Sci. Total Environ. 2019, 669, 459–470. [Google Scholar] [CrossRef]
  4. Sutton, P.C.; Costanza, R. Global estimates of market and non-market values derived from nighttime satellite imagery, land cover, and ecosystem service valuation. Ecol. Econ. 2002, 41, 509–527. [Google Scholar] [CrossRef]
  5. Metzger, M.J.; Schröter, D. Towards a spatially explicit and quantitative vulnerability assessment of environmental change in Europe. Reg. Environ. Change 2006, 6, 201–216. [Google Scholar] [CrossRef]
  6. Zagonari, F. Using ecosystem services in decision-making to support sustainable development: Critiques, model development, a case study, and perspectives. Sci. Total Environ. 2016, 548–549, 25–32. [Google Scholar] [CrossRef] [PubMed]
  7. Cortinovis, C.; Geneletti, D. A performance-based planning approach integrating supply and demand of urban ecosystem services. Landsc. Urban Plan. 2020, 201, 103842. [Google Scholar] [CrossRef]
  8. De Groot, R.S.; Alkemade, R.; Braat, L.; Hein, L.; Willemen, L. Challenges in integrating the concept of ecosystem services and values in landscape planning, management and decision making. Ecol. Complex. 2010, 7, 260–272. [Google Scholar] [CrossRef]
  9. Bennett, E.M.; Peterson, G.D.; Gordon, L.J. Understanding relationships among multiple ecosystem services. Ecol. Lett. 2010, 12, 1394–1404. [Google Scholar] [CrossRef]
  10. Burkhard, B.; Kroll, F.; Nedkov, S.; Muller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic. 2012, 21, 17–29. [Google Scholar] [CrossRef]
  11. Villamagna, A.M.; Angermeier, P.L.; Bennett, E.M. Capacity, pressure, demand, and flow: A conceptual framework for analyzing ecosystem service provision and delivery. Ecol. Complex. 2013, 15, 114–121. [Google Scholar] [CrossRef]
  12. Wolff, S.; Schulp, C.J.E.; Verburg, P.H. Mapping ecosystem services demand: A review of current research and future perspectives. Ecol. Indic. 2015, 55, 159–171. [Google Scholar] [CrossRef]
  13. Castillo-Eguskitza, N.; Martín-López, B.; Onaindia, M. A comprehensive assessment of ecosystem services: Integrating supply, demand and interest in the Urdaibai Biosphere Reserve. Ecol. Indic. 2018, 93, 1176–1189. [Google Scholar] [CrossRef]
  14. Fu, B.J.; Tao, T.; Liu, Y.X.; Zhao, W.W. New developments and perspectives in physical geography in China. Chin. Geogr. Sci. 2019, 29, 363–371. [Google Scholar] [CrossRef] [Green Version]
  15. Lorilla, R.S.; Kalogirou, S.; Poirazidis, K.; Kefalas, G. Identifying spatial mismatches between the supply and demand of ecosystem services to achieve a sustainable management regime in the Ionian Islands (Western Greece). Land Use Policy 2019, 88, 104171. [Google Scholar] [CrossRef]
  16. Wu, X.; Liu, S.; Zhao, S.; Hou, X.; Xu, J.; Dong, S.; Liu, G. Quantification and driving force analysis of ecosystem services supply, demand and balance in China. Sci. Total Environ. 2019, 652, 1375–1386. [Google Scholar] [CrossRef]
  17. Quintas-Soriano, C.; Castro, A.J.; García-Llorente, M.; Cabello, J.; Castro, H. From supply to social demand: A landscape-scale analysis of the water regulation service. Landsc. Ecol. 2014, 29, 1069–1082. [Google Scholar] [CrossRef]
  18. Stürck, J.; Poortinga, A.; Verburg, P.H. Mapping ecosystem services: The supply and demand of flood regulation services in Europe. Ecol. Indic. 2014, 38, 198–211. [Google Scholar] [CrossRef]
  19. Baró, F.; Haase, D.; Gómez-Baggethun, E.; Frantzeskaki, N. Mismatches between ecosystem services supply and demand in urban areas: A quantitative assessment in five European cities. Ecol. Indic. 2015, 55, 146–158. [Google Scholar] [CrossRef] [Green Version]
  20. Castro, A.J.; Verburg, P.H.; Martín-López, B.; Garcia-Llorente, M.; Cabello, J.; Vaughn, C.C.; López, E. Ecosystem service trade-offs from supply to social demand: A landscape-scale spatial analysis. Landsc. Urban Plan. 2014, 132, 102–110. [Google Scholar] [CrossRef]
  21. Blanco, V.; Holzhauer, S.; Brown, C.; Lagergren, F.; Vulturius, J.; Lindeskog, M.; Rounsevell, M.D.A. The effect of forest owner decision-making, climatic change and societal demands on land-use change and ecosystem service provision in Sweden. Ecosyst. Serv. 2017, 23, 174–208. [Google Scholar] [CrossRef] [Green Version]
  22. Schulp, C.J.E.; Lautenbach, S.; Verburg, P.H. Quantifying and mapping ecosystem services: Demand and supply of pollination in the European Union. Ecol. Indic. 2014, 36, 131–141. [Google Scholar] [CrossRef]
  23. Larondelle, N.; Lauf, S. Balancing demand and supply of multiple urban ecosystem services on different spatial scales. Ecosyst. Serv. 2016, 22, 18–31. [Google Scholar] [CrossRef]
  24. Turpie, J.K.; Forsythe, K.J.; Knowles, A.; Blignaut, J.; Letley, G. Mapping and valuation of South Africa’s ecosystem services: A local perspective. Ecosyst. Serv. 2017, 27, 179–192. [Google Scholar] [CrossRef]
  25. Hitzhusen, F.J. Economic Valuation of River Systems; Edward Elgar: Cheltenham, UK, 2007. [Google Scholar]
  26. Peng, J.; Yang, Y.; Xie, P.; Liu, Y.X. Zoning for the construction of green space ecological networks in Guangdong Province based on the supply and demand of ecosystem services. Acta Ecol. Sin. 2017, 37, 4562–4572. [Google Scholar]
  27. Nedkov, S.; Burkhard, B. Flood regulating ecosystem services—Mapping supply and demand, in the Etropole municipality, Bulgaria. Ecol. Indic. 2012, 21, 67–79. [Google Scholar] [CrossRef]
  28. Peña, L.; Casado-Arzuaga, I.; Onaindia, M. Mapping recreation supply and demand using an ecological and a social evaluation approach. Ecosyst. Serv. 2015, 13, 108–118. [Google Scholar] [CrossRef]
  29. Wang, J.; Zhai, T.; Lin, Y.F.; Kong, X.; He, T. Spatial imbalance and changes in supply and demand of ecosystem services in China. Sci. Total Environ. 2019, 657, 781–791. [Google Scholar] [CrossRef]
  30. Wei, H.; Fan, W.; Wang, X.; Lu, N.; Dong, X.; Zhao, Y.; Zhao, Y. Integrating supply and social demand in ecosystem services assessment: A review. Ecosyst. Serv. 2017, 25, 15–27. [Google Scholar] [CrossRef]
  31. Xing, L.; Xue, M.; Hu, M. Dynamic simulation and assessment of the coupling coordination degree of the economy–resource–environment system: Case of Wuhan City in China. Environ. Manag. 2019, 230, 474–487. [Google Scholar] [CrossRef]
  32. Li, W.; Wang, Y.; Xie, S.; Cheng, X. Coupling coordination analysis and spatiotemporal heterogeneity between urbanization and ecosystem health in Chongqing municipality, China. Sci. Total Environ. 2021, 791, 148311. [Google Scholar] [CrossRef]
  33. Liu, W.; Zhan, J.; Zhao, F.; Wei, X.; Zhang, F. Exploring the coupling relationship between urbanization and energy eco-efficiency: A case study of 281 prefecture-level cities in China. Sustain. Cities Soc. 2020, 64, 102563. [Google Scholar] [CrossRef]
  34. Guan, Q.; Hao, J.; Ren, G.; Li, M.u.; Chen, A.; Duan, W.; Chen, H. Ecological indexes for the analysis of the spatial–temporal characteristics of ecosystem service supply and demand: A case study of the major grain producing regions in Quzhou, China. Ecol. Indic. 2020, 108, 105748. [Google Scholar] [CrossRef]
  35. Li, P.; Liu, C.; Liu, L.; Wang, W. Dynamic Analysis of Supply and Demand Coupling of Ecosystem Services in Loess Hilly Region: A Case Study of Lanzhou, China. Chin. Geogr. Sci. 2021, 31, 276–296. [Google Scholar] [CrossRef]
  36. Yang, Z.; Zhan, J.; Wang, C.; Twumasi-Ankrah, M.J. Coupling coordination analysis and spatiotemporal heterogeneity between sustainable development and ecosystem services in Shanxi Province, China. Sci. Total Environ. 2022, 836, 155625. [Google Scholar] [CrossRef]
  37. Huang, M.D.; Xiao, Y.; Xu, J.; Liu, J.Y.; Wang, Y.Y.; Gan, S.; Lv, S.X.; Xie, G.D. A Review on the Supply-Demand Relationship and Spatial Flows of Ecosystem Services. J. Resour. Ecol. 2022, 13, 925–935. [Google Scholar]
  38. Sun, X.; Tang, H.; Yang, P.; Hu, G.; Liu, Z.; Wu, J. Spatiotemporal patterns and drivers of ecosystem service supply and demand across the conterminous United States: A multiscale analysis. Sci. Total Environ. 2020, 703, 135005. [Google Scholar] [CrossRef]
  39. Peng, J.; Tian, L.; Liu, Y.X.; Zhao, M.Y.; Hu, Y.N.; Wu, J.S. Ecosystem services response to urbanization in metropolitan areas: Thresholds identification. Sci. Total Environ. 2017, 607-608, 706–714. [Google Scholar] [CrossRef]
  40. Wang, L.; Gong, J.; Ma, S.; Wu, S.; Zhang, X.; Jiang, J. Ecosystem service supply–demand and socioecological drivers at different spatial scales in Zhejiang Province, China. Ecol. Indic. 2022, 140, 109058. [Google Scholar] [CrossRef]
  41. Goldenberg, R.; Kalantari, Z.; Cvetkovic, V.; Mörtberg, U.; Deal, B.; Destouni, G. Distinction, quantification and mapping of potential and realized supplydemand of flow-dependent ecosystem services. Sci. Total Environ. 2017, 593–594, 599–609. [Google Scholar] [CrossRef] [PubMed]
  42. Wilkerson, M.L.; Mitchell, M.G.E.; Shanahan, D.; Wilson, K.A.; Ives, C.D.; Lovelock, C.E.; Rhodes, J.R. The role of socio-economic factors in planning and managing urban ecosystem services. Ecosyst. Serv. 2018, 31, 102–110. [Google Scholar] [CrossRef]
  43. Wang, S.F.; Zhuang, Y.N.; Cao, Y.G.; Yang, K. Ecosystem Service Assessment and Sensitivity Analysis of a Typical Mine–Agriculture–Urban Compound Area in North Shanxi, China. Land 2022, 11, 1378. [Google Scholar] [CrossRef]
  44. Xu, M.J.; Feng, Q.; Zhang, S.R.; Lv, M.; Duan, B.L. Ecosystem Services Supply–Demand Matching and Its Driving Factors: A Case Study of the Shanxi Section of the Yellow River Basin, China. Sustainability. 2023, 15, 11016. [Google Scholar] [CrossRef]
  45. Pan, H.H.; Wang, J.Q.; Du, Z.Q.; Wu, Z.T.; Zhang, H.; Ma, K.M. Spatiotemporal evolution of ecosystem services and its potential drivers in coalfields of Shanxi Province, China. Ecol. Indic. 2023, 148, 110109. [Google Scholar] [CrossRef]
  46. Hu, B.A.; Kang, F.F.; Han, H.R.; Cheng, X.Q.; Li, Z.Z. Exploring drivers of ecosystem services variation from a geospatial perspective: Insights from China’s Shanxi Province. Ecol. Indic. 2021, 131, 108188. [Google Scholar] [CrossRef]
  47. Wang, J.F.; Li, Y.; Wang, S.; Li, Q.; Liu, X.L. Assessment of Multiple Ecosystem Services and Ecological Security Pattern in Shanxi Province, China. Int. J. Environ. Res. Public Health 2023, 20, 4819. [Google Scholar] [CrossRef]
  48. National Bureau of Statistics of China. China Statistical Yearbook (2022); China Statistics Press: Beijing, China, 2023. [Google Scholar]
  49. Guo, S.; Ma, Y. Comprehensive evaluation for sustainable development capacity of resource-based region. Chin. Popul. Resour. Environ. 2017, 27, 72–79. [Google Scholar]
  50. Yang, G.; Ge, Y.; Xue, H.; Yang, W.; Shi, Y.; Peng, C.; Du, Y.; Fan, X.; Ren, Y.; Chang, J. Using ecosystem service bundles to detect trade-offs and synergies across urban–rural complexes. Landsc. Urban Plan. 2015, 136, 110–121. [Google Scholar] [CrossRef]
  51. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses—A Guide to Conservation Planning; Agricultural Handbook No. 537; US Department of Agriculture Science and Education Administation: Washington, DC, USA, 1978. [Google Scholar]
  52. Peng, J.; Chen, X.; Liu, Y.X.; Lu, H.L.; Hu, X.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]
  53. Lyu, R.; Clarke, K.; Zhang, J.M.; Feng, J.L.; Jia, X.H.; Li, J.J. Spatial correlations among ecosystem services and their socio-ecological driving factors_ A case study in the city belt along the Yellow River in Ningxia, China. Appl. Geogr. 2019, 108, 64–73. [Google Scholar] [CrossRef] [Green Version]
  54. Potter, C.S.; Randerson, J.T.; Field, C.B.; Matson, P.A.; 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]
  55. Crabtree, R.; Potter, C.; Mullen, R.; Sheldon, J.; Huang, S.; Harmsen, J.; Rodman, A.; Jean, C. A modeling and spatio-temporal analysis framework for monitoring environmental change using NPP as an ecosystem indicator. Remote Sens. Environ. 2009, 113, 1486–1496. [Google Scholar] [CrossRef] [Green Version]
  56. Raudsepp-Hearne, C.; Peterson, G.D.; Bennett, E.M. Ecosystem service bundles for analyzing tradeoffs in diverse landscapes. Proc. Natl. Acad. Sci. USA 2010, 107, 5242–5247. [Google Scholar] [CrossRef] [PubMed]
  57. Zhao, T.; Ouyang, Z.; Zheng, H.; Wang, X.; Hong, M. Forest ecosystem services and their valuation in China. Nat. Resour. 2004, 19, 480–491. [Google Scholar]
  58. Williams, J.R. The erosion-productivity impact calculator (EPIC) model: A case history. Philos. Trans. R. Soc. B 1990, 329, 421–428. [Google Scholar]
  59. Mccool, D.K.; Foster, G.R.; Mutchler, C.K.; Meyer, L.D. Revised Slope Length Factor for the Universal Soil Loss Equation. Trans. ASAE 1989, 32, 1571–1576. [Google Scholar] [CrossRef]
  60. Liu, B.; Nearing, M.; Shi, P.; Jia, Z. Slope Length Effects on Soil Loss for Steep Slopes. Soil. Sci. Soc. Am. J. 2000, 64, 1759–1763. [Google Scholar] [CrossRef] [Green Version]
  61. Cai, C.; Ding, S.; Shi, Z.; Huang, L.; Zhang, G. Study of applying USLE and geographical information system IDRISI to predict soil erosion in small watershed. Soil Water Conserv. 2000, 14, 19–24. [Google Scholar]
  62. Kumar, A.; Devi, M.; Deshmukh, B. Integrated remote sensing and geographic information system based RUSLE modelling for estimation of soil loss in Western Himalaya, India. Water Resour. Manag. 2014, 28, 3307–3317. [Google Scholar] [CrossRef]
  63. Mouchet, M.A.; Paracchini, M.L.; Schulp, C.J.E.; Stürck, J.; Verkerk, P.J.; Verburg, P.H.; Lavorel, S. Bundles of ecosystem(dis)services and multifunctionality across European landscapes. Ecol. Indic. 2017, 73, 23–28. [Google Scholar] [CrossRef] [Green Version]
  64. Zhang, K.; Lü, Y.; Fu, B.; Yin, L.; Yu, D. The effects of vegetation coverage changes on ecosystem service and their threshold in the Loess Plateau. Acta Geogr. Sinica 2020, 75, 949–960. [Google Scholar]
  65. Wang, M.; Bai, Z.; Dong, X. Land Consolidation Zoning in Shaanxi Province based on the Supply and Demand of Ecosystem Services. China Land Sci. 2018, 32, 73–80. [Google Scholar] [CrossRef] [Green Version]
  66. Zhai, T.; Wang, J.; Jin, Z.; Qi, Y.; Fang, Y.; Liu, J. Did improvements of ecosystem services supply-demand imbalance change environmental spatial injustices? Ecol. Indic. 2020, 111, 106068. [Google Scholar] [CrossRef]
  67. Ouyang, X.; Wang, Z.; Zhu, X. Construction of the ecological security pattern of urban agglomeration under the framework of supply and demand of ecosystem services using Bayesian network machine learning: Case study of the Changsha–Zhuzhou–Xiangtan urban agglomeration. China. Sustainability 2019, 11, 4616. [Google Scholar] [CrossRef] [Green Version]
  68. Zhou, S.; Yang, L.; Gao, R.; Wang, X.; Gao, X.; Nie, W.; Xu, P.; Zhang, Q.; Wang, W. A comparison study of carbonaceous aerosols in a typical North China Plain urban atmosphere: Seasonal variability, sources and implications to haze formation. Atmos. Envion. 2017, 149, 95–103. [Google Scholar] [CrossRef]
  69. Shi, T.; Yang, S.; Zhang, W.; Zhou, Q. Coupling coordination degree measurement and spatiotemporal heterogeneity between economic development and ecological environment—Empirical evidence from tropical and subtropical regions of China. J. Clean. Prod. 2020, 244, 118739. [Google Scholar] [CrossRef]
  70. Fu, B.; Wang, S.; Liu, Y.u.; Liu, J.; Liang, W.; Miao, C. Hydrogeomorphic Ecosystem Responses to Natural and Anthropogenic Changes in the Loess Plateau of China. Annu. Rev. Earth Planet. Sci. 2017, 45, 223–243. [Google Scholar] [CrossRef]
  71. Hu, M.; Li, Z.; Wang, Y.; Jiao, M.; Li, M.; Xia, B. Spatio-temporal changes in ecosystem service value in response to land-use/cover changes in the Pearl River Delta. Resour. Conserv. Recycl. 2019, 149, 106–114. [Google Scholar] [CrossRef]
  72. Chen, W.; Chi, G.; Li, J. The spatial aspect of ecosystem services balance and its determinants. Land Use Policy 2020, 90, 104263. [Google Scholar] [CrossRef]
  73. Wang, J.F.; Li, X.H.; Christakos, G.; Liao, Y.L.; Zhang, T.; Gu, X.; Zheng, X.Y. Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun region, China. Int. J. Geogr. Inf. Sci. 2010, 24, 107–127. [Google Scholar] [CrossRef]
  74. Wang, J.; Zhang, T.; Fu, B. A measure of spatial stratified heterogeneity. Ecol. Indic. 2016, 67, 250–256. [Google Scholar] [CrossRef]
  75. Zhang, Y.; Lu, X.; Liu, B.; Wu, D.; Fu, G.; Zhao, Y.; Sun, P. Spatial relationships between ecosystem services and socioecological drivers across a large-scale region: A case study in the Yellow River Basin. Sci. Total Environ. 2021, 766, 142480. [Google Scholar] [CrossRef]
  76. Liu, H.; Tang, F.; Ding, Y.; Zhang, Y.; Guo, X.; Tan, J.; Cheng, Y. Temporal and spatial evolution characteristics of the coupling between county high-quality development and ecosystem services in Shanxi Province. Arid. Zone Res. 2022, 39, 1234–1245. [Google Scholar]
  77. Ma, Y.; Qian, J.; Su, Z. Climate Change and Its Impact on Land Desertificatiin in Northwestern Shanxi Province. Desert Res. 2011, 31, 1585–1589. [Google Scholar]
  78. Seppelt, R.; Dormann, C.F.; Eppink, F.V.; Lautenbach, S.; Schmid, S. A quantitative review of ecosystem service studies: Approaches, shortcomings and the road ahead. Appl. Ecol. 2011, 48, 630–636. [Google Scholar] [CrossRef]
  79. Yang, M.; Zhao, X.; Wu, P.; Hu, P.; Gao, X. Quantification and spatially explicit driving forces of the incoordination between ecosystem service supply and social demand at a regional scale. Ecol. Indic. 2022, 137, 108764. [Google Scholar] [CrossRef]
  80. Xin, R.; Skov-Petersen, H.; Zeng, J.; Zhou, J.; Li, K.; Hu, J.; Liu, X.; Kong, J.; Wang, Q. Identifying key areas of imbalanced supply and demand of ecosystem services at the urban agglomeration scale: A case study of the Fujian Delta in China. Sci. Total Environ. 2021, 791, 148173. [Google Scholar] [CrossRef]
  81. Hou, W.; Zhou, W.; Li, J.; Li, C. Simulation of the potential impact of urban expansion on regional ecological corridors: A case study of Taiyuan, China. Sustain. Cities Soc. 2022, 83, 103933. [Google Scholar] [CrossRef]
  82. Bryan, B.A.; Gao, L.; Ye, Y.; Sun, X.; Connor, J.D.; Crossman, N.D.; Stafford-Smith, M.; Wu, J.; He, C.; Yu, D.; et al. China’s response to a national land-system sustainability emergency. Nature 2018, 559, 193–204. [Google Scholar] [CrossRef]
  83. Ouyang, Z.; Zheng, H.; Xiao, Y.; Polasky, S.; Liu, J.; Xu, W.; Wang, Q.; Zhang, L.; Xiao, Y.; Rao, E.; et al. Improvements in ecosystem services from investments in natural capital. Science 2016, 352, 1455–1459. [Google Scholar] [CrossRef]
  84. Schneider, A.; Mertes, C.M. Expansion and growth in Chinese cities, 1978–2010. Environ. Res. Lett. 2014, 9, 024008. [Google Scholar] [CrossRef]
  85. Kim, J.H.; Jobbagy, E.G.; Jackson, R.B. Trade-offs in water and carbon ecosystem services with land-use changes in grasslands. Ecol. Appl. 2016, 26, 1633–1644. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Shanxi Forestry and Grassland Bureau. Available online: http://lcj.shanxi.gov.cn/lcfm/lcgk/ (accessed on 16 June 2023).
  87. Shen, J.; Li, S.; Wang, H.; Wu, S.; Liang, Z.; Zhang, Y.; Wei, F.; Li, S.; Ma, L.; Wang, Y.; et al. Understanding the spatial relationships and drivers of ecosystem service supply-demand mismatches towards spatially-targeted management of social-ecological system. J. Clean. Prod. 2023, 406, 136882. [Google Scholar] [CrossRef]
  88. Schirpke, U.; Kohler, M.; Leitinger, G.; Fontana, V.; Tasser, E.; Tappeiner, U. Futureimpacts of changing land-use and climate on ecosystem services of mountain grassland and their resilience. Ecosyst. Serv. 2017, 26, 79–94. [Google Scholar] [CrossRef] [PubMed]
  89. Liu, Y.; Bi, J.; Lv, J.S.; Ma, Z.W.; Wang, C. Spatial multi-scale relationships of ecosystem services: A case study using a geostatistical methodology. Sci. Rep. 2017, 7, 9486. [Google Scholar] [CrossRef] [Green Version]
  90. Liu, Y.; Li, T.; Zhao, W.W.; Wang, S.; Fu, B. Landscape functional zoning at a county level based on ecosystem services bundle: Methods comparison and management indication. J. Environ. Manag. 2019, 249, 109315. [Google Scholar] [CrossRef]
  91. Raudsepp-Hearne, C.; Peterson, G.D. Scale and ecosystem services: How do observation, management, and analysis shift with scale-lessons from Québec. Ecol. Soc. 2016, 21, 16. [Google Scholar] [CrossRef] [Green Version]
  92. Cui, F.; Tang, H.; Zhang, Q.; Wang, B.; Dai, L. Integrating ecosystem services supply and demand into optimized management at different scales: A case study in Hulunbuir, China. Ecosyst. Serv. 2019, 39, 1200984. [Google Scholar] [CrossRef]
  93. Ding, T.; Chen, J.; Fang, L.; Ji, J.; Fang, Z. Urban ecosystem services supply-demand assessment from the perspective of the water-energy-food nexus. Sustain. Cities Soc. 2023, 90, 104401. [Google Scholar] [CrossRef]
Figure 1. Study area location in China (a), elevation (b), and land use/land cover type of Shanxi Province (c).
Figure 1. Study area location in China (a), elevation (b), and land use/land cover type of Shanxi Province (c).
Applsci 13 09262 g001
Figure 2. Flowchart of the proposed methodology.
Figure 2. Flowchart of the proposed methodology.
Applsci 13 09262 g002
Figure 3. Spatial distributions of the key five ESs and ESSI across the 107 counties of Shanxi Province in 2000 and 2020.
Figure 3. Spatial distributions of the key five ESs and ESSI across the 107 counties of Shanxi Province in 2000 and 2020.
Applsci 13 09262 g003
Figure 4. Spatial distributions of ES demand and the ESDI across the 107 counties of Shanxi Province in 2000 and 2020.
Figure 4. Spatial distributions of ES demand and the ESDI across the 107 counties of Shanxi Province in 2000 and 2020.
Applsci 13 09262 g004
Figure 5. The coupling coordination degree between ES supply and demand of Shanxi Province in 2000 and 2020.
Figure 5. The coupling coordination degree between ES supply and demand of Shanxi Province in 2000 and 2020.
Applsci 13 09262 g005
Figure 6. Factor-detected results of socio-ecological variables of CCD using GDM in 2000 (a) and 2020 (b). “***” p < 0.001, “**” p < 0.01, “*” p < 0.05. Panels (a) and (b) display the prominent influencing factors of CCD in 2000 and 2020, respectively. All factors are organized in descending order based on their influencing coefficients (q-statistic values).
Figure 6. Factor-detected results of socio-ecological variables of CCD using GDM in 2000 (a) and 2020 (b). “***” p < 0.001, “**” p < 0.01, “*” p < 0.05. Panels (a) and (b) display the prominent influencing factors of CCD in 2000 and 2020, respectively. All factors are organized in descending order based on their influencing coefficients (q-statistic values).
Applsci 13 09262 g006
Figure 7. Interaction-detected results of socio-ecological variables of CCD using GDM for 2000 (a) and 2020 (b). “*” indicates nonlinear enhancement: q (X1∩X2)> q (D1) + q (D2); “+” indicates bivariate enhancement: q (D1 ∩ D2) > q (D1) and q (D2). A deeper shade of blue indicates a smaller interaction coefficient between X1 and X2 on CCD, and a darker shade of red signifies a larger interaction coefficient.
Figure 7. Interaction-detected results of socio-ecological variables of CCD using GDM for 2000 (a) and 2020 (b). “*” indicates nonlinear enhancement: q (X1∩X2)> q (D1) + q (D2); “+” indicates bivariate enhancement: q (D1 ∩ D2) > q (D1) and q (D2). A deeper shade of blue indicates a smaller interaction coefficient between X1 and X2 on CCD, and a darker shade of red signifies a larger interaction coefficient.
Applsci 13 09262 g007
Table 1. Datasets used in the study.
Table 1. Datasets used in the study.
ES VariableData TypeSpatial ResolutionData Source
Land use/land coverRaster30 mNational Geomatics Center of China (http://www.globallandcover.com/GLC30Download/index.aspx, accessed on 11 December 2020)
NDVIRaster250 mNational Aeronautics and Space Administration and United States Geological Survey
(http://e4ftl01.cr.usgs.gov/MOLT/MOD13Q1.006/, accessed on 7 December 2020)
DEMRaster90 mGeospatial Data Cloud (https://www.gscloud.cn/#page 1, accessed on 20 December 2020)
Meteorological dataNumericSitesChina Meteorological Data Sharing Service System
(http://www.escience.gov.cn/metdata/page/index.html, accessed on 5 November 2020)
Soil databaseRaster30 arc-secondHarmonized World Soil Database
(http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonied-world-soil-datebse-v12/en/, accessed on 1 December 2020)
Administrative mapVectorCountyNational Geomatics Center of China (http://ngcc.sbsm.gov.cn/ngcc/, accessed on 11 December 2020)
Crop yield, GDP and populationNumericCountyStatistical Yearbook
Table 2. ES indicators from the MA categories and their quantitative methods, units, and ES variable requirements.
Table 2. ES indicators from the MA categories and their quantitative methods, units, and ES variable requirements.
ES IndicatorCodeModel or ProxyUnitAssessment Process
Crop productionCroCrop yield per square kilometerton/km2Crop production was calculated by dividing the crop yield of each county by its territory to illustrate per-unit provision service.
Water retentionWretWater balance equationm3/km2 T Q = i = 1 j ( P i R i E T i )
TQ is water conservation, Pi is the annual average rainfall (mm), Ri is the annual average surface runoff (mm), and ETi is the annual evapotranspiration (mm).
Soil conservationSconUSLE (Universal Soil Loss Equation)t/(hm2·a) Δ A = R × K × L × S ( 1 C × P )
ΔA = soil conservation (t/ (hm2·a)), R = rainfall erosivity index (MJ·mm/(hm2·h·a)), K = soil erodibility factor (t·hm2·h/ (MJ·mm·hm2)), L S = slope length and steepness factor (unitless), C = cover and management factor (unitless), P = conservation practice factor (unitless). The parameters R were from Wischmeier and Smith [57], K from Williams [58], L S from McCool et al. [59] and Liu et al. [60], C from Cai et al. [61], and P from Kumar et al. [62].
Carbon sequestrationCseqCASA (Carnegie–Ames–Stanford approach)kg C/km2 N P P = A P A R × ξ
NPP = net primary productivity (g C/m2), APAR = absorbed photosynthetic active radiation (MJ/m), ξ = the utilization rate of light energy (g C/MJ).
Outdoor recreationRecTourists per square kilometerpersons/km2Outdoor recreation was calculated via the area of forest land in each county.
Table 3. Details of the driving variables for coupling coordination between ESSI and ESDI in this study.
Table 3. Details of the driving variables for coupling coordination between ESSI and ESDI in this study.
VariableCodeDescriptionUnit
ElevationDEMDerived from the SRTM3 global digital elevation modelMeter
SlopeSLOPEDerived from the SRTM3 global digital elevation modelDegree
PrecipitationPREAnnual trends of precipitation for the period 1956–2017mm
TemperatureTEMAnnual trends of temperature for the period 1956–2017°C
Normalized Difference Vegetation IndexNDVIVegetation cover%
CroplandCROPCounty land area that is occupied by area that is classified as cropland%
ForestlandFORESTCounty land area that is occupied by area that is classified as forest%
GrasslandGRASSCounty land area that is occupied by area that is classified as grassland%
Construction landCONCounty land area that is occupied by area that is classified as construction land%
PopulationPOPAnnual total populationperson
Economic levelGDPGross domestic productyuan
Urbanization rateURBANUrban population proportion%
Distance to the nearest countyCOUNTYDistance to the nearest county centerkm
Table 4. The statistical values of standardized values of ES supply, ES demand, and coupling coordination degree in 2000 and 2020.
Table 4. The statistical values of standardized values of ES supply, ES demand, and coupling coordination degree in 2000 and 2020.
ESSIESDICCD
Year200020202000202020002020
Minimum value0.5280.512 0.649 0.638 0.002 0.059
Maximum value0.859 0.870 0.979 0.980 0.566 0.427
Mean value0.670 0.649 0.044 0.067 0.2120.196
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Liu, B.; Sui, R. Evaluation and Driving Determinants of the Coordination between Ecosystem Service Supply and Demand: A Case Study in Shanxi Province. Appl. Sci. 2023, 13, 9262. https://doi.org/10.3390/app13169262

AMA Style

Zhang Y, Liu B, Sui R. Evaluation and Driving Determinants of the Coordination between Ecosystem Service Supply and Demand: A Case Study in Shanxi Province. Applied Sciences. 2023; 13(16):9262. https://doi.org/10.3390/app13169262

Chicago/Turabian Style

Zhang, Yushuo, Boyu Liu, and Renjing Sui. 2023. "Evaluation and Driving Determinants of the Coordination between Ecosystem Service Supply and Demand: A Case Study in Shanxi Province" Applied Sciences 13, no. 16: 9262. https://doi.org/10.3390/app13169262

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