Next Article in Journal
Pathologic Complete Response (pCR) and Survival of Women with Inflammatory Breast Cancer (IBC): An Analysis Based on Biologic Subtypes and Demographic Characteristics
Next Article in Special Issue
An Evaluation Model for Urban Comprehensive Carrying Capacity: An Empirical Case from Harbin City
Previous Article in Journal
A Toxicogenomic Approach Reveals a Novel Gene Regulatory Network Active in In Vitro and In Vivo Models of Thyroid Carcinogenesis
Previous Article in Special Issue
Evolution of the Pattern of Spatial Expansion of Urban Land Use in the Poyang Lake Ecological Economic Zone
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Land Use/Cover Change on Yangtze River Delta Urban Agglomeration Ecosystem Services Value: Temporal-Spatial Patterns and Cold/Hot Spots Ecosystem Services Value Change Brought by Urbanization

1
College of Management Science, Chengdu University of Technology, No.1 Dongsan Road, Erxian Bridge, Chenghua District, Chengdu 610059, China
2
Key Laboratory of GeoSpatial Information Technology of Ministry of Land and Resources, Chengdu University of Technology, No.1 Dongsan Road, Erxian Bridge, Chenghua District, Chengdu 610059, China
3
College of Earth Sciences, Chengdu University of Technology, No.1 Dongsan Road, Erxian Bridge, Chenghua District, Chengdu 610059, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(1), 123; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010123
Submission received: 21 November 2018 / Revised: 18 December 2018 / Accepted: 29 December 2018 / Published: 4 January 2019

Abstract

:
Land use/cover change (LUCC) from increased urbanization significantly impacts regional ecosystem services. Based on a cold/hot spots analysis, this paper used grain yield, food prices, price index statistics, and a land use thematic map to study the impact of LUCC on four ecosystem services values (ESVs) in the Yangtze River Delta urban agglomeration, and determine the spatial differences resulting from the rapid urbanization LUCC. The correlation between the four ecosystem services was then studied and sensitivity analyses conducted to investigate whether any changes in the ESVs could lead to unacceptable unit value transfer uncertainties. It was found that most urban land was converted from farmland, and that before 2000, the total ESVs and the regulating services values (RSVs) increased significantly, after which it declined, the provisioning services values (PSVs) declined year on year, the habitat services value (HSV) and cultural and amenity services value (CSV) declined sharply after 2000, and the spatial distribution of the four ESVs were significantly different. Over time, it was found that the hot spots were shrinking and the cold spots were spreading. The provisioning services were found to be negatively correlated with habitat services and cultural and amenity services, the regulating services were weakly positively correlated with the provisioning services and significantly positively correlated with the habitat services and cultural and amenity services, and the habitat services were significantly positively correlated with cultural and amenity services. In the Yangtze River Delta urban agglomeration, the water area is the most important for the total ESVs, followed by non-bush forest. Paddy field is ranked third. Dryland, bush, grassland, and wetland are less important. The importance of barren land is almost zero. This research provides the government with a scientific basis from which to formulate spatial planning and environmental protection policies for ecological sustainable development in the Yangtze River Delta urban agglomeration.

1. Introduction

