Next Article in Journal
Treating the Symptoms as Well as the Root Causes: How the Digital Economy Can Mitigate the Negative Impacts of Land Resource Mismatches on Urban Ecological Resilience
Previous Article in Journal
Defining and Verifying New Local Climate Zones with Three-Dimensional Built Environments and Urban Metabolism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies in the Yangtze River Delta

1
College of Geography and Environment, Shandong Normal University, Jinan 250358, China
2
National Institute of Natural Hazards, Ministry of Emergency Management of China, Beijing 100085, China
3
School of Resources and Environmental Engineering, Ludong University, Yantai 264025, China
4
Shandong Institute of Territorial and Spatial Planning, Jinan 250014, China
5
State Key Laboratory of Water Environment Simulation, School of Environment, Beijing Normal University, Beijing 100875, China
*
Authors to whom correspondence should be addressed.
Land 2024, 13(9), 1462; https://doi.org/10.3390/land13091462
Submission received: 3 August 2024 / Revised: 30 August 2024 / Accepted: 6 September 2024 / Published: 9 September 2024

Abstract

:
The transformation of ecosystem structure leads to changes in ecosystem services (ESs) and their relationship. However, most research in this area has focused on particular scales and timeframes, often overlooking the significance of spatial and temporal variations. Therefore, we used the equivalent value method to evaluate seven typical ESs in the Yangtze River Delta (YRD) between 2000 and 2020: food production (FP), water supply (WS), climate regulation (CR), environmental purification (EP), soil conservation (SC), biodiversity maintenance (BM), and aesthetic landscape (AL). We further employed the Spearman correlation coefficient and bivariate Moran’s I to evaluate the relationship of ESs and their spatial heterogeneity at grid, township, county and city scales. Our results show that (1) All ESs except AL exhibited a fluctuating upward trend from 2000 to 2020, resulting in a total increase in ecosystem service (ES) value of RMB 650.63 billion. (2) Approximately 70% of the ES pairs demonstrated a synergistic relationship, with the exception of FP and other ESs, which mainly showed a trade-off. (3) With the scale increased from grid to city level, the degree of trade-off between FP and other ESs strengthened at different levels, while the synergy degree of among other ESs gradually decreased. (4) The relationship between ESs demonstrated strong spatial heterogeneity, with FP and other ESs exhibiting trade-offs primarily in the northern and southern YRD, while other ES pairs exhibited mostly synergy in these regions. This study provides scientific information for governments to optimize land use distribution and improve ESs.

1. Introduction

Nature provides fundamental resources for human survival. Ecosystem services (ESs) encompass various benefits that humans derive directly or indirectly from nature [1]. As the concept of sustainable development and ES gains more attention, people are increasingly recognizing that the ultimate research goal on ESs is incorporating this concept into urban planning and management [2,3], thereby promoting the coordinated development of natural ecological systems and socio-economic systems. However, due to human activities and global climate change, significant changes in ecosystem structure have negatively impacted the provisioning of ecosystem services (ESs) [4,5]. Relevant studies have demonstrated that 15 out of the 24 global ESs have been degraded to varying degrees [6]. Thus, different countries have taken different measures aimed at improving the level of ESs and human well-being based on the related studies.
The evaluation of ES values (ESVs) primarily relies on the service function price per unit area and value equivalent factor per unit area [7]. The former one is primarily concerned with value accounting that is predicated on the evaluation of the physical quantity of ESs, which needs many ecological parameters. The latter one obtains the equivalent value based on meta-analysis and the area of each ecosystem, which is more effective and widely used in research. On the basis of evaluation of ES material or value, scholars have conducted extensive research across a multitude of domains, such as driving mechanisms, ES relationships, human welfare, landscape planning, and ecosystem management [8]. A large number of studies have shown that the impacts of drivers vary across different ESs [9]. The fluctuating dynamics of ESs leads to evolving trade-offs and synergies among them, which are crucial for sustainable ecosystem management [10]. Based on their different changing trends, the relationship can be classified into synergy, trade-off, and neutrality, which means that one ES may have positive, negative, and negligible effects on another ES [11]. Trade-off/synergy analysis methods mainly include correlation analysis [12], scenario analysis [13], and ES bundles [14]. Compared with other methods, correlation analysis has the advantages of being simple and intuitive, requiring a small amount of data and having high flexibility. Based on the previous studies, this study will use correlation analysis to analyze the ES trade-offs/synergies based on the ES value which is calculated by value equivalent factor per unit area.
The relationship between ESs is not static but evolves with changes in time, spatial scale, and surrounding environment [15]. Different studies have analyzed trade-offs/synergies at global [16], national [17], regional, and landscape scales [18], where ESs exhibit complex interactions [19]. Drivers, including land use, socio-economic conditions, and meteorological factors, evolve over time and vary across spatial scales [20], potentially leading to diverse change trends among various ESs and affecting the relationship between ESs. However, most studies have researched the relationship between ESs at a single scale, neglecting their scale effects. The relationship between ESs also changes over time due to changes in key factors [21]. Studies in the Baota Mountain area by Hou [22] found that the trade-off and synergy degrees between grain and aquatic production varied with time at grid scale. However, most existing research has only examined their relationship at a single point in time. In addition, spatial heterogeneity of the relationship between ESs (spatial trade-offs) arises from different ecological process depending on the diverse natural environment and socio-economic situation. Therefore, the relationships between the same pairs of ESs may vary due to regional environmental differences. Many studies have used Moran’s Index to analyze the spatial distribution of individual ESs, which does not reflect the relationships between ESs as it pertains to only one ES. Furthermore, few studies have explored the spatial trade-offs in the relationships between ESs. Consequently, we have introduced the bivariate Moran’s Index, which considers two ESs, to assess the spatial trade-offs and synergies between them. A limited number of studies have explored the temporal changes and scale effects in the relationships between ESs, leading to an incomplete understanding of them [22,23,24,25,26]. Current studies have only revealed the average situation of ES trade-offs/synergies in a specific research area [27], without identifying specific locations. Understanding the temporal and spatial characteristics of the relationships between ESs is fundamental for urban development planning and is increasingly incorporated into decision-making by governments and environmental organizations [28]. Therefore, there has been a growing body of research on the scale effects of ES trade-offs/synergies, and some progress has been achieved; however, a considerable gap still exists in the depth and breadth of the related research. Firstly, the range of ESs addressed in current research is not comprehensive enough, and the temporal scope of the studies is also inadequate. Secondly, the prevailing research has a strong emphasis on the scale effects of the relationship between ESs, frequently neglecting the examination of their spatial distribution patterns.
Benefiting from unique terrain, climate, and other natural factors, the Yangtze River Delta (YRD) has formed a unique ecosystem which has a great influence on the ecological security of China. Additionally, the YRD is one of China’s fastest developing and most developed regions, which is very significant to China’s economic growth [29]. But, a large amount of ecological land has been replaced by gradually expanding cities, and the area of ecological land has decreased significantly. Additionally, after years of efficient and rapid development, the YRD is facing many ecological issues, such as inadequate water supply, air pollution, and water and soil loss. These issues have made the ecosystem structure more vulnerable and altered the trade-offs/synergies among ESs [30]. Therefore, we selected the YRD as our study area with the objectives to (1) examine the spatial–temporal change trends for seven key ESs (Food Production (FP), Water Regulation (WR), Climate Regulation (CR), Soil Conservation (SC), Environmental Purification (EP), Biodiversity Maintenance (BM), and Aesthetic Landscape (AL)) from the year 2000 to 2020; (2) elucidate the dynamics of the relationship between ESs over time and their spatial scale effects, ranging from grid to township, county, and city levels; and (3) ascertain the spatial heterogeneity in the relationships between ESs. This research is intended to inform decision-makers in the development of strategies for land management and the promotion of sustainable ecosystem development.

