Next Article in Journal
The Effect of Surgical Masks on the Featural and Configural Processing of Emotions
Previous Article in Journal
Thoron Interference on Performance of Continuous Radon Monitors: An Experimental Study on Four Devices and a Proposal of an Indirect Method to Estimate Thoron Concentration
Previous Article in Special Issue
Analysis for the Interaction Relationship between Urbanization and Ecological Security: A Case Study in Wuhan City Circle of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatio-Temporal Evolution, Prediction and Optimization of LUCC Based on CA-Markov and InVEST Models: A Case Study of Mentougou District, Beijing

1
Key Laboratory of National Forestry and Grassland Administration on Ecological Landscaping of Challenging Urban Sites, National Innovation Alliance of National Forestry and Grassland Administration on Afforestation and Landscaping of Challenging Urban Sites, Shanghai Engineering Research Center of Landscaping on Challenging Urban Sites, Shanghai Academy of Landscape Architecture Science and Planning, Shanghai 200232, China
2
School of Design, Shanghai Jiao Tong University, Shanghai 200240, China
3
Shanghai Foundation Ding Environmental Technology Company, Shanghai 200063, China
4
Jiangxi Institute of Ecological Civilization, School of Resources, Environmental & Chemical Engineering, Nanchang University, Nanchang 330031, China
5
Taihu Basin Monitoring Central Station for Soil and Water Conservation, Taihu Basin Authority of Ministry of Water Resources, Shanghai 200434, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(4), 2432; https://doi.org/10.3390/ijerph19042432
Submission received: 26 November 2021 / Revised: 9 February 2022 / Accepted: 9 February 2022 / Published: 19 February 2022

Abstract

:
With the rapid advancement of urbanization and industrialization, the contradiction between the social economy and resources and the environment has become increasingly prominent. On the basis of limited land resources, the way to promote multi-objective comprehensive development such as economic, social development and ecological and environmental protection through structure and layout regulation, so as to maximize regional comprehensive benefits, is an important task of current land spatial planning. Our aim is to obtain land-use-change data in the study area using remote-sensing data inversion and multiple-model simulation. Based on land suitability evaluation, we predict and optimize the land use structure of the study area in 2030 and evaluate and compare ecosystem services. Based on remote-sensing images and eco-environmental data from 1985 to 2014 in the study area, land use/land cover change (LUCC) and future simulation data were obtained by using supervised classification, landscape metrics and the CA-Markov model. The ecosystem services were evaluated by the InVEST model. The analytic hierarchy process (AHP) method was used to evaluate the land suitability for LUCC. Finally, the LUCC in 2030 under two different scenarios, Scenario_1 (prediction) and Scenario_2 (optimization), were evaluated, and the ecosystem service functions were compared. In the last 30 years, the landscape in the study area has gradually fragmented, and the built-up land has expanded rapidly, increased by one-third, mainly at the cost of cropland, orchards and wasteland. According to the suitability evaluation, giving priority to the land use types with higher environmental requirements will ensure the study area has a higher ecosystem service value. The rapid development of urbanization has a far-reaching impact on regional LUCC. Intensive land resources need reasonable and scientific land use planning, and land use planning should be based on the suitability evaluation of land resources, which can improve the regional ecosystem service function.

1. Introduction

Under the influence of natural and human factors, land use/land cover change (LUCC) has had a latent influence on the material cycle and energy conversion in the ecological chain with the continuous advancement of urbanization. These impacts alter habitat quality and ecological change processes, and ultimately affect the ecosystem structure and function [1,2,3,4,5]. With the rapid economic development in recent years, the process of global urbanization has been accelerated, and the world’s urban population has reached half of the total population by 2020 [6]. The Chinese urbanization rate has increased from 36.22% in 2000 to 63.89% in 2020 [7]. By 2030, The Chinese urban population will account for about three-quarters of the total population [8]. Rapid urbanization has brought about a series of problems related to urban expansion, which might lead to inefficient LUCC, the heat island effect, soil erosion and human settlement deterioration [9,10,11]. According to the Millennium Ecosystem Assessment (MA), more than 60% of global Ecosystem Services (ES) were degraded or overutilized, and LUCC change caused by human activities was one of the important driving factors leading to the degradation of ES [12]. Therefore, monitoring, evaluating, forecasting and optimizing the evolution of LUCC and the resulting evolution of habitat and ES, so as to improve the quality of the urban ecological environment and the human settlement environment, was of great practical significance to maintain regional sustainable development [13,14].
In recent years, spatial-temporal forecasting has been used to help optimize the development direction of urban planning, so as to obtain the urban development dynamics and ES capacity in advance and provide timely feedback and adjustment, which has achieved good results [15,16]. Cellular Automata (CA), as one of the most important and widely used methods across many models, was the basis of many models, such as the Logistic-CA model [17], ANN-CA model [18], CLUE-S model [19] and FLUS model [20]. The CA-Markov model was a relatively successful simulation method, which combined the ability of the CA model to simulate the spatial changes of complex systems and the advantage of the Markov model in long-term prediction [21,22,23]. It not only improved the prediction accuracy of LUCC transformation, but also effectively simulated the spatial change of LUCC, which had good scientific and practical characteristics [24]. Jenerette et al. [25] used the CA-Markov model to study the impact of urbanization and population growth on landscape changes in Phoenix, Arizona (USA). Nourgolipour et al. [26] simulated the spatial pattern distribution of local palm plantations in the Kuala Lengyue region of Malaysia by combining CA-Markov and MCE models. Alcamo et al. [27] set four different scenarios and found that the tradeoff between different ES will be increasingly intensified in the 21st century. These studies not only provided a good research paradigm for LUCC management and planning at different scales, but also provide important information about the implications for the ecological processes and ES changed by LUCC.
The concept of ES was originally derived from the concept of the environmental service function proposed in the research on key environmental issues, and then gradually developed by Ehrlich et al. (1987) [28]. Since then, the value evaluation of global ES and natural capital published by Constanza et al. (1997) has caused the research on ES to gradually become a hot topic in global ecological research [29]. In 2000, the launch of the MA once again promoted the development of ES research [30]. The project divided ES into four functions: Supply, culture, support and regulation [31]. It not only gave the scope and mode of ES value, but also pointed out the relationship between human well-being and ES [32]. The project defined ES as human benefits from natural ecosystems [33]. The evaluation methods of ES could be divided into two types: Value quantity and material quality evaluation [34].
At present, the use of remote sensing to assess ES has become a mainstream trend. It mainly carries out the key characteristic parameters of ecological assets, monitoring the changing trend at different times and the application of ES evaluation with the help of mature remote-sensing means and spatial analysis methods [35,36,37]. The ES assessment models that attracted scholars' attention include Multiscale Integrated Models of Ecosystem Services (MIMES) [38], the Global Unified Metamodel of the Biosphere (GUMBO) [39] and Artificial Intelligence for Ecosystem Services (ARIES) [40] developed by the University of Vermont, the Integrated Biosphere Simulator (IBIS) [41] developed by the University of Wisconsin, CITYgreen [42] developed by the U.S. Forestry Administration and the InVEST model developed by Stanford University [43]. Among them, the InVEST model has been more widely used because of its accurate quantification, visualization of results and low cost [44]. The InVEST model was mainly divided into marine, land and freshwater modules, and each module was specifically divided into several assessment projects. The assessment mainly included three aspects: The change in ES caused by the change in ecological environment dominated by the change in LUCC and landscape patterns [45]; crop yield evaluation, grassland resource evaluation and forest volume estimation aiming at evaluation and application [46]; and ES driven by ecosystem characteristic parameters, such as the net primary productivity, vegetation coverage and leaf area index [47].
The Mentougou District of Beijing is the western mountainous area of Beijing, the capital of China. As the western green barrier of the capital, it was also the only pure mountainous area in Beijing, providing important ES functions [48,49]. In the mountainous areas of the Chinese capital, land resources have economic, social, natural and ecological attributes [50,51]. This study took the Mentougou District of Beijing as the research object and took seven issues of remote-sensing images and other data (including terrain, landform, climate, soil, etc.) from 1985 to 2014 as the basic data source. Using the methods of supervised classification, landscape metrics, spatial statistical analysis, scenario simulation, grey-linear model, analytic hierarchy process and ES evaluation, the LUCC and ES function in the study area were used, evaluated and predicted. The purpose of this study was to (1) analyze the temporal and spatial transfer characteristics and landscape patterns of land use in the study area from 1985 to 2014; (2) carry out land suitability evaluation and scenario prediction and simulation of future land use; (3) evaluate the differences of ES functions in the study area in 2030 under different scenarios. The intention of the present study was to provide a scientific reference for the formulation of regional ecological environment construction and sustainable development.

2. Materials and Methods

2.1. Study Area

