GISdevelopment.net ---> AARS ---> ACRS 2000 ---> Landuse



A Land Use Change Study Using Cellular Automata

Jinn-Guey Lay
Associate Professor, Department of Geography
National Taiwan University, PO Box, 23-76, Taipei, Taiwan
Tel: (886)-2-23621499 Fax: (886)-2-23622911
Email: jglay@ccms.ntu.edu.tw

Key words Spatial modeling, Land Use, Cellular Automata (CA)

Abstract
Land use change is a major issue of global environment change. The modeling and projecting of land use change is essential to the assessment of consequent environmental impacts. Recent development of cellular automata (CA) provides a powerful tool for the dynamic modeling of land use change. This research adopts the spatial evolution concept embedded in CA and applies it to land use change study in Tansui Watershed. Digital land use data of two separate years were complied and analyzed using GIS software. A computer program was coded to analyze the neighborhood condition of each target cell. Summary of neighborhood conditions may reveal the dynamic process of land use change and thus enhance our understanding on transition rules, the heart of a CA.

Introduction
Land use activity is a linkage between human and environment. It is like a seesaw carrying human activities and environment on two opposite ends. Human activity of development is a primary driving force for global environmental changes while the changes of environment in return affects land use types. Envisioning the consequent effects of land use changes, IGBP (International Geosphere and Biosphere Programme) and IHDP (International Human Dimension Programme) co-organized a working group to set up research agenda and promote research activity for land use and land cover changes (LUCC). The working group suggested three core subjects for LUCC research, such as: situation assessment, modeling and projecting, and conceptual scaling. The ultimate goal of global change study is to assess the impacts under each possible scenario and suggest preventive actions. The modeling and projecting of land use change is essential for scenario analysis and the assessment of LUCC. Consequently, issues related to data, information, and modeling have attracted many research interests ranging from local authorities to global organizations. The importance of land use change has long been recognized in Taiwan with many research projects undergoing. This paper presents partial result of an undergoing research, supported by National Science Council, for the study of land use change in Taiwan.

Cellular Automata and Land Use Change
Recent development of GIS technology enhances the analytical power needed for the study of land use and land cover change. Existing land use data are being digitized while much more data are readily created in digital format, i.e., remotely sensed imagery. The development and applications of digital land use database have been successful as documented in many reports (Lay, 2000). With the help of GIS, these data can be easily overlaid and provide information needed for situational assessment. On the other hand, the methodology for predicting and modeling land use change is relatively immature. Land use change is a dynamic spatial process involves complex interactions between many factors at various spatial extents. The complexity of this dynamic process makes the creation of a comprehensive model very challenging. Recent development of cellular automata (CA) brings in a new perspective for land use change. There are increasing interests in using CA for land use study (Batty, 1998; Clarke et al., 1997; White and Engelen, 1997) The concept of cellular Automata is originated from research for a self-replicating machine, a robot type machine equipped with visions and encoded Turing machine that can assemble a copy of itself from component parts (Firebaugh, 1988). Although such machine has not been successfully constructed, an abstract model of this self-replicating system was constructed as two-dimensional cellular automata. Cellular automata may be represented by a set of simple production rules while its outcome may mimic a very complex system (Firebaugh, 1988, p320-323). From the application side, cellular automata are dynamic model that inherently integrates spatial and temporal dimension. CA are composed of four elements as described below (White and Engelen, 2000).
  • Cell space: The cell space is composed of individual cell. Theoretically, these cells may be in any geometric shape. Yet, most CA adopt regular grids to represent such space, which make CA very similar to a raster GIS.
  • Cell states: The states of each cell may represent any spatial variable, e.g., the various types of land use.
  • Time steps: A CA will evolve at a sequence of discrete time steps. At each step, the cells will be updated simultaneously based on transition rules.
  • Transition rules: These rules are the heart of a CA that guide its dynamic evolution. A transition rule normally specifies the states of cell before and after updating based on its neighborhood conditions.