LUCC (Land use/cover change) programs involve the study of land surface changes, and were inspired by the NASA (National Aeronautics and Space Administration) Earth Science Program [1]. As land is the most essential ecosystem element, LUCC can have a significant impact on regional ESVs (ecosystem services values) [2,3,4,5]. The physical aspects of land (biological), such as the rocks, soil, vegetation, and water, are referred to as land cover, and the human uses, such as arrangements, activities, and inputs, are referred to as land use [6,7]. Land use changes inevitably affect the ecosystem, biodiversity, as well as communities [8]. Land cover can change because of natural phenomena, such as weathering, glaciers, and vegetation succession; natural disasters, such as earthquakes and mudslides; and human activities [9]. When these changes occur, the land can return to equilibrium; however, some land restoration efforts, such as recovering from desertification, can take a long time [10] and require significant resources. Land use/land cover changes can cause significant environmental impacts, such as soil degradation, deforestation, a loss of biodiversity, and declines in both water quantity and water quality [11], all of which affect the ability of the ecosystems to provide the natural products and services critical to human survival and well-being [12]. Consequently, the ESVs has become an important indicator to assess whether a region is able to attain sustainable development [13]. It has, therefore, become important to identify, quantify, and evaluate all ecosystem services changes [14], with quantitative assessments of land-use changes being the main approach [15,16]. Some scholars have been advocated for integrating ecosystem services into land use decisions [17,18,19], and it can help to promote interaction between researchers and policymakers to develop more effective policies to improve the ecological environment, but there are few methods and demonstrations in practice. So, Liang et al. have developed a method, integrating ESVs to balance future ecosystem-service benefit and risk to optimize investment in land for ecological conservation in land use planning [19].
The ESVs concept, which assesses the ecosystem-service benefits that humans receive from natural systems, was first proposed by Costanza et al. The ecosystem services economic value is determined from the interactions between the ecosystem supply and the human/social needs [15]. As natural systems are vitally important to economic well-being, the ecosystem value includes an assessment of the many ways human and natural systems interact, with the goal being to ensure that “nature” is fully reflected in government decision-making processes [20]. The benefit transfer assessment method uses existing valuation studies or data to estimate the ESVs at a site and then transfers these results to ESVs at similar locations [21]. In 1997, Constanza et al. took the lead and estimated the system functional value of the global biosphere ecosystem services. The publication of the results in the Nature magazine attracted widespread attention and prompted increased research interest in ecosystem services economic values assessments. Since that time, there has been significant research published on the impact of land use change on the ESVs on global scales [22], national scales, [23], and regional scales [24], with many watershed scale studies [25,26,27] as well as multi-scale comparative studies [28]. Research has also focused on assessing individual ESVs for forests [29,30,31], grasslands [32], farmland [33,34], wetlands [35,36], and marine environments [37]. Therefore, because of these extensive studies, there have also been many assessment methods developed, such as cost-benefit analysis [38], ecosystem services willingness to pay assessments [39], value transfers [21,40], public participation in GIS (PPGIS) methods [41], selection experiments [42], conditional valuations [43], and multiple or mixed methods (group reviews, semi-structured interviews, and analytical techniques, such as factor analysis) [44]. Li et al. estimated the impact of land-based land changes on the ESVs using the product of the area and a coefficient to reflect the land use impact [45]. Fei et al. revised the ESVs calculation method based on biophysical and socio-economic factors, and reflected the land use change impact on the ESVs in terms of space and quantity [46]. Wang et al. studied the driving factors behind land use change and the ESVs responses in Hengduan Mountains along the elevation gradient [47]. Cao et al. analyzed the ESVs gains and losses resulting from the rapid urbanization of the coastal areas of Zhejiang Province from a land use change perspective [48]. Arowolo et al. assessed the ESVs changes based on the land use/land cover dynamics in different administrative regions in Nigeria [49]. Yu et al. quantitatively evaluated the ESVs and the influence of human activities in the upper reaches of the Yangtze River, and proposed related planning methods [50]. Most previous research has quantitatively assessed land use change impacts on the ESVs. However, other studies have examined ESVs changes from the perspective of LUCC [48], have integrated ecosystem services into land use impact assessment methods to differentiate the different land use types [51], and have analyzed regional spatial differences [23]. However, few studies have considered the spatial heterogeneity of the LUCC impact on the various ESVs types and have offered a visual representation of these spatial changes.
Rapid economic development and urban expansion has resulted in significant changes in land use/cover, which has in turn resulted in changes in the ecosystem structures around cities, reducing the ESVs [52]. In particular, the growth in the Yangtze River Delta urban agglomeration in China has been at the expense of farmland, with the urban ecological function significantly weakening and threatening the overall stability of the ecological environment [53]. In recent years, governments at all levels have sought to control the decline in the Yangtze River Delta urban agglomeration ecological environment. However, as no quantitative assessments have been conducted to assess the advantages and disadvantages of the ecological environmental changes or the accompanying regional changes due to the growing urban sprawl, the ecosystem services trends and spatial patterns remain unclear.
This paper takes the Yangtze River Delta urban agglomeration as a research area to explore the potential impact of the LUCC on four ESVs; provisioning services, regulating services, habitat services, and cultural and amenity services; and to assess their spatial differences and correlations to determine whether the spatial land use/land cover layout is balanced. The supplier valuation method was used to improve Xie’s [54] ESVs coefficient and to evaluate the four ESVs, after which a sensitivity analysis was conducted. Using the cold/hot spots analysis method in the ArcGIS spatial statistical tool (ESRI, Redlands CA, USA, available online: http://www.esri.com), a cold/hot spots spatial distribution map for the four ESVs was generated. From a regional perspective, changes in the ESVs in 26 cities in the Yangtze River Delta were analyzed, after which they were clustered into high and low ESVs regions to clearly show the regional imbalances in the spatial land use/land cover patterns. The research results provide decision support for the government to balance urban development and ecological environmental protection, and provide guidance for the development of ecological environmental spatial management and protection plans for the Yangtze River Delta.

2. Materials and Methods

2.1. Study Area and Data Source

The Yangtze River Delta urban agglomeration (Figure 1), located in the Yangtze River Economic Belt, is part of China’s latest stage of “reform and opening up”: The transformation and implementation of the new regional open development strategy. The area belongs to the alluvial plain and has a long history of agriculture. It is one of the most developed and densely populated industrial areas in China [55]. According to the “Changjiang Delta Urban Agglomeration Development Plan” approved by the State Council, the administrative scope of the Yangtze River Delta urban agglomeration includes Shanghai, Jiangsu Province (Nanjing City, Wuxi City, Changzhou City, Suzhou City, Nantong City, Yancheng City, Yangzhou City, Zhenjiang City, Taizhou City), Zhejiang Province (Hangzhou City, Ningbo City, Jiaxing City, Huzhou City, Shaoxing City, Jinhua City, Zhoushan City, Taizhou City), and Anhui Province (Hefei City, Wuhu City, Chaohu City, Maanshan City, Tongling City, Anqing City, Zhangzhou City, Chizhou City, Xuancheng City). It covers an area of 210,000 km2 with a population of 139.2 million and an average population density of 581 people per square kilometer. Shanghai is the core city of the Yangtze River Delta urban agglomeration.
The statistics on grain yield and price index used in this study are derived from the China Statistical Yearbook. The grain price data are from the government website of the National Development and Reform Commission. The land use data (Figure 2) are from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences [56].

2.2. The Assessment Method of ESVs

To date, the methods used to evaluate the ESVs in studies within the academic community have not been consistent [15]. The broader ESVs valuation method of “supplier valuation” was applied in the present study. The analysis presented in this paper draws on Xie’s [54] supplier valuation method and uses the product of the standard equivalence factor and the equivalent coefficient to determine the ESVs coefficient. The standard equivalence factor, as a benchmark for other ecosystem services, refers to the value provided by a natural ecological component of an ecosystem per unit of land area [54]. The standard equivalence factor is defined as the economic value of the grain yield of 1 hectare of farmland, which is obtained by multiplying the national grain yield and the national food price average during the study period. Adopting the classification of ecosystem and biodiversity economics [57], the ecosystem services are divided into four categories: Provisioning services, regulating services, habitat services, and cultural and amenity services [38]. Considering the representativeness of ecosystem services and its particular importance to the research field, according to present data, the existing value coefficient of different ecological land types in China further categorizes ecosystem services into food supply, water supply, material supply, climate regulation, air quality regulation, wastewater treatment, regulation of water flow, erosion protection, soil fertility maintenance, biodiversity, and aesthetic landscape. Eleven aspects were assessed for the ESVs generated in the Yangtze River Delta urban agglomeration. Since the research area of this study is an urban agglomeration, the research period extends from 1980 to 2015, so the timespan is long. Therefore, national-level statistical data were used for calculations, and the grain price used in the analysis is the average of the values during the study period. Taking into account the differences in food prices, a normalized food price index was used to correct the standard equivalence factor (Table 1):
S c = N P I × Y × P
C i j = S c × E i j
ESV = i = 1 n j = 1 m C i j × A i
where S c is the standard equivalence factor; N P I is the normalized price index; Y is the average grain yield per unit (kg/ha) at the national level for the research period; P is the national average grain price per unit in CNY/kg; C i j is the value coefficient of the i-th land use type’s j-th ecosystem service; E i j is the equivalent coefficient of the i-th land use type’s j-th ecosystem services; A i is the area of the i-th land use type in hectares.
According to the Classification Standards for Land Use Types, land use types are divided into seven first-class land types, and there are six types of ESVs per unit area. According to the available land use data, farmland and forests are divided into subcategories: Paddy field and dry land for farmland; bush and non-bush for forests.
In this study, we introduced a simplified version of Aschonitis et al.’s method to establish the importance of each land use type based on its contribution to the total ESVs [58]; the importance values were then ranked:
C S i = C b i × A i ESV b
where C S i is the sensitivity factor of the i-th land use type; ESV b is the total ESVs; C b i is the value coefficient of the i-th land type; A i is the area of the i-th land use type in hectares.

2.3. Global Spatial Autocorrelation

Global Moran’s I measures spatial similarity and was used to describe the ESVs’ distribution pattern of different units in the study area, i.e., spatial autocorrelation [59]. The spatial relevance of the ESVs of the Yangtze River Delta urban agglomeration can be measured by the Global Moran’s I Index [60]. The calculation formula for Global Moran’s I is as follows:
I t = n i j w i , j z i z j S i z i 2
where z i and z j are the deviations of the k-th ecosystem service value (from the space unit i, j, respectively) from its mean, where k = 1, 2, 3, 4; w i , j is the spatial weight between spatial units i and j; S is the sum of all the elements of the weight matrix; and n is the total number of spatial units. The range for Global Moran’s I index is [−1, 1]. The exponent with this variable greater than 0 presents an overall positive spatial correlation; if it is less than 0, it means there is a negative correlation; if it is equal to 0, it means that there is no correlation. The importance of each calculated Global Moran’s I index can be tested by Z statistics [60].

2.4. Cold/Hot Spots Analysis

In order to detect the high and low values and the degree of clustering of potential ESVs in different regions, the Getis–Ord G i * statistic was calculated for different ESVs using the spatial statistical tool in ArcGIS, the cold/hot spots analysis tool. The higher the Z-score of the G i * index, the closer the aggregation of high values (hot spots) for a given ecosystem service [61] and the higher the values of the attributes around the unit. The lower the score, the tighter the aggregation of low values (cold spots), the lower the value of the attributes around the unit, and the lower the ecosystem services supply value [62]. The G i * index calculation formula is as follows:
G i * = j = 1 n w i , j x j X ¯ j = 1 n w i , j S [ n j = 1 n w i , j 2 ( j = 1 n w i , j ) 2 ] n 1
where x j is the value of the k-th ecosystem service of space unit j; w i , j is the spatial weight between spatial units i and j; n is the total number of spatial units.

3. Results and Analysis

3.1. Analysis of the Changes in the ESVs Resulting from LUCC

During the study period, land use/land cover underwent major changes (Figure 3). From 1980 to 2015, farmland (paddy field and dry land), grassland, and wetland were gradually decreasing; the trend of decreasing paddy field area gradually increased and the trend of decreasing grassland area gradually weakened. The area of water and towns gradually increased, and both bush and non-bush forests first increased and then decreased. From 1980 to 2000, the area of bush and non-bush forests increased greatly, while wasteland first decreased and then increased.
Figure 4 shows the ESVs trends in the Yangtze River Delta urban agglomeration at three time points in the period between 1980 and 2015. The PSVs (provisioning services values) were gradually decreasing. It can be seen from the interval between the two time points that the ESVs decreases most during the period from 2000 to 2015. The RSVs (regulating services values) were gradually increasing, the HSV (habitat services value) and the CSV (cultural and amenity services value) first increased and then decreased. Habitat services gradually started decreasing after 2000, which is an indirect indication that the quality of habitats is declining. The total ESVs first increased and then decreased. Before 2000, the ESVs increased significantly, and then it decreased. This is an early warning sign of the destruction of ecological land, and it needs attention.
Comparing Figure 3 and Figure 4, the temporal change trend of forest area (bush and non-bush) is similar to the change trend of the total ESVs, indicating that the change in forest area is the main factor affecting the change in the total ESVs [63]. The PSVs in the decreasing farmland (paddy field and dry land) also decreased, indicating that changes in farmland have a great impact on provisioning services. Regulating services, habitat services, and cultural and amenity services may be affected by a combination of these land use/cover types.
During the period from 1980 to 2000, the farmland (paddy field and dry land) and grassland areas decreased the most: The farmland area decreased by 456,600 hectares, and the ecological service value provided by farmland decreased by about 7.437 billion yuan. Although the wetland only decreased by 55,300 hectares, the value of the ecological services it provides decreased by about 7250.9 million yuan. Although the water area only increased by 55,300 hectares, the value of the ecological services it provides increased by 175.083 billion yuan. The disproportionately large effects of the changes in the wetland and water area are because these land types are important ecological protection barriers and are associated with higher value factors. The area of forests (bush forests and non-bush forests) increased by 261,600 hectares, and the value of its ecological services increased by 14.481 billion yuan. The urban area increased by 366,000 hectares. The total ESVs increased by 11.765 billion yuan, which was mainly due to the increase in forest and water area. From 2000 to 2015, the farmland area decreased by 766,600 hectares, and the ESVs decreased by 14.251 billion yuan, about double that of the previous period. The water area increased by 61,400 hectares, and the urban area increased by 801,400 hectares. During the study period, except for a small increase in the area of wasteland, the area of the remaining land use/cover types was reduced. The value coefficient and the change in the area of the wasteland were small, and the ESVs provided was low. So, the total ESVs decreased by 2.818 billion yuan. During the period of 1980–2015, a large amount of farmland was transformed into towns, causing serious degradation of farmland and severely affecting provisioning services (Table 2).
According to the radar chart, the regional variation in the ESVs is significant (Figure 5). From 1980 to 2000, the largest values’ variation of ecosystem services was in Jinhua City, Hangzhou City, Ningbo City, and Taizhou City, and it was also positive. The ESVs of Huzhou City, Chizhou City, Anqing City, and Zhoushan City increased slightly, and other cities experienced a decrease. From 2000 to 2015, the values’ variation in ecosystem services was the largest in Hangzhou, Xuancheng, and Anqing. The values’ variation of ecosystem services was positive due to a small increase in ESVs in Taizhou, Zhoushan, Shaoxing, Jinhua, Huzhou, Zhangzhou, Chizhou, and Chaohu, while the ESVs in the other cities decreased. There were 14 cities with reduced ecosystem services during the two study periods: Shanghai, Zhenjiang, Yangzhou, Yancheng, Wuxi, Taizhou, Suzhou, Nanjing, Nantong, Changzhou, Wuhu, Tongling, Maanshan, and Hefei. Among them, Shanghai had the largest reduction, and its ESVs continue to decrease.

3.2. Spatial Autocorrelation Analysis and Correlation Analysis of Four Ecosystem Services

It can be seen from Table 3 that the Moran I index of the ESVs in the study area is greater than 0, and the Z-score is higher than the standard 1% significance level value of 1.65. This indicates that the spatial distribution of ESVs in the study area has a strong spatial autocorrelation—it is aggregated in space. The PSVs index I and the Z-score showed a downward trend. The RSVs index I and the Z-score showed an upward trend, and the CSV index I and the Z-score showed an upward trend.
The Pearson correlation coefficients (5% significance level) for the four potential ecosystem services are presented in Table 4. Provisioning services were negatively correlated with habitat services and cultural and amenity services. Regulating services were weakly positively correlated with provisioning services and were significantly positively correlated with habitat services and cultural and amenity services. The correlation coefficients were 0.219, 0.857, and 0.871, respectively. Habitat services were significantly positively correlated with cultural and amenity services; these two services had the highest correlation (0.998).

3.3. Evolution of the Spatial Pattern Of ESVs Based on Cold/Hot Spots

The Getis–Ord G i * statistic [64] is one of the Local Indicators of Spatial Association (LISA) measures [65]. Calculating Getis–Ord G i * statistics can identify hot spots and cold spots of the potential supply of each ecosystem service. Cold/hot spots are areas with spatial clustering features and high/low ecosystem service supply values [62]. The cold/hot spots analysis results enable spatial and visual comparisons among ecosystem services [66]. The cold/hot spots analysis tool of ArcGIS was used to calculate the spatial distribution of the high and low clusters of ESVs of the four different ecological service types (Figure 6). Higher Z-scores indicate a higher intensity of ESVs clustering. Lower Z-scores indicate a lower intensity of ESVs clustering.
From the overall spatial layout, the hot spots of provisioning services mainly aggregated in the northern part of the Yangtze River Delta urban agglomeration, and the cold spots accumulated in the southern part of the Yangtze River Delta urban agglomeration and along the coast. In the time dimension, from 1980 to 2015, provisioning services hot spots showed a trend of shrinking and aggregating, and the high-ESVs area was decreasing. The hot spots of regulating services were mainly concentrated in the southern part of the Yangtze River Delta urban agglomeration. There were a few hot spots in the western and central regions. The cold spots were mainly concentrated in the north and east of the Yangtze River Delta urban agglomeration along the coastal area. From 1980 to 2015, the hot spots of regulating services showed a trend of shrinking and aggregating, the cold spots area was spreading, and the low-ESVs area was increasing. The value of habitat services was divided into a low-ESVs north and a high-ESVs north. There are almost very large cold spots in the north. Cultural and recreational services have distributions similar to that of habitat services. From 1980 to 2015, the hot spot area of habitat services and cultural and amenity services gradually decreased, and the cold spots area gradually expanded.
From the analysis of the local spatial layout, provisioning services hot spots gathered in the administrative areas of Anhui Province, Jiangsu Province, and Shanghai, indicating that they contributed a high PSVs. Where the farmland was concentrated in these areas, the forest area was small. There were also large cold spots in Jiangsu Province, mainly in coastal cities. The area of service value in these cities decreased, and the trend of cold spots spread in Shanghai was the most obvious. According to the land use/cover change thematic map (Figure 2), Shanghai’s urban expansion has been the most obvious over time. In the past 35 years, the urban area has been encroaching on and occupying the farmland in Shanghai, and farmland degradation is the most serious in this location. Several concentrated cold spots areas in Zhejiang Province have gradually disappeared. The hot spots of regulating services were mainly concentrated in Zhejiang Province, while the proportion of forest coverage in Zhejiang Province was relatively high. The large cold spots area in Jinhua City disappeared and was replaced by hot spots. The hot spots of habitat services and cultural and amenity services were mainly concentrated in Zhejiang Province, and a small number were gathered in Anhui Province.

3.4. Analysis of the Importance of Land Use Types on the ESVs

To reflect the contribution from each land use type to the ESVs, the importance values of different land use types based on their contribution to the total ESVs were calculated. The results are listed in Table 5, and their importance values are ranked. In the study area, the water area is the most important for the total ESVs, followed by non-shrub forests, with paddy fields ranked third. Drylands, shrubs, grasslands, and wetlands are less important, and the importance of barren land is almost zero. This is because the value coefficient of the water area is much higher than that of the other land types, while the area of the wasteland is smaller. For the time points, 1980, 2000, and 2015, the sensitivity coefficient of the water area increased, the sensitivity coefficient of the paddy field decreased, and the sensitivity coefficient of the non-shrub forest first increased and then decreased. These changes in the sensitivity coefficients reflect the changes in area: The water area increased, the paddy field area decreased, and the non-shrub forest area first increased and then decreased.
In the Yangtze River Delta urban agglomeration, water, non-shrub forest, and paddy fields are important indicators when considering the total ESVs. In general, when planning land use in the region, incorporating the importance of different land use types based on their contribution to the total ESVs is necessary.

4. Discussion

4.1. Data Method Restrictions

The Yangtze River Delta urban agglomeration is located on the East coast of China. Because of the difficulties in acquiring all the needed data, the research results did not include the ocean ESVs. Further, to better unify the study area, the impact of sea reclamation on ecological services was also not included; however, these will be included in future research. The research method was based on Constance and Xie’s supplier valuation method, for which the value coefficient was improved to increase the accuracy of the results. This method has been widely used to assess the monetary costs of ecosystem services, and has proven to be a more effective method than assessing the benefits [67]. It has also been found that more extensive and more accurate field survey data is needed, if scholars do more accurately assess each type of ESVs in terms of ecosystem. This survey data not only eliminates some of the biases inherent in the overgeneralization and simplification of LULC [68], but also reduces the uncertainty in the analysis [69].

4.2. Causes of ESVs Changes in the Study Area

In the study area, changes in the ESVs were caused by changes in the type of land use. For example, in order to improve the ecological environment, the Zhejiang Provincial Government proposed in 1989 a strategic plan of “two years of preparation, five years to eliminate the barren hills, and ten years to green Zhejiang” [70]. This is the largest natural transformation activity in the history of Zhejiang Province. In 2000, the Zhejiang Provincial People’s Government announced that the strategic goal of “Greening Zhejiang for Ten Years” was achieved on schedule. Afforestation amounted to about 582,000 hectares, which is equivalent to more than 1/10 of the province’s existing forest area. All the barren hills suitable for planting forests in the province have been eliminated. The forest land area, forest stock volume, and forest coverage rate of the province all met the plan’s targets. In the same year, the Anhui Provincial Government proposed the “five eight” greening planning target of “Five-year elimination of barren hills and greening of Anhui in eight years” [71]. During the “five eight” afforestation period, a total of 324,000 ha of economic forestry was completed, accounting for 39.8% of the area of artificial afforestation, an increase of more than 20% above the quantity before the implementation. The following year, the State Council promulgated the “Outline of National Afforestation and Greening Plan from 1989 to 2000” to promote greening and afforestation throughout the country. This led to a sharp increase in the area of forestland within the study area before 2000. From 1980 to 2000, the state adjusted the industrial structure of the Yangtze River Delta region, especially by increasing the aquaculture area, which led to an increase in the water area. Before 2000, the increase of the forest and water area was the main cause of the total ESVs growth. After 2000, a series of national economic and social development, as well as urban and rural construction policies, such as the “Ninth Five-Year Plan”, resulted in the Yangtze River Delta and the Yangtze River Delta urban agglomeration ranking among the top seven economic zones, leading to faster urban residential construction and a rapidly increasing construction land. As the process of economic development and urbanization accelerates, a growing number of people are moving to cities, causing the towns to expand. The enhancement of urbanization has taken up a large amount of farmland, forest, and grassland in the Yangtze River Delta urban agglomeration, which has reduced the total ESVs.
The method used in this study, cold/hot spots analysis, can well represent the spatial difference in ecosystem services, and many scholars apply it to the spatial analysis of ecosystem services [72,73,74]. In the study area, the spatial differences between provisioning services and regulating services may have been caused by the uneven distribution of land use/land cover. The hot spots of provisioning services were in the northern part of the Yangtze River Delta because the provisioning services were mainly provided by farmland as a crop production area [75], and farmland is distributed in the north of the Yangtze River Delta urban agglomeration (Figure 1) and covers Jiangsu Province, Anhui Province, and Shanghai. The regulating services were mainly provided by forests, which are distributed in the southern part of the Yangtze River Delta urban agglomeration, mainly in Zhejiang Province. The habitat services and cultural and amenity services were found to have similar spatial distributions, which reflected the synergy between their ecological services [76]. Some scholars have found that in the middle and lower reaches of the Yangtze River, there are hot spots of species [77], and more species may live in the middle and lower reaches of the Yangtze River. The forest resources in the middle and lower reaches of the Yangtze River are abundant, indicating that the habitat quality is good, which is likely because good quality habitats are beneficial for the tourist culture. Nowadays, urbanization is becoming ever more serious. In some areas, the habitat area is reduced, biodiversity is threatened, land use changes are leading to a fragmented habitat distribution, and species are threatened with extinction [77,78,79]. In response to these environmental threats, scholars have developed appropriate decision support systems to incorporate biodiversity, ecosystem services, and land use change [19].
Some scholars have found that reductions in farmland and forest areas can lead to a decline in the ESVs. As urban land tends to be taken from high-quality farmland, forests, and grasslands, these land use changes result in ESVs losses. Additionally, it was found that the increase in forest and water areas from 1980 to 2000 resulted in an increase in the ESVs. Therefore, these changes in the ESVs were strongly related to changes in the farmland, forests, and water bodies. Further, the results of this study are consistent with those of Xin and Meng et al. [48,80]. Correlation analyses have been commonly used to explore ecosystem services trade-offs and synergies [81,82]; for example, Roces-Díaz used correlation analysis to determine the synergy between cultural and amenity services and regulating services. Therefore, the findings validated past research as they also indicated that there was a synergy between the cultural and amenity services, the regulating services, and the habitat services.

5. Conclusions and Policy Recommendations

LUCC maps were used to analyze the potential supply of four ecosystem services and evaluate for four ESVs in the Yangtze River Delta urban agglomeration. The cold/hot spots analysis tool in ArcGIS was used to calculate the spatial distributions for the high and low ESVs clusters for the four ecosystem services types. The spatial difference in four ESVs brought about by LUCC under rapid urbanization was analyzed. In most areas, and especially in the northern part of the Yangtze River Delta, there has been a significant increase in urban land use, which has led to a decline in high-quality farmland, forests, and grasslands, all of which have resulted in ESVs losses. The ESVs changes had obvious regional differences. The ESVs in Hangzhou and Jinhua City had noticeable increases; however, Shanghai had the largest ESVs declines. The spatial distributions for the four ESVs differed significantly. The provisioning services hot spots were primarily in the northern part of the Yangtze River Delta urban agglomeration, and the cold spots were primarily in the southern part. The regulating services hot spots were mainly concentrated in the southern part, there were a few hot spots in the western and central regions, and the cold spots were in the north and east along the coast. The habitat services value also show that the north is low and the south is high. The cultural and amenity services had a similar distribution to the habitat services. However, the analysis showed that over time, the hot spots were gradually shrinking and the cold spots were increasing. The Pearson’s correlation analysis showed the interrelationships between the four ecosystem services, and reflected their trade-offs and synergies. The provisioning services were negatively correlated with the habitat services and the cultural and amenity services. There was a weak positive correlation between the regulating services and the provisioning services, and a significant positive correlation with habitat services and the cultural and amenity services. Finally, the habitat services were significantly positively correlated with the cultural and amenity services. In the Yangtze River Delta urban agglomeration, the water area was the most important for the total ESVs, followed by non-bush forest. Paddy field was ranked third. Dryland, bush, grassland, and wetland were less important. The importance of barren land was almost zero. Therefore, water, non-bush forest, and paddy fields are important indicators when considering the total ESVs.
In recent years, the government has built a shelterbelt system on the banks of the Yangtze River, the Huaihe River, the Beijing-Hangzhou Grand Canal, and other lakes and has established a soil erosion control project in the small watershed area, which has improved the RSVs. Additionally, the Chinese government has formulated a policy of returning farmland to forests and returning farmland to grassland. These policies have helped ecological restoration. However, due to rapid urbanization, the current ecological environment is still deteriorating. It is still important to balance the protection of cultivated land and forest protection during the process of regional land development and utilization. The ESVs in the Yangtze River Delta urban agglomeration are declining and shows differences in space, which also reflects the disproportionate spatial pattern of land use structure. When formulating regional land use planning policies, ecological protection should be given top priority, and the ESVs should be incorporated into decision making. The importance of the contribution of different land use types to the total ESVs in the region should be considered. Therefore, within the study area, the following two suggestions are proposed based on the research results: (1) According to the research results, the contribution of the water area to the total ESVs is the most important. However, in recent years, due to rapid urbanization, agricultural pollution has increased due to, for example, the use of chemical fertilizers and pesticides and the expansion of livestock and poultry farming. An excessive amount of industrial sewage discharge exacerbates the seriousness of water pollution. Although the government has formulated a “transboundary water body joint guarantee operation” as a response to the insufficiencies in environmental protection laws, the punishment for illegal sewage disposal is not severe enough, so the challenge of protecting water from pollution still exists. Therefore, for the illegal discharge of pollutants, a high penalty system would be conducive to the protection of water. (2) The expansion of the urban scape has occupied a large amount of farmland, which has reduced the ESVs contributed by farmland. Shanghai, the largest urban area in the Yangtze River Delta, had an average annual expansion of about 3225 hectares from 1980 to 2000; however, from 2000 to 2015, its annual expansion was about 5907 hectares, and the urban annual expansion area increased by 83.16%. Therefore, it is very necessary to strictly control the construction scale of megacities and large cities by designating areas of permanent farmland as boundaries to urban development. (3) In terms of space, the western and southern parts of Zhejiang Province are rich in forest resources, and hot spots for regulating services, habitat services, and cultural and amenity services are concentrated in these areas. It is necessary to list these areas as ecological protection barriers and to define ecological corridor protection areas in the Yangtze River Basin to protect the water resources along the Yangtze River basin.

Author Contributions

Formal Analysis, Z.L.; Investigation, Z.S. and Y.T.; Methodology, Z.S. and J.Z.; Funding acquisition, W.Y. and Z.L.; Writing—original draft, Z.S.; Writing—review & editing, Z.S. and Z.L.

Funding

This work was supported by the National Natural Science Foundation of China: “Remote Sensing Dynamic Monitoring Method of Vegetation Eco-Water and Water Stress in West Sichuan Plateau” (No. 41671432) and The Center Project of Urban Development Research of Resource-Based of Sichuan (No. ZYZX-YB-1405).

Acknowledgments

The authors thank the editors and the anonymous reviewers for providing valuable comments and suggestions to improve the paper and thank MDPI and LetPub for its linguistic assistance during the preparation of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Justice, C.; Gutman, G.; Vadrevu, K.P. NASA Land Cover and Land Use Change (LCLUC): An interdisciplinary research program. J. Environ. Manag. 2015, 148, 4–9. [Google Scholar] [CrossRef] [PubMed]
  2. Estoque, R.C.; Murayama, Y. Examining the potential impact of land use/cover changes on the ecosystem services of Baguio city, the Philippines: A scenario-based analysis. Appl. Geogr. 2012, 35, 316–326. [Google Scholar] [CrossRef]
  3. Deng, X.; Li, Z.; Huang, J.; Shi, Q.; Li, Y. A revisit to the impacts of land use changes on the human wellbeing via altering the ecosystem provisioning services. Adv. Meteorol. 2013, 2013, 907367. [Google Scholar] [CrossRef]
  4. Gao, J.; Li, F.; Gao, H.; Zhou, C.; Zhang, X. The impact of land-use change on water-related ecosystem services: A study of the Guishui River Basin, Beijing, China. J. Clean. Prod. 2017, 163, S148–S155. [Google Scholar] [CrossRef]
  5. Quintas-Soriano, C.; Castro, A.J.; Castro, H.; García-Llorente, M. Impacts of land use change on ecosystem services and implications for human well-being in Spanish drylands. Land Use Policy 2016, 54, 534–548. [Google Scholar] [CrossRef]
  6. De Gregorio, A.; Jansen, L.J.M. Land Cover Classification System (LCCS) Classification Concepts and User Manual; Environment and Natural Resources Service (SDRN), Food and Agricultural Organisation of the United Nations: Rome, Italy, 2005. [Google Scholar]
  7. Knoester, L. Linking land use and ecosystem services. Development of a life cycle impact assessment method to improve evaluation of biobased products. Habitat Int. 2015, 45, 163–168. [Google Scholar]
  8. Liu, H.R.; Lin, Y.P.; Lien, W.Y. Impacts of Climate Change and Land Use Change on Ecosystem Services—A Case Study in Taiwan. Available online: http://adsabs.harvard.edu/abs/2018EGUGA..20.7449L (accessed on 21 November 2018).
  9. USGS. Land Change Science Program. 2013. Available online: http://www.usgs.gov/climate_landuse/lcs/pdfs/LCSinfosheetMarch2013.pdf (accessed on 10 August 2018).
  10. Koellner, T.; Scholz, R.W. Assessment of land use impacts on the natural environment: Part 1: An analytical framework for pure land occupation and land use change. Int. J. Life Cycle Assess. 2007, 12, 16–23. [Google Scholar]
  11. Krausmann, F.; Gingrich, S.; Eisenmenger, N.; Erb, K.H.; Haberl, H.; Fischer-Kowalski, M. Growth in global materials use, GDP and population during the 20th century. Ecol. Econ. 2009, 68, 2696–2705. [Google Scholar] [CrossRef]
  12. Millennium Ecosystem Assessment. Ecosystems and Human Well-being: Synthesis. 2005. Available online: https://www.pik-potsdam.de/ateam/avec/ma_pdf/ash.pdf (accessed on 21 November 2018).
  13. Shi, C.; Yao, S.; Chen, X. Economic Valuation of Ecosystem Services Based on Choice Experiments: A Case Study of Wenjiang in Sichuan Province. J. Nat. Resour. 2016, 31, 767–778. [Google Scholar]
  14. Braat, L.C.; Groot, R.D. The ecosystem services agenda: Bridging the worlds of natural science and economics, conservation and development, and public and private policy. Ecosyst. Serv. 2012, 1, 4–15. [Google Scholar] [CrossRef]
  15. Costanza, R.; de Groot, R.; Sutton, P.; van der Ploeg, S.; Anderson, S.J.; Kubiszewski, I.; Farber, S.; Turner, R.K. Changes in the global value of ecosystem services. Glob. Environ. Chang. 2014, 26, 152–158. [Google Scholar] [CrossRef]
  16. Song, W.; Deng, X. Land-use/land-cover change and ecosystem service provision in China. Sci. Total Environ. 2017, 576, 705–719. [Google Scholar] [CrossRef] [PubMed]
  17. Goldstein, J.H.; Caldarone, G.; Duarte, T.K.; Ennaanay, D.; Hannahs, N.; Mendoza, G.; Polasky, S.; Wolny, S.; Daily, G.C. Integrating ecosystem-service tradeoffs into land-use decisions. Proc. Natl. Acad. Sci. USA 2012, 109, 7565–7570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Mark, S.; Erica, G.; Bartuska, A.M.; Sutton-Grier, A.; Lubchenco, J. Nature as capital: Advancing and incorporating ecosystem services in United States federal policies and programs. Proc. Natl. Acad. Sci. USA 2015, 112, 7383. [Google Scholar] [CrossRef]
  19. Liang, J.; Zhong, M.; Zeng, G.; Chen, G.; Hua, S.; Li, X.; Yuan, Y.; Wu, H.; Gao, X. Risk management for optimal land use planning integrating ecosystem services values: A case study in Changsha, Middle China. Sci. Total Environ. 2017, 579, 1675–1682. [Google Scholar] [CrossRef] [PubMed]
  20. Pritchard, L.; Folke, C.; Gunderson, L. Valuation of Ecosystem Services in Institutional Context. Ecosystems 2000, 3, 36–40. [Google Scholar] [CrossRef]
  21. Costanza, R.; D’Arge, R.; Groot, R.D.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.O.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 25, 3–15. [Google Scholar] [CrossRef]
  22. Alexander, A.M.; List, J.A.; Margolis, M.; d’Argea, R.C. A method for valuing global ecosystem services. Ecol. Econ. 1998, 27, 161–170. [Google Scholar] [CrossRef]
  23. Wang, W.; Guo, H.; Chuai, X.; Dai, C.; Lai, L.; Zhang, M. The impact of land use change on the temporospatial variations of ecosystems services value in China and an optimized land use solution. Environ. Sci. Policy 2014, 44, 62–72. [Google Scholar] [CrossRef]
  24. Nahuelhual, L.; Carmona, A.; Aguayo, M.; Echeverria, C. Land use change and ecosystem services provision: A case study of recreation and ecotourism opportunities in southern Chile. Landsc. Ecol. 2014, 29, 329–344. [Google Scholar] [CrossRef]
  25. Zheng, J.K.; Yu, X.X.; Jia, G.D.; Xia, B. Dynamic evolution of the ecological service value based on LUCC in Miyun Reservoir Catchment. Trans. Chin. Soc. Agric. Eng. 2010, 26, 315–320. [Google Scholar]
  26. Fu, Y.; Du, X.; Peng, W.; Dong, F. Agro-ecosystem value connotation based on watershed land use. Trans. Chin. Soc. Agric. Eng. 2015, 31, 243–250. [Google Scholar] [CrossRef]
  27. Zhang, D.; Sun, X.; Yuan, X.; Liu, F.; Guo, H.; Xu, Y.; Li, B. Land use change and its impact on habitat quality in Lake Nansi Basin from 1980 to 2015. J. Lake Sci. 2018, 30, 349–357. [Google Scholar] [CrossRef]
  28. Ran, S.; Li, X.; Lu, C. Multi-scale Modeling of Land-cover Change and Ecosystem Service Values: A Case of the Yuzixi Catchment in Sichuan. Acta Geogr. Sin. 2006, 61, 1113–1120. [Google Scholar]
  29. Xie, G.; Li, W.; Xiao, Y.; Zhang, B.; Lu, C.; An, K.; Wang, J.; Xu, K.; Wang, J. Forest Ecosystem Services and Their Values in Beijing. Chin. Geogr. Sci. 2010, 20, 51–58. [Google Scholar] [CrossRef]
  30. Yu, X. The assessment of the forest ecosystem services valuation in China. Acta Ecol. Sin. 2005, 25, 2096–2102. [Google Scholar]
  31. Li, T.; Cui, Y.; Liu, A. Spatiotemporal dynamic analysis of forest ecosystem services using “big data”: A case study of Anhui province, central-eastern China. J. Clean. Prod. 2016, 142, 589–599. [Google Scholar] [CrossRef]
  32. Homolová, L.; Schaepman, M.E.; Lamarque, P.; Clevers, J.G.P.W.; de Bello, F.; Thuiller, W.; Lavorel, S. Comparison of remote sensing and plant trait-based modelling to predict ecosystem services in subalpine grasslands. Ecosphere 2016, 5, 1–29. [Google Scholar] [CrossRef]
  33. Sandhu, H.S.; Wratten, S.D.; Cullen, R.; Case, B. The future of farming: The value of ecosystem services in conventional and organic arable land. An experimental approach. Ecol. Econ. 2008, 64, 835–848. [Google Scholar] [CrossRef] [Green Version]
  34. Kroeger, T.; Casey, F. An assessment of market-based approaches to providing ecosystem services on agricultural lands. Ecol. Econ. 2007, 64, 321–332. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, Z.; Zhang, B.; Zhang, S.; Li, X.; Liu, D.; Song, K.; Li, J.; Li, F.; Duan, H. Changes of Land Use and of Ecosystem Service Values in Sanjiang Plain, Northeast China. Environ. Monit. Assess. 2006, 112, 69–91. [Google Scholar] [CrossRef] [PubMed]
  36. Gunderson, L.H.; Cosens, B.; Garmestani, A.S. Adaptive governance of riverine and wetland ecosystem goods and services. J. Environ. Manag. 2016, 183 Pt 2, 353–360. [Google Scholar] [CrossRef]
  37. Barbier, E.B.; Hacker, S.D.; Kennedy, C.; Koch, E.W.; Stier, A.C.; Silliman, B.R. The value of estuarine and coastal ecosystem services. Ecol. Monogr. 2011, 81, 169–193. [Google Scholar] [CrossRef] [Green Version]
  38. De Groot, R.; Brander, L.; van der Ploeg, S.; Costanza, R.; Bernard, F.; Braat, L.; Christie, M.; Crossman, N.; Ghermandi, A.; Hein, L.; et al. Global estimates of the value of ecosystems and their services in monetary units. Ecosyst. Serv. 2012, 1, 50–61. [Google Scholar] [CrossRef] [Green Version]
  39. Rawlins, M.A.; Westby, L. Community participation in payment for ecosystem services design and implementation: An example from Trinidad. Ecosyst. Serv. 2013, 6, 117–121. [Google Scholar] [CrossRef]
  40. 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]
  41. Brown, G. The relationship between social values for ecosystem services and global land cover: An empirical analysis. Ecosyst. Serv. 2013, 5, 58–68. [Google Scholar] [CrossRef]
  42. Jobstvogt, N.; Watson, V.; Kenter, J.O. Looking below the surface: The cultural ecosystem service values of UK marine protected areas (MPAs). Ecosyst. Serv. 2014, 10, 97–110. [Google Scholar] [CrossRef] [Green Version]
  43. Brander, L.M.; Wagtendonk, A.J.; Hussain, S.S.; McVittie, A.; Verburg, P.H.; de Groot, R.S.; van der Ploeg, S. Ecosystem service values for mangroves in Southeast Asia: A meta-analysis and value transfer application. Ecosyst. Serv. 2012, 1, 62–69. [Google Scholar] [CrossRef] [Green Version]
  44. Bullock, C.; Joyce, D.; Collier, M. An exploration of the relationships between cultural ecosystem services, socio-cultural values and well-being. Ecosyst. Serv. 2018, 31, 142–152. [Google Scholar] [CrossRef]
  45. Li, J.C.; Wang, W.L.; Hu, G.Y.; Wei, Z. Changes in ecosystem service values in Zoige Plateau, China. Agric. Ecosyst. Environ. 2010, 139, 766–770. [Google Scholar] [CrossRef]
  46. Fei, L.; Zhang, S.; Yang, J.; Chang, L.; Yang, H.; Bu, K. Effects of land use change on ecosystem services value in West Jilin since the reform and opening of China. Ecosyst. Serv. 2018, 31, 12–20. [Google Scholar] [CrossRef]
  47. Wang, Y.; Dai, E.; Yin, L.; Ma, L. Land use/land cover change and the effects on ecosystem services in the Hengduan Mountain Region, China. Ecosyst. Serv. 2018, 34, 55–67. [Google Scholar] [CrossRef]
  48. Cao, L.; Li, J.; Ye, M.; Pu, R.; Liu, Y.; Guo, Q.; Feng, B.; Song, X. Changes of Ecosystem Service Value in a Coastal Zone of Zhejiang Province, China, during Rapid Urbanization. Int. J. Environ. Res. Public Health 2018, 15, 1301. [Google Scholar] [CrossRef] [PubMed]
  49. Arowolo, A.O.; Deng, X.; Olatunji, O.A.; Obayelu, A.E. Assessing changes in the value of ecosystem services in response to land-use/land-cover dynamics in Nigeria. Sci. Total Environ. 2018, 636, 597–609. [Google Scholar] [CrossRef] [PubMed]
  50. Yu, Z.; Qin, T.; Yan, D.; Yang, M.; Yu, H.; Shi, W. The Impact on the Ecosystem Services Value of the Ecological Shelter Zone Reconstruction in the Upper Reaches Basin of the Yangtze River in China. Int. J. Environ. Res. Public Health 2018, 15, 2273. [Google Scholar] [CrossRef] [PubMed]
  51. Koellner, T.; de Baan, L.; Beck, T.; Brandão, M.; Civit, B.; Margni, M.; Canals, L.M.; Saad, R.; de Souza, D.M.; et al. UNEP/SETAC guideline on global land use impact assessment on biodiversity and ecosystem services in LCA. Int. J. Life Cycle Assess. 2013, 18, 1188–1202. [Google Scholar] [CrossRef]
  52. 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. Landsc. Urban Plan. 2013, 116, 60–72. [Google Scholar] [CrossRef]
  53. Rong, Y.; Li, C.; Xu, C.; Yan, Y. Ecosystem service values and spatial differentiation changes during urbanization: A case study of Huanghua City. Chin. J. Ecol. 2017, 36, 1374–1381. [Google Scholar]
  54. Xie, G.; Zhang, C.; Zhen, L.; Zhang, L. Dynamic changes in the value of China’s ecosystem services. Ecosyst. Serv. 2017, 26, 146–154. [Google Scholar] [CrossRef]
  55. Du, H.; Wang, D.; Wang, Y.; Zhao, X.; Qin, F.; Jiang, H.; Cai, Y. Influences of land cover types, meteorological conditions, anthropogenic heat and urban area on surface urban heat island in the Yangtze River Delta Urban Agglomeration. Sci. Total Environ. 2016, 571, 461–470. [Google Scholar] [CrossRef] [PubMed]
  56. Xu, X.; Liu, J.; Zhang, S.; Li, R.; Yan, C.; Wu, S. China’s Multi-Period Land Use Land Cover Remote Sensing Monitoring Data Set (CNLUCC); Resource and Environment Data Cloud Platform: Beijing, China, 2018. [Google Scholar]
  57. Kumar, P. (Ed.) The Economics of Ecosystems and Biodiversity (TEEB): Ecological and Economic Foundations; Earthscan: London, UK, 2010. [Google Scholar]
  58. Aschonitis, V.G.; Gaglio, M.; Castaldelli, G.; Fano, E.A. Criticism on elasticity-sensitivity coefficient for assessing the robustness and sensitivity of ecosystem services values. Ecosyst. Serv. 2016, 20, 66–68. [Google Scholar] [CrossRef] [Green Version]
  59. Xiong, Y.; Zhang, F.; Gong, C.; Luo, P. Spatial-temporal Evolvement of Ecosystem Service Value in Hunan Province Based on LUCC. Resour. Environ. Yangtze Basin 2018, 27, 1397–1408. [Google Scholar]
  60. Yu, H. The influential factors of China’s regional energy intensity and its spatial linkages: 1988–2007. Energy Policy 2012, 45, 583–593. [Google Scholar] [CrossRef]
  61. Gimona, A.; van der Horst, D. Mapping hotspots of multiple landscape functions: A case study on farmland afforestation in Scotland. Landsc. Ecol. 2007, 22, 1255–1264. [Google Scholar] [CrossRef]
  62. Roces-Díaz, J.V.; Burkhard, B.; Kruse, M.; Müller, F.; Díaz-Varela, E.R.; Álvarez-Álvarez, P. Use of ecosystem information derived from forest thematic maps for spatial analysis of ecosystem services in northwestern Spain. Landsc. Ecol. Eng. 2017, 13, 45–57. [Google Scholar] [CrossRef]
  63. Ferreira, L.M.R.; Esteves, L.S.; Souza, E.P.D.; Costa dos Santos, C.A. Impact of the Urbanisation Process in the Availability of Ecosystem Services in a Tropical Ecotone Area. Ecosystems 2018, 1–17. [Google Scholar] [CrossRef] [Green Version]
  64. Getis, A.; Ord, J.K. The Analysis of Spatial Association by Use of Distance Statistics. Geogr. Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  65. Timilsina, N.; Escobedo, F.J.; Cropper, W.P.; Abd-Elrahman, A.; Brandeis, T.J.; Delphin, S.; Lambert, S. A framework for identifying carbon hotspots and forest management drivers. J. Environ. Manag. 2013, 114, 293–302. [Google Scholar] [CrossRef]
  66. Bagstad, K.J.; Reed, J.M.; Semmens, D.J.; Sherrouse, B.S.; Troy, A. Linking biophysical models and public preferences for ecosystem service assessments: A case study for the Southern Rocky Mountains. Reg. Environ. Chang. 2015, 16, 1–14. [Google Scholar] [CrossRef]
  67. Zhan, J.; Chu, X.; Li, Z.; Jia, S.; Wang, G. Incorporating ecosystem services into agricultural management based on land use/cover change in Northeastern China. Technol. Forecast. Soc. Chang. 2018. [Google Scholar] [CrossRef]
  68. Eigenbrod, F.; Armsworth, P.R.; Anderson, B.J.; Heinemeyer, A.; Gillings, S.; Roy, D.B.; Thomas, C.D.; Gaston, K.J. The impact of proxy-based methods on mapping the distribution of ecosystem services. J. Appl. Ecol. 2010, 47, 377–385. [Google Scholar] [CrossRef] [Green Version]
  69. Hou, Y.; Burkhard, B.; Müller, F. Uncertainties in landscape analysis and ecosystem service assessment. J. Environ. Manag. 2013, 127 (Suppl. 3), 117–131. [Google Scholar] [CrossRef] [PubMed]
  70. Zhejiang News. Available online: https://zj.zjol.com.cn/news/492259.html (accessed on 8 December 2018).
  71. State Forestry Administration. Available online: http://www.forestry.gov.cn/portal/main/s/1017/content-266121.html (accessed on 8 December 2018).
  72. Hu, X.; Hong, W.; Qiu, R.; Hong, T.; Chen, C.; Wu, C. Geographic variations of ecosystem service intensity in Fuzhou City, China. Sci. Total Environ. 2015, 512–513, 215–226. [Google Scholar] [CrossRef]
  73. Dvarskas, A. Mapping ecosystem services supply chains for coastal Long Island communities: Implications for resilience planning. Ecosyst. Serv. 2018, 30, 14–26. [Google Scholar] [CrossRef]
  74. 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]
  75. Lin, S.; Wu, R.; Yang, F.; Wang, J.; Wu, W. Spatial trade-offs and synergies among ecosystem services within a global biodiversity hotspot. Ecol. Indic. 2018, 84, 371–381. [Google Scholar] [CrossRef]
  76. Kong, L.; Zheng, H.; Xiao, Y.; Ouyang, Z.; Li, C.; Zhang, J.; Huang, B. Mapping Ecosystem Service Bundles to Detect Distinct Types of Multifunctionality within the Diverse Landscape of the Yangtze River Basin, China. Sustainability 2018, 10, 857. [Google Scholar] [CrossRef]
  77. Liang, L.; Xing, W.; Zeng, G.; Li, X.; Peng, Y.; Li, X.; Gao, X.; He, X. Where will threatened migratory birds go under climate change? Implications for China’s national nature reserves. Sci. Total Environ. 2018, 645, 1040–1047. [Google Scholar] [CrossRef]
  78. Liang, J.; He, X.; Zeng, G.; Zhong, M.; Gao, X.; Li, X.; Li, X.D.; Wu, H.; Feng, C.; Xing, W.; et al. Integrating priority areas and ecological corridors into national network for conservation planning in China. Sci. Total Environ. 2018, 626, 22–29. [Google Scholar] [CrossRef]
  79. Shackelford, N.; Standish, R.J.; Ripple, W.; Starzomski, B.M. Threats to biodiversity from cumulative human impacts in one of North America’s last wildlife frontiers. Conserv. Biol. 2018, 32, 672–684. [Google Scholar] [CrossRef] [PubMed]
  80. Xin, Z.; Li, C.; Liu, H.; Shang, H.; Ye, L.; Li, Y.; Zhang, C. Evaluation of Temporal and Spatial Ecosystem Services in Dalian, China: Implications for Urban Planning. Sustainability 2018, 10, 1247. [Google Scholar] [CrossRef]
  81. Benis, E.; Belinda, R.; Mathieu, R.; Bode, M.; Richardson, D.M. Spatial congruence between biodiversity and ecosystem services in South Africa. Biol. Conserv. 2009, 142, 553–562. [Google Scholar] [CrossRef]
  82. Mouchet, M.A.; Lamarque, P.; Martín-López, B.; Crouzat, E.; Gos, P.; Byczek, C.; Lavorel, S. An interdisciplinary methodological guide for quantifying associations between ecosystem services. Glob. Environ. Chang. 2014, 28, 298–308. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Ijerph 16 00123 g001
