Journal of Resources and Ecology >
Damage or Recovery? Assessing Ecological Land Change and Its Driving Factors: A Case of the Yangtze River Economic Belt, China
Received date: 2020-09-24
Accepted date: 2020-11-30
Online published: 2021-05-30
Supported by
The National Social Science Fund of China(19BGL283)
The National Natural Science Foundation of China(41301619)
Ecological land can provide people with ecological products and ecological services; and it plays an important role in maintaining the health and safety of the ecosystem. With China’s rapid urbanization development, ecological land has been invaded in large quantities, and damaged seriously, even resulting in loses of its ecological function. Based on land use data from 1995 to 2015, our study explores the spatial and temporal evolution of the damage or recovery of ecological land in the Yangtze River Economic Belt (YREB). Two spatial models, geographic detector and geographic weighted regression (GWR), were employed to assess the global effects and the local effects of the driving factors for ecological land change, respectively. Our study divided the ecological land change into five types based on the degree of change as severe damage, slight damage, unchanged, slight recovery, and obvious recovery. The results show that from 1995 to 2015, the total area of ecological land in the YREB increased initially and then decreased, but the overall trend was decreasing. The total damaged area was larger than the recovered area. Arable land and woodland both showed downward trends. In terms of ecological land change over the past 20 years, the type of unchanged had the largest area, followed by slight damage and slight recovery. Our study further revealed that ecological land change was the net result of the interaction of many factors, and the explanatory power between any two driving factors was greater than that of any individual driving factor. In addition, driving factors have different impacts on ecological land change in different geographical locations. This knowledge should help land managers and policymakers to be better informed when developing pertinent land use policies at the regional and local levels. The lessons can also be extended to other regions for better management of their ecological land for sustainable use.
ZHOU Ting , QI Jialing , XU Zhihan , ZHOU De . Damage or Recovery? Assessing Ecological Land Change and Its Driving Factors: A Case of the Yangtze River Economic Belt, China[J]. Journal of Resources and Ecology, 2021 , 12(2) : 175 -191 . DOI: 10.5814/j.issn.1674-764x.2021.02.005
Fig. 1 The technical flowchart of this study |
Fig. 2 Geographical location of the YREB |
Table 1 Connotation of ecological land |
Types | Connotation | Literatures |
---|---|---|
Type 1 | Ecological land is defined as the space of ecological elements | Dong et al., 1999 |
Ecological land is one of the key resources and conditions for the survival of humans, and the total amount of the environment required or occupied by a species in a stable state | Xie et al., 2013; Zhu et al., 2015; Li et al., 2016; Li et al., 2020 | |
Ecological land’s main function is to provide ecological products and ecological services. It plays an important role in regulating, maintaining and ensuring regional ecological security | Deng et al., 2009; Zhang et al., 2015; Hu et al., 2020; Huang et al., 2017; Peng et al., 2017; Liu et al., 2018; Gao et al., 2020 ; Li et al., 2020; Zhang et al., 2020 | |
Ecological land includes green ecological land (grassland and unused land) and water ecological land | Chen et al., 2015 | |
Type 2 | Ecological land includes unproductive forest land, grassland, water and unused land. Production- ecological land includes cultivated land, garden land, productive forest land, grassland and water. Living-ecological land includes parks and green space, and land for scenic spots | Dang et al., 2014 |
Ecological land is relatively less used by humans. Ecological-production land has the dual functions of ecology and agricultural production, but the ecological function is stronger than the production function; while production-ecological land is mainly aimed at obtaining agricultural products, so the production function is stronger than ecological function | Zhang et al., 2015; Yu et al., 2017 | |
Ecological land includes complete ecological land, semi-ecological land, weak ecological land and non-ecological land | Liu et al., 2017a; Lin and Feng, 2018; Zou et al., 2018 |
Table 2 Types of ecological land change |
Levels | Transformation | Land use changes |
---|---|---|
Severe damage | Complete ecological land→Non-ecological land | Woodland, grassland, water, and other ecological land→Construction land |
Slight damage | Complete ecological land→Semi-ecological land | Woodland, grassland, water, and other ecological land→Arable land |
Semi-ecological land→Non-ecological land | Arable land→Construction land | |
Unchanged | Complete ecological land→Complete ecological land | Woodland, grassland, water, and other ecological land→Woodland, grassland, water, and other ecological land |
Semi-ecological land→Semi-ecological land | Arable land→Arable land | |
Slight recovery | Non-ecological land→Semi-ecological land | Construction land→Arable land |
Semi-ecological land→Complete ecological land | Arable land→Woodland, grassland, water, and other ecological land | |
Obvious recovery | Non-ecological land→Complete ecological land | Construction land→Woodland, grassland, water, and other ecological land |
Table 3 Index system for the driving factors of ecological land change |
Classes | First-level indicators | Basic-level indicators | References | |
---|---|---|---|---|
Natural factors (A) | Topography | A1 | Elevation (m) | Xie, 2011; Wang, 2012; Zhou, 2019 |
A2 | Slope (°) | Xie, 2011; Peng et al., 2017; Zhou, 2019 | ||
Climate | A3 | Annual average precipitation (mm) | Long, 2015; Zhou, 2019 | |
A4 | Annual average temperature (℃) | Long, 2015; Zhou, 2019 | ||
Social factors (B) | Urbanization level | B1 | Proportion of non-agricultural population to total population (%) | Zhang et al., 2007; Xie, 2011; Wang, 2012; Long, 2015; Wang, 2018 |
Population | B2 | Population density (person km-2) | Xie, 2011; Wang, 2012; Long, 2015; Tang et al., 2016; Wang and Chen, 2016; Wang, 2018 | |
Development scale | B3 | Proportion of construction land to the total land area (%) | Zhang et al., 2007; Tang et al., 2016 | |
Infrastructure | B4 | Urban road density (m2 person-1) | Zhang et al., 2007; Wang, 2018 | |
Savings level (Social Stability Index) | B5 | Year-end balance of savings of urban and rural residents (104 yuan) | Sun, 2005; Zhang and Wang, 2017 | |
Economic factors (C) | Economic development level | C1 | GDP (104 yuan) | Xie, 2011; Wang, 2012; Cui, 2015; Long, 2015; Zhang and Wang, 2017; Wang, 2018 |
C2 | Total investment in fixed assets (104 yuan) | Zhang et al., 2007; Cui, 2015 | ||
Consumption level | C3 | Total retail sales of social consumer goods (104 yuan) | Cui, 2015; Yu, 2016 | |
Industrial structure | C4 | Proportion of the secondary industry in the regional GDP (%) | Wang, 2012; Wang and Chen, 2016; Zhang and Wang, 2017; Wang, 2018 | |
C5 | Proportion of the tertiary industry in the regional GDP (%) | Wang, 2012; Wang and Chen, 2016; Zhang and Wang, 2017; Wang, 2018 |
Table 4 Types of interactions between two covariates |
Description | Interaction |
---|---|
q(X1∩X2) < min(q(X1), q(X2)) | Weaken, nonlinear |
min(q(X1), q(X2)) < q(X1∩X2) <max (q(X1), q(X2)) | Weaken, univariate |
q(X1∩X2) > max(q(X1), q(X2)) | Enhance, bivariate |
q(X1∩X2) = q(X1) + q(X2) | Independent |
q(X1∩X2) > q(X1) + q(X2) | Enhance, nonlinear |
Notes: The nonlinear-weaken effect means a smaller interactive effect of the driving factors X1 and X2 than each of their separate effects, indicating that the two driving factors weaken each other; the uni-weaken effect means that the less active factor reduces the effect of the other one, or a mild weakening effect; the bi-enhance effect means a greater interactive effect of the driving factors X1 and X2 than each of their separate effects; the independent effect means that the interactive effect amounts to the sum of the separate effects of X1and X2, implying the two driving factors are independent of each other; and the nonlinear-enhance effect means the strongest interactive effect of the driving factors X1 and X2 over the sum of their separate effects, a strong enhancement that does not show a simple (linear) proportional relationship. |
Fig. 3 Changes of ecological lands in the YREB from 1995 to 2015 |
Fig. 4 Spatial pattern of ecological lands in the YREB from 1995 to 2015 |
Table 5 Areas of ecological land changes in the YREB (km2) |
Stages | Severe damage | Slight damage | Unchanged | Slight recovery | Obvious recovery |
---|---|---|---|---|---|
1995‒2015 | 9053 | 161490 | 1713018 | 148307 | 3513 |
1995‒2000 | 3805 | 149134 | 1731878 | 146935 | 3511 |
2000‒2005 | 881 | 6525 | 2138318 | 3400 | 65 |
2005‒2010 | 677 | 4119 | 2043259 | 1289 | 26 |
2010‒2015 | 3433 | 8474 | 2035341 | 2006 | 220 |
Table 6 Areas of ecological land changes in different regions in the YREB from 1995 to 2015 (km2) |
Regions | Severe damage | Slight damage | Unchanged | Slight recovery | Obvious recovery |
---|---|---|---|---|---|
Upstream | 2477 | 87028 | 942775 | 84836 | 943 |
Midstream | 3303 | 35974 | 480717 | 33635 | 1142 |
Downstream | 3253 | 38475 | 259397 | 29839 | 1428 |
Fig. 5 Spatial distribution of types of ecological land damage in the YREB from 1995 to 2015 Note: (a), (b), (c) and (d) mainly show severe damage, slight damage, slight recovery and obvious recovery, respectively. |
Table 7 The area transition matrix of ecological land in the YREB from 1995 to 2015 (km2) |
Land types | Semi-ecological land | Complete ecological land | Non-ecological land | |||
---|---|---|---|---|---|---|
2015 1995 | Arable land | Woodland | Grassland | Water | Other ecological land | Construction land |
Arable land | - | 98872 | 27475 | 9664 | 312 | 29177 |
Woodland | 95973 | - | 57989 | 3973 | 780 | 5790 |
Grassland | 28806 | 55389 | - | 1587 | 671 | 1507 |
Water | 7308 | 3242 | 1032 | - | 414 | 1737 |
Other ecological land | 226 | 565 | 707 | 672 | - | 19 |
Construction land | 11984 | 2103 | 532 | 853 | 25 | - |
Fig. 6 Spatial characteristics of ecological land change for 130 cities from 1995 to 2015 Note: Red, green, and white indicate the percentages of the damaged area, the recovered area, and the unchanged area of a land type in the total administrative area, respectively. |
Table 8 Results of the factor detector |
Classes | First-level indicators | Basic-level indicators | Rank | q value | P value | |
---|---|---|---|---|---|---|
Natural factors (A) | Topography | A1 | Elevation | 14 | 0.0568 | 0.2379 |
A2 | Slope | 11 | 0.1338 | 0.0176* | ||
Climate | A3 | Annual average precipitation | 8 | 0.1598 | 0.0031** | |
A4 | Annual average temperature | 7 | 0.1599 | 0.0046** | ||
Social factors (B) | Urbanization level | B1 | Proportion of non-agricultural population to total population | 5 | 0.1650 | 0.0024** |
Population | B2 | Population density | 13 | 0.1302 | 0.0109* | |
Development scale | B3 | Proportion of construction land to the total area of the city | 6 | 0.1635 | 0.0110* | |
Infrastructure | B4 | Urban road density | 10 | 0.1572 | 0.0057** | |
Savings level (Social Stability Index) | B5 | Year-end balance of savings of urban and rural residents | 1 | 0.2700 | 0.0470* | |
Economic factors (C) | The level of economic development | C1 | GDP | 12 | 0.1316 | 0.0550 |
C2 | Total investment in fixed assets | 2 | 0.2078 | 0.1998 | ||
Consumption level | C3 | Total retail sales of social consumer goods | 9 | 0.1582 | 0.0594 | |
Industrial structure | C4 | Proportion of the secondary industry in the regional GDP | 4 | 0.1754 | 0.0025** | |
C5 | Proportion of the tertiary industry in the regional GDP | 3 | 0.1826 | 0.0021** |
Note: * P<0.05, indicates statistical significance at the 5% level; **P<0.01, indicates statistical significance at the 1% level. |
Table 9 Factor detector results of the five land type |
Basic-level indicators | Arable land | Woodland | Grassland | Water | Other ecological land |
---|---|---|---|---|---|
Elevation (A1) | 0.0675 | 0.0088 | 0.0243 | 0.0754 | 0.0176 |
Slope (A2) | 0.1326 | 0.0611 | 0.0383 | 0.0436 | 0.0229 |
Annual average precipitation (A3) | 0.2396*** | 0.0213 | 0.0068 | 0.0128 | 0.0189 |
Annual average temperature (A4) | 0.1696* | 0.0089 | 0.0259 | 0.0085 | 0.0099 |
Proportion of non-agricultural population to total population (B1) | 0.2700*** | 0.0132 | 0.0092 | 0.0537 | 0.01095 |
Population density (B2) | 0.2838* | 0.0007 | 0.0140 | 0.0408 | 0.0072 |
Proportion of construction land to the total area of the city (B3) | 0.1116* | 0.0377 | 0.0307 | 0.0414 | 0.0176 |
Urban road density (B4) | 0.1821** | 0.0268 | 0.0011 | 0.0152 | 0.0236 |
Year-end balance of savings of urban and rural residents (B5) | 0.0597 | 0.0518 | 0.0712 | 0.0924 | 0.0074 |
GDP (C1) | 0.0667 | 0.0640 | 0.0876 | 0.0661 | 0.0071 |
Total investment in fixed assets (C2) | 0.0719 | 0.0070 | 0.0318 | 0.0299 | 0.0196 |
Total retail sales of social consumer goods (C3) | 0.0569 | 0.0687 | 0.0517 | 0.0646 | 0.0107 |
Proportion of the secondary industry in the regional GDP (C4) | 0.2355*** | 0.0499 | 0.0032 | 0.0660 | 0.0424 |
Proportion of the tertiary industry in the regional GDP (C5) | 0.2106*** | 0.0232 | 0.0113 | 0.0508 | 0.0221 |
Note: *P<0.05, indicates statistical significance at the 5 % level; ** P<0.01, indicates statistical significance at the 1 % level; ***P<0.001, the independent variable is significantly related to the dependent variable. |
Fig. 7 The risk detector of the driving factors for the ecological land change Note: The X-axis is the subregion of driving factors; the Y-axis is the area of ecological land change (km2). A1-Elevation, A2-Slope, A3-Annual average precipitation, A4-Annual average temperature, B1-Proportion of non-agricultural population to total population, B2-Population density, B3-Proportion of construction land to the total area of the city, B4-Urban road density, B5-Year-end balance of savings of urban and rural residents, C1-GDP, C2-Total investment in fixed assets, C3-Total retail sales of social consumer goods, C4-Proportion of the secondary industry in the regional GDP, C5-Proportion of the tertiary industry in the regional GDP. |
Table 10 The explanatory power between any two driving factors |
Indicators | A1 | A2 | A3 | A4 | B1 | B2 | B3 | B4 | B5 | C1 | C2 | C3 | C4 | C5 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
A1 | 0.057 | |||||||||||||
A2 | 0.368 | 0.134 | ||||||||||||
A3 | 0.424 | 0.178 | 0.160 | |||||||||||
A4 | 0.413 | 0.186 | 0.179 | 0.160 | ||||||||||
B1 | 0.408 | 0.180 | 0.176 | 0.179 | 0.165 | |||||||||
B2 | 0.263 | 0.149 | 0.173 | 0.173 | 0.175 | 0.130 | ||||||||
B3 | 0.482 | 0.189 | 0.176 | 0.178 | 0.176 | 0.175 | 0.163 | |||||||
B4 | 0.392 | 0.169 | 0.196 | 0.204 | 0.191 | 0.169 | 0.214 | 0.157 | ||||||
B5 | 0.441 | 0.302 | 0.291 | 0.300 | 0.293 | 0.286 | 0.301 | 0.326 | 0.270 | |||||
C1 | 0.384 | 0.160 | 0.174 | 0.185 | 0.178 | 0.152 | 0.195 | 0.185 | 0.295 | 0.132 | ||||
C2 | 0.377 | 0.232 | 0.239 | 0.258 | 0.251 | 0.224 | 0.247 | 0.241 | 0.286 | 0.249 | 0.208 | |||
C3 | 0.477 | 0.196 | 0.174 | 0.179 | 0.178 | 0.173 | 0.176 | 0.204 | 0.290 | 0.229 | 0.266 | 0.158 | ||
C4 | 0.464 | 0.215 | 0.192 | 0.186 | 0.188 | 0.185 | 0.188 | 0.242 | 0.308 | 0.220 | 0.282 | 0.194 | 0.175 | |
C5 | 0.477 | 0.239 | 0.197 | 0.193 | 0.194 | 0.192 | 0.201 | 0.273 | 0.311 | 0.229 | 0.289 | 0.201 | 0.194 | 0.183 |
Note: The diagonal q values represent the explanatory power of the individual driving factors, whereas the off-diagonal q values represent the interactions between pairs of the driving factors. |
Table 11 Descriptive statistics for the regression coefficients of the geographically weighted regression model |
Indicators | Min | Max | Mean | Q1 | Median | Q3 | STD |
---|---|---|---|---|---|---|---|
Elevation (A1) | ‒0.07 | ‒0.04 | ‒0.05 | ‒0.06 | ‒0.06 | ‒0.05 | 0.01 |
Slope (A2) | 0.04 | 0.07 | 0.05 | 0.05 | 0.05 | 0.06 | 0.01 |
Annual average precipitation (A3) | ‒0.05 | 0.04 | ‒0.01 | ‒0.04 | ‒0.02 | 0.01 | 0.02 |
Annual average temperature (A4) | 0.01 | 0.09 | 0.05 | 0.03 | 0.05 | 0.07 | 0.02 |
Proportion of non-agricultural population to total population (B1) | 0.01 | 0.07 | 0.04 | 0.02 | 0.04 | 0.05 | 0.01 |
Population density (B2) | ‒0.10 | ‒0.05 | ‒0.07 | ‒0.09 | ‒0.08 | ‒0.06 | 0.01 |
Proportion of construction land to the total area of the city (B3) | ‒0.33 | ‒0.14 | ‒0.21 | ‒0.25 | ‒0.20 | ‒0.17 | 0.05 |
Urban road density (B4) | ‒0.06 | ‒0.05 | ‒0.05 | ‒0.05 | ‒0.05 | ‒0.05 | 0 |
Year-end balance of savings of urban and rural residents (B5) | ‒0.53 | ‒0.13 | ‒0.30 | ‒0.41 | ‒0.28 | ‒0.20 | 0.11 |
GDP (C1) | ‒0.65 | 0.04 | ‒0.39 | ‒0.55 | ‒0.43 | ‒0.22 | 0.18 |
Total investment in fixed assets (C2) | ‒0.59 | ‒0.51 | ‒0.54 | ‒0.56 | ‒0.54 | ‒0.52 | 0.02 |
Total retail sales of social consumer goods (C3) | 0.31 | 0.48 | 0.43 | 0.40 | 0.45 | 0.47 | 0.04 |
Proportion of the secondary industry in the regional GDP (C4) | ‒0.04 | 0.00 | ‒0.02 | ‒0.03 | ‒0.02 | ‒0.01 | 0.01 |
Proportion of the tertiary industry in the regional GDP (C5) | ‒0.01 | 0.01 | 0 | 0 | 0 | 0 | 0 |
Note: Q1 represents the Lower Quartile; Q3 represents the Upper Quartile; STD represents the Standard Deviation. |
Fig. 8 Spatial distribution of regression coefficients of driving factors Note: A1-Elevation, A2-Slope, A3-Annual average precipitation, A4-Annual average temperature, B1-Proportion of non-agricultural population to total population, B2-Population density, B3-Proportion of construction land to the total area of the city, B4-Urban road density, B5-Year-end balance of savings of urban and rural residents, C1-GDP, C2-Total investment in fixed assets, C3-Total retail sales of social consumer goods, C4-Proportion of the secondary industry in the regional GDP, C5-Proportion of the tertiary industry in the regional GDP. |
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|
17 |
|
18 |
|
19 |
|
20 |
|
21 |
|
22 |
|
23 |
|
24 |
|
25 |
|
26 |
|
27 |
|
28 |
|
29 |
|
30 |
|
31 |
|
32 |
|
33 |
|
34 |
|
35 |
|
36 |
|
37 |
|
38 |
|
39 |
|
40 |
|
41 |
|
42 |
|
43 |
|
44 |
|
45 |
|
46 |
|
47 |
|
48 |
|
49 |
|
50 |
|
51 |
|
52 |
|
53 |
|
54 |
|
55 |
|
56 |
|
57 |
|
58 |
|
59 |
|
60 |
|
61 |
|
62 |
|
63 |
|
64 |
|
65 |
|
66 |
|
67 |
|
68 |
|
69 |
|
70 |
|
71 |
|
/
〈 |
|
〉 |