Next Article in Journal
The Outward Extension of an Ecological Footprint in City Expansion: The Case of Beijing
Next Article in Special Issue
Using Urban-Carrying Capacity as a Benchmark for Sustainable Urban Development: An Empirical Study of Beijing
Previous Article in Journal
Towards Life Cycle Sustainability Assessment of Alternative Passenger Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Projection of the Relationship between Industrial Structure and Land Use Structure in China

1
School of Mathematics and Physics, China University of Geosciences (Wuhan), Wuhan 430074, China
2
Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
3
Center for Chinese Agricultural Policy, Chinese Academy of Sciences, Beijing 100101, China
4
State Key Laboratory of Water Environment Simulation, School of Environment, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Sustainability 2014, 6(12), 9343-9370; https://0-doi-org.brum.beds.ac.uk/10.3390/su6129343
Submission received: 3 September 2014 / Revised: 14 November 2014 / Accepted: 17 November 2014 / Published: 16 December 2014

Abstract

:
Based on the computable general equilibrium (CGE) modelling method, this research analyzes the relationship between industrial structure and land use structure in China. The results show that our model is feasible, and the simulation results are of a certain stability. Under the scenario analysis and projection of the relationship between the industrial structure and land use structure of the thirty-one provinces in China from 2010 to 2020, the proportions of secondary and tertiary industry in each province have been increasing; correspondingly, the proportion of agriculture has been decreasing. This means that the industrial structure of China is changing. As for land use, in general, the trend is similar to the industrial structure changes. The transformation of the structure of industrial development and land use has driven economic structure changes in China. The economic structure has an inclination to transform from agriculture to both secondary and tertiary industry. Along with industrial transformation, the cultivated land in China shows a trend of continuous decline. Empirical analysis results indicate that a decrease of cultivated land is acceptable under the scenario of economic growth in the next ten years. This shows a possibility that the economic efficiency of land use for cultivation and business services will decline, and more attention ought to be paid to increasing the economic efficiency of land use.

1. Introduction

Land is the foundation of human survival, as well as the basic material for social production [1]. Land use is a process in which humans transform natural land to meet the demands of socio-economic development. Various types of land use are formed in that structural transformation process. It directly affects the human activities of using and planning their living environment [2]. Land use structure is the proportional land use that reflects the pattern of socioeconomic development in an economic structure [3]. The industrial structure represents the interrelationships and the interdependence of sectoral production. Their proportional relationships among sectors systematically demonstrate the integrity of the interaction of the social production process and its proportional changes in the reformation of social production [4]. In other words, it shows the composition of each sector in the economic structure [5]. Grasping the trend of land use change and its impacts on the improvement of sectoral production efficiency for economic development is of unusual significance.
Industrial structure and land use structure are closely related. Industrial structure is the structure of the proportions of various industries in a specific region. Land use structure is the structure of the shares of different types of land use for various industries in a specific region. Because land is the fundamental factor for production activities, the adjustment of the economic industrial structure will inevitably change the demand for land use and, thus, alter the land use structure [6]. Kuznets (1955) said that the proportion for the industrial and service sectors would increase until the economy developed to a certain level. The production factors of agriculture would be transferred to secondary and tertiary industry [7]. The investment-oriented and policy-oriented adjustment of the industrial structure can drive the reallocation of production inputs. In return, land use structure interacts with the structural changes of industrial production. Meanwhile, economic structure is the structure based on the industrial structure under the designed scenarios. Government strategic policy on sustainable development somehow reshapes this interaction process for strategically guaranteeing the balance between economic growth and resource utilization, which guides the adjustments of industrial structure and land use structure.
The computable general equilibrium (CGE) model can quantitatively describe the relationships and interdependencies among components of various economies within a complex social system [8,9,10]. Theoretically, this model reflects the structure of a bundle of systematic equilibriums through a market mechanism to depict both the supply and demand of commodities and services in a unit pricing system. Moreover, by giving a well-structured depiction of physical land use, the application of the CGE model has been extended to analyze land use analysis based on traditional market economy theory.
In this paper, based on an analysis of the current pattern of land use in China, we examined the relationship between industrial structure and land use structure in the thirty-one provinces (Note: Taiwan, Hong Kong and Macao have not been included, due to a lack of statistical data). This will help us research future decision making in the optimization of land use structure and contribute to further reformation of industrial development under advanced land management planning. This model and data description will be presented in the following Section 2. Thereafter, empirical analysis and simulation results will be given in Section 3. At the end of this paper, we provide the conclusions of this research and a discussion of further research directions.

2. Model and Data

2.1. Multi-Regional CGE Model

Based on multi-regional CGE modelling, the constraints of land use were introduced into the potential productivity process. Through a mathematic operation research approach by using time series data of land use and its demand changes, a robust land use prediction model was built. In the following sections, three nested production functions are combined with the coordinated environmental sustainability (CES) production function and Leontief production function, which will be introduced into the CGE model, as in Figure 1.
Figure 1. Nested production function structure. CES, coordinated environmental sustainability.
Figure 1. Nested production function structure. CES, coordinated environmental sustainability.
Sustainability 06 09343 g001
The third level of production functions includes capital and labor, which represents an integrated value. The second level includes two production levels. One is the integration of labor and capital and land, which are converted to added value through the transaction of the CES production function. The other one is the combination of intermediate inputs through the transaction of the Leontief production function. The function at each level of production functions as follows:
The first level of nested CES final demand:
Q i = A i ( η A V i ρ i + ( 1 η ) M i ρ i ) 1 ρ i
The second level of nested CES added value:
A V i = A i A V ( δ i A V I V i ρ i A V + ( 1 δ i A V ) L N D i ρ i A V ) 1 ρ i I V
The second level of nested CES intermediate inputs:
M i = min { a 1 M 1 , a 2 M 2 , L , a n M n } = a j M j
The third level of nested CES intermediate inputs of labor and capital:
I V i = A i I V ( δ i I V L i ρ i I V + ( 1 δ i I V ) K i ρ i I V ) 1 ρ i I V
The production process is a process to pursue maximum profit (or the minimum of the production cost). When profits reach the maximum level (or the cost reaches the minimum level), each factor input is the factor demand. Except the above input-output production function, solving the objective function also requires introducing the price relationships. Under this condition of the nest CES structure, if the objective function is the cost minimum function, an objective function of a certain good will be represented as follows:
M i n .   C i = P i A V A V i + P i M M i
The relevant price relationships are:
p i Q Q i = P i A V A V i + p i M M i
p i A V A V i = P i I V I V i + p i L N D L N D i
p i I V I V = p i L L i + p i K K i
p i M M i = i j p X i M i
where and below,
  • i is the sector of crop, wool, other agriculture, secondary industry, construction and service industry;
  • j is the other sectors of i;
  • t, is the year from 2000–2020;
  • r is a region of the31 provinces in China;
  • l is the other regions of r;
  • d is the land use type of crop, wool, other agriculture, secondary industry, construction and service industry.

2.1.1. Constraints of Land Use

According to the restrictions of land use and land planning, the constraints of land use are shown in the following equations.
The sum of the total area of land use for each sector in each province is smaller or equal to the total area of land in each province:
i = 1 6 X L N D i , r X L N D _ I I r
where i is the sector and r is the region.
Land supply is not smaller than the inventory in the planning period:
X L N D 1 , r X L N D P c r o p , r
where X L N D P c r o p , r is the inventory in the planning period.
Land supply in each region is not beyond the upper bound of the inventory in the planning period:
X L N D 4 , r + X L N D 5 , r + X L N D 6 , r X L N D P c o n s , r
where X L N D P c o n s , r is the upper bound of the inventory in the planning period.
In addition, in the simulation process, with respect to other land use types, built-up land is special. In the short term, the built-up area will not shrink; in other words, the current built-up area will not be smaller than it was in the last year, no matter how much the outputs increase in the area.
X c o n s L N D t X c o n s L N D t 1