Figure 2. Spatial distribution of land use and land cover from 1980 to 2015.
Figure 2. Spatial distribution of land use and land cover from 1980 to 2015.
Ijerph 16 00123 g002
Figure 3. Time trends of different land use/cover types.
Figure 3. Time trends of different land use/cover types.
Ijerph 16 00123 g003
Figure 4. Time-varying trends of ESVs in the Yangtze River Delta urban agglomeration. (a) Regulating services values (RSVs); (b) provisioning services values (PSVs); (c) habitat services value (HSV); (d) cultural and amenity services value (CSV); and (e) total ecosystem services values (ESVs).
Figure 4. Time-varying trends of ESVs in the Yangtze River Delta urban agglomeration. (a) Regulating services values (RSVs); (b) provisioning services values (PSVs); (c) habitat services value (HSV); (d) cultural and amenity services value (CSV); and (e) total ecosystem services values (ESVs).
Ijerph 16 00123 g004
Figure 5. Changes in ESVs in different urban areas of the Yangtze River Delta urban agglomeration.
Figure 5. Changes in ESVs in different urban areas of the Yangtze River Delta urban agglomeration.
Ijerph 16 00123 g005
Figure 6. Potential ecosystem services supply value based on the G i * index Z score.
Figure 6. Potential ecosystem services supply value based on the G i * index Z score.
Ijerph 16 00123 g006
Table 1. The ecosystem services values (ESVs) per unit area of terrestrial ecosystems, yuan/ha.
Table 1. The ecosystem services values (ESVs) per unit area of terrestrial ecosystems, yuan/ha.
FarmlandForestGrasslandWater AreaWetlandBarren Land
Paddy FieldDry LandBushNon-BushGrasslandWater AreaWetlandBarren Land
Provisioning servicesFood3427.9342142.459478.9026781.3674554.51882016.4321285.4750
Material226.84861008.2161083.8321789.583831.7782579.72421260.270
Water2974.23750.4108554.5188932.5998453.697220,895.286528.1990
Regulating servicesAir quality regulation2797.7991688.76183553.9615923.2692873.4161940.8164789.02650.4108
Climate regulation1436.708907.394410,661.8817,719.47612.0315772.0379073.9440
Wastewater treatment428.4918252.0543226.2915015.8752520.5413,9899073.944252.054
Regulation of water flow6855.869680.54588443.8098847.0955570.393257,70061,072.6875.6162
Erosion protection25.20542596.15624335.3297208.7443503.5512344.1025822.44750.4108
Soil fertility maintenance478.9026302.4648327.6702554.5188277.2594176.4378453.69720
Habitat servicesBiodiversity529.3134327.67023957.2486553.4043201.0866427.37719,836.6550.4108
Cultural and amenity servicesAesthetic landscape226.8486151.23241739.1732873.4161411.5024763.82111,922.1525.2054
Table 2. Changes in land use/cover area and ESVs in the Yangtze River Delta urban agglomeration from 1980 to 2015.
Table 2. Changes in land use/cover area and ESVs in the Yangtze River Delta urban agglomeration from 1980 to 2015.
Land Use/Cover Type1980–2000 Change2000–2015 Change
Area (10,000 ha)%Value (100 Million Yuan)%Area (10,000 ha)%Value (100 Million Yuan)%
Paddy field−33.71−3.565−65.425−3.565−69.96−7.673−135.779−7.673
Dry land−8.85−4.178−8.945−4.178−6.7−3.301−6.772−3.301
Bush3.736.93414.3096.934−0.22−0.382−0.844−0.382
Non Bush22.434.524130.5054.523−5.37−1.036−31.260−1.036
Grassland−20.1−20.643−55.359−19.915−0.71−0.919−2.045−0.918
Water5.533.984175.0833.9836.144.254194.3964.253
Wetland−5.53−20.044−72.509−20.043−3.5−15.866−45.892−15.865
Barren land−0.1−22.727−0.005−22.7270.1852.9410.00952.941
Urban36.629.508990080.1449.8910
Total value change117.6541.156Total value change−28.188−0.274
total area2095.06total area2095.06
Table 3. ESVs in the Yangtze River Delta urban agglomeration Moran I spatial autocorrelation index.
Table 3. ESVs in the Yangtze River Delta urban agglomeration Moran I spatial autocorrelation index.
Spatial Autocorrelation of ESVs
Provisioning ServicesRegulating ServicesHabitat ServicesCultural and Amenity Services
Index IZ ScoreIndex IZ ScoreIndex IZ ScoreIndex IZ Score
19800.32585915.768063510.34131513.544840100.44389921.479966420.43389121.3177653
20000.29666314.354844760.34493516.690576570.47195322.836650610.45632521.7954661
20150.26461512.807529280.34594216.7435760.47847122.974222870.47462322.9187627
Table 4. Correlation coefficients of four potential ecosystem services.
Table 4. Correlation coefficients of four potential ecosystem services.
Ecosystem ServicesRegulating ServicesProvisioning ServicesHabitat ServicesCultural and Amenity Services
Regulating services1---
Provisioning services0.2191--
Habitat services0.857−0.1351-
Cultural and amenity services0.871−0.1500.9981
Table 5. Sensitivity coefficient of ESVs in the Yangtze River Delta urban agglomeration.
Table 5. Sensitivity coefficient of ESVs in the Yangtze River Delta urban agglomeration.
198020002015
Sensitivity CoefficientRankSensitivity CoefficientRankSensitivity CoefficientRank
Paddy field0.178530.172130.15893
Dry land0.020860.019970.01937
Bush0.020170.021560.02146
Non-bush0.280120.292720.28972
Grassland0.027350.021750.02155
Water area0.427510.444510.46341
Barren land080808
Wetland0.035240.028140.02374