This research adopts the concept of CA by analyzing the neighborhood conditions for each target cell of change. A computer program is coded to analyze the states of neighborhood of each center cell. This study used the data of Tansui Watershed as sample data for testing and analysis. The watershed is located in northern Taiwan and the heartland of Taiwan both economically and politically. The target area of research is the plain region of this watershed, selected with elevation below 100 meters, which is about 32400 hectares of acreage. The change of land use in this area was analyzed based on the spatial condition of neighborhood. Summary of neighborhoods states is fundamental to the development of transition rules. The purpose of this analysis is to identify the relationship between land use change and surrounding environment. Findings from this analysis may reveal the dynamic process of land use change and thus enhance our understanding on transition rules, the heart of a CA.

Land Use Change in Tansui Watershed
With a population close to 6 millions, Tansui watershed has experienced a dramatic change during the past decades with huge amount of land use and land cover change. Land use data of the research area are available in various formats, such as: historical documents, paper maps, remotely sensed imagery, and digital data. For convenience reason, digital data are chosen as source data for this research. These data are available separately in 84 data files for each year of 1971 and 1977 in Arc/Info vector formats. The data files of each year are merged together into a big one to cover the whole research area. Since vector files tempt to result in erroneous sliver polygons thus false changes, this research converted the vector files into raster format before further process. Original data is composed of 99 categories of land use types. For the purpose of environmental change, they are reclassified into seven categories: settlement, farmland, bare land, forest, river and reservoirs, other water body, and roads. The data of two years are overlaid with each other to identify the land use change.

The overlay function is conducted using Arc View's Tabulate Area function. The change is shown in Table 1 and Table 2. As shown in Table 1, the increase of settlement is most significant with a figure near 4% or 1250 hectares. Each figure in Table 2 represents amount of change from 1971 to 1977 for each category. Based on these figures, this research therefore concentrated on investigating the change from farmland to settlement. A computer program is coded for investigating the neighborhood states of the targets cells.

Table 1: Summay of Land Use Change from 1971 to 1977
Settlement Farmland Forest River and reservoir Other water body Roads Bare land
1971 acreage (hectare) 16328.97 7842.66 782.65 4194.23 827.37 590.34 1904.20
1977 acreage (hectare) 17583.53 6554.05 917.52 3532.54 733.33 800.20 2349.24

1971 percentage (%) 50.29 24.15 2.41 12.92 2.55 1.82 5.86
1977 percentage (%) 54.15 20.18 2.83 10.88 2.26 2.46 7.24


Table 2?Summary of Land Use Change
1977 1971 Settlement Farmland Forest River and reservoir Other water body Roads Bare land
Settlement 46.92 1.36 0.29 0.12 0.13 0.42 1.05
Farmland 3.58 17.13 0.49 0.34 0.64 0.66 1.33
Forest 0.30 0.22 1.48 0.05 0.01 0.04 0.31
River and reservoir 0.83 0.20 0.15 9.75 0.28 0.11 1.59
Other water body 0.49 0.37 0.05 0.17 0.98 0.07 0.42
Roads 0.47 0.14 0.02 0.04 0.06 1.04 0.05
Bare land 1.56 0.76 0.34 0.42 0.16 0.13 2.49


Neighborhood
Analysis The states of neighborhood is summarized based on the pattern of its 8 neighboring cell regardless its orientation or sequences. Taking the example shown in Figure 1, the land use of a central cell changed from farmland to settlement, and its neighborhood are composed of one settlement cell, six farmland cells, and one forest cell. This neighborhood condition is coded as: 1610000, with each digit indicating the number of cells under each land use sequentially. T1 T2 Farmland Farmland Farmland Settlement Farmland Forest =

T1                                                   T2
Farmland Farmland Farmland

Settlement Farmland Forest =>

Settlement
Farmland Farmland Farmland

Figure1 Illustration of Neighborhood Condition The
Computer program consists of five types of query as below.
  1. One to many query: showing the neighborhood statistics of target cells from one land use type to any type;
  2. Many to one query: showing the neighborhood statistics of target cells from any type of land use to a particular one;
  3. One to one query: showing the neighborhood statistics of target cells from one land use type to another;
  4. One to many neighborhood analysis?given a neighborhood condition, showing the percentage of land use change from one type to any other;
  5. Many to one neighborhood analysis?given a neighborhood state, showing the percentage of land use change from any type to a particular one.
