A Local Land Use Competition Cellular Automata Model and Its Application

Cellular automaton (CA) is an important method in land use and cover change studies, however, the majority of research focuses on the discovery of macroscopic factors affecting LUCC, which results in ignoring the local effects within the neighborhoods. This paper introduces a Local Land Use Competition Cellular Automata (LLUC-CA) model, based on local land use competition, land suitability evaluation, demand analysis of the different land use types, and multi-target land use competition allocation algorithm to simulate land use change at a micro level. The model is applied to simulate land use changes at Jinshitan National Tourist Holiday Resort from 1988 to 2012. The results show that the simulation accuracies were 64.46%, 77.21%, 85.30% and 99.14% for the agricultural land, construction land, forestland and water, respectively. In addition, comparing the simulation results of the LLUC-CA and CA-Markov model with the real land use data, their overall spatial accuracies were found to be 88.74% and 86.82%, respectively. In conclusion, the results from this study indicated that the model was an acceptable method for the simulation of large-scale land use changes, and the approach used here is applicable to analyzing the land use change driven forces and assist in decision-making.


Introduction
Land resources are the most fundamental and important resources in human life and production.They are not only related to the structure and function of ecosystems, but also affect the global interaction between the ecosystems and land surface environment [1].Different land use types have different social, economic, environmental, and ecological development functions, and exist in a state of dynamic equilibrium.Land use change is a complex process caused by the interactions between nature and social systems at several spatial and temporal scales [2,3].The diverse demand of society for limited land resources is the root of land use competition occurrences, and in essence it is a type of competing interest, which eventually creates mutual conversions between land use types [4].Dynamic changes in land resources are embodied in the competitive process between the various types of land use.It is a process that maximizes the satisfaction of the new demand of the land use by weighing and changing the land use patterns following the changes of the human demand for the land use [5].
Models characterized by simplifying and abstracting play a very important role in understanding and predicting the patterns and evolution process of LUCC.Currently, there are many valuable land use change simulation and prediction models and related research work, which typically include the system dynamics model (SD) [6], Markov model [7], conversion of land use and its effects at small regional extent (CLUE-S) model, multi-agent systems (MAS) model [8,9], cellular automata (CA) model [10][11][12][13], and the integration of these models [14].Many studies emphasize in particular a macro perspective of the driving factors to simulate land use changes.However, in fact, land use changes are often the results of the integrated effects of the natural and human macroscopic and microscopic factors.The question remains as to how to begin from the point of view of the micros of local land use competition, and then combine the macro driving factors and the micro effects for the purpose of building a reliable land use change simulation model, which has currently become a breakthrough worthy of studying.
Miller suggests that relationships among near entities do not signify a simple and sterile geography, and complex geographic processes and structures can emerge from local interactions [15].CA is a grid dynamics model, in which time, space, and state are all discrete, and the spatial interaction and temporal causality are both local.This model has been able to simulate the spatiotemporal evolution process of complicated systems [16][17][18][19].The state of each cell changes over time, and is determined by the local transition rules, where a large number of cells comprise the evolution of the dynamic systems through simple interaction [20].Wolfram proved that a CA model could simulate the complex natural phenomena, and the local behavior between individuals can evolve the global change patterns in time and space [21].Recently, CA models have been increasingly applied in the simulation and prediction of urban expansion and land use change research, and have obtained significant results, which indicate that a CA model can effectively reflect the complex features of land use evolution at the microscopic pattern [22][23][24][25].
It is a basic feature of a CA that global dynamic behavior can evolve from local rules.In land use models, CA can typically model the transition of a cell from one land use to another, depending on the land use within the neighborhood of the cell [26].Spatial variations in the relationships between driving forces and land use change are ignored in conventional land use CA models, which are based on assumption of spatial homogeneity [27].The existing CA models greatly weakened the fundamental characteristics of the original scheme of CA with the introduction of intelligent methods, and the continuous relaxation of the model constraints [28].The majority of the previous studies have focused on constraint design, influencing factors, and their intelligent acquisition, in order to improve the CA simulation results from a macro perspective.However, one of the most important characteristics of a CA is the local effect within a neighborhood.The lack of research on local effects and principles has lost the advantage of global change simulation through local interaction.One important method of solving this problem is to pay more attention to the study of the spatial relationship of the local scale, while introducing a number of macro geographical restrictions to the CA.Yang et al., (2014) designed a land use change simulation CA model, which focused on the cell type of land use conversion in multiple directions, and indicated that the transition probability of each cell could be influenced by not only the same land use cells, but also the different land use cells [29].Sirakoulis et al., (2012) designed a CA model based on an irregular neighborhood, which reflected a very good consideration of local interaction, and the transition potential for the cell could be defined by using the intensities of the signals received from the other cells within the neighborhood [30].
Land use change study not only requires both global and macro-scales, but also needs local and micro-scale to be cooperative [31,32].CA provides the spatial-temporal dynamic simulation computing framework by using a bottom-up approach, and it can achieve a land use change simulation under the comprehensive impact of macroscopic natural and social factors, as well as the interaction of the local neighborhood of cells [33][34][35].Based on the achievements listed above, in this study, a Local Land Use Competition Cellular Automata (LLUC-CA) model is introduced, which is based on local land use competition, land suitability evaluations, demand analysis of the different land use types, and a multi-target land use competition allocation algorithm, for the purpose of simulating the land use changes at a micro level, which focuses on the understanding of the neighborhood effects in the simulation process.In order to discuss the land use change simulation method based on a LLUC-CA, and to verify the reliability of the model, as well as to provide theoretical and technical support for the simulation and prediction of land use changes, a case study was conducted at the Jinshitan National Tourist Holiday Resort (hereinafter referred to as Jinshitan), which is a very local extent.