2.1.2. Agricultural Production Impact Function

Apparently, each sector can be influenced by a combination of factor inputs and land use changes. Relatively speaking, agricultural production can be influenced by land use changes, which are more than the proportional changes in sectoral production. Usually, agricultural production is directly and fundamentally affected by land use changes, and this kind of influence will easily transfer to the secondary and tertiary sectors, sequentially affecting the whole economic system. Therefore, we introduce the agricultural production impact function for comprehensively assessing the impact factors of the agriculture sector. We assume this might influence the secondary and tertiary sectors by using the dynamic multi-regional CGE model.
Agricultural production usually is influenced by many impact factors, such as natural hazards, human factor intervention, production planning reformation and technology innovation. Among those factors, natural disasters have the most impact on agricultural production. In 2013, agricultural production had a loss of 125.6 billion CNY due to natural disasters. Therefore, we focus on how much impact natural hazards have on land use, because this mainly induces agricultural production loss, in order to analyze the possible impacts on the secondary and tertiary industries in the next ten years.
The impact factor of natural disasters is represented by the ratio of the production loss area over the total area of land use in the following equation:
λ i , t = D Q i , t L N D i , t
where λ i , t is the coefficient of the i type of land suffering from natural hazards, D Q i , t is the area of the i type of land with production loss in the t year and L N D i , t is the total area of the i type of land in the t year.
It is easy to figure out that natural hazards mainly influence cultivated land productivity and, then, influence the cultivated land supply. Moreover, natural hazards have less impact on other types of land use, such as built-up land, unused land, and so forth, and could be ignored. Thereby, when cultivated land converts to other types of land use, we consider this as D Q i , t = 0 . Thus, the natural hazard impact factor on each type of land use is as follows:
λ i , t = { D Q i , t L N D i , t , i = 1 , 2 , 3 0 , i = 4 , 5 , 6
Natural hazards mainly have impacts on productive cultivated land, woodland and grassland. Thereafter, there are three different levels to describe the impact of natural hazards on agricultural lands, which represent disaster, harmed and lost. Disaster represents less than 10 percent of the expected agricultural productivity being lost. Harmed land represents about 30 percent, but less than 70 percent of the expected agricultural productivity being lost. Lost land represents more than 70 percent of the expected agricultural productivity being lost. Therefore, the total amount of agricultural productivity loss due to natural hazards is computable in the following equations:
Total area of agricultural productivity loss = flood impacted land + hailstorm impacted land + frost impacted land + drought impacted land + hurricane impacted land = (flood impacted land *0.1 + flood impacted land * 0.3 + flood impacted land * 0.7) + (hailstorm impacted land * 0.1 + hailstorm impacted land * 0.3 + hailstorm impacted land * 0.7) + (frost impacted land * 0.1 + frost impacted land * 0.3 + frost impacted land * 0.7) + (drought impacted land * 0.1 + drought impacted land * 0.3 + drought impacted land * 0.7) + (hurricane impacted land * 0.1 + hurricane impacted land * 0.3 + hurricane impacted land * 0.7).
The above conceptual equation is represented by the mathematical equation as follows:
D Q i , t = f l Q i , t + w i Q i , t + l t Q i , t + d r Q i , t + t y Q i , t
where,
f l Q i , t = f l S Q i , t 0.1 + f l C Q i , t 0.3 + f l J Q i , t 0.7
w i Q i , t = w i S Q i , t 0.1 + w i C Q i , t 0.3 + w i J Q i , t 0.7
l t Q i , t = l t S Q i , t 0.1 + l t C Q i , t 0.3 + l t J Q i , t 0.7
d r Q i , t = d r S Q i , t 0.1 + d r C Q i , t 0.3 + d r J Q i , t 0.7
t y Q i , t = t y S Q i , t 0.1 + t y C Q i , t 0.3 + t y J Q i , t 0.7
A natural disaster happens suddenly, which may bring about huge, unpredictable loss. Based on the 2000–2011 historical data, which is produced by Jin Qing [11], we use the average of that value to predict the impacts of a natural disaster on land use. According to Equation (16), we obtain the data for the total area of land suffering from a natural disaster in each province of China from 2000 to 2011. Sequentially, according to Equation (15), the impact factor of natural disasters on land use and changes in each province of China from 2000 to 2011 can be calculated. Finally, the average of the impact factor of natural disasters on land use and changes in each province of China from 2000 to 2011 are listed in Table 1, which can be used for forecasting natural disaster impacts on land use and changes in each province of China.
Table 1. Impact factors of natural disasters on land use in each province of China * (104 hm2).
Table 1. Impact factors of natural disasters on land use in each province of China * (104 hm2).
RegionDisaster AreaAgricultural AreaNatural Disaster Impact FactorRegionDisaster AreaAgricultural AreaNatural Disaster Impact Factor
China1,749.8915,463.580.1132Henan92.191,335.980.0690
Beijing6.3034.200.1842Hubei100.55735.500.1367
Tianjing8.6052.280.1645Hunan133.04778.990.1708
Hebei125.33893.510.1403Guangdong59.60493.240.1208
Shanxi58.22390.050.1493Guangxi46.76629.990.0742
Inner Mongolia96.55588.700.164Hainan11.1086.330.1286
Liaoning62.06380.920.1629Chongqing49.31346.710.1422
Jilin52.14468.770.1112Sichuan71.56956.460.0748
Heilongjiang128.92985.840.1308Guizhou44.31464.540.0954
Shanghai1.2947.670.027Yunnan51.82581.310.0891
Jiangsu24.30779.790.0312Tibet1.3123.290.0562
Zhejiang23.68306.450.0773Shaanxi74.20419.830.1767
Anhui43.45899.760.0483Gansu47.47364.990.1301
Fujian30.85266.140.1159Qinghai12.6549.430.2559
Jiangxi43.02535.510.0803Ningxia16.44114.790.1432
Shandong209.131,104.780.1893Xinjiang23.73347.830.0682
* Note: natural disasters include: flood, hailstorm, frost, drought and hurricanes; data from China Statistic Yearbook 2001–2012 and China Agricultural Statistic Yearbook 2001–2012.

2.1.3. Land Demand Prediction Model

Here, the growth rate of land is equal to the growth rate of production; then, the growth of land demand can be represented as follows:
A _ L N D i , d = A _ Q i , d A _ T F P i , d + 1
where A _ L N D i , d is the growth rate of land demand, A _ Q i , d is the growth rate of production and A _ T F P i , d is the growth rate of other factor inputs.
Under the constraints of the total number of each type of land, the minimum deviance of the land demand and its satisfaction coefficient in a certain period are set up in the objective function.
min i ( R D L i , d , Δ t M D L i , d , Δ t ξ i , d , Δ t 1 ) 2
s . t . i R D L i , d , Δ t = 0 ,   i = 1 , 2 , , 6 ; d = 1 , 2 , , 31
where R D L i , d , Δ t is actual changes of land at time t, M D L i , d , Δ t is land demand at time t and ξ i , d , Δ t 1 is the satisfaction coefficient in the previous time t − 1, which is equal to the ratio of land use changes divided by the land demand in previous time t − 1, where the previous land use changes are calculated by the LaGrange multiplier method within the prediction model as follows:
(1)
Auxiliary function:
Φ = i = 1 6 ( R D L i , d , Δ t M D L i , d , Δ t ξ i , d , Δ t 1 ) 2 + λ i = 1 6 R D L i , d , Δ t
(2)
Take the derivative of R D L i , d , Δ t , and solve for the static point of Φ :
Φ R D L i , d , Δ t = 2 M D L i , d , Δ t ( R D L i , d , Δ t M D L i , d , Δ t ξ i , d , Δ t 1 ) + λ = 0
(3)
According to Equation (21), the relationship of R D L i , d , Δ t and λ can be derived:
R D L i , d , Δ t = M D L i , d , Δ t 2 2 ( ξ i , d , Δ t M D L i , d , Δ t λ )
(4)
Plug Equation (22) into the equation i = 1 6 R D L i , d , Δ t = 0 to get λ as the following equation:
λ = i = 1 6 ξ i , d , Δ t 1 M D L i , d , Δ t i = 1 6 M D L 2 i , d , Δ t
(5)
According to Equations (22) and (23), the equation is reached:
R D L i , d , Δ t = M D L i , d , Δ t 2 2 ( ξ i , d , Δ t M D L i , d , Δ t + i = 1 6 ξ i , d , Δ t 1 M D L i , d , Δ t i = 1 6 M D L 2 i , d , Δ t )

2.2. CES Production Function Elastic Parameter Estimation

Bayesian estimation of Markov chain Monte Carlo methods (MCMC) [12] is used for estimating the elasticity of the CES production function, to renew the original data of the parameters and to implement dynamic parameter inputs. The parameter estimation is in Table 2 and Table 3. All of the numbers for the constant elasticity of substitution are generated by the Bayesian estimation method, which are all inside of the range of the CES production function from 0.1 to 2.5 [13]. Therefore, the elasticities used in this research are acceptable as parameter inputs to the multi-regional CGE model (see the Appendix).
Table 2. The elasticity of the substitution between the land and the integrated value of labor and capital in each sector.
Table 2. The elasticity of the substitution between the land and the integrated value of labor and capital in each sector.
SectorElasticity of SubstitutionSectorElasticity of Substitution
Planting1.99Manufacturing0.86
Dairy1.99Construction0.26
Other Agriculture1.99Tertiary0.26
Table 3. The elasticity of substitution between added value and intermediate inputs in each sector.
Table 3. The elasticity of substitution between added value and intermediate inputs in each sector.
SectorElasticity of SubstitutionSectorElasticity of Substitution
Planting0.34Manufacturing2.01
Dairy0.34Construction2.02
Other Agriculture0.34Tertiary1.01

2.3. Sensitivity Analysis

Assessing the CGE model usually requires examining whether the tiny changes of some generally key parameters have impacts on the optimum solution. Therefore, sensitivity analysis is used to test the robustness of the optimum solution.
Assume the key parameter as θ ; the optimum solution of the function is F ( X , θ ) = 0 is X * . With respect to one of the neighborhoods Ο ( θ Δ θ , θ + Δ θ ) of θ , if θ Ο ( θ Δ θ , θ + Δ θ ) , which let F ( X , θ ) = 0 have the optimum solution at X * ˈ , then ε > 0 and | X * ˈ X * | < ε . Here, we call the model robust. Theoretically, a smaller ε is better. However, actually, as long as | X * ˈ X * | stays in a certain acceptable range, we can also claim that the model is robust.
It is necessary to perform numerous extractions in the neighborhoods to statistically describe the average and variance of the model or the given confidence interval [14], so that we reach a reliable conclusion.
The elasticity of substitution of the CES production function is distributed in (0.1, 2.5). Thereby, when doing sensitivity analysis on that, we extract the elasticity of substitution from the average distribution in (0.1, 2.5). Therefore, in the simulation process, the elasticity of substitution is randomly selected from (0.1, 2.5). Computing program of the model run 1000-times under this scenario of two percent of the growth in all factor inputs, and then, the elasticity of substitution of the endogenous variables in the CGE model within the 95% confidence interval is derived in Table 4 (See Appendix).
Table 4. Sensitivity analysis results *.
Table 4. Sensitivity analysis results *.
VariableAverageSDLower Bound of CIUpper Bound of CI
ES11.9190.0321.8922.033
ES20.8420.0130.6780.902
ES30.2580.2010.2960.287
ES40.2620.4500.2640.313
ES50.3410.0230.320.381
ES61.9880.0561.9322.012
ES72.1010.0231.9982.174
ES81.0010.2200.9991.035
GDP10.0300.0440.0270.035
GDP20.0250.0220.0240.027
GDP30.0110.0510.0120.017
GDP40.0150.0330.0140.018
GDP50.0120.3100.010.016
GDP60.0160.2560.0150.018
GDP70.0130.3450.0120.017
GDP80.0150.6170.0090.02
GDP90.0310.2400.0260.033
GDP100.0280.0250.0260.031
GDP110.0270.0520.0260.034
GDP120.0140.0380.0130.016
GDP130.0200.0290.0160.024
GDP140.0130.5150.0140.017
GDP150.0150.3620.0140.018
GDP160.0120.0280.0130.016
GDP170.0130.0530.0120.018
GDP180.0150.0360.0170.021
GDP190.0330.050.0320.038
GDP200.0130.20.0090.015
GDP210.0080.0450.0090.012
GDP220.0120.030.0110.019
GDP230.0120.030.010.018
GDP240.0070.010.0070.01
GDP250.0070.020.0060.012
GDP260.0040.030.0030.008
GDP270.0120.030.0110.015
GDP280.0090.030.0090.011
GDP290.0100.010.0090.013
GDP300.0090.020.0080.014
GDP310.0110.030.0090.016
* Note: ES1 to ES8 sequentially represent the elasticity of substitution at the first level and second level of the CES production function, respectively, in agriculture, manufacturing sectors, construction sectors and tertiary industry; GDP1 to GDP31 represent changes in the rate of GDP in the 31 provinces of China, respectively.
From the simulation results, under the condition of two percent growth in all production factor inputs, the growth rate of GPD in the 31 provinces of China is about 0.7 to 3.4 percent, with the confidence interval at 0.5 to 4 percent. This illustrates that the fluctuation of the elasticity of substitution has a small impact on the simulation results; therefore, the simulation results are robust.

3. Empirical Analysis and Results

3.1. Scenario Design

According to the draft of the 13th Five-Year Plan of China Development reported by the National Development and Reform Commission, three scenarios were designed for analyzing the relationships between industrial structure and land use structure, which were named the business-as-usual (BAU) scenario, the rapid economic growth (REG) scenario and the coordinated environmental sustainability (CES) scenario.
Here, we chose the population growth rate, investment growth rate and total factor productivity (TFP) growth rate as the variables to set the scenarios for future socio-economic development (Table 5). Firstly, the BAU scenario is set according to the forecast of the economic development in China, which could act as the baseline to be compared with other scenarios. In this scenario, the rate of labor force growth is obtained based on the projection of population and workforce age. The investment growth rate is set to decrease with the increase of economic development. We assume that the growth rate of TFP will continue this pattern and remain around 2% in 2011–2020 [15]. On this basis, we further develop the CES and REG scenarios according to the opportunities and challenges that China might be facing, together with the environmental conservation policies in the process of development. In the CES scenario, we conventionally assume that the population growth pressure and investment growth rate are both slightly lower than that in the baseline scenario [16]. In the REG scenario, we conventionally assume that the population growth rate is relatively high, that the investment growth rate is higher than that in the BAU scenario and that the TFP growth rate is higher than that in the BAU scenario.
Table 5. Socio-economic development scenarios in China, 2011–2020. BAU, business-as-usual; TFP, total factor productivity; REG, rapid economic growth.
Table 5. Socio-economic development scenarios in China, 2011–2020. BAU, business-as-usual; TFP, total factor productivity; REG, rapid economic growth.
ScenarioIndicatorGrowth Rate (%)
BAULabor0.28
Capital8.40
TFP2.00
CESLabor0.30
Capital8.05
TFP1.90
REGLabor0.25
Capital8.70
TFP2.10

3.2. Industrial Structure Change Projection

3.2.1. Industrial Structure Change in China

By applying the CGE model, the change rate of the output value of various industries in the thirty-one provinces was simulated under three scenarios in China. Based on the output value and the simulated change rate in 2010, the output values in 2010–2020 were generated. Additionally, the output values’ contribution was calculated as shown in Figure 2.
Figure 2. Analysis results of industrial structure change in China from 2010 to 2020 under scenarios of BAU, CES, and REG.
Figure 2. Analysis results of industrial structure change in China from 2010 to 2020 under scenarios of BAU, CES, and REG.
Sustainability 06 09343 g002aSustainability 06 09343 g002b
The simulated results showed that, in the three scenarios in 2010–2020, the secondary and service industries take the largest proportion. The changing trend of the industrial structure is similar (Table 2), with agriculture and the service industry transferring to secondary industry. The proportion of agriculture is shrinking, while the trend is stable. The proportion of the secondary industry shows an obvious increase. The proportions of the construction and service industry are decreasing.
The change rate of the industrial structure in the next ten years in China was calculated based on the proportion of various industries from 2010 to 2020 (Table 6). The proportion of the secondary industry under all three scenarios shows an increasing trend, while other industries are decreasing, with the service industry decreasing the most. Among the three scenarios, in REG, all industries, other than the secondary one, decreased the most. The range of BAU is smaller than that in REG, while the CES possesses the lowest range.
Under the BAU scenario, the proportion of the crop planting industry will decrease by 1.63% from 2010 to 2020, while animal husbandry will decrease by 1.89%; the other agricultural sectors will decrease by 0.73%; secondary industry will decrease by 23.63%; the construction and service industries will respectively decrease by 3.67% and 15.71%. Under the CES scenario, crop planting, animal husbandry and other agriculture will respectively decrease by 1.52%, 1.81% and 0.73%; while the secondary industry will increase by 22.63%. The construction and service industries will respectively decrease by 3.6% and 14.98%. Under the REG scenario, crop planting, animal husbandry and other agriculture will respectively decrease by 1.71%, 1.95% and 0.76%; while the secondary industry will increase by 24.54%, and the construction and service industry will respectively decrease by 3.76% and 16.37% from 2010 to 2020.
Table 6. Change rate of industrial structure in China under three scenarios, 2010–2020 (%).
Table 6. Change rate of industrial structure in China under three scenarios, 2010–2020 (%).
BAU
TimeCropWoolOther AgricultureSecondary IndustryConstructionService Industry
2010–2011−0.08−0.13−0.061.87−0.44−1.16
2011–2012−0.10−0.15−0.062.00−0.43−1.26
2012–2013−0.12−0.16−0.072.12−0.42−1.35
2013–2014−0.13−0.18−0.072.24−0.41−1.45
2014–2015−0.15−0.19−0.082.36−0.39−1.55
2015–2016−0.17−0.20−0.082.47−0.37−1.64
2016–2017−0.19−0.21−0.082.56−0.35−1.73
2017–2018−0.21−0.22−0.082.63−0.32−1.80
2018–2019−0.23−0.22−0.082.68−0.29−1.86
2019–2020−0.25−0.22−0.082.71−0.25−1.92
CES
TimeCropWoolOther AgricultureSecondary IndustryConstructionService Industry
2010–2011−0.08−0.12−0.061.80−0.42−1.12
2011–2012−0.09−0.14−0.061.92−0.42−1.20
2012–2013−0.11−0.16−0.072.03−0.41−1.29
2013–2014−0.13−0.17−0.072.14−0.40−1.38
2014–2015−0.14−0.18−0.072.25−0.38−1.47
2015–2016−0.16−0.20−0.072.35−0.36−1.56
2016–2017−0.18−0.20−0.082.44−0.34−1.64
2017–2018−0.20−0.21−0.082.51−0.31−1.72
2018–2019−0.21−0.21−0.082.56−0.28−1.78
2019–2020−0.22−0.22−0.092.63−0.28−1.82
REG
TimeCropWoolOther AgricultureSecondary IndustryConstructionService Industry
2010–2011−0.09−0.13−0.061.94−0.45−1.20
2011–2012−0.10−0.15−0.072.07−0.45−1.30
2012–2013−0.12−0.17−0.072.20−0.44−1.40
2013–2014−0.14−0.19−0.072.33−0.42−1.51
2014–2015−0.16−0.20−0.082.46−0.41−1.61
2015–2016−0.18−0.21−0.082.57−0.38−1.71
2016–2017−0.20−0.22−0.082.66−0.35−1.80
2017–2018−0.22−0.23−0.082.74−0.32−1.89
2018–2019−0.24−0.23−0.082.78−0.29−1.94
2019–2020−0.26−0.22−0.082.80−0.25−1.99

3.2.2. Industrial Structure Change in Each Province

Based on the analysis of the industrial structure changes in China, we further select Beijing, Guangdong, Heilongjiang, Hubei, Qinghai and Jiangxi as the representatives of the thirty-one provinces, to analyze the industrial structure. (Figure 3).
In 2010, under all three scenarios, the secondary and service industry are in an absolutely advantageous state. Among the six provinces, the proportion of service industry in Beijing is the largest. The secondary and service industries in Guangdong account for more than 80%. In Heilongjiang, the crop planting industry takes the largest proportion. The proportion of crop planting industry in Hubei is relatively higher, while the secondary industry proportion is lower, and the proportion for the service industry is higher than that of other provinces other than Beijing. In Qinghai province, the proportion for animal husbandry is relatively higher, while in Jiangxi province, the proportion of other agriculture and construction industries are higher.
Figure 3. Scenario analysis results of industrial structure change in the six typical provinces from 2010 to 2020 under scenarios of BAU, CES, and REG.
Figure 3. Scenario analysis results of industrial structure change in the six typical provinces from 2010 to 2020 under scenarios of BAU, CES, and REG.
Sustainability 06 09343 g003
In 2020, under the three scenarios, the secondary and service industries are still in an absolutely advantageous state. Compared with 2010, in every province, there are obvious increases in the proportion of the secondary industry, obvious decreases in the service industry, and slight increases in crop planting and slight decreases in animal husbandry, other agriculture and the construction industry.
During 2010 to 2020, the change rates of the industrial structure are different under different scenarios. (Table 7). Under the BAU and REG scenarios, the proportions of secondary industry are all increasing in the six provinces, but the proportions of other industries are decreasing. Under the CES scenario, the proportions for crop planting and secondary industry are increasing, while the proportions for other industries are decreasing.
Table 7. Change rate of industrial structure in six provinces of China under three scenarios, 2010–2020 (%).
Table 7. Change rate of industrial structure in six provinces of China under three scenarios, 2010–2020 (%).
BAU
RegionCropWoolOther AgricultureSecondary IndustryConstructionService Industry
Beijing−0.06−0.21−0.0525.16−3.06−21.78
Guangdong−0.91−0.85−0.6921.11−2.16−16.49
Heilongjiang−2.50−2.65−0.4423.05−3.53−13.93
Hubei−2.13−2.56−1.2026.34−3.73−16.72
Qinghai−1.75−3.68−0.1029.95−6.23−18.20
Jiangxi−2.39−2.78−1.9928.05−6.78−14.11
CES
RegionCropWoolOther AgricultureSecondary IndustryConstructionService Industry
Beijing1.71−0.33−0.0837.91−4.05−35.17
Guangdong4.56−1.19−0.9324.34−2.67−24.10
Heilongjiang11.70−3.85−0.6118.73−4.40−21.57
Hubei14.89−4.02−1.7722.57−4.76−26.92
Qinghai5.77−4.98−0.1331.30−7.45−24.52
Jiangxi11.39−4.11−2.7925.07−8.43−21.14
REG
RegionCropWoolOther agricultureSecondary IndustryConstructionService Industry
Beijing−0.06−0.22−0.0526.48−3.14−23.00
Guangdong−0.95−0.88−0.7121.93−2.21−17.18
Heilongjiang−2.61−2.74−0.4523.91−3.61−14.49
Hubei−2.25−2.66−1.2427.42−3.82−17.45
Qinghai−1.83−3.80−0.1030.92−6.35−18.83
Jiangxi−2.51−2.88−2.0529.05−6.93−14.67
Under the BAU scenario, the growth rate of secondary industry in Qinghai is the highest by 29.95 percent among all provinces. The recession rates for crop planting, animal husbandry and other agriculture in Beijing are the lowest. While in Heilongjiang, Qinghai and Jiangxi, crop planting, animal husbandry and other agriculture will decrease the most. In Jiangxi, the construction industry will decrease the most at 6.78%; while in Beijing, the service industry will decrease the most at 21.78%. Under the CES scenario, Hubei and Beijing will have the largest increase in crop planting and secondary industry. Qinghai and Jiangxi will have the largest decrease in animal husbandry and other agriculture. Jiangxi and Beijing will have the largest increase in the construction and service industries. Under the REG scenario, the secondary industry will increase the most in Qinghai at 30.92%, and the least increase will be in Guangdong at 21.93%. The provinces that will have the largest decrease in crop planting, animal husbandry and other agriculture remain in a similar situation under the BAU scenario, and those provinces that have the largest decrease in the construction and service industry remain in a similar situation under the CES scenario.
From the above analysis, we can see that the distinguished structures of the advantageous industries in each province have been influenced by different scenarios. Under the CES scenario, the secondary industry will increase relatively faster, but will decrease faster in the advantageous industries. Under the REG scenario, the myth of the advantageous industries will fade with an increasing trend in the secondary industry.

3.3. Land Use Change Projection

3.3.1. Land Use Changes in China

In this land use prediction model, we calculated land use changes in 31 provinces in China. By adding up the total amount of changes in each province, we analyzed the characteristics of the trend of land use changes in the next ten years.
The three scenarios show a similar trend with respectively changing variance. Land use under the REG scenario shows the largest variance compared to the other two scenarios; the second largest is under the BAU scenario, and the lowest variance is under the CES scenario. From 2010 to 2020, the growth rate of land use for manufacturing and urban construction is decreasing, while the growth rate for cultivated land, grassland and other land use for agriculture is decreasing. From Figure 4, until 2020, the land use structure will remain balanced, since increasing land uses can compensate for decreasing land uses.
Figure 4. Land use changes in China from 2010 to 2020 based on three scenarios of BAU, CES and REG.
Figure 4. Land use changes in China from 2010 to 2020 based on three scenarios of BAU, CES and REG.
Sustainability 06 09343 g004
From the Figure 5, we see that cultivated land, grassland and land use for other agriculture are increasing, while land use for manufacturing and urban construction are decreasing year by year. However, each type of land use is decreasing in different amounts.
Figure 5. The trends of six types of land use under the three scenarios in China 2010–2020.
Figure 5. The trends of six types of land use under the three scenarios in China 2010–2020.
Sustainability 06 09343 g005
From 2010 to 2020, cultivated land use will shrink at a decreasing rate. Under the economic growth-oriented scenario, cultivated land will shrink relatively more than under the other two scenarios, which is about 119.3 thousand ha. Next is under BAU, which is about 113.9 thousand ha. Under the CES scenario, cultivated land use will shrink 108.3 thousand ha. Grassland changes with the same trend as cultivated land. Under the REG scenario, it will shrink the most: about 147.6 thousand ha. Under the BAU scenario, this number is about 139.1 thousand ha. Under the CES scenario, it will shrink 129.7 thousand ha. Other land uses for agriculture are also decreasing. Under the same sequence of scenarios, it will decrease 14.3 thousand ha, 13.6 thousand ha and 12.9 thousand ha, respectively. However, the land use for manufacturing is annually increasing, but with a decreasing growth rate. It will increase 54.2 thousand ha, 51.3 thousand ha and 48.3 thousand ha, respectively, under the REG scenario, BAU scenario and CES scenario. The land use for urban construction will also increase 137 thousand ha, 129.5 thousand ha and 121.4 thousand ha under the same sequence of scenarios. The land use for other built-up areas will also increase, even more than the land use for manufacturing, about 89.9 thousand ha, 85.8 thousand ha and 81.3 thousand ha, respectively.
According to the national land use plan in 2006–2020, until 2010, total cultivated land in China was about 121 million ha, total land use for urban construction was about 0.3374 million ha and, until 2020, the total amount of cultivated land will not be less than 120.33 million ha; the total land use for urban construction will reach 37.24 million ha [17]. Based on this modified CGE model, we predict the total land use for cultivation and urban construction under the three scenarios from 2010 to 2020 in Table 8. Based on Table 8 and the constraint data of land use for cultivation and urban construction in 2010, we predict the total land use for cultivation and urban construction in 2020 in Table 9.
Table 8. Changes of land use for cultivation and urban construction under the three scenarios in China 2010–2020 (10 thousand ha).
Table 8. Changes of land use for cultivation and urban construction under the three scenarios in China 2010–2020 (10 thousand ha).
ScenariosCultivationUrban Construction
BAU−11.3926.66
CES−10.8325.09
REG−11.9328.11
Table 9. Estimated land use for cultivation and urban construction under the three scenarios in China 2010–2020 (10 thousand ha).
Table 9. Estimated land use for cultivation and urban construction under the three scenarios in China 2010–2020 (10 thousand ha).
ScenariosCultivationUrban Construction
BAU12,108.613400.66
CES12,109.173399.09
REG12,108.073402.11
Until 2020, cultivated land will be over 120.33 million ha, so that the target for cultivated land protection will be reached; while the target of 37.24 million ha for urban construction will not be reached under the three scenarios. Therefore, although cultivated land will decrease, it is in the acceptable range. Until 2020, food security will be protected under the total area of cultivated land at 1.8 billion mu. With economic growth, urban construction will require more land, but its expansion will be acceptable, as well, and still have the possibility of rising.

3.3.2. Land Use in Each Province

We selected Beijing, Guangdong, Heilongjiang, Hubei, Qinghai and Jiangxi to analyze the trend of land use changes from 2010 to 2020 in China.
Figure 6 shows that the six types of land use changes in the six provinces will decrease during 2010–2020. Until 2020, land use in China will tend to be stable. Sequentially, under the REG scenario, the fundamental scenario and the CES scenario, the variance of the changes of land use is sequentially smaller than the previous scenario. From 2010 to 2020, six types of land use changes will decrease, until 2020, becoming relatively stable. Specifically, cultivated land in Heilongjiang, Qinghai and Jiangxi is increasing, while it will decrease in Beijing, Guangdong and Hubei. Meanwhile, Guangdong will be the one with the greatest decreases. Grassland in Heilongjiang will increase, while it will decrease in Beijing, Guangdong, Hubei, Qinghai and Jiangxi. Meanwhile, Guangdong will be the one with the greatest decrease again. Other land use for agriculture in Jiangxi will increase, while it will decrease in the other five provinces, and Heilongjiang will be the one with the greatest decrease. Land use for manufacturing will increase in all six provinces, but with a decreasing growth rate. Meanwhile, Guangdong will be the one with the greatest increase, but Qinghai will be the one with the lowest increase. Land use for urban construction will increase in the provinces, except Jiangxi, and Guangdong will be the one with the greatest increase. Other land for built-up areas will increase in all six provinces; Guangdong will be the one with the greatest increase, and Heilongjiang will be the one with the lowest increase.
Figure 6. Trends of six types of land use changes in six provinces in China 2010–2020.
Figure 6. Trends of six types of land use changes in six provinces in China 2010–2020.
Sustainability 06 09343 g006aSustainability 06 09343 g006b
From 2010 to 2020, as shown in Figure 7, land use changes in six provinces will be most influenced under the REG scenario, but the least under the CES scenario. Specifically, cultivated land will increase with the greatest amount in Jiangxi province, about 4360.86 ha, 4130.36 ha and 3869.73 ha, which are respectively under the REG scenario, fundamental scenario and CES scenario. Cultivated land will decrease the greatest amount in Guangdong province, about 10,138.04 ha, 9649.81 ha and 9122.97 ha, which are respectively under REG, BAU and CES. Grassland will only increase in Heilongjiang, about 697.11 ha, 663.85 ha and 629.22 ha, which are respectively under REG, BAU and CES. Grassland will decrease the greatest amount in Guangdong, about 2725.29 ha, 2594.73 ha and 2458.49 ha, which are respectively under REG, BAU and CES. Other land use for agriculture will only increase in Jiangxi, about 12,972.94 ha, 12,344.04 ha and 11,666.72 ha, which are respectively under REG, BAU and CES. Other land use for agriculture will decrease the greatest in Heilongjiang, about 4072.94 ha, 3886.75 ha and 3688.24 ha, which are respectively under REG, BAU and CES. Land use for manufacturing will increase in all six province. Guangdong will increase the greatest, about 6098.39 ha, 5777.86 ha and 5435.54 ha, which are respectively under the economic growth-oriented scenario, the fundamental scenario and environmentally-oriented scenario, while Qinghai will increase least, about 55.09 ha, 52.09 ha and 48.9 ha, which are respectively under REG, BAU and CES. Land use for urban construction will increase the most in Guangdong province, about 3543.28 ha, 3401.9 ha and 3245.29 ha, which are respectively under REG, BAU and CES. Comparing the other five provinces, other land use for built-up areas will increase the most in Guangdong, as well, about 6972.64 ha, 6647.32 ha and 6293.85 ha, which are respectively under REG, BAU and CES.
Figure 7. Trends of six types of land use changes in six provinces in China under three scenarios 2010–2020.
Figure 7. Trends of six types of land use changes in six provinces in China under three scenarios 2010–2020.
Sustainability 06 09343 g007aSustainability 06 09343 g007b

3.4. Relationship Analysis of Future Industrial Structure and Land Use Structure

Based on the land use amount in 2010 and the land use change from 2010 to 2020, the land use amount for 2020 was calculated. Then, based on the land use amount for 2010 and 2020, the land use structure of 2010 and 2020 and the structural change during 2010–2020 were calculated. Based on the data in the three scenarios, the change in the industrial structure and land use structure in China from 2010 to 2020 was analyzed. (Table 10).
Table 10. Industrial structure and land use structure change under three scenarios in China, 2010–2020 (Unit: %).
Table 10. Industrial structure and land use structure change under three scenarios in China, 2010–2020 (Unit: %).
Industrial Structure Change
ScenarioCropWoolOther AgricultureSecondary IndustryConstructionService Industry
BAU−1.63−1.89−0.7323.63−3.67−15.71
CES−1.54−1.73−0.7222.73−3.68−15.06
REG−1.71−1.95−0.7624.54−3.76−16.37
Land Use Structure Change
ScenarioCropWoolOther AgricultureSecondary IndustryConstructionOther Construction
BAU8.51−2.65−6.430.69−0.170.05
CES8.01−2.53−6.000.63−0.160.04
REG8.94−2.80−6.760.74−0.170.05
Among the three scenarios, REG has the largest range of industrial structure and land use change, followed by BAU, and CES has the lowest range. From the perspective of industrial structure adjustment, only the proportion of the industrial sector is increasing, while the proportions of other industries are all decreasing. Other industries have a tendency to transform to the secondary industry. From the perspective of land use structure change, the proportions of cultivated land, industrial and mining area and land for other construction are all increasing, while other types of land use are decreasing. Among them, the proportion of cultivated land increased the most, while the proportion of other farmland decreased the most. From the perspective of the relationship between the industrial structure and land use structure, though the proportion of the cultivated land is increasing, the output value of the planting industry, which is based on cultivated land, is decreasing. The proportion of grassland and animal husbandry are both decreasing, besides the reduction of the grassland proportion being larger than for animal husbandry. Land use for agriculture is decreasing in both structures, while the reduction of the proportion of land for primary agriculture is much larger than that for other agriculture. The proportion of the secondary industry is increasing in both structures, while the proportion of the secondary industry increased much more than that for land structure. The proportion of construction is decreasing in both structures, while the decrease in the industrial structure is much more than that for land use structure. The proportion of the service industry is decreasing in both structures, while it is decreasing for the land use structure, which is much less than for the industrial structure.
Through comparison analysis of the industrial structure and land use structure change, we found that the change of cultivated land and the crop planting industry, the service industry and other construction land is inconsistent. The proportions of cultivated land and other construction land are increasing, while the proportions of the output value in the crop planting industry and the service industry are decreasing, which to a certain extent shows that the land production efficiency of cultivated land and other construction land (or land for the service industry) may have the possibility of decreasing.

4. Discussion and Conclusions

In this paper, a multi-regional dynamic CGE model was built to analyze and project the relationship between industrial structure and land use structure in China. The constraints of land supply and the impact of natural hazards have been designed with mathematical modelling. By nesting the mathematical model into a multi-regional dynamic CGE model, this research makes an improvement to the standard CGE model. Furthermore, from the perspective of the transformation of economic development under the BAU scenario, REG scenario and CES scenario, possible outcomes in the future are simulated.
The results shows that: (1) From 2010 to 2020, the industrial proportion in each province and the whole of China has been increasing. The agriculture and service industries gradually shift to the secondary industry. The rate of economic transformation has a significant impact on the industrial structure. The lower the rate, the faster the increase of the proportion of the secondary industry. On the other hand, if the economic development changes at a high rate, the increase of the proportion of the secondary industry will be impeded to a certain extent. Meanwhile, the transformation of the economic structure has various impacts on the structures of the preponderant industries in each province. (2) From the general trend, the land use of the secondary and tertiary industries has been increasing, and the land use of agriculture has been decreasing during the projection period, both in China as a whole and in each province. However, no matter if it is increasing or decreasing, the range of the change is shrinking, until 2020, when the land use structure in both China and each province will reach a stable stage. (3) The transformation of economic development modes, however the rate, is beneficial for raising the industrial development and output value. The industrial structure showed an inclination toward the secondary and tertiary industries; in response, the land use pattern converted from agriculture to the secondary and tertiary industries. The higher the rate of transformation of the economic growth rate, the higher the rate of industrial transformation, and correspondingly, the secondary and tertiary industries expanded at a much higher rate. (4) The transformation of the economic structure will cause cultivated land to decrease in China; however, the amount of the decline is acceptable in the next ten years. By 2020, the amount of cultivated land will not exceed 1.8 billion mu, which is enough to achieve food security. In the next ten years, the transformation of the economic structure will promote the expansion of the built-up area within a certain limit. The ultimate limit of 2020 has not been reached, which means there is still space for built-up area expansion. (5) The changing trends of land use for agriculture and the service industry are inconsistent with those of the other corresponding industries. There is a possibility that the productivity of the cultivated and service industry land will decrease. More attention ought to be paid to increasing cultivated land and service industry land use efficiency to promote the production of the corresponding land use industries.
Based on the conclusion of this research, we provide the following policy suggestions:
(1)
Promote a smooth transformation of the economic structure during the rapid development.
Researchers found that, at represent, the secondary and service industries are the leading modes of economic growth in China. However, from the perspective of the development tendency, the increasing trend for other industries is higher than the service industry; besides, the service industry shows a slight hint of transformation to the secondary industry. Industrialization-oriented development is able to promote rapid economic growth, but too much emphasis on the secondary industry will create great pressure on the natural environment. In this circumstance, in order to achieve economic sustainable development, the Chinese government ought to actively adjust and accelerate the transformation of the economic structure toward the service industry as the leading industry; meanwhile, the secondary industry should be the major driving force.
(2)
Adjust the land management pattern, and boost the market-oriented distribution of land factor.
With the current land use management with land policy regulation as the major method, guaranteeing the stable supply of land resources, in the long run, with the increasing market economy, land marketization is an inevitable inclination. In the process of promoting land marketization, the open land market to the international market could be used as a reference. Noting that there are specific characteristics of China, land marketization will not be lost in the international market due to governmental behavior being highly involved. Since the government has interfered too much in land use, it is good to mitigate conflicts of allocation such that there is the appropriate release of rights for land use and management of free trade in the market within a well-regulated exchange system.
(3)
Strengthen the regional surveys of the changes in the industrial structure and land use structure.
Industrial structure and land use are interacting with each other; to do research on industrial structure that is compatible with economic development, we have to strengthen the study on land use structure and vice versa. The close relationship between the two remind us that in the future process of socio-economic development, regional survey ought to be strengthened. Besides, there are less land use data, and the updating of the data is slow, which influences the industrial structure analysis to a certain extent. Therefore, strengthening regional surveys of industrial structure and land use structure change is particularly important.

Acknowledgments

This research was sponsored by the National Natural Science Foundation of China (Grant No. 41171434) the National Natural Science Funds of China for Distinguished Young Scholar (Grant No. 71225005). The authors would like to thank Feng Wu from Beijing Normal University, Yingzhi Lin from China University of Geosciences (Wuhan), and Gui Jin from Hubei University for their contributions on revision suggestions and data clearing in this paper.

Author Contributions

All authors read and approved the final manuscript.

Appendix

A1. Bayesian Estimation

Bayesian estimation addresses the minimization of posterior expected value, which is derived from Bayesian probability statistics, which can be divided into discrete and continuous forms [18]:
Discrete Bayesian estimation: suppose a set of events A 1 , A 2 , , A n , is complete, to any latent event B, P ( B ) > 0,
P ( A | i B ) = P ( A i B ) P ( B ) = P ( A i ) P ( B | A i ) i = 1 n P ( A i ) P ( B | A i ) ,   i = 1 , 2 ,
P ( A | i B ) is the probability of Event A happening under Event B, being informative a priori, which is a posteriori probability and P ( A i ) is the probability of Event A happening whenever Event B is informative a priori, which is a priori probability.
The Bayesian theorem of continuous random variable x: set θ as a parameter to be estimated; the a priori probability distribution of θ is g(θ), X is a finite set of a priori information of the random variable x; the probability density function of X is f ( X | θ ) ; thus, g ( θ | X ) is defined as the a posteriori probability distribution of θ as follows:
g ( θ | X ) = f ( X | θ ) g ( θ ) f ( X ) = f ( X | θ ) g ( θ ) θ f ( X | θ ) g ( θ ) d θ
From Equations (25) and (26), the Bayesian theorem is a method to reach estimated a posteriori probability distribution g ( θ ) of estimated parameter θ derived from the prior probability distribution g ( θ | X ) .
From Equation (26), by taking the integral of θ f ( X | θ ) g ( θ ) d θ to θ , under a certain range, θ is a constant, if i.i.d. f ( X | θ ) is similar as the likelihood function L ( θ | X ) of the estimator θ ; the Equation (26) can be rewritten in a reduced form as follows:
g ( θ | X ) f ( x | θ ) g ( θ ) = L ( θ | X ) g ( θ )
If the i.i.d. prior probability distribution of the estimator θ is given, the a posteriori probability distribution g ( θ ) can be derived by Equation (27) [19]. Generally speaking, θ is a vector of multiple parameters. Thus, taking the integration of the probability density function will require calculating the high-dimensional integral problem. Additionally, it is usually much more complex to reach a solution. Therefore, to solve Bayesian estimation, we have to calculate the high-dimensional integral problem. The Markov chain Monte Carlo method (MCMC) is a commonly used method for solving this problem.

A2. Markov Chain Monte Carlo Methods (MCMC)

The MCMC method can be traced to the “Buffon needle problem” in the 18th century. It is a quantitative method by using probability theory based on a large number of trials.
For example, to solve a simple unary integral problem, such as x 0 x 1 f ( x ) d x , it can be interpreted as a problem to solve the geometrical area of a trapezoid with the frontier of f ( x ) in a two-dimensional plane that can be divided as an N small area by taking the derivative with respect to x at the N − 1 differential point; therefore, the value of every differential point at f ( x ) can be represented by the area of each tiny rectangle. When N is continuously increased, the number of N small tiny areas is accumulated to approach the total area of the trapezoid with the frontier of f ( x ) .
The Monte Carlo method is based on the method of taking the integral, but randomly selecting the differential point inside of the integral interval and reaching an average of the N value of the accumulative value of the total area of the trapezoid with the frontier of f ( x ) as follows:
A ( x ) = i = 1 N f ( x i ) N
By increasing randomized selection points, the average of the total area of the trapezoid with the frontier of f ( x ) is approaching the actual value.
If the randomized selection points i.i.d. obey a certain statistical probability distribution, the expectation of those points can be given by the following equation:
A ( x ) = x 0 x 1 f ( x ) φ ( x ) d x
where φ ( x ) represents the probability density function of x.
Overall, randomization of the selection point is a key issue in the MCMC method.
The MCMC can be used for solving the integration of a complex high-dimensional function, which is based on a stationary distribution π ( x ) derived by the Markov chain sampling process.
The MCMC method follows the steps [20]:
(1)
Select Markov chain sampling with a well-defined stationary probability distribution of p ( * , * ) , which is the transition probability (or transition kernel) [19];
(2)
Select an initial point of x ( 0 ) , and generate a series of x ( 1 ) , x ( 2 ) , , x ( n ) [20] based on Equation (1);
(3)
With respect to m times and enough n points, the expectation of any f(x) is as follows:
E ˜ π f = 1 n m t = m + 1 n f ( x ( t ) )
where, Step (1) is the key to the success of the MCMC method [18].
Definition 1: Set a set of random variables { x ( 0 ) , x ( 1 ) , x ( 2 ) , } , at a time of j ( j 0 ) ; x ( j + 1 ) at j + 1 only depends upon the current status of j, but is unrelated to priori { x ( 0 ) , x ( 1 ) , , x ( j 1 ) } , where p ( x ( j + 1 ) | x ( 0 ) , x ( 1 ) , , x ( j ) ) = p ( x ( j + 1 ) | x ( j ) ) is held, and this stochastic sequence is the so-called Markov Chain [18], p ( * , * ) , which is the transition probability (or transition kernel):
p ( x , x ) = p ( x x ) = p ( x ( j + 1 ) = x | x ( j ) = x )
Definition 2: If p ( x , x ) π ( x ) d x = π ( x ) , x Θ , π ( x ) is the stationary distribution of the transition probability, p ( * , * ) (or transition kernel). In other words, x ( 0 ) can be any value, after iteration; the marginal distribution of Markov chain { x ( j ) } is converging to the stationary distribution of π ( x ) [19].
Gibbs sampling method is one of the most commonly used methods to reach a stationary distribution of the Markov chain with the transition probability. The basic idea of the Gibbs sampling method is to construct an irreducible normal return of the Markov chain { x ( j ) } from a high-dimensional complex sampling to reach a stationary distribution of π through the conditional probability distribution of π [20]. When n is large enough, the distribution of x ( j ) approaches π, which is approximately the distribution of the total samples of π [19].
Definition 3: μ i = n = 1 n f i i ( n ) represents the average time to go back and forth to i, if μ i < ; it is called the positive return to the normal of i, and if μ i = , it is called the zero return to the normal of i.
Definition 4: If in a conditional distribution of p ( x J | x J ) , all variables x J = { x i , i J } , x J = { x i , i J } are involved as variables or conditions, we call this condition distribution a full conditional distribution [21].
The Gibbs sampling method is as follows:
Define the initial value of the given arbitrary variable, x ( 0 ) = ( x 1 ( 0 ) , x 2 ( ( 0 ) ) , , x m ( 0 ) ) , at observation time t = ( t 1 , t 2 , , t n ) in a testing time period of M; the following statements are:
(1)
Sampling selection from the full conditional distribution:
  • select x 1 ( j + 1 ) from π ( x 1 ( j + 1 ) | x 2 ( j ) , , x m ( j ) , t ) ;
  • then, select x i ( j + 1 ) from π ( x i ( j + 1 ) | x 1 ( j + 1 ) , , x i 1 ( j + 1 ) , x i + 1 ( j ) , , x m ( j ) , t ) ;
  • select x m ( j + 1 ) from π ( x m ( j + 1 ) | x 1 ( j + 1 ) , , x m 1 ( j + 1 ) , t ) ;
(2)
Iterate (1) M times, and let ( x 1 ( 0 ) , x 2 ( 0 ) , x m ( 0 ) ) = ( x 1 ( M ) , x 2 ( M ) , x m ( M ) ) ;
(3)
Until the convergence of ( x ( j + 1 ) | x ( j + 1 ) = ( x 1 ( j + 1 ) , , x m ( j + 1 ) ) ) , and the average is derived.
x ^ 1 = x 1 ( M + 1 ) + x 1 ( M + 2 ) + + x 1 ( M + K ) K
x ^ m = x m ( M + 1 ) + x m ( M + 2 ) + + x m ( M + K ) K
Thus, a Markov chain is generated from a stochastic series of { x ( j ) } with a stationary probability distribution of π ( x | t ) , where x = ( x 1 , x 2 , , x m ) .

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xiubin, L. A review of the international researches on land use/land cover change. Acta Geogr. Sin. 1996, 51, 558–565. [Google Scholar]
  2. Liu, J.; Zhan, J.; Deng, X. Spatio-temporal patterns and driving forces of urban land expansion in China during the economic reform era. AMBIO A J. Hum. Environ. 2005, 34, 450–455. [Google Scholar]
  3. Chen, Y.; Li, X.; Tian, Y.; Tan, M. Structural change of agricultural land use intensity and its regional disparity in China. J. Geogr. Sci. 2009, 19, 545–556. [Google Scholar] [CrossRef]
  4. Bao, C.; Fang, C.; Chen, F. Mutual optimization of water utilization structure and industrial structure in arid inland river basins of Northwest China. J. Geogr. Sci. 2006, 16, 87–98. [Google Scholar] [CrossRef]
  5. Deng, X.; Zhang, F.; Wang, Z.; Li, X.; Zhang, T. An Extended Input Output Table Compiled for Analyzing Water Demand and Consumption at County Level in China. Sustainability 2014, 6, 3301–3320. [Google Scholar] [CrossRef]
  6. Hubacek, K.; Sun, L. Economic and Societal Changes in China and their Effects on Water Use a Scenario Analysis. J. Ind. Ecol. 2005, 9, 187–200. [Google Scholar] [CrossRef]
  7. Kuznets, S. Economic growth and income inequality. Am. Econ. Rev. 1955, 45, 1–28. [Google Scholar]
  8. Rose, A.; Liao, S.Y. Modeling regional economic resilience to disasters: A computable general equilibrium analysis of water service disruptions. J. Reg. Sci. 2005, 45, 75–112. [Google Scholar] [CrossRef]
  9. Gowdy, J.; Erickson, J.D. The approach of ecological economics. Camb. J. Econ. 2005, 29, 207–222. [Google Scholar] [CrossRef]
  10. Seung, C.K.; Harris, T.R.; Englin, J.E.; Netusil, N.R. Impacts of water reallocation: A combined computable general equilibrium and recreation demand model approach. Ann. Reg. Sci. 2000, 34, 473–487. [Google Scholar] [CrossRef]
  11. Jin, Q.; Jiang, Q.; Wu, F.; Li, X.; Deng, X. Ecological Risk Assessment of Benzo (a) pyrene in Yellow River Delta. CLEAN–Soil Air Water 2013, 41, 370–376. [Google Scholar] [CrossRef]
  12. Adkins, L.C.; Rickman, D.S.; Hameed, A. Bayesian estimation of regional production for CGE modeling. J. Reg. Sci. 2003, 43, 641–661. [Google Scholar] [CrossRef]
  13. Ramsey, J.B.; Zarembka, P. Specification error tests and alternative functional forms of the aggregate production function. J. Am. Stat. Assoc. 1971, 66, 471–477. [Google Scholar] [CrossRef]
  14. Sali, E.; Wolfson, H. Texture classification in aerial photographs and satellite data. Int. J. Remote Sens. 1992, 13, 3395–3408. [Google Scholar] [CrossRef]
  15. Chao, X.; Ren, B. The Fluctuation and Regional Difference of Quality of Economic Growth in China. Econ. Res. J. 2011, 4, 26–40. [Google Scholar]
  16. Liu, R.; An, T. Trend and Factor Analysis of Chinese Economic Growth Performance under Restrictions of Resource and Environment—A Research Based on a New Method of Productivity Index’s Construction and Decomposition. Econ. Res. J. 2012, 11, 34–47. [Google Scholar]
  17. The Central People’s Government of the People’s Republic of China. National Land Use Plan in 2006–2020; Xinhua News Agency: Beijing, China, 2008. [Google Scholar]
  18. Zellner, A. Bayesian estimation and prediction using asymmetric loss functions. J. Am. Statist. Assoc. 1986, 81, 446–451. [Google Scholar] [CrossRef]
  19. Tierney, L. Markov chains for exploring posterior distributions. Ann. Statist. 1994, 22, 1701–1728. [Google Scholar] [CrossRef]
  20. Gelfand, A.E.; Smith, A.F. Sampling-based approaches to calculating marginal densities. J. Am. Statist. Assoc. 1990, 85, 398–409. [Google Scholar] [CrossRef]
  21. Gilks, W. R. Markov Chain Monte Carlo. In Encyclopedia of Biostatistics; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2005. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Jin, Q.; Deng, X.; Wang, Z.; Shi, C.; Li, X. Analysis and Projection of the Relationship between Industrial Structure and Land Use Structure in China. Sustainability 2014, 6, 9343-9370. https://0-doi-org.brum.beds.ac.uk/10.3390/su6129343

AMA Style

Jin Q, Deng X, Wang Z, Shi C, Li X. Analysis and Projection of the Relationship between Industrial Structure and Land Use Structure in China. Sustainability. 2014; 6(12):9343-9370. https://0-doi-org.brum.beds.ac.uk/10.3390/su6129343

Chicago/Turabian Style

Jin, Qin, Xiangzheng Deng, Zhan Wang, Chenchen Shi, and Xing Li. 2014. "Analysis and Projection of the Relationship between Industrial Structure and Land Use Structure in China" Sustainability 6, no. 12: 9343-9370. https://0-doi-org.brum.beds.ac.uk/10.3390/su6129343

Article Metrics

Back to TopTop