Journal of Geographical Sciences >
The response of key ecosystem services to land use and climate change in Chongqing: Time, space, and altitude
Gao Jie (1986-), Associate Professor, specialized in the ecosystem services, land use change. E-mail: gaojieswu@swu.edu.cn |
Received date: 2021-02-21
Accepted date: 2021-09-17
Online published: 2022-04-25
Supported by
National Natural Science Foundation of China(41701611)
National Natural Science Foundation of China(41830648)
General Program of Social Science and Planning of Chongqing(2020YBZX15)
Fundamental Research Funds for the Central Universities(XDJK2019C090)
Mountainous landscapes are particularly vulnerable and sensitive to climate change and human activities, and a clear understanding of how ecosystem services (ES) and their relationships continuously change over time, across space, and along altitude is therefore essential for ecosystem management. Chongqing, a typical mountainous region, was selected to assess the long-term changes in its key ES and their relationships. From 1992 to 2018, the temporal variation in water yield (WY) revealed that the maximum and minimum WYs occurred in 1998 and 2006, which coincided with El Niño-Southern Oscillation and severe drought events, respectively. Soil export (SE) and WY were consistent with precipitation, which reached their highest values in 1998. During this period, carbon storage (CS) and habitat quality (HQ) both decreased significantly. ES in Chongqing showed large variations in altitude. Generally, WY and SE decreased with increasing altitude, while CS and HQ increased. For spatial distribution, WY and SE showed positive trends in the west and negative trends in the east. In regard to CS and HQ, negative trends dominated the area. Persistent tradeoffs between WY and soil conservation (SC) were found at all altitude gradients. The strong synergies between CS and HQ were maintained over time.
Key words: ecosystem services; InVEST model; Mann-Kendall; tradeoffs; synergies; mountainous region
GAO Jie , BIAN Hongyan , ZHU Chongjing , TANG Shuang . The response of key ecosystem services to land use and climate change in Chongqing: Time, space, and altitude[J]. Journal of Geographical Sciences, 2022 , 32(2) : 317 -332 . DOI: 10.1007/s11442-022-1949-x
Figure 1 Elevation, meteorological stations, main rivers, and the mountainous landscape of Chongqing |
Table 1 Description of the datasets adopted in this study |
Dataset | Source | Description | Spatial resolution or distribution | Temporal resolution |
---|---|---|---|---|
Land-use data | European Space Agency http://maps.elie.ucl.ac.be/CCI/viewer/ | Land use and land cover (LULC) maps on an annual basis from 1992 to 2018 | 300 m×300 m | Yearly |
Climate data | China Meteorological Data Network http://data.cma.cn/ | Maximum, minimum and average temperatures, precipitation, and solar radiation | 12 stations | Daily/ monthly |
Soil data | Harmonized World Soil Database (HWSD) http://webarchive.iiasa.ac.at/Research/LUC/External-World-soil-database/HTML/ | Soil properties including the texture, organic matter content, and root depth | 30 arc-second | - |
Table 2 Methods for ecosystem service assessment |
Ecosystem service indicators | Methods | Unit | Equation | References |
---|---|---|---|---|
Water yield | Water balance equation | m3. ha-1 | ${{Y}_{x}}=\left( 1-\frac{AE{{T}_{x}}}{{{P}_{x}}} \right)\times {{P}_{x}}$ Yx: the water yield for pixel x; Px: the annual precipitation for pixel x; and ATEx: the annual actual evapotranspiration for pixel x. | (Budyko, 1974; Sharp et al., 2014) |
Soil export (inverse indicator of soil conservation) | Universal soil loss equation | t.ha-1.y-1 | USRLx= Rx×Kx×LSx×Cx×Px USRLx: the average annual soil loss for pixel x; Rx: the rainfall factor for pixel x; Kx: the soil erodibility factor for pixel x; LSx: the field topography factor for pixel x; Cx: the cropping and management factor for pixel x; and Px: the factor describing the supporting conservation practices for pixel x. | (Wischmeier and Smith, 1978) |
Carbon storage | Sum of the carbon stored in vegetation, litter and soil | Mg. ha-1 | Cx =Cxabove+Cxbelow+Cxsoil+Cxdead Cxabove, Cxbelow, Cxsoil and Cxdead are the aboveground carbon density, belowground carbon density, soil organic carbon density, and dead organic matter, respectively, for pixel x. | (Tappeiner et al., 2008) |
Habitat quality | Habitat quality model of InVEST | Dimensionless index (0-1) | ${{Q}_{xj}}={{H}_{j}}\times \left( 1-D_{xj}^{z}/(D_{xj}^{z}+{{k}^{z}}) \right)$ Qx: the HQ for pixel x in land-use j; Hj: the habitat suitability of land-use type j; Dxjz: the total threat level in grid cell x of land-use type j; k: the half-saturation value; and z: normalized constant. | (Sharp et al., 2014) |
Figure 2 Trends in cropland, built-up land, grassland, and woodland changes in Chongqing from 1992 to 2018 |
Figure 3 Temporal variation in the ecosystem services in Chongqing from 1992 to 2018 (a. water yield; b. soil export; c. habitat quality; d. carbon storage) |
Table 3 The average values (n=27, mean+SD) and change trends (Z values) of ecosystem services in different altitude gradients. SD represents the standard deviation. |
DEM | |||||
---|---|---|---|---|---|
0-500 m | 500-1000 m | 1000-1500 m | 1500-2000 m | >2000 m | |
Water yield | 673.67±155(0.91)1 | 665.23 ±138(1.08) | 588.81±135(0.70) | 546.59±137(0.041) | 560.50±144(-0.29) |
Soil export | 88.97±13(0.96) | 88.51±12(0.91) | 42.06±6(0) | 23.06±3(0.29) | 28.87±4(0.75) |
Carbon storage | 1315.17±41(-6.56) | 1766.79±53(3.95) | 2268.65±68(1.03) | 2434.68±74(-6.59) | 2380.98±72(-5.67) |
Habitat quality | 0.46±0.0035(-6.71) | 0.64±0.0032(4.34) | 0.87±0.0024(1.18) | 0.95±0.0026(-6.73) | 0.92±0.0034(-5.73) |
1 The values in parentheses indicate the Mann-Kendall test results. The bold numbers indicate significant trends at α=0.05. |
Figure 4 Spatial distributions of average annual ecosystem services in Chongqing from 1992 to 2018 (WY: water yield; SE: soil export, CS: carbon storage; HQ: habitat quality) |
Figure 5 Spatial distribution of the annual ecosystem service trends in Chongqing from 1992 to 2018 (a. water yield; b. soil export; c. carbon storage; d. habitat quality). The black circles indicate significant trends at α=0.05. MK refers to the Mann-Kendall test. |
Figure 6 Variations in the relationships between the ecosystem services in Chongqing from 1992 to 2018 (WY: water yield; SE: soil export; CS: carbon storage; HQ: habitat quality). WY vs. HQ indicates the relationship between water yield and habitat quality. Y-axis represents the Pearson correlation coefficient. |
Figure 7 Variations in the relationships between the ES at different attitude gradients in Chongqing. WY vs. HQ indicates the relationship between water yield and habitat quality. Y-axis represents the Pearson correlation coefficient. |
Table 4 Basic statistical properties of precipitation and different land use types in Chongqing |
Annual average (mm) | SD2 (mm) | CV2 (%) | |
---|---|---|---|
Precipitation | 1122 | 144 | 12.8 |
Annual average (km2) | SD (km2) | CV (%) | |
Cropland | 46067 | 486 | 1.0 |
Woodland | 30724 | 197 | 0.6 |
Grassland | 4083 | 55 | 1.3 |
Built-up land | 718 | 359 | 50.0 |
Water bodies | 804 | 7 | 0.7 |
2 SD represents the standard deviation, and CV represents the coefficient of variation. |
Figure 8 Major socioeconomic, policy, and climate events influencing ES in Chongqing during 1992-2018 (a, b). The urbanization rate from 1992 to 2018 is shown in the middle graph (c). Trends of the built-up land change at low altitudes (d). |
Figure A1 The observed water yield versus simulated water yield (1010 m³/a) |
Figure A2 The observed soil export versus simulated soil export (105 t/a) |
Figure A3 The observed soil organic carbon versus simulated soil organic carbon (107 t/a) |
Table A1 Threats and their maximum distance of influence and weights |
Threats | MAX_DIST | WEIGHT | DECAY |
---|---|---|---|
Express way | 8.0 | 0.5 | Exponential |
Highway | 10 | 0.5 | Exponential |
Railway | 7 | 0.7 | Exponential |
Urban area | 10 | 1 | Exponential |
MAX_DIST: the maximum distance over which each threat affects habitat quality WEIGHT: the impact of each threat on habitat quality, relative to other threats. DECAY: the type of decay over space for the threat. |
Table A2 The sensitivity of habitat types to each threat |
NAME | HABITAT | Express way | Highway | Railway | Urban area |
---|---|---|---|---|---|
Woodland | 1.00 | 0.7 | 0.6 | 0.7 | 0.7 |
Grassland | 0.60 | 0.7 | 0.9 | 0.8 | 0.75 |
Built-up land | 0 | 0 | 0 | 0 | 0 |
Cropland | 0.40 | 0.4 | 0.3 | 0.6 | 0.5 |
Water bodies | 0.7 | 0.5 | 0.6 | 0.5 | 0.7 |
NAME: the name of each land use/land cover type. HABITAT: assign each land use/land cover type a relative habitat suitability score from 0 to 1 where 1 indicates the highest habitat suitability. |
Table B1 Parameters for water yield and soil conservation in Chongqing |
LULC_desc | Kc | root_depth | usle_c | usle_p | LULC_veg |
---|---|---|---|---|---|
Cropland | 0.65 | 2100 | 0.25 | 0.4 | 1 |
Woodland | 1.1 | 7000 | 0.003 | 1 | 1 |
Grassland | 0.58 | 2600 | 0.04 | 0.8 | 1 |
Urban area | 0.3 | 500 | 0 | 0 | 0 |
Water bodies | 1 | 1000 | 0 | 0 | 0 |
LULC_desc: Descriptive name of land use/land cover class. Kc: The plant evapotranspiration coefficient for each LULC class. root_depth (mm): The maximum root depth for vegetated land use classes. usle_c: Cover-management factor for the USLE. This factor is used to measure the inhibitory effect of vegetation cover and management on soil erosion. usle_p: Support practice factor for the USLE. It is the ratio of soil loss with a support practice, such as contouring, stripping, or terracing to soil loss with straight-row farming up and down the slope. LULC_veg: Values should be 1 for vegetated land use types except wetlands, and 0 for all other land use types, including urban, water bodies, etc. |
Table B2 Parameters for carbon storage in Chongqing |
Items | Cropland | Woodland | Grassland | Built-up land | Water bodies |
---|---|---|---|---|---|
Aboveground (Mg/ha) | 19.5 | 35 | 15.1 | 0 | 0 |
Belowground (Mg/ha) | 45.5 | 85.1 | 67.5 | 0 | 0 |
Soil organic (Mg/ha) | 72.4 | 140.6 | 100 | 1.8 | 0 |
Dead organic (Mg/ha) | 2.5 | 30.5 | 1 | 0 | 0 |
Table B3 Parameters for habitat quality in Chongqing |
Threats | ||||||
---|---|---|---|---|---|---|
Habitat | Express way | Highway | Railway | Urban area | ||
Properties of threats | Weight | - | 0.5 | 0.5 | 0.7 | 1 |
Maximum distance of influence | - | 8 | 10 | 7 | 10 | |
Sensitivity of different land cover types | Woodland | 1 | 0.7 | 0.6 | 0.7 | 0.7 |
Grassland | 0.5 | 0.7 | 0.9 | 0.8 | 0.75 | |
Built-up land | 0 | 0 | 0 | 0 | 0 | |
Cropland | 0.4 | 0.4 | 0.3 | 0.6 | 0.5 | |
Water bodies | 0.7 | 0.5 | 0.6 | 0.5 | 0.7 |
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
MEA, 2005. Ecosystems and Human Well-being. Washington, DC, USA: Island Press.
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
/
〈 |
|
〉 |