2. Materials and Methods

2.1. Study Area

The YRD is situated in the lower reaches of the Yangtze River and proximal to the East China Sea and encompasses 26 cities (Figure 1a) [31]. The YRD spans an area of 210 thousand km2, accounting for 2.3% of China’s total area. Its topography chiefly entails lowlands in the center and highlands around the periphery. Notably, the altitude is mostly below ten meters (Figure 1b). The YRD has a subtropical monsoon climate and concurrent rainy and warm periods. The average temperature remains above 0 °C, and subtropical evergreen broad-leaved forests dominate the vegetation. Furthermore, the YRD is the Chinese region with the densest river network.
The YRD had a population of 227 million and a high level of urbanization in 2020. During the study period, the Gross Domestic Product (GDP) of the YRD grew significantly from RMB 1.86 trillion to RMB 24.47 trillion (in comparison to 2000), with an average annual growth rate of 12.8%, which was higher than the national average. As such, it is an essential contributor to China’s overall economic development. However, the expanding cities have replaced a substantial amount of arable and forested land, resulting in significant damage to the ecosystem. Remarkably, the area of arable land (irrigated cropland and dry land), forest land (broad-leaved forest and shrubbery), grassland (grassland, open and dense), and wetlands in the YRD declined by 12,800 km2, 650 km2, 330 km2, and 573 km2, respectively, between 2000 and 2020 (Figure 1c,d)).

2.2. Data Sources and Processing

The primary datasets utilized in this study encompass both land use and socio-economic information. We obtained land use data with a spatial resolution of 1000 × 1000 m from the Resource and Environmental Science Data Platform of the Chinese Academy of Sciences and refined the land use classification system (Table 1) in alignment with the methodologies established by Xie et al. [32], taking into account the unique regional characteristics of the YRD. Due to the minimal ESs provided by construction land or barren or sparsely vegetated regions, we assigned a value of zero to these areas in this study. Xie et al. [32] introduced the initial equivalent value for ESs per unit area, as detailed in Table 2, which represented the average ES levels across China. However, the intensity of ESs can vary considerably depending on the ecological conditions of different regions. This variability necessitates the adjustment of the equivalent value per unit area when assessing the ESV for a particular area. The detailed exposition of the methods and processes for equivalent value correction is shown in Section 3.1.
In this study, socio-economic data are sourced from the Statistical Yearbooks of various cities and provinces in the YRD and China from 2000 to 2020. Key indicators mainly included the agricultural, forestry, animal husbandry, fishery production value, forest coverage, tourism income, and demographic data for both China and the YRD. Additionally, we considered the prices, yields, cultivation areas, and producer price indexes (PPI) of the three principal cereal crops: wheat, corn, and rice. The PPI for agricultural products is also derived from the National Bureau of Statistics.

2.3. Methods

2.3.1. Equivalent Value Correction

The correction coefficient was introduced in this study to carry out the first and second corrections to the equivalent values per unit area. First, we revised the equivalent value of ESs according to the economic characteristics in the YRD. The adjustment coefficients for FP were recalibrated for various land use types, based on the per-unit-area output values of agriculture, forestry, animal husbandry, and fisheries within the YRD and the national averages [33]. This recalibration is detailed in Table 3. The equivalent value for CR, EP, SC, and BM was modified based on forest coverage [34]. Since tourism can indirectly reflect the aesthetic value, the value equivalent of AL can be revised according to the total tourism revenue and the total number of people. We have corrected the equivalent value of AL based on the tourism revenue per person [33]. The functional adjustment index of equivalent value factors of ESs in the YRD between 2000 and 2020 is shown in Table 3.
The magnitude of various ESs is intrinsically linked to the ecosystem’s biomass; an increase in biomass generally corresponds to a higher capacity for ES provision [35]. In accordance with this principle, we conducted a secondary correction of the equivalent value factors, which was based on the biomass factors and the respective proportions of the provinces and municipality within the YRD, namely Jiangsu, Zhejiang, Anhui, and Shanghai. The biomass factor for the YRD was determined to be 1.52 and was utilized as the secondary correction coefficient. The ES equivalent value underwent both primary and secondary revisions in accordance with the method outlined in Formula 1.
b i j = n f × B i j
where b i j represents the corrected unit equivalent value, i = 1, 2… 6 represents different types of ESs, and j = 1, 2… 8 represents different land use types; B i j represents the equivalent value prior to correction; n f represents the first and second correction coefficients.

2.3.2. Evaluation of Ecosystem Services

The equivalent value for ESs per unit area was ascertained by calculating one-seventh of the income generated from grain production per unit area (Formula (2)).
C t = 1 7 i = 1 i S i × P i × Q i S
where C t represents ESV of unit value equivalent (RMB 10,000/km2); i represents wheat, corn and rice. S i and S refer to the sown area of the ith grain and total sown area of three grains, respectively (km2); P i and Q i represent the price (RMB 10,000/kg) and the yield per unit land use area of the ith grain (kg/km2), respectively.
To ensure the comparability of ESVs across various time periods, we adjusted for price fluctuations by employing the PPI. This allowed us to calculate the equivalent prices of crops for different years, normalized to the year 2000. Between 2000 and 2020, the equivalent value factor for ESs in the YRD was RMB 7.12, 10.90, 17.27, 22.68, and 21.30, respectively, as detailed in Table A1.

2.3.3. Ecosystem Services Trade-Off/Synergy Analysis

This study employed GIS technology and correlation analysis to assess the trade-off and synergy intensity of seven key ESs across different spatial scales, including 1000 × 1000 m grid, township, county, and city levels, for the years 2000, 2005, 2010, 2015, and 2020. The correlation analysis methods mainly include Pearson and Spearman correlations. The selection between these methods was contingent upon the data’s adherence to a normal distribution. The Pearson correlation was applied in cases where data exhibited normality, while the Spearman correlation was utilized for non-normally distributed data [36]. The Kolmogorov–Smirnov test was utilized to assess the normality of the data distribution, which indicated that 64% of the ESs exhibited non-normal distributions. Therefore, Spearman correlation analysis was selected to examine the relationships among the ESs in this study. The calculation of the Spearman correlation coefficient is presented in Formula (3).
ρ = 1 6 i = 1 n d i 2 n ( n 2 1 )
where ρ represents the correlation coefficient; d i denotes the rank difference in variable values; i is the rank index of the variable value; n is the total number of variable values.

