Scenario-Based Analysis of Land Use Competition and Sustainable Land Development in Zhangye of the Heihe River Basin, China

Rapid economic growth has a significant impact on land use change, which would threaten the natural ecology. Zhangye city of the Heihe River Basin, China is an ecologically vulnerable region where land use changes significantly due to socioeconomic development and population increases. The study employed a computable general equilibrium of land use change (CGELUC) model to simulate land use change and then used a dynamic land system (DLS) model to spatialize land use change during 2015–2030 under three development scenarios in Zhangye city. The three development scenarios are the baseline scenario (BAU), the resource consumption scenario (RCS) and the green development scenario (GDS). We found that economic growth would lead to land demand increases in high value-added industries and decreases in low value-added industries. The cultivated land would decrease while the built-up area would increase. By 2030, the cultivated land will decrease by 8.16%, 10.89% and 4.16%, respectively, under BAU, RCS and GDS, while the built-up area will increase by 8.61%, 10.39% and 4.75%, respectively. The expansion of built-up area under RCS presents spatial characteristics of centralized distribution, while spatial characteristics of uniform discrete distributions are presented under GDS. The expansion of ecological land under GDS would be considerable, especially in the north of Sunan County and Gaotai County, and around the natural reserve of Ganzhou County. This paper provides a scientific reference for coordinating economic development and ecological protection in the rapidly developing urbanized areas in western China.


Introduction
Land use change has both social and natural attributes, and is influenced by the combination of climate change and human activities [1]. Since natural factors are relatively stable with cumulative effects over short time scales, human activities have played a dominant role in regional environmental change at the scale of decades to centuries [2]. With the rapid development of social economy, the role of human socioeconomic activities in driving the evolution of land use systems has become more obvious [3,4]. Since the reform and opening up, China has experienced rapid urbanization [5]. The rapid growth of urbanization has had a significant impact on land use change. The construction of socioeconomic development scenarios to analyze the impact of human activities on land use change has indicated the direction for land resource management in the context of sustainable development.
The Heihe River Basin is the second largest inland river basin in northwest China, located in the middle of the Hexi Corridor, China. The midstream area of Heihe River Basin accounts for 88.47% of the population size and 87.93% of the GDP of the whole basin [6]. erosion and degradation, agricultural land and forestry conservation policies and climate change policies related to land use change, using a CGE model embedding land resource accounts (Table 1). However, currently, many CGE models do not consider the land supply constraint caused by the heterogeneity of different land use types. Therefore, they cannot simulate land conversion among different uses, and the land supply response to price change would be exaggerated. The characterization of land heterogeneity is an important but challenging issue in CGE modeling. Many models are employed to simulate the spatial distribution of land use, including the CLUE-S model [42][43][44], the agent-based model [45][46][47], the cellular automata (CA) model and the dynamic of land system (DLS) model [41,48]. Zhang et al. [49] adopted the models of CLUE-S and SWAT to simulate the pollution load under land use scenarios in the upstream basin of the Miyun Reservoir in Beijing, China. Villamor et al. [50] used an agent-based model to explore the trade-offs between the provision of ecosystem services and rural income under a series of land use dynamic scenarios in lowland Sumatra. Guan et al. [51] adopted a method of CA and the Markov model to simulate the spatial and temporal distribution of land use in Saga City from 2015 to 2042. Compared with the traditional models, the DLS model constructs a spatial statistical model of land system spatial distribution and its driving factors in a spatially explicit way, which can better reflect the complex self-organization and the competition process among land use types in a land system [52,53]. Based on the principle of spatial allocation and the balance of land supply and demand, the DLS model considers the driving mechanism of natural factors and social-economic factors on land use change, and realizes the dynamic simulation of a land system on a fine grid scale. Under the integration of the CGE and DLS models, the dynamic macro-analysis of land systems is closely combined with micro-econometric analyses on the grid scale. The comprehensive consideration of regional socioeconomic development characteristics, natural conditions, policy conditions and historical trends of land use change is realized to make the simulation results more scientific and reasonable. Therefore, this paper developed a computable general equilibrium of land use change (CGELUC) model to simulate land use change in Zhangye city of Heihe River Basin, China in 2015-2030 under three development scenarios: baseline, resource consumption, and green development. The dynamic land system (DLS) model was further applied to spatialize land use change on the spatial grid scale. This paper aims to provide policy decision information on regional land resource allocation and management, regional economic sustainable development, and ecological environment protection in rapid urbanization areas in western China.

