Relationship Study on Land Use Spatial Distribution Structure and Energy-Related Carbon Emission Intensity in Different Land Use Types of Guangdong, China, 1996–2008

This study attempts to discuss the relationship between land use spatial distribution structure and energy-related carbon emission intensity in Guangdong during 1996–2008. We quantized the spatial distribution structure of five land use types including agricultural land, industrial land, residential and commercial land, traffic land, and other land through applying spatial Lorenz curve and Gini coefficient. Then the corresponding energy-related carbon emissions in each type of land were calculated in the study period. Through building the reasonable regression models, we found that the concentration degree of industrial land is negatively correlated with carbon emission intensity in the long term, whereas the concentration degree is positively correlated with carbon emission intensity in agricultural land, residential and commercial land, traffic land, and other land. The results also indicate that land use spatial distribution structure affects carbon emission intensity more intensively than energy efficiency and production efficiency do. These conclusions provide valuable reference to develop comprehensive policies for energy conservation and carbon emission reduction in a new perspective.


Introduction
Since energy-related carbon emission was recognized as the most important anthropogenic factor that leads to global warming [1,2], numerous studies have focused on the mechanism of how human activities affect energy consumption. Economic development, as well as industrial structure and energy consumption structure are typically considered as primary influencing factors for energy consumption [3][4][5][6][7]. Actually, human activity affects the regional carbon emission through changing the land use patterns, which in turn changes the energy consumption patterns and finally influences the amount and rate of carbon emission [8][9][10]. However, studies on land use change-induced carbon emission are usually directed toward calculating carbon storage in soil and vegetation cover given that soil respiration also considerably accounts for carbon dioxide emission into the atmosphere [11][12][13][14][15][16][17][18][19][20]. The social attribution of land use, which affects energy consumption via the spatial distribution of land use, is always neglected.
Related research focuses on the relationship between urban form and energy consumption [21,22]. The majority of such research suggests a strong negative correlation between urban density and energy consumption [23][24][25][26][27][28]. Concentration distribution reduces energy consumption, resulting in the emission of fewer pollutants and maintenance of sustainable urban development, particularly for traffic and residential lands [29][30][31][32][33][34][35][36][37]. This view is accepted by policy makers in many developed countries, and is considered as an effective way to reduce energy consumption [26,27]. Opposing perspectives, however, also exist. For example, Mindali et al. (2004) argues that the relationship between urban form and energy consumption is uncertain [26]. Holden and Norland (2005) indicated that decentralized concentration may reduce energy consumption in households [38]. Ma and Jin (2011) found a weak correlation between urban density and energy consumption under high population density in China [39]. Urban form, defined on the basis of the spatial distribution of land use, encompasses urban size, population density, and other related elements [40]. Except for traffic land and residential land distribution, other types of land use are seldom regarded as influencing factors for energy consumption in existing research. Thus, the existed urban form perspective is an unsuitable approach to determining the relationship between the distribution of different land use types and energy-related carbon emission. A comprehensive relation study between the spatial distribution structure of different land use types and energy-related carbon emission should be considered at a regional scale, specific to domestic conditions in terms of the large population, defective infrastructure, and binary structures of urban and rural areas.
Some researchers analyzed the energy-related carbon emission caused by land use patterns change in China, but none of them studied on their relationships.  found that the energy-related carbon emission per unit land area in Jiangsu Province increased from 8.24 t/hm 2 to 15.53 t/hm 2 from 2003 to 2007, and that the highest carbon emission intensities originated from residential and industrial land [8].  studied the energyrelated carbon emissions of different industrial spaces in China [41]. The results indicate that living and industrialcommercial spaces contribute the highest carbon emissions, and that the carbon emission intensities of each industrial space in Shanghai are the highest among the investigated 30 provinces. Using data envelopment analysis, You and Wu (2010) measured the provincial energy-related carbon emission efficiency of land use [42]. The results show that carbon emission per unit land area is effective in only three provinces. Zhang (2011) analyzed the land use constitution structure in Yunnan Province through information entropy method, and confirmed that the changes in land use structure cause energy-related carbon emission in the province [43].
Guangdong, which located in south China and comprises 21 cities (Figure 1), is one of the most developed provinces in the country. Its rapid development is attributed to the reform and opening policy implemented in 1978. Statistical data (taken from the Guangdong Statistical Yearbook and the Guangdong Agricultural Statistical Yearbook) show that Guangdong's GDP increased from 18.58 billion Yuan to 4601.31 billion Yuan from 1978 to 2010, and that its contribution to China's GDP increased from 5.1% to 11.47% during the same period. The proportions of output value in three major industrial sectors changed from 29.8% : 46.6% : 23.6% to 5% : 50% : 45%. The permanent population increased from 50.64 million to 104.40 million. In addition, land use spatial distribution and energy consumption significantly changed as human activities increased. Crop land decreased from 6.6416 million hectares in 1978 to 4.5245 million hectares in 2010, while terminal energy consumption increased from 18.13 million tons in 1980 to 263.4485 million in 2010. Guangdong was identified as one of the pilot provinces for national low-carbon economic development in 2010. To this end, the Guangdong government set a goal of reducing carbon dioxide emission per unit of GDP by 35% in 2015 relative to 2005 levels [44]. Land use pattern, a factor that significantly affects energy consumption, is also considered in the government's macro control policy. Therefore, ascertaining the relationship between land use spatial distribution structure and energyrelated carbon emission is important for formulating land use policies, which promote efficient energy conservation and emission reduction.
This study analyzes the relationship between land use spatial distribution structure and energy-related carbon emission intensity in Guangdong from 1996 to 2008. The core components of the research are as follows: (1) The To accurately calculate the energy-related carbon emission in each type of land use, we attempt to establish the corresponding relationships between industrial energy consumption and land use types according to the content of land use status classification in China [45]. Based on the research of Li (2009) and , land uses in this study are classified into five different types: agricultural land, industrial land, residential and commercial land, traffic land, and other land (includes lakes, rivers, and special designated land;

Spatial Lorenz Curve and Gini Coefficient.
Numerous methods are used to study the spatial distribution structure of land use. Such approaches include information entropy theory [47][48][49], geographic information system method combined with remote sensing imagery [50][51][52][53][54], landscape ecology method [55][56][57], and mathematical methods, such as Markov chain modeling and linear programming [58][59][60][61][62]. However, most of these methods are complicated and provide abstract results [63]. Spatial Lorenz curve analysis combines location entropy with the Lorenz curve concept in economics, making it an intuitive method for studying land use distribution structures [63][64][65]. The spatial Lorenz curve is quantified by the Gini coefficient, whose value indicates the concentration degree of land use spatial distribution. Thus, on the basis of relevant research, we use the spatial Lorenz curve and Gini coefficient to measure the concentration degree of the spatial distribution structures of different land use types.

Land Use Spatial Lorenz
Curve. The spatial Lorenz curve of land use can be established by citing two economic concepts, which are Lorenz curve and location quotient. The original Lorenz curve, which comprises population percentage and the corresponding income percentage, is an extensively used tool for examining national wealth distribution [66]. Lorenz curve bends indicate the concentration degree of national wealth and a straight line with no bending translates to equitably distributed wealth ( Figure 2) [67]. For intuitive distribution analysis, Lorenz curve is applied to many fields by using relevant percentage variables [68][69][70][71]. To measure the spatial distribution of different land use types in this study, the spatial Lorenz curve of land use should comprise specific land percentage and regional area percentage. The specific land percentage refers to the ratio of specific land use area in each city to the whole specific land use area in the province, while the regional area percentage denotes the ratio of each city area to the province area. Actually, the relationship between specific land percentage and regional area percentage can be defined by location quotient in economics [64,65]. location quotient is the ratio between a special factor percentage and the entire factor percentage. In this study, the special factor percentage is the specific land percentage, and the entire factor percent is the regional area percentage. These percentages are expressed as follows: where and are the indices for land use types and cities in Guangdong, respectively; is the location quotient of land use type ( = 1, 2, 3, 4, 5) in city ( = 1, 2, 3, . . . , 21) (Table 1 and Figure 1); represents the area occupied by land use type in city ; ∑ denotes the area occupied by land use type in Guangdong; represents the area occupied by city ; and ∑ is the area occupied by Guangdong Province. The numerator represents the area ratio of land use type in city , and the denominator represents the regional area percentage.
On the basis of the implications of Lorenz curve, we rank from the smallest to the largest; consider an ordinate value as the accumulative area ratio of different land use types and an abscissa value as the accumulative regional area percentage; and then plot the spatial Lorenz curve (Figure 3). The curve bends indicate the concentration degree of land use spatial distribution in the studied region. A straight line with no bending corresponds to equitably distributed land use.  The accumulative area ratio of land use type

Gini Coefficient.
The Gini coefficient measures the land use spatial distribution based on quantifies the bends of spatial Lorenz curve [72]. It is defined as where is the Gini coefficient, A represents the area between the equitable distribution curve and spatial Lorenz curve ( Figure 3, area A), A + B denotes the area under the equitable distribution curve (Figure 3, areas A and B).
The Gini coefficient reflects deviation from the equitable distribution curve. In the spatial Lorenz curve, therefore, a high Gini coefficient represents high degree of land use concentration.
We calculate the Gini coefficient in (2) by the curve regression method [73,74]. When the vertical and horizontal axes are both set to one unit, (2) is converted into where and are the indices for the abscissa value and ordinate value of the spatial Lorenz curve, respectively; thus, 0 < ≤ 1. According to the commonly used division standard of the Gini coefficient, land use distribution characteristics can be estimated as follows: when 0 < ≤ 0.2 stands for absolute decentralization, 0.2 < ≤ 0.3 represents decentralization; 0.3 < ≤ 0.4 indicates appropriate concentration; 0.4 < ≤ 0.5 means concentration; and 0.5 < ≤ 1 corresponds to absolute concentration [73].

Calculation of Energy-Related Carbon Emission.
Fossil energy accounts for more than 90% of total energy consumption. It is the main source of energy-related carbon emission in China [75]. Energy-related carbon emissions, which include the carbon content of both carbon dioxide and methane, are released from the consumption of fossil and biomass energy. To calculate the fossil energy-related carbon emission both from final consumption and transformational consumption, we separately calculate the carbon emission from final fossil energy consumption and that from terminal thermal power and heating consumption because the majority of fossil energy consumption in the energy conversion process is used to produce thermal power and heating (refer to the energy balance table for Guangdong in the China Energy Statistics Yearbook). The energy-related carbon emissions are calculated as follows [8,76]: where CT represents the total energy-related carbon emission of land use type (classified in Table 1); CR is the carbon emission from final fossil energy consumption, excluding the thermal power and heating supply for land use type ; CE and CH are the indices for the fossil energy-related carbon emissions during the supply of thermal power and heating, respectively, in land use type ; CB denotes the carbon emission from noncommercial energy, which is used by the residents of residential and commercial land.
The detailed calculation methods are as follows: The Scientific World Journal 5 In (5), CR (10 4 t) is the carbon emission from final fossil energy consumption, excluding thermal power and heating for land use type ; QRN (TJ) represents the combusted quantity of energy type ( refers to the fossil energy types in Table 2) in land use type ; QR (10 4 t (10 7 m 3 )) is the consumed fossil energy type in land use type ; NCV (TJ/Gg) denotes the net calorific value of fossil energy ; QRN is the product of QR and NCV ; EF (t/TJ) is the default carbon emission factor of fossil energy ; and (t/TJ) represent the default carbon contents in carbon dioxide and methane, respectively; and denotes the indices for default oxidation carbon factors; and 10 −3 is the coefficient of unit conversion.
The carbon emission from consumed biomass energy (CB ), which includes biogas, crop stalk, and firewood, can also be calculated by (5). The carbon emission factors are taken from the Provincial Greenhouse Gas List Preparation Guide (Table 2) [77].
In (6) and (7), CE (10 4 t) and CH (10 4 t) are the indices for the carbon emissions from thermal power consumption and heating supply, respectively, in land use type . The carbon emissions from thermal power and heating can be calculated by using standard coal consumption given that coal is the main source of thermal power and heating supply in China. QE (10 10 kWh) and QH (10 13 kJ) represent the quantities of thermal power consumption and heating supply, respectively, in land use type ; EE (g/kWh) and EH (kg/GJ) are the standard coal consumption levels per unit of thermal power and heating supply, respectively; and CF is the carbon emission coefficient of unit standard coal consumption. EE and EH are based on the China Power Statistical Yearbook (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009) (Table 3) [78]. CF equals 0.67 kg C/kgce [79].

Energy-Related Carbon Emission Intensity in Different
Land Use Types. The energy-related carbon emission intensity of land use pertains to the energy-related carbon emission in the unit area of land use [8]. It is calculated as follows: where CI (t/hm 2 ) represents the energy-related carbon emission intensity in land use type , CT (10 4 t) denotes the energy-related carbon emission in land use type ; and (10 4 hm 2 ) is the area occupied by land use type .

Regression Analysis and Cointegration Test.
Regression analysis is commonly used in research on relationships among time series variables. Because of the relevance of time series data, the unit root always exists in time series. Time series with unit roots are called nonstationary time series, and may considerably influence the results of a regression model [80]. Although some time series are nonstationary, certain linear combinations among such series are stationary. This linear combination reflects a stable long-term relationship among variables, and suggest that cointegration among the variables exist. If the unit roots of time series variables are uniformly integrated, the cointegration relationship can be verified by the Engle-Granger testing method [81,82]. Engle and Granger argue that the residual series of regression models should be stationary if cointegration among variables exists [83]. Thus, the cointegration relationship can be evaluated by determining the stationarity of residual series. Eliminating inaccurate regression in relationship research entails the following steps for regression analysis and cointegration testing: first, the stationarity of variables is determined by using Augmented Dickey-Fuller (ADF) values. The variables in a uniform integration order are then chosen to build a regression model. Finally, in accordance with the Engle-Granger method, ADF values are used to test the stationarity of the residual series in the regression model. If the residual series exhibit integration of order 1, then cointegration exists in the regression variables and the regression model is reasonable [83].
Studies on the energy-related carbon emission factors in China show that energy-related carbon emission is influenced by multiple factors, including industrial structure, energy consumption structure, energy efficiency, and production efficiency [84][85][86][87]. Different land use types are classified by industries and are investigated as separate research objects. These attributes indicate that the influence of structures, including industrial structure and energy consumption structure, are insignificant. Thus, GDP per unit energy consumption (denoted as GE) and GDP per unit capital investment (denoted as GC), as well as land use distribution structure (denoted as LS), are chosen as the independent variables in this study. GE and GC represent energy efficiency and production efficiency, respectively. Considering the effect of area changes on energy-related carbon emission, we choose 6 The Scientific World Journal

Spatial Distribution Structures of Different Land Use Types.
The spatial Lorenz curves of different land use types are shown in Figure 4. The spatial Lorenz curves of agricultural land are always the closest ones to the equitable distribution curve, whereas the industrial land curves are always the farthest. These characteristics indicate that agricultural land exhibits the most decentralized distribution among the five land use types, whereas industrial land shows the most concentrated distribution. Traffic land is more concentrated than residential and commercial land, as well as other land, over 1996 to 2008. The Gini coefficients are calculated based on the spatial Lorenz curve to quantify the distribution characteristics of land use. The results are shown in Figure 5 and Table 4. The average Gini coefficient of agriculture land was 0.0491, the smallest among the Gini coefficients of the five land use types. Although the Gini coefficient of agricultural land did not visibly increase in terms of absolute value, its average annual growth rate was higher than that observed for the other land use types. These results indicate that during the study period, agricultural land was absolutely decentralized, with a trend of concentrated use. Industrial land had the highest average Gini coefficient at 0.6163, indicating absolute concentration. Although minimal fluctuation occurs, the increasing trend of industrial land concentration was clearly observable, especially from 2002 to 2004. The average Gini coefficients of residential and commercial land and other land were 0.2999 and 0.2874, respectively, indicating that these two types of land use were distributed in decentralized manners, with different change trends. Residential and commercial land tended toward appropriate concentration, whereas other land tended toward decentralization. Traffic land exhibited appropriate concentration, with an average Gini coefficient of 0.3679; however, this value increased, indicating a concentrated trend for traffic land.

Energy-Related Carbon Emission Intensity in Different
Land Use Types. Using the methods discussed in Section 2, we calculate the energy-related carbon emissions in different land use types (Table 5, Figures 6 and 7). During the studied period, the total amount of energy-related carbon emission increased from 55.53 million tons to 159.33 million tons, with the annual increased rate of 9.18%. The energy-related carbon emissions in all land use types significantly increased in varying degrees. The carbon emission in industrial land increased from 32.79 million tons to 101.43 million tons. The proportion of total carbon emission in this land use type also increased from 59.06% to 63.66%, the fastest increase among all the land use types. On the basis of the increase in carbon emission, we rank the rest of the land use types as follows: residential and commercial land > traffic land > other land > agricultural land (Table 5). Although the carbon emissions in residential and commercial land and agricultural land increased 19.29 million tons and 0.49 million tons, respectively, their proportions of total carbon emission decreased by 5.86% and 2.47%, respectively. The carbon emissions in traffic land and other land increased 9.83 million tons and 5.56 million tons, respectively, and their proportions of total carbon emission increased by 2.02% and 1.71%, respectively. These results show that during the studied period, energy-related carbon emission intensively increased in all the land use types. The anthropogenic activities in the other three land use types more considerably increased than those in agricultural land and residential and commercial land.
The Scientific World Journal Energy-related carbon emission intensities are calculated according to (8). Results are shown in Table 6. For the period 1996 to 2008, the average energy-related carbon emission intensity in Guangdong increased from 3.2326 t/hm 2 to 9.2525 t/hm 2 , with an annual increase rate of 9.16%. The carbon emission intensities in industrial land, traffic land, and residential and commercial land were higher than the provincial level. In addition, the carbon emission intensities in agricultural land and other land were lower than the average values, and the carbon emission intensity in agricultural land was the lowest among all the land use types. The carbon emission intensities in various land use types increased with different ratios. On the basis of annual increase rates, we rank the land use types as follows: other land > traffic land > residential and commercial land > industrial land > agricultural land. The results indicate that the anthropogenic activities per unit area in each land use type increased, and that the continued development in other land caused a high annual increase rate of energy-related carbon emission in this land use type.

Regression Analysis and Cointegration Testing of Different
Land Use Types. Implementing the steps in Section 2.3, we test the stationarity of each variable using ADF values. Table 7 shows that the variables in the models of agricultural land, industrial land, residential and commercial land, and 8 The Scientific World Journal  other land are integrated in the same order. This integration indicates that stable long-term relationships exist among the variables in each of the four models. Regression models can therefore be constructed using these variables. In the traffic land model, the integration order of lnGC differs from those of the other three variables, indicating that a stable long-term relationship exists among the other three variables. Thus, the other three variables are considered in constructing the traffic land model. Five regression models are constructed, with their coefficients shown in Table 8. All of the regression models exhibit good fit, as indicated by the large R-square values. The models also pass the F-test at highly significant levels. Table 9 shows that the residual series of the regression models pass the stationarity test. These results indicate that cointegration exists among the variables in each model, and confirm the rationality of the five regressions models.
As indicated by the five models, land use distribution structure influences energy-related carbon emission intensities at different directions and extents. Aside from industrial land, the spatial distribution Gini coefficients of the other four land use types show a significant positive correlation with The Scientific World Journal 9 Notes: * , * * , * * * denote significance at the 10%, 5%, and 1% levels, respectively.  carbon emission intensities. These results indicate that concentrated land use distribution increases the carbon emission intensities in agricultural land, residential and commercial land, traffic land, and other land. Although the land structure variable is significant only at 20% level in the statistical t-test ( value, 0.1722), the model of industrial land also indicates that the Gini coefficient of spatial distribution is negatively correlated with carbon emission intensity in the long term. With respect to extent of influence, every 1% change in the spatial distribution Gini coefficients causes changes of 1.4794%, -1.4262%, 2.4352%, 3.8939%, and 9.7493% in the agricultural land, industrial land, residential and commercial land, traffic land, and other land, respectively. Among the land use types, other land is the type most strongly affected by spatial distribution structure. This result indicates that the other four types of land use, where intensive anthropogenic activities occur, are affected by other factors, such as energy efficiency and production efficiency. Meanwhile, the regression model of other land shows that the spatial distribution structure of land use strongly influences carbon emission intensity when no obvious industrial development and social activities occur in this land use type.
On the basis of the integrative effects of all the independent variables in the regression models (except for the other land model), we can analyze measures for reducing carbon emission intensity as follows. In the regression model of agricultural land, energy efficiency and product efficiency are negatively correlated with energy-related carbon emission intensity, whereas the concentration degree of spatial distribution is positively correlated with the dependent variable.

10
The Scientific World Journal Notes: * , * * , * * * denote significance at the 10%, 5%, and 1% levels, respectively.   Thus, improving the efficiency of energy use and production, as well as decentralize the agricultural land distribution can reduce the carbon emission intensity in agricultural land. The industrial land model shows that improving energy efficiency increases the carbon emission intensity in industrial land. This phenomenon also occurs in the residential and commercial land model, in which improving the efficiency of energy utilization and production increases carbon emission intensity (Table 8). Carbon emission per unit land use may increase when production per unit land use rises as a result of improved energy efficiency and production efficiency. Improving energy efficiency and production is the development trend driven by technological progress. Thus, controlling energy-related carbon emission through  land use distribution structures is crucial. The coefficients of the models indicate that centralize the industrial land and decentralize the residential and commercial land can reduce the carbon emission intensity, respectively. In traffic land, energy efficiency improvement and land use decentralization both can reduce energy-related carbon emission intensity.
The analyses above indicate that adjusting the spatial distribution structure of land use is important in reducing the energy-related carbon emission intensities in all land use types. The coefficients shown in Table 8 indicate that among all the variables in each model, the spatial distribution structure of land use has the largest effect on carbon emission intensity. Thus, focusing on efficient spatial distribution structure of land use can effectively controls carbon emission.

Conclusion
The spatial distribution structure of land use influences energy consumption patterns and is one of the most important factors for regional energy-related carbon emission. However, studies on the relationship between land use spatial distribution structure and energy-related carbon emission are rare. The current work discusses this relationship for different land use types in Guangdong, an economically advanced region in China and has been chosen as a pilot province for national low-carbon development. On the basis of land use pattern and industry types, land use is classified into five types which contain agricultural land, industrial land, residential and commercial land, traffic land, and other land. The spatial distribution structures of different land use types are analyzed by using the spatial Lorenz curve and Gini coefficient. The corresponding energy-related carbon emissions are calculated with factors taken from the Provincial Greenhouse Gas List Preparation Guide [77]. In constructing regression models, the energy-related carbon emission intensities in different land use types are set as dependent variables. Energy efficiency, production efficiency, and land use spatial distribution structure are chosen as independent variables. According to the results form regression analysis and cointegration testing, we draw the following conclusions.
(1) Aside from other land, the other four types of land use displayed different concentration trends during 1996 to 2008. Among the five types of land use, industrial land was the most concentrated, whereas agricultural land was the most decentralized but exhibits the fastest change rate.
(2) During the study period, energy-related carbon emission in each type of land use considerably increased. The anthropogenic activities in industrial land still account for the highest energy consumption. (3) Spatial distribution concentration is negatively correlated with energy-related carbon emission intensity in industrial land in the long term; by contrast, these two variables show a positive correlation in agricultural land, residential and commercial land, traffic land, and other land. These results indicate that the centralized distribution of industrial land may reduce the corresponding carbon emission intensity, but that the centralized distribution of the other four land use types can increase carbon emission intensity. (4) Land use spatial distribution structure has a more intensive effect on carbon emission intensity than energy efficiency and production efficiency do. Thus, adjusting the spatial distribution structure of land use is an effective way to control carbon emission. The findings provide valuable reference to develop comprehensive policies for regional energy conservation and carbon emission reduction at a new perspective. Meanwhile, some of the related suggestion can be summarized as follows.
(1) Strengthening the intensive use of agricultural land by improving the input-output efficiency of energy, labor, and capital. Ensure the spatial distribution of agricultural land keep on decentralizing when the land is occupied and compensated by other land types or adjusted in internal structure. (2) Promoting industrial agglomeration and unify the layout of energy-intensive industries coordinately. Meanwhile, technological innovation application, infrastructure construction improvement, and recycling economy should be encouraged in industrial development to make sure that resource utilization and scale economy effects are sufficiently developed. (3) Guiding the residents and commercial industry expand or transfer from population intensive region to the relative sparse region, avoiding the unnecessary energy waste from intensive energy consumption and the low commute efficiency caused by excessive concentration of residential and commercial land. (4) Increasing the density of road network in underdeveloped regions of Guangdong. Adjusting the transport instrument structure and industrial management standards, improve the traffic efficiency and reduce the unnecessary repeatedly energy cost by passengers and freight transport in the region of sparse road network. (5) Improving the development efficiency in other land by formulating specification management system. Trying to develop the other land evenly in spatial distribution instead of concentrated develop.

12
The Scientific World Journal  Although the results in this study are credible and representative, the study has certain limitations. Aside from the different spatial distributions of various land use types, the land use factors that affect regional energy consumption patterns also include land use composition. Thus, the result of this empirical study does not represent the common mechanism in other regions. In the future study, coordinating the internal relationship among land use factors and discussing their integrated effects on regional energy-related carbon emission are more comprehensive ways to define the mechanism of how the changes in land use influence regional energy-related carbon emission. Another limitation is the difficult access to land use classification data. Collecting land use area data on earlier years by land use classification survey is difficult, and accurately identifying different land use types through Landsat TM images with a 30 m resolution is very hard [27]. We anticipate the collection of more relevant land use data on earlier years, so that the argument in this study can be validated for longer time series.