2.3.4. Spatial Heterogeneity of Ecosystem Service Trade-Offs/Synergies

Spatial autocorrelation refers to the possible interconnections between variables within a specific area, which is typically assessed using Moran’s I index. This statistical measure has been widely employed to investigate the spatial relationships of ESs [37]. The bivariate global Moran’s I (BGM) is a tool that evaluates the overall spatial correlation between two variables, reflecting the extent to which their spatial distributions are interrelated. On the other hand, the bivariate local Moran’s I (BLM) delves into the local spatial correlation between the two variables, highlighting the degree of clustering within proximate regions. The BLM identifies four distinct clustering patterns: High–High (H-H), Low–Low (L-L), High–Low (H-L), and Low–High (L-H), each representing different spatial associations. In this study, we employed the queen adjacency approach to calculate the BGM, which allowed us to evaluate the overall spatial correlation between ESs. Additionally, we utilized the BLM to quantify the local spatial agglomeration among these services. The specific methodologies for calculating BGM and BLM are delineated in Formulas 4 and 5, respectively. This analytical framework provides a comprehensive understanding of the spatial dynamics within the ESs under investigation.
I k l = i = 1 n j = 1 n W i j ( x k i x ¯ k ) ( x l i x ¯ l ) S 2 i = 1 n j = 1 n W i j
I k l = x k i x ¯ k S k 2 j = 1 n W i j x l j x ¯ l S l 2
where I k l and I k l represent the BGM and BLM, respectively; x k i ( x l j ) is the observed value of k (l) variable on space element I (j); x l i is the observed value of l variable on space element i; x ¯ k , x ¯ l are the average of the observed values of k and l variables; W i j is the spatial weight matrix; S k 2 , S l 2 are the variance of the observed values of k and l variables, respectively; n is the total number of geographical units in the study area.

3. Results

3.1. Spatial and Temporal Evolution of Ecosystem Services

3.1.1. Temporal Evolution of Ecosystem Services

The ESV per unit area generally displayed a trend of fluctuating ascent, characterized by an initial increase followed by a decline in the YRD (Table A1). By synthesizing the various land use areas within the YRD across different time periods, we have generated the ESV fluctuations for the seven distinct ESs over five distinct periods, as illustrated in Figure 2. Upon a comparative analysis of Figure 2 and Table 1, it was evident that the ESVs exhibited a trend of fluctuating increase. Generally, the ESVs witnessed a substantial increase from 2000 to 2020, amounting to a total rise of RMB 650.625 billion. The zenith of this growth was reached in 2015, with the total ESV peaking at RMB 1017.05 billion. This was followed by the year 2020, with the total ESV stood at RMB 982.31 billion, a slight decrease from the peak but still significantly higher than the initial value in 2000, which was recorded at RMB 331.68 billion. The most remarkable surge in the ESV occurred between the years 2005 and 2010, during which the value increased by RMB 304.469 trillion, marking a growth rate of 59%. When examining the change rates of the seven different ESVs from 2000 to 2020, they were ranked from highest to lowest as follows: CR, WR, EP, BM, SC, FP, and AL. Notably, the value of CR consistently emerged as the most significant contributor to the ESV, with its value peaking annually and constituting nearly one-third of the total ESV.

3.1.2. Spatial Evolution of Ecosystem Services

The ESV density across various land use types showed a distinct geographical pattern within the study area. It was observed that the ESV in the northern regions was relatively lower when compared to the southern regions. This disparity could be attributed to differences in land use practices, environmental conditions, and the inherent capacity of the ecosystems in these areas to provide services. Furthermore, the ESV associated with aquatic ecosystems, including water bodies, rivers, lakes, wetlands, and estuaries of the Yangtze River and Qiantang River, was notably higher than in their adjacent areas. This elevated ESV could be attributed to the significant ecological functions performed by these water bodies, such as WR, EP, BM, and AL. The higher ESV in these water-related land use types underscores the importance of preserving and enhancing the ecological integrity of these areas, as they contribute substantially to the overall environmental health and economic value of the region. Furthermore, the northeastern inland region of the YRD demonstrated a markedly lower ESV in comparison to the coastal zones of the YRD.
FP exhibited a spatial distribution that was relatively low in the southern part of the YRD but notably high in the northern areas (Figure 3). This region was characterized by its flat topography, abundance of rivers and lakes, and extensive tracts of arable land, which contributed to its status as the principal area for grain production within the YRD. The favorable conditions for agriculture in the northern plains have made it a vital contributor to the region’s food security. WR in the northern region of the YRD was comparatively lower than in the south. This difference could be attributed to the predominance of shrubbery and broad-leaved forests in the southern region, which were instrumental in water conservation due to their capacity to retain moisture and regulate water flow. Conversely, the northern YRD was characterized by extensive arable land, which necessitates significant water usage for irrigation and other agricultural practices. The high demand for water in agriculture can lead to a relative depletion of available water resources. Furthermore, the northeast coastal areas of the YRD also showed high values of WR. The distribution of CR mirrored that of SC, with lower values observed in the northern YRD and higher values in the south. This pattern primarily stemmed from the predominance of shrubbery and broad-leaved forests in the south, which were highly effective in regulating various environmental factors such as soil quality, temperature, humidity, wind speed, and precipitation. These forests played a crucial role in enhancing physical and mental comfort and conserving soil, thereby contributing significantly to the high values of CR and SC in the southern areas. The distribution pattern of EP, BM, and AL exhibited similarities across the region. Broadly speaking, the southern forested areas demonstrated higher values for these services compared to the northern ones, which were predominantly characterized by arable land.

3.2. Temporal Evolution of Ecosystem Services Trade-Offs/Synergies

The visual representation of these relationships is conveyed through a series of “bubbles”, where the size of each bubble corresponds to the magnitude of the correlation coefficient, indicating the strength of the relationship between the ESs (Figure 4). Additionally, the color of the bubbles signifies the direction of the correlation: the darker the color, the more pronounced the correlation coefficient. Red and blue bubbles indicate a positive and negative correlation, respectively. In the YRD, our findings revealed that a substantial majority, 70%, of the ES pairs exhibited synergistic interactions.

3.2.1. Spatial Scale Effects

The relationships between ESs observed at the grid, township, county, and city scales exhibited a notable degree of similarity. Although the correlation between ESs at grid scale was roughly similar to that of the three administrative regional scales, the correlation value and change trend are still quite different due to different research scales. Compared with the other three scales, the trade-off between FP and other ESs at the grid scale was weak, and most of them even showed weak synergy. The synergy between WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL was higher than that of the other three scales. CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL also showed medium- to high-intensity synergy.
At the township and county scales, the relationships between WR and the other ESs (CR, EP, SC, BM, and AL) exhibited a consistent and moderate level of synergy. This stability suggested a balanced and harmonious interaction among these services within these administrative scales, contributing to a resilient and sustainable ecological framework. However, the degree of synergy observed was more pronounced at the township scale compared to the county and city scales. At the city scale, the interactions between FP and the other ESs (CR, SC, and BM) were characterized by strong trade-off relationships, while other ES pairs demonstrated synergy.
As the scale expanded from the township to the city level, the synergistic relationships of CR-EP, CR-BM, CR-Al, CR-SC, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL exhibited a gradual decline. Conversely, the trade-offs of FP-CR, FP-EP, FP-SC, FP-BM, and FP-AL were most pronounced at the city scale. Interestingly, the strongest trade-off between FP and WR was observed at the township scale. From the township to city scales, the relationship of FP-WR underwent a transition from a weak trade-off to a weak synergy. The synergy of WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL experienced a rapid decline as the scale expanded. This reduction was particularly notable, with the initially weak synergies of WR-CR and WR-SC eventually turning into weak trade-offs at larger scales. Nevertheless, it persisted, with a consistent degree of robustness despite fluctuations in its intensity. Despite the shifts observed in other relationships between ESs, the synergy of CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL remained resilient. These relationships demonstrated a consistent level of robustness, even though their intensity may have varied.