The Mentougou District is located in the northwest of Beijing, northern China (115°25′ E–116°10′ E, 39°48′ N–40°10′ N), with a total area of 1447.85 km2 (Figure 1). The region has a mid-latitude continental monsoon climate, and 98.5% of the region is mountainous [52]. In 2019, the average annual precipitation in Mentougou District was 405.7 mm and the average temperature was 13.8 °C [53]. Mentougou district has 4 streets and 9 towns, and in 2019, 254,000 people registered for residence and 344,000 permanent residents, of whom 42,000 belonged to the agricultural population and 211,000 to the non-agricultural population [54]. The regional GDP is 391.33 million dollars, and the primary, secondary and tertiary industries accounted for 1.24%, 26.94% and 71.84%, respectively [55].

2.2. LUCC Classification and Quantification of Landscape Patterns

In this study, we used remote-sensing image data of the Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM) and Operational Land Imager (OLI) for nearly 30 years (one image per year from 1985, 1990, 1995, 2000, 2005, 2010 and 2014) [56,57]. Based on the Chinese land use classification system and the characteristics of the study area, land use is divided into nine types through supervised classification [58,59]. LUCC includes cropland (CL), orchard (OC), forestland (FL), shrubland (SL), grassland (GL), bareland (BL), waterbodies (WB), wasteland (WL) and built-up land (BUL), forming a grid image with a spatial resolution of 30 m (Table A1). Using the method of stratified random sampling, we randomly selected f 270 points (30 points for each land use type) in the study area, using higher-resolution remote-sensing images and related planning data (before 2005) and field survey (after 2005), and tested the accuracy of interpretation, and the overall accuracy of data interpretation for each year was over 85%.

2.2.1. LUCC Transfer Matrix

Through the analysis of the land use transfer matrix, two different periods of transformation could be obtained. In order to clearly express the data of land use transfer from 1990 to 2015, we made a graph called circos and listed the data as an Appendix (Figure 2 and Table A2, Table A3, Table A4, Table A5, Table A6, Table A7 and Table A8). The mapping method using the template is from Canada’s Michael Smith Genome Sciences Centre (CMSGSC) (http://mkweb.bcgsc.ca/tableviewer/ accessed on 14 August 2021) [60,61].

2.2.2. Statistical Analysis

To assess changes in the structural characteristics of landscape patterns at the sub-watershed level, we selected the number of patches (NP), mean patch size (MPS), patch density (PD) and patch cohesion index (COHESION), to characterize landscape patterns (Table 1) [62,63]. All calculations were extracted from the FRAGSTATS 4.2 software (Developed by the Clark Labs, Clark University).

2.3. AHP and CA-Markov

According to ES demand and natural climate factors, a land suitability evaluation map was constructed (Figure 3). The analytic hierarchy process of determining the importance of each indicator was conducted by more than 20 experts in ecology-related fields. All results were tested for consistency using Yaahp software (developed by the Shanxi Yuan Decision Software Technology Co., LTD in China). Spatial simulation is based on the CA-Markov module in IDRIS 17.0 software (developed by the Clark Laboratory of Clark University in the United States) to simulate LUCC changes in the study area [64,65,66].
The process was as follows (Figure 4): (1) Determine the transition area matrix and transition probability matrix of LUCC in the study area from 1985 to 2014, which were used to test the accuracy and predict the results with the CA-Markov model. (2) Then, we established the atlas of LUCC change suitability rules; the spatial distribution probability atlas of each LUCC was obtained by analyzing the limiting factors and driving factors as independent variables. (3) Next, the CA filter was constructed; in this study, a 5 × 5 mole neighborhood [21] (four adjacent cells above and below a cell are domains) was used as the filtering parameter of the CA-Markov model. (4) We then determined the starting time and the number of iterations. Firstly, the period from 1990 to 2010 was taken as the starting time of prediction, and the number of CA iterations was set as 10 to simulate the spatial distribution of land use types in the study area in 2010. Compared with the real results, kappa coefficients of all types of land use were > 0.70%, which was feasible. (5) Taking 1990 to 2014 as the starting point of prediction, the number of iterations was set as 16 to predict the LUCC in the study area in 2030. The prediction result was obtained (Scenario_1). (6) According to the land quantity structure and land suitability evaluation results, the land use transformation was restrained and controlled. The gray-linear model was used to solve the quantity of land use with equal emphasis on ecology and economy (Table A9). According to the results of the suitability analysis, priority should be given to the land use types requiring higher environmental conditions, and these types with higher land-use requirements should be allocated to the appropriate space, so as to obtain the optimization of the LUCC spatial distribution in 2030 (Scenario_2, regarding which land-use types are allocated to more suitable areas). (7) Lastly, the climate factors in 2014 were brought into the InVEST model to calculate and compare the ES of the predicted and optimized LUCC in 2030.
The kappa coefficient is calculated as follows [62]:
Kappa = P o P c P p P c
where P o is the observed consistency rate between the reference data and the simulation data, P c is the expected correct simulation proportion in the case of randomness and P p is the correct simulation proportion in the case of ideal classification, which is generally 1.

2.4. Methods to Assessment ES

2.4.1. Water Yield

The water yield (WY) module of the InVEST model can estimate the WY function of different LUCC. The core algorithm is used to calculate the WY of each grid by using the water balance method combined with climate, terrain and LUCC types. The WY is the precipitation of each grid in the region minus the actual evapotranspiration. The calculation formula of WY is as follows [67].
Y x = P x T A E x
T A E x / P x = ( 1 + w x R x ) / ( 1 + w x R x + 1 / R x )
R x = T P E x / P x
C K x = T P E x / T E 0 x
T A E x = min ( C K x × T E 0 x , P x )
w x = Z × ( C A w x / P x ) + 1.25
Y x is the annual water yield of grid x (mm); T A E x is the actual annual evapotranspiration of grid x (mm); P x is the precipitation (mm) of grid x (mm); T A E x is approximated from the Budyko curve [67]; R x is the Budyko dryness index of grid x, which is the ratio of potential evapotranspiration ( T P E x ) to precipitation ( P x ); C K x is the evapotranspiration coefficient of vegetation, which is different in different vegetation types, and represents the ratio of potential evapotranspiration to reference evapotranspiration ( T E 0 x ) of plants at different growth stages; T A E x is evaluated directly by T E 0 x , and its value is determined jointly by T E 0 x and P x ; C A w x is the ratio of soil water availability to precipitation; C A w x is the soil available water content of grid x (mm); Z is the seasonal constant.

2.4.2. Soil Conservation and Soil Loss

The Sediment Delivery Ratio model of InVEST is used to calculate soil loss and soil conservation based on the soil loss equation. The original soil loss equation is added to the model to intercept the sediment and grid upstream sediment, and on this basis, the soil conservation amount is calculated [68].
R U S L E n = R n × K n × L S n × C n × P n
U S L E n = R n × K n × L S n
S E D R E T n = U S L E n R U S L E n + S E D n
where R U S L E n is the actual amount of soil erosion, U S L E n is the potential soil loss under no vegetation coverage, S E D R E T n is the soil conservation, R n is the rainfall erosivity, K n is the soil erodibility factor, L S n is the slope length–gradient factor, C n is the crop management factor, P n is the support practice factor, and S E D n is the amount of sediment and sediment intercepted upstream.

2.4.3. Carbon Stocks

The Carbon stocks model of InVEST quantifies carbon storage and sequestration based on four carbon pools: Aboveground biomass, underground biomass, dead organic matter and soil organic matter [69,70]. Carbon storage in vegetation is estimated by multiplying the evegetation carbon density by the area of each LUCC based on local research results [71].
The carbon density was cropland (30.96 tC·ha−1), orchard (30.96 tC·ha−1), forested land (150.40 tC·ha−1), scrubland (118.42 tC·ha−1), grassland (96.68 tC·ha−1), bare land (10.00 tC·ha−1), water bodies (0 tC·ha−1), wasteland (96.68 tC·ha−1) and built-up land (0 tC·ha−1), respectively, according to the related reference [49].
C t o t a l = C a b o v e + C b e l o w + C s o i l + C d e a d
where C t o t a l is the total carbon stocks (t/hm2), C a b o v e is the aboveground carbon stocks t/hm2), C b e l o w is the underground carbon stocks (t/hm2), C s o i l is the soil carbon stocks (t/hm2), and C d e a d is the dead organic carbon stocks (t/hm2).

3. Results

3.1. Transfer Characteristics of LUCC and Landscape Patterns

3.1.1. Temporal and Spatial Transfer Characteristics of LUCC

During the study period, the total BUL increased by 21.11 km2 (32.40%). This increased BUL was mainly the cost of CL (15.49 km2), OC (2.58 km2) and WL (2.42 km2). During the study period, the CL decreased rapidly, with a total decrease of 66.97% (20.13 km2). Before 2000, the reduced area of CL was mainly transformed into BUL and WL. After 2000, CL was mainly transformed into OC and BUL. Before 2000, the transferred CL area was 18.80 km2, of which 53.14% (9.99 km2) was converted to BUL and 43.24% (8.15 km2) was converted to WL. From 2000 to 2010, the total area of CL transferred out was 14.94 km2. Of the reduced CL, 32.20% (4.81 km2) was transformed into BUL and 42.24% (6.31 km2) into OC. In general, in the past 30 years, the urbanization development of the study area was mainly manifested in the transformation of agricultural land (CL and OC) into BUL. CL accounts for 73.38% of the total land transformed into BUL and OC accounts for 12.22%. In addition, the WL also decreased significantly by 80.87%, mainly transformed into BUL, FL, SL and OC, of which the area converted to FL accounted for 2.58%, BUL accounted for 2.36%, OC accounted for 10.85% and SL accounted for 82.03%. The transfer into SL comprised 224.94 km2 and the transfer out SL was 74.82 km2. The area transferred into FL comprised 77.43 km2 and the area of FL transferred out was 120.46 km2 (Figure 5).