Study Area
Jinshitan, the back garden of Dalian, Liaoning Province, China (39 ˝11 49 2 -39 ˝81 10 2 N, 121 ˝56 1 11 2 -122 ˝41 42 2 E), is located in the Liaodong Peninsula in the northeastern area of Dalian (Figure 1 simulating the land use changes at a micro level, which focuses on the understanding of the neighborhood effects in the simulation process.In order to discuss the land use change simulation method based on a LLUC-CA, and to verify the reliability of the model, as well as to provide theoretical and technical support for the simulation and prediction of land use changes, a case study was conducted at the Jinshitan National Tourist Holiday Resort (hereinafter referred to as Jinshitan), which is a very local extent.

Study Area
Jinshitan, the back garden of Dalian, Liaoning Province, China (39°1′49″-39°8′10″N, 121°56′11″-122°4′42″E), is located in the Liaodong Peninsula in the northeastern area of Dalian (Figure 1).It faces the Yellow Sea on the east, and is composed of the eastern peninsula, western peninsula, open heartland, and the beach between the two peninsulas.It covers a land area of approximately 58 km 2 .Located adjacent to the Dalian urban area, it is a typical area with a rural-urban zone and a coastal region.Beginning as a rural town dominated by traditional agriculture, it was identified as a National Scenic Spot in 1988, National Geopark of China in 2004, and a 5A Tourist Area in 2010.Its construction and development into a tourist resort in China had a strong representation.For nearly 25 years, the study area has experienced rapid economic development and land use change, and therefore possesses high research value.

Data Sources and Processing
The years 1988, 2004 and 2010 are three time nodes that have important significance to the development of Jinshitan.In view of the data accessibility, this study based its research on the land use data of 1988, 2003 and 2012.The real land use maps were mainly obtained from the 1:10,000 map scale vector land use database.The administrative map, road, town center and coastline data were mainly obtained from the comprehensive underlying database.All of these data were provided by

Data Sources and Processing
The years 1988, 2004 and 2010 are three time nodes that have important significance to the development of Jinshitan.In view of the data accessibility, this study based its research on the land use data of 1988, 2003 and 2012.The real land use maps were mainly obtained from the 1:10,000 map scale vector land use database.The administrative map, road, town center and coastline data were mainly obtained from the comprehensive underlying database.All of these data were provided by the Dalian Municipal Bureau of Land Resources and Housing, including the topographic maps and remote sensing images.
By combining the land use vector data with the high spatial resolution remote sensing images, the respective maps of land use of Jinshitan in 1988, 2003 and 2012 were acquired by a series of processing, such as projection transformation, spatial adjustment, topology error checks, cutting, image interpretation, and data normalization.ArcGIS10.2and ERDAS Imaging 2013 software were used.In order to simplify simulation process, the land use types were reclassified into agriculture land, construction land, forestland, and water by ArcGIS.The first three land use types obtained a large proportion and changed more obviously, and water was an important natural resource that had an outstanding effect on Jinshitan (Table 1).All of the spatial data were standardized to the Xi'an 80 geographic coordinate system and Gauss Kruger zone dividing projected coordinate system, and the raster data was resampled with a grid size of 5 m ˆ5 m.

The Model
Through local interactions, a CA can simulate global changes, while the global computing capacity is achieved by local computation.In global models, it is assumed that all variables are the same.However, one possible situation is that global models cannot be expressed very effectively in a specific location.Therefore, an LLUC-CA model uses cells in different land use types within a neighborhood as the neighborhood scenario, and forms local transition rules through the combination of the local land use competition and the macro driving factors.With these rules in place, it can evolve the global and complex dynamic behavior, and then simulate all of the land use changes.

Local Land Use Competition
At the present time, the land use science community has not yet clearly defined the meaning of local land use competition.Before further research is conducted, a clear definition should be made regarding this scientific question.In physical geography, the term "local" may be some area over which a particular process has an obvious effect.However, in terms of spatial data analysis, a "local space" is often expressed as a certain distance of space starting from a specific point or an area, for example, the locality or neighborhood of that point or area [36].Local interactions among entities are capable of generating complex global behaviors, as well as intricate structures in space and time [15].In this study, "local land use competition (LLUC)" is defined as a mutual competition between the different land use types in the neighboring regions.Local land use competition materializes as the mutual competition between the different state of cellular within a neighborhood, with the interaction intensity of the cell decreasing with the increase of the distance between the cells, and the scope of the workspace embodied in the locality.Therefore, the competition result is the land use state change of the central cell.
Land use changes occur due to competition within and among different land use types.Spatial autocorrelation not only acts on the same land use type cells, but also on the different types of land use cells.This study assumed that there were three land use types, namely agricultural land (A), construction land (C) and forestland (F), in a 5 ˆ5 Moore neighborhood (radius r = 2), and the different land use type cells were denoted with different colors.Under the effects of local land use competition and macro driving factors, the cells in a certain state can transition to multiple directions, and the competing relationship between the land use types is illustrated with Figure 2. As shown in the figure, Num (A) = 7 expressed that the number of cells of the agricultural land in the neighborhood scenario is 7, Num (C) = 8 expressed that the number of cells of construction land in the neighborhood scenario is 8, and so on.Due to the fact that the land type transformation of the central cell was based on its independent neighborhood scenario, the interactions between the land use types were not equal.The competition relationship between the central and neighbor cells is expressed by the land use competition intensity.The land use competition intensity was constantly changing, it was a dynamic parameter, and was inspired by spatial interaction model, and the local land use competition intensity of the cell (i,j), which state (or land use type) is k, at time t can be described as: where N is the total number of cells in the central cell's neighborhood (in this study it is 25); J(S t (i+x,j+y)) is used for judging the state of cell (i+x,j+y), whether it is k, 1 for true, and 0 for false; d(i+x,j+y) is convolution kernel, and it represents the distance from the central cell (i,j) to neighborhood cell (i+x,j+y), if (i+x,j+y) is the central cell, d(i+x,j+y)=1; otherwise if (i+x,j+y) is just close to (i,j), then d(i+x,j+y)=2; or if (i+x,j+y) and (i,j) are separated by one cell, then d(i+x,j+y)=3.If the neighborhood cell is closer to the central cell, the impact on it will be greater.The competition relationship between the central and neighbor cells is expressed by the land use competition intensity.The land use competition intensity was constantly changing, it was a dynamic parameter, and was inspired by spatial interaction model, and the local land use competition intensity of the cell (i,j), which state (or land use type) is k, at time t can be described as:

Land Suitability and Transition Potential
rJ `st pi `x, j `yq ˘{dpi `x, j `yqs where N is the total number of cells in the central cell's neighborhood (in this study it is 25); J(S t (i+x,j+y)) is used for judging the state of cell (i+x,j+y), whether it is k, 1 for true, and 0 for false; d(i+x,j+y) is convolution kernel, and it represents the distance from the central cell (i,j) to neighborhood cell (i+x,j+y), if (i+x,j+y) is the central cell, d(i+x,j+y)=1; otherwise if (i+x,j+y) is just close to (i,j), then d(i+x,j+y)=2; or if (i+x,j+y) and (i,j) are separated by one cell, then d(i+x,j+y)=3.If the neighborhood cell is closer to the central cell, the impact on it will be greater.

Land Suitability and Transition Potential
The process of LUCC is very complex, and the suitability conditions play a very important role in the land use transition, in addition to the local neighborhood effect.In order to make the land suitability analysis feasible both biophysically and socioeconomically, a comprehensive consideration of the social, economic, political, environmental, and other driving factors, as well as a reflection of the different influences of these factors are needed.Therefore, in this model a series of suitability index layers L K pi, jq for each type of land use was built using the method of multi-criteria evaluation (MCE): where F k,f refers to the scoring of criteria f to land type k, and this score is normalized to 0~1 by the maximum value method.When 1 ő m ő α, Fk,f indicates that the factors composed of a group of spatial distance variables and topography.When 1 ő m ő β, F k,f indicates the constraints composed by protected land, topography, and so forth, and the value of this score is 1 or 0. The weight of factor f for land use type k is represented by W k,f , and was acquired by an AHP method [37].
Each cell corresponded to a certain type of land use.The historical trends of the land use changes, land suitability, and related policies and economics were used to construct the transition rules.They controlled the transition possibility of each cell together, and the state of the cell changed when its probability exceeded the control threshold.The cells within the cellular space synchronize update state according to these local transition rules, and the entire cellular space was shown as the change in the discrete time dimension.The definition of the transition rules is the simulate core part of a standard CA, and they can guide a dynamic evolution.The LLUC-CA model's transition rules can be indicated as the land use transition potential.Due to the complexity of nature, some unforeseen characteristics of the process of land use change cannot be explained by independent variables.This can cause the simulation results to be closer to reality through the introduction of random variable to express the uncertainty of the nature and geography of a CA.In each simulation, the final transition probabilities to land use k of cell (i,j) at time t were calculated by the local land use competition intensity, land suitability index, and random variable as follows: where γ is a random variable with a range from 0 to 1, and α is a parameter for the control of the random variable range.The dynamic iterative calculation of the neighborhood change is required in the land use change simulation process, the neighborhood cell state is dynamic, so that the land use transition potential is also dynamic, and at each time step the transition rules of the model are updated.

Demand Analysis of Different Land Types
Aside from the competition between the different land use types and the land suitability, the land use change pattern also depended on the demand for the land use.The Markov model, as a raster-based spatial probability model, is often used to forecast geographical events that have no after-effect characteristics, and has been widely applied to LUCC simulations.Under the assumption that the current social economic mode is under the condition of constant development, the Markov model has the capacity to predict the quantity change of the land use types in the future [38,39].The transition areas file from a Markov Chain analysis of two prior land use maps established the quantity of the expected land use change from each existing category to each other category during the next time period.Assuming that there are m kinds of land use types, from time t to time t+1 the land use state i transform to j(i,j = 1, 2, ..., m), then the transition probability can be described as: where A t i is the area of land use type i at time t, and ∆A ij is the changed area from time t to time t+1 , then the land use state i transform to j.Since there are m kinds of possible states during the land use change process, the transition probability matrix is represented as follows: where 0 ő Pij ő 1 and ř m j"1 P ij " 1.Therefore, at time t + 1, the area of land use type k can be described as:

Multi-Target Land Use Competition Allocation Algorithm
The land use transition potential determined the possible location and direction of the land use state change.The different types of land use demands determined the quantity of the land use changes, and the two worked together to determine the final land use pattern.In most models, a rule-based system is used to allocate the actual land use changes based on the suitability map, including using a threshold to select the locations with the highest suitability, or simulate the land use competition based on the land use type specific characteristics [26].Assuming that there are m kinds of land use types, from Figure 3 and Equation (3) it can be determined that there are m kinds of transition possibilities and corresponding potentials for cell (i,j) at time t.A simple method is to compare the m transition potentials, and the land use types with the highest potential as the land use change direction for cell (i,j).However, this will lead to the explosion phenomenon of the dominant land use type.In this study, a multi-target land use competition allocation algorithm was built, with the land use transition potential as the basis, and the land use demand as a constraint.It was able to meet the demands of the different land use types on the spatial locations and quantities in the process of the simulation.The algorithm was an allocation mechanism based on a variety of land categories as the initial simulation, which could simulate the multiple land use changes.In Figure 3 k is exceed or equal to A t`1 k , go to step four; otherwise change the land use state S t pi, jq of cell (i,j) to k, set the potential p t pi, jq of cell (i,j) as zero, discard the other potential of cell (i,j), then go back to step one; (4) If all kinds of land use types demands are satisfied, or all the cells are changed, then stop the allocation, we will acquire the simulated land use map at time t+1; otherwise, we set all cells' transition potential P t k to land use k as zero, give up the possibility of other cells transform to land use k, then go back to step one.

LLUC-CA Model Structure
Figure 4 is the model structure of the local land use competition cellular automata (LLUC-CA).The CA was combined with the Markov model to control the transition quantity when coupled with the GIS, and the local rules were used to simulate the global and complex patterns of the land use.The model adopted a loose structure, and the GIS, MCE, Markov, and CA programs run independently, exchanging information via data files.The spatial data's input, conversion, and spatial analysis were performed in the GIS software, the data were extracted from the GIS database, and then input into the MCE, Markov, and CA program for processing.The computation results were then transmitted to the GIS for display and further processing.The LLUC-CA model consisted of three stages, as follows: data preprocessing and input, parameter setting and production, and decision-making of the land use type changes and output.In the first stage, ArcGIS, ERDAS, IDRISI, as well as other software, were used to process the land use map, remote sensing images, topography, and planning data, as well as to create a comprehensive geographic information database, where the data structure was consistent, and the spatial reference was unified.In the second stage, the required parameters for the LLUC-CA were acquired, which included: the land use competition intensity index under the influence of local land use competition at the scale of a cellular neighborhood; the land suitability index computed by MCE under the influence of land use macro driving factors; and the area demand for the different land use types predicted by the Markov model.In the last stage, based on the land use transition potential, which was computed by the local land use competition intensity, land suitability and random variables, and different land use area restriction demands, every cell's land use type change was determined with a multi-target land use competition allocation algorithm.The prototype of the LLUC-CA was implemented in Visual Studio 2010 integrated development environment, and the program language was Python.The GDAL and Numpy library functions were also used.
three stages, as follows: data preprocessing and input, parameter setting and production, and decision-making of the land use type changes and output.In the first stage, ArcGIS, ERDAS, IDRISI, as well as other software, were used to process the land use map, remote sensing images, topography, and planning data, as well as to create a comprehensive geographic information database, where the data structure was consistent, and the spatial reference was unified.In the second stage, the required parameters for the LLUC-CA were acquired, which included: the land use competition intensity index under the influence of local land use competition at the scale of a cellular neighborhood; the land suitability index computed by MCE under the influence of land use macro driving factors; and the area demand for the different land use types predicted by the Markov model.In the last stage, based on the land use transition potential, which was computed by the local land use competition intensity, land suitability and random variables, and different land use area restriction demands, every cell's land use type change was determined with a multi-target land use competition allocation algorithm.The prototype of the LLUC-CA was implemented in Visual Studio 2010 integrated development environment, and the program language was Python.The GDAL and Numpy library functions were also used.

Model Parameters Prepare and Implementation
In addition to the basic map of the actual land use, the following parameters also needed to be prepared in order to implement the LLUC-CA model in the study area: the computation of the land suitability maps based on the MCE method, and the computation of the transition matrix using the Markov chain analysis.This study considered the effects of the current land use situation, town centers, main roads, light rail, coastline, and the topographic and spatial regulation of the construction land, based on the biophysical and socioeconomic characteristics, on the land use changes, for the purpose of computing the suitability maps for the agricultural land, construction land, forestland, and water, respectively.Among these, the scores of the first five factors to the four land use types were evaluated by distance, the topographic data were expressed by slope, and the results were standardized (values between 0 and 1).The areas with a slope greater than 25% were not allowed for cultivation and construction, the transportation land was not suitable for cultivation, and the spatial regulation of the construction land did not allow construction.These three criteria were the constraints for the land suitability evaluation.An integrated method combining the domain experts and the AHP method was adopted to obtain the relative significance values of the factors contributing to the suitability of the land for development.A satisfactory Consistency Ratio was obtained for the four land use types (Table 2).Subsequently, according to Equation (3), and the standardized criteria images, the suitability maps for 2003 were calculated using the MCE method (Figure 5).

Model Parameters Prepare and Implementation
In addition to the basic map of the actual land use, the following parameters also needed to be prepared in order to implement the LLUC-CA model in the study area: the computation of the land suitability maps based on the MCE method, and the computation of the transition matrix using the Markov chain analysis.
This study considered the effects of the current land use situation, town centers, main roads, light rail, coastline, and the topographic and spatial regulation of the construction land, based on the biophysical and socioeconomic characteristics, on the land use changes, for the purpose of computing the suitability maps for the agricultural land, construction land, forestland, and water, respectively.Among these, the scores of the first five factors to the four land use types were evaluated by distance, the topographic data were expressed by slope, and the results were standardized (values between 0 and 1).The areas with a slope greater than 25% were not allowed for cultivation and construction, the transportation land was not suitable for cultivation, and the spatial regulation of the construction land did not allow construction.These three criteria were the constraints for the land suitability evaluation.An integrated method combining the domain experts and the AHP method was adopted to obtain the relative significance values of the factors contributing to the suitability of the land for development.A satisfactory Consistency Ratio was obtained for the four land use types (Table 2).Subsequently, according to Equation (3), and the standardized criteria images, the suitability maps for 2003 were calculated using the MCE method (Figure 5).The Markov chain analysis was used to compute the transition probabilities matrix based on the actual land use maps from 1988 to 2003, and then Equation ( 6) was used to compute the different land use types demand for 2012 (Table 3).(3) the demand area for 2012) were input into the LLUC-CA model before running the program.After all the required basic parameters were calibrated (About the details of the model calibration, please see reference [40]) and configured, the LLUC-CA simulation program was begun.During the simulation process, the program first produced the local land use competition intensity index, and then generated the land use transition potential for every cell, which had been integrated with the land suitability index by Equation (3).The simulated map from 2012 was computed by the multi-target land use competition allocation algorithm, according to the land use demands.The Markov chain analysis was used to compute the transition probabilities matrix based on the actual land use maps from 1988 to 2003, and then Equation ( 6) was used to compute the different land use types demand for 2012 (Table 3).(3) the demand area for 2012) were input into the LLUC-CA model before running the program.After all the required basic parameters were calibrated (About the details of the model calibration, please see reference [40]) and configured, the LLUC-CA simulation program was begun.During the simulation process, the program first produced the local land use competition intensity index, and then generated the land use transition potential for every cell, which had been integrated with the land suitability index by Equation (3).The simulated map from 2012 was computed by the multi-target land use competition allocation algorithm, according to the land use demands.