3.2.2. Temporal Evolution

Compared with the other three scales, from 2000 to 2020, the trade-off between FP and other ESs at the grid scale gradually weakened and most of them turned into synergy. Similar to other scales, the synergy of WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL also showed a decreasing trend from 2000 to 2020. The synergistic relationship of CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL showed an increasing trend
At the township, county, and city scales, the trade-off degree of FP-WR, FP-CR, FP-EP, FP-SC, FP-BM, and FP-AL experienced a decline, and FP-WR even turned into weak synergy from 2000 to 2020. Moreover, the synergistic relationships of WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL also declined. Although there was a strong and significant synergistic effect of CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL, they generally showed a weakening trend at the township scale and an increasing trend at the county and city scales, but neither the strengthening nor the weakening trend is obvious.
Over time, the trade-off degrees of FP-WR, FP-CR, FP-SC, FP-BM, and FP-AL gradually diminished, and even partly turned into a weak synergy. Similarly, the level of synergy and its significance for WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL declined. Furthermore, the relationships of WR-CR and WR-SC evolved from synergy into weak trade-offs. In contrast, the ES pairs of CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL showed a significant synergy.

3.3. Spatial Heterogeneity of Ecosystem Services Trade-Offs/Synergies

One purpose of our study was to show the spatial distribution characteristics of the interaction between ESs in the YRD so as to better formulate and adjust the regional development plan. At the grid scale, the correlation between FP and other ESs was relatively weak. Therefore, we did not consider this particular scale in this part. Among the other three scales, the township scale had the smallest area and could more finely show the spatial heterogeneity of ESs. Attached Figure A1 and Figure A2 are the spatial distribution of the interaction between ESs at the county and city scale in 2020, respectively. By comparing the research results at city, county, and township scales, we found that the spatial distribution characteristics of the interaction between the ESs they reflect were very similar. Compared with the township scale, the county and city scales were the simplified or even more simplified versions of their characteristics. Only the township scale can reflect the spatial distribution characteristics of the interaction between ESs in the most detailed way.
Finally, considering time validity, we no longer focused on the spatial distribution characteristics of the interaction between ESs before 2020 and selected the township scale in 2020 to conduct spatial heterogeneity analysis of ES relationships and obtained the Local Indicator of Spatial Association (LISA) diagram between ESs that passed the significance test at the level of 0.05. Statistical analysis confirmed that there was a significant correlation between different ESs (0.05 level), as shown in Figure 5. The results of BGM are consistent with the previous analysis introduced in Section 3.2.1 and Section 3.2.2.
Utilizing the spatial distribution patterns of the relationships between ESs, we categorized the ES pairs into three distinct groups to better understand their interactions and spatial dynamics. The first group included pairs that involve FP with other ESs: FP-WR, FP-CR, FP-EP, FP-SC, FP-BM, and FP-AL. The second group consisted of pairs that involve WR with other services, including WR-CR, WR-EP, WR-SC, WR-BM, and WR-AL. The final group encompassed ES pairs that did not include FP or WR: CR-EP, CR-SC, CR-BM, CR-AL, EP-SC, EP-BM, EP-AL, SC-BM, SC-AL, and BM-AL.
The first group was characterized by a distinctive spatial distribution pattern: H-H and L-L associations were sparse, indicating that areas of high or low ES provision tended not to cluster together. Conversely, L-H and H-L associations were widespread, suggesting ES pairs in the first group mostly exhibited trade-off relationships in the YRD (Figure 6). The areas exhibiting trade-offs (H-L and L-H) were predominantly situated in the northern and southern parts, while the central region remained largely insignificant. L-L areas were observed to be dispersed around the areas of H-L, while H-H areas tended to be distributed around the L-L areas. Notably, the FP-WR pair did not exhibit a significant relationship in the northern of the YRD.
In the second group, the relationship between WR and other ESs showed two different patterns in the north and south of the YRD. The northern part was dominated by the synergistic area, which was mainly characterized by the distribution characteristics of L-L with a small amount of H-H, indicating that the supply of water and other ESs in this area was complementary. The south was mainly dominated by trade-offs, with a small number of synergistic areas. The distribution of L-H areas surrounded H-H areas. The northern and southern parts of the YRD presented a spatial distribution pattern in which high (low) service areas and low (high) service areas were interspersed with each other. In contrast, the central part of the YRD showed an insignificant presence of these synergistic patterns.
In the third group, the relationships between ESs in the study area predominantly demonstrated synergistic patterns, with only minor variations observed among different ES pairs. Areas characterized by L-L and H-H were notably concentrated in the northern and southern parts of the YRD, respectively. The central part of the YRD, however, remained an area of insignificance in terms of significant ES relationships. Regions with H-L and L-H patterns were found widespread around the L-L and H-H areas. Notably, the SC-BM and SC-AL exhibited large areas of H-L in the northern YRD. In contrast, the EP- SC had extensive areas of H-L distributed around the L-L regions, which may suggest a more complex interplay between these services.

4. Discussion

4.1. Spatio-Temporal Evolution and Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies

Land use is a pivotal determinant in shaping the dynamics of trade-offs and synergies between ESs [38]. The existence of trade-offs is primarily due to the inherent incompatibility between certain ESs and particular types of land use. Grasping the ES relationship is crucial for enhancing total ES level and human well-being. In the YRD, the majority of ES pairs demonstrated synergistic relationships, similar to ESs in other regions [39]. The relationship between FP and other ESs exhibited a trade-off situation, which was predominantly attributed to the spatial competition among different land uses such as arable land, forests, grasslands, and water bodies. The need to allocate space for agricultural production, which is essential for FP, often leads to the conversion of ecological land, thereby affecting other ESs like WR, CR, EP, SC, and BM. The majority of ES pairs between WR and other ESs showed moderate synergies; however, not all exhibited synergistic relationships. The synergies were typically observed among regulation, support, and cultural services, which is consistent with the findings reported in the study by Li [40].
When integrating ESs into regional planning, it is imperative to take into account the interrelationships of ESs across spatial and temporal scales. This multi-dimensional perspective is essential for a comprehensive understanding of how ESs interact within and across different landscapes and over time, enabling more effective and sustainable land use strategies. Previous research has predominantly focused on examining the relationships between ESs at a particular scale, often neglecting the influence of spatial scale and temporal dynamics on these relationships [41]. This study revealed that the synergies of most ES pairs, except CR with other ESs and SC with other ESs, weakened with increasing scale. In 2020, a notable transition occurred in the dynamics of ES interactions as the spatial scale expanded from grid to city level. Specifically, the relationships of WR-CR and WR-SC shifted from weak synergies to weak trade-offs. Concurrently, the trade-offs between FP and other ESs, with the exception of the FP-WR, became more pronounced. Additionally, this study noted that the significance of the relationships between ESs tended to diminish as the spatial scale increased. This observation is in line with the findings reported in the research conducted by Pan [27].
The response of certain ESs to environmental changes has been identified as being relatively slow [37]. Consequently, it is imperative to take a long-term perspective when considering the evolution of the relationship between ESs. For our research, we selected the time frame from 2000 to 2020 for analysis. Within this two-decade span, we observed that the trade-offs of most ES pairs experienced a weakening trend, while synergies among them showed an increase. This suggests a more harmonious coexistence of these services over time.