Using function 3 of this program, the change from farmland to settlement was investigated first. As shown in Table 3, 58% of the farmland surrounded by farmland remains farmland while 8.9% of them changed to settlement.

Table3?Settlement Change?Percentage > 1??10M grids?
1971 1977 Settlement Farm Land Forest River and reservoir Other Water body Roads Bare land No. of cells Percen- tages (%)
Farmland Farmland 0 8 0 0 0 0 0 451584 58.0
Farmland Settlement 0 8 0 0 0 0 0 69116 8.9
Farmland Bare land 0 8 0 0 0 0 0 26662 3.4
Farmland Roads 0 8 0 0 0 0 0 14687 1.9
Farmland Farmland 1 7 0 0 0 0 0 13478 1.7
Farmland Other water 0 8 0 0 0 0 0 12925 1.7
Farmland Farmland 2 6 0 0 0 0 0 10917 1.4
Farmland Farmland 3 5 0 0 0 0 0 9482 1.2
Since an initial type of land use may result in several different land use changes, this study analyzes the various neighborhood of a particular land use change pattern type. As shown in Table 4, a farmland cell in 1971 may change to settlement in 1977 while the neighborhoods of these farmland cells may be very different. For example, sixty percent of the changed cells are surrounded by 8 farmland cells followed by the category of cells with 3 settlement cells and 5 farmland cells. These finding suggests that a land use change pattern may occur in various situations.

Table 4: Neighborhood Conditions of 1971 Farmland to 1977 Settlement (10M grids)
Neighborhood conditions Number of cells Percen- tages (%)
Settlement Farm Land Forest River and reservoir Other Water body Roads Bare land
0 8 0 0 0 0 0 69158 60
3 5 0 0 0 0 0 9120 8
1 7 0 0 0 0 0 6014 5.3
4 4 0 0 0 0 0 5785 5.1
2 6 0 0 0 0 0 5374 4.7
0 5 0 0 0 3 0 2154 1.9
0 7 0 0 0 1 0 1862 1.6
5 3 0 0 0 0 0 1491 1.3
0 6 0 0 0 2 0 1357 1.2


Conclusion
This research adopts concepts of CA to investigate land use change from a spatial evolution perspective. Factor of concern is limited to neighborhood conditions in this research. Although this approach is rather simple, yet the underneath merit lies in its ability in revealing the spatial dynamic of land use evolution. From a quantitative perspective, the prediction of future land use change may include factor of: amount of changes, rate of changes, type of changes, and location of changes. CA alone may be too simplified to predict all characteristics of changes, yet there is great possibility for CA to integrate with other methods. The creation of transition rules is but a fundamental step, yet the most challenging one, in building a comprehensive model of land use change. Findings from this research help reveal the complex process of land use change, as shown in the various changes resulted from a same pattern of neighborhood. Such finding indicates a deterministic approach of CA may not fit with real world situation. The calibration of such deficiency lends itself to many possible topics for further research.

Acknowledgement
This research is funded by National Science Council, Taiwan, Republic of China (NSC89-2750-P-002-009). The test data are provided by the Council of Agriculture, Taiwan, Republic of China.

Reference
  • Batty, M. 1998, Urban evolution on the desktop: simulation with the use of extended cellular automata, Environment and Planning A, 30:1943-1967.
  • Clarke, K.C., S. Hoppen, and L. Gaydos ,1997:A self-modifying cellular automata model of historical urbanization in the San Francisco Bay area, Environment and Planning B: Planning and Design, 24:247-261.
  • Firebaugh, M.W. 1988, Artificial Intelligence A Knowledge-Based Approach, Boyd and Fraser Publishing Company, Boston, 740p.
  • Lay, J.G., 2000, A primary study on the integration of cellular automata and GIS, Bulletin of the Chinese Geographical Society, 28:109-126.
  • White R., and G.Engelen,1997:Cellular automata as the basis of integrated dynamic regional modelling, Environment and Planning B, 24:235-246.
  • White, R., and G. Engelen, 2000, High-resolution integrated modeling of the spatial dynamics of urban and regional systems, Computer, Environment and Urban System 24:383-400.