Analysis of Land Use/Cover Changes
The land use cover of Jinshitan developed greatly from 1988 to 2012, from the National Scenic Spot to the National Geopark, and then to the 5A Tourist Area.The general trend was the rapid growth of construction land, and the continuing decline of agricultural land and forestland, with water decreasing slightly, and then leveling off (Figures 6 and 7a-c).It can be seen from these figures that the agricultural land area in 1988 was 21.29 km 2 , and had decreased to 8.97 km 2 in 2003, a 57% decrease in 15 years, and then further decreased to 4.08 km 2 in 2012.Meanwhile, the construction land experienced a rapid growth, from just 3.62 km 2 to 18.16 km 2 , a four-fold increase in 15 years, resulting in a total land area of 27.97 km 2 in 2012.At the same time, the decrease of the forestland was relatively small, decreasing from 26.32 km 2 to 20.28 km 2 over the 25-year period, and the area of water remained stable.The results indicated that the increase of the construction land was at the cost of a transition from a large agricultural land area, as well as some forest area.This transition is reflected in the probability matrix in Table 3 from 2003 to 2012, which indicates that there may be 37.83% of agricultural land converted to construction land, and 12.71% will be converted to forestland; 15.74% of forestland will be converted to construction land; and 8.61% of the water will be converted to construction land.At the same time, 9.43% of the construction land will be converted to forestland, and the mutual conversion between the other land use types is less significant.