4.2. Spatial Heterogeneity of Ecosystem Service Trade-Offs/Synergies for Policy Implications

In terms of the spatial heterogeneity of ES relationships, FP was found to be negatively correlated with other ESs. In contrast, WR demonstrated a weak positive correlation with other ESs, and the remaining ESs showed a moderate positive correlation. It was particularly noteworthy that the value of BGM for the relationship between FP and other ESs closely mirrored the results of the Spearman correlation coefficient. This concurrence lends robust support to the findings of our study, underscoring the reliability of the observed spatial patterns in ES interactions.
When devising a management plan of land use and ecosystems, pinpointing the spatial distribution of trade-offs and synergies between ESs is of paramount importance. This precise location identification of ES relationships is essential for targeted interventions and the optimization of land use strategies. Despite the significance of this task, a dearth of research exists that systematically analyzes the spatial heterogeneity of these ESs interactions [36]. Previous research on trade-offs and synergies in the YRD predominantly concentrated on average levels, neglecting the distinctive spatial distribution patterns of these interactions. This paper underscored that the dynamics of trade-offs and synergies for identical pairs of ESs were subject to geographical variation, a conclusion that aligns with the findings of Xue et al. [42]. For example, while there was a pronounced trade-off between FP and SC, the LISA diagram revealed instances of synergistic co-existence in specific regions (central areas), characterized by H-H and L-L associations. This highlighted the importance of considering the spatial dimension when assessing ES interactions to capture the full complexity of their relationships. By comparing the results related to ES relationships across three scales, we found that the interactions were quite similar in spatial distribution. As the research scale expands from the township to the city level, the spatial distribution becomes increasingly simplified, and hence the information reflected at the township scale is more abundant.
Various ESs exhibited distinctive spatial clustering characteristics, which were significantly correlated with the distribution of different land use types. Specifically, the YRD’s wide distribution of arable land in the north and forests in the south resulted in spatial aggregation of ES relationships exhibiting “upper, middle, and lower” characteristics. This spatial feature was manifested in that ESs in the northern and southern YRD showed different correlation aggregation characteristics, and most of the ESs in the north and south of the YRD were dominated by one correlation feature. However, in the central region of the YRD, the correlation of ESs showed insignificant characteristics.
The LISA maps of various ES pairs involving FP-WR, FP-CR, FP-EP, FP-SC, FP-BM, and FP-AL displayed extensive areas of L-H and H-L (trade-offs) relationships, consistent with the findings of the Spearman correlation coefficient and the BGM in this study. The results regarding the relationship between FP-WR and FP-AL at the township level were also consistent with those of Chen [15]. Furthermore, each ES pair was dominated by two aggregation characteristics, while the other two were dispersed around them, similar to the findings of research conducted in the source area of the three rivers [37].

4.3. Research Limitations and Prospects

This study has delivered notable insights into the scale effects, temporal variations, and spatial heterogeneity of ES trade-offs/synergies; however, certain limitations persist. Firstly, the method of value equivalent factor per unit area simplifies land into various types while ignoring the internal differences of land use types under different natural backgrounds. Although primary and secondary corrections were made in this study, differences in climate, vegetation, and soil quality within the same land use type could not be detailed. Additionally, this study roughly assumed that the ESV of construction land was zero, neglecting the effects of street greening and small urban parks in cities and towns. Secondly, while the Spearman correlation coefficient and bivariate Moran’s index were utilized to quantify trade-offs/synergies and spatial clustering of ESs, the driving mechanisms of interaction between these ESs were not explored [43]. Therefore, there is a need to explore more driving factors of the dependencies and link them to government decisions [44]. Thirdly, the limitation of the analysis is its concentration on the subjective valuation of ecosystems, overlooking their essential objective existence. It is essential to implement an objective assessment that recognizes the intrinsic significance of nature for human life and well-being and effectively conveys the intrinsic and functional value of nature. Lastly, future research should include more spatial scales to provide more comprehensive and scientific information for decision-makers to regulate ecosystems and land use management.

5. Conclusions

This study evaluated the values of seven typical ESs in the YRD using the equivalent value method based on land use data from 2000 to 2020. We examined the time variation, spatial scale effects, and spatial heterogeneity of the ES relationships. The findings revealed that the value of the seven ESs and the total ESV exhibited a fluctuating upward trend. The total ESV increased by RMB 650.625 billion, with the highest growth rate (59%) occurring between 2005 and 2010. Approximately 70% of the ES pairs showed synergies. The direction and degree of trade-offs/synergies between ESs were similar at the grid, township, county, and city scales, whereas the grid scale revealed slight differences. Moreover, the relationship between FP and other ESs showed trade-offs mainly in the northern and southern areas of the YRD, while the other ES pairs exhibited mainly synergy. The spatial distribution results of the relationships between ESs were similar to those of the Spearman correlation analysis and the spatial aggregation of the relationship between ESs exhibiting “upper, middle, and lower” characteristics. This study provides scientific evidence to support local governments in formulating sustainable land and ecosystem management policies.

Author Contributions

Conceptualization, W.L. and F.Z.; methodology, Y.C.; software, Q.Z.; validation, Y.C., W.L. and F.Z.; formal analysis, Z.X.; investigation, W.L. and F.Z.; resources, F.Z.; data curation, Y.C.; writing—original draft preparation, Y.C.; writing—review and editing, W.L. and M.A.K.; visualization, Y.C.; supervision, W.L. and F.Z.; project administration, W.L.; funding acquisition, W.L. and F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the research funds from the National Natural Science Foundation of China, grant number 72104130, 42201046, U2039201; Jinan City-School Integration Project, grant number JNSX2023036; 2024 Experimental Teaching Reform Research Project of Shandong Normal University, grant number 2024MS22.

Data Availability Statement

Some or all data, models, or code generated or used during the study are available from the corresponding author by request.

Acknowledgments