3.1.2. Evolution Characteristics of Landscape Patterns of LUCC

During the study period, BUL became more fragmented, and landscape heterogeneity increased. The landscape heterogeneity of OC first increased and then decreased. Before 2000, the landscape patterns of OC showed a trend of fragmentation, and after 2000, landscape connectivity increased. CL and FL became fragmented, reducing the connectivity. There were no obvious changes in WB, GL and BL during the study period. The PD of BUL increased from 1.35 patches/km2 to 1.59 patches/km2. The NP and MPS of CL decreased from 1985 to 2014, reaching 340 patches (8.84 × 10−2 km2) and 224 patches (4.43 × 10−2 km2), respectively (Figure 3). The NP of WL increased first and then decreased, from 2379 in 1985 to 2873 in 2000, and then decreased to 1639 in 2014. The MPS and PD of WL decreased, and the COHESION decreased. The NP of SL increased from 1491 patches in 1985 to 2189 patches in 2014. The NP of SL decreased first and then increased, from 2409 patches in 1985 to 2128 patches in 2000, and then increased to 2460 patches in 2014 (Figure 6).

3.2. Spatial Distribution of ES and Eco-Environmental Suitability of LUCC

3.2.1. ES of Mentougou District in 2014

The total amounts of SC, SLO, WY and CS in 2014 were 0.50 × 108 t, 0.66 × 106 t, 2.02 × 108 m3 and 1713.54 × 104 t, respectively. In terms of spatial characteristics, CS, SC and SLO were mainly distributed in some hilly areas in Mentougou District, followed by hilly areas with relatively fragmented LUCC along the river. This part of LUCC was mainly composed of CL, OC and WL. The smallest area was mainly concentrated in the southeast Mentougou District plain area and part of the BL. In general, the CS and WY of different LUCC vary greatly. FL and SL had the strongest ES capacity, with a total CS of 1.73 × 108 m3 (85%), SC of 0.43 × 108 t (80%), SLO of 0.46 × 106 t (70%), and WY of 0.16 × 108 t (96.48%) (Figure 7).

3.2.2. Spatial Distribution and Quantitative Structure of Eco-Environment

According to geomorphology, the study area is mainly divided into three types: Plain area (≤200 m), hilly area (200–500 m) and mountainous area (≥500 m), covering 85.00 km2 (5.84%), 333.15 km2 (22.90%) and 85.00 km2 (71.26%), respectively. The main soil types were cinnamon and brown, accounting for 1136.32 km2 and 298.93 km2 or 78.10% and 20.55%, respectively. The area of bare land and water area was 14.48 km2, covering an area of 1.00% (Figure 8).
Mentougou district has a total area of 107.59 km2 with a slope between 0° and 5°, accounting for 7.40% of the whole region. An area of 386.62 km2 has a slope between 5° and 15°. An area of 475.07 km2 has a slope between 15° and 25°, accounting for 32.65%. The area with a slope of >= 25° was 485.60 km2, accounting for 33.38% of the county area. The sunny slope accounts for 25.01% of the total area of the district, and the shady slope for 38.62%. Soil pH was nearly neutral with values between 6.5 and 7.5, covering a total of 1027.39 km2, accounting for 70.62% of the whole region. The acidic soil (pH 5.5~6.5) covers an area of 395.5 km2, covering 27.18%. The area of alkaline soil (pH 7.5~8.5) was 31.99 km2, covering 2.20%.
The area with SLO less than or equal to 10 t/hm2 had a total of 1405.77 km2 (96.63%). The area with SC less than or equal to 200 t/hm2 was 140.06 km2, accounting for 9.63% of the total area of the whole region. The area of SC between 200 t/hm2 and 600 t/hm2 was 1281.19 km2, covering 88.06%. The area with WY less than or equal to 500 m3/hm2, 500~1000 m3/hm2, 1000~1500 m3/hm2 and greater than or equal to 1500 m3/hm2 were 84.63 km2 (5.82%), 361.03 km2 (24.82%), 375.98 km2 (25.84%) and 633.24 km2 (43.53%), respectively. The evaluation was rated on a scale of 0 to 2, 2 to 3 and 3 to 4, which were suitable, relatively suitable and unsuitable, respectively (Table A7 and Figure A1).

3.2.3. Spatial Distribution Characteristics of Suitable LUCC

The results of the land suitability grade evaluation showed that the area suitable for CL was 514.36 km2 (35.35%), mainly distributed in the gully and hilly regions (Figure A1 and Table A10). The area unsuitable for CL was 940.52 km2, accounting for 64.65% of the total area of the whole region (Figure 9 and Table 2). The area suitable for OC was 1130.64 km2 (77.71%), mainly distributed in the eastern plain and part of the central hilly region, mainly in most plain and hilly areas, and some mountainous areas. The area of suitable FL and relatively suitable FL was 1350.94 km2 (92.86%). The area of suitable GL and relatively suitable GL was 1350.94 km2 (92.85%). The area unsuitable for SL was 103.94 km2 (7.14%).

3.3. Comparative Results of LUCC and ES under the two Scenarios

3.3.1. LUCC Prediction and Optimization Results in 2030

The comparison between Scenario_1 in 2030 and the situation in 2014 showed that BUL will mainly expand into adjacent areas in the future. The difference was that the patch size and expansion rate of each BUL were different. The spread degree of BUL in plain areas was higher than that in the hilly areas, and the spread degree of BUL was the slowest in mountainous areas (Figure 10). Compared with the LUCC situation in 2014, all the WL in the Scenario_2 (Optimization of 2030) was changed into other LUCC with higher comprehensive benefits, and the WL in the prediction results increased by 1.22 km2 (5.07%). The area of FL and WB increased by 22.91 km2 (4.12%) and 5.88 km2 (35.74%), respectively, while the area of CL and SL had little change. The FL, OC and WB decreased by 53.47 km2 (9.61%), 11.78 km2 (22.14%) and 0.28 km2 (1.70%), respectively, and the BUL increased in both scenarios. In terms of spatial distribution, Cl and OC migrate from an unsuitable area to a suitable area (Figure 10 and Table 3).

3.3.2. Comparison of the ES of Present, Prediction and Optimization

The predicted and optimized land use map of 2030 was used to predict the model, together with the meteorological data and soil data of 2014 (precipitation, potential evapotranspiration, water available to vegetation, soil depth and rainfall erosivity, etc.). The predicted and optimized land use ES was calculated and compared with the ES in 2014.
From the modelled data, in 2030. all ecosystem functions were declining. WY, CS and SC decreased by 0.16 × 108 m3 (7.92%), 6.15 × 104 t (0.36%) and 0.08 × 108 t (16.00%), while SLO increased by 0.08 × 106 t (12.12%) (Figure 11 and Table 4). After land-use optimization, all ES except CS would be improved. WY would be increased by 31.19% (0.63 × 108 m3) and SC would be increased by 12.00% (0.06 × 108 t). SLO would be decreased by 9.09% (0.06 × 106 t). However, CS would be decreased slightly, by 0.35% (Figure 11 and Table 4). In general, at the optimized land use, WY and SC capacity would be improved, and ES would be improved (Figure 11 and Table 4).

4. Discussion

4.1. Impacts of Human Activities and Policies on LUCC

The landscape patterns of Mentougou District have undergone great changes. Since the promulgation of the Land Management Law in 1986, the market economy system (an economy in which social resources were allocated through the market) have gradually occupied a dominant position [72]. Mentougou District lost more than half of its CL in the nearly 30 years from 1985 to 2014. People tend to choose more fertile land for farming and thus the productivity of the land decreased distinctly during this period [73,74,75]. Meanwhile, BUL increased steadily during the study period. Such changes are common in rapidly urbanizing suburbs, and the Los Angeles and California suburbs of Orange County and San Bernardino County have witnessed a dramatic loss of prime farmland and orange groves to suburban development [76]. The study area implemented the policy “Grain to Green” in 1999, prompting further losses of cropland [49]. Future economic development should take into account the unique natural conditions of mountainous areas and their landscape patterns. The government should adopt development and land management policies so that new urban development can be sustainable.

4.2. The Impact of Landscape Patterns Changes