Analysis of Land Use/Cover Changes
The land use cover of Jinshitan developed greatly from 1988 to 2012, from the National Scenic Spot to the National Geopark, and then to the 5A Tourist Area.The general trend was the rapid growth of construction land, and the continuing decline of agricultural land and forestland, with water decreasing slightly, and then leveling off (Figures 6 and 7a-c).It can be seen from these figures that the agricultural land area in 1988 was 21.29 km 2 , and had decreased to 8.97 km 2 in 2003, a 57% decrease in 15 years, and then further decreased to 4.08 km 2 in 2012.Meanwhile, the construction land experienced a rapid growth, from just 3.62 km 2 to 18.16 km 2 , a four-fold increase in 15 years, resulting in a total land area of 27.97 km 2 in 2012.At the same time, the decrease of the forestland was relatively small, decreasing from 26.32 km 2 to 20.28 km 2 over the 25-year period, and the area of water remained stable.The results indicated that the increase of the construction land was at the cost of a transition from a large agricultural land area, as well as some forest area.This transition is reflected in the probability matrix in Table 3 from 2003 to 2012, which indicates that there may be 37.83% of agricultural land converted to construction land, and 12.71% will be converted to forestland; 15.74% of forestland will be converted to construction land; and 8.61% of the water will be converted to construction land.At the same time, 9.43% of the construction land will be converted to forestland, and the mutual conversion between the other land use types is less significant.