The land use dataset is provided by Resource and Environmental Science Dara Platform (https://www.resdc.cn/).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Ecosystem service value per unit land area of YRD (unit: RMB 10,000/km2).
Table A1. Ecosystem service value per unit land area of YRD (unit: RMB 10,000/km2).
ES TypeFPWRCREPSCBMALTotal
2000Dry land51.51−99.6211.473.420.24.231.14−27.65
Irrigated cropland32.20.767.252.0120.732.620.7666.33
Broad-leaved forest10.6712.51130.8438.8553.3448.5113.42308.14
Shrubbery6.998.185.1525.7734.6231.68.74200.97
Grassland, dense12.510.2104.8834.6248.3143.8812.16266.55
Grassland, open7.245.9260.7920.1327.9825.567.09154.71
Wetland5.5228.0372.4772.4746.5158.4259.89443.30
Water bodies56.19582.2746.1111.7218.7251.3323.93890.26
2005Dry land65.57−126.817.475.210.316.442.49−29.31
Irrigated cropland40.980.9611.033.0731.573.981.6693.25
Broad-leaved forest15.3818.03199.2359.1681.2273.8729.33476.22
Shrubbery10.0711.66129.6539.2352.7248.1219.09310.54
Grassland, dense14.6711.97159.6952.7273.5666.8226.56405.99
Grassland, open8.496.9592.5730.6542.638.9315.49235.68
Wetland8.4542.91110.34110.3470.8241.22130.87714.93
Water bodies77.67804.8670.19170.1128.5178.1652.291281.79
2010Dry land95.32−184.3327.538.210.4810.145.06−37.59
Irrigated cropland59.581.417.394.8349.756.283.37142.60
Broad-leaved forest22.0825.88313.9593.22128116.459.55759.08
Shrubbery14.4616.75204.3161.8283.0875.8338.76495.01
Grassland, dense22.5418.39251.6583.08115.92105.353.93650.81
Grassland, open13.0510.68145.8748.367.1461.3431.46377.84
Wetland13.3967.99173.88173.88111.57380.13265.711186.55
Water bodies116.341205.59110.61268.0744.92123.17106.171974.87
2015Dry land119.55−231.236.1610.780.6313.323.07−47.69
Irrigated cropland74.721.7622.846.3465.338.252.05181.29
Broad-leaved forest27.7932.58412.3122.42168.09152.8736.18952.23
Shrubbery18.2121.08268.3181.19109.199.5923.55621.03
Grassland, dense28.323.08330.48109.1152.24138.2832.76814.24
Grassland, open16.3813.4191.5663.4388.1780.5619.11472.61
Wetland17.5889.29228.35228.35146.53499.21161.431370.74
Water bodies145.341506.09145.26352.0458.99161.7564.52433.97
2020Dry land97.75−189.0333.5910.020.5912.372.54−32.17
Irrigated cropland61.091.4421.215.8960.697.661.69159.67
Broad-leaved forest23.7527.85383.01113.72156.15142.0129.86876.35
Shrubbery15.5618.02249.2575.42101.3592.5119.44571.55
Grassland, dense20.316.56307101.35141.42128.4627.04742.13
Grassland, open11.759.62177.9558.9281.974.8315.77430.74
Wetland16.5183.85212.13212.13136.12463.73133.231257.70
Water bodies144.531497.66134.94327.0354.8150.2653.242362.46
Note: FP: Food Production, WR: Water Resources Supply, CR: Climate Regulation, EP: Environmental Purification, SC: Soil Conservation, BM: Biodiversity Maintenance, AL: Aesthetic Landscape.
Figure A1. LISA cluster map of trade-offs/synergies between Ecosystem Services at county level in the Yangtze River Delta in 2020.
Figure A1. LISA cluster map of trade-offs/synergies between Ecosystem Services at county level in the Yangtze River Delta in 2020.
Land 13 01462 g0a1aLand 13 01462 g0a1b
Figure A2. LISA cluster map of trade-offs/synergies between Ecosystem services at city level in the Yangtze River Delta in 2020.
Figure A2. LISA cluster map of trade-offs/synergies between Ecosystem services at city level in the Yangtze River Delta in 2020.
Land 13 01462 g0a2aLand 13 01462 g0a2b

References

  1. Qiu, S.; Peng, J.; Dong, J.; Wang, X.; Meersmans, J. Understanding the relationships between ecosystem services and associated social-ecological drivers in a karst region: A case study of Guizhou province, China. Prog. Phys. Geogr. 2020, 45, 98–114. [Google Scholar] [CrossRef]
  2. Magdalena, U.R.; Gonçalves De Souza, G.B.; Amorim, R.R. Spatial analysis guiding decision making in environmental conservation: Systematic conservation planning and ecosystem services. Prog. Phys. Geogr. 2022, 47, 123–139. [Google Scholar] [CrossRef]
  3. 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. Landscape Urban Plan. 2014, 132, 102–110. [Google Scholar] [CrossRef]
  4. Dunford, R.W.; Smith, A.C.; Harrison, P.A.; Hanganu, D. Ecosystem service provision in a changing europe: Adapting to the impacts of combined climate and socio-economic change. Landscape Ecol. 2015, 30, 443–461. [Google Scholar] [CrossRef] [PubMed]
  5. Estoque, R.C.; Murayama, Y. Landscape pattern and ecosystem service value changes: Implications for environmental sustainability planning for the rapidly urbanizing summer capital of the Philippines. Landscape Urban Plan. 2013, 116, 60–72. [Google Scholar] [CrossRef]
  6. Geng, W.; Li, Y.; Zhang, P.; Yang, D.; Jing, W.; Rong, T. Analyzing spatio-temporal changes and trade-offs/synergies among ecosystem services in the Yellow river basin, China. Ecol. Indic. 2022, 138, 108825. [Google Scholar] [CrossRef]
  7. Zhao, X.; Wang, J.; Su, J.; Sun, W. Ecosystem service value evaluation method in a complex ecological environment: A case study of Gansu Province, China. PLoS ONE 2021, 16, e240272. [Google Scholar] [CrossRef]
  8. Torres, A.V.; Tiwari, C.; Atkinson, S.F. Progress in ecosystem services research: A guide for scholars and practitioners. Ecosyst. Serv. 2021, 49, 101267. [Google Scholar] [CrossRef]
  9. Wang, X.; Peng, J.; Luo, Y.; Qiu, S.; Dong, J.; Zhang, Z.; Vercruysse, K.; Grabowski, R.C.; Meersmans, J.; Cleveland, C.J. Exploring social-ecological impacts on trade-offs and synergies among ecosystem services. Ecol. Econ. 2022, 197, 107438. [Google Scholar] [CrossRef]
  10. Gong, J.; Xu, C.; Yan, L.; Zhu, Y.; Zhang, Y.; Jin, T. Multi-scale analysis of ecosystem services trade-offs in an ecotone in the eastern margin of the Qinghai-tibetan Plateau. J. Mt. Sci.-Engl. 2021, 18, 2803–2819. [Google Scholar] [CrossRef]
  11. Aryal, K.; Maraseni, T.; Apan, A. How much do we know about trade-offs in ecosystem services? A systematic review of empirical research observations. Sci. Total Environ. 2022, 806, 151229. [Google Scholar] [CrossRef]
  12. Huang, L.; Du, Y.; Tang, Y. Ecosystem service trade-offs and synergies and their drivers in severely affected areas of the Wenchuan earthquake, China. Land. Degrad. Dev. 2024, 35, 3881–3896. [Google Scholar] [CrossRef]
  13. Liu, J.; Pei, X.; Zhu, W.; Jiao, J. Scenario modeling of ecosystem service trade-offs and bundles in a semi-arid valley basin. Sci. Total Environ. 2023, 896, 166413. [Google Scholar] [CrossRef]
  14. Qiao, X.; Gu, Y.; Zou, C.; Xu, D.; Wang, L.; Ye, X.; Yang, Y.; Huang, X. Temporal variation and spatial scale dependency of the trade-offs and synergies among multiple ecosystem services in the Taihu Lake Basin of China. Sci. Total Environ. 2019, 651, 218–229. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, H.; Yan, W.; Li, Z.; Wende, W.; Xiao, S.; Wan, S.; Li, S. Spatial patterns of associations among ecosystem services across different spatial scales in metropolitan areas: A case study of Shanghai, China. Ecol. Indic. 2022, 136, 108682. [Google Scholar] [CrossRef]
  16. Petz, K.; Alkemade, R.; Bakkenes, M.; Schulp, C.J.E.; van der Velde, M.; Leemans, R. Mapping and modelling trade-offs and synergies between grazing intensity and ecosystem services in rangelands using global-scale datasets and models. Glob. Environ. Chang. 2014, 29, 223–234. [Google Scholar] [CrossRef]
  17. Armatas, C.A.; Campbell, R.M.; Watson, A.E.; Borrie, W.T.; Neal, C.; Venn, T.J. An integrated approach to valuation and tradeoff analysis of ecosystem services for national forest decision-making. Ecosyst. Serv. 2018, 33, 1–18. [Google Scholar] [CrossRef]
  18. Karimi, J.D.; Harris, J.A.; Corstanje, R. Using Bayesian Belief Networks to assess the influence of landscape connectivity on ecosystem service trade-offs and synergies in urban landscapes in the UK. Landscape Ecol. 2021, 36, 3345–3363. [Google Scholar] [CrossRef]
  19. Su, C.; Dong, M.; Fu, B.; Liu, G. Scale effects of sediment retention, water yield, and net primary production: A case-study of the Chinese Loess Plateau. Land Degrad. Dev. 2020, 31, 1408–1421. [Google Scholar] [CrossRef]
  20. Xiong, L.; Li, R. Assessing and decoupling ecosystem services evolution in karst areas: A multi-model approach to support land management decision-making. J. Environ. Manag. 2024, 350, 119632. [Google Scholar] [CrossRef]
  21. Deng, X.; Xiong, K.; Yu, Y.; Zhang, S.; Kong, L.; Zhang, Y. A review of ecosystem service trade-offs/synergies: Enlightenment for the optimization of forest ecosystem functions in karst desertification control. Forests 2023, 14, 88. [Google Scholar] [CrossRef]
  22. Hou, Y.; Lü, Y.; Chen, W.; Fu, B. Temporal variation and spatial scale dependency of ecosystem service interactions: A case study on the central loess plateau of China. Landscape Ecol. 2017, 32, 1201–1217. [Google Scholar] [CrossRef]
  23. Rodríguez, J.P.; Beard, J.T.D.; Bennett, E.M.; Cumming, G.S.; Cork, S.J.; Agard, J.; Dobson, A.P.; Peterson, G.D. Trade-offs across space, time, and ecosystem services. Ecol. Soc. 2006, 11, 28. [Google Scholar] [CrossRef]
  24. Li, B.; Wang, W. Trade-offs and synergies in ecosystem services for the Yinchuan basin in China. Ecol. Indic. 2018, 84, 837–846. [Google Scholar] [CrossRef]
  25. Zhang, B.; Li, W.; Xie, G. Ecosystem services research in China: Progress and perspective. Ecol. Econ. 2010, 69, 1389–1395. [Google Scholar] [CrossRef]
  26. Yang, M.; Gao, X.; Zhao, X.; Wu, P. Scale effect and spatially explicit drivers of interactions between ecosystem services—A case study from the Loess Plateau. Sci. Total Environ. 2021, 785, 147389. [Google Scholar] [CrossRef]
  27. Pan, J.; Wei, S.; Li, Z. Spatiotemporal pattern of trade-offs and synergistic relationships among multiple ecosystem services in an arid inland river basin in NW China. Ecol. Indic. 2020, 114, 106345. [Google Scholar] [CrossRef]
  28. Li, R.; Kong, L.; Yang, Y.; Wang, Y.; Zheng, H.; Liang, M. Dynamic bundles to detect the spatiotemporal characteristics and impact factors of ecosystem services in northern China. Prog. Phys. Geogr. 2023, 47, 687–701. [Google Scholar] [CrossRef]
  29. Tao, Y.; Tao, Q.; Sun, X.; Qiu, J.; Pueppke, S.G.; Ou, W.; Guo, J.; Qi, J. Mapping ecosystem service supply and demand dynamics under rapid urban expansion: A case study in the Yangtze River Delta of China. Ecosyst. Serv. 2022, 56, 101448. [Google Scholar] [CrossRef]
  30. Cai, W.; Gibbs, D.; Zhang, L.; Ferrier, G.; Cai, Y. Identifying hotspots and management of critical ecosystem services in rapidly urbanizing Yangtze river delta region, China. J. Environ. Manag. 2017, 191, 258–267. [Google Scholar] [CrossRef]
  31. Shu, H.; Xiong, P. Reallocation planning of urban industrial land for structure optimization and emission reduction: A practical analysis of urban agglomeration in China’s Rangtze River Delta. Land. Use Policy 2019, 81, 604–623. [Google Scholar] [CrossRef]
  32. Xie, G.; Zhang, C.; Zhang, L.; Chen, W.; Li, S. Improvement of the evaluation method for ecosystem service value based on per unit area. J. Nat. Resour. 2015, 30, 1243–1254. [Google Scholar]
  33. Luo, Q.; Zhang, X.; Li, Z.; Yang, M.; Lin, Y. The effects of China’s ecological control line policy on ecosystem services: The case of Wuhan city. Ecol. Indic. 2018, 93, 292–301. [Google Scholar] [CrossRef]
  34. Wu, X.; Wang, S.; Fu, B.; Liu, Y.; Zhu, Y. Land use optimization based on ecosystem service assessment: A case study in the Yanhe watershed. Land. Use Policy 2018, 72, 303–312. [Google Scholar] [CrossRef]
  35. Xie, G.; Zhen, L.; Lu, C.; Xiao, Y.; Li, W. Applying value transfer method for eco-service valuation in China. J. Resour. Ecol. 2010, 1, 51–59. [Google Scholar]
  36. Xu, S.; Liu, Y. Associations among ecosystem services from local perspectives. Sci. Total Environ. 2019, 690, 790–798. [Google Scholar] [CrossRef]
  37. Zheng, D.; Wang, Y.; Hao, S.; Xu, W.; Lv, L.; Yu, S. Spatial-temporal variation and tradeoffs/synergies analysis on multiple ecosystem services: A case study in the three-river headwaters region of China. Ecol. Indic. 2020, 116, 106494. [Google Scholar] [CrossRef]
  38. Tian, Y.; Jiang, G.; Zhou, D.; Li, G. Systematically addressing the heterogeneity in the response of ecosystem services to agricultural modernization, industrialization and urbanization in the Qinghai-Tibetan plateau from 2000 to 2018. J. Clean. Prod. 2021, 285, 125323. [Google Scholar] [CrossRef]
  39. Felipe-Lucia, M.R.; Comín, F.A.; Bennett, E.M. Interactions among ecosystem services across land uses in a floodplain agroecosystem. Ecol. Soc. 2014, 19, 20. [Google Scholar] [CrossRef]
  40. Li, S.; Zhao, Y.; Xiao, W.; Yellishetty, M.; Yang, D. Identifying ecosystem service bundles and the spatiotemporal characteristics of trade-offs and synergies in coal mining areas with a high groundwater table. Sci. Total Environ. 2022, 807, 151036. [Google Scholar] [CrossRef]
  41. Yang, Y.; Li, M.; Feng, X.; Yan, H.; Su, M.; Wu, M. Spatiotemporal variation of essential ecosystem services and their trade-off/synergy along with rapid urbanization in the lower pearl river basin, China. Ecol. Indic. 2021, 133, 108439. [Google Scholar] [CrossRef]
  42. Xue, C.; Chen, X.; Xue, L.; Zhang, H.; Chen, J.; Li, D. Modeling the spatially heterogeneous relationships between tradeoffs and synergies among ecosystem services and potential drivers considering geographic scale in bairin left banner, China. Sci. Total Environ. 2023, 855, 158834. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, W.; Zhan, J.; Zhao, F.; Wang, C.; Zhang, F.; Teng, Y.; Chu, X.; Kumi, M.A. Spatio-temporal variations of ecosystem services and their drivers in the pearl river delta, China. J. Clean. Prod. 2022, 337, 130466. [Google Scholar] [CrossRef]
  44. Wang, H.; Liu, L.; Yin, L.; Shen, J.; Li, S. Exploring the complex relationships and drivers of ecosystem services across different geomorphological types in the Beijing-Tianjin-Hebei region, China (2000–2018). Ecol. Indic. 2021, 121, 107116. [Google Scholar] [CrossRef]
Figure 1. Location and zoning (a), DEM model (b), and land use in 2000 (c) and 2020 (d) of Yangtze River Delta.
Figure 1. Location and zoning (a), DEM model (b), and land use in 2000 (c) and 2020 (d) of Yangtze River Delta.
Land 13 01462 g001
Figure 2. Change trends in the value of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Figure 2. Change trends in the value of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Land 13 01462 g002
Figure 3. Spatial distribution of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Figure 3. Spatial distribution of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Land 13 01462 g003aLand 13 01462 g003b
Figure 4. Trade-offs/synergies between ecosystem services under different spatial scales in Yangtze River Delta between 2000 and 2020 (Note: ** and * denote statistical significance at the 0.01 and 0.05 probability levels, respectively).
Figure 4. Trade-offs/synergies between ecosystem services under different spatial scales in Yangtze River Delta between 2000 and 2020 (Note: ** and * denote statistical significance at the 0.01 and 0.05 probability levels, respectively).
Land 13 01462 g004
Figure 5. Bivariate global Moran’s I of trade-offs/synergies between ecosystem services at township level in the Yangtze River Delta in 2020.
Figure 5. Bivariate global Moran’s I of trade-offs/synergies between ecosystem services at township level in the Yangtze River Delta in 2020.
Land 13 01462 g005
Figure 6. LISA cluster map of trade-offs/synergies between ESs at township level in the Yangtze River Delta in 2020.
Figure 6. LISA cluster map of trade-offs/synergies between ESs at township level in the Yangtze River Delta in 2020.
Land 13 01462 g006
Table 1. The adjustment of land use classification in this study.
Table 1. The adjustment of land use classification in this study.
Classification in the Land Use DataLand Use Classification
Xie et al. [32]This Study
Paddy fieldPaddy fieldIrrigated cropland
Dry landDry landDry land
Marsh, saline–alkaline land, beach, bottom landWetlandWetland
Canals, lakes, reservoirs and pondsRiver systemWater bodies
Woodland, sparse woodland, other woodlandsBroad-leaved forestBroad-leaved forest
ShrubberyShrubberyShrubbery
High0cover grasslandBushGrassland, dense
Low- and medium-cover grasslandMeadowGrassland, open
Sandy land, bare land, bare rock stony landBare landBarren or sparsely vegetated
Table 2. Ecosystem service equivalent value per unit area of land use in China.
Table 2. Ecosystem service equivalent value per unit area of land use in China.
Land UseFPWRCREPSCBMAL
Arable landDry land0.850.020.360.11.030.130.06
Irrigated cropland1.36−2.630.570.170.010.210.09
Forest landBroad-leaved forest0.290.346.51.932.652.411.06
Shrubbery0.190.224.231.281.721.570.69
GrasslandDense grassland 0.380.315.211.722.42.180.96
Open grassland0.220.183.0211.391.270.56
Water bodiesWetland0.512.593.63.62.317.874.73
Water bodies0.88.292.295.550.932.551.89
Note: FP: Food Production, WR: Water Resources Supply, CR: Climate Regulation, EP: Environmental Purification, SC: Soil Conservation, BM: Biodiversity Maintenance, AL: Aesthetic Landscape. Source: Xie et al. [32].
Table 3. The functional adjustment index of equivalent value factors of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Table 3. The functional adjustment index of equivalent value factors of different ecosystem services in Yangtze River Delta between 2000 and 2020.
Correction Factor20002005201020152020
FP of arable land3.502.912.672.552.22
FP of forest land3.403.202.902.782.53
FP of grassland 3.042.332.262.161.65
FP of water bodies6.495.865.545.275.58
CR, EP, SC, BM1.861.851.841.841.82
AL1.171.672.140.990.87
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

Chen, Y.; Liu, W.; Zhao, F.; Zhao, Q.; Xu, Z.; Asiedu Kumi, M. Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies in the Yangtze River Delta. Land 2024, 13, 1462. https://doi.org/10.3390/land13091462

AMA Style

Chen Y, Liu W, Zhao F, Zhao Q, Xu Z, Asiedu Kumi M. Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies in the Yangtze River Delta. Land. 2024; 13(9):1462. https://doi.org/10.3390/land13091462

Chicago/Turabian Style

Chen, Yongqi, Wei Liu, Fen Zhao, Qing Zhao, Zhiwei Xu, and Michael Asiedu Kumi. 2024. "Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies in the Yangtze River Delta" Land 13, no. 9: 1462. https://doi.org/10.3390/land13091462

APA Style

Chen, Y., Liu, W., Zhao, F., Zhao, Q., Xu, Z., & Asiedu Kumi, M. (2024). Multi-Scale Analysis of Ecosystem Service Trade-Offs/Synergies in the Yangtze River Delta. Land, 13(9), 1462. https://doi.org/10.3390/land13091462

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