Next Article in Journal
Investigating the Relationship between Construction Supply Chain Integration and Sustainable Use of Material: Evidence from China
Next Article in Special Issue
Postseismic Restoration of the Ecological Environment in the Wenchuan Region Using Satellite Data
Previous Article in Journal
A Comparative Analysis of the Visual Comfort Performance between a PCM Glazing and a Conventional Selective Double Glazed Unit
Previous Article in Special Issue
Does the Exhaustion of Resources Drive Land Use Changes? Evidence from the Influence of Coal Resources-Exhaustion on Coal Resources–Based Industry Land Use Changes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Land-Use/Land-Cover Changes and Their Influence on the Ecosystem in Chengdu City, China during the Period of 1992–2018

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100094, China
3
School of Life and Environmental Sciences, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA
4
School of Tourism and Geography Science, Qingdao University, Qingdao 266071, Shandong, China
*
Authors to whom correspondence should be addressed.
These authors contributed to this work equally.
Sustainability 2018, 10(10), 3580; https://0-doi-org.brum.beds.ac.uk/10.3390/su10103580
Submission received: 12 September 2018 / Revised: 30 September 2018 / Accepted: 2 October 2018 / Published: 8 October 2018

Abstract

:
Due to urban expansion, economic development, and rapid population growth, land use/land cover (LULC) is changing in major cities around the globe. Quantitative analysis of LULC change is important for studying the corresponding impact on the ecosystem service value (ESV) that helps in decision-making and ecosystem conservation. Based on LULC data retrieved from remote-sensing interpretation, we computed the changes of ESV associated with the LULC dynamics using the benefits transfer method and geographic information system (GIS) technologies during the period of 1992–2018 following self-modified coefficients which were corrected by net primary productivity (NPP). This improved approach aimed to establish a regional value coefficients table for facilitating the reliable evaluation of ESV. The main objective of this research was to clarify the trend and spatial patterns of LULC changes and their influence on ecosystem service values and functions. Our results show a continuous reduction in total ESV from United States (US) $1476.25 million in 1992, to US $1410.17, $1335.10, and $1190.56 million in 2001, 2009, and 2018, respectively; such changes are attributed to a notable loss of farmland and forest land from 1992–2018. The elasticity of ESV in response to changes in LULC shows that 1% of land transition may have caused average changes of 0.28%, 0.34%, and 0.50% during the periods of 1992–2001, 2001–2009, and 2009–2018, respectively. This study provides important information useful for land resource management and for developing strategies to address the reduction of ESV.

1. Introduction

Ecosystems are important for living on Earth and also for survival of the human population [1,2,3] and livelihood [4,5,6,7,8]. A good ecosystem covers food, water, and other raw materials, as well as environmental (e.g., hydrological system and climate) and cultural services (e.g., recreation and culture values), which have a direct impact on population, while also supporting services (e.g., pollination and soil formation) that have a relatively indirect impact [4,6,9,10]. Economic assessment of these services can quantify the benefits derived from the ecosystem [11,12,13], an important tool for raising citizen awareness, helping develop management knowledge of natural capital, improving decision-making for limited resources in competitive demand, and providing incentives to protect the ecosystems [14,15,16,17,18,19].
The evaluation method of ecosystem service value (ESV) can be divided into four categories: (1) cost-based approaches, such as avoided cost and replacement cost; (2) revealed preference approaches, such as market prices and travel cost; (3) stated preference approaches, such as contingent valuation and choice experiments; and (4) the benefits transfer method (BTM) [10]. Most notably, in 1997, Costanza et al. used BTM to estimate global ESV with an average of United States (US) $33 trillion (in 1995) per year [4]. Since then, many academics utilized BTM for ESV evaluation to examine various natural resources including farmland [20,21], forest [22], grassland [23], and wetland [24,25], and although criticized in some cases for uncertainty, this ESV assessment method proposed by Xie et al. [4] is still widely used because of its feasibility, especially in areas where data are scarce [26,27].
The ecosystem is directly affected by changes in land use/land cover (LULC). However, due to the development of society and the rapid increase in population, the speed, degree, and intensity of LULC changes are now faster compared to the past, and a large number of landscapes on Earth are getting disturbed [28]. For instance, in the tropics, more than 55% of new agricultural land was at the expense of intact forests, while 28% was associated with disturbed forests from 1980–2000 [29]. Changes in LULC influence ecosystem services by increasing the availability of certain services while reducing other services that influence the ability of the biosystem to support human needs, further impacting ecological degradation [21]. Studies showed that 1% of land conversion led to an average change in ESV of 0.10% during the period of 2000–2008 in China [30].
Efforts are being made to understand, model, evaluate, and manage ecosystem services and natural capital [31]. For example, based on the proxy-based method and geographically weighted regression (GWR), Su et al. described the ESV changes in Shanghai’s surrounding areas (China) and their relationship with urbanization, and the results showed that the total ESV dropped from 4718.1 to 3263.6 million yuan from 1994–2006 [32]. Taking the Munessa Shashemene landscape on the Ethiopian plateau as an example, Kindu et al. used revised conservative value coefficients to estimate the response of ESV changes to LULC dynamics during the period of 1973–2012 [33]. The decline in ESV reflected the impact of ecological degradation in the study area, and further research was recommended to explore future options and develop intervention strategies. Xie et al. used the equivalent value factor table of the Chinese terrestrial ecosystem to study the response of the ecosystem to LULC changes at the national, basin, and regional levels [34].
As an important commodity as a grain production base, the contradiction between food security and ecological civilization in Chengdu Plain is becoming more prominent due to the rapid growth of cities and the associated environment [35]. A study of LULC changes and their influence on the value of ecosystem services in Chengdu has important practical implications for ensuring food security while protecting ecological civilization for the evaluation of ecosystem services in southwest China [36]. However, to date, efforts in the quantitative analysis of the impact of LULC changes on ESV and the spatiotemporal variability and regional differences of ESV changes in Chengdu are very limited. This is the first study of its kind in Chengdu (the core city in southwest China) that contributes to the qualitative and quantitative analysis of LULC changes and their influence on the ecosystem based on regional value coefficients. The main aims of the present study were (1) to unfold LULC changes and spatial patterns in Chengdu for the period of 1992–2018; (2) to assess total ESV response to LULC changes and the spatial contribution of individual ecosystem service functions changes during this period based on modified coefficients; (3) to analyze ESV changes in different ecosystems through locally modified coefficients, compared with two kinds of global coefficients proposed by Costanza et al. in 1997 and 2014 [4,31]; and (4) to determine the elasticity of ecosystem service changes with respect to LULC changes.

2. Materials and Methods

2.1. Description of the Study Area

Chengdu, located in the central part of Sichuan Province, is the capital of the Sichuan Province and one of the sub-provincial cities in China. The province is the science and technology center, trade center, financial center, and transportation and communication hub of southwest China [37]. Geographically, Chengdu is located between 30°05′ and 31°26′ north and 102°54′ and 104°53′ east (Figure 1). The land area covers 14,335 km2, accounting for 2.95% of the province’s total area, of which the built-up area is 837.27 km2. The resident population of Chengdu reached 15.92 million with an urbanization rate of 70.62%, and the gross domestic product (GDP) of Chengdu ranked first among China’s sub-provincial cities in 2016, at US $183.22 billion [38]. Chengdu has 11 districts, four counties, and four county-level cities.
The climate of the study area is subtropical monsoon, with an average annual temperature of about 18 °C and an average annual precipitation of more than 1000 mm [39]. Natural resource conditions make Chengdu an important agriculturally producing area [40]. As a significant ecological buffer zone, the ecological environment of Chengdu has an important impact on the ecological security of the Yangtze River Basin and the Three Gorges reservoir area [36]. However, the rapid increase in related resource consumption resulted in resource inefficiency, which led to an increasingly prominent contradiction between Chengdu’s economic development and environmental protection. Assessing land-use changes and their influence on ESV is important for improving resource utilization and promoting sustainable development.