Model Validation
For model validation, the actual land use map for 2012 and the LLUC-CA simulated land use map were compared (Figure 7c,d).The overall simulation land use map was very similar to the actual one by visual comparison, especially in the middle and southern sections of the study area.The simulation results of the forestland and water were very close to the actual land use map both in spatial layout and quantity, and the actual water area was 5.63 km 2 , while the simulation area reached 5.32 km 2 .These results indicated that the model effectively reproduced the area for each land use type in the land use change simulation.The main deviation was located in the northern areas, where

Model Validation
For model validation, the actual land use map for 2012 and the LLUC-CA simulated land use map were compared (Figure 7c,d).The overall simulation land use map was very similar to the actual one by visual comparison, especially in the middle and southern sections of the study area.The simulation results of the forestland and water were very close to the actual land use map both in spatial layout and quantity, and the actual water area was 5.63 km 2 , while the simulation area reached 5.32 km 2 .These results indicated that the model effectively reproduced the area for each land use type in the land use change simulation.The main deviation was located in the northern areas, where the construction land of the northwest and northeast areas were scattered around the agricultural land, and the actual agricultural land area was 4.08 km 2 , while the simulation area reached 5.56 km 2 .Fan et al. ( 2013) conducted related research in Jinshitan from the perspective of a tourist land use landscape pattern, which also had a high consistency with this study from the perspective of land use patterns [41].For the purpose of further verifying the improvement effects of the LLUC-CA model, a comparative study was carried out using a CA-Markov model, and the suitability maps used were the same.The simulation results of the CA-Markov are shown in Figure 7e, and are very similar to those shown in Figure 7d.However, the construction and agricultural land simulation was found to be poorer in the northeast and southwest, where the growth area was too large and the simulation position was inaccurate.For the quantitative validation of the model's accuracy, the actual land use map and the simulated land use map were verified based on a Kappa coefficient.The LLUC-CA's Kappa coefficients were 64.46%, 77.21%, 85.30%, and 99.14% for the agricultural land, construction land and forestland, and water, respectively.The overall simulation success was 88.74%, which meant that for most of the spatial locations of the land use types, the LLUC-CA model performed correct simulations, especially for the simulation of the forestland and water.On the other hand, the CA-Markov's Kappa coefficient was 61.54%, 74.62%, 85.12%, and 99.13% for the agricultural land, construction land and forestlands, and water, respectively, and the overall simulation success was 86.82%, which demonstrated that the spatial accuracy of the simulated land use map using the LLUC-CA model was higher than that of the CA-Markov model.In conclusion, the LLUC-CA model can potentially be used to simulate land use cover change on a large scale.

