Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Coastal Acoustic Tomography of the Neko-Seto Channel with a Focus on the Generation of Nonlinear Tidal Currents—Revisiting the First Experiment
Next Article in Special Issue
Evaluation of Land Surface Phenology for Autumn Leaf Color Change Based on Citizen Reports across Japan
Previous Article in Journal
Performance Evaluation of Feature Matching Techniques for Detecting Reinforced Soil Retaining Wall Displacement
Previous Article in Special Issue
Using Eco-Geographical Zoning Data and Crowdsourcing to Improve the Detection of Spurious Land Cover Changes
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 14
Issue 7
10.3390/rs14071698
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating Future LUCC by Coupling Climate Change and Human Effects Based on Multi-Phase Remote Sensing Data

by
Zihao Huang
1,2,3,
Xuejian Li
1,2,3,
Huaqiang Du
1,2,3,*,
Fangjie Mao
1,2,3,
Ning Han
1,2,3,
Weiliang Fan
1,2,3,
Yanxin Xu
1,2,3 and
Xin Luo
1,2,3
1
State Key Laboratory of Subtropical Silviculture, Zhejiang A & F University, Hangzhou 311300, China
2
Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon Sequestration of Zhejiang Province, Zhejiang A & F University, Hangzhou 311300, China
3
School of Environmental and Resources Science, Zhejiang A & F University, Hangzhou 311300, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(7), 1698; https://doi.org/10.3390/rs14071698
Submission received: 18 February 2022 / Revised: 22 March 2022 / Accepted: 29 March 2022 / Published: 31 March 2022
(This article belongs to the Special Issue Advanced Phenology, and Land Cover and Land Use Change Studies)
Browse Figures
Figure 1
<p>Zhejiang Province study area: (<b>a</b>) location in China, (<b>b</b>) sample plots in 2014, (<b>c</b>) forest coverage rate from 2004 to 2020, and (<b>d</b>) land-use patterns in 2014.</p> "> Figure 2
<p>(<b>a</b>) Land-use patterns and (<b>b</b>) accuracy evaluation of land-use classification in Zhejiang Province from 1984 to 2014.</p> "> Figure 3
<p>Topographic data in 2014: (<b>a</b>) DEM, (<b>b</b>) slope, and (<b>c</b>) aspect; soil data in 2008: (<b>d</b>) clay fraction, (<b>e</b>) sand fraction, (<b>f</b>) silt fraction, (<b>g</b>) soil available water content, (<b>h</b>) soil bulk density, and (<b>i</b>) soil wilt point; climate data in 2014: (<b>j</b>) total precipitation, (<b>k</b>) average temperature, (<b>l</b>) average radiation, and (<b>m</b>) average relative humidity; socioeconomic data in 2015: (<b>n</b>) population density, (<b>o</b>) GDP; distance data in 2014: (<b>p</b>) distance to road, (<b>q</b>) distance to railway, and (<b>r</b>) distance to water.</p> "> Figure 4
<p>Macro statistics: (<b>a</b>) total population; (<b>b</b>) GDP; (<b>c</b>) grain yield; (<b>d</b>) aquatic product yield.</p> "> Figure 5
<p>Configurations of four development scenarios concerning human and natural effects: (<b>a</b>) annual total precipitation; (<b>b</b>) annual average temperature; (<b>c</b>) annual average radiation; (<b>d</b>) annual average relative humidity under four scenarios.</p> "> Figure 6
<p>Flow chart used in this study.</p> "> Figure 7
<p>The structure of the SD model. Variables in red indicate the inputs of the model, which varied in different scenarios. Variables in purple represent the main outputs, including the areas of different LULC types. Purple ones were used to validate the model.</p> "> Figure 8
<p>The interactive coupling mechanism of the SD model and BPNN_CA model.</p> "> Figure 9
<p>The actual land-use pattern in 2008 and 2014 and the simulated land-use pattern in 2014. (<b>A</b>,<b>B</b>) are the two sub-areas shown in magnification.</p> "> Figure 10
<p>Verification: (<b>a</b>) comparison between the classified and simulated areas using the SD model; (<b>b</b>) ROC curves and AUC values fitted by the BPNN; (<b>c</b>) normalized confusion matrix of the simulated result and the actual classification in 2014; (<b>d</b>) the FOM from 2008 to 2014.</p> "> Figure 11
<p>Areas of each land-use type in Zhejiang Province over the past 30 years and that over the next 70 years under different scenarios: (<b>a</b>) urban land area; (<b>b</b>) water body area; (<b>c</b>) cultivated land area; (<b>d</b>) bamboo forest area; (<b>e</b>) broad-leaved forest area; (<b>f</b>) coniferous forest area.</p> "> Figure 12
<p>Spatial distribution of land use in different regions of Zhejiang Province in 2084 under different scenarios. (<b>A</b>,<b>B</b>) are the two sub-areas shown in magnification.</p> "> Figure 13
<p>Predicted land-use conversion in 2014–2084 under 4 scenarios: (<b>a</b>) SD_Scenario; (<b>b</b>) HD_Scenario; (<b>c</b>) BD_Scenario; (<b>d</b>) FD_Scenario. Numbers: years; UL: urban land; WB: water body; CL: cultivated land; BLF: broad-leaved forest; CF: coniferous forest; BF: bamboo forest. For example, 14BF represents bamboo forest in 2014.</p> "> Figure 14
<p>Land use/cover change amplitude at the municipal administrative level in Zhejiang Province from 2014 to 2084 under different future scenarios.</p> ">
Review Reports Versions Notes

Abstract

:
Future land use and cover change (LUCC) simulations play an important role in providing fundamental data to reveal the carbon cycle response of forest ecosystems to LUCC. Subtropical forests have great potential for carbon sequestration, yet their future dynamics under natural and human influences are unclear. Zhejiang Province in China is an important distribution area for subtropical forests. For forest management, it is of great significance to explore the future dynamic changes of subtropical forests in Zhejiang. As a popular LUCC spatial simulation model, the cellular automata (CA) model coupled with machine learning and LUCC quantitative demand models such as system dynamics (SD) can achieve effective LUCC simulation. Therefore, we first integrated a back propagation neural network (BPNN), a CA, and a SD model as a BPNN_CA_SD (BCS) coupled model for future LUCC simulation and then designed a slow development scenario (SD_Scenario), a harmonious development scenario (HD_Scenario), a baseline development scenario (BD_Scenario), and a fast development scenario (FD_Scenario), combining climate change and human disturbance. Thirdly, we obtained future land-use patterns in Zhejiang Province from 2014 to 2084 under multiple scenarios, and finally, we analyzed the temporal and spatial changes of land use and discussed the subtropical forest dynamics of the future. The results showed the following: (1) The overall accuracy was approximately 0.8, the kappa coefficient was 0.75, and the figure of merit (FOM) value was over 28% when using the BCS model to predict LUCC, indicating that the model could predict the consistent change of LUCC accurately. (2) The future evolution of the LUCC under different scenarios varied, with the growth of bamboo forests and the decline of coniferous forests in the FD_Scenario being prominent among the forest dynamics changes. Compared with 2014, the bamboo forest in 2084 will increase by 37%, while the coniferous forest will decrease by 25%. (3) Comparing the area and spatial change of the subtropical forests, the SD_Scenario was found to be beneficial for the forest ecology. These results can provide an important decision-making reference for land-use planning and sustainable forest development in Zhejiang Province.

1. Introduction