2.2. Data Collection and LULC Dataset

Landsat-5 Thematic Mapper (1992, 2001, and 2009) and Landsat-8 Operation Land Imager/Thermal Infrared Sensor (2018) satellite images were used for the LULC change analysis with 30-m spatial resolution, which were taken from the United States Geological Survey (USGS; https://glovis.usgs.gov/). The administrative division vector border data of Chengdu were obtained from the National Science and Technology Infrastructure of China (http://nnu.geodata.cn:8008/). Additionally, some characteristic information about LULC type for remote-sensing interpretation was obtained through field survey.
Satellite images were preprocessed for geometric correction, radiometric calibration, atmospheric correction, mosaicking, and cropping before classification [41,42,43]. The Fast Line-of-Sight Atmospheric Analysis of Hypercubes (FLAASH) module was applied to eliminate atmospheric effects using the ENVI 5.3 Software [44]. For extracting useful LULC information, the maximum-likelihood classification (MLC) algorithm, a supervised classification method, was applied to land-use classification, which is considered to be a simple and powerful method when using sufficiently accurate samples in training [45,46]. Chengdu was divided into six classes: farmland (paddy fields and dry land), forest land (dense trees, shrubs, sparse woodlands, and other wooded land), grassland (various types of grasses with coverage above 5%), water area (river canal, lake, pond, and reservoir), construction land (urban and rural residential areas, transportation land, industrial and mining land, and water conservancy facilities), and unused land. Training samples were collected via the method of on-screen selection of accurate polygons based on high-resolution Google Earth images and existing LULC maps from the Chinese Academy of Sciences Resource and Environmental Data Cloud Platform (http://www.resdc.cn/), which was previously field checked [47]. The uniformity of training samples and their representation throughout the entire image were evaluated by Jeffries-Matusita (JM) distance. The training samples were evenly distributed in the study area. A total of 450 training samples (75 polygons in each LULC) were collected throughout the whole image in each period. In this study, the JM distance for each period was greater than 1.9, indicating good separability of training samples [48]. Furthermore, there was no probability threshold so that each pixel was classified to the category with highest likelihood [48]. After classifying, the LULC datasets for the four periods were obtained.
The accuracy of LULC classification was evaluated based on ground-truth samples which were collected from reference data including Google Earth satellite images, field survey data, and existing LULC maps. Four groups of random sampling points were created using the ArcGIS v.10.1 software for four periods, and each group consisted of 450 points. A similar approach was taken by Gashaw et al. [43] and Tilahun et al. [49]. We got total accuracies of 86.33, 87.33, 87.00, and 87.66% for the 1992, 2001, 2009, and 2018 classification results, respectively, while corresponding kappa coefficients were 0.76, 0.76, 0.79, and 0.83. Therefore, the LULC dataset used for the present changes in LULC in Chengdu was reliable [50].

2.3. LULC Changes

LULC changes have an impact on the ecosystem and its function [51]. Here, we carried out quantitative analysis of LULC changes of Chengdu City using two indices to reveal regional differences and changing trends [51,52]. The two indices are defined as follows:
  K = A f i n a l A i n i t i a l A i n i t i a l   ×   1 T   ×   100 % ,  
where K refers to the land-use dynamic index for a single land-use category, Afinal and Ainitial are the areas of a certain land use at the final and initial years of a period, respectively, and T is the study period. If T is one year, K refers to the annual change rate. The land-use dynamic degree of the study period (T) is defined as
  S = ( i , j n ( Δ A i j A i ) )   ×   1 T   ×   100 % ,  
where Ai is the area of land-use category i at the initial year of the period, ΔAi→j is the total area of land-use category i converted into land-use category j, and n represents the types of land use (farmland, forest land, grassland, water area, construction land, and unused land).

2.4. Assignment of Ecosystem Services Values

Costanza et al. [4] estimated the global economic value of 17 ecosystem services for 16 biomes using BTM, based on existing studies and original calculations. Since then, studies on global and regional ecosystem service value assessments developed rapidly [20,33,53]. However, owing to the low estimate of the unit value of farmland and uncertainties, this work received some critiques [34]. Costanza et al. provided an updated unit ecosystem service value to estimate the global value of ecosystem services in 2011, which was US $125 trillion/year in 2007 (Table 1), which caused large deviations in China [31]. As a result, on the basis of a questionnaire survey of 700 ecological experts, Xie et al. developed a table of ecosystem service value equivalent factors for China, which divided ecosystem services into four types and nine sub-types (Table 1) [54]. This method assumed that the equivalent value per unit area of food production of farmland is 1; thus, the equivalent value of other ecosystems can be quantified by comparing their utility to food production on farm land [30]. The table of ecosystem service value equivalent factors was considered as a modification of the method proposed by Costanza et al. (2014) [10].
In order to revise unit ecosystem service values from a national scale to a regional scale, this study employed ecosystem service value developed for China by Xie et al. in 2008 [54], but proposed a correction coefficient based on the net primary productivity (NPP) for Chengdu. The correction coefficient formula is as follows:
  C C i = N P P l i N P P i ,  
where CCi denotes the correction coefficient of a certain i-th land ecosystem, NPPi is the NPP of a certain i-th land ecosystem in China, and NPPli is the NPP of a certain i-th land ecosystem in Chengdu [55,56,57,58,59]. The ecosystem service value coefficients of each land ecosystem in Chengdu were obtained (Table 1).

2.5. Calculation of ESV

Based on ecosystem service value coefficients in Table 1, this study calculated the total ESV in Chengdu and each district, as employed in other researches [24,60,61,62], as follows:
  E S V t = i = 1 n ( A i t   ×   V C i ) ,  
where ESVt is the assessed total ESV at time t, Ait is the area (hm2) of land-use type i at time t, and VCi is the ecosystem service value coefficient of land-use type i in Table 1.
The change rate of ESV through the study period was estimated using the following formula:
  E S V c r = ( E S V f i n a l   y e a r E S V i n i t i a l   y e a r E S V i n i t i a l   y e a r )   ×   100 % ,  
where ESVcr is the change rate of the initial year to the final year, and ESVinitial year and ESVfinal year refer to the total ESV at the beginning and end of the year, respectively.
The value of individual ecosystem service functions was calculated using the following expression:
  E S V f t = i = 1 n ( A i t   ×   V C f i ) ,  
where ESVft is the assessed ESV of individual function f at time t, Ait is the area (hm2) of land-use type i at time t, and VCfi is the ecosystem service value coefficient of individual function f for land-use type i.

2.6. Elasticity for the Response of ESV to LULC Changes

In economics, elasticity is an indicator of how well an economic variable responds to another economic variable [30]. In order to discuss the response of ESV to LULC changes, this paper used an elasticity indicator, which can be applied to measure the percentage change of ESV due to the percentage change of LULC.
  E E L = | ( E S V f i n a l   y e a r E S V i n i t i a l   y e a r E S V i n i t i a l   y e a r )   ×   1 T   ×   100 % L T P | ;  
  L T P =   i = 1 n Δ L C A i i = 1 n L C A i   ×   1 T   ×   100 % .  
In the above expressions, EEL is the elasticity for the response of ESV to LULC changes, ESVfinal year and ESVinitial year have the same meanings as in Equation (5), LTP is the land transition percentage, which can display the transition speed and degree of LULC changes, ΔLCAi is the converted area of LULC category i, LCAi is the area of LULC category i, and T represents the research period.

3. Results

3.1. LULC Patterns in Chengdu

3.1.1. LULC of Chengdu from 1992–2018

Based on remote-sensing images, four maps of LULC patterns of Chengdu during the period of 1992–2018 were obtained (Figure 2). In terms of spatial distribution of LULC in Chengdu, Figure 2 shows that the farmland was mainly distributed in flat areas in the east and central regions of Chengdu, while the construction land was mainly distributed in the downtown area of Chengdu and the centers of various districts and counties. Moreover, the forest land and grassland were mainly distributed in the southwest and northwest regions of the study area. Farmland continuously represented the most dominant LULC type, with proportions of 59.56% in 1992, 62.00% in 2001, 58.71% in 2009, and 51.86% in 2018 relative to the total area, while unused land comprised the smallest LULC proportion (Table 2).
LULC in Chengdu changed dramatically during the period 1992 to 2018, which was characterized by a decline in farmland, forest land, and water area, and an increase in grassland, construction land, and unused land (Figure 2; Table 2). The farmland decreased from 59.56% (717,033.15 hm2) in 1992 to 51.86% (624,366.81 hm2) in 2018, forest land from 34.40% (414,105.48 hm2) to 26.47% (318,693.69 hm2), and water area from 1.37% (16,509.15 hm2) to 1.01% (12,126.96 hm2). On the contrary, grassland increased from 1.85% (22,228.11 hm2) in 1992 to 2.5% (30,102.30 hm2) in 2018, and construction land from 2.83% (34,057.17 hm2) to 18.16% (218,631.96 hm2), while the area of unused land was extremely limited with an increase of 31.95 ha. Moreover, in order quantify LULC changes, a change matrix of LULC conversion in Chengdu for the period of 1992–2018 was obtained (Table 3). The percentage of reservation of unchanged forest land was 68.26%, of which about 23.62% and 4.70% of the forest land was converted to farmland and grassland, respectively, while the spread of construction land was mainly due to the reduction in farmland and water area, accounting for 24.10% and 23.94%, respectively (Table 3).
The land-use dynamic index K of a single land-use category (Table 2) and the land-use dynamic degree S of study period T was calculated using Equations (1) and (2). The K of the forest land was negative, which means its area decreased continuously during the period of 1992 to 2018, while the area of construction land increased continuously in that K > 0. Although the K of the unused land varied greatly, its area change performance was not obvious in LULC maps because of its small absolute area. Furthermore, the K of unused land reached its maximum value in 2001–2009, which indicated the remarkable area change during this period. It is suggested that the reason for the significant change in the area of unused land was caused by the 2008 Wenchuan earthquake [38]. The intense earthquakes caused geological changes, causing the landscape to become fragmented, and secondary disasters such as landslides, collapses, and mudslides led to an increase in unused land. Table 2 shows the largest land-use dynamic degree S showing largest changes in LULC during the period of 1992–2001.

3.1.2. Spatial Patterns of LULC Changes during the Period of 1992–2018

To demonstrate the spatial patterns of LULC change hotspots for each LULC type during the period of 1992–2018, the grid-based variation intensity was calculated based on the ArcGIS 10.3 software neighborhood analysis tool. Here, the transfer map of each LULC type was re-divided into non-overlapping grids 1.8 km × 1.8 km, and the LULC change information in each grid was counted [10]. Figure 3 shows the spatial patterns of the change hotspots.
During this period, the most notable features show a massive expansion of urban areas and the shrinkage of farmland around the original urban area, as well as the significant reduction in forest area in the eastern hilly areas (Figure 3). Among them, the annual growth rate of construction land area was 20.84%, causing a total growth of 184,574.79 hm2 in the period of 1992–2018, which was 4.8 times the total area of grassland, water, and unused land in 1992. Similarly, the grassland area enhanced by 7874.19 hm2 with an annual growth rate of 1.36%. Conversely, compared to the initial area in 1992, the loss of forest land in 2018 was the highest 7.92%, followed by the loss of farmland by 7.70%, while water area showed a minimal decline of 0.36%.
The reduced farmland area was mainly in Shuangliu (16,815.15 hm2), Xindu (12,430.71 hm2), Pidu (11,053.62 hm2), Qionglai (9622.80 hm2), and Pengzhou (8035.56 hm2), which are all distributed around initial urban areas, accounting for 62.54% of the total converted farmland. The hotspots map shows the primary procedure of LULC changes was the transition between farmland, forest land, and construction land. The hotspots of farmland were dominated by the expansion of construction land (Figure 3). During the study period, the ecologically important types of LULC in Chengdu decreased due to the reduction in natural vegetation and water areas, which means that the impact of human activities on natural ecosystems was more serious.

3.2. Assessing Changes in Ecosystem Services in Chengdu

3.2.1. Changes in Total ESV

The ecosystem service values (ESV) and changes in total ESV (Table 4) were evaluated in Chengdu for the years 1992, 2001, 2009, and 2018 using modified coefficients. On a whole, the total ESVs of the entire study landscape were about US $1476.3, $1410.2, $1335.1, and $1190.6 million in 1992, 2001, 2009, and 2018, respectively (Table 4). As a result, the value of ecosystem services in Chengdu lost US $285.70 million, accounting for 19.35% of the total value of ecosystem services in 1992. Although the ecosystem service value of each LULC type varied across the entire study landscape for different reference years, the values showed a similar sequence (Figure 4). The forest land showed the highest proportion, while the percentage of unused land was lowest. For example, at the end of the study period (2018), forest land accounted for the largest section, i.e., US $744.15 million (62.50%), while farmland, water area, and grassland accounted for about US $389.39 million (32.71%), $32.16 million (2.70%), and $24.86 million (2.09%), respectively, of the total ESV across the entire landscape.
The changes in ESV demonstrated an obvious reduction during the first (1992–2001), second (2001–2009), third (2009–2018), and entire (1992–2018) periods. The loss in total ESV over the first study period was about US $66.08 million, corresponding to a change rate of 4.48%. The total ESV dropped further, amounting to approximately US $75.07 million and $144.50 million of the total ESV in 2001 and 2009 over the second and third periods, respectively. During the period of 1992–2018, the total ESV lost approximately US $285.70 million, with a change rate of 19.35%. Moreover, the ESV of each LULC type changed, with notable reductions in the ESVs of forest and farmland, which were reduced by US $222.80 and $57.79 million, with change rates of 23.04% and 12.92%, respectively, during the period of 1992–2018. Although the loss in ESV in the water area was less than that of the farmland and forest land, amounting to US $11.62 million, its change rate reached 26.48%. In contrast, the ESV of grassland increased by US $6.51 million, equivalent to 35.33% of the grassland ESV in 1992.

3.2.2. Spatial Patterns of Individual Ecosystem Service Function Changes

In order to reveal the regional differences in the ecosystem service functions, the spatial patterns of the change rates of ecosystem service function values in different districts of Chengdu are shown in Figure 5, using the method of natural breaks. The value of all ecosystem service functions, which included food production, raw material, gas regulation, climate regulation, hydrological regulation, waste treatment, soil conservation, biodiversity protection, and recreation and culture (Table 1), were reduced significantly throughout the entire study area (Figure 2). The change rate was highest (−82.35%) in Wuhou for total ESV in different districts, while Qionglai shows the lowest ESV reduction rate (−1.94%). Across the ecosystem service function values of the entire study area, the change rate varied from −14.22% to −21.17% for the ESVs of food production and recreation and culture, respectively. The decline in the value of ecosystem services was mainly attributed to the reduction in the area of farmland and forest land. However, the increase in construction land did not increase the value of ecosystem services, because its value coefficient was zero, which ultimately led to a dramatic decrease in the value of individual ecosystem service functions in each district.
The spatial distribution of ecosystem service function values presented an obvious circular structure, demonstrating a drastic reduction zone of the central city (in Jinniu, Wuhou, Jinjiang, Chenghua, and Qingyang) as the first circle, a reduction zone around the urban zone as the second circle, and the outermost zone as the third circle. For example, in terms of the change rate of waste treatment (Figure 5f), located in the third circles of Pengzhou, Dujiangyan, Chongzhou, Dayi, Qionglai, and Pujiang in western Chengdu, and Jintang in the east, the rates of change were smaller than those of the first and second circles. Except for the increase in food production in Pujiang (8.45%), raw material in Qionglai (0.91%), and recreation and culture in Qionglai (0.09%), the value of ecosystem service functions in other districts decreased, among which the highest reduction rate was in Wuhou for raw material (−86.68%). Furthermore, the ESV gap in different districts was obvious (Figure 6). In 1992, the ESV ranged from US $3.15 million (in Jinjiang) to US $210.08 million (in Dayi), compared to the range of US $1.01 million (in Qingyang) to US $188.25 million (in Dujiangyan) in 2018.

3.2.3. ESV Change in Different Ecosystem Types from 1992 to 2018

To show the ESV change in different ecosystem types and the impact of ecosystem service value coefficients on ESV, we compared the differences in ESV for different LULC types using three coefficients for the period of 1992–2018 (Table 5, Figure 7), based on coefficients modified by NPP (marked as ESVa), and proposed by Costanza et al. in 1997 [4] (marked as ESVb) and Costanza et al. in 2014 [31] (marked as ESVc). Due to the same study area, the total ESV calculated by the three coefficients showed the same trend. From the ESVb, the ecosystem service value of forest land continued decreasing during the period of 1992–2018, reducing from US $966.93 million to $744.15 million, while the ESVs of farmland and water area generally decreased in the corresponding period, reducing from US $447.18 million to $389.39 million, and $43.78 million to $32.16 million, respectively, although there was a slight increase from 1992–2001. Meanwhile, the ESV of unused land increased sharply over the period of 2001–2009, although there was almost no change during the period of 1992–2001, while the ESV of grassland decreased over the period of 1992–2001 and increased during the period of 2001–2018 with an overall increase of US $6.5 million. The ESVc was always found to be greater than ESVa and ESVb, except for the value of the ecosystem services of the unused land (Figure 7). The ESVb values of farmland, forest land, and grassland were less than ESVa and ESVc values, while the ESVb of water area was between that of ESVa and ESVc (Figure 7).
The ESV showed variations during the period of 1992 to 2018. The total ESVs were about US $845.16 million, $803.78 million, $742.68 million, and $656.96 million in 1992, 2001, 2009, and 2018, respectively, using the global coefficients of Costanza et al. in 1997 [4], whereas the total ESV changed to US $6091.34 million, $6182.27 million, $6238.89 million, and $6420.33 million in corresponding years when using global coefficients developed by Costanza et al. in 2014 [31]. These estimates were 0.5 and four times that of the estimates obtained using the modified coefficients.

3.2.4. Elasticity of Ecosystem Service Value Change in Response to LULC Changes

The elasticity of the ecosystem service value change responding to LULC changes during the periods of 1992–2001, 2001–2009, and 2009–2018 were 0.28, 0.34, and 0.50, respectively, indicating that a transition of 1% of land area would cause average changes of 0.28%, 0.34%, and 0.50% in ESV, respectively (Figure 8). Using coefficients proposed by Costanza et al. in 1997 [4] and in 2014 [31], the elasticity changed to 0.31 and 0.09 during the period of 1992–2001, 0.49 and 0.06 for the period of 2001–2009, and 0.54 and 0.14 for the period of 2009–2018, respectively. Some studies indicated that the elasticity of ESV with respect to LULC in China for the periods of 1988–2000 and 2000–2008 were 0.12 and 0.33, respectively, using the unit value from Costanza et al. in 2014 [30]. Although the elasticity increased notably, the response of ESV to LULC changes was still not obvious in China [30]. High elasticity, where ESV changes were sensitive to LULC, mainly situated in areas with low ESV values, were concentrated in the central part of the North China Plain, northeastern China, and central northwestern China [30].
In general, elasticity continued increasing, indicating that ecosystem services were increasingly responsive to LULC changes in Chengdu. The higher the elasticity, the greater the likelihood that a slight change in LULC can result in a considerable change in the ESV.

4. Discussion

The present study quantified and mapped the LULC changes and their influence on the ecosystem in Chengdu from 1992 to 2018. Results demonstrated that urban expansion is the main reason for the high level of total ESV loss at the expense of farmland and forest land coverage. The trend of LULC changes is in conformity with other findings in Chengdu, but the trend and absolute value of total ESV were significantly different, even resulting in opposite conclusions [36,38]. According to Li et al. [38], the ESV of Chengdu increased by 75.46% from 2000–2015, which was calculated using a corrected equivalent factor of ESV with a marginal value correction factor and CPI (Consumer Price Index) accumulation coefficient. This method obtained a table of ESV coefficients that increased over time; thus, despite the notable reduction in the area of farmland and forest land, the increased coefficient compensated for the decrease in ESV due to the reduction in area, and finally concluded that the ESV increased significantly. Peng et al. showed that the ESV of Chengdu was mainly composed of farmland and forest land, and it continuously declined substantially from 1978 to 2010, which was consistent with the findings of this paper [36]. However, the above studies did not consider the spatial heterogeneity of the LULC changes and the spatial pattern of individual ecosystem service functions in Chengdu, which were the focus of this paper.
The loss in ESV caused by an excessive pursuit of economic benefits is a very common phenomenon in China, especially in the process of regional development and urban expansion. A number of studies on the impact of land-use change on ESV were conducted elsewhere in China and yielded meaningful results [24,55,60,63]. Hu et al. reported that the huge loss of ecosystem services in Xishuangbanna was attributed to the sudden shift in land use from tropical forests and swidden fields to large rubber plantations from 1988–2006, and they pointed out that the cost of a GDP increase of US $1 corresponds to at least a US $1.39 reduction in ESV in the agricultural economy [22]. Wang et al. indicated that, due to population growth and accompanying food requirements, the expansion of cultivated areas through the zealous exploitation of grassland and marginal forest land caused tremendous damage to the ecological environment, such as the decline of ecosystem services and land degradation [61]. Zhang et al. explored the commonalities and differences in the influence of LULC on ESV between three coastal urban agglomerations and demonstrated that, if a nature reserve is established and protection laws and policies are applied, it is not inevitable that urban expansion will result in a net reduction in ESV [63].
Globally, due to historical and policy reasons, there are completely different land-use structures and trends [64,65,66,67,68,69]. Many studies in European countries indicated that the growth of forest land and abandonment of agricultural areas are common issues across Europe, which had a great impact on biodiversity and landscape [70,71]. For example, in Poland and Hungary from 2002–2016, the trend of LULC changes could be depicted with a quadratic function using a statistical method, and the trend of agricultural land was negative while the trend of forest land was positive [66,68]. This is contrary to the changes of forest land area for some cities mentioned above in China. However, based on a case of a great agricultural region in the Czech Republic from 1845–2010, a study revealed that the spread of arable land and agricultural intensification led to an increase in provisioning services, while cultural and regulation services were remarkably reduced by analyzing long-term LULC data. However, in most developing countries such as Ethiopia, the loss rate of ecosystem services is high, which was interpreted as an increase in farmland and settlements, but a decrease in forest land [64].
There are several uncertainties and limitations in the ESV assessment. For instance, the benefits transfer method which was adopted in this study assumes the homogeneity of ESV and no change within each LULC class [30,33]. However, the heterogeneity and complexity of human–environmental systems can cause inevitable errors by value generalization and transfer [30]. Also, this approach is not reliable and valid until building the empirical relationships required between ecosystem characteristics and services [72,73]. Despite many researches commenting on the impact of accurate value coefficients for ESV evaluation, it is often rarely questioned for dynamic analyses as opposed to cross-sectional analyses, because the accuracy of value coefficients affects the ESV at certain points in time more than the estimates of directional change of ESV for a time series [33,63]. In addition, a number of studies calculated the coefficient of sensitivity (CS) to validate the reliability of research results and the value coefficients (VCs) of ESV were adjusted by ±50% in most cases [22,33,47,64]. However, using the CS method to test the trustworthiness of evaluation results was criticized, because the value of CS is often less than 1 even when VCs were adjusted by ±25%, which indicates erroneously robust coefficients [74]. Hence, this study did not carry out CS analysis.

5. Conclusions

This study emphasizes the important link between LULC changes and influences on the function and structure of ecosystem services. It is vital for Chengdu to explore how ecosystem services respond to LULC transformation and to ensure sustainable urban development. We believe that this study can serve as a reference and basis for improving decision-making involving the management of land resources, and contribute to a trade-off between urban expansion and the reduction in ecosystem services.
In this study, we used net primary productivity to calibrate the ecosystem service value equivalents per unit area of China’s terrestrial ecosystems, and thus, established an ESV assessment model for Chengdu, and analyzed the impact of LULC changes on ESV during the period of 1992 to 2018. Remote sensing, map visualization, and geographic information system (GIS) technologies were applied to obtain LULC datasets and spatial patterns for the changes in LULC and ESV. Our results indicate that a decline in farmland, woodland, and water area, and a rapid expansion of construction land over for the period of 1992–2018 resulted in a continuous loss of total ESV to the tune of US $285.7 million using modified coefficients, and US $188.20 million and $328.99 million using global value coefficients in 1997 and 2014, respectively. The reduction in farmland which was mainly distributed around urban areas and forest land of the entire study area was the main contributor to the noteworthy loss of ESV. Furthermore, during the entire study period, the value of individual ecosystem functions also declined, with the maximum and minimum change rates being for food production (−14.22%) and entertainment culture (−21.17%). We also calculated the elasticity, and results showed that elasticity continued increasing significantly from 0.28 to 0.50, indicating that ecosystem services were increasingly responsive to LULC changes in Chengdu.
The decline in the value of ecosystem services reflects, to some extent, the impact of ecological degradation in Chengdu. Under the current rapid urban development model, urban land use should focus on improving land-intensive use efficiency, rather than sacrificing farmland and forest land to expand urban areas. Moreover, it is necessary to adopt appropriate ecological protection measures to achieve the balance between economic development and ecological environment health. The experience of Liaodong Peninsula which maintained ecosystem services by building a nature reserve can be used as a model [63]. In China, it is becoming increasingly urgent to realize the utmost importance of ecosystem conservation in the current process of rapid urbanization. The ecosystem service value assessment at the regional level should be regarded as a conservation strategy to promote decision-making for resource conservation and sustainable utilization. Beyond the completion of these findings, we recommend that a further challenge will be to continuously put forward specific alternative strategies and future planning for improving the ecological environment and ecological services, and for reducing the impact of human activities in the study area and elsewhere in China.

Author Contributions

Design and investigation of this research, X.L., M.X., C.C., and R.S.; collecting and analyzing data, X.L.; manuscript writing, X.L. and M.X.; manuscript revision, C.C., R.S., W.C., and H.J.

Acknowledgments

The authors greatly acknowledge the financial support from the National Key Research and Development Program of China (Grant No. 2017YFD0600903) and Technological Basic Special Funds (No. 2013FY1118063). The author is very grateful to reviewers for their helpful and constructive comments on the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. MEA. Ecosystems and Human Well-being: Current State and Trends; Island Press: Washington, DC, USA, 2003. [Google Scholar]
  2. De Groot, R.S.; Wilson, M.A.; Boumans, R.M.J. A typology for the classification, description and valuation of ecosystems functions, goods and services. Ecol. Econ. 2002, 41, 393–408. [Google Scholar] [CrossRef]
  3. Faber, J.H.; Van, W.J. Elaborations on the use of the ecosystem services concept for application in ecological risk assessment for soils. Sci. Total. Environ. 2012, 415, 3–8. [Google Scholar] [CrossRef] [PubMed]
  4. Costanza, R.; D’Arge, R.; Groot, R.D.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  5. Gómez-Baggethun, E.; Barton, D.N. Classifying and valuing ecosystem services for urban planning. Ecol. Econ. 2013, 86, 235–245. [Google Scholar] [CrossRef]
  6. MEA. Ecosystems and Human Well-Being: Synthesis; Island Press: Washington, DC, USA, 2005. [Google Scholar]
  7. Schägner, J.P.; Brander, L.; Maes, J.; Hartje, V. Mapping ecosystem services’ values: Current practice and future prospects. Ecosyst. Servi. 2013, 4, 33–46. [Google Scholar] [CrossRef]
  8. TEEB. The Economics of Ecosystems and Biodiversity: Ecological and Economic Foundations; Earthscan: London, UK, 2010. [Google Scholar]
  9. De Groot, R.S.; Alkemade, R.; Braat, L.; Hein, L.; Willemen, L. Challenges in integrating the concept of ecosystem services and values in landscape planning, management and decision making. Ecol. Complex. 2010, 7, 260–272. [Google Scholar] [CrossRef]
  10. 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]
  11. Chen, B.; Chen, G.Q. Modified ecological footprint accounting and analysis based on embodied exergy—A case study of the Chinese society 1981–2001. Ecol. Econ. 2007, 61, 355–376. [Google Scholar] [CrossRef]
  12. Kang, H.; Seely, B.; Wang, G.; Innes, J.; Zheng, D.; Chen, P.; Wang, T.; Li, Q. Evaluating management tradeoffs between economic fiber production and other ecosystem services in a Chinese-fir dominated forest plantation in Fujian Province. Sci. Total. Environ. 2016, 557, 80–90. [Google Scholar] [CrossRef] [PubMed]
  13. Troy, A.; Wilson, M.A. Mapping ecosystem services: Practical challenges and opportunities in linking GIS and value transfer. Ecol. Econ. 2007, 60, 435–449. [Google Scholar] [CrossRef]
  14. Barral, M.P.; Oscar, M.N. Land-use planning based on ecosystem service assessment: A case study in the Southeast Pampas of Argentina. Agr. Ecosyst. Environ. 2012, 154, 34–43. [Google Scholar] [CrossRef]
  15. Bateman, I.J.; Harwood, A.R.; Mace, G.M.; Watson, R.T.; Abson, D.J.; Andrews, B.; Binner, A.; Crowe, A.; Day, B.H.; Dugdale, S. Bringing ecosystem services into economic decision-making: Land use in the United Kingdom. Science 2013, 341, 45–50. [Google Scholar] [CrossRef] [PubMed]
  16. Frélichová, J.; Vačkář, D.; Pártl, A.; Loučková, B.; Harmáčková, Z.V.; Lorencová, E. Integrated assessment of ecosystem services in the Czech Republic. Ecosyst. Serv. 2014, 8, 110–117. [Google Scholar] [CrossRef]
  17. Guo, Z.; Xiao, X.; Gan, Y.; Zheng, Y. Ecosystem functions, services and their values—A case study in Xingshan County of China. Ecol. Econ. 2001, 38, 141–154. [Google Scholar] [CrossRef]
  18. Konarska, K.M.; Sutton, P.C.; Castellon, M. Evaluating scale dependence of ecosystem service valuation: A comparison of NOAA-AVHRR and Landsat TM datasets. Ecol. Econ. 2002, 41, 491–507. [Google Scholar] [CrossRef]
  19. Liu, S.; Costanza, R.; Troy, A.; D’Aagostino, J.; Mates, W. Valuing New Jersey’s ecosystem services and natural capital: A spatially explicit benefit transfer approach. Environ. Manage. 2010, 45, 1271–1285. [Google Scholar] [CrossRef] [PubMed]
  20. Li, R.Q.; Dong, M.; Cui, J.Y.; Zhang, L.L.; Cui, Q.G.; He, W.M. Quantification of the impact of land-use changes on ecosystem services: A case study in Pingbian County, China. Environ. Monit. Assess. 2007, 128, 503–510. [Google Scholar] [CrossRef] [PubMed]
  21. Polasky, S.; Nelson, E.; Pennington, D.; Johnson, K.A. The Impact of Land-Use Change on Ecosystem Services, Biodiversity and Returns to Landowners: A Case Study in the State of Minnesota. Environ. Resour. Econ. 2011, 48, 219–242. [Google Scholar] [CrossRef]
  22. Hu, H.; Liu, W.; Cao, M. Impact of land use and land cover changes on ecosystem services in Menglun, Xishuangbanna, Southwest China. Environ. Monit. Assess. 2008, 146, 147–156. [Google Scholar] [CrossRef] [PubMed]
  23. Zang, S.; Wu, C.; Liu, H.; Na, X. Impact of urbanization on natural ecosystem service values: A comparative study. Environ. Monit. Assess. 2011, 179, 575–588. [Google Scholar] [CrossRef] [PubMed]
  24. Zhao, B.; Kreuter, U.; Li, B.; Ma, Z.; Chen, J.; Nakagoshi, N. An ecosystem service value assessment of land-use change on Chongming Island, China. Land Use Policy 2004, 21, 139–148. [Google Scholar] [CrossRef]
  25. Tong, C.; Feagin, R.A.; Lu, J.; Zhang, X.; Zhu, X.; Wang, W.; He, W. Ecosystem service values and restoration in the urban Sanyang wetland of Wenzhou, China. Ecol. Eng. 2007, 29, 249–258. [Google Scholar] [CrossRef]
  26. Kreuter, U.P.; Harris, H.G.; Matlock, M.D.; Lacey, R.E. Change in ecosystem service values in the San Antonio area, Texas. Ecol. Econ. 2001, 39, 333–346. [Google Scholar]
  27. Kubiszewski, I.; Costanza, R.; Dorji, L.; Thoennes, P.; Tshering, K. An initial estimate of the value of ecosystem services in Bhutan. Ecosyst. Serv. 2013, 3, e11–e21. [Google Scholar] [CrossRef]
  28. Lambin, E.F.; Meyfroidt, P. Global land use change, economic globalization, and the looming land scarcity. Proc. Natl. Acad. Sci. USA 2011, 108, 3465–3472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Gibbs, H.K.; Ruesch, A.S.; Achard, F.; Clayton, M.K.; Holmgren, P.; Ramankutty, N.; Foley, J.A. Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proc. Natl. Acad. Sci. USA 2010, 107, 16732–16737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. 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]
  31. 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]
  32. Su, S.; Li, D.; Hu, Y.; Xiao, R.; Zhang, Y. Spatially non-stationary response of ecosystem service value changes to urbanization in Shanghai, China. Ecologi Indic. 2014, 45, 332–339. [Google Scholar] [CrossRef]
  33. Kindu, M.; Schneider, T.; Teketay, D.; Knoke, T. Changes of ecosystem service values in response to land use/land cover dynamics in Munessa-Shashemene landscape of the Ethiopian highlands. Sci. Total. Environ. 2016, 547, 137–147. [Google Scholar] [CrossRef] [PubMed]
  34. Xie, G.; Lu, C.X.; Leng, Y.F.; Zheng, D.U.; Li, S.C. Ecological assets valuation of the Tibetan Plateau. J. Nat. Resour. 2003, 18, 189–196. [Google Scholar]
  35. Lu, X.; Gan, Q.; Shi, L. Research on the Chengdu Electronic Supervision System for Food Quality Safety Based on PDA and Mobile GIS Technology; Springer-Verlag: Berlin/Heidelberg, Germany, 2012; pp. 391–397. [Google Scholar]
  36. Peng, W.F.; Zhou, J.M.; Fan, S.Y.; Yang, C.J. Effects of the Land Use Change on Ecosystem Service Value in Chengdu, Western China from 1978 to 2010. J. Indian Soc. Remote Sens. 2016, 44, 197–206. [Google Scholar] [CrossRef]
  37. Qin, B. City profile: Chengdu. Cities 2015, 43, 18–27. [Google Scholar] [CrossRef]
  38. Li, Y.; Zhan, J.; Liu, Y.; Zhang, F.; Zhang, M. Response of ecosystem services to land use and cover change: A case study in Chengdu City. Resour. Conserv. Recy. 2018, 132, 291–300. [Google Scholar] [CrossRef]
  39. Peng, W.F.; Zhou, J.M.; Xinliang, X.U.; Luo, H.L.; Zhao, J.F.; Yang, C.J. Effect of land use changes on carbon emission and its spatial patterns in Chengdu Plain and its surrounding area, western China, from 1990 to 2010. Ecol. Sci. 2017, 36, 105–114. [Google Scholar]
  40. Ma, W.L.; Ren, P.; Zeng, Y.Q. The Characteristics for the Spatial-Temporal Evolution of Water Conservation Value in Farmland Ecosystem in Chengdu Plain Economic Zone. China Land Sci. 2015, 10, 85–94. [Google Scholar]
  41. Teferi, E.; Uhlenbrook, S.; Bewket, W.; Wenninger, J.; Simane, B. The use of remote sensing to quantify wetland loss in the Choke Mountain range, Upper Blue Nile basin, Ethiopia. Hydrol. Earth Syst. Sci. Discuss. 2010, 14, 2415–2428. [Google Scholar] [CrossRef] [Green Version]
  42. Gebrehiwot, S.G.; Bewket, W.; Gärdenäs, A.I.; Bishop, K. Forest cover change over four decades in the Blue Nile Basin, Ethiopia: comparison of three watersheds. Reg. Environ. Chang. 2014, 14, 253–266. [Google Scholar] [CrossRef]
  43. Gashaw, T.; Tulu, T.; Argaw, M.; Worqlul, A.W. Evaluation and prediction of land use/land cover changes in the Andassa watershed, Blue Nile Basin, Ethiopia. Environ. Syst. Res. 2017, 6, 17. [Google Scholar] [CrossRef]
  44. Lin, C.; Wu, C.C.; Tsogt, K.; Ouyang, Y.C.; Chang, C.I. Effects of atmospheric correction and pansharpening on LULC classification accuracy using WorldView-2 imagery. Inf. Process. Agric. 2015, 2, 25–36. [Google Scholar] [CrossRef]
  45. Tehrany, M.S.; Pradhan, B.; Jebuv, M.N. A comparative assessment between object and pixel-based classification approaches for land use/land cover mapping using SPOT 5 imagery. Geocarto Int. 2014, 29, 351–369. [Google Scholar] [CrossRef]
  46. Al-Sharif, A.A.A.; Pradhan, B. Monitoring and predicting land use change in Tripoli Metropolitan City using an integrated Markov chain and cellular automata models in GIS. Arab. J. Geosci. 2014, 7, 4291–4301. [Google Scholar] [CrossRef]
  47. 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]
  48. Ahmad, A.; Quegan, S. Analysis of maximum likelihood classification on multispectral data. Appl. Math. Sci. 2012, 6, 6425–6436. [Google Scholar]
  49. Tilahun, A. Accuracy Assessment of Land Use Land Cover Classification using Google Earth. J. Environ. Prot. 2015, 4, 193–198. [Google Scholar] [CrossRef]
  50. Peng, W.; Wang, G.; Zhou, J.; Zhao, J.; Yang, C. Studies on the temporal and spatial variations of urban expansion in Chengdu, western China, from 1978 to 2010. Sustain. Cities Soc. 2015, 17, 141–150. [Google Scholar] [CrossRef]
  51. Hao, F.; Lai, X.; Ouyang, W.; Xu, Y.; Wei, X.; Song, K. Effects of land use changes on the ecosystem service values of a reclamation farm in northeast China. Environ. Manage. 2012, 50, 888–899. [Google Scholar] [CrossRef] [PubMed]
  52. Wang, X.; Bao, Y. Study on the methods of land use dynamic change research. Prog. Geogr. 1999, 18, 81–87. [Google Scholar]
  53. 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]
  54. Xie, G.D.; Zhen, L.; Lu, C.X.; Xiao, Y.; Chen, C. Expert Knowledge Based Valuation Method of Ecosystem Services in China. J. Nat. Res. 2008, 23, 911–919. [Google Scholar]
  55. Li, F.; Zhang, S.; Yang, J.; Bu, K.; Wang, Q.; Tang, J.; Chang, L. The effects of population density changes on ecosystem services value: A case study in Western Jilin, China. Ecol. Indic. 2016, 61, 328–337. [Google Scholar]
  56. Zhao, T.; Ouyang, Z.; Zheng, H.; Wang, X.; Miao, H. Forest ecosystem services and their valuation in China. J. Nat. Res. 2004, 4, 480–491. [Google Scholar]
  57. Zhang, M.L. An improved CASA model based on the grassland comprehensive and Sequential Classification System and its Application to the Grassland NPP in China; Gansu Agricultural University: Gansu, China, 2012. [Google Scholar]
  58. Xie, G.D.; Xiao, Y.; Zhen, L. Study on ecosystem services value of food production in China. Chin. J. Eco-Agric. 2005, 13, 10–13. [Google Scholar]
  59. Ai, Z.H. Impact of Land-Use Changes on the Ecosystem Service Value in Qianjiang City; Central China Normal University: Wuhan, China, 2013. [Google Scholar]
  60. Zhang, Z.; Gao, J.; Gao, Y. The influences of land use changes on the value of ecosystem services in Chaohu Lake Basin, China. Environ. Earth Sci. 2015, 74, 385–395. [Google Scholar] [CrossRef]
  61. Wang, Z.L.; Wang, Z.M.; Zhang, B.; Lu, C.; Ren, C. Impact of land use/land cover changes on ecosystem services in the Nenjiang River Basin, Northeast China. Ecol. Process. 2015, 4, 11. [Google Scholar] [CrossRef]
  62. Wang, Y.; Gao, J.; Wang, J.; Qiu, J. Value assessment of ecosystem services in nature reserves in Ningxia, China: Aresponse to ecological restoration. Plos ONE 2014, 9, e89174. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, Y.; Zhao, L.; Liu, J.; Liu, Y.; Li, C. The Impact of Land Cover Change on Ecosystem Service Values in Urban Agglomerations along the Coast of the Bohai Rim, China. Sustainability 2015, 7, 10365–10387. [Google Scholar] [CrossRef] [Green Version]
  64. Tolessa, T.; Senbeta, F.; Kidane, M. The impact of land use/land cover change on ecosystem services in the central highlands of Ethiopia. Ecosyst. Serv. 2017, 23, 47–54. [Google Scholar] [CrossRef]
  65. Sharma, R.; Nehren, U.; Rahman, S.A.; Meyer, M.; Rimal, B.; Seta, G.A.; Baral, H. Modeling Land Use and Land Cover Changes and Their Effects on Biodiversity in Central Kalimantan, Indonesia. Land 2018, 7, 57. [Google Scholar] [CrossRef]
  66. Noszczyk, T.; Rutkowska, A.; Hernik, J. Determining Changes in Land Use Structure in Małopolska Using Statistical Methods. Pol. J. Environ. Stud. 2017, 26, 211–220. [Google Scholar] [CrossRef] [Green Version]
  67. Frélichová, J.; Fanta, J. Ecosystem service availability in view of long-term land-use changes: A regional case study in the Czech Republic. Ecosyst. Health Sustain. 2016, 1, 1–15. [Google Scholar] [CrossRef]
  68. Cegielska, K.; Noszczyk, T.; Kukulska, A.; Szylar, M.; Hernik, J.; Dixon-Gough, R.; Jombach, S.; Valánszki, I.; Kovács, K.F. Land use and land cover changes in post-socialist countries: Some observations from Hungary and Poland. Land Use Policy 2018, 78, 1–18. [Google Scholar] [CrossRef]
  69. Cabral, P.; Feger, C.; Levrel, H.; Chambolle, M.; Basque, D. Assessing the impact of land-cover changes on ecosystem services: A first step toward integrative planning in Bordeaux, France. Ecosyst. Serv. 2016, 22, 318–327. [Google Scholar] [CrossRef]
  70. Muchová, Z.; Tárníková, M. Land cover change and its influence on the assessment of the ecological stability. Appl. Ecol. Environ. Res. 2018, 16, 2169–2182. [Google Scholar] [CrossRef]
  71. Baumann, M.; Kuemmerle, T.; Elbakidze, M.; Ozdogan, M.; Radeloff, V.C.; Keuler, N.S.; Prishchepov, A.V.; Kruhlov, I.; Hostert, P. Patterns and drivers of post-socialist farmland abandonment in Western Ukraine. Land Use Policy 2011, 28, 552–562. [Google Scholar] [CrossRef]
  72. Wong, C.P.; Jiang, B.; Kinzig, A.P.; Lee, K.N.; Ouyang, Z. Linking ecosystem characteristics to final ecosystem services for public policy. Ecol. Lett. 2015, 18, 108–118. [Google Scholar] [CrossRef] [PubMed]
  73. Richardson, L.; Loomis, J.; Kroeger, T.; Casey, F. The role of benefit transfer in ecosystem service valuation. Ecol. Econ. 2015, 115, 51–58. [Google Scholar] [CrossRef]
  74. 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]