The WL in Mentougou District was gradually fragmented, especially after 2005, and landscape heterogeneity increased rapidly. The BUL, CL and BL showed a trend of fragmentation in this period, and the degree of landscape aggregation decreased. The plain sub-region of Mentougou District accounted for less than 10% of the region’s area but contained more than 60% of the population [50]. People preferred to build houses in the plains, near the existing residential development [17]. Population pressure in the plains decreased the area of cropland as it was converted to built-up land. Over time, small patches of FL increased, and SL formed large patches, which became the main LUCC in the study area. The landscape heterogeneity of OC and WB increased first and then decreased. Before 2000, the landscape patterns of OC were fragmented, and after 2000, landscape connectivity increased. The habitat was divided into fragmentary patches, which reduced ecosystem stability and accelerated the invasion of alien species due to the edge effect. Since 2005, invasive alien species have been found in Mentougou District, including Solanum rostratum Dunal. and Pueraria phaseoloides (Roxb.) Benth. [77]. These alien species, with their strong ability to survive, will take over the resources of other plants and animals and cause the extinction of local species [78,79].

4.3. Spatial Suitability of Various Types of LUCC in the Region

Through the classification and evaluation of Mentougou’s basic eco-environmental factors (including slope, aspect, elevation, soil type, soil pH and ES), the suitable spatial distribution positions and quantities of various LUCC in the study area were selected. According to the Chinese policy of Grain to Green, areas with a slope greater than 25° were strictly not allowed to be CL, so the suitable area for CL accounted for only one-third of the total area [80,81]. The suitable area of CL was mainly distributed in the eastern plain area and some parts of the central hilly area. These places were usually rich in water and experienced less soil loss [82]. The CL-suitable area was mostly brown soil and tidal soil, with a small slope, complete water conservancy facilities, convenient irrigation and high management level. It was a basic farmland area with high and stable output in Mentougou District. The unsuitable area of OC accounted for 22.29%. Cinnamon accounts for more than two-thirds of the study area, and most areas were suitable for FL, GL and SL. After returning CL to FL, the FL in the study area would further increase. The FL could improve the regional microclimate, reduce soil erosion and water evaporation, prevent wind and consolidate soil and improve soil physical and chemical properties. However, our study did not compare different suitability evaluation methods. Future studies should add other evaluation methods, such as machine learning and sensitivity analysis, to increase the objectivity of suitability evaluation [83,84].

4.4. Impacts of Future LUCC on ES

Based on the existing spatial transfer matrix of LUCC, the CA-Markov model was used to predict LUCC in 2030, the control constraints were set and the Scenario_2 spatial distribution map of LUCC in 2030 was obtained. Compared with the LUCC in 2014, Scenario_2 was more in line with the protection needs of the local government for CL. From the goal of common development of ecological benefits (ecological space) and economic benefits (built-up land space), Science_2 is more suitable than Science_1. Linking the InVEST and CA-Markov models has promising applications for guiding ecological engineering. First, the CA-Markov model was applied to simulate LUCC by combing eco-environmetat, which provides a method for simulating LUCC in other areas [17,26]. More importantly, the predicted land cover map aided in the identification of areas where ecological space is more likely to occur, which has great significance for site selection for investment in ecological engineering. Second, the InVEST model can quantify a number of ecosystem services [44,46].
Against the background of the south-to-north water transfer and the government's encouragement of ecological restoration of mountains, rivers, forests, farmland and lakes, Scenario 2 would provide important guidance for planning. In addition, the OC of Scenario_1 was reduced by 11.78 km2 (22.14%) and the OC of Scenario_2 remained unchanged. In Scenario_1, SL increased, while in Scenario_2, SL remained unchanged. Because the economic and ecological benefits of the GL in the study area were not high, the area was reduced in both Scenario_1 and Scenario_2 (Table 3). In terms of spatial distribution, Scenario_2 was more conducive to the future development of the research area. ES (WY, CS, SC and SL) was evaluated using environmental and climatic factors in 2014, and each index proved that Scenario_2 had higher ES capacity.

5. Conclusions

From the perspective of landscape patterns, during the past 30 years (1985–2014), LUCC in Mentougou District had undergone drastic changes, with rapid urbanization. BUL increased by one-third, CL decreased by more than half, and WL decreased by more than two-thirds. At the same time, the landscape patterns of various LUCC also changed, generally showing a trend of gradual fragmentation, decreasing landscape connectivity and increasing heterogeneity. In terms of ES and spatial suitability, the total CS, SC, SLO and WY in the study area in 2014 were 1713.54 × 104 t, 0.50 × 108 t, 0.66 × 106 t and 2.02 × 108 m3. The area suitable for CL and OC in the study area accounted for one-third and two-thirds of the whole study area, and the area suitable for FL, GL and SL accounted for more than 90%. In 2030, the areas of FL, WB, OC and CL in the study area would decrease in Scenario_1, and CL would especially decrease by about half. In Scenario_2, FL and WB would increase by about 40% in total, and BUL also would develop well. In scenario_2, ES would be also better than scenario_1, WY and SC would increase by more than one-third and SLO would decrease. According to the land-use configuration of Scenario_2, the vegetation would grow under more suitable environmental conditions, and the spatial distribution of land use would be more reasonable. This study is helpful for policymakers, planners and landscape designers to determine urban land use schemes and promote the rational allocation of urban land.

Author Contributions