Land use and land cover change (LUCC) is a direct driving factor of the carbon balance in terrestrial ecosystems, and its impact on global warming is second only to that of fossil fuels and industrial emissions [1,2,3,4]. Therefore, the impact of LUCC and climate change on forest spatiotemporal dynamics has been widely appreciated. However, limited LUCC data may lead to significant underestimations within the impact of LUCC on carbon emissions [5], and the absence of future LUCC data under future climate change is a major limitation for revealing the response of forest ecosystems’ carbon cycles to future climate change [6]. Therefore, it is of great scientific significance to obtain future LUCC data through spatiotemporal LUCC simulations and to explore the LUCC evolution in order to reveal the impact of LUCC on the forest ecosystem’s carbon cycle.
There are the following three types of models to simulate LUCC: quantitative models, spatial models, and coupled models of the first two. Quantitative demand models focus on predicting the area transfer among different land-use types, including system dynamics (SD), grey prediction, and Markov chain models [7,8,9]. Most of these models are statistical models. The SD model is a simulation method in which inventories and flows, with corresponding feedback loops, are used to simulate large-scale, complex socioeconomic systems [10]. SD models seek to understand how physical processes, information flows, and management policies interact [11] so as to link feedback between different variables in the system and, therefore, it can be used to deal with dynamic processes [12]. Previous studies have been conducted to show that the SD model for LUCC demand prediction is the most effective of the quantitative models, mainly because it considers comprehensive factors such as the impact of markets, policies, and climate adaptation strategies [13]. However, these quantitative models lack the ability to project the spatial patterns of land use; thus, spatial models that reveal LUCC at a spatial scale are required. The cellular automata (CA) model and the conversion of land use and its effects at a small regional extent (CLUE-S) model are popular methods to simulate the LUCC spatial evolution [12,14,15,16,17]. These models with spatial function are all scale-dependent. Among these models, the CA model has simple preconditions and rules expressed by a matrix where a state defines each cell [10,15]. It can simulate spatial dynamics from a bottom-up perspective and is capable of establishing interconnections between LUCC and driving forces [18,19]. Its main advantages are in expressing the driving forces of LUCC by transition rules and expressing spatial externalities by neighborhood effects [20,21]. In addition, it is easily combined with regression models to identify locations with high suitability for LUCC based on weighted overlay suitability factors [17,22,23]. This has become an excellent method for predicting future LUCC according to different scenarios [24,25].
However, the single model described above could not fully consider the internal mechanisms of the ecosystem when simulating LUCC. Therefore, the coupled model integrated by a top-down quantitative model and a bottom-up spatial model has been popular to improve the accuracy of LUCC simulation [7]. The FLUS model could be applied to the effective Chinese LUCC simulation considering various socioeconomic and natural climatic factors, which are integrated by an artificial neural network, a CA model, and an SD model [26]. Considering that the performance of the back propagation neural network (BPNN) is better than some artificial neural networks [27], the CA model integrated with the BPNN (BPNN_CA) can be designed to control the spatial pattern changes more precisely. Furthermore, the interactive coupling of the CA model and the SD model is more effective and popular than the loose coupling of the two. Therefore, it is achievable to interactively couple the three models of BPNN, CA, and SD together as a BPNN_CA_SD (BCS) model for simulating LUCC precisely. Gradually, the refinement of coupled models meets the accuracy requirements of LUCC spatiotemporal simulation.
Subtropical forests, with their diverse types, high photosynthetic capacity, and four-season growth, have been of global concern for their carbon sequestration capacity, which accounts for approximately 40% of the world’s gross primary productivity together with tropical forests [28,29]. China is an important distribution area for subtropical forests. The vegetation carbon storage of the evergreen broad-leaved forest and coniferous and broad-leaved mixed forest in subtropical China is 2.527 billion tons, accounting for approximately 30% of the country [30]. The carbon sink here has great potential, but it is very sensitive to LUCC caused by global climate change and anthropogenic disturbances. Regarding the impact of LUCC on the subtropical carbon cycle, its past and current impacts have been studied [31,32], while the impact of LUCC is still unknown under future scenarios [33,34]. Simultaneously, future spatiotemporal LUCC modeling studies mainly focus on urban expansion [35,36], while less research has been conducted on the future distribution of subtropical forests, and the complex competitive relationships between subtropical forests and other land-use types are not yet clear. Given this situation, we posed our research questions as follows: (1) How to establish a coupled model that takes into account natural and human influences that can accurately simulate the quantitative and spatial changes of subtropical forests? (2) How will subtropical forests change in the future under natural and human effects, and which scenarios are beneficial to subtropical forests?
As a typical region of subtropical forests, Zhejiang Province has one of the highest forest coverage rates in China, and its forest changes have been of great concern. Simultaneously, the Coupled Model Intercomparison Project 5 (CMIP5) stage could provide reasonable scenario data for projecting future LUCC dynamics by combining representative concentration pathways (RCPs) 2.6, 4.5, 6.0, and 8.5 scenarios [37,38]. RCPs are used for scenario modeling under natural influences such as greenhouse gas emissions and radiative forcing. The combination of RCPs and socioeconomic development scenarios into integrated scenarios will make future LUCC projections more reasonable and accurate [39]. Thus, in this study, we take Zhejiang Province as an example to develop a BPNN_CA_SD (BCS) coupled model for future LUCC simulation under multiple development scenarios. The objectives consisted of the following four specific parts: (1) to build an SD model to project the land-use area demands, (2) to establish a BPNN_CA model to deal with the simulation of the land-use spatial patterns, (3) to interactively couple BPNN_CA with SD as the BCS model to obtain future spatiotemporal land-use patterns under the multiple natural and socioeconomic development scenarios we set, and (4) to analyze the LUCC of Zhejiang Province under different scenarios and provide a scenario reference for managing subtropical forests in response to the impacts of climate change and human disturbances.

2. Materials and Methods

2.1. Study Area

Zhejiang Province (Figure 1a) is located in the south of Yangtze River Delta on the southeast coast of China, with a total area of 1.055 × 105 km2. Zhejiang Province includes 11 administrative cities (Figure 1b), including Hangzhou, Ningbo, and Wenzhou. The province belongs to a subtropical monsoon climate. It has rich forest resources, and the forest coverage rate of the whole province has been increasing since 2004 and reached 61.15% in 2020 (Figure 1c). The major forest types are broad-leaved forest, coniferous forest, and bamboo forest (Figure 1d).

2.2. Data Sets and Processing