Study Area
Zhangye city of Heihe River Basin, China extends 97 • 20 -102 • 12 E and 37 • 28 -39 • 57 N, and is located at the northwest of Gansu Province in the middle section of the Hexi Corridor and the middle and upper reaches of Heihe River, the second largest inland river in China. The cultivated land covers an area of 4204 km 2 , and still about 2000 km 2 land is available for reclamation. There are 26 rivers distributed in this region, with an annual runoff of 2.66 billion m 3 . The groundwater resources are also abundant. This area is characterized by an arid continental climate, the mean annual temperature is 6 • C, and the annual sunshine is 3000 h.
Currently, Zhangye has 6 counties, including Linze county, Gaotai county, Shandan county, Minle county, Sunan county and Ganzhou district ( Figure 1). According to the statistical data, the GDP reached CNY 37.70 billion in 2017, which is nearly 3.4 times higher than that in 2005.The industrial structure also shows obvious characteristics of an increasing trend in the proportion of industry and the service industry. With the continuous optimization and adjustment of the industrial structure, the proportion of the service industry shows an increasing trend, reaching 54.75% in 2017 ( Figure 2).  The land use change in Zhangye from 1990 to 2015 indicates that forest area, grassland, water area and unused land showed a decreasing trend, while cultivated land and builtup area showed an increasing trend (Table 2). With the continuous acceleration of the urbanization process, the demand for built-up area is large. However, the areas suitable for construction are distributed in high-quality cultivated land, and the contradiction between major infrastructure construction and cultivated land protection is prominent. On the other hand, Zhangye city is located in the intersection of the Qinghai Tibet Plateau and the Loess Plateau, belonging to the western ecological security barrier. However, in recent years, with a large area of cropland having been reclaimed, cropland and forest compete for water, leading to more dying ecological forests and aggravated desertification.

Methodology and Data
This paper developed a computable general equilibrium of land use change (CGELUC) model to simulate the changes of land use structure in Zhangye city of Heihe River Basin, China, in 2015-2030 under three development scenarios: baseline, resource consumption, and green development. The dynamic land system (DLS) model was further applied to spatialize land use change on a spatial grid scale. The research framework is shown as follows ( Figure 3):