Conclusions
The LUCC is a giant complex system, and CA is an important tool that is suitable for complex system simulation.On the one hand, in recent years, research has been mainly concentrated on using artificial intelligence methods to obtain transition rules in order to improve the simulation quality.However, artificial intelligence methods belong to the black box model, and it is not easy to determine For the purpose of further verifying the improvement effects of the LLUC-CA model, a comparative study was carried out using a CA-Markov model, and the suitability maps used were the same.The simulation results of the CA-Markov are shown in Figure 7e, and are very similar to those shown in Figure 7d.However, the construction and agricultural land simulation was found to be poorer in the northeast and southwest, where the growth area was too large and the simulation position was inaccurate.For the quantitative validation of the model's accuracy, the actual land use map and the simulated land use map were verified based on a Kappa coefficient.The LLUC-CA's Kappa coefficients were 64.46%, 77.21%, 85.30%, and 99.14% for the agricultural land, construction land and forestland, and water, respectively.The overall simulation success was 88.74%, which meant that for most of the spatial locations of the land use types, the LLUC-CA model performed correct simulations, especially for the simulation of the forestland and water.On the other hand, the CA-Markov's Kappa coefficient was 61.54%, 74.62%, 85.12%, and 99.13% for the agricultural land, construction land and forestlands, and water, respectively, and the overall simulation success was 86.82%, which demonstrated that the spatial accuracy of the simulated land use map using the LLUC-CA model was higher than that of the CA-Markov model.In conclusion, the LLUC-CA model can potentially be used to simulate land use cover change on a large scale.