Conceptualization, Y.Y. and C.Z.; methodology, H.K.; software, Y.Y.; validation, Y.Z.; formal analysis, Y.Y.; investigation, Y.Y.; resources, C.Z., H.S. and Y.Y.; data curation, J.Z.; writing—original draft preparation, Y.Y.; writing—review and editing, Y.Y. and H.K.; visualization, Y.Y.; supervision, H.K.; project administration, C.Z.; funding acquisition, H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (31770746), the Yangfan Special Project of Shanghai Qimingxing Program (22YF1444000), the Youth initiation project of the Shanghai Academy of landscape planning (KT00262) and the Special Project of the Shanghai Municipal Economy and Information Technology Commission (201901024).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Abbreviations of land use and land coverage changes (LUCC), landscape metrics (LM) and ecosystem service (ES).
Table A1. Abbreviations of land use and land coverage changes (LUCC), landscape metrics (LM) and ecosystem service (ES).
TypesAbbreviationContent
Land use and land coverage changes (LUCC)BULBuilt-up land
OCOrchard
WBWater bodies
CLCropland
GLGrassland
WLWasteland
BLBareland
FLForestland
SLShrubland
Landscape metrics (LM)NPNumber patches
PDPatch density
MPSMean Patch Area
COHESIONPatch Cohesion Index
Ecosystem service (ES)WYWater yield
SLOSoil Loss
SCSoil conservation
CSCarbon stocks
Table A2. Land use transfer matrix in Mentougou District from 1985 to 1990 (km2).
Table A2. Land use transfer matrix in Mentougou District from 1985 to 1990 (km2).
BULOCWBCLGLWLBLFLSL
BUL65.16
OC1.3744.99 8.41 2.4
WB 18.2 0.45
CL3.07 23.01 3.98
GL 17.07
WL0.43 2.067.1882.021.9332.09
BL 0.98 1.11
FL 599.58
SL 0.8538.59
NOTE: The data in each row add up to the total area of land use in the same category in 1985, and the data in each column equal the total area of land use in the same category in 1990. The data in each cell represent the area of the row land-use type transferred to the column land-use type from 1985 to 1990. The total area of the study area is 64,500 km2.
Table A3. Land use transfer matrix in Mentougou District from 1990 to 1995 (km2).
Table A3. Land use transfer matrix in Mentougou District from 1990 to 1995 (km2).
BULOCWBCLGLWLBLFLSL
BUL70.03
OC0.8334.13 0.87 7.86 1.3
WB 18.2
CL0.4 20.52 4.15
GL 18.01 0.286.94
WL0.145.29 2.460.6268.6 8.369.39
BL 0.45 2.59
FL 0.54 600.733.63
SL 1.84 3.76532.99
NOTE: The data in each row add up to the total area of land use in the same category in 1990, and the data in each column equal the total area of land use in the same category in 1995. The data in each cell represent the area of the row land-use type transferred to column land-use type from 1990 to 1995. The total area of the study area is 64,500 km2.
Table A4. Land use transfer matrix in Mentougou District from 1995 to 2000 (km2).
Table A4. Land use transfer matrix in Mentougou District from 1995 to 2000 (km2).
BULOCWBCLGLWLBLFLSL
BUL71.39
OC 32.01 2.75 4.66
WB 20.04
CL6.520.17 17.17 0.51
GL 12.596.49
WL1.07 5.14 68.95 1.084.38
BL 0.06 2.53
FL 563.6549.46
SL 2.29 0.01 9.92572.04
NOTE: The data in each row add up to the total area of land use in the same category in 1995, and the data in each column equal the total area of land use in the same category in 2000. The data in each cell represent the area of the row land-use type transferred to column land-use type from 1995 to 2000. The total area of the study area is 64,500 km2.
Table A5. Land use transfer matrix in Mentougou District from 2000 to 2005.
Table A5. Land use transfer matrix in Mentougou District from 2000 to 2005.
BULOCWBCLGLWLBLFLSL
BUL78.98
OC 32.18
WB 18.76 3.57
CL2.494.02 15.8
GL 12.65
WL0.486.31 0.3958.13 12.89
BL 2.21 0.32
FL 554.0120.63
SL 0.91630.15
NOTE: The data in each row add up to the total area of land use in the same category in 2000, and the data in each column equal the total area of land use in the same category in 2005. The data in each cell represent the area of the row land-use type transferred to column land-use type from 2000 to 2005. The total area of the study area is 64,500 km2.
Table A6. Land use transfer matrix in Mentougou District from 2005 to 2010.
Table A6. Land use transfer matrix in Mentougou District from 2005 to 2010.
BULOCWBCLGLWLBLFLSL
BUL81.96
OC 38.03 4.48
WB 16.341.051.37
CL2.322.29 7.37 3.82
GL 15.78 0.82
WL0.3916.77 0.77 36.11 4.09
BL 0.85 1.36
FL 502.3152.61
SL 51.01612.98
NOTE: The data in each row add up to the total area of land use in the same category in 2005, and the data in each column equal the total area of land use in the same category in 2010. The data in each cell represent the area of the row land-use type transferred to column land-use type from 2005 to 2010. The total area of the study area is 64,500 km2.
Table A7. Land use transfer matrix in Mentougou District from 2010 to 2015.
Table A7. Land use transfer matrix in Mentougou District from 2010 to 2015.
BULOCWBCLGLWLBLFLSL
BUL84.67
OC 53.19 0.74 3.16
WB 16.17 0.17
CL1.6 7.59
GL 16.89 1.11
WL 0.020.281.60.6623.88 9.67
BL 0.14 1.22
FL 553.32
SL 0.29 3.22675.29
NOTE: The data in each row add up to the total area of land use in the same category in 2010, and the data in each column equal the total area of land use in the same category in 2015. The data in each cell represent the area of the row land-use type transferred to column land-use type from 2010 to 2015. The total area of the study area is 64,500 km2.
Table A8. Land use transfer matrix in Mentougou District from 1985 to 2015.
Table A8. Land use transfer matrix in Mentougou District from 1985 to 2015.
BULOCWBCLGLWLBLFLSL
BUL65.16
OC2.5839.59 1.77 13.23
WB0.62 16.420.050.940.62
CL15.492.51 7.98 0.12 3.96
GL 13.79 3.28
WL2.4211.11 0.131.7323.30.372.6484.01
BL 1.24 0.85
FL 479.12120.46
SL 0.03 74.79464.57
NOTE: The data in each row add up to the total area of land use in the same category in 1985, and the data in each column equal the total area of land use in the same category in 2015. The data in each cell represent the area of the row land-use type transferred to column land-use type from 1985 to 2015. The total area of the study area is 64,500 km2.
Table A9. The LUCC quantitative structure obtained using grey-linear model.
Table A9. The LUCC quantitative structure obtained using grey-linear model.
Scenario_2CLOCFLSHLGLWBBULWLBL
Quantitative of LUCC9.9353.21579.45689.230.1922.3399.320.001.22
Table A10. The structure table of evaluating indicators in Mentougou District.
Table A10. The structure table of evaluating indicators in Mentougou District.
FactorsTypesArea (km2)Percentage (%)Suitability Evaluation
CroplandOrchardForestlandGrasslandShrublandLegend
Slope(°)≤5107.597.444444High
5~15386.6226.5733444
15~25475.0732.6502334
≥25485.633.3800202
Aspect(°)6.480.4544444Low
0°–45°194.2113.3521434
45°–135°392.6426.9933343
135°–225°363.8225.0144444
225°–315°324.1322.2833343
315°–360°173.6111.9321434
Elevation(m)≤200855.8444444
200~500333.1522.933444
≥5001036.7371.2612444
Soil typeSoilless area14.48100000
Cinnamon soil1136.3278.133444
Brown soil298.9320.5544444
Moisture soil0.880.0644444
Meadow soil4.270.2903242
Soil Conservation(t/hm2)≤200140.069.6333444
200~400738.1950.7422444
400~60054337.3211434
≥60033.642.3100424
Water Yield(m3/hm2)≤50084.635.8200021
500~1000361.0324.8212232
1000~1500375.9825.8443333
≥1500633.2443.5344444
Soil loss(t/hm2)≤5981.6867.4844444
5~10424.0929.1533434
10~1522.911.5712424
≥1526.21.800414
pH of soil5.0~6.5395.527.1800333
6.5~7.51027.3970.6244444
7.5~8.531.992.222333
Figure A1. Map of land use type suitability evaluation in Mentougou District; (ae) the spatial distribution of suitability for cropland, orchard, forestland, grassland and shrubland; (f) areas that are not included in model prediction.
Figure A1. Map of land use type suitability evaluation in Mentougou District; (ae) the spatial distribution of suitability for cropland, orchard, forestland, grassland and shrubland; (f) areas that are not included in model prediction.
Ijerph 19 02432 g0a1