CGELUC Model
In this study, we developed a computable general equilibrium of land use change (CGELUC) model to simulate the changes in land use structure by using the ORANI-G model, which is a single-country CGE model and is widely applied to the analysis of economic policies and academic research around the world (https://www.copsmodels. com/oranig.htm) (Accessed on 1 August 2018) [35]. In this study, we expanded the ORANI-G model by embedding a land resources account. Similar to the general CGE model, this model includes production, consumption, trade, tax, income distribution, markets and other modules. The dynamic recursive mechanism is established in the model, which mainly aid in scenario prediction and the analysis of future land system structure dynamics. The model uses Gempack software (Centre of Policy Studies, Victoria University, Melbourne, Australia) as the model running environment.

Production Module
The following figure describes the structure of the production module of the CGELUC model. The module encompasses multi-level nested production functions (Figure 4), each describing an optimal combination to ensure the minimization of production costs. The top-level is the optimal combinations of goods, production factors, and other costs using Leontief functions. The second level is the combination of domestic and imported commodities using the CES function. Land resources are embedded in the production module as a production factor and nested with other production factors, labor and capital, by the CES production function.
Considering the heterogeneity of land resources and the constraints of land conversion, total land input is divided into different uses according to conversion elasticity and relative rent by applying the CET function. In this study, we adopted the method proposed by Horridge [54] to amend the CET function in order to achieve the allocation of areas of different land use types under the total fixed land area, which avoids the errors in the traditional CGE model when it is used to explain changes in physical units.
Description of variables in equations are listed in Appendix A.

Construction of Dynamic Mechanism
The static CGE model describes the behavior and characteristics of the main economic system in a given period. In this study, the model introduced a dynamic module to simulate the impact of socioeconomic factors on land use change under different future scenarios. The model introduces recursive dynamic structures of capital accumulation and employment to depict the evolution of the economic system, in order to track the trend of economic change over time under continuous policy and economic shocks and predict dynamic long-run macroeconomic change. Long-term forecasting is carried out on the basis of the dynamic accumulation of capital, which affects the total capital of the next round of production equal to the capital stock of the previous period minus depreciation plus total investment. In addition, the model takes future changes of labor employment into account, constructing a dynamic adjustment mechanism for wages and employment.
Other modules and specific explanations can be found in Bai et al. [41].

The Structure of DLS
The dynamic of land system (DLS) model is based on the theory of equilibrium in regional land use structure change and the constraint of land use type distribution at the grid scale. Based on a comprehensive consideration of the driving mechanism of natural factors and socioeconomic factors in land use changes, a dynamic simulation of spatial pattern evolution in regional land use change is realized.
The DLS model consists of the following modules: (1) Driving analysis module of land use distribution-this module is used to measure the spatial relationship between natural environmental factors, climatic factors, location factors, socioeconomic factors and land use distribution, and depict the control and transformation effect of these driving factors on land use; (2) Spatial allocation module-based on the probability of land use conversion, land use conversion rules and land use conversion area restrictions, the spatial allocation of land use change is effectively simulated through the land supply and demand balance at the grid scale.

Construction of DLS Model
The DLS model fully considers the driving characteristics of regional land use type distribution and the characteristics of land supply and demand balance, reflecting the systems, mechanisms, equilibrium and convergence of regional land systems.
The changes in land supply and demand in the DLS model are a dynamic convergent process in which the total amount remains in equilibrium. The raster-scale land supply and demand balance requires a highly accurate raster-scale land supply and demand balance model, which is an empirical statistical model with spatially explicit typical characteristics. By considering the effects of local drivers within each raster and the neighborhood effects of drivers around the raster, the competition process between land use types is finely depicted. The core is the process of the spatial allocation of land demand at the micro-scale ( Figure 5). It is assumed that the probability of land use type k occurring, based on the local raster factor i, is p k i,Local , based on the neighborhood raster factor i being p k i,NBH . The combination probability of land use, p k i = f p k i,Local , p k i,NBH , reflects the self-organization process of land use system evolution. p k i,Local and p k i,NBH are obtained from the nonlinear constraint model of land use type distribution at the raster scale: where p k i is the probability of the raster i displaying the land use type of k, and X il is the driving factor l in the raster i, which can be a local factor or neighborhood factor. a k i is the regression coefficient and r is the influence coefficient of the spatial autocorrelation factor. The spatial autocorrelation factorŷ k i is also defined.

Natural Data
Land use data and other natural environment data, including landform type, DEM, slope, aspect, etc., are provided by the Resource and Environment Science and Data Center, Chinese Academy of Science (http://www.resdc.cn/dataResource/dataResource.asp) (accessed on 1 August 2018). Climate data, including temperature, precipitation, sunshine duration (sun), ≥0 • C accumulated temperature, and ≥10 • C accumulated temperature, are provided by the Chinese Meteorological Data Center (http://data.cma.cn/) (accessed on 1 August 2018). Soil data including soil organic, soil depth and soil PH value are provided by the Harmonized World Soil Database (HWSD). Social economic data such as population and GDP at 1 km*1 km resolution and date related to distance are also derived from the Resource and Environment Science and Data Center, Chinese Academy of Science (http://www.resdc.cn/dataResource/dataResource.asp) (accessed on 1 August 2018).

Value Accounting of Land Resources
The input-output table of 48 sectors in Zhangye city in 2012 is a value input-output table. To embed the land resources account into it, it is necessary to carry out the value accounting of land resources. Referring to the System of Integrated Environmental and Economic Accounting (SEEA), the value of agricultural land is multiplied by land rent, and the value of industrial and service land is multiplied by land transfer price divided by service life.
Based on the sectoral classification of the input and output tables of 48 sectors in Zhangye city in 2012, the land areas of seven agricultural sectors, including corn, wheat, oil crops, vegetables, cotton, fruit and other agriculture, were obtained, and the land price was calculated according to the average contracted price of local agricultural land, which is CNY 0.75/m 2 . Through remote sensing interpretation technology combined with field investigation, the areas of land functional types of 41 secondary and tertiary industrial sectors were determined, and the land prices were calculated according to the average price of the local land auction divided by the service life, which are respectively CNY 56.85/m 2 and CNY 13.66/m 2 .

Scenario Design
There are many factors influencing the change of land use demand in social economic system. In national "Silk Road" economic belt planning in China, as the core zone of the Hexi Corridor, the rapid development of urbanization in Zhangye city is inevitable. Regional population, economic aggregate, industrial structure, consumption pattern, technology and management level all affect the land use demand of social and economic system. With the growth of the economy and population, the improvement of urbanization level and the transformation of consumption pattern, the demand for land use in social and economic systems would continue to rise, while the optimization and transformation of economic structure, the improvement of scientific and technological level and the progress of management measures would reduce the land use demand and improve the efficiency of land resource use. Therefore, it is of great practical significance to clarify the change trend of driving factors of socioeconomic development scenarios and simulate the land use demand of socioeconomic systems for guiding sustainable land management in Zhangye. Based on the literature [56,57] and regional development planning, the projection of population, investment, consumption, economic growth and trade in Zhangye city in 2030 was carried out. Based on the CGELUC model, long-term closure with employment growth as the main factor was set. Meanwhile, the trend of future social and economic development is simulated according to the prediction index of future social and economic development (Table 3), and the baseline scenario (BAU) of future social and economic development is constructed. This study simulates land use dynamics in Zhangye city of Heihe River Basin, China in 2015-2030 under three development scenarios: baseline (BAU), resource consumption (RCS), and green development (GDS). Referring to the BAU scenario, the trends of future social and economic development are also simulated in RCS and GDS scenarios ( Table 3). The baseline scenario simulates the trend of future socioeconomic development and the demand for land use in the socioeconomic system with the same changes in investment, exports, resident consumption, capital input, population employment and labor technological progress from the past to the present. The resource consumption scenario gives priority to rapid economic growth and social development, and has a stronger impact on the demand for land use in the socioeconomic system by being shocked by higher investment, exports, residents' consumption, capital, population employment and labor technological progress. The green development scenario makes more consideration of regional ecological environment construction, assuming that ecological security and ecological restoration will be mainly considered in future urban development. Compared with the basic scenario, investment, export, resident consumption and capital investment have a weaker impact on the demand for land use in the socioeconomic system, while the assumptions of population employment and labor technological progress remain unchanged.

Land Use Change in Social-Economic System under Different Scenarios
According to the setting of the change rate of the shock variables, it is predicted that by 2030, the GDP under the baseline scenario will reach CNY 77 billion, with an annual growth rate of about 6%. This is consistent with regional development planning, showing that the scenario's design is relatively reasonable.
This paper simulates the land demand of the socioeconomic system in Zhangye city from 2015 to 2030 under the three future scenarios. Figure 6 shows the percentage of land use change in different industries predicted up to 2030. The changes under the resource consumption scenario are the most intense, followed by the baseline scenario, and finally the green development scenario. From the perspective of industrial structure, the land use demands for wheat, corn and cotton in agricultural sectors decrease by more than 10%. The reason for this result is, on the one hand, with the development of modern agriculture, the efficiency of land resource utilization in traditional crop production is constantly improving, intensive management is strengthening, and the quality of cultivated land is improving. On the other hand, the internal structure of agricultural sectors in Zhangye city is optimized and adjusted. Agriculture has gradually changed from planations to the dual structure of grain and industrial crops, and even the triadic structure of grain, industrial crop and grass, which promotes a virtuous circle in the internal structure of agricultural sectors. In the secondary industry, the steady development of the social economy has stimulated the demand for land in mining, smelting and processing, food manufacturing, electric heat supply, construction and other industries. In particular, the demand for land in the general and special equipment manufacturing industry has increased by more than 50%. The land demand for most of the tertiary industries also showed an increasing trend, especially in the land demand of the real estate industry, which increases by more than 130% year-on-year. The land demands of other modern service industries, such as environment and public facilities management, information software technology, finance and health, social security and social welfare, are also stimulated by the socioeconomic system.  Figure 7 shows the percentage change of land demand for the three major industries under different scenarios in the future. From the perspective of the industrial structure, the land demand of the primary industry shows a downward trend on the whole, the land demand of the secondary and tertiary industries shows an increasing trend, and the land demand of the tertiary industry is higher than that of the secondary industry. The analysis shows that this is the result of the continuous optimization and upgrading of the industrial structure and the outward transfer of rural labor. The simulation results show that the land demand of the three major industries shows a development trend in the same direction under the three scenarios, but the change degree is different. The analysis shows that the continuous growth of investment, consumption and export will promote the rapid development of the economy; the increase in population employment and the technological progress of labor force will promote industrial transformation and upgrading; the economic structure will be adjusted and optimized, so that the land demand of the agricultural sector is reduced; the land demand of the secondary and tertiary industries will increase, and within a certain level of development, the greater impacts of the social and economic systems are affected by the coupling of investment, consumption, export, capital quantity, population employment and technological progress of labor force. The lower the land demand of the agricultural sector, the greater the land demand of the secondary and tertiary industries.

Comparative Analysis of Land Use Structure Changes under Different Scenarios
Based on the analysis of the changes in land demand in different industrial sectors under different future scenarios, 7 agricultural sectors and 41 industrial and service sectors are merged (Figure 8). It is found that the same land use types present the same direction variation trend under the three scenarios, but the change degree is different.

Spatial Patterns of Land Use Change under Different Scenarios
The DLS model is further applied to simulate the land use spatial pattern of Zhangye city. Based on the land use data interpreted by remote sensing in 2015, the simulated results obtained are compared with the actual land use spatial pattern to test the accuracy of the DLS model. The consistency of the simulation results can be quantitatively reflected by using the Kappa index. In this study, the Kappa index of the simulation results was calculated as 0.9036, indicating the robustness and reliability of the DLS model.
The land use spatial pattern of Zhangye city is shown in Figure 9. It is found that the built-up area expands within a certain scale under all three scenarios, but the expansion is most pronounced under the resources consumption scenario, with large-scale concentrated and contiguous distribution occurring in northeastern Ganzhou District and southwestern Gaotai County (Figure 9b); the second is under the baseline scenario, wherein the builtup area extends mainly in the northeastern industrial park of Ganzhou District and its surrounding Linze County, with a small concentration in the southwestern Gaotai County (Figure 9a). In contrast, built-up area expansion is slowest in the green development scenario, with a relatively even and discrete spatial layout (Figure 9c  Simulations of the spatial pattern of land use under different future scenarios provide reference information for the optimal allocation and management of regional land resources, sustainable economic development and ecological protection policy decisions. The simulation results show that cultivated land has a significant trend of contraction in the central part of Shandan and Minle Counties. However, these two counties mainly undertake the function of agricultural development in Zhangye. The government should pay attention to these sensitive areas of land use change, and take the relevant measures to strengthen the protection of cultivated land and basic farmland in these areas while strictly controlling the unreasonable conversion of land use. At the same time, to ensure ecological safety, the government should formulate a reasonable land resource development plan and industrial restructuring plan, considering the socioeconomic and ecological benefits of urbanization development. In the meantime, the transformation and leapfrogging of Zhangye will be promoted to achieve the goal of building "Ecological Zhangye, Green Zhangye", so that the land use pattern of Zhangye can develop in a benign and sustainable direction that is conducive to the harmony between humans and land.

Industrial Structure Adjustment and Sustainable Land Management
Combined with the current development in Zhangye city and some relevant policies, the above simulation results are reasonable. Against the background of China's new round of Western development, the State Council, the Ministry of Natural Resources and relevant departments have issued a series of policies and measures to vigorously support Gansu Province in speeding up urbanization construction. Zhangye is in a new stage of rapidly promoting industrialization and urbanization development. With the rapid development of urbanization, the urban population increases, much of the labor force is transferred to the city, and the transfer of the industrial structure from the primary industry to the secondary and tertiary industries is becoming increasingly obvious, which leads to a continuous increase in construction land area and a continuous decrease in cultivated land area. The continuous growth of investment, consumption and exports will promote rapid economic development, and the increases in population employment and the technological progress of the labor force will promote industrial transformation, upgrading and optimization. Therefore, among the three scenarios in this study, exogenous impact factors such as investment, consumption, export, labor force and technological progress have the most significant impact on land use change in resource consumption scenarios, followed by the baseline scenario, and have the least impact on the green development scenario. According to the research on land use change in the socioeconomic system under different scenarios in the future, this paper puts forward relevant suggestions regarding urbanization development and land use planning in Zhangye city. The government should vigorously promote industrial upgrading, on the one hand, to realize the optimization of industrial structure in agriculture, and on the other hand, to actively develop the modern service industry and the tourism industry, encourage industrial sectors to strengthen technological innovation, and drive the rapid development of the regional economy. At the same time, the simulation results show that the cultivated land is reduced sharply under the resources consumption scenario. In order to ensure food security, the local government should take relevant measures to strengthen the protection of cultivated land and basic farmland, reasonably and responsibly develop and utilize land resources, and effectively supplement the cultivated land occupied by construction land.

Urban Development Model Based on Land Use Change
In the accelerating process of urbanization, the construction of Zhangye will manifest a great change in the pattern of land use. Zhangye is still an underdeveloped region. Economic development is still the main goal of urbanization construction. Therefore, it is easy for Zhangye to be directed by the current policy towards the one-sided pursuit of the economic benefits of the development model. The ecosystem structure of Zhangye is fragile. The natural background conditions are not superior. As an oasis city with high development intensity in the Heihe River Basin, the contradiction in Zhangye city between people and land is becoming more and more prominent in the process of rapid development. Against the strategic background of "sustainable development" and the "ecocity" concept, the rational co-ordination of ecological and economic benefits has become an unavoidable problem in the future development model of Zhangye. This study confirms that under the BAU, RCS and GDS scenarios, cultivated land is reduced to a certain extent due to the impact of the large-scale occupation of built-up areas. Since the area suitable for construction coincides with the area of high-quality cultivated land, it is easy to cause great tension between urban infrastructure construction and farmland protection. To ensure food security, relevant departments should take relevant measures to strengthen the protection of cultivated land and basic farmland, stick to the red line of cultivated land, and ensure that farmland occupied by construction will be effectively supplemented. At the same time, the ecological policy of returning farmland will also lead to the contraction of cultivated land. In the GDS scenario, the conversion of ecological land mainly occurs in Sunan County, the northern part of Gaotai County, and the Heihe Wetland in Ganzhou District and around the Dongdashan Nature Reserve. Therefore, it can also be combined with the simulation results of land use change under the GDS to explore a reasonable model for converting a large quantity of unused land into ecological land, and reduce the scale of cultivated land contraction, which can activate both the socioeconomic and ecological benefits of urbanization development.

Future Prospect
(1) The value accounting method of embedding land resources into the value inputoutput table has a great influence on the accuracy of the model simulation results. Due to the temporal and spatial differences of regional land rent and land transfer price, the value of the three major industrial land types in Zhangye city is relatively simplified by using the unified average price, so it is necessary to carry out more accurate accounting of land resources. (2) This study addresses three scenarios constrained by socioeconomic development to simulate future land use change in Zhangye city. Combined with the actual situation of the Heihe River Basin, the next step for promoting the research is to carry out scenario design under different environmental development policies, such as water resource constraint scenarios, so as to provide a more comprehensive reference for the optimal allocation of land use for sustainable development. (3) On the basis of static CGE, the model introduces capital accumulation and future trend of employment to realize a dynamic recursive model for long-term scale simulation and prediction. However, due to the partial delay of capital and investment transformation, the prediction results are unreasonable in the first few years. The model needs to be improved. Besides, we will carry out further studies on land use change and natural resource management on the watershed scale in the Heihe River Basin. Therefore, the CGELUC model will also be improved in cross-regional studies.

Conclusions
Taking Zhangye city of the Heihe River Basin, China, as an example, this paper constructs three scenarios-the baseline scenario, the resource consumption scenario and the green development scenario-to simulate land use dynamics from 2015 to 2030 by developing the CGELUC model. The DLS model is selected to dynamically simulate the spatial patterns of land use under three different future scenarios. The main conclusions are as follows: (1) Analyzing land use demand in the socioeconomic system under future scenarios, it is found that the land use demand in the three scenarios develops in the same development direction, but the degree of change is different-from the least to the greatest, they rank as the green development scenario, the baseline scenario and the resource consumption scenario. From the perspective of industrial structure, the land demand for primary industry is decreasing, the land demand for the secondary and tertiary industries is increasing, and the land demand for the tertiary industry is higher than that of the secondary industry. From the perspective of the internal industrial structure, the land demand for wheat and corn in the agricultural sector is showing a decreasing trend, but the land demand for vegetables and other agricultural sectors is increasing. For the secondary and tertiary industries, except for a few traditional industrial sectors, the land use demand for most sectors is increasing; (2) Analyzing the structure of land use change under future scenarios, it is found that land use change in the three scenarios develops in the same development direction, but the degree of change is different-from the least to the greatest, these are ranked as the green development scenario, the baseline scenario and the resource consumption scenario. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A   Table A1. Description of variables in equations.  Other agriculture 8

PPRI M i
Mining and washing of coal 9 Extraction of petroleum and natural gas 10 Mining of metal ores 11 Mining and processing of nonmetal ores and other ores 12 Manufacture of foods and tobacco 13 Manufacture of textiles 14 Manufacture of textile wearing apparel, footwear, caps, leather, fur, feather (Down) and its products 15 Processing of timbers and manufacturing of furniture 16 Papermaking, printing and manufacturing of articles for culture, education and sports activities 17 Processing of petroleum, coking, processing of nuclear fuel 18 Chemical industry 19 Manufacturing of nonmetallic mineral products 20 Smelting and rolling of metals 21 Manufacture of metal products 22 Manufacture of general-purpose and special-purpose machinery 23 Manufacture of transport equipment 24 Manufacture of electrical machinery and equipment 25 Manufacture of communication equipment, computer and other electronic equipment 26 Manufacture of measuring instrument and machinery for cultural activities and office work 27 Manufacture of artwork, other manufacture 28 Scrap and waste 29 Production and supply of electric power and heat power 30 Production and distribution of gas 31 Production and distribution of water 32 Construction 33 Traffic, transport and storage 34 Post 35 Information transmission, computer services and software 36 Wholesale and retail trades 37 Hotels and catering services 38 Financial intermediation 39 Real estate 40 Leasing and business services 41 Research and experimental development 42 Comprehensive technical services 43 Management of water conservancy, environment and public facilities Health, social security and social welfare 47 Culture, sports and entertainment 48 Public management and social organization