Does the Exhaustion of Resources Drive Land Use Changes ? Evidence from the Influence of Coal Resources-Exhaustion on Coal Resources – Based Industry Land Use Changes

Analyzing the spatial-temporal changes of resources–based industrial land is essential to the transformation and development of resources–exhausted cities. In this paper, we studied coal resources–based industrial land use changes and their driving factors in a typical coal resources–exhausted city, Anyuan District, Pingxiang city. The changes between coal resources–based industrial land and other land-use types were analyzed. The logistic regression models were applied to identify the main driving factors and quantify their contributions to coal resources–based industrial land-use changes during the two periods of 2003–2008 and 2008–2013. The results show that coal resources–based industrial land declined by 34.37% during the period 2008–2013 as coal resources were being exhausted. Altitude, distance to roads, distance to town, population density change, fixed-asset investment per area change, and GDP per capita change drove coal resources–based industrial land-use changes. However, the patterns of the driving effects differed, and even the same factors had different influences on coal resources–based industrial land-use changes during the two periods. The changes in the driving factors can be seen as responses to socioeconomic transformation and development in the city, which is experiencing the exhaustion of coal resources. As a result of the comprehensive effects of these driving factors, coal resources–based industrial land use has changed in complex ways.


Introduction
Nowadays resources-exhausted cities whose economy used to rely on development of the natural resources such as mineral resources, forest resources and so on are facing the phenomenon of resource exhaustion.The cumulative extraction of resources in these cities has reached more than 70% of recoverable reserves and resource exploitation has come to the end stage [1].In recent years, resources extraction and initial processing enterprises in these cities had gradually closed, causing a dramatic change of industrial land.On the one hand, enormous abandoned industrial land brings environmental problems [2,3].On the other hand, there is an urgent need for construction land in economic transformation and sustainable development of these cities [4].However, few studies have explored coal resources-based industrial land-use changes and their driving forces.
In recent decades, dramatic land-use changes have occurred around the world, especially in rapidly developing regions [5,6].As an important component of global environmental change, land use and cover change (LUCC) has become a hot topic in sustainable development studies [7,8].Increasing numbers of studies have investigated land-use changes and their driving forces [9][10][11], including comprehensive land-use changes [12,13], construction land-use changes [14][15][16], farmland changes [17][18][19], and ecological land changes [20,21].Generally speaking, researchers believe that natural environment changes, industrialization and urbanization are the main driving forces in the conversion of farmland to nonagricultural use [22][23][24].Socioeconomic developments and land-use policies drive urban and industrial land expansion for future urban land-use spaces [25,26].On the other hand, the wetland and forest has increased and environment has improved with the ecological restoration projects are implemented [27].Most of these studies are focused on large-scale or mesoscale land use changes and their driving forces [28][29][30].At different spatial scales, the interactions between land use change and its associated factors are also different [31][32][33].An in-depth study of small-scale land-use changes can further reveal the micro-processes and internal mechanisms of regional spatial-temporal changes in land use [34].Moreover, as a special and important land use type in resources-exhausted cities, resources-based industrial land use changes are related to the phases of coal resources-based cities, which include the two different phases of prosperity and exhaustion.Analyzing the spatial-temporal changes of resources-based industrial land is essential to provide important information for its reasonable use and reclamation, and is beneficial to the transformation and development of these cities.Therefore, it is necessary to expand the research to coal resources-based industrial land-use changes and their driving forces at a small scale.
Pingxiang is a typical coal resources-exhausted city located in Jiangxi province in China.Anyuan District in Pingxiang was selected in this study to explore coal resources-based industrial land-use changes and their driving forces.By the using of high-resolution images from 2003, 2008 and 2013, changes in coal resources-based industrial land areas and their directions of flow were considered to determine spatiotemporal characteristics.The possible driving forces and their impact on coal resources-based industrial land-use changes were spatially identified and quantified using logistic regression models.