Conclusions
The LUCC is a giant complex system, and CA is an important tool that is suitable for complex system simulation.On the one hand, in recent years, research has been mainly concentrated on using artificial intelligence methods to obtain transition rules in order to improve the simulation quality.However, artificial intelligence methods belong to the black box model, and it is not easy to determine the rule hidden under the spatial pattern change.On the other hand, in all types of the extended CA models, without using artificial intelligence algorithms, most can only simulate one target land use type, without the ability to multi-target land use simulation at the same time, and the discovery comprehensive analysis of the relationship between the land use types is weak.The existing regional land use models, as well as the related research, are mainly used for the scope of the study area [14], while in this study, the object of study was not only aimed at a small scale area, but also emphasized the research method of the LLUC-CA cell transition rules within the neighborhood.
In this study, a new CA model of local land use competition was introduced, which integrated the advantages of the models, on the basis of summarizing previous research for land use change simulation.In order to explore the interaction mechanism between the various land use types, and to meet the needs of more land types comprehensive simulation, a cellular neighborhood was taken as a local unit.The interaction mechanism within a cellular neighborhood was analyzed, and research of the local land use competition analysis method influenced by the conditions of the cell itself, and land use type distribution within the neighborhood were studied.This study focused on researching this new method of using local land use competition to develop the cell transition rules, explored a multi land use type CA model that had clear physical meaning, and established a complete experiment method.In addition, the land use transition potential generated by combining the local land use competition intensity index, which was computed by local land use competition analysis, the land suitability index computed by the MCE method, and random variables and different area demands that were predicted by the Markov model, were integrated by a multi-target land use competition allocation algorithm, followed by implementing the multi-target land use change simulation.The results indicated that the local land use competition could provide a good expression of the spatial autocorrelation, and the Markov chain analysis helped to control the quantity of land use change.The MCE generated suitability maps in order to help control the spatial distribution of the land use change, and the multi-target land use competition allocation algorithm determined the last change location and quantity.The prototype of the model was programmed by Python and GDAL, which has been successfully applied to the land use change simulation of Jinshitan, using land use maps from 1988, 2003, and 2012.When compared with the LLUC-CA simulated map, the CA-Markov simulated map and the actual land use map, the results indicated that the model had a high simulation accuracy.In the simulation, at the same time, the Kappa coefficient was 64.46%, 77.21%, 85.30%, and 99.14% for the agricultural, construction and forestland, and the water, respectively.The proposed model showed a high accuracy of quantity and spatial distribution in the simulation of the land use changes.The model had a strong spatial self-organization ability for the simulation of the land use changes, and had a certain advantage which made it an effective method.
Land use change is a very complicated geographical process.It is affected by natural conditions and cultural influences, as well as political, planning, technological, and many other socioeconomic factors.Consequently, an accurate simulation and prediction of land use change is always difficult.The LLUC-CA model proposed in this study can achieve the simulation of multi land use type comprehensive evolution, and monitor the status and trends of the land use changes.It can provide a method for the land use change driven forces analysis, the analysis and prediction of urban land use expansion, and the development and protection of agricultural land, forestland, and water.It can also provide important data to assist the planners, government, and public in future decision-making.Although the model considers cell interaction within the local neighborhood and the influence of macro natural driving and planning factors, the continuity and variability of the transition rules in different time intervals, as well as the influence of policy, economic, and demographic factors on land use change simulation, still require further research.Additionally, in the follow-up study, the demand analysis of the land use area should consider adding planning of the land use areas, in addition to using the previous years' data of the land use changes.The land suitability evaluation can add a great deal of economic and demographic data, and an artificial intelligence method could be introduced to generate the local land use competition intensity indexes.
).It faces the Yellow Sea on the east, and is composed of the eastern peninsula, western peninsula, open heartland, and the beach between the two peninsulas.It covers a land area of approximately 58 km 2 .Located adjacent to the Dalian urban area, it is a typical area with a rural-urban zone and a coastal region.Beginning as a rural town dominated by traditional agriculture, it was identified as a National Scenic Spot in 1988, National Geopark of China in 2004, and a 5A Tourist Area in 2010.Its construction and development into a tourist resort in China had a strong representation.For nearly 25 years, the study area has experienced rapid economic development and land use change, and therefore possesses high research value.ISPRS Int.J. Geo-Inf.2016, 5, 106 3 of 15