Share and Cite

MDPI and ACS Style

Li, Z.; Sun, Z.; Tian, Y.; Zhong, J.; Yang, W. Impact of Land Use/Cover Change on Yangtze River Delta Urban Agglomeration Ecosystem Services Value: Temporal-Spatial Patterns and Cold/Hot Spots Ecosystem Services Value Change Brought by Urbanization. Int. J. Environ. Res. Public Health 2019, 16, 123. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010123

AMA Style

Li Z, Sun Z, Tian Y, Zhong J, Yang W. Impact of Land Use/Cover Change on Yangtze River Delta Urban Agglomeration Ecosystem Services Value: Temporal-Spatial Patterns and Cold/Hot Spots Ecosystem Services Value Change Brought by Urbanization. International Journal of Environmental Research and Public Health. 2019; 16(1):123. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010123

Chicago/Turabian Style

Li, Zhigang, Zishu Sun, Yangjie Tian, Jialong Zhong, and Wunian Yang. 2019. "Impact of Land Use/Cover Change on Yangtze River Delta Urban Agglomeration Ecosystem Services Value: Temporal-Spatial Patterns and Cold/Hot Spots Ecosystem Services Value Change Brought by Urbanization" International Journal of Environmental Research and Public Health 16, no. 1: 123. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16010123

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