Overview of the Study Region
Anyuan District is situated in the hilly and mountainous area of central China, lying between 113 • 45 -113 • 59 E and 27 • 33 -27 • 44 N (Figure 1).It is the political, economic, and cultural center of Pingxiang in Jiangxi province.The landform is plain in the northwest part of the study area and mountainous in the southeast.The east-west width is about 22.75 km, and the north-south length is about 20.33 km, covering an area of approximately 212.58 km 2 .The area is located in a humid subtropical monsoon climate zone and is characterized by four distinct seasons with an average annual precipitation of 1630 mm and a mean annual temperature of 17.2 • C. The annual rainfall is higher than the annual evaporation.Anyuan District is among the earliest recorded coal mining areas in China and can be traced back to the Tang Dynasty.In the past, it was called the "Coal Capital" south of the Yangtze River because of its abundant coal resources.Coal resources-based industrial land is the main component of industrial land.As they are located in hilly and mountainous areas, coal mines are widely distributed at relatively high altitudes.The exploitation of coal resources has greatly contributed to the modernization of the region, as well as the whole country.As a result of high-intensity coal mining over a long period, coal resources in Anyuan District since 2000 have become nearly exhausted, resulting in a reduction of coal mining and coal processing activities.The coal economy is collapsing.In March 2008, the State Development and Reform Commission identified Pingxiang as a resources-exhausted city.As the most important coal production area of Pingxiang, Anyuan District is experiencing dilemmas in both its socioeconomic and ecological development.2).Slope and altitude data were extracted from ASTER GDEM data (30 m resolution) obtained from the Geospatial Data Cloud, Chinese Academy of Sciences (http://www.gscloud.cn).The data for distances to roads and towns were calculated using the ArcGIS.The socioeconomic data including urbanization rates, financial revenues, population, and fixed-asset investment were collected from Anyuan District's yearbooks for 2004-2014.These data were assigned to the sub-districts of Anyuan District's vector layer, and converted to raster layers with a resolution of 30 m [36].It is important to note that we take the changes of the data from 2003-2008, and from 2008-2013 into consideration, therefore, the study is actually based on cross-sectional data for two periods.In the study, coal resources-based industrial land refers to land used for coal mining, coal product processing and coal storage.The land use data of the corresponding periods and high-resolution images from Google Earth were used as reference data.At the same time, the coal resources-based industrial land and the changes were revised according to the investigation of mine management bureau and the village cadres of the area.The accuracy of the comprehensive evaluation of results was more than 86% [35].All images were converted into a grid map with 2.5 m resolution to reduce errors in counting the data.To obtain gains and losses in coal resources-based industrial land during the periods of 2003-2008 and 2008-2013, the land-use maps for 2003, 2008 and 2013 were overlaid using ArcGIS 10.3 software (Figure 2).Slope and altitude data were extracted from ASTER GDEM data (30 m resolution) obtained from the Geospatial Data Cloud, Chinese Academy of Sciences (http://www.gscloud.cn).The data for distances to roads and towns were calculated using the ArcGIS.The socioeconomic data including urbanization rates, financial revenues, population, and fixed-asset investment were collected from Anyuan District's yearbooks for 2004-2014.These data were assigned to the sub-districts of Anyuan District's vector layer, and converted to raster layers with a resolution of 30 m [36].It is important to note that we take the changes of the data from 2003-2008, and from 2008-2013 into consideration, therefore, the study is actually based on cross-sectional data for two periods.

Driving Factors Analysis of Industrial Land-Use Change
At any point on the earth's surface, land-use change can be divided into two types: Changed or not changed.Regression model was popular in the quantitative analysis of land-use change driving factors.However, the traditional linear regression model is not applicable for solving this problem with a dichotomous dependent variable [37][38][39].In this situation, the logistic regression model has the advantage of addressing this problem properly as it is a nonlinear model for binary or multivariate analysis [40].Therefore, logistic regression models are often used to analyze the driving forces in land-use changes [36,41,42].In this study, coal resources-based industrial land-use change was considered the dependent variable (Y) with a binary value (i.e., changed or not changed), linked to a series of possible driving factors (X), with the following logistic regression equation [43]: where ( ) Λ ⋅ represents the standard logistic cumulative distribution function.P is the occurrence probability of coal resources-based industrial land-use change, and β is a vector of parameters to be estimated.While the coefficients of logistic model only provide the direction and significance of explanatory variables, marginal effects which shows how much the event probability changes when a given explanatory variable is changed by one unit usually need to be computed.When one of the driving factors Xk is a continuous variable, its partial effect on coal resources-based industrial land-use change is obtained from the partial derivative: The explanatory variable will only pass the significance test if the associated p-value of the z-statistic is less than the given significance level (0.05 in our study).In addition, socioeconomic data including urbanization rates, financial revenues, population, and fixed-asset investment are based on Anyuan District's yearbooks for 2004-2014, which is measured at an aggregation level.However, the dependent variable is measured as a detailed spatial unit.The different measurement level of dependent variable and explanatory variables may lead to a bias in standard errors.Therefore, we use clustered adjusted standard errors to alleviate this problem [44].

Driving Factors Analysis of Industrial Land-Use Change
At any point on the earth's surface, land-use change can be divided into two types: Changed or not changed.Regression model was popular in the quantitative analysis of land-use change driving factors.However, the traditional linear regression model is not applicable for solving this problem with a dichotomous dependent variable [37][38][39].In this situation, the logistic regression model has the advantage of addressing this problem properly as it is a nonlinear model for binary or multivariate analysis [40].Therefore, logistic regression models are often used to analyze the driving forces in land-use changes [36,41,42].In this study, coal resources-based industrial land-use change was considered the dependent variable (Y) with a binary value (i.e., changed or not changed), linked to a series of possible driving factors (X), with the following logistic regression equation [43]: where Λ(•) represents the standard logistic cumulative distribution function.P is the occurrence probability of coal resources-based industrial land-use change, and β is a vector of parameters to be estimated.While the coefficients of logistic model only provide the direction and significance of explanatory variables, marginal effects which shows how much the event probability changes when a given explanatory variable is changed by one unit usually need to be computed.When one of the driving factors X k is a continuous variable, its partial effect on coal resources-based industrial land-use change is obtained from the partial derivative: The explanatory variable will only pass the significance test if the associated p-value of the z-statistic is less than the given significance level (0.05 in our study).In addition, socioeconomic data including urbanization rates, financial revenues, population, and fixed-asset investment are based on Anyuan District's yearbooks for 2004-2014, which is measured at an aggregation level.However, the dependent variable is measured as a detailed spatial unit.The different measurement level of dependent variable and explanatory variables may lead to a bias in standard errors.Therefore, we use clustered adjusted standard errors to alleviate this problem [44].

Selection of Driving Factors
Previous studies have shown that land-use change is a complex process that is affected by the comprehensive effects of natural environment changes and human activity [45][46][47].As a special land-use type, coal resources-based industrial land-use change is also affected by changes in the natural environment and human activity.After investigating the status of coal resources-based industrial land-use changes, physical geography conditions, location conditions, and socioeconomic conditions in Anyuan District-and combining them with principles of driving-factor selection such as representativeness, comprehensiveness, and data availability-9 possible driving factors were selected.However, financial revenue per area, which had a strong collinear relationship with other driving factors, was excluded.Thus, 8 independent variables were finally selected as possible driving factors (Table 1).In these driving factors, slope and altitude represent natural driving factors; distances to roads and towns represent location driving factors; and population density change, fixed-asset investment per area change, urbanization rate change, and GDP per capita change represent socioeconomic driving factors.Changes in socioeconomic factors were calculated in terms of the difference between the two periods of 2003-2008 and 2008-2013.The socioeconomic data were assigned to the towns of Anyuan District as attribute data in vector layers using ArcGIS 10.3.Finally, maps of the driving factors were converted to 30 m resolution grid maps to facilitate sampling.

Spatial Sampling
To avoid spatial autocorrelation, stratified random sampling was used by creating a random point tool in ArcGIS, and n observations were selected randomly throughout the study area.For every observation, the data for dependent and independent variables were extracted from coal resources-based industrial land-use changes and driving factor layers, and 200 sampling points were finally retained.For each sampling point, 1 was assigned if the changes between coal resources-based industrial land and other land-use types occurred during the two periods of 2003-2008 and 2008-2013; otherwise, 0 was assigned.Using the same observations, the data were extracted from the driving factor maps as the independent variable.To avoid influence on the constant term, the number of observations with values of 0 and 1 was equal in every model [48,49].Since the independent variables were selected from natural conditions, socioeconomic conditions, location, and so on, the sampled data have different dimensions and dimensional units.To eliminate the influence of the dimension, the sampled data were standardized using the standard scores standardization method using the following expression [36]: where X ' i,j stands for the standardization result; X i,j is the raw value of the ith possible driving factor in the jth sampling point; X i is the mean of the ith possible driving factor for all sampling points.Once all the data were ready, Stata 12 was used to construct the logistic regression models and analyze the driving forces of coal resources-based industrial land-use change.where X ' i,j stands for the standardization result; Xi,j is the raw value of the ith possible driving factor in the jth sampling point; i X is the mean of the ith possible driving factor for all sampling points.

Coal Resources-Based Industrial Land-Use Change
Once all the data were ready, Stata 12 was used to construct the logistic regression models and analyze the driving forces of coal resources-based industrial land-use change.During the period 2003-2013, the area of farmland, forest, and coal resources-based industrial land decreased.The area of general construction land and other land types increased, and the area of water remained mostly steady over the same period.Anyuan District is the urban area of Pingxiang City whose residential space is the main construction land.Although coal resources-based industrial land accounts for a small proportion of construction land, it accounts for a large proportion of industrial land.The total area of coal resources-based industrial land was 743.99 hm 2 in 2003, 498.55 hm 2  Of the total reduced coal resources-based industrial land, 38.15% was converted to general construction land (area: 179.60 hm 2 ), 32.37% was converted to forest (area: 152.36 hm 2 ), and 26.98% was converted to other land-use types (area: 127.03 hm 2 ).On the other hand, some general construction land, forest, and farmland were converted to coal resources-based industrial land.General construction land increased 25.64% from 5556.01 hm 2 to 6980.61 hm 2 whereas 109.70 hm 2 of general construction land was converted to coal resources-based industrial land.Farmland and forest were reduced by 29.31% (from 4975.91 hm 2 to 3517.61 hm 2 ) and 2.07% (from 9499.20 hm 2 to 9302.16 hm 2 ), respectively.Of the total reduced area in these two types, 66.82 hm 2 of forest and 36.61hm 2 of farmland were converted to coal resources-based industrial land.In model 1, altitude, distance to roads, distance to town, population density change, and fixed-asset investment per area change have significant effect on the loss of coal resources-based industrial land.Altitude has positive impact on coal resources-based industrial land loss: One unit increase in altitude will produce 0.111 units increase in the probability of coal resources-based industrial land loss.Given the difficulty of coal mining at higher altitudes, industrial land with the exhaustion of coal resources is more likely to be transferred.Distance to roads and distance to town as location factors positively affect coal resources-based industrial land loss: With one unit increase in these two factors, the probabilities of coal resources-based industrial land loss increase by 0.231 units and 0.095 units, respectively.Generally speaking, most of the coal mines are located in the hilly and mountainous regions far away from the town.The probability of coal resources-based industrial land loss increases indicates that more coal resources-based industrial land changes to other land use types.This is mainly because coal declined in quality during the period of coal-resources exhaustion and could not be sold at good prices, and long-distance transportation increases production costs.Population density change is another important driving factor that positively affects coal resources-based industrial land loss.Loss is more likely to occur in regions with higher population density.With a one unit increase in population density, the probability of coal resources-based industrial land loss increases by 0.153 units.Population agglomeration drives socioeconomic development.As a consequence, when the coal industry declines, labor forces and capitals that originally worked or invested in small and micro coal enterprises transfer to other industries for higher profits, which caused changes in coal resources-based industrial land use.It was found that many coal enterprises, especially small and micro private coal enterprises, shut down between 2003 and 2008.Further, the variable of fixed-asset investment per area change has negative impact on coal resources-based industrial land loss.With one unit increase in fixed-asset investment per area, the probabilities of coal resources-based industrial land loss decreased by 0.228 units.The reason for this could be that to ensure industrial output value, the government increased investment in the mining industry, since coal is the leading industry in the area.

Factors Driving Changes in
In model 2, only the variable of distance to town passes the significance test.Distance to town has a negative impact on coal resources-based industrial land gain.One unit increase in distance to town will produce a 0.122 units decrease in the probability of coal resources-based industrial land gain.The probability of coal resources-based industrial land gain declines indicates that it is more difficult to make the other land use types change to coal resources-based industrial land.With the coal-resources exhaustion, the difficulties of coal exploitation and the costs of production increased and as a result, the expanding of land for coal mining became more and more difficult.

Factors Influencing Industrial Land-Use Change, 2008-2013
Table 4 shows the results of the logistic regression models for the driving factors of coal resources-based industrial land losses and gains during the period 2008-2013.In model 3, altitude, distance to roads, distance to town, population density change, fixed-asset investment per area change and GDP per capita change significantly affect the loss of coal resources-based industrial land.Altitude has positive impact on coal resources-based industrial land loss: With a one unit increase in altitude, the probability of coal resources-based industrial land loss increases by 0.125 units.This suggests that coal resources-based industrial land at higher altitudes is more likely to be transferred.It was found that with the aggravation of coal-resource exhaustion during this period, coal shortages caused the exit of state-owned coal enterprises.Distance to roads has a negative impact on coal resources-based industrial land loss: With a one unit increase in distance to roads, the probability of coal resources-based industrial land loss decreases by 0.186 units.Meanwhile, distance to town positively affects coal resources-based industrial land loss: With one unit increase in distance to town, the probability of coal resources-based industrial land loss increases by 0.128 units.This is mainly because coal mining, processes, and storage cannot generate satisfactory profits in advantageous locations, and the coal resources-based industrial lands near roads are likely to change to other land types.At the same time, other coal mining enterprises that are far from towns are in continuous recession as coal resources become exhausted.As for socioeconomic factors, population density change, fixed-asset investment per area change, and GDP per capita change positively affect coal resources-based industrial land loss.With one unit increase in these three factors, the probabilities of coal resources-based industrial land losses increase 0.235 units, 0.242 units and 0.118 units, respectively.Along with socioeconomic transformation and development, the coal industry-previously the area's leading industry-is replaced by other industries, and capital, labor, and other production factors flow into the new leading industries with higher outputs.Therefore, these factors accelerate coal resources-based industrial land loss.
Similarly to Model 2, only the variables of distance to town passed the significance test at the 1% significance level in model 4. As with 2003-2008, distance to town negatively affected coal resources-based industrial land increase during this period.A one unit increase in distance to town will produce a 0.145 units decrease in the probability of coal resources-based industrial land gain.Note: *, **, and *** denote variables significant at the 10%, 5% and 1% levels, respectively.

Exhaustion of Coal Resources Improves the Land Environment Indirectly
The prosperity of mining activities will lead to a rapid reduction in the scale of ecological land such as forests and result in some ecological environment problems [50,51].After high-intensity coal mining over a long period, the coal resources had been exhausted from 2000 in the Anyuan District.At the first stage of 2003 to 2008, small and micro coal enterprises were rapidly shut down because of the lack of capital and labor, causing a dramatic decrease of coal resources-based industrial land.Although the state-owned coal enterprises possessing advanced production equipment and technology remained, the exhaustion of coal resources had influenced their production capacity, which led to a long-term decrease of coal resources-based industrial land.With the implementation of the policy of environmental improvement for abandoned mine, the reduced coal resources-based industrial land was converted to forest with the area of 152.36 hm 2 from 2003 to 2013, which went through a process from abandonment to recovery (Figure 4).With the depletion of coal resources and demand for urban transformation and development, the structure and function of the ecosystem have been restored under the help of self-repair and artificial restoration.Although the exhaustion of coal resources does not directly promote ecological restoration, the land's environmental conditions improve indirectly.Note: *, **, and *** denote variables significant at the 10%, 5% and 1% levels, respectively.

Exhaustion of Coal Resources Improves the Land Environment Indirectly
The prosperity of mining activities will lead to a rapid reduction in the scale of ecological land such as forests and result in some ecological environment problems [50,51].After high-intensity coal mining over a long period, the coal resources had been exhausted from 2000 in the Anyuan District.At the first stage of 2003 to 2008, small and micro coal enterprises were rapidly shut down because of the lack of capital and labor, causing a dramatic decrease of coal resources-based industrial land.Although the state-owned coal enterprises possessing advanced production equipment and technology remained, the exhaustion of coal resources had influenced their production capacity, which led to a long-term decrease of coal resources-based industrial land.With the implementation of the policy of environmental improvement for abandoned mine, the reduced coal resources-based industrial land was converted to forest with the area of 152.36 hm 2 from 2003 to 2013, which went through a process from abandonment to recovery (Figure 4).With the depletion of coal resources and demand for urban transformation and development, the structure and function of the ecosystem have been restored under the help of self-repair and artificial restoration.Although the exhaustion of coal resources does not directly promote ecological restoration, the land's environmental conditions improve indirectly.

Complexity of Coal Resources-Based Industrial Land Changes Driving Factors
It has been shown that the processes of land-use changes are closely related to environmental changes, population increases, socioeconomic developments, and land-use policies [24,52,53].

Complexity of Coal Resources-Based Industrial Land Changes Driving Factors
It has been shown that the processes of land-use changes are closely related to environmental changes, population increases, socioeconomic developments, and land-use policies [24,52,53].Natural factors control long-term and large scale land use change while socioeconomic factors play a leading role in short-term and small scale land use change [12].For different land use types, the driving factors and the driving mechanism of factors are different.As a special type of land use, coal resources-based industrial land changes have the both directions of increase and decrease in space in the depletion stage of coal resources in this area.The driving forces of social and economic development can also be seen as responses to socioeconomic transformation and development in the city, which is eager to get rid of the development dilemma caused by the exhaustion of coal resources.Apart from the social and economic development driving factors, the altitude is still an important factor influencing coal resources-based industrial land change in the hilly and mountainous areas.This may because particularity that coal production difficulty will increase with altitude.As a result of the comprehensive effects of these driving factors, the coal resources-based industrial land use changes in complex ways.

Bidirectional Process of Coal Resources-Based Industrial Land Changes
Land use change is a complex process in spatial and temporal dimensions, which could be drastic in a short time under the influence of human activities.In a period of time, the land use change of any point in space is bidirectional [23].Coal resources-based industrial land refers to land used for coal mining, coal product processing and coal storage, and the role of the driving factors of coal resources-based industrial land gain and loss may be different.Study on the spatiotemporal bidirectional process combining with investigation can help to explore the detailed driving mechanisms of land use change.

Influence of Policies on the Coal Resources-Based Industrial Land Changes
Taking various factors of physical geography, location and socioeconomic into consideration, the driving mechanism of coal resources-based industrial land changes is explained to some extent in this research.In addition to these driving factors, the policies of environment protection, ecological restoration, treatment of abandoned mines and transformation and development of resource exhausted cities can drive the coal resources-based industrial land changes in a short time.Pilot projects addressing waste land in industrial and mining have begun to be carried out in this area and obtained achievements in the improvement of the ecological environment.However, the policy impacts have not been taken into consideration, because the impacts of policies on coal resources-based industrial land changes are still hard to quantify by the logistic regression.The influence of policies on the coal resources-based industrial land changes should be examined in future studies if the data and methods are available.

Conclusions
Taking Anyuan District as a typical coal resource-exhausted city, this study analyzed coal resources-based industrial land-use changes and quantified the driving factors of physical geography, location, and socioeconomics using logistic regression models.Some conclusions can be drawn as follows.
SPOT images from 2003 and aerial photographs from 2008 and 2013 obtained from the Pingxiang Municipal Bureau of Land and Resources were chosen as the base images.The resolution of the SPOT images is 2.5 m, and the resolution of the aerial photography images is 1 m.The images from 2003, 2008 and 2013 were geo-referenced and imported into a Gauss-Kruger coordinate system.An accuracy of less than 0.5 pixel root mean square error was achieved using image rectification in ERDAS Imagine 9.2 software.The land-use categories were divided into six kinds: Coal resourcesbased industrial land, general construction land, farmland, forestland, and other use types by the objected-based classification method.In the study, coal resources-based industrial land refers to land used for coal mining, coal product processing and coal storage.The land use data of the corresponding periods and high-resolution images from Google Earth were used as reference data.At the same time, the coal resources-based industrial land and the changes were revised according to the investigation of mine management bureau and the village cadres of the area.The accuracy of the comprehensive evaluation of results was more than 86% [35].All images were converted into a grid map with 2.5 m resolution to reduce errors in counting the data.To obtain gains and losses in coal resources-based industrial land during the periods of 2003-2008 and 2008-2013, the land-use maps for 2003, 2008 and 2013 were overlaid using ArcGIS 10.3 software (Figure
SPOT images from 2003 and aerial photographs from 2008 and 2013 obtained from the Pingxiang Municipal Bureau of Land and Resources were chosen as the base images.The resolution of the SPOT images is 2.5 m, and the resolution of the aerial photography images is 1 m.The images from 2003, 2008 and 2013 were geo-referenced and imported into a Gauss-Kruger coordinate system.An accuracy of less than 0.5 pixel root mean square error was achieved using image rectification in ERDAS Imagine 9.2 software.The land-use categories were divided into six kinds: Coal resources-based industrial land, general construction land, farmland, forestland, and other use types by the objected-based classification method.

Figure 2 .
Figure 2. Changes in coal resources-based industrial land in Anyuan District.

Figure 2 .
Figure 2. Changes in coal resources-based industrial land in Anyuan District.

Figure 3
Figure 3 shows the area for each type of land use in Anyuan District from 2003 to 2013.To analyze coal resources-based industrial land-use changes, the land-use transition between coal resources-based industrial land and non-coal resources-based industrial land during 2003-2013 was calculated using a land-use change matrix.Table 2 shows the results of the land-use change matrix for 2003-2013 in Anyuan District.

Figure 3
Figure 3 shows the area for each type of land use in Anyuan District from 2003 to 2013.To analyze coal resources-based industrial land-use changes, the land-use transition between coal resources-based industrial land and non-coal resources-based industrial land during 2003-2013 was calculated using a land-use change matrix.Table 2 shows the results of the land-use change matrix for 2003-2013 in Anyuan District.

Figure 3 .
Figure 3. Area for each type of land use in Anyuan District during the period 2003-2013.

Figure 3 .
Figure 3. Area for each type of land use in Anyuan District during the period 2003-2013.
in 2008, and 488.30hm 2 in 2013, which represents a dramatic decrease during 2003 to 2008.
Industrial Land-Use 4.2.1.Driving Factors of Industrial Land-Use Change, 2003-2008 Table 3 shows the coefficients, clustered adjusted standard errors and marginal effects of the two logistic regression for the driving factors of coal resources-based industrial land-use losses and gains from 2003 to 2008.

( 1 )
From 2003 to 2013, the area of coal resources-based industrial land was significantly reduced with the exhaustion of coal resources.Coal resources-based industrial land was mainly converted to general construction land and forest.And some general construction land, forest, and farmland were converted to coal resources-based industrial land.(2) Factors of physical geography, location, and socioeconomics exerted varying degrees of impact on coal resources-based industrial land gains and losses during the periods of 2003-2008 and 2008-2013.From 2003 to 2008, distance to town was the main factor affecting coal resources-based industrial land gains, while altitude, distance to roads, distance to town, population density change, and fixed-asset investment per area change were the main factors affecting the loss of coal resources-based industrial land.From 2008 to 2013, altitude and distance to town drove coal resources-based industrial land gains, while altitude, distance to roads, distance to town, population density change, fixed-asset investment per area change, and GDP per capita change drove coal resources-based industrial land decreases.(3) Generally speaking, altitude, distance to roads, distance to town, population density change, and fixed-asset investment per area change were the main factors affecting the change of coal resources-based industrial land.Although the driving factors of coal resources-based industrial land gains and losses shared some similarities, the patterns of driving effects were different, and even the same factors had different influences on coal resources-based industrial land-use changes during different periods.(4) In the exhaustion process of coal resources, the transformation from coal resources-based industrial land into other types of land is the main trend.The land administration department and mine management department should make relevant plans in advance to ensure the orderly transformation of coal resources-based industrial land into other land use types.At the same time, the interests of residents around the mining area, grass-roots governments and other stakeholders need to be considered for sustainable development to occur.Although the study is based on evidence for county scale, it provides a reference for the study of land use change in hilly and mountainous coal resources-exhausted cities.

Table 1 .
Description and value assignment of variables.

Table 3 .
Results of the logistic regression models for coal resources-based industrial land-use changes, 2003-2008.

Table 4 .
Results of the logistic regression models for coal resources-based industrial land-use changes, 2008-2013.

Table 4 .
Results of the logistic regression models for coal resources-based industrial land-use changes, 2008-2013.