Figure 1 .
Figure 1.Location of study area, Dalian, Liaoning, China.The remote sensing image was generated based on 10 m multispectral and 2.5 m panchromatic spatial resolution SPOT5 image.

Figure 1 .
Figure 1.Location of study area, Dalian, Liaoning, China.The remote sensing image was generated based on 10 m multispectral and 2.5 m panchromatic spatial resolution SPOT5 image.
ISPRS Int.J. Geo-Inf.2016, 5, 106 5 of 15 competing relationship between the land use types is illustrated with Figure2.As shown in the figure, Num (A) = 7 expressed that the number of cells of the agricultural land in the neighborhood scenario is 7, Num (C) = 8 expressed that the number of cells of construction land in the neighborhood scenario is 8, and so on.Due to the fact that the land type transformation of the central cell was based on its independent neighborhood scenario, the interactions between the land use types were not equal.

Figure 2 .
Figure 2. Sketch of mutual competition between different land use type cells in the neighborhood.

Figure 2 .
Figure 2. Sketch of mutual competition between different land use type cells in the neighborhood.
, the flow chart of the multi-target land use competition allocation algorithm can be seen, which was divided into four steps: (1) Search the cell which has the maximum value of transition potential as Maxpp t k pi, jqq of the land use map, keep a record of the position (i,j) and the land use change direction k; (2) Compute the area of land use type k as A t k ; (3) Compare A t k and the area of land use demand A t`1 k predicted by Markov model, if A t

Figure 3 .
Figure 3. Flow chart of multi-target land use competition allocation algorithm.Figure 3. Flow chart of multi-target land use competition allocation algorithm.

Figure 3 .
Figure 3. Flow chart of multi-target land use competition allocation algorithm.Figure 3. Flow chart of multi-target land use competition allocation algorithm.

Figure 6 .
Figure 6.Area of land use classes in the study area.

Figure 6 .
Figure 6.Area of land use classes in the study area.
ISPRS Int.J. Geo-Inf.2016, 5, 106 12 of 15 Fan et al. (2013) conducted related research in Jinshitan from the perspective of a tourist land use landscape pattern, which also had a high consistency with this study from the perspective of land use patterns [41].

Figure 7 .
Figure 7. Land use maps: (a) actual land use map in 1988; (b) actual land use map in 2003; (c) actual land use map in 2012; (d) LLUC-CA simulated land use map in 2012; and (e) CA-Markov simulated land use map in 2012.

Figure 7 .
Figure 7. Land use maps: (a) actual land use map in 1988; (b) actual land use map in 2003; (c) actual land use map in 2012; (d) LLUC-CA simulated land use map in 2012; and (e) CA-Markov simulated land use map in 2012.

Table 1 .
Classification of land use types.

Table 2 .
Weights of influencing factors of land suitability.

Table 3 .
Land use transition probability matrix and area demand.

Table 2 .
Weights of influencing factors of land suitability.

Table 3 .
Land use transition probability matrix and area demand.