References

  1. Worm, B.; Barbier, E.B.; Beaumont, N.; Duffy, J.E.; Folke, C.; Halpern, B.S.; Jackson, J.B.C.; Lotze, H.K.; Micheli, F.; Palumbi, S.R.; et al. Impacts of Biodiversity Loss on Ocean Ecosystem Services. Science 2006, 314, 787–790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fisher, B.; Turner, R.K.; Morling, P. Defining and classifying ecosystem services for decision making. Ecol. Econ. 2009, 68, 643–653. [Google Scholar] [CrossRef] [Green Version]
  3. Yuan, Z.; Xu, J.; Wang, Y.; Yan, B. Analyzing the influence of land use/land cover change on landscape pattern and ecosystem services in the Poyang Lake Region, China. Environ. Sci. Pollut. Res. 2021, 28, 27193–27206. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, L.; Pei, S.; Liu, X.; Qiao, Q.; Liu, C. Mapping and analysing tradeoffs, synergies and losses among multiple ecosystem services across a transitional area in Beijing, China. Ecol. Indic. 2021, 123, 107329. [Google Scholar] [CrossRef]
  5. Mu, L.; Fang, L.; Dou, W.; Wang, C.; Qu, X.; Yu, Y. Urbanization-induced spatio-temporal variation of water resources utilization in northwestern China: A spatial panel model based approach. Ecol. Indic. 2021, 125, 107457. [Google Scholar] [CrossRef]
  6. Ulucak, R.; Erdogan, F.; Bostanci, S.H. A STIRPAT-based investigation on the role of economic growth, urbanization, and energy consumption in shaping a sustainable environment in the Mediterranean region. Environ. Sci. Pollut. Res. 2021, 28, 55290–55301. [Google Scholar] [CrossRef] [PubMed]
  7. National Bureau of Statistics of China. China Statistical Yearbook; China Statistics Press: Beijing, China, 2021. [Google Scholar]
  8. Wu, Z.X.; Zhang, Q.; Song, C.Q.; Zhang, F.; Zhu, X.D.; Sun, P.; Fan, K.K.; Yu, H.Q.; Shen, Z.X. Impacts of urbanization on spatio-temporal variations of temperature over the Pearl River Delta. Acta Geogr. Sin. 2019, 74, 2342–2357. [Google Scholar] [CrossRef]
  9. Zięba-Kulawik, K.; Hawryło, P.; Wężyk, P.; Matczak, P.; Przewoźna, P.; Inglot, A.; Mączka, K. Improving methods to calculate the loss of ecosystem services provided by urban trees using LiDAR and aerial orthophotos. Urban For. Urban Green. 2021, 63, 127195. [Google Scholar] [CrossRef]
  10. Lyu, R.; Clarke, K.C.; Zhang, J.; Feng, J.; Jia, X.; Li, J. Dynamics of spatial relationships among ecosystem services and their determinants: Implications for land use system reform in Northwestern China. Land Use Policy 2021, 102, 105231. [Google Scholar] [CrossRef]
  11. Elhacham, E.; Alpert, P. Temperature patterns along an arid coastline experiencing extreme and rapid urbanization, case study: Dubai. Sci. Total Environ. 2021, 784, 147168. [Google Scholar] [CrossRef] [PubMed]
  12. Millennium Ecosystem Assessment. In Ecosystems and Human Well-Being. Synthesis; Island Press: Washington DC, USA, 2005.
  13. Kim, I.; Lee, J.; Kwon, H. Participatory ecosystem service assessment to enhance environmental decision-making in a border city of South Korea. Ecosyst. Serv. 2021, 51, 101337. [Google Scholar] [CrossRef]
  14. Shah, M.; Cummings, A.R. An analysis of the influence of the human presence on the distribution of provisioning ecosystem services: A Guyana case study. Ecol. Indic. 2021, 122, 107255. [Google Scholar] [CrossRef]
  15. Dong, X.; Wang, X.; Wei, H.; Fu, B.; Wang, J.; Uriarte-Ruiz, M. Trade-offs between local farmers’ demand for ecosystem services and ecological restoration of the Loess Plateau, China. Ecosyst. Serv. 2021, 49, 101295. [Google Scholar] [CrossRef]
  16. Giefer, M.M.; An, L.; Chen, X. Normative, livelihood, and demographic influences on enrollment in a payment for ecosystem services program. Land Use Policy 2021, 108, 105525. [Google Scholar] [CrossRef]
  17. Guan, D.; Zhao, Z.; Tan, J. Dynamic simulation of land use change based on logistic-CA-Markov and WLC-CA-Markov models: A case study in three gorges reservoir area of Chongqing, China. Environ. Sci. Pollut. Res. 2019, 26, 20669–20688. [Google Scholar] [CrossRef]
  18. Wei, Z.; Liu, Y. Construction of super-resolution model of remote sensing image based on deep convolutional neural network. Comput. Commun. 2021, 178, 191–200. [Google Scholar] [CrossRef]
  19. Mamanis, G.; Vrahnakis, M.; Chouvardas, D.; Nasiakou, S.; Kleftoyanni, V. Land Use Demands for the CLUE-S Spatiotemporal Model in an Agroforestry Perspective. Land 2021, 10, 1097. [Google Scholar] [CrossRef]
  20. Hu, S.; Chen, L.; Li, L.; Zhang, T.; Yuan, L.; Cheng, L.; Wang, J.; Wen, M. Simulation of Land Use Change and Ecosystem Service Value Dynamics under Ecological Constraints in Anhui Province, China. Int. J. Environ. Res. Public Health 2020, 17, 4228. [Google Scholar] [CrossRef]
  21. Motlagh, Z.K.; Lotfi, A.; Pourmanafi, S.; Ahmadizadeh, S.; Soffianian, A. Spatial modeling of land-use change in a rapidly urbanizing landscape in central Iran: Integration of remote sensing, CA-Markov, and landscape metrics. Environ. Monit. Assess. 2020, 192, 695. [Google Scholar] [CrossRef]
  22. Huang, S.; Xi, F.; Chen, Y.; Gao, M.; Pan, X.; Ren, C. Land Use Optimization and Simulation of Low-Carbon-Oriented—A Case Study of Jinhua, China. Land 2021, 10, 1020. [Google Scholar] [CrossRef]
  23. Manca, F.; Robinson, E.; Dillon, J.F.; Boyd, K.A. Eradicating hepatitis C: Are novel screening strategies for people who inject drugs cost-effective? Int. J. Drug Policy 2020, 82, 102811. [Google Scholar] [CrossRef] [PubMed]
  24. Gao, Y.; Chen, J.; Luo, H.; Wang, H. Prediction of hydrological responses to land use change. Sci. Total Environ. 2020, 708, 134998. [Google Scholar] [CrossRef] [PubMed]
  25. Jenerette, G.D.; Wu, J. Analysis and simulation of land-use change in the central Arizona—Phoenix region, USA. Landsc. Ecol. 2001, 16, 616–626. [Google Scholar] [CrossRef]
  26. Nourqolipour, R.; Mohamed Shariff, A.R.B.; Balasundram, S.K.; Ahmad, N.B.; Sood, A.M.; Buyong, T.; Amiri, F. A GIS-based model to analyze the spatial and temporal development of oil palm land use in Kuala Langat district, Malaysia. Environ. Earth Sci. 2015, 73, 1687–1700. [Google Scholar] [CrossRef] [Green Version]
  27. Alcamo, J.; van Vuuren, D.; Ringler, C.; Cramer, W.; Masui, T.; Alder, J.; Schulze, K. Changes in Nature’s Balance Sheet: Model-based Estimates of Future Worldwide Ecosystem Services. Ecol. Soc. 2005, 10, art19. [Google Scholar] [CrossRef] [Green Version]
  28. Ehrlich, P.R.; Ehrlich, A.H. Environmental Problem Solving. Ecology 1987, 68, 2067–2068. [Google Scholar] [CrossRef]
  29. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  30. Long, H.; Liu, Y.; Hou, X.; Li, T.; Li, Y. Effects of land use transitions due to rapid urbanization on ecosystem services: Implications for urban planning in the new developing area of China. Habitat Int. 2014, 44, 536–544. [Google Scholar] [CrossRef]
  31. Raskin, P.D. Global Scenarios: Background Review for the Millennium Ecosystem Assessment. Ecosystems 2005, 8, 133–142. [Google Scholar] [CrossRef]
  32. Heydinger, J.M. Reinforcing the Ecosystem Services Perspective: The Temporal Component. Ecosystems 2016, 19, 661–673. [Google Scholar] [CrossRef]
  33. Foudi, S.; Osés-Eraso, N.; Tamayo, I. Integrated spatial flood risk assessment: The case of Zaragoza. Land Use Policy 2015, 42, 278–292. [Google Scholar] [CrossRef]
  34. FU, B.J. The integrated studies of geography: Coupling of patterns and processes. Acta Geogr. Sin. 2014, 69, 1052–1059. [Google Scholar] [CrossRef]
  35. Ożgo, M.; Urbańska, M.; Marzec, M.; Kamocki, A.; Andrzejewski, W.; Golski, J.; Lewandowski, K.; Geist, J. Lake-stream transition zones support hotspots of freshwater ecosystem services: Evidence from a 35-year study on unionid mussels. Sci. Total Environ. 2021, 774, 145114. [Google Scholar] [CrossRef] [PubMed]
  36. Li, R.; Shi, Y.; Feng, C.; Guo, L. The spatial relationship between ecosystem service scarcity value and urbanization from the perspective of heterogeneity in typical arid and semiarid regions of China. Ecol. Indic. 2021, 132, 108299. [Google Scholar] [CrossRef]
  37. Lee, L.S.H.; Zhang, H.; Jim, C.Y. Serviceable tree volume: An alternative tool to assess ecosystem services provided by ornamental trees in urban forests. Urban For. Urban Green. 2021, 59, 127003. [Google Scholar] [CrossRef]
  38. Grunhut, J.H.; Wade, G.A.; Marcolino, W.L.F.; Petit, V. A MiMeS analysis of the magnetic field and circumstellar environment of the weak-wind O9 sub-giant star HD 57682. Proc. Int. Astron. Union 2010, 6, 188–189. [Google Scholar] [CrossRef] [Green Version]
  39. Boumans, R.; Costanza, R.; Farley, J.; Wilson, M.A.; Portela, R.; Rotmans, J.; Villa, F.; Grasso, M. Modeling the dynamics of the integrated earth system and the value of global ecosystem services using the GUMBO model. Ecol. Econ. 2002, 41, 529–560. [Google Scholar] [CrossRef]
  40. Bagstad, K.J.; Johnson, G.W.; Voigt, B.; Villa, F. Spatial dynamics of ecosystem service flows: A comprehensive approach to quantifying actual services. Ecosyst. Serv. 2013, 4, 117–125. [Google Scholar] [CrossRef]
  41. Zhu, Q.; Jiang, H.; Liu, J.; Wei, X.; Peng, C.; Fang, X.; Liu, S.; Zhou, G.; Yu, S.; Ju, W. Evaluating the spatiotemporal variations of water budget across China over 1951-2006 using IBIS model. Hydrol. Process. 2010, 24, 429–445. [Google Scholar] [CrossRef]
  42. Carver, A.D.; Unger, D.R.; Parks, C.L. Modeling Energy Savings from Urban Shade Trees: An Assessment of the CITYgreen® Energy Conservation Module. Environ. Manag. 2004, 34, 650–655. [Google Scholar] [CrossRef]
  43. Lahiji, R.N.; Dinan, N.M.; Liaghati, H.; Ghaffarzadeh, H.; Vafaeinejad, A. Scenario-based estimation of catchment carbon storage: Linking multi-objective land allocation with InVEST model in a mixed agriculture-forest landscape. Front. Earth Sci. 2020, 14, 637–646. [Google Scholar] [CrossRef]
  44. Caro, C.; Marques, J.C.; Cunha, P.P.; Teixeira, Z. Ecosystem services as a resilience descriptor in habitat risk assessment using the InVEST model. Ecol. Indic. 2020, 115, 106426. [Google Scholar] [CrossRef]
  45. Yang, D.; Liu, W.; Tang, L.; Chen, L.; Li, X.; Xu, X. Estimation of water provision service for monsoon catchments of South China: Applicability of the InVEST model. Landsc. Urban Plan. 2019, 182, 133–143. [Google Scholar] [CrossRef]
  46. Cong, W.; Sun, X.; Guo, H.; Shan, R. Comparison of the SWAT and InVEST models to determine hydrological ecosystem service spatial patterns, priorities and trade-offs in a complex basin. Ecol. Indic. 2020, 112, 106089. [Google Scholar] [CrossRef]
  47. Sun, X.; Yang, P.; Tao, Y.; Bian, H. Improving ecosystem services supply provides insights for sustainable landscape planning: A case study in Beijing, China. Sci. Total Environ. 2022, 802, 149849. [Google Scholar] [CrossRef]
  48. Redhead, J.W.; Stratford, C.; Sharps, K.; Jones, L.; Ziv, G.; Clarke, D.; Oliver, T.H.; Bullock, J.M. Empirical validation of the InVEST water yield ecosystem service model at a national scale. Sci. Total Environ. 2016, 569–570, 1418–1426. [Google Scholar] [CrossRef] [Green Version]
  49. Yi, Y.; Shi, M.; Liu, C.; Wang, B.; Kang, H.; Hu, X. Changes of Ecosystem Services and Landscape Patterns in Mountainous Areas: A Case Study in the Mentougou District in Beijing. Sustainability 2018, 10, 3689. [Google Scholar] [CrossRef] [Green Version]
  50. Yi, Y.; Zhao, Y.; Ding, G.; Gao, G.; Shi, M.; Cao, Y. Effects of Urbanization on Landscape Patterns in a Mountainous Area: A Case Study in the Mentougou District, Beijing, China. Sustainability 2016, 8, 1190. [Google Scholar] [CrossRef] [Green Version]
  51. Li, F.; Liu, X.; Zhao, D.; Wang, B.; Jin, J.; Hu, D. Evaluating and modeling ecosystem service loss of coal mining: A case study of Mentougou district of Beijing, China. Ecol. Complex. 2011, 8, 139–143. [Google Scholar] [CrossRef]
  52. Liu, L.; Qin, M.; Tian, N.; Zhou, C.; Wang, D.; Basinger, J.F.; Xue, J. Belowground rhizomes and roots in waterlogged paleosols: Examples from the Middle Jurassic of Beijing, China. Geobios 2018, 51, 419–433. [Google Scholar] [CrossRef]
  53. Sheng, W.; Zhen, L.; Xie, G.; Xiao, Y. Determining eco-compensation standards based on the ecosystem services value of the mountain ecological forests in Beijing, China. Ecosyst. Serv. 2017, 26, 422–430. [Google Scholar] [CrossRef]
  54. Mentougou District Local Chronicles Compilation Committee of Beijing Municipality. Beijing Mentougou Yearbook; Chinese Communist Party History Publishing House: Beijing, China, 2015. [Google Scholar]
  55. Mentougou District People's Government of Beijing Municipality. Overall Land Use Planning of Mentougou District (2006-2020). 2010. Available online: http://www.mnr.gov.cn/gk/ghjh/201811/t20181101_2324823.html (accessed on 7 October 2021).
  56. People’s Government of Beijing Municipality. In Beijing urban master plan (2004-2020); 2007. Available online: http://ghzrzyw.beijing.gov.cn/zhengwuxinxi/zxzt/bjcsztgh2004/202201/t20220110_2587452.html (accessed on 7 October 2021).
  57. United States Geological Survey. Remote Sensing Images. Available online: https://www.usgs.gov/ (accessed on 19 October 2021).
  58. Geospatial Data Cloud. Remote Sensing Images. Available online: https://www.gscloud.cn/ (accessed on 11 October 2021).
  59. Chinese Academy of Sciences. National Earth System Science Data Sharing Infrastructure. Available online: http://www.geodata.cn/ (accessed on 27 October 2021).
  60. Jia, J.L.; Zhang, Y.; Wang, C.; Li, D.; Liu, B.W.; Liu, Y.; Zhao, L.; Yang, S.Q. Soil organic pollution characteristics and microbial properties in coal mining areas of mentougou. Huanjing Kexue/Environ. Sci. 2011, 32, 875–879. [Google Scholar]
  61. Krzywinski, M.I.; Schein, J.E.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Yi, Y.; Zhang, C.; Zhang, G.; Xing, L.; Zhong, Q.; Liu, J.; Lin, Y.; Zheng, X.; Yang, N.; Sun, H.; et al. Effects of Urbanization on Landscape Patterns in the Middle Reaches of the Yangtze River Region. Land 2021, 10, 1025. [Google Scholar] [CrossRef]
  63. Wu, H.; Hu, X.; Sun, S.; Dai, J.; Ye, S.; Du, C.; Chen, H.; Yu, G.; Zhou, L.; Chen, J. Effect of increasing of water level during the middle of dry season on landscape pattern of the two largest freshwater lakes of China. Ecol. Indic. 2020, 113, 106283. [Google Scholar] [CrossRef]
  64. Yi, Y.; Shi, M.; Liu, C.; Kang, H.; Wang, B. On Landscape Patterns in Typical Mountainous Counties Middle Reaches of the Yangtze River in China. Int. J. Environ. Res. Public Health 2021, 18, 4000. [Google Scholar] [CrossRef]
  65. Feng, X.; Xiu, C.; Bai, L.; Zhong, Y.; Wei, Y. Comprehensive evaluation of urban resilience based on the perspective of landscape pattern: A case study of Shenyang city. Cities 2020, 104, 102722. [Google Scholar] [CrossRef]
  66. Kurttila, M.; Pesonen, M.; Kangas, J.; Kajanus, M. Utilizing the analytic hierarchy process (AHP) in SWOT analysis—A hybrid method and its application to a forest-certification case. For. Policy Econ. 2000, 1, 41–52. [Google Scholar] [CrossRef]
  67. Adhikari, S.; Southworth, J. Simulating Forest Cover Changes of Bannerghatta National Park Based on a CA-Markov Model: A Remote Sensing Approach. Remote Sens. 2012, 4, 3215–3243. [Google Scholar] [CrossRef] [Green Version]
  68. Etemadi, H.; Smoak, J.M.; Karami, J. Land use change assessment in coastal mangrove forests of Iran utilizing satellite imagery and CA–Markov algorithms to monitor and predict future change. Environ. Earth Sci. 2018, 77, 208. [Google Scholar] [CrossRef]
  69. Guo, H.W.; Sun, X.Y.; Lian, L.S.; Zhang, D.Z.; Xu, Y. Response of water yield function of ecosystem to land use change in Nansi Lake Basin based on CLUE-S model and InVEST-model. Chin. J. Appl. Ecol. 2016, 27, 2899–2906. [Google Scholar] [CrossRef]
  70. Ma, S.; Wang, L.J.; Zhu, D.; Zhang, J. Spatiotemporal changes in ecosystem services in the conservation priorities of the southern hill and mountain belt, China. Ecol. Indic. 2021, 122, 107225. [Google Scholar] [CrossRef]
  71. Zhang, C.; Ju, W.; Chen, J.M.; Zan, M.; Li, D.; Zhou, Y.; Wang, X. China’s forest biomass carbon sink based on seven inventories from 1973 to 2008. Clim. Change. 2013, 118, 933–948. [Google Scholar] [CrossRef]
  72. Edmondson, J.L.; Davies, Z.G.; McCormack, S.A.; Gaston, K.J.; Leake, J.R. Land-cover effects on soil organic carbon stocks in a European city. Sci. Total Environ. 2014, 472, 444–453. [Google Scholar] [CrossRef] [Green Version]
  73. Feng, Z.; Jin, X.; Chen, T.; Wu, J. Understanding trade-offs and synergies of ecosystem services to support the decision-making in the Beijing–Tianjin–Hebei region. Land Use Policy 2021, 106, 105446. [Google Scholar] [CrossRef]
  74. Jin, J.; Wang, R.S.; Li, F.; Huang, J.L.; Zhou, C.B.; Zhang, H.T.; Yang, W.R. Conjugate ecological restoration approach with a case study in Mentougou district, Beijing. Ecol. Complex. 2011, 8, 161–170. [Google Scholar] [CrossRef]
  75. Du, J.; Thill, J.C.; Peiser, R.B.; Feng, C. Urban land market and land-use changes in post-reform China: A case study of Beijing. Landsc. Urban Plan. 2014, 124, 118–128. [Google Scholar] [CrossRef]
  76. Du, J.; Thill, J.C.; Peiser, R.B. Land pricing and its impact on land use efficiency in post-land-reform China: A case study of Beijing. Cities 2016, 50, 68–74. [Google Scholar] [CrossRef]
  77. Zhou, D.C.; Zhao, S.Q.; Zhu, C. The Impact of the Grain for Green Project on the Land Use/Cover Change in the Northern Farming-pastoral Ecotone, China—A Case Study of Kezuohouqi County. Sci. Geogr. Sin. 2012, 32, 442–449. [Google Scholar] [CrossRef]
  78. Historic Land Use Conversion 1984–2010. Available online: http://www.conservation.ca.gov/dlrp/fmmp/products/Pages/DownloadGISdata.aspx (accessed on 7 October 2021).
  79. Duan, M.; Liu, Y.; Yu, Z.; Li, L.; Wang, C.; Axmacher, J.C. Environmental factors acting at multiple scales determine assemblages of insects and plants in agricultural mountain landscapes of northern China. Agric. Ecosyst. Environ. 2016, 224, 86–94. [Google Scholar] [CrossRef]
  80. Qin, D.T.; Cai, B.F. Risk Assessment of Biological Invasion in Beijing. Environ. Prot. 2004, 44–47. [Google Scholar] [CrossRef]
  81. Wang, H.H.; Liu, J.L.; Zhang, R.; Liu, J.K.; Zou, Y.Q.; Zou, D.L.; Nan, H.L.; Zhang, Z.M. The Study of Invasive Alien Plants in Beijing. J. Agric. 2014, 4, 49–52. [Google Scholar] [CrossRef]
  82. Bai, Y.L.; Wang, R.S.; Jin, J.S. Water eco-service assessment and compensation in a coal mining region—A case study in the Mentougou District in Beijing. Ecol. Complex. 2011, 8, 144–152. [Google Scholar] [CrossRef]
  83. Merghadi, A.; Yunus, A.P.; Dou, J.; Whiteley, J.; Pham, B.T.; Bui, D.T.; Avtar, R.; Abderrahmane, B. Machine learning methods for landslide susceptibility studies: A comparative overview of algorithm performance. Earth Sci. Rev. 2020, 207, 103225. [Google Scholar] [CrossRef]
  84. Abraham, M.T.; Satyam, N.; Lokesh, R.; Pradhan, B.; Alamri, A. Factors Affecting Landslide Susceptibility Mapping: Assessing the Influence of Different Machine Learning Approaches, Sampling Strategies and Data Splitting. Land 2021, 10, 989. [Google Scholar] [CrossRef]