Natural environmental factors and socioeconomic factors have been proposed as the two main driving forces of land use/cover change by the IGBP and IHDP [40]. To simulate how climate change and human activities affect LUCC, these two types of driving forces (Table 1) were considered and consisted of the following three components: (1) geospatial data were used for spatial simulation. All the geospatial data were projected onto a WGS_1984 coordinate system and resampled to the same resolution (30 × 30 m) of land-use patterns. (2) Macro statistical data were used for quantitative simulation. (3) Sample plotted data were the classification verification data of each land-use type, which were used to verify the accuracy.
Table
Table 1. List of data used in this study.
Table 1. List of data used in this study.
Data SetCategoryDataYearResolutionDiagramData Resource
Geo-spatial dataLand useLand-use patterns1984–201430 mFigure 2a,bThe data were based on the Satellite 30 m multispectral data of Landsat-5 TM (1984–2008) and Landsat-8 OLI (2014). After radiation correction, atmospheric correction, and geometric correction, the maximum likelihood classification method was used to obtain land-use patterns.
TerrainDEM201430 mFigure 3a–cDownloaded from the Geospatial Data Cloud site (http://www.gscloud.cn, accessed on 24 December 2021).
SlopeCalculated from DEM.
Aspect
SoilSilt fraction20081 kmFigure 3d–iSilt fraction, clay fraction, sand fraction, and available water content were derived from the Harmonized World Soil Database (HWSD 1.2).
The bulk density and soil wilt point were calculated by the silt and clay fraction [41].
Clay fraction
Sand fraction
Available water content
Bulk density
Wilt point
ClimateTotal precipitation1984–20141 kmFigure 3j–mThe annual data were calculated from the averages or sums of the daily data. The daily data were interpolated from observations at 410 meteorological stations in Zhejiang Province and its surrounding provinces using the inverse distance weighted method [42].
Average temperature
Average radiation
Average relative humidity
Human influencePopulation20151 kmFigure 3n,oObtained from the Resource and Environmental Science and Data Center of the Chinese Academy of Sciences (http://www.resdc.cn, accessed on 24 December 2021).
Gross domestic product (GDP)
Distance to roads201430 mFigure 3p–rCalculated from the vector maps of the roads, the railways, and the water systems, which were downloaded from the Open Street map (https://www.openstreetmap.org/, accessed on 12 October 2020).
Distance to railways
Distance to water
Macro statistics dataLand-use area1984–2014--Calculated from the land-use patterns.
Total precipitation statistics-Calculated from the total precipitation.
Average temperature statistics-Calculated from the average temperature.
Population statisticsFigure 4a–dCollected from the Zhejiang Statistical Yearbook (http://tjj.zj.gov.cn/, accessed on 12 October 2020).
GDP statistics
Grain yield
Aquatic product yield
Forest coverage rate2004–2020Figure 1cCollected from the Announcement of Forest Resources and Its Ecological Function Value of Zhejiang Province (http://lyj.zj.gov.cn/index.html, accessed on 24 December 2021).
Sample plots dataClassification verification plots1984–2014-Figure 1b and
Table 2		Classification verification plots of BLF, CF, and BF were derived from the data of the National Forest Inventory. Verification plots of other land-use types were based on field investigation and image visual interpretation.
The land use of Zhejiang Province (Figure 2a) was classified into the following six types: urban land (UL), water body (WB), cultivated land (CL), broad-leaved forest (BLF), coniferous forest (CF), and bamboo forest (BF) [43,44]. The sample plots data (Table 1 and Table 2) were used to calculate the overall accuracy (OA) and kappa coefficient (Kappa) to evaluate the accuracy of the classification results. Figure 2b shows the normalized confusion matrices from 1984 to 2014. Each value of the main diagonal in the normalized confusion matrix corresponds to the producer’s accuracy (PA) value for each land-use type. It can be seen that the overall accuracy (OA) and kappa coefficients (Kappa) in different years were all higher than 0.78 and 0.73, respectively, indicating that the classification results were highly consistent with the actual situation. Moreover, most of the PA values were higher than 76%, which provides an important guarantee for spatiotemporal LUCC simulation in the study area.

2.3. Future Scenario Description

Four development scenarios (Figure 5) were designed based on the CMIP5 while considering climate variations together with different socioeconomic developments in Zhejiang Province.
Future climate data were derived from BCC-CSM1-1 climate change modeling data under the RCP 2.6, 4.5, 6.0, and 8.5 scenarios proposed by the IPCC AR5 [45]. The data under each scenario included annual total precipitation, annual average temperature, annual average radiation, and annual average relative humidity, which correspond to past climate data so that the BCS model can alternate input data in the prediction phase. Future socioeconomic scenario settings (Table 3) were mainly based on past macro-statistical data (Figure 4).
The slow development scenario (SD_Scenario) was constructed to predict the demand for land use under the influence of slow socioeconomic growth and mild climate change. Corresponding to the RCP2.6 scenario, it is the lowest scenario for greenhouse gas emissions and radiative forcing. Radiative forcing will increase first and decrease after 2054, reaching 2.6 W · m−2 by 2100, with global average warming limited to 2 °C [47]. In this scenario, the population, GDP, and technological innovation are considered to be at their lowest levels, while the annual temperature and precipitation change slowly. It promotes the use of biomass energy and advocates the restoration of forests, so the forests will be well restored and the forest coverage rate will be the highest.
The harmonious development scenario (HD_Scenario) is a sustainable development model. Moderate growth of the population and GDP are taken into account, and the proportion of investment is assumed to have more input into the productivity of agriculture and fisheries. Additionally, the natural environment will undergo moderate changes (RCP4.5 scenario), forming a model of sustainable development with the economy and society. This scenario limits greenhouse gas emissions through low-emission energy technologies so that increasing radiative forcing will only reach 4.5 W · m−2 by 2100.
The baseline development scenario (BD_Scenario) is established based on the past and current development trends of Zhejiang Province. Under this scenario, the population, economy, and technological level are recognized as stable and advanced. Moreover, the climate is assumed to maintain its current temperature and precipitation rates, which is consistent with the RCP6.0 scenario [48].
Contrary to the SD_Scenario, the fast development scenario (FD_Scenario) aims to maximize the social and economic benefits of Zhejiang Province. The economy and population increase at a high speed, and science and technology develop rapidly. At the same time, massive-scale human activities will accelerate greenhouse gas emissions and increase atmospheric radiation, leading to a sharp increase in temperature (RCP 8.5) and, thus, resulting in drastic climate change [49,50].

2.4. Methodology

In the paper, we present the BCS model for multiple LUCC scenarios for future land use by coupling human and natural effects. The proposed model consists of an SD model and a BPNN_CA model. The SD model was used to project the land-use demand at the macro-level, and the BPNN_CA model was used to allocate land use spatially at the micro-level.
The flow chart used is shown in Figure 6. Firstly, based on the geospatial data in the past, sample sets were randomly sampled to train the BPNN, and then the overall probability was calculated, and the simulated land-use types were determined according to the overall probability and roulette-wheel selection mechanism. Secondly, macro-statistical data were used to establish the relationship between land-use area and macro-factors so as to capture the feedback loop between them. Thirdly, the SD model was interactively coupled with the CA model, and the simulated land-use pattern was output only when the spatially allocated area of the CA model reached the requirements of the SD model. Finally, when the model accuracy was valid, the inputs of the two sub-models were updated using future data to obtain the land-use pattern in the future.

2.4.1. SD Model

The SD model for LUCC was mainly used to describe the quantitative conversion between different land-use types to achieve quantitative demand forecasting [10,51]. The relationship of quantitative transformation requires consideration of the market, investment, policies, and climate adaptation strategies from a macro perspective [8]. By understanding the impact of the internal factors of the macro system on LUCC, a top-down quantitative model can be formed [12].
The SD model was built using the Vensim PLE software (Figure 7). The simulation period of the SD model was 1984–2084 and the time step of the model was one year. It included the following two stages: (1) 1984–2014 was the model test stage, and the historical data were used to set the parameters, adjust the model, and validate the model; (2) 2014–2084 was the prediction stage, and the future land-use demand under different scenarios could be simulated.
To simulate the quantitative demand by considering both human and natural effects, the SD model in our study consisted of the following five subsystems: population, economy, productivity, climate, and land use. As the main social factor, the population influences the land-use system in many aspects. The economy has strong influences on population and land use such as gross domestic product (GDP), which affects the change in fixed-asset investments, thereby driving economic investment in various land-use types. The productivity mainly includes food and aquatic products; the amount of productivity indirectly determines the increase and decrease in cultivated land and water bodies. The climate subsystem involves temperature and precipitation, which have various influences on the growth and regeneration capacity of vegetation, leading to changes in cultivated land and forest land. The land-use subsystem includes the increment or decrement of each land-use demand. The change in each land-use type is constrained by the integrated influences of socioeconomic and climate conditions as well as by the interactions among the various land-use types.
The values of the level variables (e.g., total population, GDP) for each year are obtained based on the annual growth rate for that year compared to the previous year. Under different development scenarios, the annual growth rates are set differently, and the land-use demand outputs obtained by relevant formulas will also be greatly different.

2.4.2. BPNN_CA Model

The CA model is the most common and widely used model for LUCC spatial simulations. Based on the Python language programming, the CA model can be integrated with BPNN as the BPNN_CA model. Whether a cell state changes is not only related to external factors (land suitability) but also depends on the interrelationships between various land uses (e.g., neighborhood effects, land inheritance) and stochastic occurrences [26]. Therefore, the transition rules of the model can be expressed as follows:
S t a t e i , t + 1 = f P r o i , k t
where S t a t e i , t + 1 is the land-use type of a cell i at time t + 1 . f is the land-use transition rule that decides the change of a cell i from the state ( S t a t e i , t ) at time t to the state ( S t a t e i , t + 1 ) at time t + 1 . P r o i , k t is the overall probability of the cell i being occupied by a specific land-use type k at time t , which is defined as follows:
P r o i , k t = L s i , k × Ω i , k t × I k t × 1 − C s → k
(1) L s i , k is the land suitability probability of a cell i being occupied by land-use type k , which represents the impact of the driving factors on land-use transition. The collected driving factors (Figure 3) normalized to the range of [0,1] were selected to establish the BPNN model. The BPNN had the following four layers: an input layer, two hidden layers, and an output layer, with the corresponding numbers of neurons set to (19, 30, 14, 6). Simultaneously, the sigmoid function was selected as the activation function to ensure that L s i , k fell within [0,1]. A total of 6 million randomly selected samples were divided into a training set and a test set in a ratio of 7:3. The training set was used to train the neural network, while the test set was used to evaluate its performance. The minibatch gradient descent method (Algorithm 1) was adopted to input the training set into the BPNN module in batches to train and optimize the parameters of the network. In the algorithm, the number of iterations ( K ) was 200, and the initial learning rate was 0.0005. In each iteration, the loss function of the mean square error ( M S E ) was used to compute the gradients ( g r a d ), and the adaptive optimization of Adam was used to automatically adjust the learning rate ( l r ) after training a batch of images. The training of the BPNN stopped after K passed through the data set.
Algorithm 1: Train BPNN with the minibatch Adam optimization algorithm.
initialize ( n e t )
for e p o c h = 1, …, K do
  for b a t c h = 1, …, #   i m a g e s / b do
     i m a g e s   ← uniformly sample b a t c h − s i z e images
     X , y   ← preprocess(images)
     L s i , k   ← forward (net, X )
     M S E   ← loss ( L s i , k , y )
     l r , g r a d   ← backpropagation ( M S E )
    update ( n e t , l r , g r a d )
  end for
end for
(2) Ω i , k t is the neighborhood effect of land-use type k on grid cell i at iteration time t affected by the surrounding neighborhood at iteration time t − 1 , which is given as follows:
Ω i , k t = ∑ N × N S t a t e i , t − 1 = k N × N − 1
where ∑ N × N S t a t e i , t − 1 = k is the total number of grid cells occupied by land-use type k at iteration time t − 1 within the N × N window. In this study, the 5 × 5 Moore neighborhood was adopted.
(3) I k t is the inertia coefficient of land-use type k at iteration time t . If the development of a specific land-use type contradicts the future quantitative demands, the inertia coefficient would dynamically control the inheritance of the land-use type to increase or decrease to rectify the changing trend in the next iteration, which is expressed as follows:
I k t = I k t − 1 , i f   D k t − 1 ≦ D k t − 2 I k t − 1 × D k t − 2 D k t − 1 , i f D k t − 1 < D k t − 2 < 0 I k t − 1 × D k t − 1 D k t − 2 , i f 0 < D k t − 2 < D k t − 1
where D k t − 1 、 D k t − 2 are the differences between the quantitative demand and the allocated amount of land-use type k until the iteration times of t − 1 and t − 2 , respectively. Based on the ratio of the two values, the model updates I k t in real time.
(4) C s → k is the conversion cost, which is the difficulty of converting a specific grid pixel from land-use type s to the target land-use type k . It is a constant parameter, the value of which is fixed within [0,1] (Table 4). Additionally, larger values indicate more difficult conversions.
(5) The f for this study was a roulette-wheel selection mechanism [52], which assumes that the probability of selection is proportional to the fitness of a sector. Let us consider N sectors, each characterized by its fitness p i , k t > 0 ( n   = 1, …, N ). As shown in Algorithm 2, one constructs a line segment of length 1 out of consecutive sectors of length p i , k t , generates a uniformly distributed random number r i t such that r i t   ≤ q i , n t , and locates the corresponding sector, thus selecting the respective label T y p e n to assign to the state ( S t a t e i , t + 1 ) at time t + 1 . In this paper, n   = k = 1, …, 6 and T y p e n = (UL, WB, CL, BLF, CF, BF). Spatial allocation using this mechanism not only ensures the possibility of all events occurring but also ensures the dominance of high probability events, therefore ensuring the randomness and fairness in the allocation and effectively reducing the uncertainty of LUCC.
Algorithm 2: Using a roulette-wheel selection mechanism to allocate the probability.
input: P r o i , k t
p i , k t   ←   P r o i , k t / ∑ n = 1 N P r o i , k t
q i , n t   ←   ∑ k = 1 n p i , k t
r i t   ← a uniformly distributed random number ranging from 0 to 1
for n   = 1, …, N do
  if r i t   ≤ q i , n t then
     S t a t e i , t + 1 ←   T y p e n
    break
  else
    continue
end for

2.4.3. Interactive Integration of the BCS Model

Many integrated models are loosely coupled with the quantitative model and the spatial model based on the final land-use demands so that the accuracy is not as high as with the interactive coupling model [9]. Thus, to strengthen the mutual feedback between the SD and BPNN_CA sub-models, the two sub-models in the BCS model are interactively coupled during the study time series.
The schematic diagram of the coupling mechanism is shown in Figure 8. The future 70 year period from 2014 to 2084 was divided into the following 14 segments at 5 year intervals: 2014–2019, 2019–2024, …, and 2079–2084. Using the actual classification in 2014 as input data, the process was as follows: (1) the land-use pattern was input to the BPNN_CA to calculate land suitability and neighborhood effects, etc., while the area of each land-use type was fed into the SD model; (2) the demand of each land-use type at the next time was predicted through the configured SD model using the area at the previous time combined with the effects of human and natural factors during this time interval; (3) the demand was input to the BPNN_CA model as a constraint at the end of the model iteration, and the CA model explored the local competition and interactions by adjusting the I k t through multiple iterations; (4) the spatial distribution of the land-use pattern for the next moment was generated until the area allocated by the BPNN_CA model met the demand of the SD model. Looping steps (1)–(4) above, finally, the SD and BPNN_CA models would exchange input/output information to generate the land-use pattern in 2019, 2024, …, 2084 sequentially.

2.4.4. Assessment Methods of the BCS Model

The performance evaluation of the BCS coupled model needs to verify its ability to simulate the area and spatial pattern of land use.
For the area simulation, the determination coefficient R 2 and root mean square error (RMSE) were calculated using predicted area and actual area to account for the fitness of the SD model. The higher the R 2 and the smaller the RMSE, the better the result. The R 2 and RMSE are defined as follows:
R 2 = 1 − ∑ i = 1 N y p − y a i 2 ∑ i = 1 N y p − y ¯ a i 2
R M S E = 1 N ∑ i = 1 N y p − y ¯ a i 2
where y p represents the simulated area, y a represents the actual area, y ¯ a represents the average actual area, and N represents the number of land-use types, N = 6 .
For the spatial simulation, not only were OA, Kappa, and PA used for the evaluation but also a figure of merit (FOM) was introduced. The reason was that OA can only show the overall consistency of the simulation results and the actual results, but it cannot show the consistency of cell state changes, and the F o M index can directly show the ability of the model to simulate changes [53]. The OA, Kappa, PA, and FOM are expressed as follows:
O A = a 0 = 1 n ∑ k = 1 N a k k  
K a p p a = a 0 − a c 1 − a c ,   a c = ∑ k = 1 N ∑ i = 1 N a k i * t k n * n  
P A k = a k k / t k
F O M = R i g h t s / R i g h t s + M i s s e s + F a u l t s + E x t r a s
where n represents the total number of pixels; k and i   represents the land-use type, k , i = 1, 2, 3 …   N ; a k k   refers to the number of the category k correctly simulated; a k i is the number of type i simulated as k ; t k refers to the actual total number of category k ; M i s s e s is an area of error due to the actual changes simulated as no change; R i g h t s is an area of correctness due to the actual changes consistent with the simulated changes; F a u l t s is an area of error due to the actual change simulated as changing to an incorrect category; E x t r a s is an area of error due to the fact of no actual change simulated as change.

3. Results

3.1. Model Validations

The simulation for 2014 obtained through the BCS model is shown in Figure 9. Two partial enlargements not only show that the result was consistent with the actual pattern but also indicate that the spatial LUCC from 2008–2014 was also well simulated. However, there were still some errors in local details, and not all changes could be simulated; therefore, we quantitatively evaluated them to determine the performance of the BCS model.
Based on the simulation results and the actual classification in 2014, we evaluated the BCS model quantitatively, which consisted of the following four components: (1) the R 2 and RMSE were used to evaluate the SD model for predicting land-use area. The high R 2 (0.99) with the small RMSE (707 km 2 ) demonstrates that the calibrated SD model provides an important basis for the accurate prediction of future land-use area demand (Figure 10a); (2) The receiver operating characteristic (ROC) curve and the area under the ROC curve (AUC) values were used to quantify the BPNN’s performance. A larger AUC value corresponds to a better BPNN fitting performance [17]. The AUC values were all above 0.8 (Figure 10b), indicating that the land suitability fit for each land use can be strongly explained by the selected driving factors; (3) The normalized confusion matrix (Figure 10c) was used to evaluate the spatial consistency of the two maps in 2014. The high OA (79.61%) and Kappa (75.53%) show that the BCS model simulation results in 2014 were in excellent spatial consistency with the actual classification; (4) The FOM calculated from the simulation results in 2014 and the actual classifications in 2008 and 2014 was 28.21%, representing a large proportion of the area where the actual change was consistent with the simulated change (Figure 10d). Pontius et al. [54] reported FOM ranging from 1 to 59%, most of which were lower than 30%. Considering that our study area was a large region with complex climatic conditions and significant regional differences, the simulation accuracy was quite acceptable.
In summary, the BCS coupled model can not only be used to simulate great land-use spatial distributions but can also be used to show consistent LUCC between the simulated results and the actual classification, guaranteeing an accurate simulation of the future spatial distribution of land use.

3.2. Future Land Use Demand Projection

According to the four scenarios described in Section 2.3, the SD model with the interaction mechanism was used to predict the land-use demands for the next 70 years (Figure 11).
Figure 11a shows that the UL area will increase constantly from 2014 to 2084 under all scenarios due to the population increase. In the SD_Scenario and FD_Scenario, it will grow to 1.28 × 104 km2 and 1.63 × 104 km2 in 2084. Compared to the area in 2014, it will increase by 10% and 15% (+10% and +15%), respectively. The area in the BD_Scenario and HD_Scenario is in between that in the first two scenarios.
Figure 11b shows that there will be slight WB area changes in the future. It is between 2250 and 2750 km2 in the next 70 years under each scenario.
Figure 11c shows that the CL area will gradually decrease in the future. In the SD_Scenario and FD_Scenario, it will be reduced to 2.08 × 104 and 2.37 × 104 km2 in 2084 (–28% and –18%), respectively. The area in the BD_Scenario and HD_Scenario was in between that in the first two scenarios.
Figure 11d shows the BF area has slowly increased since 1984. Its differences between different scenarios were small in 2014–2054, while the differences started to become apparent after 2054. By 2084, in the FD_Scenario, it will reach 1.25 × 104 km2 (+37%), which is 1217 km2 more than that in the SD_Scenario.
Figure 11e shows the BLF areas have increased obviously over the past 30 years, while in the future 70 years, they will fluctuate between 2.25 × 104 km2 and 3 × 104 km2, with minor differences and changes among different scenarios.
Figure 11f shows the CF area varied clearly under different scenarios. In the FD_Scenario, it will have a decreasing trend after 2019, decreasing to 2.04 × 104 km2 in 2084, indicating that the urbanization and other development in Zhejiang Province under this scenario will have a strong negative effect on the development of the CF area. While in the SD_Scenario, it will reach 2.72 × 104 km2 (+13%), revealing that the ecological focus and better protection of CF under this scenario has led to an increase rather than a decrease in CF area.

3.3. Future Spatiotemporal Land-Use Pattern

Based on the BCS coupled model, future land-use patterns in 2084 were simulated. To demonstrate the differences between different scenarios, two enlargements of the simulation results are shown in Figure 12.

3.4. Analysis of Future Land-Use Conversion

With a total of six land-use types in this study, there were 36 possibilities for future land conversion. The Sankey diagrams (Figure 13) were used to show the land-use conversion from 2014 to 2084 under four scenarios, and the top five highest change amplitudes of all land uses are marked on the diagrams.
In both the SD_Scenario and HD_Scenario, the land-use types with the largest area outflow is BLF, with 4.8 × 103 km2 (18.61%) and 6 × 103 km2 (23.37%) of BLF will convert to CF in 2084 as a result of different climate changes and anthropogenic disturbances, respectively, and therefore the CF area will increase under both scenarios (Figure 13a,b). However, in the BD_Scenario and FD_Scenario, the percentage of CL outflow was both the highest due to the rapid urban development, and 5.2 × 103 km2 (18.11%) and 5.7 × 103 km2 (19.74%) of CL will convert to UL, respectively. These demonstrate a contrast between the ecological conservation scenarios and the urban development scenarios.

3.5. Analysis of Land-Use Change Amplitude at the Administrative Level

Based on the land-use pattern between 2014 and 2084, under different scenarios, the change amplitude of each land-use type at the Zhejiang Province administrative level was mapped as shown in Figure 14. The following three points are worth noting: (1) The BF area is in a state of growth in most of prefectures (+20% to +60%) under each scenario; (2) under the FD_Scenario, the CF area in northeast Zhejiang will be greatly reduced (–20% to –95%); (3) the BLF area in Jiaxing will decrease under either scenario (–24% to –90%).

4. Discussion

4.1. Future Enhancements of the BCS Model

The reasons for the great simulation results based on the BCS coupled model are as follows: (1) the SD model considers complex factors such as social, economic, policy, and planning, and the uncertainty in predicting the quantitative area is reduced, thus achieving more accurate predictions than other quantitative models; (2) the CA model integrates a BPNN to extract extremely accurate land suitability considering more comprehensive driving data.
However, the model simulation still leaves some errors, mainly in the simulation of three forest types. For BLF, CF, and BF, the accuracy was 0.75, 0.69, and 0.73, respectively. Forest succession needs to take into account not only climate change and human disturbance but also its own growth conditions [55]. How to incorporate the spatial distribution of these factors for more accurate modeling and prediction is a challenge to overcome for our future research. Furthermore, parameters, such as land suitability and conversion costs, in the BCS coupled model are assumed to remain constant during the simulation, whereas these parameters will change over time in realistic situations [26]. Studies have already been conducted to predict the future spatial distribution of the population and economy on the basis of which future land suitability can be predicted [56]. Therefore, we need to invest more effort in improving the model based on the appropriate approach in future works.

4.2. Future Strategy for Land-Use Management

Different scenarios have different levels of ecological protection and different shares of investment in forestry; therefore, the forest coverage rate is higher in scenarios that focus on ecological protection. Consequently, the areas of BLF and CF are higher in the SD_Scenario and HD_Scenario than in the BD_Scenario and FD_Scenario after 70 years (Figure 11e,f). However, the reasons why the area of BF is higher in the FD_Scenario than in others are for the following two reasons: (1) The growth of BF is affected by the climate and has certain requirements for precipitation and temperature. The annual total precipitation required for BF growth is 1200–2500 mm, and the average daily temperature conducive to BF growth ranges from 15–25 °C [45]. Despite the increasing radiation intensity after the mid-21st century, in the FD_Scenario it results in increasing temperature and decreasing precipitation, and the temperature and precipitation in Zhejiang Province eventually ranged from 16.8–22.5 °C and 858.2–1693 mm, respectively (Figure 5a,b). Therefore, it can be recognized that Zhejiang Province will still be a natural environment suitable for bamboo forest growth in this future scenario. (2) In recent years, with the adjustment of the industrial structure of Zhejiang Province, the economic benefits of BF have become prominent. Moso bamboo and lei bamboo, with higher economic value, have become new growth points in the regional economy. The management intensity of BF has increased, and the area of BF has continued to increase. In the future FD_Scenario, the rapid development of the social economy will expand the management and investment of BF to a certain extent. The SD model feedback system takes this man-made activity into account; thus, the BF area obtained by simulation was larger than that simulated in other scenarios.
In the four scenarios, the simulated future LUCC evolutions were consistent with realistic changes. In the future, the proportion of each land-use type converted ranges from 0 to 23.27% (Figure 13), with the highest values being the conversion of BLF to CF in the HD_Scenario (23.27%) and the conversion of CL to UL in the FD_Scenario (19.74%). In the HD_Scenario, due to the impact of policies, in order to prevent the massive reduction of CF, the planned goals are achieved by converting part of the BLF into CF. In the FD_Scenario, the man-land contradiction will intensify, and CL will be mainly converted into UL, making it difficult to achieve the goal of protecting CL and posing a serious challenge to the ecological environment. LUCC is caused by complex and diverse factors, and the factors we can consider are limited. Therefore, it is likely that the future change will exceed our expectations; that is, the LUCC amplitudes simulated by the model may be smaller than the actual future change.
As seen in the future area change (Figure 11), the future spatial distribution (Figure 12), the LUCC conversion (Figure 13), and the change at the prefecture level (Figure 14), the SD_Scenario has less land-use change than the other scenarios, and this scenario is the more ecologically friendly one. Without considering the promotion of social and economic development, to maximize the restoration of forest ecology, and promote ecological protection, various indicators established in the SD_Scenario can be referred to for future land planning and forest resource management.

5. Conclusions

In summary, the aim of this study was to develop a BPNN_CA_SD (BCS) coupled model for future LUCC simulation and to analyze the LUCC of Zhejiang Province under different scenarios from 2014 to 2084. The BCS coupled model consisted of the BPNN_CA model and the SD model. The top-down SD model and the bottom-up BPNN_CA model were interactively integrated during the simulation. The simulation results in 2014 showed a great OA (0.8), Kappa (0.75), and relatively high FOM (>28%) value, indicating that the proposed model can simulate LUCC accurately. Under different scenarios, the future evolution of the LUCC simulated by the BCS model varied due to the different natural and human effects. By 2084, bamboo forests would increase by 37% under the FD_Scenario, while coniferous forests would decline by 25%. A comparison of the simulated subtropical forest area and spatial variation in the four scenarios revealed that the SD_Scenario was favorable to forest ecology. We also analyzed the future transfer area between land uses and the changes in each prefecture. These study results could provide an effective reference for decision makers regarding sustainable forest development and land-use planning under future climate conditions in Zhejiang Province.

Author Contributions

Conceptualization, H.D.; methodology, Z.H.; software, Z.H.; validation, Z.H., X.L. (Xuejian Li) and F.M.; formal analysis, Z.H., N.H. and W.F.; investigation, Z.H., Y.X. and X.L. (Xin Luo); data curation, Z.H.; writing—original draft preparation, Z.H.; writing—review and editing, H.D.; supervision, H.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation (No. 32171785, U1809208, 31901310), the State Key Laboratory of Subtropical Silviculture (No. ZY20180201), the Zhejiang Provincial Collaborative Innovation Center for Bamboo Resources and High-Efficiency Utilization (No. S2017011), and the Key Research and Development Program of Zhejiang Province (No. 2021C02005).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the support of various foundations. The authors are grateful to the editor and anonymous reviewers whose comments have contributed to improving the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationDescription
LUCCLand Use and Land Cover Change
ULUrban Land
WBWater Body
CLCultivated Land
BLFBroad-Leaved Forest
CFConiferous Forest
BFBamboo Forest
CACellular Automata
SDSystem Dynamics
BPNNBack Propagation Neural Network
BPNN_CACA Model Integrated with the BPNN
BCSBPNN_CA Model Integrated with the SD
CLUE-SThe Conversion of Land Use and Its Effects at the Small Regional Extent
OAOverall Accuracy
KappaKappa Coefficients
PAProducer’s Accuracy
ROCReceiver Operating Characteristic
AUCArea under ROC Curve
F o M Figure of Merit
SD_ScenarioSlow Development Scenario
HD_ScenarioHarmonious Development Scenario
BD_ScenarioBase Development Scenario
FD_ScenarioFast Development Scenario
RCPRepresentative Concentration Pathway
CMIP5Coupled Model Intercomparison Project 5
IPCCIntergovernmental Panel on Climate Change

References

  1. Jiao, J.-G.; Yang, L.-Z.; Wu, J.-X.; Wang, H.-Q.; Li, H.-X.; Ellis, E.C. Land use and soil organic carbon in China’s village landscapes. Pedosphere 2010, 20, 1–14. [Google Scholar] [CrossRef] [Green Version]
  2. Matthews, H.D.; Graham, T.L.; Keverian, S.; Lamontagne, C.; Seto, D.; Smith, T.J. National contributions to observed global warming. Environ. Res. Lett. 2014, 9, 468–475. [Google Scholar] [CrossRef] [Green Version]
  3. Su, Z.-Y.; Xiong, Y.-M.; Zhu, J.-Y.; Ye, Y.-C.; Ye, M. Soil organic carbon content and distribution in a small landscape of Dongguan, South China. Pedosphere 2006, 16, 10–17. [Google Scholar] [CrossRef]
  4. Yue, C.; Ciais, P.; Houghton, R.A.; Nassikas, A.A. Contribution of land use to the interannual variability of the land carbon cycle. Nat. Commun. 2020, 11, 1–11. [Google Scholar] [CrossRef] [PubMed]
  5. Yu, Z.; Lu, C.; Tian, H.; Canadell, J.G. Largely underestimated carbon emission from land use and land cover change in the conterminous United States. Glob. Chang. Biol. 2019, 25, 3741–3752. [Google Scholar] [CrossRef] [PubMed]
  6. Brovkin, V.; Boysen, L.; Arora, V.K.; Boisier, J.P.; Cadule, P.; Chini, L.; Claussen, M.; Friedlingstein, P.; Gayler, V.; van den Hurk, B.J.J.M.; et al. Effect of anthropogenic land-use and land-cover changes on climate and land carbon storage in CMIP5 projections for the twenty-first century. J. Clim. 2013, 26, 6859–6881. [Google Scholar] [CrossRef]
  7. Khoshnood Motlagh, S.; Sadoddin, A.; Haghnegahdar, A.; Razavi, S.; Salmanmahiny, A.; Ghorbani, K. Analysis and prediction of land cover changes using the land change modeler (LCM) in a semiarid river basin, Iran. Land. Degrad. Dev. 2021, 32, 3092–3105. [Google Scholar] [CrossRef]
  8. Rasmussen, L.V.; Rasmussen, K.; Reenberg, A.; Proud, S. A system dynamics approach to land use changes in agro-pastoral systems on the desert margins of Sahel. Agric. Syst. 2012, 107, 56–64. [Google Scholar] [CrossRef]
  9. Liu, D.; Zheng, X.; Zhang, C.; Wang, H. A new temporal–spatial dynamics method of simulating land-use change. Ecol. Model. 2017, 350, 1–10. [Google Scholar] [CrossRef]
  10. Liu, X.; Ou, J.; Li, X.; Ai, B. Combining system dynamics and hybrid particle swarm optimization for land use allocation. Ecol. Model. 2013, 257, 11–24. [Google Scholar] [CrossRef]
  11. Zheng, X.-Q.; Zhao, L.; Xiang, W.-N.; Li, N.; Lv, L.-N.; Yang, X. A coupled model for simulating spatio-temporal dynamics of land-use change: A case study in Changqing, Jinan, China. Landsc. Urban Plan. 2012, 106, 51–61. [Google Scholar] [CrossRef]
  12. Luo, G.; Yin, C.; Chen, X.; Xu, W.; Lu, L. Combining system dynamic model and CLUE-S model to improve land use scenario analyses at regional scale: A case study of Sangong watershed in Xinjiang, China. Ecol. Complex. 2010, 7, 198–207. [Google Scholar] [CrossRef]
  13. Jiao, M.; Hu, M.; Xia, B. Spatiotemporal dynamic simulation of land-use and landscape-pattern in the Pearl River Delta, China. Sustain. Cities Soc. 2019, 49, 101581. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Li, Y.; Lv, J.; Wang, J.; Wu, Y. Scenario simulation of ecological risk based on land use/cover change—A case study of the Jinghe county, China. Ecol. Indic. 2021, 131, 108176. [Google Scholar] [CrossRef]
  15. Gomes, E.; Inácio, M.; Bogdzevič, K.; Kalinauskas, M.; Karnauskaitė, D.; Pereira, P. Future land-use changes and its impacts on terrestrial ecosystem services: A review. Sci. Total Environ. 2021, 781, 146716. [Google Scholar] [CrossRef]
  16. He, M.; Chen, C.; Zheng, F.; Chen, Q.; Zhang, J.; Yan, H.; Lin, Y. An efficient dynamic route optimization for urban flooding evacuation based on Cellular Automata. Comput. Environ. Urban Syst. 2021, 87, 101622. [Google Scholar] [CrossRef]
  17. Priem, F.; Canters, F. Modelling transitions in sealed surface cover fraction with Quantitative State Cellular Automata. Landsc. Urban Plan. 2021, 211, 104081. [Google Scholar] [CrossRef]
  18. Aguejdad, R. The Influence of the Calibration Interval on Simulating Non-Stationary Urban Growth Dynamic Using CA-Markov Model. Remote Sens. 2021, 13, 468. [Google Scholar] [CrossRef]
  19. Li, K.; Feng, M.; Biswas, A.; Su, H.; Niu, Y.; Cao, J. Driving Factors and Future Prediction of Land Use and Cover Change Based on Satellite Remote Sensing Data by the LCM Model: A Case Study from Gansu Province, China. Sensors 2020, 20, 2757. [Google Scholar] [CrossRef]
  20. Surabuddin Mondal, M.; Sharma, N.; Kappas, M.; Garg, P.K. Modeling of spatio-temporal dynamics of land use and land cover in a part of Brahmaputra River basin using Geoinformatic techniques. Geocarto. Int. 2013, 28, 632–656. [Google Scholar] [CrossRef]
  21. Jokar Arsanjani, J.; Helbich, M.; Kainz, W.; Darvishi Boloorani, A. Integration of logistic regression, Markov chain and cellular automata models to simulate urban expansion. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 265–275. [Google Scholar] [CrossRef]
  22. Mallick, J.; AlQadhi, S.; Talukdar, S.; Pradhan, B.; Bindajam, A.A.; Islam, A.R.M.T.; Dajam, A.S. A Novel Technique for Modeling Ecosystem Health Condition: A Case Study in Saudi Arabia. Remote Sens. 2021, 13, 2632. [Google Scholar] [CrossRef]
  23. Roodposhti, M.S.; Aryal, J.; Bryan, B.A. A novel algorithm for calculating transition potential in cellular automata models of land-use/cover change. Environ. Model. Software 2019, 112, 70–81. [Google Scholar] [CrossRef]
  24. Chen, Z.; Huang, M.; Zhu, D.; Altan, O. Integrating Remote Sensing and a Markov-FLUS Model to Simulate Future Land Use Changes in Hokkaido, Japan. Remote Sens. 2021, 13, 2621. [Google Scholar] [CrossRef]
  25. Marondedze, A.K.; Schütt, B. Predicting the Impact of Future Land Use and Climate Change on Potential Soil Erosion Risk in an Urban District of the Harare Metropolitan Province, Zimbabwe. Remote Sens. 2021, 13, 4360. [Google Scholar] [CrossRef]
  26. Liu, X.; Liang, X.; Li, X.; Xu, X.; Ou, J.; Chen, Y.; Li, S.; Wang, S.; Pei, F. A future land use simulation model (FLUS) for simulating multiple land use scenarios by coupling human and natural effects. Landsc. Urban Plan. 2017, 168, 94–116. [Google Scholar] [CrossRef]
  27. Liu, Y.; Jing, W.; Xu, L. Parallelizing backpropagation neural network using MapReduce and cascading model. Comput. Intell. Neurosci. 2016, 2016, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Zhou, G.; Liu, S.; Li, Z.; Zhang, D.; Tang, X.; Zhou, C.; Yan, J.; Mo, J. Old-growth forests can accumulate carbon in soils. Science 2006, 314, 1417. [Google Scholar] [CrossRef] [Green Version]
  29. Tan, Z.-H.; Zhang, Y.-P.; Schaefer, D.; Yu, G.-R.; Liang, N.; Song, Q.-H. An old-growth subtropical Asian evergreen forest as a large carbon sink. Atmos. Environ. 2011, 45, 1548–1554. [Google Scholar] [CrossRef]
  30. Fang, J.; Yu, G.; Liu, L.; Hu, S.; Chapin, F.S. Climate change, human impacts, and carbon sequestration in China. Proc. Natl. Acad. Sci. USA 2018, 115, 4015–4020. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, F.-h.; Xu, C.-Y.; Yang, X.-x.; Ye, X.-c. Controls of climate and land-use change on terrestrial net primary productivity variation in a subtropical humid basin. Remote Sens. 2020, 12, 3525. [Google Scholar] [CrossRef]
  32. Guo, X.; Chen, H.Y.H.; Meng, M.; Biswas, S.R.; Ye, L.; Zhang, J. Effects of land use change on the composition of soil microbial communities in a managed subtropical forest. For. Ecol. Manag. 2016, 373, 93–99. [Google Scholar] [CrossRef]
  33. Zhu, K.; Zhang, J.; Niu, S.; Chu, C.; Luo, Y. Limits to growth of forest biomass carbon sink under climate change. Nat. Commun. 2018, 9, 2709. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Yao, Y.; Piao, S.; Wang, T. Future biomass carbon sequestration capacity of Chinese forests. Sci. Bull. 2018, 63, 1108–1117. [Google Scholar] [CrossRef] [Green Version]
  35. Poelmans, L.; Van Rompaey, A. Complexity and performance of urban expansion models. Comput. Environ. Urban Syst. 2010, 34, 17–27. [Google Scholar] [CrossRef]
  36. Feng, Y.; Tong, X. Incorporation of spatial heterogeneity-weighted neighborhood into cellular automata for dynamic urban growth simulation. GISci. Remote Sens. 2019, 56, 1024–1045. [Google Scholar] [CrossRef]
  37. Van Vuuren, D.P.; Edmonds, J.; Kainuma, M.; Riahi, K.; Thomson, A.; Hibbard, K.; Hurtt, G.C.; Kram, T.; Krey, V.; Lamarque, J.-F.; et al. The representative concentration pathways: An overview. Clim. Chang. 2011, 109, 5–31. [Google Scholar] [CrossRef]
  38. Jefferson, M. IPCC fifth assessment synthesis report: “Climate change 2014: Longer report”: Critical analysis. Technol. Forecast. Soc. Chang. 2015, 92, 362–363. [Google Scholar] [CrossRef]
  39. Li, J.; Chen, X.; Kurban, A.; Van de Voorde, T.; De Maeyer, P.; Zhang, C. Coupled SSPs-RCPs scenarios to project the future dynamic variations of water-soil-carbon-biodiversity services in Central Asia. Ecol. Indic. 2021, 129, 107936. [Google Scholar] [CrossRef]
  40. Dubovyk, O.; Sliuzas, R.; Flacke, J. Spatio-temporal modelling of informal settlement development in Sancaktepe district, Istanbul, Turkey. ISPRS J. Photogramm. Remote Sens. 2011, 66, 235–246. [Google Scholar] [CrossRef]
  41. Zheng, J.; Mao, F.; Du, H.; Li, X.; Zhou, G.; Dong, L.; Zhang, M.; Han, N.; Liu, T.; Xing, L. Spatiotemporal simulation of net ecosystem productivity and its response to climate change in subtropical forests. Forests 2019, 10, 708. [Google Scholar] [CrossRef] [Green Version]
  42. Li, X.; Du, H.; Mao, F.; Zhou, G.; Xing, L.; Liu, T.; Han, N.; Liu, E.; Ge, H.; Liu, Y.; et al. Mapping spatiotemporal decisions for sustainable productivity of bamboo forest land. Land. Degrad. Dev. 2020, 31, 939–958. [Google Scholar] [CrossRef]
  43. Li, Y.; Han, N.; Li, X.; Du, H.; Mao, F.; Cui, L.; Liu, T.; Xing, L. Spatiotemporal estimation of bamboo forest aboveground carbon storage based on Landsat data in Zhejiang, China. Remote Sensing 2018, 10, 898. [Google Scholar] [CrossRef] [Green Version]
  44. Zhang, M.; Du, H.; Mao, F.; Zhou, G.; Li, X.; Dong, L.; Zheng, J.; Zhu, D.e.; Liu, H.; Huang, Z.; et al. Spatiotemporal evolution of urban expansion using Landsat time series data and assessment of its influences on forests. ISPRS Int. J. Geoinf. 2020, 9, 64. [Google Scholar] [CrossRef] [Green Version]
  45. Li, X.; Mao, F.; Du, H.; Zhou, G.; Xing, L.; Liu, T.; Han, N.; Liu, Y.; Zhu, D.e.; Zheng, J.; et al. Spatiotemporal evolution and impacts of climate change on bamboo distribution in China. J. Environ. Manag. 2019, 248, 109265. [Google Scholar] [CrossRef]
  46. Tian, H.; Liang, X.; Li, X.; Liu, X.; Ou, J.; Hong, Y.; He, Z. Simulating Multiple Land Use Scenarios in China during 2010–2050 Based on System Dynamic Model. Trop. Geogr. 2017, 37, 547–561. [Google Scholar]
  47. Vuuren, D.; Stehfest, E.; Elzen, M.; Kram, T.; Vliet, J.V.; Deetman, S.; Isaac, M.; Goldewijk, K.K.; Hof, A.; Beltran, A.M. RCP2.6: Exploring the possibility to keep global mean temperature increase below 2°C. Clim. Chang. 2011, 109, 95–116. [Google Scholar] [CrossRef]
  48. Chaturvedi, R.K.; Joshi, J.; Jayaraman, M.; Bala, G.; Ravindranath, N. Multi-model climate change projections for India under representative concentration pathways. Curr. Sci. 2012, 103, 791–802. [Google Scholar]
  49. Riahi, K.; Rao, S.; Krey, V.; Cho, C.; Chirkov, V.; Fischer, G.; Kindermann, G.; Nakicenovic, N.; Rafaj, P. RCP 8.5—A scenario of comparatively high greenhouse gas emissions. Clim. Chang. 2011, 109, 33–57. [Google Scholar] [CrossRef] [Green Version]
  50. Zou, L.; Zhou, T. Near future (2016–40) summer precipitation changes over China as projected by a regional climate model (RCM) under the RCP8.5 emissions scenario: Comparison between RCM downscaling and the driving GCM. Adv. Atmos. Sci. 2013, 30, 806–818. [Google Scholar] [CrossRef]
  51. Xu, X.; Du, Z.; Zhang, H. Integrating the system dynamic and cellular automata models to predict land use and land cover change. Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 568–579. [Google Scholar] [CrossRef]
  52. Liang, X.; Guan, Q.; Clarke, K.C.; Liu, S.; Wang, B.; Yao, Y. Understanding the drivers of sustainable land expansion using a patch-generating land use simulation (PLUS) model: A case study in Wuhan, China. Comput. Environ. Urban Syst. 2021, 85, 101569. [Google Scholar] [CrossRef]
  53. Xu, H.; Song, Y.; Tian, Y. Simulation of land-use pattern evolution in hilly mountainous areas of North China: A case study in Jincheng. Land Use Policy 2022, 112, 105826. [Google Scholar] [CrossRef]
  54. Pontius, R.G.; Boersma, W.; Castella, J.; Clarke, K.C.; De Nijs, T.; Dietzel, C.; Duan, Z.; Fotsing, E.; Goldstein, N.; Kok, K. Comparing the input, output, and validation maps for several models of land change. Ann. Regional Sci. 2008, 42, 11–37. [Google Scholar] [CrossRef] [Green Version]
  55. Wu, Z.; Dai, E.; Wu, Z.; Lin, M. Future forest dynamics under climate change, land use change, and harvest in subtropical forests in Southern China. Landsc. Ecol. 2019, 34, 843–863. [Google Scholar] [CrossRef]
  56. Reimann, L.; Jones, B.; Nikoletopoulos, T.; Vafeidis, A.T. Accounting for internal migration in spatial population projections—A gravity-based modeling approach using the Shared Socioeconomic Pathways. Environ. Res. Lett. 2021, 16, 074025. [Google Scholar] [CrossRef]
Remotesensing 14 01698 g001 550
Figure 1. Zhejiang Province study area: (a) location in China, (b) sample plots in 2014, (c) forest coverage rate from 2004 to 2020, and (d) land-use patterns in 2014.
Figure 1. Zhejiang Province study area: (a) location in China, (b) sample plots in 2014, (c) forest coverage rate from 2004 to 2020, and (d) land-use patterns in 2014.
Remotesensing 14 01698 g001
Remotesensing 14 01698 g002 550
Figure 2. (a) Land-use patterns and (b) accuracy evaluation of land-use classification in Zhejiang Province from 1984 to 2014.
Figure 2. (a) Land-use patterns and (b) accuracy evaluation of land-use classification in Zhejiang Province from 1984 to 2014.
Remotesensing 14 01698 g002
Remotesensing 14 01698 g003 550
Figure 3. Topographic data in 2014: (a) DEM, (b) slope, and (c) aspect; soil data in 2008: (d) clay fraction, (e) sand fraction, (f) silt fraction, (g) soil available water content, (h) soil bulk density, and (i) soil wilt point; climate data in 2014: (j) total precipitation, (k) average temperature, (l) average radiation, and (m) average relative humidity; socioeconomic data in 2015: (n) population density, (o) GDP; distance data in 2014: (p) distance to road, (q) distance to railway, and (r) distance to water.
Figure 3. Topographic data in 2014: (a) DEM, (b) slope, and (c) aspect; soil data in 2008: (d) clay fraction, (e) sand fraction, (f) silt fraction, (g) soil available water content, (h) soil bulk density, and (i) soil wilt point; climate data in 2014: (j) total precipitation, (k) average temperature, (l) average radiation, and (m) average relative humidity; socioeconomic data in 2015: (n) population density, (o) GDP; distance data in 2014: (p) distance to road, (q) distance to railway, and (r) distance to water.
Remotesensing 14 01698 g003
Remotesensing 14 01698 g004 550
Figure 4. Macro statistics: (a) total population; (b) GDP; (c) grain yield; (d) aquatic product yield.
Figure 4. Macro statistics: (a) total population; (b) GDP; (c) grain yield; (d) aquatic product yield.
Remotesensing 14 01698 g004
Remotesensing 14 01698 g005 550
Figure 5. Configurations of four development scenarios concerning human and natural effects: (a) annual total precipitation; (b) annual average temperature; (c) annual average radiation; (d) annual average relative humidity under four scenarios.
Figure 5. Configurations of four development scenarios concerning human and natural effects: (a) annual total precipitation; (b) annual average temperature; (c) annual average radiation; (d) annual average relative humidity under four scenarios.
Remotesensing 14 01698 g005
Remotesensing 14 01698 g006 550
Figure 6. Flow chart used in this study.
Figure 6. Flow chart used in this study.
Remotesensing 14 01698 g006
Remotesensing 14 01698 g007 550
Figure 7. The structure of the SD model. Variables in red indicate the inputs of the model, which varied in different scenarios. Variables in purple represent the main outputs, including the areas of different LULC types. Purple ones were used to validate the model.
Figure 7. The structure of the SD model. Variables in red indicate the inputs of the model, which varied in different scenarios. Variables in purple represent the main outputs, including the areas of different LULC types. Purple ones were used to validate the model.
Remotesensing 14 01698 g007
Remotesensing 14 01698 g008 550
Figure 8. The interactive coupling mechanism of the SD model and BPNN_CA model.
Figure 8. The interactive coupling mechanism of the SD model and BPNN_CA model.
Remotesensing 14 01698 g008
Remotesensing 14 01698 g009 550
Figure 9. The actual land-use pattern in 2008 and 2014 and the simulated land-use pattern in 2014. (A,B) are the two sub-areas shown in magnification.
Figure 9. The actual land-use pattern in 2008 and 2014 and the simulated land-use pattern in 2014. (A,B) are the two sub-areas shown in magnification.
Remotesensing 14 01698 g009
Remotesensing 14 01698 g010 550
Figure 10. Verification: (a) comparison between the classified and simulated areas using the SD model; (b) ROC curves and AUC values fitted by the BPNN; (c) normalized confusion matrix of the simulated result and the actual classification in 2014; (d) the FOM from 2008 to 2014.
Figure 10. Verification: (a) comparison between the classified and simulated areas using the SD model; (b) ROC curves and AUC values fitted by the BPNN; (c) normalized confusion matrix of the simulated result and the actual classification in 2014; (d) the FOM from 2008 to 2014.
Remotesensing 14 01698 g010
Remotesensing 14 01698 g011 550
Figure 11. Areas of each land-use type in Zhejiang Province over the past 30 years and that over the next 70 years under different scenarios: (a) urban land area; (b) water body area; (c) cultivated land area; (d) bamboo forest area; (e) broad-leaved forest area; (f) coniferous forest area.
Figure 11. Areas of each land-use type in Zhejiang Province over the past 30 years and that over the next 70 years under different scenarios: (a) urban land area; (b) water body area; (c) cultivated land area; (d) bamboo forest area; (e) broad-leaved forest area; (f) coniferous forest area.
Remotesensing 14 01698 g011
Remotesensing 14 01698 g012 550
Figure 12. Spatial distribution of land use in different regions of Zhejiang Province in 2084 under different scenarios. (A,B) are the two sub-areas shown in magnification.
Figure 12. Spatial distribution of land use in different regions of Zhejiang Province in 2084 under different scenarios. (A,B) are the two sub-areas shown in magnification.
Remotesensing 14 01698 g012
Remotesensing 14 01698 g013 550
Figure 13. Predicted land-use conversion in 2014–2084 under 4 scenarios: (a) SD_Scenario; (b) HD_Scenario; (c) BD_Scenario; (d) FD_Scenario. Numbers: years; UL: urban land; WB: water body; CL: cultivated land; BLF: broad-leaved forest; CF: coniferous forest; BF: bamboo forest. For example, 14BF represents bamboo forest in 2014.
Figure 13. Predicted land-use conversion in 2014–2084 under 4 scenarios: (a) SD_Scenario; (b) HD_Scenario; (c) BD_Scenario; (d) FD_Scenario. Numbers: years; UL: urban land; WB: water body; CL: cultivated land; BLF: broad-leaved forest; CF: coniferous forest; BF: bamboo forest. For example, 14BF represents bamboo forest in 2014.
Remotesensing 14 01698 g013
Remotesensing 14 01698 g014 550
Figure 14. Land use/cover change amplitude at the municipal administrative level in Zhejiang Province from 2014 to 2084 under different future scenarios.
Figure 14. Land use/cover change amplitude at the municipal administrative level in Zhejiang Province from 2014 to 2084 under different future scenarios.
Remotesensing 14 01698 g014
Table
Table 2. The number of verification plots in different years.
Table 2. The number of verification plots in different years.
YearULWBCLBLFCFBFTotal
19841511413022043851141297
19881571042871643171491178
19921631122371962661821156
19961771342671442151701107
2000128146139159165232969
2004128142139142152215918
2008123128127127127246878
2014123132138147154139833
Table
Table 3. Parameter settings of different socioeconomic development scenarios.
Table 3. Parameter settings of different socioeconomic development scenarios.
FactorsPatternsAnnual Growth Rate Settings from 2014 to 2084
PopulationHigh growth (P1)7.2‰ average from 2004 to 2014
Steady growth (P2)Growth rate simulated by logistic population retardation growth model
Moderate growth (P3)0.85× growth rate simulated by logistic population retardation growth model
Slow growth (P4)7.2‰ linearly down to 3.4‰
GDPHigh growth (G1)14% average from 2004 to 2014
Steady growth (G2)14% linearly down to 10.5%
Moderate growth (G3)14% linearly down to 8%
Slow growth (G4)14% linearly down to 6.5%
TechnologyRapid innovation (T1)Grain yieldMaintain 5‰ in 2014
Aquatic yieldMaintain 8% in 2014
Steady innovation (T2)Grain yield5‰ linearly down to 3‰
Aquatic yield8% linearly down to 5%
Moderate innovation (T3)Grain yield5% linearly down to 1%
Aquatic yield8% linearly down to 2%
No innovation (T4)Grain yield0%
Aquatic yield0%
EcologyHigher forest coverage rate (E1)
High forest coverage rate (E2)
Medium forest coverage rate (E3)
Low forest coverage rate (E4)		60.89% linearly up to 65%
60.89% linearly up to 63%
60.89% linearly up to 61%
60.89% linearly down to 60%
The logistic population retardation growth model [46] was constructed to calculate p m and r . p m is the maximum population of 4,328,457, and r is the population growth rate of 6.097%.
Table
Table 4. The conversion cost of different land-use types from 2004 to 2084.
Table 4. The conversion cost of different land-use types from 2004 to 2084.
TypesULWBCLBLFCFBF
UL00.850.70.990.990.99
WB0.800.80.90.90.9
CL0.30.700.50.40.5
BLF0.90.90.700.60.5
CF0.90.90.60.500.6
BF0.90.90.70.60.60
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, Z.; Li, X.; Du, H.; Mao, F.; Han, N.; Fan, W.; Xu, Y.; Luo, X. Simulating Future LUCC by Coupling Climate Change and Human Effects Based on Multi-Phase Remote Sensing Data. Remote Sens. 2022, 14, 1698. https://doi.org/10.3390/rs14071698

AMA Style

Huang Z, Li X, Du H, Mao F, Han N, Fan W, Xu Y, Luo X. Simulating Future LUCC by Coupling Climate Change and Human Effects Based on Multi-Phase Remote Sensing Data. Remote Sensing. 2022; 14(7):1698. https://doi.org/10.3390/rs14071698

Chicago/Turabian Style

Huang, Zihao, Xuejian Li, Huaqiang Du, Fangjie Mao, Ning Han, Weiliang Fan, Yanxin Xu, and Xin Luo. 2022. "Simulating Future LUCC by Coupling Climate Change and Human Effects Based on Multi-Phase Remote Sensing Data" Remote Sensing 14, no. 7: 1698. https://doi.org/10.3390/rs14071698

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