Figure 1. Location map of Chengdu City. The letters in the figure represent the abbreviations of the district names.
Figure 1. Location map of Chengdu City. The letters in the figure represent the abbreviations of the district names.
Sustainability 10 03580 g001
Figure 2. Land-use/land-cover (LULC) patterns of Chengdu in the years 1992, 2001, 2009, and 2018.
Figure 2. Land-use/land-cover (LULC) patterns of Chengdu in the years 1992, 2001, 2009, and 2018.
Sustainability 10 03580 g002
Figure 3. Hotspots of LULC changes for each type from 1992 to 2018 for (a) farmland, (b) forest land, (c) grassland, (d) water area, (e) construction land, and (f) unused land.
Figure 3. Hotspots of LULC changes for each type from 1992 to 2018 for (a) farmland, (b) forest land, (c) grassland, (d) water area, (e) construction land, and (f) unused land.
Sustainability 10 03580 g003
Figure 4. Proportions of LULC changes (%) and ecosystem service values (ESVs; %) during the period of 1992 to 2018 in Chengdu.
Figure 4. Proportions of LULC changes (%) and ecosystem service values (ESVs; %) during the period of 1992 to 2018 in Chengdu.
Sustainability 10 03580 g004
Figure 5. ESV change rates (%) for individual ecosystem functions of (a) food production, (b) raw material, (c) gas regulation, (d) climate regulation, (e) hydrological regulation, (f) waste treatment, (g) soil conservation, (h) biodiversity protection, and (i) recreation and culture. The letters in the figure represent the abbreviations of the district names.
Figure 5. ESV change rates (%) for individual ecosystem functions of (a) food production, (b) raw material, (c) gas regulation, (d) climate regulation, (e) hydrological regulation, (f) waste treatment, (g) soil conservation, (h) biodiversity protection, and (i) recreation and culture. The letters in the figure represent the abbreviations of the district names.
Sustainability 10 03580 g005
Figure 6. (a) Ecosystem service values (US$ in millions) for each district, and (b) ESV change rates (%) from 1992 to 2018.
Figure 6. (a) Ecosystem service values (US$ in millions) for each district, and (b) ESV change rates (%) from 1992 to 2018.
Sustainability 10 03580 g006
Figure 7. ESVs of different ecosystems of (a) farmland, (b) forest land, (c) grassland, (d) water area, (e) construction land, and (f) unused land (ESVa calculated using own modified coefficients; ESVb calculated using the coefficient proposed by Costanza et al. in 1997 [4]; ESVc calculated using the coefficient proposed by Costanza et al. in 2014 [31]).
Figure 7. ESVs of different ecosystems of (a) farmland, (b) forest land, (c) grassland, (d) water area, (e) construction land, and (f) unused land (ESVa calculated using own modified coefficients; ESVb calculated using the coefficient proposed by Costanza et al. in 1997 [4]; ESVc calculated using the coefficient proposed by Costanza et al. in 2014 [31]).
Sustainability 10 03580 g007
Figure 8. Elasticity of ecosystem service value changes in response to LULC changes based on different unit values.
Figure 8. Elasticity of ecosystem service value changes in response to LULC changes based on different unit values.
Sustainability 10 03580 g008
Table 1. Details of corrected ecosystem service value per unit area for ecosystem service functions of each land ecosystem type. LULC—land use/land cover; US—United States.
Table 1. Details of corrected ecosystem service value per unit area for ecosystem service functions of each land ecosystem type. LULC—land use/land cover; US—United States.
Ecosystem ServicesSub-typesEach LULC Type Ecosystem Service Value Coefficients (US$/hm2/year)
FarmlandForest LandGrasslandWater AreaConstruction LandUnused Land
Provisioning services Food production78.94 27.40 30.42 30.99 01.17
Raw material30.79 247.45 25.47 20.47 02.34
Regulation servicesGas regulation56.84 358.72 106.13 29.82 0 3.51
Climate regulation76.58 337.96 110.38 120.46 0 7.60
Hydrological regulation60.79 339.62 107.55 1097.61 0 4.09
Waste treatment109.73 142.82 93.40 868.38 0 15.20
Supporting servicesSoil conservation116.05 333.81 158.49 23.98 0 9.94
Biodiversity protection80.52 374.50 132.32 200.57 0 23.39
Cultural servicesRecreation and culture13.42 172.72 61.56 259.64 0 14.03
Total623.65 2334.99 825.73 2651.91 0 81.28
Costanza et al. in 1997 [4] 12613383211172700
Costanza et al. in 2014 [31] 5567380041661251266610
Table 2. Changes in LULC in Chengdu during the period of 1992–2018.
Table 2. Changes in LULC in Chengdu during the period of 1992–2018.
Land-Use TypePercentage of Land-Use Type (%)Land-Use Dynamic Index K (%)
19922001200920181992–20012001–20092009–20181992–2018
Farmland59.5662.0058.7151.860.46−0.66−1.30−0.50
Forest land34.4031.7429.7426.47−0.86−0.79−1.22−0.89
Grassland1.850.832.052.50−6.1118.372.441.36
Water area1.371.381.181.010.09−1.84−1.62−1.02
Construction land2.834.058.2618.164.8013.0013.3020.84
Unused land0.0010 0.0026 0.0559 0.0037 16.43260.71−10.389.89
Land use dynamic degree S (%) 29.2029.1225.7011.37
Table 3. Percentage of LULC conversion in Chengdu from 1992 to 2018.
Table 3. Percentage of LULC conversion in Chengdu from 1992 to 2018.
2018 (%)1992 (%)
FarmlandForest LandGrasslandWater AreaConstruction LandUnused Land
Farmland71.0623.6228.1028.9417.520
Forest land4.2268.2623.791.380.7563.04
Grassland0.234.7040.060.140.1236.96
Water area0.380.240.6045.592.140
Construction land24.103.177.4223.9479.470
Unused land0.0010.0080.0040.00200
Direction of change
Table 4. Assessed ecosystem service values (ESVs; US$ in millions) for each LULC type of Chengdu using modified coefficients. C—change; CR—change rate (calculated using Equation (5)).
Table 4. Assessed ecosystem service values (ESVs; US$ in millions) for each LULC type of Chengdu using modified coefficients. C—change; CR—change rate (calculated using Equation (5)).
Land Use Type19922001200920181992–20012001–20092009–20181992–2018
ESVESVESVESVCCR (%)CCR (%)CCR (%)CCR (%)
Farmland447.18465.52440.8389.3918.344.09−24.72−5.31−51.41−11.66−57.79−12.92
Forest land966.93892.24836.2744.15−74.69−7.73−56.04−6.28−92.05−11.01−222.80−23.04
Grassland18.358.2520.3924.86−10.1−54.8912.14145.784.4722.066.5135.33
Water area43.7844.1537.6632.160.370.91−6.49−14.71−5.5−14.59−11.62−26.48
Construction land00000-0-0-0-
Unused land0.0010.00250.05470.00360.00152000.05221733.3−0.051−92.730.0026300
Total1476.31410.21335.11190.6−66.08−4.48−75.07−5.33−144.50−10.82−285.70−19.35
Table 5. Ecosystem service value coefficients of different ecosystem types (US$ in millions).
Table 5. Ecosystem service value coefficients of different ecosystem types (US$ in millions).
LULC TypesFarmlandForest LandGrasslandWater AreaConstruction LandUnused Land
Modified coefficients623.65 2334.99 825.73 2651.91 0 81.28
Costanza et al. in 1997 [4]12613383211172700
Costanza et al. in 2014 [31]556738004166125126s6610

Share and Cite

MDPI and ACS Style

Lin, X.; Xu, M.; Cao, C.; P. Singh, R.; Chen, W.; Ju, H. Land-Use/Land-Cover Changes and Their Influence on the Ecosystem in Chengdu City, China during the Period of 1992–2018. Sustainability 2018, 10, 3580. https://0-doi-org.brum.beds.ac.uk/10.3390/su10103580

AMA Style

Lin X, Xu M, Cao C, P. Singh R, Chen W, Ju H. Land-Use/Land-Cover Changes and Their Influence on the Ecosystem in Chengdu City, China during the Period of 1992–2018. Sustainability. 2018; 10(10):3580. https://0-doi-org.brum.beds.ac.uk/10.3390/su10103580

Chicago/Turabian Style

Lin, Xiaojuan, Min Xu, Chunxiang Cao, Ramesh P. Singh, Wei Chen, and Hongrun Ju. 2018. "Land-Use/Land-Cover Changes and Their Influence on the Ecosystem in Chengdu City, China during the Period of 1992–2018" Sustainability 10, no. 10: 3580. https://0-doi-org.brum.beds.ac.uk/10.3390/su10103580

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