Figure 1. Location of Mentougou district and topography.
Figure 1. Location of Mentougou district and topography.
Ijerph 19 02432 g001
Figure 2. An explanation of the circos chart. Each stripe represents the occurrence of a transformation process at two different times, and the color represents the land use types. (a) The total amount of A1 transferred in and transferred out; (b) the proportion of each land-use type transferred into A1; (c) the proportion of each land-use type in transferred from A1; (d) the total area of A1 transformed out and transformed into other types; (e) the types of land use transformed from A1; (f) the area of A1 transformed into other land-use types; (g) the types of land use transformed into A1; (h) the area of A1 transformed from other land-use types.
Figure 2. An explanation of the circos chart. Each stripe represents the occurrence of a transformation process at two different times, and the color represents the land use types. (a) The total amount of A1 transferred in and transferred out; (b) the proportion of each land-use type transferred into A1; (c) the proportion of each land-use type in transferred from A1; (d) the total area of A1 transformed out and transformed into other types; (e) the types of land use transformed from A1; (f) the area of A1 transformed into other land-use types; (g) the types of land use transformed into A1; (h) the area of A1 transformed from other land-use types.
Ijerph 19 02432 g002
Figure 3. LUCC suitability evaluation structure map of Mentougou District.
Figure 3. LUCC suitability evaluation structure map of Mentougou District.
Ijerph 19 02432 g003
Figure 4. Optimization structure of LUCC.
Figure 4. Optimization structure of LUCC.
Ijerph 19 02432 g004
Figure 5. LUCC transfer flow charts in the Mentougou District from 1985 to 2014.
Figure 5. LUCC transfer flow charts in the Mentougou District from 1985 to 2014.
Ijerph 19 02432 g005
Figure 6. Landscape pattern changes of LUCC in Mentougou District. (a) NP; (b) MPS; (c) PD; (d) CHOESION. The types of LUCC were built-up land (BUL), orchard (OC), waterbodies (WB), cropland (CL), shrubland (SL), grassland (GL), wasteland (WL), bare land (BL) and forestland (FL).
Figure 6. Landscape pattern changes of LUCC in Mentougou District. (a) NP; (b) MPS; (c) PD; (d) CHOESION. The types of LUCC were built-up land (BUL), orchard (OC), waterbodies (WB), cropland (CL), shrubland (SL), grassland (GL), wasteland (WL), bare land (BL) and forestland (FL).
Ijerph 19 02432 g006
Figure 7. Spatial distribution of soil conservation (SC) (a), soil loss (SLO) (b), water yield (WY) (c) and (d) carbon stocks (CS) by InVEST in 2014 of Mentougou District.
Figure 7. Spatial distribution of soil conservation (SC) (a), soil loss (SLO) (b), water yield (WY) (c) and (d) carbon stocks (CS) by InVEST in 2014 of Mentougou District.
Ijerph 19 02432 g007
Figure 8. Map of evaluating indicator distribution in Mentougou District; (ad) the slope, aspect, elevation and soil type; (eh) the soil loss, soil conservation, water yield and soil pH.
Figure 8. Map of evaluating indicator distribution in Mentougou District; (ad) the slope, aspect, elevation and soil type; (eh) the soil loss, soil conservation, water yield and soil pH.
Ijerph 19 02432 g008
Figure 9. Land suitability evaluation map in Mentougou District. Note: The five digits of each legend, from left to right, represent cropland, orchard, forestland, grassland and shrubland. The score of each digit is 4, 3, 2 and 1; 4 represents the suitable distribution area of a certain land use type; 3 represents the relatively suitable distribution area of a certain land use type; 2 and 1 represent the unsuitable distribution area of a certain land-use type.
Figure 9. Land suitability evaluation map in Mentougou District. Note: The five digits of each legend, from left to right, represent cropland, orchard, forestland, grassland and shrubland. The score of each digit is 4, 3, 2 and 1; 4 represents the suitable distribution area of a certain land use type; 3 represents the relatively suitable distribution area of a certain land use type; 2 and 1 represent the unsuitable distribution area of a certain land-use type.
Ijerph 19 02432 g009
Figure 10. LUCC of prediction (Scenario_1: Prediction of 2030) (a) and optimization (Scenario_2: Optimization of 2030) (b) of Mentougou in 2030.
Figure 10. LUCC of prediction (Scenario_1: Prediction of 2030) (a) and optimization (Scenario_2: Optimization of 2030) (b) of Mentougou in 2030.
Ijerph 19 02432 g010
Figure 11. Spatial distribution of ES after prediction and optimization by CA-Markov and InVEST; (ad) the ES as predicted for 2030; (eh) the ES as optimization for 2030.
Figure 11. Spatial distribution of ES after prediction and optimization by CA-Markov and InVEST; (ad) the ES as predicted for 2030; (eh) the ES as optimization for 2030.
Ijerph 19 02432 g011
Table 1. Description of landscape metrics.
Table 1. Description of landscape metrics.
TypesModel NameCalculation MethodConnotation
Area/ DensityNP N P = n i The number of patches in landscape.
PD P D = n i a i j ( 10000 ) ( 100 ) The density of landscape
MPS M P S = j = 1 n a i j n i ( 1 10000 ) It can be used to characterize landscape fragmentation.
AggregationCOHESION C O H E S I O N = 1 j = k n p i k j = 1 n p i j a i j n 1 1 n i ( 100 % ) It reflects the aggregation degree of patches in the landscape.
NOTE: a i j is the area of patch i j , p i j is the common border length of patch ij, and n i is the number of patches in the landscape of patch type (lass) i .
Table 2. The main LUCC suitability structure in Mentougou District (km2).
Table 2. The main LUCC suitability structure in Mentougou District (km2).
TypeCLOCFLGLSL
Appropriate16.4527.791341.781331.261346.02
Relatively appropriate497.911102.869.1619.684.92
Inappropriate940.52324.23103.94103.94103.94
Table 3. Area of LUCC in present situation, prediction and optimization in Mentougou (km2).
Table 3. Area of LUCC in present situation, prediction and optimization in Mentougou (km2).
LUCC2014Prediction of 2030 (Scenario_1)Optimization of 2030 (Scenario_2)
CL9.935.079.93
OC53.2141.4353.21
FL556.54503.07579.45
SL689.23759.39689.23
GL17.9813.050.19
WB16.4516.1722.33
BUL86.2790.8199.32
WL24.0525.270.00
BL1.220.611.22
Table 4. The ES of present situation, prediction and optimization.
Table 4. The ES of present situation, prediction and optimization.
ES2014Prediction for 2030 (Scenario_1)Optimization for 2030 (Scenario_2)
WY(m3)2.02 × 1081.86 × 1082.65 × 108
CS(t)1713.54 × 1041707.39 × 1041707.53 × 104
SC(t)0.50 × 1080.42 × 1080.56 × 108
SLO(t)0.66 × 1060.74 × 1060.60 × 106
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yi, Y.; Zhang, C.; Zhu, J.; Zhang, Y.; Sun, H.; Kang, H. Spatio-Temporal Evolution, Prediction and Optimization of LUCC Based on CA-Markov and InVEST Models: A Case Study of Mentougou District, Beijing. Int. J. Environ. Res. Public Health 2022, 19, 2432. https://doi.org/10.3390/ijerph19042432

AMA Style

Yi Y, Zhang C, Zhu J, Zhang Y, Sun H, Kang H. Spatio-Temporal Evolution, Prediction and Optimization of LUCC Based on CA-Markov and InVEST Models: A Case Study of Mentougou District, Beijing. International Journal of Environmental Research and Public Health. 2022; 19(4):2432. https://doi.org/10.3390/ijerph19042432

Chicago/Turabian Style

Yi, Yang, Chen Zhang, Jinqi Zhu, Yugang Zhang, Hao Sun, and Hongzhang Kang. 2022. "Spatio-Temporal Evolution, Prediction and Optimization of LUCC Based on CA-Markov and InVEST Models: A Case Study of Mentougou District, Beijing" International Journal of Environmental Research and Public Health 19, no. 4: 2432. https://doi.org/10.3390/ijerph19042432

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