02 Whole
02 Whole
02 Whole
Asep Karsidi
May 2004
Abstract
This study is concerned with land use/land cover change detection, identification,
analysis and prediction using remote sensing and GIS techniques in the downstream of
'West
the Citarum watershed and its surroundings in Java, Indonesia. Supervised
Maximum Likelihood classification of PCA and NDVI transformed images are used to
approach was used to detect land use/land cover changes, and a Markov Cellular
automata model is then used to predict possible future land use/land cover patterns in
the study area."Leafon" and "leaf off'phenomena of the broad leaf vegetation cover
have been recognised related to dry and wet season as well as rice field (planted) and
rice field (unplanted) related to growing season in the study area. Forest and plantation
area were extensive in wet season and less in dry season. Rice field þlanted) area was
large in harvesting time and less in planting time. Settlement has increased
continuously and is not influenced by season or weather. Overall' the KIA of the
classification was 0.89. Settlement and rice field are the main land use/land cover types
that have been changed and this is related to factors such as proximity to roads and to
urban and semi-urban centres. There is an indication that land use/land cover in the
study area was converted from intensive agriculture land such as rice field to
settlement, rather than from less intensive uses such as open/dry land, plantation or
forest. Discriminant analysis as well as overþ and simple linear analysis support
factors such as proximity to roads, urban and semi-trban centres, as well as slope, as
being most influential in land use/land cover change in the study area. The Markov
Cellular automata model affords a powerful descriptive and predictive model for land
11
use/land cover change and for future land use/land cover distribution in the study area,
but it needs some adjustment in order to obtain suitable results. Markov transition, as
well as suitabilit¡ maps of each land use/land cover category are created'
ll1
Abbreviations and Glossary
BAKOSURTA}IAL The Coordinating Agency for National Survey and
Mapping
BAPPEDA Regional Development Planning Boards
BAPPENAS The National Development Planning Agency
BPN The National Land Agency
BPPT The Agency for Assessment and Application of
Technology
CA Cellular Automata
'Water
DGWRD Directorate General of Resources Development
lV
Gontents
ii
Abstract
iv
Abbreviations and Glo ssarY
vlll
List of Figures xi
List of Tables
xiv
Declaration xv
Acknowledgement
CIIAPTER ONE I
l.l Introduction
7
1.2 Aims and Objectives 9
1.3 The Context ofthe StudY
18
1.4 Research Framework
2t
1.5 Thesis Outline
CHAPTER TWO 22
2.1 Introduction 22
2.2 Defining Land use/land cover change 25
2.3 The dynamic of Land use/land cover change
change 30
2.4 Some aspects related to Land use/land cover
pressur 30
2. .LLanduse and Demo graphic 32
Z.4.ZLanduse conflict and Regional food security
35
2.4.3 Landuse, Industry and Manufactruing 37
2.4.4 Landuse and Agriculture 39
2.4.5 SpatiaÍdimension of Land use/land cover change
detection
z.+.0 rrre important of Land use/land cover change
42
in sPatial Planning
43
2.4.TLanduseevolutionanddemographicpressrrreinlndonesia
planning in
2.+.t yalduse/land cover change dãtection and spatial 47
Indonesia
50
2.5 Conclusion
v
IIAPTER THREE
3.1 Introduction 52
3.2 Selection of Study area 53
3.3 Socio-EconomicCharacteristics 55
3.3.1 PoPulation and demograPhY 56
3.3.2 Employment 57
3.3.3 The Land Tenure 58
3.3.4 Existing Land use 6l
3.3.5 Infrastructure 63
3.4 BiophysicalCharacteristics 65
3.4.1 The Physiography and Soil type 66
3.4.3 Climate, Irrigation and Growing season 68
3.5 The availability of spatial data 7l
3.5.1 Existing sPatial data 7l
3.5.2 Satellite imagery 74
3.6 Conclusion 75
CTIAPTER F'OUR
4.1 Introduction 77
4.2 Land use/land cover change detection techniques 78
4.2.1Digital change detection of land use/land cover 83
4.2.1.1 Pre-classification 88
4.2.1.2 Po st-Classification Comparison 96
4.2.1.3 Hybrid 98
4.2.2 Selection of the appropriate detection and identification
Techniques 99
4.3 Method 101
4.3.1 Digital change detection techniques 103
4.3.1.1 Image Pre-Processing 103
4.3.1.2 Image Enhancement 104
4.3 .2 Land use/land cover Classification 106
4.3.3 Landuse/land cover change detection r07
4.3.4 Accuracy Assessment 108
4.4 Results digital change detection and identification in the study area 108
4. 4 . I lrlrttal reco gnitio n
108
4.4.2 Image Enhancement tt2
4.4.3 Larà use/land cover Classification 11s
4.4.4Theland use/land cover change detection t2l
4.4.4.1Annual Land use/land cover change from 1989 to 1993 122
4.4.4.2land use/land cover change in 1993 to 1995 and
1995 to 1997 125
4.4.4.3 Land use/land cover change 8- year time interval
1989-1997 t26
4.4.5 Landuse/land cover change identification and re-classification t28
4.5 Discussion 131
4.6 Conclusion t34
V1
CHAPTER FIVE
5.1 Introduction 137
CIIAPTER SIX
6.1 Introduction t7l
6.2 Land use/land cover change prediction mode 172
174
6.2.1 Markov chain model
177
6.2.2 Cellular automata
6.2.3 Prediction based on Markov-Cellular automata
181
Model within GIS
6.3 Method 183
184
6.3 . 1 Markov transition probability calculation
6.3.2 Multi criteria evaluation approach to create land suitability 185
6.3.3 Prediction process based on Markov-Cellular automata model 186
CHAPTER SEVEN
CHAPTER ONE 7
1.1 The Study area of the Citarum Watershed
11
1.2 Province ofthe Republic of Indonesia
Indonesian Country t2
1.3 Percentage of Forest cover in the main Island of
Bekasi and
l.A Housingãnd Industrial estate development area in Bogor,
I7
Karawang
20
1.5 Research framework
CHAPTER TWO 30
2.1 The mutual interaction of land use change dynamic
37
2,.2 Housing and Industrial estate development area in JABOTABEK
38
2.3 Land use in West Java1993
2.4 Overall model spatial planning structure related to
43
Land use
44
2.5 The conc Indonesia
45
2.6 Land use ounding West Java
49
2.7 Integrated Planning System for Regional Development
CHAPTER THREE 52
3.1 Study area as part of the Citarum'Watershed
54
3.2 The irve Kabupaten and Jakarta within study area
in Jakarta and
3.3 The new rown and Industrial estate development
54
Surrounding area
55
3.4 iand rrse tõq¡ in five Kabupaten within study area and Jakarta 62
3.5 Main Land use in West Java: 1989-2000
63
3.6 Roads Networks in the studY area
66
3.7 The Relief of West Java
67
3.8 The Physiography in five Kabupaten of study area
68
3.9 Soiltype itfive Kabupatenwithin study area
69
3.10 The Rainfall Regimes in West Java
70
3.11 The monthly rainfall regimes in the study area
7l
3.12 The Citarum river and inigation in the study area
CHAPTER FOUR
analysis of
4.1 The integrated model of detection, identification and
102
Land use/land cover change
4.2 post-classification approach of land use/land cover
"o-p*irot, r07
Change detection and identification
109
4.3a TM Composite 432 t989
110
4.3b TM Composite 4321990
110
4.3c TM Composíte 4321995
111
4.3d TM Composire 4321997
v111
4.4a The Spectral features of the Principal Component 5 (PC5) tt4
The Spectral features of Band Ratio (NDVD 115
4.4b
The Final classified image of land use/land cover 1989 116
4.5a
4.sb The Final classified image of land use/land cover 1990 l17
4.5c The Final classified image of land use/land cover 1991 t17
The Final classified image of land use/land cover 1992 118
4.5d
The Final classified image of land use/land cover 1993 118
4.5e
4.sf The Final classified image of land use/land cover 1995 t19
119
4.5g The final classified image of land use/land cover 1997
4.6 The Annual land use/land cover change in the study area l2l
4.7 Land use/land cover in the study area 1989 and 1997 r29
130
4.8 Land use/land cover change in the study area
CHAPTER FIVE
5.1a The Classification of population density in the study area r52
153
5.lb The Settlement distribution inthe study area
5.2 The increasing of Settlement and decreasing Rice field within
urban and semi-wban area 155
158
5.3a
159
5.3b
159
5.3c
s.3d
160
161
5.3e
161
s.3f
5.3g
t62
CHAPTER SIX
6.1 The Neighbourhood of the (ij) cell is formed by the (Ij) cell itself
And theiight adjacent cells 182
6.2a Suitability-map for settlement without proximity to Jakarta city 191
6.2b Suitabitity -up for settlement with proximity to Jakarta city 192
6.3a prediction t; óf g-iteration,3-kernel filter, suitability map of settlement
Without proximþ to Jakarta cþ 194
6.3b Original^Classifred Inage 1997 195
x
List of Tables
CHAPTER ONE
I .1 Distribution of Indonesian Population and Growth Rate I 980-2000 1 0
1.2 Forest cover and Population Density in the main Island in Indonesia
Year 1996 13
1.3 The Change of Agriculture area 1981 to 1999 (change fromrice fie1d
to other land use) in Java andBali 14
1.4 Area of rice field and settlement within Kabupaten in the study area 15
1.5 Decreasing rice field and increasing settlement in V/est Java 1989-2000 16
1.6 Number ofNational lndustrial and trading companies 16
1.7 Number and area of Industrial Estate in Kabupaten Bekasi year 2000 18
CHAPTER TWO
2.1 Driving force behind land use change 27
2.2 Change in Agricultural Land between 1981-1999 in Indonesian
Provinces 28
2.3 Changes in Agricultural land in Indonesia7983-1993 55
2.4 The average rice productivity in Java and offJava (Ton/halyr) 34
2.5 Conversion from rice field to settlement and industrial area in the
main island of Indonesia 1981-1999 34
2.6 Number of large and medium manufacturing 1997-2000 35
2.7 The land use classification system in Indonesia 40
2.8 Land use/land cover classification level and their source 4t
2.9 Scale map at different level of Spatial planning application 48
CHAPTER THREE
3.1 Population and Demographic Indicator 56
3.2 Labour Force and Occupation in five kabupaten and Kodya Bogor 58
3.3 Changes in Agricultural Land and Households 1983-1993 59
3.4 The change of Agricultural land to other land use year 1981-1999 60
3.5 Existing Land use in the study area (ha) 62
3.6 Infrastructures in five Kabupaten within the study area 64
3.7 The Growing season schedule 70
3.8 Partial list of Land Resources Inventory and Mapping Agencies in
Indonesia 73
3.9 Landsat TM data 74
CHAPTER FOUR
4.I The Comparison of the Terrestrial, Air photos and Remote-
Sensing techniques 80
4.2 The various techniques to identiS and monitor land use/land cover
change 82
4.3 Comparison of digital detection techniques categorised by
Independentl¡ Simulation and Hybrid 86
x1
4.4 The type of Image enhancement and function 105
4.5 1989
P principal component analysis transformation TM 113
4.6 Land use/land cover classificationscheme 116
4.7 assessment
Summery of Kappa index of agreement (KIA) 120
4.8a The Summary major change in land use/land cover category 198911990122
4.8b The Summary major change in land use/land cover category l990ll99l 123
4.8c The Summary major change in land use/land cover category 199211993 124
4.8d The Summary major change in land use/land cover category l99lll992l24
4.9a The Summary major change in land use/land cover category 199311995125
4.9b The Summary major change in land use/land cover category 199511997 126
4.10 The Summary major change in land use/land cover category 198911997 I27
4.11 Land use/land cover in the study areayear 1989-1997 (ha,) 130
CHAPTER FIVE
5.1 The Method to create area of possible influence or driver factors 144
5.2 Variable used in the discriminant analysis r46
5.3 Summery of land use/land cover change versus slope (1989-1997) 148
5.4 Summery of land use/land cover change versus Physiography
(te9e-tee7) t49
5.5 Summery land use/land cover change versus Land system
(1989-tee7) 150
5.6 Land use/land cover and Population density Classes 152
5.7 Land use/land cover change within 10 km of urban area (ha) 154
5.8 Land use/land cover change within 10 km of semi-urban area (ha) r54
5.9 Percentage of decreasing and increasing rice field and settlement
within urban and semi-urban in the study area r54
5.10a Land use/land cover in 1 km buffer of To11 rcad/way 157
5.10b Land use/land cover in 1 km buffer of State road 157
5.10c Land use/land cover in I km bufÊer of other road 158
5.11 R2 of Linear regression between population and land use/land cover 160
5.12 CovarianceMatrices t63
5.13 Eigenvalues 163
5 .l 4a Standardized Canonical Discriminant Function Coefficients 164
5.14b Structure Matrix 164
5.15 Classifrcation function Coeffi.cients 16s
CHAPTER SD(
6.1 Land suitability factors of each category 185
6.2 The Markov transition of Land use/land cover 1989-1997 189
6.3 The weight value of suitability factors 190
6.4a The Kappa Index of Agreement (KIA) value between
original 1997 -1997 predicted I and original 1997 -1997 predicted 2 t96
6.4b The Kappa Index of Agreement (KIA) value between
original 1989 - original 1997 and original 1989- predictedl997 196
6.5 Summary of probability expected to stay as same category l0- year
Scenario 200
6.6a The Markov transition of land use/land cover for lO-years prediction
(1997-2007), based on 1989-1997 images 20r
xl1
6.6b The Markov transition of land use/land cover for lO-years
(scenariol 2007) assumed that change will low 202
6.6c The Markov transition of land use/land cover for l0-years
(scenario2 2007) assumed that change will very low 20r
6.7 The lO-years change characteristic between Simulation and Scenario 206
6.8a Markov transition of land use/land cover for 2O-year (Scenario 2017)
based on scenariol 208
6.8b Markov trarsition of land use/land cover for 2û-year (Scenario 2017)
based on scenario2 208
6.9 The 20-years change cha¡acteristic between Simulation and Scenario 2lt
6.10 Land use/land cover within 10 km of main roads tn 1989,1997 ,
Scenario2 2007 and Scem¡io2 2017 2ls
xl11
Declaration
This work contains no material which has been accepted for the award of any other
degree or diploma in any university or any other tertiary institution, and to the best of
another person, except where due reference has been made in the text.
I give consent to this copy of my thesis, when deposited in the University Library,
Dl:.
a
SI DATE:
(Asep Karsidi)
x1v
Acknowledgement
and Environmental
This thesis was carried out in the Department of Geographical
Professor Lesley Potter for her encouragement and enthusiasm in her role as
Financial support for this study was provided by
the
postgraduate coordinator.
scholarship Schema, and financial
University of Adelaide,s international postgraduate
was provide by Deputy of TPSA BPPT'
support for the research presented in this thesis
Enol Bamford and all staff in
I would particularly like to thank Dr. Kevin Harris, Mr'
many yeaÏs of friendship and
GISCA for their constant support provided tliloughout
xv
encouraging and truly a reminder
common word coming from my son, is very much
thatlmustcompletemystudyandgainthebesteverinmylife.
xv1
Chapter One
INTRODUCTION
1.1 lntroduction
Accelerating changes to the earth's environment are being driven by growth in human
changes in technology and socio-political organnation. Changes in land use and land
cover arc part of this transformation. Land use and land cover are closely interrelated
but are not the srrme. The term "land use" denotes the human uttltzation of the land,
while "land cover" denotes the physical state of the land. Baulies (1997) indicates that
Iand use and land cover change dynamícally due to demographic pressure, road
construction, fire, nutrient deposition, and many other factors. Land use/land cover
change may be the result of natural processes or human activities (Turner II et al.,
1994, Tucker et al., 1991, Sage, 1994). Furthermore, technological advances and
expanding population have put increasing pressure on scarce resources and have
created a variety of complex land use dilemmas that affect persons at all levels of
In the last decade, land use/land cover change has become a major issue at the global
and regional level. International attention has focussed on it because ofconcern over
1
Many studies have been conducted regarding the land use/land cover change issue
(Turner, 199I; Turner, I994;Mayer et al., 1994; Baulis, 1997; Thomas, 1997; Heilig'
1997; Rindfuss, 1998; Koning, 1999; Mc Connell, 2000; Solecki, 2001; Bicik, 2001).
Some of those studies focus on understanding the process of land use/land cover
change as well as developing predictive models of change at the global level (Turner,
I99l;Tumef, I994;Mayer et al., I994;Baulis, 1997). There have also been numerous
studies conducted at regional or country level looking at land use/land cover change
(Thomas, 1997; Heilig, 1997; Koning, 1999; Mc Connell, 2000; solecki, 2001; Bicik,
Land use/land cover change is also a major issue in Indonesia. Demographic pressure
with an annual rate of population growth between 1990 - 2000 at the national level of
I.49 % (Central Bureau of Statistics, 2000) and an uneven population distribution has
had an effect on the land use pattern in Indonesia. The 2000 National Census placed
Indonesia's population at206 million (Central Bureau of Statistics, 2000). Of this total
Java's population (Central Bureau of Statistics, 2000). However between the 1990 and
2000 census, the population in rural areas slightly decreased while the urban
population increased at the rate of 5.4%o per year over the period of 1990'2000. There
were 55.4 million urban residents in Indonesia in 1990 and 85.4 million in 2000
(Central Bureau of Statistics, 2000). Approximately one-fourth of this urban
population is concentrated in the Jakarta, Bogor, Tanggerang and Bekasi
population has had an impact on land use. Other land has been converted from rural to
2
urban use, from less intensive to more intensive agriculture, or from agricultural to
non-agricultural use. Nasution (1999) reported that about 20 000 ha per year of
agricultural land, has been converted into non-agricultural land mainly in the areas
In terms of the land use and population growth relationship, Sandy (1977) produced a
land use evolution model for Indonesia. People will start to cultivate an area at 25 m
above sea levef which has a gentle slope and is free from flooding. When population
increases, they will extend to cultivate the areas to the upper slopes and with more
population growtt¡ the area that could be conserved (such as forest in the steep slopes
and wetlands or swamp forests near the coast) will be cultivated, and as a result this
marginal land will deteriorate. This simplifred model shows that population pressure
can create ecological stress and land resources crises. Many other studies (Clark, 1977;
WolmarU 1993; Haub, 1996) have revealed that where there has been a long history of
population growth, there have been substantial effects on land and environment change
(Wolmar¡ 1993). The evidence suggests that the combination of a large population
base, a relatively rapid rate of population growtl¡ and rapid rates of technological
change has an effect of increasing the rate of change on land and in the environment.
One ofthe major impacts is associated with the structural economic changes producing
a shift of people from rural to urban areas, leading to a greater proportion of the
population living in urban centres and a greater share of land use being in urbanized
J
There have been numerous studies on the land use issue in Indonesia (Sandy, 1975,
studies have mostly focussed on making an inventory of land use conditions and the
causal relationship between social and biophysical factors. Very few studies have
sought to understand the dynamics of land use/land cover change over time and space
and to predict future changes. In this context, knowledge on the causes and eflects,
magnitude, spatial and temporal distribution of land use/land cover change dynamics,
as well as change prediction in Indonesia is far from complete. One of the problems is
the limited availability of spatial information regarding land use/land cover conditions.
There are a range of techniques for collecting information on land use/land cover
conditions, including field surveys and mapping from aerial photography and satellite
remote sensing. These techniques have various advantages and disadvantages. Field
surveys can provide extremely accurate data but are very time consuming and usually
impractical for large areas. Accurate mapping can be achieved from aeÅaI
scale (Lunetta et at., 1999). Satellite remote sensing with various spatial resolutions
analysis (Foster, 1935). The high frequency of satellite overpass makes it relatively
easy to detect changes in land cover such as those resulting from land use changes.
Various studies have been conducted regarding land use/land cover identification using
remote sensing techniques (GordofL 1986; Jadhav, 1993; RanU 1993; Dobson, 1996;
4
Karsidi et al., 1997; Henderson, 1997 Foody, 1999; Morisette, 1999; and V/ard;
2000).
The availability of accurate spatial information of land use/land cover, integrated with
vector data in a GIS, can facilitate frirther spatial analysis and model development.
Spatial analyses of land use/land cover change, such as where and what type of land
use/land cover has changed, and to what extent this change relates to social and
biophysical factors including why, is important for spatial planning purposes and to
ensgre scarce high qualþ agricultural land is not taken out of production prematurely.
manipulate, analyse and dispþ spatialty referenced data (Burrough, 1986). A common
approach of GIS is map overþ, which is capable of integrating many kinds of social
and biophysical data þers. This approach has been widely used in many spatial
analysis studies (Pathan et al., 1993; Schmid et al., 1995; Zeidle4 1997; Li, and Yeh"
1998; Cox, 1998; Debinski, 1999; and by Frederick et a1.,2001) but has beenrarely
done in Indonesia. However, this traditional GIS technique cannot effectively handle
the dynamic processes involved in land use/land cover change. Other techniques such
as Cellular Automata (CA) are a potential solution to handle dynamic process when
coupled with GIS techniques (Shi and Matthew, 2000). Cellular Automata are capable
of representing the spatial interaction among different systems in the real world and
5
This study explores land use/land cover change dynamics in an Indonesian region
using detection and spatial analysis via remote sensing and GIS techniques. It uses this
as a basis to predict future land use/land cover using a GIS based Cellular Automata
'West
model in the downstream areas of the Ci Tarum watershed, Java, Indonesia
(Figgre 1.1). The area includes the large-scale, multipurpose Jatiluhur reservoir. The
construction of the reservoir has created a dynamic response in the land use/land cover
change. Equalty important, this area has been heavily influenced by the expansion of
the mega-city of Jakarta and surroundings (Jakarta, Bogor, Tanggerang and Bekasi or
JABOTABEK). The development of new housing and other urban land use associated
with urbanization is a powerful force of change in the study area. Onthe eastern side of
Jakarta, where the study area is located, there has been dynamic e4pansion of housing
and industrial estate development. This is especially evident in the Bekasi kabupater¡
which separates the study area from Jakarta where there are many new housing estates
adjacent to industrial areas such as the Tambun industrial zone, Cikarang industrial
6
INDONESIA
c
F
Ð gtæ¡
Ë
gF î
A
WEST JAVA
Serang
Bogor
Garut Tasik
Kabupaten boundary
Study area 30 0 30 60 Kllomelers
Wateshed
use/land cover change using remote sensing and GIS techniques in Indonesia through a
case study in the downstream area, and surroundings of the Ci Tarum watershed in
West Java.
7
The study consists of,
It is hoped that by understanding more about the spatial phenomenon of land use - land
cover change dynamics, and developing and applying techniques to detect and predict
land use/land cover change using remote sensing and GIS, it will be possible to reduce
the misallocation of resources and mismanagement of land use in this area. Finall¡ it
policies that will improve land use planning and land resources mrinagement in
o To detect and determine land use - land cover change dynamics using
o To develop and apply methods that can be used to detect and predict land
8
1.3 The context of the studY
(Tomascik et al., lggT), has experienced significant land use/land cover change. The
issue of change becomes more critical, as the population increases, leading to the
increase in the demand for space. At the 2000 Census the population of Indonesia was
206,265,000 but it was unevenly distributed between the nation's 31 provinces (Table
l.l and Figure 1.2). Of the population about 60 Yo are located in Java and Bali, which
together have about 9 Vo of Indonesia's land. The average annual population growth
Land use/land cover change, including from deforestatior¡ has been of significant
magnitude in the last 50 years (Fredeick et al., 2001). Forest cover has decreased and
this is related to increasing population density. In Java and Nusa Tenggara, which
have high population densities, forest cover was less than25 Yo of the total area, while
in Kalimantan and lrian Jaya where population density is lower, forest covers more
Recent studies (World Bank, 2001) indicate that because of land utilization policies,
coupled with the impact of climate variability, a much larger magnitude of change in
forest cover has been experienced in Indonesia than elsewhere. According to a remote
sensing study, estimates showed that between the years 1985-1997 the rate of forest
loss was 19,700,500 hectares or 1.5 million hectare every year. It was 119,700,500
9
Table l.l:Distribution of Indonesian Population and Growth Rate, 1980 - 2000
10
UTAEA
Ê
\ Ê
KEPULÂUAH Êtau
t
)t' U-ALUKU ÜTNNA
.r
T &luß
\ KALIMAI¡TAI¡ JÂYÂ S
BA.RÂT
u
AJi
s4-t|G Êt"r q
TE¡{GAH
.l
\ O,ffi.
BEt¡G untu
uþ \
S€LATÂN *{ 'ql
A
ßAH 9ULASESI
çl'- SELAT¡H ¡t
Ê
?
8A¡fT EH æ '-o . ói {
10Êf
JAJflATT¡Uß
.ty'-
5dl 0
11
Sumatera
I
Sulawesi
F
ryv'fl*"
Ë
ava &.lvladura
"
þ
Bali-
oy'o Fo¡ect covet
À
fl10 0
o/o Non Forest
t2
Table 1.2: Forest cover and Population density in the main island of Indonesia
Year 1996
In the more developed islands of Java and Sumatra, land use/land cover changes
because of land use conversion have been especially pronounced. Nasution (1999)
estimated that between 1985 and 1995, the rate of land use conversion from rice field
to non-rice field in these islands reached 20,000 hectare per year. In the islands of Java
and Bali alone, more than 223,000 hectares of rice fields were converted to other uses
between 1981 and 1999 (Table 1.3), with nearly 35%o of thts involved in conversion to
and National Land Agency, 2000). Other estimates indicate that in the areas
surrounding the National Capital of Jakarta, the conversion of prime rice fields to
housing estates amounted to 2,000 out of 23,000 hectares in 1986 alone. Between 1989
and 2000, settlement in kabupaten Bekasi and Karawang increased, while rice fields
decreased (Table 1.4). The increasing density of settlement in this area includes an
13
Table 1.3: The of the Iture area 1981 to1999 rice field to other land in
Percentage 27.62 7 27
Source: Tabutation from Directorate Extensivication of Agricufture area Dep. Of Agriculture and National Land Ag
t4
Table 1.4: Area of rice field and settlement within kabupatens in the study area
Land use/land cover in the study area changed dramatically not just because of the
influence of the expansion of the mega city of Jakarta but also the influence of the
created a dynamic response in the land use/land cover change pattern. The flood
control function of the reservoir has made more land located in the downstream area
being freed from annual floods. It makes areas suitable for permanent use such as
settlement. The inigation. function of the reservoir has made water more readily
available all year round, making many of the downstream areas more suitable for
intensive agriculture. Since Jatiluhur was constructed, rice fields with technical
irrigation have increased due to the efÊect of the reservoir. The increase in rice fields
with technical irrigation continued until 1994, but it has been decreasing slightly since
then. The other categories of rice field such as rain fed, semi-technical and non-
technical irrigation, have decreased since 1989 (Table 1.5). Hence the decrease of rain
fed, semi-technical and non-technical irrigation are partly due to conversion into
technical irrigation rice field, but since 1995 alltypes of rice field including technical
irrigation have decreased in area due to conversion into settlement. Settlement has
l5
Table 1.5: Decreasing rice field and increasing settlement in V/est Java 1989-2000
The decrease in technical irrigation rice fields since 1994 is related to increasing
industry and housing development. This is clearly seen in the kabupatens within the
study area as shown in Table 1.6. The installation of hydropower plants has provided
electricity supply for industrial development along the watershed. The increase of
companies. Table 1.6 shows that in four kabupaten within the study area, the number
of national industrial and trading companies between 1989 and 2000 increased, with
I6
Industrial development has also involved establishing industrial estates such as those in
Cikampek, Karawang, Cikarang, and Bekasi. From Figure 1.4 it can be seen that new
developing housing and industrial estates were concentrated between Bekasi, Cikarang
and Karawang along the main road of Jakarta-Purwakarta a¡d toll way Jakarta-
Cikampek. The Karawang integrated industrial estate alone covered 2000 ha, and6234
thiln Fl( I
¡¡Ëffinl.d
llmnH*l ¿ffi1¡
låitr?É¡
hr
@ New Housing
o Indushial Estate
N
À
@ Industrial zone Iln Scaled
New Town/industrial estate
Figurel.4. Housing and Industrial estate developmertt area in Bogor, Bekasi and
Karawang (Adoptedfrom Jayadinata, 1999 pp- 185)
17
Table 1.7: Number and area of Industrial Estates in Kabupaten Bekasi year 2000
method to detect and predict land use/land cover change dynamics using
l8
This study includes three main elements: remote sensing, existing data collections and
geographical information systems to detect, analyse and predict land use/land cover
change (Figure 1.5). Figwe 1.5 presents a general overview of the conceptual
framework employed in this study. Remote sensing techniques are a way to detect land
use/land cover change on a digital basis. The integration of remote sensing and GIS
combined with existing spatial data is used to analyse land use/land cover change and
the possible driver factors as well as to operationalise a land use/land cover change
prediction model.
Identification and change detection of land use/land cover was based on Landsat TM
cover type. Landsat TM images were available from 1989 to 1997 anrrually. These
images enable us to detect land use/land cover change on an annual basis, and the
influence of weather or season variability within this time interval can be recognised.
The spatial analysis of land use/land covers change, and its relationship with the
possible driver factors, was based on examining both static and dynamic driver factors.
Static factors such as slope, elevation and physiography, and dynamic factors such as
population densit¡ proximity to urban and semi-urban centres, proximity to roads and
proximity to toll \¡/ay, were the main focus of this analysis. The relationship of these
driver factors with land use/land cover was assessed in order to recognise the factors
that have a strong influence on land use/land cover change in the study area.
l9
Land use/land cover prediction is another focus in this study. Future land use/land
cover information is important since it can be used to identiff areas that require
priority attention or to anticipate mismanagement or misarrangement of land resources.
The Markov Cellular Automata (MCA) model was selected to predict future land
use/land cover in the study area. This model was based on a combination of cellular
automata and Markov chain models that can represent the spatial distribution of land
The ¿n't anil sc¡unLte af'collecting Mcr¡ss ttrttl tùrtct c¿illec:Íed by ather
b),¿t ¡levice, which is not in conÍücl i n¿{ ivi ¿lt rct I or agenc i es
i n tu grü ri on, rnani pu Ia Í i o n, an a lys i s' ctnd d i sp I cty aJ' spa ti a I d ata
a qualitative ¿pl)'sis
Overla¡' and (or,erla1,, buffer &
cross tabulatiorr intersectioÐ Markov basecl
a quantitative analy'sis Cellular Automata
(discriminant
frrnction)
Figure 1.5. Research framework
20
1.5 Thesis Outline
This thesis consists ofseven chapters. A general review ofland use/land cover change,
as well as defuritions and some aspects related to land use/land cover change
dynamics, are presented and discussed in Chapter 2. Clnpter 3 consists of a description
of the socio-economic and the biophysical characteristics of the study area and the
availability of spatial data for it. This chapter provide a description of socio-economic
and the physical characteristics within the fives kabupatens that cover the study area.
Firstly there is a description of population growtt¡ occupation, and land tenure as well
climate and hydrology are discussed. Chapter 4 discusses digital change detection
using remote sensing and GIS. Change detection techniques are intensively reviewed
detection from the satellite imagery. The dynamic change of land use/land cover
related to the season and growing cycle in the study area is discussed. Chapter 5
identifies possible driver factors in land use/land cover change. Static and dynamic
drivers using overlays as well as statistical analysis are discussed in this chapter.
Cellular Automata model that are available in ldris32 were conducted to get a suitable
prediction result. Future land use/land cover change in the study area was predicted or
simulated 10 and 20 years ahead using this model. Finally Chapter 7 presents the
conclusions of the study. It includes a sunmary as well as some major findings and
some implications for policy makers and planners in land use panning. Some
recommendations for future research are put forward in the last part of this chapter.
2I
Chapter Two
2.1 lntroduction
Human factors such as population growth and distribution, economic growth, and
physical factors such as topography, slope, soil type, climate and others strongly
influence land use/land cover changes (Skole and Tucke4 1993). Land use change is a
matter of historical process as relating to how people use the land. Nowadays, human
attachment to land is not just as habitat and living space, but also involves more
complex purposes such as industry and tourism (Mather, 1986). Land use/land cover
change is intricately linked to the dynamics of human activities. The dynamics of land
use change is interrelated between demand and supply structure of land uses. Land use
demand increases with population and elements of economic growth such as income
per capita, GNP and industrialisation. On the other side, supply of land as space to be
This chapter presents an overview of land use/land cover change dynamics, providing
definition and reviews some aspects that relate to land use/land cover change.
precisely what constitutes land use or land cover, or how to define them (McConnell
22
and Emilio, 2000). There are many definitions and descriptions of land use and land
cover, depending on the purpose ofthe application and the context ofthe study. It is
hence necessary to define the definitions and descriptions of land use, land cover, and
land use/land cover change terms that are used in this study.
In general land cover is the biophysical state of the earth's surface (Turner et al.,
1995). It may be defined as the observed physical (including the vegetatior¡ natural or
planted) and human, constructions which cover the earth's surface (Baulies, 1997).
McConnell and Emilio (2000) add that water, ice, bare rock or sand and salt flats or
similar un-vegetated surface, although strictly speaking part of the land (and water)
itself and not its coveÍ, are for practical reason often included in land cover. Moser
(1996) notes that the term "land cover" originally referred to the type ofvegetation that
covered the land surface, but has broadened subsequently to include human structures,
such as buildings or pavement, and other aspects of the physical environment, such as
Land zse refers to both the manner in which the biophysical attributes of the land are
manipulated, and the intent underlying that manipulation or the purpose of the land
used (Turner et al., 1995). Moreover, Iand use involves considerations of human
the interlevel integration of processes at one level with those at other levels of
manipulation of the land to provide food, building materials, buildings site, firewood,
and clothing.
23
A second definition of land use is the purpose for which the land is being used (Young,
1994). This definition includes the management of land and human activities which are
directly related to the land. In addition, land use constitutes a series of activities
land, carried out by humans, with the intention to obtain products andlor benefits
The term land use/land cover has commonly been used in association with land use,
particularþ in the field of remotely sensed analysis. Land use changes aÍe a major
determinant of land cover change. Land use/land cover inventories tend to be a mixture
of land use and land cover. One category of land use may correspond fairþ well to one
class of land cover. On the other hand, a single class of land cover may consist of
multiple uses. There are two categories of land use/land cover change. conversions
from one land use/land cover into another, and transformations within a given land
use/land cover type. Lambin et al. (1999'. p.37) refer to the process by which land
. The activities (or operation) and inputs that are undertaken (or
and
outputs (goods or services) that are expected, and theforces that cause
24
In terms of detection or identification, land use/land cover is the observed physical
cover at a given location and time, as might be seen on the ground or from remote
whereas information on land use requires a statement of purpose from the person who
controls or carries out the land use. Remotely sensed data (e.g. from satellite images)
can often be used to map land uselland cover, for example, by identifying multi-
spectral signatures characteristic of certain land cover types. Land use, in turn, is often
related to land cover, so that land cover may be used to infer land use (Adger, 1994).
Land use/land cover can change due to natural or human factors. Natural factors such
as flooding, drought, forest fire, and volcanic eruption result in a change of land
use/land cover. Human factors such as demographic pressure, level of poverty, and the
economic and institutional structures of resource use are an indirect cause of land
growth, human activities have occupied more land (Mather, 1986). Changes in the way
which populations live their lives will change the function of the land, and the human
factors such as population and their activities are a driving force in changes in land
use/land cover. Sage and Grubler in Mayer (1994) argue that the main driving forces of
land use - land cover change are population, income and technology. Robinson (1994)
argues that the initial driving forces to be considered are population, income,
technology and price. Hence population and their activities, such as income and
technology, are forces shaping land use. But recognition of the link between human
25
activities and land use/land cover type is not simple. Single or multiple factors may
Land use type changes from forest to agriculture or from agriculture to settlement or
from agricultural to industrial or cornmercial, for example, are the result of changes in
human activities. On the other hand, the nature of soil such as mineral content, texture,
wetness and dryness as well as topography are factors that affect land use/land cover
change. Therefore, the interaction of human and biophysical factors influences the
dynamics of land use/land cover change. This interaction is complex and needs to be
simplified and categorized in order to identify the most influential human and
biophysical causes of land use/land cover change. McNeil (1994) concludes that there
are four major driving forces - political, economic, demographic and environment
(Table 2.1). Within these broad categories, he identifies specific attributes that strongly
influence land use/land cover patterns. Land uses for settlement, for example, may be
levels and income per-capita or may act together with other factors such as economic
growth, policy or technology. In Indonesia, land use change is related to all these
factors. The policy of privatisation of forest concession companies for example forced
forest conversion in the main islands such as Kalimantan, Irian and Sumatra (World
Bank, 2001). In the islands with a high-density of population such as Java, land has
been converted from agricultural to settlement and industrial use. From Table 2.2 it can
be seen that a high rate of conversion from agricultural to settlement and industrial
areas has occurred in West Java, Central Java, East Java and Bali. This is evidence that
26
population or demographic pressure is one of many factors that influence land use
change in Indonesia.
1.Political:
Decision-making process of public participation nitary of federal structure
Local pressure, special ed, centralized/ umber of special-interest
interest and corruption
State capacity s ector pressure/influence sector expenditure/
land area
2.Economic:
Vulnerability to external vs. closed economy
pressure(economic/p olitical) concentration
sector dependence sector export/
export
of exchange rate management (Population -
active)
total
exchange rate
service ratio
3.Demography:
Population pressure on High-low land
the land in population
agriculture/PEA total
4. Environment:
Natural resource quality Scarcity lr,o.u, yield, flow
Source: Adapted from McNeil ( 1 99a.p. 60)
27
Table 2.2: Change of agriculture area to other land use between 1981-1999 in Indone
Source:Tabulation from Directorate Extensivication of Agriculture area Dep. Of Agriculture and National L
28
The forces driving land use/land cover change can be distinguished into two categories
- direct and indirect. Direct drivers are associated with activities which directly interact
with and modify the physical environment, such as deforestation, urbanization, and
agricultural expansion (Turner and Meyer, 1994). The indirect drivers, or root causes,
influence how individuals or groups interact with, and change land use/land cover
@laikie et al., 1994). These are generally more complex as they are built into the
human system underþing a land use activity (Adger and Brown,1994; Krummer and
Tunner, 1994). An illustration of driving forces related to human systems can be seen
in the demand and supply relationship of a land use change model @igure 2 l). The
demand structure of land uses is very dynamic; it means that the demands for using
land for human activities increases from time to time under the influence of population
growth, community structure and economy. On the supply side, structure is relatively
static, the surface area is constant, and elasticity from this supply side is on the aspect
of physical characteristic and function of the land. Figure 2.1 shows the mutual
dynamics of land use change. This flowchart explains the increase of land use demand
demand for goods and services, which in turn affects the increased demand for
productive land use. The flowchart also explains the opportunity of increasing on the
natural quality of land (soil type, slope, altitude and climate) are the factors that could
lead to an increase on the supply side. Institutions and policy at national, local and
community level have an influence on technology and land use decisions. The
29
technology level will affect the capability to cultivate and employ the land
(Saetulhakim, 1999).
Danographic pressure;
Population Crrowth Institution and
and shift Policy
on Change of
ofGoods socio-economic
and and
ptrowth
Increase
demand
Dynamics of
Land use
Increase of supply;
A¡ea constant
Elasticity factors
supplies are:
in land
and
and spatial
landscaoe
land
Biophysical Technology
Characteristics Achievement
Change ofsupply
Elasticity of land use
change
Demographic pressure such as population growth and shift are a driving force of land
use/land cover change which is unique and can be quantified (Meyer, 1992). However,
population growth is not the main driving force as its importance is relative to the other
30
forces generating land use/land cover change such as technology and social advances
that improve the conditions of life. Young et al.(1991) did not find a significant
correlation between the land use/land cover changes with increases in population
density at national scale across the world. Hence increase in population density may be
sufficient to cause land use change but it is not necessary. Significant correlations
between population and land use/land cover change have been found when
characteristics. Comparative studies offer statistical evidence supporting the claim that
A study that analyses the interactions of population growth and land cover change
@ilborrow, 1992) concludes that population growth is an important factor, but one
significantly modified by natural and institutional context Land use/land cover change
is closely associated with agriculture, but not clearly linked to population (Meyer,
7992'. pp. 54). Empirical evidence higtrlights the complexity of factors concerning
population pressure and its links with agricultural extensification and intensification.
These would appear to imply that whilst the initial response to increases in population
access, or only very marginal areas remain, then intensification occurs. Turner et al.
p.38a). Hence land use/land cover change is usually an outcome of population growth.
Moreover, population change has a direct association with change in the way people
lead their lives and their activities, including agriculture activities. Population growth
31
is argued by some to have exceeded the capacity of the biosphere, as managed by
society, to sustain it @hrlich. I., 1990 and Ehrlich. I. 1988 cited by Meyer, 1992:
pp 52). Growth in population and economic activity require land not only to fulfil the
basic needs but also to ensure needs in social life and leisure activities are met. Thus,
the more people, the more land needed for agriculture to produce enough food and for
Land as space remains static and demand to use the land rapidly increases as a result of
population growth. The more people the more space that is needed for their living and
population growth can bring land use conflict. Land use conflict frequently occurs
between agricultural and non-agricultural uses. This conflict occurs because of the
basic principle that each piece of land should be devoted to the use in which it would
yield the highest rent (see Hoover and Giarratani, 1999: pp.7). Therefore demand from
the non-agriculture sector is usuaþ much stronger than from agriculture due to the
much higher incomes which this sector can obtain from the same area of land. If this
condition occurs in highly productive areas for agriculture it can have an effect on
overall food supply. There is a significant correlation between the amount of cultivated
land and food supply (Li and Hong, 2000). The conversions of agriculture area into
non-agriculture area will other have an effect on regional food supply. This is because
urban areas lead to be located in the most productive agriculture area. This situation is
this time the rice fields began to decrease. In 1983, the area of agricultural land in Java
32
was about 5.42 million ha, but after 10 years decreased to 4.41 million ha or decreased
From 4.41 million ha agricultural area, Java, supported about 63Yo of domestic supply
of rice, the remaining 38Yo supported by 11.04 million ha from out of Java. The
decreasing of rice field in Java has an effect on decreasing rice production and rice
supply at the national and regional level (Table 2.4). Therefore, to compensate the
decreasing rice field in Java, the government of Indonesia has a long-term policy to
open one million ha rice field in peat area in Central Kalimantan under President Decry
JJ
Table 2.4. The average rice productivity in Java and offJava (Ton/halyear)
can be seen that the highest conversion of land from rice field to settlement and
industrial areas was in Java. It was 59125.32 ha change from rice field to settlement
and 15600.87 ha to industrial area. While in other islands such as in Irian Jaya and
Maluku there was no conversion from rice field to either settlement or industrial areas.
34
2.4.3 Land use, lndustry and Manufacturing
The demand for land by industry and other commercial purposes has strongly
expanded because of the rapidity of industrial development. Many factories have been
constructed and spread into the countryside in Java (Chapin, 1972; Benstein, 1994;
Hidding, 2002). The number of large and medium manufacturers in Java has increased
with about 80o/o located in Java (Table 2.6). Often industry requires more land than
needed at the time to provide for subsequent growth. In effect, a lot of agriculture or
6. Extractive industries.
35
Each category has a different interrelation with land use/land cover change. Office and
allocation for this category has an effect on urban land use change (Chapin, 1972;
Rhind et al., 1980; Dawson, 1984; Cloke, 1989). However, in practice, there are
nowadays, that industries look at fresh sites in the countryside to build a factory. This
condition has an effect on land use/land cover change, because construction of new
factories needs new roads, communications and other facilities. In Indonesia, for
example, new towns are often built adjacent to industrial estates @igure 2.2). Lippo
Jakarta adjacent to industrial estates and service-oriented facilities such as hotels and
Tanggerang and Bekasi). The development trend becomes clear when the land area
used for some new town projects is considered. For instance, to the west of Jakarta in
the Tanggerang district, the new towns of Serpong has taken ó000 ha and Tigaraksa
3100 ha. To the east of Jakarta in the Bekasi district, Cikarang Baru adjoining
Cikarang stretches over 5400 ha, and Bekasi'2000 an area of 2000 ha. These new
towns are adjacent to the industrial areas such as Bekasi integrated industrial estate,
Tambun industrial zone, Cikarang industrial zone and Karawang integrated industrial
36
I(AÀTA¡FERA¡E DrqJÆqwn XA¡.E(ASI r(A¡.u( Rr$ro
liler!Í¡r{t
ft!fiì¡,lffir
8h OF¿Ctt ¿IEit
Tbg{rÌ{. l.ùOln ,fn¡tÉrtl
r#h*t0
hlaooh¡
O.FJI¡.ÈÍu
n¡ry*n 0.H,6.¡F0lr
Tripr¡
$.@lf,
l{ñ FK
tþpoûreíR
anFllo?!0hr I
t l¡tÉt3Hl
I l(mùÉll t.ff!r
I EddbfÊ-Ú
I !¡fhJaÊ
lllrtna
s.ml'
eüOS4ù¡ füÇrLrôtf
fLa 'ñe, i.5@lÐ l(ntrglrDl.t¡
2mfr
@u.m;r
New Housing
@
a Industrial Estate
@ Industrial Zone A
New Town/Industrial Estate Un Scaled
therefore the behavour of the agricultural sector with respect to demand for land from
adaptation and others. Agricultural land use involves growth of food crops and other
cassav4 sweet potato, peanut and soybean. Other economically valuables crops consist
of tea, coffee, oil paln¡ tobacco and others. The extension of this agricultural land use
is constrained by climate and soil factors and also depends on the agricultural market
37
situation, changed production conditions, and government intervention such as
incentives and regulation (OECD, 1976: p. 3l). In the humid tropics, such as in
Indonesia, intensive cropping systems occupy most of the resource-rich land, fertile
soils, low slopes, and areas with adequate rainfall or irrigation. For example rice fields
are mostly located in the flat areas with good irrigation and fertile soil (Figure 2.3).
bon
Land use
üi;;i.; Mangrove
1ffi Forest
fl:!.î.Ti Seüement
Plantation estate
iliî.Wetland
+ÉÉË RicE field 2 cycles
üå+ Rice iield 1 rycle
¡p64 Water/reseruoir A
ffi
'l:i::;l+
Fishpond
30 g m Kilûm eåers
fren/dry land
i;::.;,:t frerúdry land
The other economic value crops such as coffee or tea are located in areas upland of
paddy fields. These crops mainly occupy fertile soil with gentle slopes at an altitude
within 500 m above sea level. In many cases areas with poor soils and water resources
have received far less attention in agricultural development, but increasing population
and development pressures are forcing conversion of these areas into more productive
38
and intensive land uses. These pressures continue especially in areas where there is
The issue of spatial scale is important in land use/land cover change analysis. Land
use/land cover change analysis at the world level, for example, is quite diflerent from
analysis at the national or local levels. Moreoveq the analysis of land use/land cover
change is essentially performed at the level of the land use/land cover classification
systems. These systems are tied usually to a particular spatial scale. For example, at the
world level, FAO distinguishes four or five major land use/land cover types
(Briassoulis, 2000). At the level of nations, the number of land use/land cover types
increases to around ten. At this level, each type was divided into a more detailed land
use/land cover typology- Agricultural land is further subdivided according to the type
At the national level of Indonesia, for example, the National Land Agency of Indonesia
divides land use classification system into three levels: level I, II and III (Table 2.7).
The level of this classification system relates also to the source of data, acquisition
techniques and spatial resolution or map scale. Anderson ef al. (1976) summarised the
level of land use/land cover classification based on data source and spatial resolution
into five levels (Table 2.8). Table 2.8 shows that Landsat MSS is capable of producing
land use/land cover classification to level I with map scale of 1: 250 000 and smaller.
39
Table 2.7 . The land use classification system in Indonesia
40
category such as built-up area, agricultural area, open/bare land and water body or
wetland.
Recently, satellites can produce land use/land cover maps with more detail or higher
classification level. For example, Landsat TM has a 30 m spatial resolution and can
produce a map with scale of l:80 000 or l. 50 000. This map scale can produce land
Indonesian system. New satellite such as Ikonos has a I m spatial resolution and can
produce a map with scale up to 1: 2000. Based on this spatial resolution, land use/land
cover classification system can be derived at a parcel level. At this level land use/land
These spatial scales have an influence on the level of land use/land cover change
observed. For a given time interval, land use-land cover change may not be discernible
at higher spatial levels while at lower levels (e g. at the level of a settlement), very
4l
large changes may be measured. The role of the land use/land cover classification
system is critical in this context, as are also the source and acquisition techniques of
Hence, regardless of the intended level of analysis, the results obtained and the ensuing
description of land use/land cover change refers to land use/land cover change at the
level of the spatial scale. The smaller the scale the broader the land use/land cover
spat¡al planning
In many countries there have been significant changes in land use/land cover as a
has been used for agricultural and urban purposes. The lack of collection of land
use/land cover data or information will have implications for land use planning or
spatial planning processes. Scholten et al.(1999) state that the need for land use/land
implementation Recently, with land use/land cover data available on a digital basis
from remoteþ sensed data and GIS, the lack of information on land use/land cover
use/land cover supported by a digital model of spatial data is shown in Figure 2.4. This
figure shows that when land use/land cover data are available in digital format, they
42
may be integrated with other data attributes to create suitability maps for each land use
type and finally to support the land use allocation process for each parcel. It is a
dynamic model because it deals with changes in land use taking into account the
Current land
use/land cover Selected Selected
attribute per grid spatial policy planning
cell optrons scenario
Exogenous
Grid cell Suitabilrty sector specific
attributes map regional model
specification results
Figure 2.4. Overall model spatial planning structure related to digital data of land
use/land cover. (Adopted and modifiedfrom Scholten et al., 1999)
Historically, land use inventories were initiated in 1930 when Sir Dudley Stamp for the
first time undertook a land use inventory of the British Isles (Stamp, 7931) In
Indonesia, systematic land use inventory started in 1960 (Sandy, 1995) as a follow-up
of the Master Agrarian Act (UUPA) 1960, and land use inventory has been continually
43
conducted since 1969 (Silalahi, 1999). Since that year, land use in all provrnces rn
Indonesia has been mapped. The result is a land use map consisting of information of
the pattern and spatial distribution of land use in all provinces in Indonesia.
Based on the spatial information of land use condition in all provinces, in 1977, Sandy
(1977) produced a simplified model of land use evolution in Indonesia called the
Indonesia". This schema shows the evolution of land use pattern from schema A to
schema I (see Figure 2.5). Schema A represents land use condition in area that humans
have not settled yet, and schema I represents the land use pattern on critical condition,
which is when some reserve area has been changed into marginal land.
1000
Figure 2.5
The concept of land use
evolution schema in
500
Indonesia
25
0 Í
1000 t ll Legend
ffi Forest
Shifting Cultivation
c
'*u**"t
25
ffi..E
tÉi=;:t
Swamp
Marginal land
0
44
This land use pattern evolution is clearþ recognized in Indonesia. One example can be
seen in Figure 2.6 below that shows the land use evolution in the Indramayu area, West
Java. Figure 2.6 (A) shows the land use condition as the same as schema D and E in
1857 when people start to develop their agricultural practices. Some dry agriculture
areas have been converted into rice field and rain fed areas with I cycle and improved
to 2 cycles rice field when supported by imigation. Figure 2.6 (B) illustrates land use
condition that has been changed similar to schema G and H by 1940. As an effect of
population growth and improving agricultural practices, shifting cultivation has ceased,
and some swamp areas have been drained and used for other agriculture practices.
Land use condition the same as schema I is shown in figure 2.6 (C) for a depression
area and Figure 2.6 (D) at high elevation area (mountain) in 1969. At this stage some
marginal land present in the upper and lower areas begins to show effects of overuse.
í€IERåilÊè8. ¡s¿t
ffiñr
kiÈùÈ\96ãr
e*"-a, f"til{ f.åHtliffi
c i êru
*sÈr ts49r.
t. a 3&
å +å¡q* ì.ìe,Þe3
1sìù ltre d-lÞs¡q- .î4þìêe
ÞeãÊ{4 G-11
T¡u.,ç S;i$ilô F* F
PAt{tÁl
tr
lrp
1tðc
Pã€L!¡li-ËtGAN iS Êã
q*-*--L.:*Jx{
'lAi-¡\P ÈäËiii& 3 *--+.
'West
Figure 2.6Land use evolution in Indramayu and its surrounding, Java
(Kartono et a|.,1998)
45
Silalahi (1999) conducted another study of land use change in Indonesia, which
focuses on analytic comparison between conservation and cultivation areas. The study
concludes that many factors affect land use change but the most significant factor is
population, specifically the density, growth rate and the percentage of poor people^
Land use changes in conservation areas are affected by population density and
The evolution of land use pattern that has been shown above and the study conducted
by Silalahi (1999) focuses on the general scale and does not contain actual figures of
land use conversion such as settlement and urban land use. The urban land use will
(2000) conducted a study of rural to urban land conversion in Indonesia. This study
concludes that the land conversion on the fringe areas of large cities in Indonesia in the
1990s has been largely due to the issuing of excessive land development permits
granted by the National Land Development Agency. This conversion mainly involves
land use change from agricultural into urban land use such as settlement, industrial and
other commercial areas. Competition of land use in the urban areas will have an effect
on land use in the countryside. The growth of Jakarta metropolitan city in Indonesia,
for example, has a strong influence on land use change in other areas such as Bogor,
Tanggerang, Bekasi and now through to Karawang and Cikampek as has been shown
inFigare2.2.
46
2.4.8 Land use/land cover change detection and spatial planning in
lndonesia
Documentation of land use as well as other biophysical spatial data was conducted in
Indonesia in the earþ 1960s (Sandy, 1995), but documentation in digital format began
project which was launched in the 1980s (Rais, 1997). Since this year several studies
and a pilot project dealing with land use inventory and spatial planning on a digital
basis have been conducted. Land Resources Evaluation and Planning (LREP), for
project began in 1987 and finished in 1993 and was coordinated by Bakosurtanal (The
and Planning (MREP) is another National level project to support spatial planning, but
providing spatial maps (base maps as well as bio-physical and land use/land cover
maps) in digital format. Another pilot project that deal with spatial planning is the
Land Use Planning and Mapping (LUPAM). This project, under technical aid from the
German government, was implemented in the 1990s in the East Kalimantan province
(GTL,2000). Some important information from these projects is that the land use/land
cover change detection should be pursued in tandem with a streamlining of the spatial
and the roles and required contributions from its stakeholders. The new paradigm in
the period of Reformasi year 1999 supported by law No.221I999 attempts to do the
47
opposite (bottom-up planning approach) Gru No.22l1999). Through this new
approach, the spatial planning process must be supported by adequate spatial data or
information at appropriate scale. The government regulation No. 10/2000 supports this
appropriate map that has to be used. This law guides the office or institution
responsible for spatial planning activities at any level (National, Province and District)
As can be seen in Table 2.9, different scales are to be used for different level of spatial
planning. To support spatial planning at the district level, for example, the map used
required by the law, is a map with minimal scale of 1: 50.000 and for the small
district/municipality a map of scale 1:25.000 or l:10.000 must be used. This new law
has implications for providing spatial data, include land use/land cover data, for spatial
planning process. Figure 2.7 shows the integrated planning system for regional
development in Indonesia, which uses current land use/land cover condition and other
spatial and non-spatial data as a basis of this spatial planning process. Spatial planning
covers physical planning and land use planning, which is closely related to, and should
48
provide the basis for, regional or local development planning and sectoral planning as
Regional
Cunent Land use/land
cove¡ condition Development
and others spatial and Planning
non spatial data
Sectoral
Physical Land Use
Planning
Planning Planning
Land
Resources
management
Spatial planning is the fundamental process to identiSr, analyse, and organise the
Physical planning is a more specific form of planning, which is the designing of the
railway, airports and - facilities for towns and other human settlements in anticipation
49
Land use planning is a process that facilitates the allocation of land to the uses that
provide the greatest sustainable benefits (UNEP, Agenda 21, paragraph 10.5).
This schema has been introduced by GTZ to avoid the confusion of terms that indicate
the differences in interrelations between the various plans and the inherent danger for
overlaps and hence contradictions, between the spatial plans, development plans and
2.5 Gonclusion
This chapter has provided a review of land use/land cover change dynamics as well as
land use/land cover definition and understanding, including some aspects related to
land use/land cover change. The evolution of land use pattern in Indonesia has been
presented to illustrate the relation of land use change and demographic pressure in the
study area.
Land use and land cover definition depends upon the purpose of the application and the
context of the study. Land use and land cover in this study has been defined as terms of
"land use/land cover", which is land use representing human-induced intervention or
manipulation of the land to provide food, building material, building site etc. Hence,
land cover referred to the biophysical state of the earth's surface. It may incorporate
the observed physical properties, including the vegetation (natural or planted) and
50
The interaction of human and biophysical factors influences the dynamics of land
use/land cover change. This interaction is complex and needs to be simplified and
categorized to identiff the most influential human and biophysical causes of land
use/land cover change. In relation to demographic pressure, land use/land cover change
is associated with population growth. Population has a direct association with demand
of land for their life and their activities, including agriculture activities. The initial
becomes more difücult to access or only very marginal areas remain, after which
intensiflcation occurs.
The magnitude of lands useJand cover change is different at the global, national and
regional level. It requires different scale maps or spatial resolutions. The more global
the view the more simple is the land uselland cover category to be analysed and a
The spatial information of land use/land cover is the most important part of the spatial
spatial planning must be supported by appropriate scale at each spatial planning level;
5l
Chapter Three
3.1 lntroduction
This study conducts an analysis of land use/land cover change in the downstream and
5" 54,- 7o18' S. Administratively, the Ci Tarum watershed contains all or part of the
Cianjur (see Figure 3.1). The upper watershed has a hitly and mountainous topography
Sumedang
Sukabumi
Cianjur
Garut
Å
20 0 ã] 40
æ Kllonôlars
52
This area was covered by two Landsat TM scenes 122164 and 122165 þath/row)'
Scene 122164 covered the downstream and surrounding areas and scene 122165
available
covered the upper watershed. Scene 122164 is the only satellite data that are
for 5 consecutive yeafs. Therefore, the area of study was focused especially in the
of the five districts (Kabupaten) Bekasi, Bogor, Karawang, Purwakarta, and Subang
This chapter describes the social and biophysical environmental conditions in the
selected study area aflecting land use/land cover change and explains
the availability
area of this study is only the downstream watershed and its surroundings'
The region of
settlement. This area has been a prime rice-producing region of Indonesia (World
associated with the expansion of Jakarta, a city of 9 million people (Figure 3'3)' Three
large reservoirs have been constructed in the region; one is located in study area
(namely the Jatiluhur reservoir). The other two are located in the upper Ci
Tarum
watershed - namely the Ci Rata and Saguling reservoirs. Jatiluhur is the largest
paddy/rice field
reservoir and is a multipurpose reservoir, which supports irrigation for
53
downstream and as part of a hydroelectric power generation and supply system for the
whole of Java and Bali, especially the growing mega city of Jakarta.
Bogor
Kabupatên boundary
Sludy area boundary
A
Downstream of Citarum \^Æ 20 0 20 Kilometers
Jatiluhur reservoir
Figure 3.2The five Kabupatens and Jakarta within the study area
-
Figure 3.3. The New Town and Industrial estates development in Jakarta and
surrounding areas. (Source: adaptedfromJayadinata, 1999: p 18a)
54
3.3 Socio-Economic Gharacteristics
kabupaten of the study area except kabupaten Bekasi that adjoins the special capital
city district of Jakarta. The settlements are characterized by ribbon development along
the roads. Currently, since the government of Indonesia constructed the Cikampek
highway in mid 1985, some real estate and industrial plants, including a golf course'
Java 59"
Sï41Ê roads
tæ
'l|ì1,Þj Fo re sl
SÉttìemÊnl
Plantation estete
Flicc fiùl ? cycloo
A
Figgre 3.4. Land use of 1993 in five kabupaten within the s]udy area and Jakarta city
Source: SARI Project BPPT, 1998.
55
3.3.1 Population and DemograPhY
The total population in kabupaten Bekasi, Bogor, Karawang, Purwakarta and Subang
in 1990 was 9,1 million and 11.3 million in 2000. The population growth rate in
kabupaten Subang is low, 0.93%lyr during the period 1990-2000 (Table 3.1). On the
other hand, population growth rates during the same period in kabupaten Bogor,
Karawang, Purwakarta were high; about 1.33 - 2.40 %lyt and extremely high in
kabupaten Bekasi (about 5.60 %lyr} For the whole kabupaten within study area,
population growth rates during 1990- 2000 was 2.41 %lyr compared with 2.03 % n
west Java and l.4o/o in Indonesia.
Sex Ratio rn/f t990 1.00 1.03 1.00 1.O2 0.98 0.98
2000 1.01 1.03 1.01 1.O2 1.00 1.01
Table 3.1 also reveals that the number of households (one household consists of 4 to 5
persons) increases more rapidly than the population. In kabupaten Bekasi, and also in
56
%olyr whereas the
growth rate. ln total, the number of households increased at 6'77
increases' the
population growth rate was 2.41%lyr. As the number of households
increases. Thus, the trend (i'e'
demand for housing units and settlement areas also
signals the pressures from
relative high growth rates of the number of households)
another forms of population dimension. This signal is the need for more spacious
housing units.
3.3.2 EmPloYment
populatior¡ a total of 4.4 mitlion n lggT and 4.8 million in 2000 were in the labour
high level of employment or
force (Table 3.2). This firther comes with a relatively
job'
increased in 2000 by 4.5 million working and 0.37 million seeking
At the kabupaten level, the role of the agriculture sector in offering employment
57
the municipality. This indicates that a lot of areas have been changed to industrial or
land use
other commercial use such as for hotels, restaurants and golf. Therefore,
Occupation bY Sector;% of
Workino labour force: ,l
152,67 6 3 -2
a. 147,41
3,39 -2 0,97 0,59 -0,38
M and 6,21
76,6 -11,32 25,24 22,6 -2,64
I
87
1,71 -1 ,4 3 -0,36
EI Gas & Water Su 3,11
Construction 31,95 21,28 -10,67 I 6 1 -2,82
21
oÃ
111,41 129,27 17,86
Trade Hotel & Restaurants
39,25 52,67 13,42 61 0,15
Tran
5,47 2,52 -2,95 1 I -1,02
h Financial lntermediaries
59,89 -7,08 27 27 0,2
Services 66,97
0 -0,3 1 -0,16
Others 0,3
100 1 1
Total 100
Source: West Java in Figures, 1997 and 2000.
especially in Java'
substantial decrease in the area of agricultural land in Indonesia,
number of agricultural
However, Table 3.3 shows that there was an increase in the
land in Indonesia was
households over the same period. The reduction of agricultural
58
other hand, the number of agricultural households increased by about 2.49
million or
of this
approximately 249,000 households per yeaf. The most important consequence
1993-Agrieultural Chanp
11,L08 1 488
Java
18,693 21,183 2,490
Indonesia
Source: Central Bureau Statistics
In
Apart from farm size, the distribution of land holdings needs to be considered'
less than 0.50
Indonesia, in 1993, the number of households with size of land-holdings
hectares was 10.94 million or 51.63o/oof the total agricultural households. In Java, the
59
distribution was even worse. In the same year, numbers of households
with farm size
less than 0.50 hectare was 8.10 million or 69.85%o of total households
in Java. The
by 795'000 or
number of households of this category (farm size<0.50 ha) has increase
productivity of the
The type of agriculture land use is also significant since it impties
most productive
land and, thus, farm incomes to the owÏrers and operators' One of the
facilities'
agricultural lands is wetland or irrigated rice farms supported by inigation
hectares in 1993; a reduction by 14.87% or 43,800 hectares per ye¿r' The reduction in
agricultural land, particularly in Java, gives more space for the development of
Change from
West Java Central Java* East Java Indonesia
land to:
15847.00 22405.45 20872.87 72994.12
where rice
In the case of wetland rice freld, the developments have led to a situation
such as horticulture'
farming no longer is competitive, even compared to other crops
60
crops (vegetable, flowers, etc.) and for animal husbandry (especially
broiler)' on the
people.
in the discussion
Land use and land tenure are two aspects that should not be neglected
area of inigated wetland rice field (i.e. technical, semi technical and non technical
hectare
irrigation); from349437 hectare in 1989 to 288108 hectare :lr'1997 and272604
The areas of dry land
in 2000; which confnms what has been discussed in land tenure.
compound) have
farm (i.e. dry field, grassland, temporary fallow etc' except housing
in land used for house
also decreased. The reduction of this area is due to the increase
6t
Table 3.5 Existing Land use in the five kabupaten within study area (Ha)
I 400
1200
1 000
(t
E +RlcêFieH
o 800 +Houæ comPound
8 +St.b Forott
È- +Pond & W€ü¡nd
(! 600
o
400
200
0
Year'89 Year'97 Year'98 Year'99 Year'20
Yoar
62
3.3.5 lnfrastructure
use allocation.
Bogor
¡o
Provlnclol Road
RoEd À
TollWay
way n 0 20 Kilometers
Adm lnlstratlve boundarY
63
in the
Table 3.6.below provides also some information on the infrastructures situation
study area. Road length in whole area has increased by 2334.94 kilometres
or 54'82 %ó
in an eleven-year interval (V/est Java Statistic ofüce). During the 1989 to 2000 time
Karawang
interval, there was a significant increase in road length in kabupaten Bekasi,
and Subang. This increase related to construction of the Jakana - Cikampek highway
(Toll way); this had an effect on construction of connecting roads in this surrounding
kabupaten.
ROADS (km)
Type of surface:
1219.00 3t3.45 307.32 238.41 2344.88
Asphalted 1989 266.70
1288.84 768.23 330.80 722.61 39s3.30
2000 842.82
298.20 169.00 79.97 228.37 949.54
Stoneigravel's 1989 174.00
121.70 315.32 552.09 72.20 139.90 t20l.2l
2000
113.00 s5.15 1s0.70 65.81 s80.24 964.90
Soils 1989
55.67 t204.48 0.00 63.30 1439.75
2000 116.30
Condition:
489.80 319.81 252.63 238.04 13 5 s.08
Good 1989 s4.80
r33.96 479.45 663.4r r14.70 1s0.ó6 t542.r8
2000
770.20 150.79 111.91 228.38 1547.98
Moderate r.989 286.70
404.11 77.83 204.06 342.2s 1248.3r
2000 220.06
2r3.00 20r.75 68.60 76.6s 5 80.25 It40.2s
Damaged 1989
719.6s r 139.1 26.99 32.40 276.20 2t94.34
2000
t6.20 110.60 93.9s 1191.00 0.00 I4tI.75
Heavily damaged 1989
t02.34 t756.57 51.84 2s7.50 2295.40
2000 t27.r5
CTRICITY
t53791 66001 90828 t0975 data 101
Installed
1617586 940905 5 38961 I o data 353 1069
64
Rapid increases of the road length suggest rapid development of industrial and
settlement areas located in the back-up area surrounding the highway (Karawang and
Bekasi) where the price of land is still relatively cheap compared to areas close to
Jakarta. Only with new road construction can the accessibility of those back-up areas
be enhanced meaningfullY.
The development of industry and services in the metropolitan area of Jakarta and other
cities within the study area such as Boor, Bekasi and Karawang, also places more
pressure on the demand for electricity. Factories, hotels, restaurants, resort areas'
recreational facilities and places, shopping centres, education centres, etc. demand
more electricity as compared to when such activities have not been intensiffing. Table
3.7 shows that the number of electricity customers signifrcantly increased in all
kabupaten except kabupaten Subang (no data available) from 910608 in 1997 to
climate, influence the variations in land use/land cover type both temporally and
topography and physiography are the primary constraints of agricultural expansion and
other land use practices. Therefore these parameters have to be known and recognised
65
3.4.1 The Physiography and Soil type
based coastal belt of alluvial plain, a central belt of volcanic mountains with a series of
inter-mountain basins, and a southern belt of limestone hills. The study area is a part of
north coastal belt of West Java. In more detail, this area can be divided into four types
*Rolling hill and Volcano".
of physiography: "FIat", "Rolling hi11", 'Volcarto", and
The flat area is mostly located within elevation between 0- 499 m above sea level; for
the rest, elevation varies from 500 to 1499 and > 1500 m. The volcanic area is mostly
located within elevation between 500 - 1499 mabove sea level (Figure 3.7 and 3.8).
,¿ a
s<dl.llñ I
| çl@ 50
Figure 3.7. The Relief of v/est Java (source: adapted fromHugo, 1981: p 20)
66
Java 5E"
Figure 3.9 shows the soil type in the five kabupaten within study area. Distribution of
soil type in this area was influenced by active volcanos in the south part of study area;
the type of soils are varied and underlined sedimentation of tuffvolcanic. The alluvial
soils mainly lay on the flat area in the north part, which was influenced by tuffvolcanic
sedimentation. Some andosol lay on the mountainous ateainthe southern part of study
area, while latosol and podsolic are dominant in the central of the study area.
67
Java 56"
-l
Soilfype
Kabupeten bündary
Figure 3.9. Soil type in five kabupaten within the study area
- lsources: RePPProt, 1986)
In general, .urest Java receives rainfall under the influence of the West Monsoon
(November to April) and East Monsoon (May to october). Due mainly to topographic
effects, the rainfall is lowest along the northern coast and tends to increase as one
moves inland to the more elevated areas. Figure 3.10 shows rainfall zones
of the region
68
ü
llol¡tto !
agriculture and hence density of settlement. The monthly rainfall regime throughout
the year shows the dry and wet season cycle in the study area (Figure 3.11). The wet
and dry season affect the growing season of rice. In the areas that are not irrigated, rice
has only one growing season. The growing season cycle is improved to two cycles with
the support of an irrigation system. The Jatiluhur reservoir has as its main function to
irrigate the agricultur aI area on the northern part, which is mainly wetland rice fields,
69
Wet season Wet
@
?m
È
-uÆ Dry season
I
ï0
¡m
'Eo
I
0
JFN FB [/lAR AFR iVE JN JI- ÊCÐ SF O(r IIP EB
tvtrtì
Figure 3.11. The monthly rainfall regimes in the study area
(Source: Meteorology and Geophysics Agency)
Cycle I
Bloc Inigating Planting
I October lst November lst January I st
II October 16st November 16st Ianuary 16st
tII November lst December 1st February lst
N November 16st December 16st February 16st
Cycle II
Bloc Inigating Planting
t March lst April lst Iunelst
tI Mæch 16st April l6st Iune l6st
il April lst May lst July lst
IV April l6st Mav 16st Iuly 16st
Source: fafifuhu¡ ComPønY, 1998
the rice
There are three main irrigation systems from Jatiluhur reservoir that inigate
part of
field area in the study area. One is the east Tarum, which irrigates the eastern
rice field afea $rch as in Subang and Indramayu. Another is west Tarum which
fresh
irrigates the western part of rice field area such as in Bekasi, as well as supplying
70
the rice
water to the Jakarta metropolitan city. The third is central Tarum that irrigates
Jâva c" rl
Bo0or
Irr ig'atiot t
Citarum river A
Jatiluhur regervoir l<ilem qlere
analysis and
The main aims of this study are land use/land cover change detection,
scales have an implication on the magnitude of land use/land cover change analysis'
at the spatial scale
This study is concerned with land use/land cover change analysis
7l
a spatial
constructed and converted into digital format and used as sources to construct
The spatial data were available from several government agencies such as
Energy. Each agency has its own form and content of maps. BAKOSURTÆ{AL,
maps for
which is the main agency to provide base Ílaps, for example, Topographical
the whole Indonesian authority at scales l:250 000, l: 100 000, and 1: 25 000' Other
Research, and the Department of Mining and Energy are the institutions
which provide
addition spatial data such as forest cover' soil types, and geology as
well as
1: 25 000 and l: 50 000, while BPN creates land use maps at scales of l:25 000 and 1:
at least every
50 000. These maps consists of land use type which have to be updated
5-year (silalahi, lggg). Table 3.8 shows the agencies that are
involved in spatial data
collection.
Several spatial datasets that were used in this study consists of administrative
(Kabupaten and Kecamatan) boundary, roads, rivers and irrigation of west Java
l: 250 000; soil type and physiographic condition maps from RePProt
for the year
72
1986 with a scale of l: 250 000; land use type for the yeat 1993 source from sARI
Agencies in Indonesia
Table 3.g. partial list of Land Resources Inventory and Mapping
Directorate of GeologY
Geological MaPPing Division
Ministry of Agriculture
Centre of Soil and Agro-climate Research
Directorate General of Estates
Ministry of Transmigration
Directorate õeneral for Site Selection Planning and Programming
73
3.5.2 Satellite imagerY
Seven Landsat TM images path/row (122164) taken at the Dry season and the late Wet
in this area
3.9). It is difficult to get imagery for the same times due to the cloud cover
by
- the acceptable cloud cover is < 10% (NRCT, 1997). These images were provided
This satellite data were acquired at NRCT's (National Remote Sensing Centre
of
selected sites within the study area. Landsat TM has 30x30 m spatial
resolution and 16
The satellite image for the Indonesian afea can also be acquired from the LAPAII
receiving
ground station in Pare-pare South Sulawesi. This ground station is capable of
and processing SpOT and Landsat TM data. These satellite data were useful to provide
74
information of natural fesources on spatial basis, which overpass the same area every
data is
2g days for SpOT and 16 for Landsat satellite. The limitation of these satellite
3.6 Gonclusion
The aim of this chapter was to describe the social and biophysical environmental
of the Ci
of the study. The selected study area is the downstream and its surroundings
an area of
Tarum watershed. This area was a prime rice-producing region and includes
was needed to
yeafs consecutive and free of cloud cover. The annual satellite image
land use/land
The region of interest of this study is mainly lowland, with the dominant
population growth is relatively high compared with population growth at the national
other kabupaten'
level. Population growth in Kabupaten Bekasi was higher than in the
people'
Kabupaten Bekasi was located near to Jakarta the city of 9 million
75
Several types of spatial data and non-spatial data that related to this study were
collected from different agencies. The satellite images were acquired for the
dry and
1997 '
the late wet season on annual base years: 1989,1990,1991,1992,1993,l995,and
76
ChaPter Four
4.1 lntroduction
factors'
Land use/land cover dynamically changes as a result of social and biophysical
Land use/land cover change when one land use is converted to another
is relatively
land use'
within urban land use or between food crop plantations within agricultural
There are many image processing techniques (called digital change detection
The selection ofan
techniques) that are used to identiff land use/land cover changes'
can
In Indonesia, which is a tropical area,land use/land cover is often dynamic and
change due to weather or climate conditions and seasonal growing cycles' Sharkov
such
(1998, p. 175) stated that the various components of the physical environmental
77
as climate, soil, vegetation type, weather and relief contribute significantly to land
use/land cover change. V/ithin these conditions, land use/land cover can be difficult to
short time scales. Although these techniques may achieve more detail, in terms of
spatial resolutior¡ they need more time to collect information over large areas. Land
use/land cover may change significantly while the mapping is still being undertaken.
Therefore, techniques that can quickly identiff and detect land use/land covsr, such as
The spatial accuracy of land use/land cover detection is also related to the spatial
which can differentiate the detail of urban land cover from vegetation and agricultural
This chapter presents a review of digital change detection methods as well as the
results of land use/land cover change detection in the study area. The results of the land
use/land cover change are then used as an input into the land use/land cover change
model described.
There are a number of methods for obtaining land use/land cover data, such as
mapping are direct mapping technþes that have the advantage of obtaining more
78
detailed information and can produce large-scale maps. The main disadvantage of this
technique is the time required to cover large areas. Other techniques include mapping
from aerial photography from aircraft or by remote sensing from airborne or satellite
and analysis of land use/land cover and land use/land cover change, both at alocalized
scale and for broader areas (Foster, 1985; Bryant et al., 1986). Some disadvantages of
over time when different interpreters are involved. Using remote sensing techniques,
land use/land cover can be recorded over relatively large area and with reduced time
between the completion of data collection and actual field conditions. However, the
disadvantage of remote sensing techniques is that the spatial resolution is generally less
In spite of its timitations, remote sensing ptays an important role in the acquisition of
land use/land cover data. Remote sensing techniques enable the collection of land
use/land cover data in a digitalformat and the easy integration and analysis with other
digital spatial databases. Table 4.1 shows the comparison of the terrestrial and remote
79
Table 4.1 The Comparison of the Terrestrial, Air photos and remote sensing
Techniques.
Source: From various literature: Jensen, 1996, Lillesand, 1987 and Lindgren, 1985.
Over the last two decades, remote sensing techniques have been increasingly
developed and intensively used for natural resource inventories including land use -
Iand cover monitoring (Weismiller et al., 1977; Angelici et al., 1977; Brera and
Shahrokhi, 1978; Gordon, 1980; Jensen 1986; Quarmby and Cishnie, 1989; Green ¿/
al., 1994;Liet a1.,1993). Since the first Earth Resources Technology Satellite (ERTS-
l), later called Landsat 1, was launched tn 1972, activities dealing with the natural
resources inventory and land use/land cover change detection using remote sensing
with this spatial resolution, land use/land cover types smaller than 80 m x 80 m could
not be recorded. Later, other natural resources satellites with finer spatial resolutions
were launched. Spot 1 was launched in 1985, Spot 2 in 1990, and Spot 4 in 1998, all of
which have four bands, two in the visible infrared, and one in infrared with a spatial
80
with a 30 m spatial resolution and has eight spectral bands including thermal infrared
band and panchromatic (Foster, 1985: p. A\. The resolution element (spatial and
spectral) of these later systems is more comparable with the size of features requiring
detection. V/ith this improvement, land use/land cover change detection and
monitoring using remote sensing can be more accurately conducted.
Many studies have been conducted regarding the application of satellite imagery data
Maximum Likelihood using the original Landsat TM5 band combination (4,3 and 2)
has been extensively used to identiff land use/land cover in the Southeast Asia region -
e.g. in Indonesia, Malaysi4 Philippines, and Thailand (SARCS, 1997). All of these
studies result in land use - land cover change maps at the scale of 1: 100 000. Another
conventional approach is Band Ratios, which have been used to assess the capabilities
demonstrated that finer scale maps can be produced from Landsat TM at scales
(NDVÐ and Principal Components have also been used to identiff land use/land cover
types (Byrne et al., 1980; Fung and LeDrew, 1985; Jadhav et al., 1993; Qi, 1993;
Elvidge and Chen, 1995; Salem et al., 1995;Lyon et al., 1998; Li and Yeh, 1998; Fung
al.,2002).
81
Table 4.2list some of the various techniques used to identifu and monitor land use/land
cover using remote sensing data.
Table: 4.2.The various techniques to identiff and monitor land use/land cover
condition
82
4.2.1 Digital change detection of land use/land cover
Land use/land cover change detection aims to measure different features of land
use/land cover categories at two diflerent dates. Using large-scale aerial photography
at two dates of the same geographic area enables land use/land cover change to be
detected, but visual interpretation has some limitations and it is diffrcult to undertake
subsequent analysis due to the non-digital nature of the data. The digital change
detection method uses digital imagery to correlate and compare two sets of images to
identiff changes. Four main conditions are recoÍrmended for successful change
Singh (1989: p. 990) categorizes the digital change detection techniques into two
categories. The first is the independent analysis at different dates, and the second is the
change detection are still considered the basis of digital change detection and are called
change detection approach, two different date images are independently classified and
labelled. From the result, the areas that have changed are extracted. The pre-
83
different spectral signatures over time represent a change in land surface conditions.
This technique involves the transformation of two original images to a new single-band
or multi-band image in which the areas of spectral change are highlighted (Yuan et al.,
1999: p.22). The advantage of the pre-classification approach is the ability to eliminate
The disadvantage of pre-classification is the need for image co registration and the
difficult in discriminating "change" and "no change" areas. Hence, the disadvantage of
image.
Recentl¡ many studies have used the combination of pre-classification and post-
overcoming the complexity of land use/land cover change detection (e.g. V/ard et al.,
Table 4.3 shows the summary of digital change detection categorised by pre-
approaches. All of these techniques can include three basic approaches: pixel level
approach, feature level approach and object level approach (Deep, 1999). Pixel level if
detection refers to a numerical value attached to each band of each pixel in an image
such as digital number (DN). Feature level refers to some transformation of pixel level
84
values, and object level refers to symbolic or thematic labelling or classification of the
data. From this point of view pre-classification can be classified as feature level or
pixel level differencing, whilst most of the post-classification belongs to an object level
differencing , and ahybrid approach is a combination of feature level or pixel level and
The common digital change detection techniques that were categorised as Pre-
85
Table: 4.3. Comparison of digital detection techniques Categorized by Independently, Si
86
Table 4.3 continues
Categories Method Procedure/ Required Advantage
Multi date ingle Classification rectification and Efficient method
lassification two or more dataset correction Direct change detectio
class change accurate identification
analysis change area
ect level differencing normalization or
sformation image
87
4.2.1.1 Pre-classification
In image differencing techniques, two images at different dates are spatially registered
and subtracted one from the other. The corresponding pixel values are subtracted to
produce a new image that represents the change between the two dates.
Mathematicall¡ this method produces a new change image with the same number of
bands as the input images. Pixels presenting a signiflrcant radiance change can be
expected to lie in the tails of the distributions of the different image, whereas the
remaining pixels should be grouped about the mean. Images used in this technique
Image differencing is simple to implement and has been widely used in a variety of
applications. Brera and Shahrokhi (197S) applied this technique for the analysis of
desertification in the Sahara region, and Vogelmann (19SS) investigated the detection
of change in temperate forests. Criticism of this method comss from Riordan (1980 in
existence of mixed pixels, and to radiometric differences between the input images.
This method fails to speciff the mechanism of types of information loss (Singh, 1989:
p. 993). For example a forest with a radiance value of 170 in band 4 on one band and
140 on the second date shows a change of 30. This nurnber could come from the
difference between radiance values of 70 and 40 which does not represent the radiance
value of the forest pixel. Thus, there is a potential loss of information with the use of
simple image differencing technique. The other disadvantage is deciding the threshold
boundaries between change and no-change pixels displayed in the histogram (Jansen,
88
1986: p.2a\. However, the normalized or standardized image differencing of two or
more images from different dates that have been radiometrically corrected have the
potential to be used for digital change detection. In a resent stud¡ Fung (2000) was
using image differencing. Image difÊerencing was also determined to provide the best
'Wiregrass
results for detecting change in Long life Pine ecosystems (Houhouis and
Micher, 1996).
4.2.1.1.2 lmageRatioing
This method calculates the ratio of the values of corresponding pixels of registered
images at different dates. The basic concept of this approach is that if the intensity of
reflected energy is nearly the same in each image then this indicates no change. If
change occurred in a particular pixel, the ratio is expected to be either higher or lower
than the ratio in the no-change areas. Similar to other techniques, the standardization of
This technique has not been as intensively investigated as image diflerencing (Nelson,
1983: p. 130a). Todd (1977) used this method to determine urban change in Atlanta,
Georgia. As a result 91.4% of land cover change was correctly identified (Todd,1977
in Nelsoru 1983: p. 1304). Fung (1988) tested the accuracy of this method by
comparing it with image differencing and principal components methods and found it
of over 9}Yo,but it was difhcult to detect change from cropland to residential (Fung,
1988: p. 14s3).
89
Howarth (1931 : p.287) found that the use ofband ratios using data from different dates
produced an enhanced image of the study area. He used this technique successfully to
detect flooding area in Mamawi Lake, Canada. Rachman (1997) concluded that the
band ratioing technique is suitable and efficient for the study of detecting altered rocks,
Ratios of band 3 and band 7 of Landsat TM for example, is best suited for recognizing
roads or other man-made features that appear in a lighter tone due to their relatively
high reflectance in the red band (band 3) and low reflectance in the mid-infrared
(band7) (Lillesand, 2000: p.52Ð. With a strong local knowledge and correct selection
of the band ratio, the change detection using this method can be successfully
conducted. This technique is well established and is stilt being used for purposes such
(Watson et a1.,2001).
4.2.1.1.3 VegetationlndexDifferencing
This technique also has been established for a long time and is the most common
method used to identify vegetation. In the vegetation index differencing technique, the
basis, or in combinations of two or more bands. The most commonly used band
combination techniques involve the calculation of vegetation indices. These have been
developed on the basis of the strong absorption of red and strong reflectance of near-IR
band ratios, or a combination of both. It has been found that the ratio of near-IR to red
90
reflectance is significantly correlated with green leaf biomass (Tucker, 1979). There
are many vegetation indices in use (Qi, 1993; Elvidge and Chen, 1995, Salerr¡ 1995)
but the most common is the Normalized Difference Vegetation Index (NDVI)
(Eastmar¡ 2001), given by the ratio of near-IR minus red reflectance and near-IR plus
red reflectance.
It is important to note that ratioing two spectrat bands negates the effect of any
extraneous multiplicative factors in sensor datatl:øit act equally in all bands (Lillesand
and Kiefer, tggT). Singh (1939) underlined that the band ratio technique might
enhance noise that is not correlated in diflerent bands. Qi et al.(1993) stated that it is
useful especially when satellite data are available at high temporal frequencies, wffi
fine spatial resolutior¡ and with atmosphere correction of data, they advocate a
correction to account for first-order soil background effects when using vegetation
indices. In a study on the use of vegetation indices for change detectior¡ Lyon et al.
(tggg) found that the Normalised Difflerence Vegetation Index (NDVI) produced
results superior to those of other vegetation indices. Nelson (1983) examined this
indicated that of three methods tested (i.e. image diflerencing, ratioing, and vegetation
index difterence) the vegetation index difference most accurately delineated forest
canopy change. Angelici et at. (1977) employed the difFerence in the ratio between
NIR and red band and the thresholding technique to delineate the change area. The
existence of mixed pixels. However, this technique could be used on land use/land
cover change detection combined with other techniques (i.e. principal component
91
techniques) to delineate other vegetation signatures. Recently, these techniques have
been intensively used to identiff vegetation and change detection (V/ard et al., 2000;
Fung and Siu, 2000; Lanjeri et a1.,2001). Fung and Siu (2000) conclude that NDVI
values were related to woodland, tall scrubland and high-density urban areas. They
also found the decreasing NDVI values covered a large part of rural area revealing
urban expansion. Ward (2000) found the NDVI band 5 and band 3 can successfully be
Byren (1980) stated that Principal Components Analysis (PCA) is a powerful method
orthogonal coordinate system such that the data can be represented without correlation.
Therefore the degree of correlation between bands in the original multi-spectral image
Principal Components Analysis (PCA) can also be applied to compress image datasets
from two or more dates i.e. to multi-temporal image data (e.g. Byrne et al., 1980: p.
175; Richards, 1984: p. 35; Ingebritsen and Lyon, 1985: p. 687; Fung and LeDrew,
92
1987; Ceballos and Bottino, 1997; Li and Yeh, 1998; Edge, 2001). There is a high
correlation between image data for regions that have not changed significantly and a
relatively low correlation between regions that have changed substantially. Richards
localized change, associated with bushfire damage and vegetation re-growth following
fire burns. Singh and Harrison (1985: p. 895) compared the use of a covariance matrix
and a correlation matrix. They found that there is a significant improvement in signal-
to-noise ratio and image enhancement by employing the correlation matrix rather than
justified due possibly to the differences in radiometric resolution between the spectral
bands of sensor. Fung and LeDrew (1987: p.1657) also examined image enhancement
stated that principal components analysis is a "scene dependent" technique and the
exact nature of the principal components derived is not known without an examination
of the eigenstructure and visual inspection of the image. However, it is a powerful data
characteristics of the study area to avoid drawing any faulty conclusions. Ceballos and
landscapes in Northeast Brazil. This approach also has been used to discriminate
healthy and non-healthy corals in Fiji (Holden and LeDrew, 1998). Li and Yeh (1998)
change in monitoring rapid land use/land cover change and urban expansion in the
93
programming methods was conducted to calculate fraction botmds of main seaweeds
and inter-tidal components for mixed pixels (Bajjouk et al., 1998). A recent study that
has been successfully conducted is discrimination of urban and non-urban area using
principal component analysis (Edge, 2001). This study has successfully discriminated
'West
urban and non-trban areas as a basis to predict urban growth in part of Java.
Change vector analysis is the analysis of change based on the vector difference
between the multi-band digital vectors of the pixel on two dates. A particular pixel in
an image can be plotted as a point in this vector space with co-ordinates that
correspond to its brightness values in the appropriate spectral components. The data
values associated with each pixel thus define a vector in multi-dimensional space. If a
pixel undergoes a change from time tl to time t2, avector describing the change can be
defined by the subtraction of the vector at t1 from the vector at t2. This is called the
spectral change vector. It may be calculated from either the original or transformed
(e.g. PCA or NDVI) data, and using either individual pixel or clusters formed by a
spectral change vector exceeds some specified threshold criterior¡ it may be concluded
that change has occurred. The direction of the vector contains information about the
type of change. Malila (1980) introduced this method for the fnst time for detecting
forest change in northern Idaho (Yuan, 1999: p.25).Colwell (1981) used this technique
to detect forest change in the Palouse District of the Clearwater National Forest in
94
Idaho. He used two different strategies of digital change detection: delta-classification
that using this two-change detection procedure could significantly affect the estimates
ofamount ofchange.
Lambin and Strahler (I99$ used this method to detect and categorue land cover
change processes. They concluded that change vector analysis combined with a PCA
change between time trajectories of NDVI data. Cohen et al. (1999) used this
Northwest region. They conclude that a reference image was important for accurate
was compated with image differencing techniques, the image differencing technique
performed better than change vector analysis (Cohen, 1999: p. 99). The disadvantages
of this method are similar to those of other methods such as the raw image differencing
4.2.1.1.6 DirectMulti-DateClassification
more datasets to identifr areas of change. The combined data set whether, an original
approach. In the supervised approacl¡ training sets pertaining to change and no-change
95
areas are used to derive statistics to defure sub-spaces of the feature (normally spectral)
analysis, and subsequent inspection can reveal where changes have occurred. Hoffer
and Lee (1939) used this method for forest change detection and combined it with PCA
to reduce the number of bands. They were able to detect at an accuracy of 90.6Yo for
the fuIl 12 band, and 90.4Yo using 6-band after reduction by principal components
analysis. Muchoney and Haack (1994) conduct a study of change detection for
gnder the similar condition, change areas were significantly different. A disadvantage
of this method is that it needs accurate identification of the training areas for each
change class. Other disadvantages are: the data need to be re-processed whenever a
new image is added, it is difficult to label the change class, and there is limited
information on "from-to" changes class. However, it is a simple technique because it
post - classification comparison is the most common method of detecting change' This
comparison of the corresponding pixel signature to identifu areas where change has
unsupervised approach. It reveals that change has occurred, and also reveals the precise
nature ofthe change, but requires an accurate rectification and classification processes.
96
Joyce et al. (1980) concluded that post - classification comparison "appears suitable for
detecting land cover change with Landsat MSS data in sites where large areas of
forestland are being converted to cropland". Fisher and Pathiram (1993) found that
area studied had undergone change, whereas it was known that little change had
actually occurred.
classification of aircraft MSS data for wetland change detection, but achieved limited
success. Martin (1939) concluded post - classification comparison gave better results
detection in the rural-urban fringe. Jensen et al. (1995) reported on a later study also
using post - classification comparison for wetland change detection. The result was that
cattail and cattail mixture (sawgrass) classes could be detected. This technique is
commonly used and combined with image enhancement processes to get better feature
signatures in order to classiff land use land cover change in what are called hybrid
techniques (e.g. Deep, 1999; Ward, 2000 and Madhavar¡ 2000). The disadvantages of
the method are: it needs ground truth information to get precise training areas, it is
classif,rcations. The advantages are: it reveals the area that has been changed, it
provides "from-to" change class information and bypasses the difficulties associated
97
4.2.1.3 Hybrid
classification based methods in various ways to minimize effors in land cover change
example, can be used in land use/land cover change detection. Radiometrical image
differencing may be used to identifr the area of significant spectral change, and then a
post-classification comparison applied within areas where spectral change was detected
Pilon et at. (1988) developed this method further and applied it to the detection of
using new approaches to classification such as fiizzy logic and neural networks. These
techniques can be used to classiff land use/land cover categories on digital change
detection techniques . Fvz-zy logic was successful in recognizing thresholds of land use
by using the fivzy membership principle. ANN, as well as the frrzzy membership
approacl¡ was successfully used to classiff land use categories when the spectral
Berberoglu et al. (2000) used hybrid methods of land cover mapping in the
Mediterranean basin with the integration of spectral and textural information using
neural networks. The ANN improved the accuracy of classification results compared
98
with the maximum likelihood approach. Berberoglu concluded that when using
spectral data alone and using spectral and textural data, the accuracy of ANN was
greater than that of maximum likelihood. Ward et al. (2000) used hybrid methods as a
spectral change detection to separate impervious surfaces and exposed soil by using
classification. The result improved the accuracy from 75 % to 83 Yo. Madhavan et al.
(2001) detected spatial growth of the Bangkok Metropolitan area using post-
model. The vegetation-impervious - soil model has been used to distinguish urban from
non-urban morphology.
An advantage of hybrid methods is overcoming the mixed pixel problem and providing
The focus on change detection in this study was to provide an approach that could be
used to detect and identifr land use/land cover change in the study area. Considerations
99
l. The study area is in the Tropics area that has varied physical environmental
conditions such as climate, soil vegetation type and weather as well as growing
cycle.
an original band combination constrained to get correct land use/land cover due
new transformed images and can be used to replace the original band image to
change in the study area. This change detection was based on Maximum Likelihood
and NDVI.
1980); 2) it has been successfully employed in land cover change analysis (Fung and
LeDrew, 1987; Holden and LeDrew, 1998; Ceballos and Bottino, 1997; Lie and Yeh
1998); 3) it has been successfully demonstrated to separate crop types cover without
performing data pre-processing (Jiaju, 1988) and 4) it has been successfully used to
100
NDVI transformation was the most coÍrmon vegetation index that was successfully
4.3 Methods
The main method of digital change detection using remote sensing in this study is
image classification were used to identifu land use/land cover categories. Change
sections describe the integration of remote sensing and GIS operations to detect land
101
Multitemporal
geocorrected
TM
TM Composite
432
NDVI
Training
Classification Area
I 2 3 4 5 6 7 8 9
Check
No
Yes
LU-LC
Classification
Map
Change dstection,
LU-LC identifi cation and analysis
GIS Database
Note:
1 = Forest 4 = Open land 7 = Dry land
2 = Plantation 5 = Settlement 8 = Fishpond
3 = Rice fîeld 6 = Rice field unplanted 9 = Water/reservoir
Figure 4.1. The integrated model of detectior¡ identification and analysis of land
use/land cover change.
r02
4.3.1 Digital change detection techniques
cofltmon image processing procedures to detect and identifr land use/land cover from
The satellite image could not be used directly to interpret or delineate physical
introduced during image acquisition. The ideal or perfect remote sensing system has
yet to be developed (Jenser¡ 1936). Remote sensing devices have constraints such as
pre-process the remotely sensed data before analysing it in order to remove some of
these errors.
during the imaging process (Moik, 1930) will produce a corrected image that is as
characteristics of the original scene. Radiometric and geometric errors are the most
Most digital satellite images have already had systematic error removed. Landsat
satellite data for example, which are acquired through the EOSAT Company, have
103
consists of spatial enhancement, radiometric enhancement and spectral enhancement.
and edge enhancement. Table 4.4 below shows the Image enhancement types and their
function:
These two transformation techniques enable one to enhance and provide new set
images, which have strong spectral signatures of built-up areas or settlement and
vegetation cover. Using the image-processing module in ldrisi32 that provides PCA
andNDVI transformation algorithms, a new set of PCA and NDVI transformed images
were created
105
4.3.1.2.2 Filtering
The result of image classification produces some noise due to a scattering of individual
or small groups of pixel. Filtering processes can be used to enhance the spatial features
of the image. A filter creates a new image by calculating new values using a
mathematical operation on the original cell value and its neighbours. The mode or
focal majority filter, which assigns the most coflrmon value to the centre pixel, is
commonly used to remove very small areas from a qualitative image. All of the images
classified in this study were enhancing by mode filtering to remove isolated pixels.
Multi-spectral classification is the method most often used to extract the information
from remotely sensed data. This procedure assumes that imagery of a specific
that the images have good registration (Jenser¡ 1996). Therefore the physical features
of the earth, in particularly land use/land cover types, could be extracted through the
The supervised maximum likelihood classifier was selected in this study, which has
been commonly used in many studies. Three to five training areas for each land
use/land cover (LULC) category in this study were selected based on the object feature
106
4.3.3 Land use/land cover change detection
This study used a post-classifrcation comparison approachto detect land use/land cover
changes in the study area. This approached is based on Maximum Likelihood classifier
of the PCA and NDVI transformed images. Figure 4.2 below shows the algorithm of
this approach.
1989
t *
Post
Max likelihood
classifier
+ Classification
Comparison
1997
AnalYsis
Original band Transformation classification
The first step is to transform the original image by using principal component (PCA)
and normalized different vegetation index (NDVI). The second step is to combine
these two sets of images and then classiff using supervised Maximum Likelihood
classification techniques. Finally, these two image representing difFerent dates were
compared with each other to detect and identiff any land use/land cover change. GIS
107
techniques, including overlaying and cross tabulatiorU were used to calculate the area
The Kappa index of agreement (KIA) has been selected to assess the accuracy of the
classification result. This method compares two images; one image contains the
interpreted land cover map while the second image contains the result of ground truth
classifications derived from remotely sensed data (Eastmar¡ 2001). This operation
creates an e¡¡or matrix that tabulates the different land cover classes to which ground
truth cells have been assigned. Output also includes column and row marginal totals,
effors of omission and commissior¡ an overall effor measure, confidence intervals for
that figure, and a Kappa Index of Agreement (KIA), both for all classes and on a per
category basis.
Some physical characteristic of the study area have been described and dispþed in the
previous chapter. Physiographically, the study area is flat and consists of alluvial plain
deposits. Most of the land use/land cover conditions in the study area have been
modifred through human settlement and agriculture practices. Population pressure has
108
resulted in the expansion of the agriculture area, especially for wetland rice fields as
well as settlements.
As a result of initial image analysis using simple original band combination (FCC TM
band 4 red, TM band 3 green and TM band 2 blue), the initial identification of land
cover was recognued. There were significant land cover changes as shown in
Figure 4.3a
TMComposite 432
1989
l.Kabupaten Bogor
2.Kabupaten Bekasi
3.Kabupaten Karawang
4.Kabupaten Subang
S.Kabupaten Purw akarata
r09
Figure 4.3b
TMComposite 432
1990
LKabupaten Bogor
2.Kabupaten Bekasi
3.Kabupaten Karawang
4.Kabupaten Subang
5. Kabupaten Purw akanla
Figure 4.3c
TMComposite 432
1995
l.Kabupaten Bogor
2.Kabupaten Bekasi
3.Kabupaten Karawang
4.Kabupaten Subang
5 .Kabupaten Purwakarata
110
Figure 4.3d
TMComposite 432
1997
l.Kabupaten Bogor
2.Kabupaten Bekasi
3.Kabupaten Karawang
Subang
5 .Kabupaten Purwakarata
Four images that have been taken for the same season (dry season): June 1989 Fig.
4.3a, June 1990 Fig. 4.3b, August 1995 Fig. 4.3c and June 1997 Fig.4.3d show the
different land use/land cover conditions in the study area. Red represents vegetative
cover (mostly rice field arca); light blue represents open or arable land and dark blue
represents wetland and water. As seen on those images, the light blue area in 1990 is
slightly larger than in 1989, indicating that open and arable land has increased.
Vegetation is the dominant cover in this area in both 1989 and 1990 (see in kabupaten
Karawang no. 3 and kabupaten Bekasi no. 2). This is very different when compared
with conditions in years 1995 and1997.In 1995 arable and open land (light blue area)
was larger than in L997. The other interesting feature is the vegetation (red area). This
area shifts between the north and south parts of Karawang no 3. In 1995 (Figure 4.3c)
vegetative area w¿Ìs located in the north part, while n 1997 (Figure 4.3d) vegetative
111
atea was located in the south part of Karawang. If this area is rice field, it easy to
understand the shifting of this area as it could be related to the local cycle of growing
From this initial analysis it can be seen that land cover in the study area is dynamic
and changes over time. However the investigation of this change is still at a broad
category level, for example, settlement areas are still not clearly recognised.
Table 4.5 shows the eigenvector matrix as well as the loading matrix from principal
component transformation for 1989 (see Appendix Chapter IV for other years). The
total eigenvalue indicates that the first principal component generally contains 85 - 95
per cent variance. This means the first component contains the majority of the
information. On the other hand, the loading matrix expresses the degree of correlation
between each component and the original bands. The highest correlation is between
tt2
Table 4.5: Principal components analysis transformation TM 1989
The strong correlation in this table shows that component 1 can represent the entire
original band and that the most representative band is band 4. The eigenvector value
within the matrix show that principal components 4 and 5 have a consistently high
value and correlate with eigenvector 3 and 6 for whole year, which is related to
settlement and other features that are sensitive to soil moisture content. Others
principal components did not demonstrate a consistently high value except principal
is related to coastal sedimentation and this feature is not included in land use/land
and PC 2) contain of the majority of information in the image. Therefore, all principal
113
components, except principal component 6, were to be combined with the NDVI image
Io anaLyze land use/land cover in the study area. Figure 4.5 shows spectral features of
principal component 5 (4.4 a) and spectral features ofNDVI (4.4b) from TM 1989'
4.4a
The dark black colour in Figure. 4.4a represents settlement features and the light
white
colours on Figure. 4.4b rcpresent forest and other vegetation feattrres. Through the
combination of pCA and NDVI transformed images, the settlement and vegetation
tt4
4.4b
The spectral signature information from the new two transformed images indicates that
land use/land cover types have a strong feature signature. Nine categories of land
use/land cover in the study area have been classified by using a supervised maximum
likelihood classification approach. Table 4.6 shows the description of the nine
115
Table 4.6The Land use/land cover Classification Scheme
The results of final classification fromthe image 1989 to 1997 after filtering are shown
Figure 4.5a. The final classifred image of land use/land cover 1989
116
Land use-.land cover 1990
Foresl
Plantatìon
Rice Reld
Oponrberc lend
Settlement
Rice unFlanted
Ä
Dry land 10 0 10 20 Kilometêrs
Fish Þond
Water/RêseNoir
Figure 4.5b. The final classif,red image of land use/land cover 1990
Figure 4.5c. The final classiflred image of land use/land cover 1991
tt7
Land u6eland cover 1992
Fùrosl
Plantetion
Rìce field
Openibare land
Settlemont
Ricâ unplanted
À
Dry lsnd 10 0 10 20 Ktlomqters
Flsh pqnd
Water/RosBNoir
Figure 4.5d. The final classified image of land use/land coYer 1992
Figure 4.5e. The final classified image of land use/land cover 1993
118
Land use-land çover 1995
Forest
PlåntÉtion
Rìce field
Openibare land
Settlsment
Rice unplántêd
A
Dry land 10 0 10 20 Kilometers
Fish pond
Watar/RåseNo¡r
Figure 4.5f. The final classified image of land use/land cover 1995
Figure 4.5g. The final classified image of land use/land covet 1997
119
Table 4.7 shows the accuracy assessment of years 1989 and 1997. Overall accuracy
was 90 %o for 1989 and 89Yo for 1997. From this table it can be seen that commission
enor (misidentiffing no change and changed pixels) in 1989 was an average of 0.15
and omission (misidentifying change and no change) error was an average of 0.06. In
1997 commission error was 0.20 and omission error was 0.08. The result was an
overall accuracy of 90Vo and 89%o, which is good when compared with 85% as a
standard minimum accuracy of land use/land cover classification from satellite data
(Janser¡ 1996).
t20
4.4.4.The land use/land cover change detection
Once the classifications have been generated, it is possible to carry out a theme
comparison between two classifications. This operation is performed by cross
Figure 4.6 shows the annual fluctuation of land use/land cover in the study area- From
this figure it can be seen that, in general all land use/land cover categories fluctuate
annually. Some categories increased and others are decreased. The land use/land cover
category that fluctuated the most was dry land. Settlement constantly increased in the
five consecutive years (1989-1993) and increased significantly between 1995 and
1997. On the other hand, rice fietd was tending to decrease and rice field unplanted
increased but decreased significantly from 1995 to 1997; fishpond and water/reservoir
remain constant.
300000
250000 *Forest
*Plantataion
200000
G -#Ricefield
QPg¡ l¿¡d
IU 150000
o -..#Settlement
100000 -Ë{- Rice unplanted
#Dryland/Grass
50000 *FFish pond
r
Lî: -l[-Water/Reservoir
0
1989 1990 1991 1992 1993 1995 1æ7
Year
Figure 4.6. The arurual land use/land cover change in the study area
r2t
More detailed analysis of land use/land cover change detection in annual 5-year
consecutive (1939 to1993), 2-year time interval between 1993-1995 and 1995-1997,
and 8-year time interval between 1989 and 1997 n the study area ate explained and
discussed as follow.
The cross tabulation of annual land use/land cover classiflrcation generates a land
use/land cover change matrix (in percentage and hectares) within the study area
(Appendix Chapter IV). Tables 4.8a to 4.8d show the summary of major land use/land
Table 4.8 a: The summary major changes in land use/land cover category
1989 I 1990 (dryldry season)
From Table 4.8a it can be seen which major land use/land cover category changed
between 1989 and 1990. Open land increased while dry land decreased. The increase
of open land was the result of conversion from rice field (planted), and decrease of dry
land and plantation was due to conversion to rice field (planted) and rice field
(unplanted). Both 1989 and 1990 images were acquired in the same season (dry season
at the harvesting time-Table 3.7 n Chapter 3). In the time interval 1990 to 1991 the
t22
land use/land covers that changed significantly were rice field and dry land. Rice field
Table 4.8b: The swnmafy major changes in land use/land covef category
r990lt99t
Category Change Conversion shifr
Rice field 630/o deqease (116924 ha) to dry land and rice
held unplanted Within rice field
and rice field (unplanted)
Dry land 247Vo ncrease (78565 ha) from rice field (unplanted)
and rice field (planted)
The decreasing of rice field (planted) and increasing of dry land in this time interval
was due to the different seasons. The 1990 image was acquired in the dry season at
harvesting time and the 1991 image was acquired in the wet season at planting time. In
the 1990 image (harvesting time), rice fields were being planted while in 1991, rice
fields have been harvested and are beginning to grow new plants (planting time).
Therefore the rice field þlanted) area was greater in 1990 than in 1991. Rice field
(planted) in 1990 was converted to dry land and rice field (unplanted) in 1991. This
phenomenon also was found when one analyzes change between 1992 arnd 1993;
which is 1992 was dry season or harvesting time and 1993 was wet season or planting
time. The change pattern in this time interval was the same as in time interval 1990 -
l99l rice field (ptanted) decreased and dry land increased (Table 4.8c). The decrease
of rice field was due to conversion to dry land or rice field (unplanted)'
t23
Table 4.8 c: The summary major changes in land use/land cover category
19921 1993 (dry/wet season)
In the opposite direction of the change between dry (harvesting time) to wet (planting
time), is the change between 1991 and 1992 as show in Table 4.8 d.
Table 4.8 d: The summary major changes in land use/land cover category
199I I 1992 (weldry season)
Rice field (planted) in wet season was less than rice field (planted) in dry season. In the
\Met season of 1991, most of the rice fields were under preparation to grow new plants
and many rice fields were appearing as dry land or rice field (unplanted). Therefore in
this time period (1992-1992), rice fields (planted) were seen larger than dry land or rice
freld (unplanted). The increasing and decreasing of rice field (planted) is related to
decreasing and increasing of dry land and rice field (unplanted). This is the cycle of
change among these categories between different growing seasons (harvesting and
planting times).
124
The other interesting feature was the change in forest category. In the dry season,
forest area was smaller than in the wet season (see Table 4.8c and 4.8d). This
phenomena is related to the phenomenon of the "leaf on and leaf off' of broad leaf
reflectance (Yuan et al., 1999). In the dry season lots of leaves fall from broad leaf
plants, so at this time the reflectance of leaf decreases, while in the wet season the
reflectance increased from the leaves that have re-grown. Therefore decreasing and
effect of the dry and wet seasons and could be a classification error due to a mix
4.4.4.2 Land use/land cover change in 1993 to 1995 and 1995 to 1997
Wet and dry season phenomena of land use/land cover changes can be continually seen
in the 2-year time interval (Table 4.9a). From Table 4.10U it can be seen that forest
and dry land were decreasing and rice field (planted) was increasing; this is similar to
the patterns of change from wet to dry season between l99l and 1992.
Table 4.9a: The summaly major changes in land use/land cover category
19931t995
Category Change Conversion shifr
Forest 4IVo de,crease (8933 ha) to plantation
Rice field 22%o increase (21209 ha) from rice field (unplanted) Within rice field
and dry land (unplanted)
Dry land 75%o deqease (58591 ha) to rice field (unplanted)
and rice field (planted)
Another 2-year interval detection is in the dry-to-dry season between 1995 and 1997.
In this time interval, the change pattern was different from the annual detection. In this
time interval settlement increased, while rice field (unplanted) decreased and dry land
t25
increased. The increasing of settlement areas related to decreasing of rice field
Table 4.9b: The summary major changes in land use/land cover category
1995 I 1997 (dryldry season)
In dry-to-dry season between 1989 and 1990 settlement did not increase, while in dry-
to-dry between 1995 and 1997 settlement was significantly increased. It means since
Table 4.11 shows a summary of land use/land cover change between 1989'1997 (the 8-
year interval). The change of land use/land cover in this time interval should represent
the real conversion or change ofthe land use/land cover in the study area after 8 years.
t26
Table 4.10: The suÍìmary major changes in land use/land cover category
1989 I t997 (dryldry season)
Within the 8-year interval there are many land use/land cover categories that changed
(Table 4.10). Open land and rice field þlanted) decreased and settlement increased,
while plantation decreased and forest increased. In this time interval it was also
settlement in this time interval related to the conversion from rice field (unplanted),
open land and dry land, while decreasing of rice freld þlanted) related to conversion to
rice field (unplanted). From this change trajectory, the sequence of land use/land cover
change in this aÍea caîbe recognized. Rice freld (planted) was converted to rice field
(unplanted) and open land, then rice field (unplanted) and open land was converted to
settlement.
There is no change or conversion from forest to rice field or settlement. Forest change
was related to plantations: if plantation decreases then forest will increase (confuse
classification). The mix between forest and plantation was recognized in the field when
ground truthing was conducted. In the study area, there were some broad leaf
plantation areas such as rubber and teak wood plantations that 3 or 5 year later will
t27
appear as forest on the image. These plantations were located in kabupaten Karawang
The change detection of nine categories of land use/land cover in the study area has
resulted in the various shifts and conversions among categories, especially in annual
time intervals. Change detection on the 8-year time interval shows no annual
variability except for rice field (planted) and rice field (unplanted). Rice field (planted)
and rice freld (unptanted) fluctuate due to different growing cycle. V/ith this condition
it was difficult to recognize the permanent change of rice field as well as open land
aÍeas. It was recognized in the freld that rice field þlanted) and rice held (unplanted)
have the same function as rice field. Open land and dry land categories are the same
as
land which is unused, or under preparation for construction or for dry land agriculture.
to recognize the change of rice field and open/dry land within the 8-year time interval.
Figure 4.7 shows land use/land cover after reclassificationinto seven categories. Rice
field (planted) and rice field (unplanted) were regrouped as rice fields. open land and
dry land were regrouped into open/dry land. From this figure it can be seen that rice
field area is the largest and was the dominant land cover in whole of the study area.
Settlement was seen to extend along the road corridor from the west to the east in the
middle of the study area, and open/dry land was located within rice field and settlement
area.
r28
Lond u¡r 1989
Foro¡l
Phnlrllon
Rlcr fl.ld
Òprn/dry land
,S¡ül¡nirnl
flrh pond 10 0 10 20 Kilomêters
w6ler/feso]volr
t29
Table 4.ll:Landuse/land cover in the study arcayear 1989-1997 (ha)
Water/reservoir sso2l
l.oe s4681 1.09 441 0
350000
A
300000
250000
tú
.C 200000
(ú
E 1 50000
1 00000
50000 tr 1989
E 1997
0
20000
B
1 5000 E Forest
1 0000
B Plantation
o E Rice field
! 5000
o g Open/dry land
c) 0
-5000
I Settlement
¿ Fishpond
-1 0000
E Water/resenoir
-1 5000
Land use change
Figure 4.8. The land use/land covers change 1989-1997 in study area (a) and Total
change of each category @).
130
From Table 4.11 and Figtre 4.8a and 4.8b. above, the change of land use/land cover
among seven categories between 1989 and 1997 can be seen. Rice field decreased as
well as plantatior¡ while settlement, open/dry land, and forest increased. The decrease
in rice field was related to increasing settlement and open/dry land. It is clearly shown
that there was conversion from rice field into settlement. Settlement increased by 8l %
or 18262 ha, while rice field decreased by 14866 ha or 5Yo. Plantation decreased, forest
increased, and water /reservoir and fishpond was relatively constant. As has been
mentioned, the increase in forest in the study area was due to the effect of "leaf on-leaf
off in the different seasons and confirsion between forest and broad leaf plantation.
4.5 Discussion
The focus on change detection in this study is to recognise land use/land cover change
There are three common digital change detection strategies: pre-classification, post -
classification and the hybrid approach. Pre-classification is the method that directly
differentiates between two dates' spectral signatures. This approach was based on the
assumption that different spectral signature over time represent a change in land
surface conditions. post-classification is based on images from two different dates that
131
From the initial recognition based on the false colour composite analysis of the study
area, it was recognised that land cover in the study area is dynamic and changes over
time. However, this initial recognition did not give detailed information of land cover
category due to limitations of spectral signatures on the original images of the false
colour composite. The PCA and NDVI transformation was used to enhance the spectral
were combined and used for a post-classification comparison approach to detect more
This study used post-classification comparison of two transformed (PCA and NDVI)
images for the following reasons: 1) the study area has varied physical environmental
conditions due to the weather/climate conditions; 2) the study area has varied land
cover types; 3) the land use in study area has extensively changed due to the population
and economic growth close to the Jakarta metropolitan atea; 4) PCA and NDVI
combination is an effrcient method and has potential for data processing and
This approach was successful in detecting land use/land cover change in the study area.
It was found that weather or season and growing cycle strongly influenced the land
use/land cover conditions in the study area. In the classified image 1990 (acquired on
June-the early dry season) forest cover was less than in 1991 (acquired on May - the
late wet season) or within 1992 (acquired on July - dry season) and 1993 (acquired on
Aprit - wet season). Dry and wet seasons influenced the "leaf on-leaf off'
t32
phenomenon on the classified image. This phenomenon was clearly recognised in the
change of plantation and forest categories. In the dry seasor¡ forest cover was less than
in the wet season due to lots of leaves fall from plants ("leaf off'), while in wet season
lots of leaves have re-grown ("leaf on"). Therefore in wet season, forest cover
appeared to be increasing because of a lot of "leaf onl' and also mixing with increasing
"leaf on" of some broad leaf plantation such as rubber and teak wood The decreasing
and increasing forest and plantation in this study could also be an error in classifrcation
due to mix between forest and broad leaf plantation such as ruber and teak wood.
Rice field (planted) and rice freld (unplanted) were clearly recognised and changed
according to growing cycles. In the harvesting time (the month of June - JulÐ, a lot of
rice fields were covered by rice plantation that was ready to be harvested; therefore, in
this harvesting time the identified rice field area was large. In the planting time (the
month of April - MÐ a lot of rice fields were empty, because the land was being
prepared to grow new rice. In this planting time the rice field area \ilas less than in the
month of June-July (harvest time). This phenomenon was clearly recognised in the
From the summaries of annual change detection from 1989 to t993, it can be
recognised that for 1990 to 1991 the area rice fields (planted) was larger than rice field
(unplanted). The 1990 image was acquired in June in the early dry se¿rson or harvesting
time, while the 1991 image was acquired in May in the wet season at the planting time.
Therefore, the changes of land use/land cover within this time interval dry-to-wet
season were seen as conversion from rice field þlanted) to rice field (unplanted) or dry
133
land. On the other hand, change within the time interval 1991 - 1992 wet-to-dry season
was seen as conversion from rice field (unplanted) to rice field (planted). The 1992
image was acquired in July atthe harvesting time, when rice field areas were covered
For the long-term time interval detection (the 8-year 1989-1997) as well as the 2-year
time intervals 1993 - 1995 and 1995 - 1997, seasonal variability was not an influence
on land use/land cover change. In this time interval, especially in the 8- year time
interval, it was recognized that some land use/land cover categories have permanently
changed. Settlement was significantly increased, but rice field (planted) and rice field
(unplanted) were still seen to fluctuate due to different growing cycles. Therefore,
regrouping of rice field (unplanted) and rice field (planted) into rice field as well as
open land and dry land into open/dry land, was necessary and resulted in more clearly
significantly increased as a result of the conversion from rice field. Forest increased
4.6 Conclusion
Land use/land cover was clearly differentiated and could be classified into land
Open land; 5) Settlement; 6) Rice field unplanted;7) Dry land; 8) Fishpond and 9)
Water/reservoir.
t34
The trajectory of land use/land cover change in this area \¡ras recognized as follows:
Rice field was converted to rice field (unplanted) and to open land, then rice field
(unplanted) and open land was convert to settlement. There was no conversion from
The change of forest in this study area did not represent deforestation or reforestation,
it was error in classification due to mixes between broad leaf plantation and forest that
related to dry and wet season ("leaf off- leaf on" phenomena). As found in the field, a
lot of some broad leaf plantation (such as rubber plantation and teak wood found
especially in the south part of Karawang near to the northern pafi of Jatiluhur
It was necessary, for long-term detection change analysis, to reclassifu or regtoup the
land use/land cover categories in this study area. This regrouping was important in
order to recognise which land use/land cover had permanently changed. Rice field
(planted) and rice field (unplanted) have the same function in the field as rice field
area, and open land and dry land categories have the same function as a land which is
unused or under preparation for construction or for dry land agricultwe- These
135
After regrouping, the categories that really change are settlement and rice field.
Settlement increased while rice field decreased. There was no change or conversion
136
Chapter Five
5.1 lntroduction
Land cover in the study area changed on an annual basis related
to seasonal variability
detection, it
rice field (planted) the next year. From 8-year or long-term time interval
changed
was recog nzed that the land use/land covers that have permanently
are
settlement and rice field. Settlement increased while rice field decreased' The
temporarily, related
land use/land cover categories such as, open/dry land was changed
related to the
to the weather or seasonal conditions. Forest and plantation was mixed
..leaf on and leaf off' phenomenon aS an effect of weather or Season condition (see
study alea are interesting
chapter 4). Moreover, the land use/land cover changes in this
as well as proximity to
cover in the study area. Population density, proximity to roads
factors that influence
the crty centre can represent social-economic growth as dynamic
and land system could
land use/land cover change; factors such as slope, physiography
137
strong influence on land
possible (dynamic and static) factors are the dominant or have
proximity to urban
and static influence factors such as population density,
as well as physiography
and semi-urban center, proximity to roads, slope,
Tunner, 1994).
138
change does not stem solely from population growth. Several studies argue
that human
(The
population aflects land use change more by human actions than by sheer numbers
The focus of the driver or influence factors analysis in this study was divided
into static
and dynamic factors. Static factors were defined as factors that might never change,
urban, and
factors that can be changed, such as population growtt¡ population density,
affluence
semi urban areas and transport routes, which represent impact of population,
GDp, industrialization and cultural structure, are also factors that can be included as
informatioru
dynamic drivers, but due to the diffrculties and limitation of collecting this
this analysis only includes factors such as population growttU population density'
proximity to the center of urban and semi urban areas and proximity to roads'
between
Discriminant analysis is a technique for simultaneously examining differences
introduced this technique for the first time in 1936 (Fisher, 1936). The concept
Linear combinations
underlying discriminant analysis is fairly simple (Marija, 1990)'
ofthe independent variables are formed and serve as the basis for classiffing cases into
difference will be reflected in the difference in means (for the variable height).
Therefore, the variable height allows us to discriminate between males and females
139
with a better than chance probability: if a person is tall, then they are likely to be male;
at a|,2000;Tu et al., 2003). Recently, this technique has been used for many
purposes
in geography and remote sensing (Sullivan, 1980; M*Üq 1990)' In pattern recognitior¡
(Du ef
linear discriminant analysis was successful in ctassifring a hyperspectral image
discriminant analysis (LCDA), which replaces the Fisher ratio used in linear
forced to separate.
criterion for optimality. The result was that all classes of interest are
et al'
By means of this direction" LCDA can detect and classifu similar targets. Dt¿tn
instead of
the data, such as similarities or the distances to a subset of the object,
140
and allowed the qualification and
useful tool for explaining the vegetation distribution
geology changes.
used the cartographic variables of elevatior¡ aspect and slope to predict forest cover
that a feed-forward artificial neural
type. The result was that the comparison indicated
type than did a traditional
network model more accurately predicted forest cover
t41
Himalayan river birds. They compared the perforrnance of multiple discriminant
analysis (MDA), logistic regression (LR) and artificial neural networks (Allf$. The
result was that MDA was preferable over LR or ANIN, but where there are complex or
advantageous.
5.3 Method
Overlay analysis in GIS was used to analyse the relationships between static as well as
dynamic factors and land use/land cover change in the study area. Static factors such as
slope, physiography and land system have been created as well as dynamic factor
(Table 5.1). The classes of each static factor followed the standard classification system
of slope class, physiography and land system used in Indonesia @irectorate Land use
and RePPProt, 1986), while dynamic factors were created by buffering a 10 km radius
from the centre of urban area and semi-urban area. This radius assumes that l0 km is
the longest distance people walk to work in agricultural areas. Urban and semi-urban
areas were defined based on the status of the city as a capital city.The capital city of a
kabupaten (district) was defined as urban and the capital city of kecamatan (sub-
district) was defined as semi-urban. The total population at kecamatan level in which
The buffer lkm from main roads (the state roads), others roads (the provincial and
kabupaten roads), and toll way the proximity factor. This assumes that I km is the
maximum distance of built up area located along the roads. This bufler area was
r42
overlayed with land use/land cover to identiff the relationship between roads as an
accessibility factor with land use/land cover change in the study area'
Another dynamic factor, population density, has been created based on the total
The area
total population data were collected from the statistical offrce of Kabupaten.
of each kecamatan was calculated based on administrative boundary map that was
to
derived from topographic maps (chapter 3). Table 5.1 shows the method employed
t43
factors
Table 5.1. The Method to create area of possible influence or driver
factor. while buffering of "line" features such as main roads, toll way and other roads
to create ar: areaof proximity to roads also as a possible dynamic
has been conducted
linear regression to
The quantitative analysis undertaken here consists of simple
use/land cover' and the
analyse relationships between populations and land
factors had the most influential on
discriminant analysis approach to determine which
analysis is used to build a predictive
land use/land cover in the study area. Discriminant
or independent
In this stud¡ discriminant analysis was used to analyze the predictor
r45
proximity to the city' to main roads'
Independent variables, such as population density,
Table5.2:Yaiableusedinthediscriminantanalysis
Based on this database, simple linear regression was conducted to identiff relationships
between population and land use/land cover for the years t989 and 1997 '
5.4 Results
5.4.1 Overlay analys¡s
An intersection between land use/land cover and static driver factors resulted in
cover
identification of various relationships among driver factors and land use/land
change in the study area. Tables V.1, Y.2 arrd V.3 (Appendix V) show the matrix of
land systems'
intersection between land use/land cover versus slope, physiography and
use/land cover
From these tables it can be recognized that only fow categories of land
categories such as operVdry land, fisþond and water/reservoir did not change
significantly.
cover change
Table 5.3 shows the summary of the relationship between land use/land
categories of forest,
and slope. From this table it is evident that the land use/land cover
plantation, rice field and settlement significantly changed within all slope classes'
Forest decreased mostly within slope 30-4 5% and > 45yo and settlement increased
plantation decreased within slope g-15%o and ls -30yo, while rice fields decreased
shows that
within slope 0-3 Yobulincreased within slope more ttøî3%. This change
147
(settlement
land use/land cover that represents human intervention or intensive land use
and rice field) is mostly located in the area within slope 0-3%.
Table 5.3: summary of land use/land cover change versus slope (1989-1997)
15-3Ûo/o >458/a
Tvpes\SloPe
-201
Forest 1) 6116 609 110 108 -191
-8.02 -9.16
2) 1.31 3.48 1.31 1.69
-2898 -1325 -1 049 -221 -151
Plantation 1) -5850
-1.25 -16.59 -15.85 -16.45 -9.28 -6.89
2)
345
Rice field 1) -1 9486 1677 1 166 909 434
18.23 15.71
2) -4.16 9.60 13.94 14.25
1 8043 330 -18 16 -4 -3
Settlement 1)
3.85 1.89 -0.22 o.24 -0.15 -0.13
2)
1) = unit in Ha
2)= Percentage of change area to total area of slope c/ass
the intersection
Another relationship that supports this phenomenon is also revealed by
land use/land cover and physiography (Table 5.4). Settlement and rice
fields
of
in rolling
significantly changed within flat areas, while forest and plantation changed
From both the intersection of land use/land cover with slope, as well as with
use/land cover
physiography, it was clearly recognized that, in general, the major land
located in the
in the study area such as settlement and rice field, were predominantly
flat areas within srope 0-3%, fo[owed by prantation and forest in rolling
hill and
148
Table 5.4: Summary of land use/land cover change versus Physiography 1989-1997
= unit in Ha
2)= Percentage of change area to total area of physiograPhY class
use/land
Further information was obtained from the result of intersection between land
cover
cover and land systems. From Table 5.5 it can be recognized that land use/land
such as rice field and settlement significantly changed in areas with land
systems of
forest changed in
coalescent estuary and undulating to rolling sedimentary plains, while
hillocky plains (Table 5.5). Land systems of coalescent estuary and undulating to
rolling sedimentary plains were located in the flat area, and consist of a fertile soil and
plains are located
are rich with minerals from volcanic sedimentary material. Hillocky
149
land system (1989-
Table 5.5: Summary of land use/land cover change ver$$
Land system
Forest 11
t) 4
I 0 0
)" 0 0 0.
Plantation
353 1
t) -0 0
2.6 -7 1.3
2) 1
tso
5.4.1.2 Dynamic drivers
Overlay analyses of land use/land cover change with dynamic driver factors were
conducted through examining the intersection between land use/land cover with
population density, proximity to the center of urban and semi-utban areas, and
Population densitY
Table 5.6 shows the relationship between land use/land cover change and population
density. It is evident that settlement has consistently increased in area with higher
and class III increase d (by aaal. ha (134%) (Figures 5.la and 5.1b). Areas of rice field
decreased with increasing population density. Rice field area decreased by 25277 ha (-
l l%) in areas with lowest population density (class I), while in class II it decreased by
Forest increased in area with population density class I and class II, while plantation
area decreased in low population density (class I) and increased in area with high
population (class II and III) (Table 5.6). The increasing forest and plantation in areas
with high population density could be due to areas of secondary forest, as well
as
plantations which consists of broad leaf plantation such as rubber and teak
wood which
151
Table 5.6: Land use/land cover and Population densþ Classes
Ä A
152
À
^
greatest change
Rice field and settlement land use/land cover categories underwent
also occurred
within the 10 km radius from the urban centers (Table 5.7). This change
areas (Table 5.8). In this
in the area within l0 km radius from the center of semi-urban
as Cikarang. These areas afe about 30 km fiom Jakarta (Table 5'9 and Figure 5'2)'
153
of urban area (Ha)
Table 5.7:Landuse - land cover change withinl0 km
Urban
Bekasi *)
;*J lr,,n,",,qJ
-11225 1756
,.nnr*"JF¡s¡pq!
0. -35.91 5.62 0
2) 0.
*) -1781 -795 1
Karawang 1
*) change -- unit in ha
Distance
% Change of % Change of from
Urban and
Rice field Settlement Jakarta (km)
Sub urban
* -36 29 1f
Bekasi
-19 14 30
Cikarang**
-6 10 3C
Cileunqsi**
-6 5 5C
Karawang*
** -¿ 75
GikamPek
-2 2 7E
Purwakarta*
5.8
Source Summarized fromTable 5.7 and
**
* (Jrban area; semi-urban area'
154
. Clty
[T:ïi-j Jakarta
Larid use 1997
Rice field 10 2tì Klþmders
Sdtlement
30
Clrange from 1989 to 1997
20
0
o
cD
0
(ú
o -10
s I Rice lield
-20
¡Settbnìent
-30
40
Belesi Cikarang CiÞungsi lGrawang Cikanpek R¡rwakarta
155
From Figure 5.2 it can be seen that the percentage increase of settlement and
percentage decrease ofrice freld is greater for areas that are close to Jakarta than for
areas that are further from Jakarta. It can also be recognized that status as a capital city
of kabupaten (urban) and status as a capital city of kecamatan (semi-urban) does not
purwakarta are capital cities of kabupaten (urban areas) but the large and significant
However, Cikarang and Cileungsi, which are the capital cities of kecamatan (semi-
urban areas), increased in their area of settlement (Table 5.9). This phenomenon
supports the view that the distance from the Jakarta metropolitan area strongly
more than 75 km away from Jakarta. This increase of settlement was related to
trade in
increasing industrial development in this area. Industry, manufacturing and
this area and surroundings increased from 4777 :urruts of industry and manufacturing in
19g9, to 6923 units in 1997 and 8354 units of industry and manufacturing in 2000 (he
Accessibility
cover
goods and the movement of people from one place to another. Land use/land
along roads tends to undergo rapid change due to the accessibility of the
surtounding
afeas
156
From Tables 5.10a to 5.10c it can be seen that rice fields and settlements have
undergone change along the 1 km buffer from the roads. The largest change
is along
the toll road. In this area, rice fields decreased by 3913 ha or '23 % and settlement
along
increased 3810 ha or 22o/o. The second largest change is rice field and settlement
the state road, where settlement increased by 3552 ha or 13 % and rice field was
decreased by 2t35 ha or -8.02 Yo. Alongthe state roads it was rècognised that open land
decreasedby -2346 ha or -8.81%. This decrease could be some open/dry land that has
field. During the fieldwork in1999 it was found that inthis area a lot of new
housing
70. Lu¿Ul
Land use
60 293 233 0.88
Forest
2124 2783 659 2.48
Plantation
ts28l 13146 -2t35 -8.02
Rice field
4234 887 -2346 -8.81
Ooen/drv land 1
t57
Road
Tabte 5.10c: Land use/land cover change in 1 km buffer of other
such as settlement'
Four land use/land cover categories that had changed significantly
linear regression
rice field, forest and plantation have been analysed through simple
5'3d)'
with the total population within the study area (Figures 5.3a to
200000
y=76.217x+50250
R2 = 0.3927
o 50 000
o
ê o o
tr taa' O' a
o o
(E o a
E 00 000 o
CL
o
È
o
(É
o
F
50000
o
a
o a
a
l¡
0
0 500 1000 1 500
settlement 1989
Figure 5.3a Linear regression of Total Population and
158
8.00
a y=0.0017x+0.1129
6.00
R2 = 0.1704
s a a a o
4.00 a
5 a
ìo o
Ð
tr
!¿' o
a aa a
o
.g -500 00
a
500.û l ooo.oo 1500.00 2000.00 2500.00
ê
=
o
o-
-+ où
Settlement change I 989'l 997(Ha)
-6.00
Figure 5.3d Linear regression of the Population growth and Settlement change.
From Figgre 5.3 a to 5.3d it can be seen that there is positive relationship between total
..population and total area of settlement'., "population change and settlement change"
as well as "population growth and settlement change". This means that when
population increases settlement will increase, but the correlation of these relationships
was low with R2 below 0.65, and the lowest was correlation between population
growth
160
in Figure
Relationships between population and rice field, forest and plantation shown
and plantation will decrease. This relationship was poor; with most of IÚ below 0'40,
and especially for forest category was very poor 0.0013 it could be due to
4000
o¡' . 3000
y=-338.91x-57.939
às
o 2000 R2 = O.3477
lE oa
o 1000
E o oaa
o
ED o I
tr -b 4 6
o
-4
-reûo a a
a O
.E
=ÉL
-- C'c
-2000
o a
o-
-3000
a
-4000
Rice field change (ha)
400
y=-3.1765x+3.4528
300
R2 = 0.0013
s 200 o
o
G 100 O
s
ìo
o -4 ^ô
-JU -1
-1 00
1 2 3
Ê
o -200 a
G
=e
o
-300
o-
-400
o
-500
Forest change (ha)
change'
Figure 5.3f Linear regression of Population growth and Forest
161
700
dto a
èS O
o 'soo
6 y = -34.824x+ 176
400
.C R2 = 0.1402
ìo 300
aa
o oo a
C
o a
a 100
-!! o
=
o.
o o
È -4
., 2 4 o I
-1 00
o a a
-200
Plantation change (ha)
Discriminant function analysis results in variables that discriminate between six land
use/land cover categories in the study area. Table 5.12 shows the covariance matrices
for separate group. The values in this table gave some idea of how the relationship
between variables changes from group to group. For example, in group 5 (Settlement),
all variables have negative values except variable population density. It means these
variables have negative and positive relation with settlements. Settlement areas are
large in areas where population density increases (positive), and also will be small in
areas when slope decreases, and proximity to urban, proximity to main road, and
proximity to toll way all decrease (negative si$). This relationship makes sense, in that
settlement areas are close to roads, close to urban, have little slope, and are in high
population density areas. Other groups such as forest, plantation and rice field have
different relationships. For example, group I @orest) has a positive relation with
population density and a negative relationship with slope, proximity to road and
proximity to urban. It means that forest areas are large in areas with high population
t62
density and close to road and to urban. Rice fields also have positive relationship with
population density, slope and proximity to urban areas, and a negative relationship with
proximity to main roads and proximity to toll way. Oper/dry land has the same
relationship. This is a single land use/land cover (one year 1997) analysis. Forest,
plantation and rice field as well as settlement in this study area are located in areas with
high population density and near to the roads. Therefore, rice field as well as forest has
a positive relationship with population density and proximity to the roads. Especially
for forest it could be a misinterpretation of this category due to mixes with broad leaf
From this
The initial statistics from the discriminant analysis are shown in Table 5.13.
table it can be seen that the high values, occur only in four functions
with an eigenvalue
>0.1 or percentage of variance >15 vo. This means that only of this four function that
t63
úr Table 5.14a it can be seen that all variables have a different contribution in each
function. For example, the population density has a greater contribution to function 4
than the other, slope has a greater contribution to function 2, as well as proximity to
main road. It means factors such as proximity to urban areas have a strong contribution
in function 4, proximity to toll way has strong contribution in functioî 2, slope has
road in function 3. This relation was also supported by looking at the relationship
between the dependent variable and discriminant variables in the structure matrix
(Table 5.14b), which gives the canonical correlation coefficients. These values are
'When
comparable to factor loading and indicate the substantive nature of the variables.
some dependent variable has high canonical correlations while other have low, then
the
ones with high conelations contribute most to group separation (Bargman, 1970).
164
and toll way
From Table 5.t4b it can be seen that proximity to urban areas' main road
density in function
have a high correlation in function 1, slope in firnction 2, population
and proximity to toll way have a strong contribution' It means these factors have a
overlay analYsis.
<r vuy ,
1 2 3 4 5 o
6.868 6.975 7.447 6.911 9.808 9.941
6.543 3.249 2.888 2.762 2.783 1.852
Slope
0.689 0.689 o.773 0.511 1.41
Proximity to urban 0.853
-0.267 -0.237 -0.306 -0.190 -0.225 -o.240
to main road
0.759 0.656 0.928 0.553 0.675 1.041
to toll
-r7.609 -10.708 -12.093 -10.410 -13.560 -22.673
(Constant)
Fisher's linear discriminant function
function can be
From Table 5.15 the linear equation of discriminant classifrcation
Groupl(Forest):.17.609+6.868*Popden+6.543*Slope+0.853*Prox.towban-0.267
*Prox to main road + 0.759*Prox to toll way
urban-
Group 4 (Oper/dry land): -10.410+6.911*PopderÛ2362*Slope+0'773*Prox'to
O.tgO* Prox to main road + 0.553*Prox to toll way
165
Group 5 (Settlement) : -13.560+9.808*Popden+2J-83lS1ope+0'5ll*Prox'to urban-
0.225*Prox to main r< ad + 0.675*Prox to toll way
grouped cases coÍectly classifred. It means the accruacy of this linear equation of
5.5 Discussion
shed more light on the research freld of land use/land cover changes' both
use/land cover. The land use/land cover types that represent human intervention
<
(settlement, rice field, plantation) were mostly located in the flat areas with slope 3yo'
The dynamic factors such as population density, proximity to the urban areas and
proximity to the road also had a strong relationship with land use/land cover'
population
Settlement increased and rice field decreased significantly with increasing
roads. This land use/land cover distribution conformed to the schema of land use
In urban areas such as Bekasi, Karawang and Purwakart4 decreasing rice field and
such as in
increasing settlement are clearly recognized, as well as in semi-tlrban areas
decreasing rice freld in the study area was related to increasing population density'
increased 23Yo, class II increased 64yo, and class III increased t34%. while rice field
areas decreased -llyo in area with lowest population density (class I), in class II it
decreased -3lyo, and in class III it decreased -63%. This condition indicates that the
decreasing rice
increasing population density has an effect on increasing settlement and
population
field in the study area. There is no indication of strong relation between
in areas with
density on decreasing or increasing forest and plantation. Forest increased
population density classes I and class II, while plantation area decreased in low
population density (class I) and increased in areas with high population (class II and
III). The increasing forest and plantation in areas that had high population density could
Purwakarta
be due to the type of forest in the study area, especially in the kabupaten
167
progfams (Local
and Karawang, where there is secondary forest as part of reforestation
field in the study atea, is the distance from the Jakarta metropolitan area
The
km away from Jakarta, settlement increased the greatest and rice field
decreased'
afeas.
decreasing rice field and increasing settlement' This proximity factor is easy to
effect, areas close to the road will have the greatest chance to be changed into
168
and negative relation with, slope, proximity to urban, proximity to main road, and
proximþ to toll way. It means that open/dry land and settlement will increase in areas
with high population densþ, little slope, close to urban, close to main road, and close
to toll way. This relationship makes sense because most settlements are located and
extended in the flatateaalong the road and near to the urban afeas'
The other interesting feature is the positive relationship between group l(forest) and2
with slope, proximþ to road and proximity to urban. It means large forests in areas
with high population density and close to roads and urban areas. A similar result was
found from the overlay analysis in that forest and plantation areas were mostþ located
in the areas with high population density and close to the main roads, especially in
kabupaten purwakarta and Karawang. Rice fields had a positive relationship with
population density, slope and proximity to urban, and negative relationship with
proximþ to main road and toll way. It represents the real condition that rice fields are
mostly located in areas vvith high population density, near to main road and toll way
and structure matrix tables indicate the contribution of relationship between dependent
variable and discriminant variables. The canonical conelation is the correlation of two
other a
canonical (atent) variables, one representing a set of independent variables, the
set of dependent variables. These values are comparable to factor loading and indicate
the substantive nature of the variables. The high value of canonical correlation
169
density, proximity to main road and proximity to urban are the variables that strongly
contribute to separation of each group of land use/land cover in the study area.
5.6 Gonclusion
This chapter has analysed the influence of land use/land cover change drivers in the
study area. As the study area is mostly flat in the downstream area of the Citarum
watershed, land use/land cover predominantly consists of irrigated rice fields. Static
factors such as slope and physiography have a strong influence on land use/land cover
in the study area. The land use/land cover types that represent human intervention
(settlement, rice field, plantation) were mostly located in the region of little slope or flat
areas. Dynamic factors such as population densit¡ proximity to the urban areas and
proximity to the roads also had a strong relationship with land use/land cover. This
density, proximity to the main roads and proximity to urban were the variables that had
the strongest contribution to discrimination of each group of land use/land cover in the
study area. Distance to the Jakarta metropolitan area strongly influenced the increasing
amount of settlement and the decreasing of rice field in the surrounding areas such as
170
Chapter Six
6.1 lntroduction
information that can be used by planners to identiff areas that require priority attention
models to simulate and e4plore land use/land cover changes. These models are various
and relate the purpose and the scale of studies (Berry, 7996; Clarke, 1997; Verburg e/
al., 1999;Irwin, 2001;Pontius et ø1.,2001; Fischer, 2001 and Gobin" 2002). Pontius e/
al. (2001), for example, developed a model to simulate the location of land use/land
cover change specifrcally for forest distr¡rbance. They present a GlS-based model,
GEOMOD2,whichqualifies factors associated with land use, and simulates the spatial
pattern of land use forward and backward in time. GEOMOD2 was able to classify
forest cover conectly to between 84 and 86 % with Kappa about 0.34 and 0.45. Gobin
et at. (2002) developed a model to predict the probability of private local agricultural
land use. The result was a binary logistic model for estimating probabilities of private
agricultural land use that correctly predicted 95.7% of the 300 sample plots. Significant
progress in spatial modelling of land use/land cover change has been made when
spatial data sets were available from remotely sensed data (Irwin et al., 2001).
t7t
conceptualizngthe basic geographic and environmental processes associated with land
use/land cover change, and developing spatial models that fit the spatial process of
Cellular automata (CA) have been used to simulate land use/land cover change (White,
principles, is a way to predict land use/land cover on a spatial basis, based on transition
This chapter presents the land use/land cover change prediction as well as simulation in
the study area based on Markov cellular automata models. The objectives are:
2. To simulate the future land use/land cover in study area based on scenarios
3. To analyse the past and the future land use/land cover in the study area.
Many studies have investigated land use/land cover change analysis, prediction and
modelling. These vary from those that focus on the concepts of land use/land cover
172
use/land cover change have progressed since spatial data sets from remotely sensed
data have been available. Time series remotely-sensed data with multi-temporal data
acquisition allow prediction of the timing of changes, opening new avenues better to
perturbations.
Serneels et al. (2007) assessed the driving forces of land use changes through a spatial
statistical analysis. They developed spatial statistical models of the proximate causes of
different processes of land use change. The descriptive spatial models developed in
their study suggest some important factors driving the land use that can be related to
proximate causes of the changes. It did not address the location of land use change
taking place. Pontius et al. (2001) studied land use/land cover changes related to
location and quantity issues. This study aimed at predicting the spatial patterr¡ location
and rate of tand use/land cover change. Other models that address the location and rate
of land use/land cover change include the Markov Cellular automata approach
(Veldkamp et a1.,2001). This is a spatial model that combines the Markov transitional
probability approach with cellular automata to predict land use/land cover change.
Through cellular automata techniques, behavioural models of land use can be made
The Markov Cellular automata model was selected in this study for three reasons.
Firstly, the Markov chain principle has the simple rule that transition probability can be
t73
used to predict future land use/land cover. Secondly, cellular automata can represent
the spatial dimension of the process, in this case land use/land cover states. And lastly,
Markov and cellular automata approaches complement each other to predict land
use/land cover change-addressing where and at what rate land use/land cover changes
Landscape change and spatial diffusion processes can be simulated using linear,
are appropriate for land use/land cover change processes given the complexity of
relationships between interacting variables, and the poor understanding of the driving
forces behind land use/land cover change. This approach has become common in
A Markov chain is a mathematical model for describing a certain type of process that
applied to land use/land cover change, is that past land use/land cover is helpful in
predicting future land use/land cover, given the present land use/land cover (Lambin,
1994). Bel et al. (1977) stated 'the conditional probability of land use/land cover at
any time, given all previous uses at earlier times, depends at most upon the most recent
use and not upon any earlier ones". The central mechanism of a Markov chain is a
probability pi¡ whichrefers to the likelihood of transition or movement from state i to7
in a given time interval, where i and j are either location-to model spatial diffusion
174
processes - or locationally relevant classes - to model landscape change processes
(Brown, 1970). For land use/land cover change studies, the states of the system are
defined as the amount of land covered by various land use/land cover, measured as
The Markov chain model will describe land use/land cover change from one period to
another and than use this as the basis to project future changes. This is accomplished
by developing a transition probability matrix of land use/land cover change from time
(t) to time (t+1), which will be the basis for projecting to a later time period. In
probabilistic terms for the sequence X (t), X (t+1) the Markov probability is expressed
as follows:
Plx(t+1)=ilx(t)=il
'Where j ofthe Markov series ofthe system (Bell, t974).
i and are indices
The Markov chain has been used in many studies. In the 1970s, for example, it was
used for analysing the location pattem of firms, migrations, and land use change (e.g.
Collin, 1973;Bell,1974; and Bell, 1977). Aaviksoon (1995) and Howard et al' (1995)
conducted further studies. Aaviksoon (1995) used a Markov chain model to simulate
vegetation dynamics; he concluded that the development of plant cover and land use
t)¡pes can be modelled well both on first and second order levels. The general trends
are well reproduced and the predictions are reliable to a few steps. Over longer time
periods (more than three or four steps), both approximations may fail to give conect
results. Howard et al. (1995) argued that a stationary Markov chain model allows the
future consequences of a given pattem of land use changes to be studied, and can be
175
used as an analytical tool. Logsdon et al. (1997) conducted another study on
Markov
probability mapping of land use change. They successfully showed that the
provide timely
In terms of predictioru the Markov chain models have the potential to
the Markov chain model does not help to address the "'Why?"
question. It can predict
that integration
change from the integration of remote sensing and GIS. He concludes
in
of remote sensing and GIS with Markov modelling \¡/as found to be beneficial
t76
6.2.2 Gellular automata
Spatial proximity underlies many dynamic processes in the landscape. The Game of
Life, for example, is a game based on the contiguous neighbour cell. Cells will live if 3
or 4 neighbour cells are alive and will die if 3 or 4 neighbouring cells are dead. This
Cellular automata (CA) are dynamic models that are discrete in time, space and state
(Batty, l99S). Irwin (2001) stated that a cellular automaton is a cellular entity that
independently varies its state based on its previous st¿te and that of its immediate
process. The only difference is application of a transition rule that depends not only
upon the previous state, but also upon the state of the local neighbourhood'
Cellular automata have been extensively used in many applications (White et al', 1993;
Butty, 1994; Cecchini, 1996;Xie, 1996; B"tty, 1997; V/u, 1998; Besussi et al', 1998;
Balzter et al., 1998; Candau et a|.,2000; Bryan, 2000; Sirakoulis, 2000)' White et al'
(1993), for example, simulate the spatial structure of urban land use using a cellular
automata model. Cecchini (1996) used cellular automata to build and test urban
development models, and Xie (1996) developed a generic model to simulate growth
dynamics for multi-sector land use in a sub-urban area. White et al. (1997) conducted
another study simulating the urban land use pattem using constrained cellular automata
for high-resolution modelling of urban land use dynamics. Xie (1996) integrated
cellular automata with GIS. He found that the integration a of cellular automata model
with a GIS would be advantageous because in many cases land use data that could be
177
used are already available in a GIS. Since that time, various studies have been
conducted which are related to integration of cellular automata with GIS (Takeyama et
maniputations within GIS. V/u (199S) developed Simland, to simulate land conversion
through the integrated GIS and cellular automata. He uses multi-criteria evaluation
process method. A preliminary result of this model shows the interesting properties of
macro dynamics. Batty et al. (1999) developed a model of urban dynamics through
GIS based cellular automata. He presents ways in which land uses are structured
through their life cycles, and ways in which urban activities spawn locations for new
activities.
On the urban growth application, Ward et al. (2000) developed a cellular automata
model integrated with a stochastic constraint to eliminate the broad scale factors of
constrained urban growth within GIS. The model can potentially simulate a wide
variety of urban forms depending on parameter values associated with the constraint
variables. All these studies show that cellular automata provide explicit handling of
interaction between the human systems and the natural environment within a GIS
178
Recently, studies have been conducted on developing integrated cellular automata and
GIS models that focus on spatial dynamics. Shi ¿r al. (2000) developed extended
spatial objects. The Voronoi based cellular automata can model local interactions
among spatial objects to generate complex global patterns. The Voronoi spatial model
dynamically. White et at. (2000) states that urban and regional models based on
cellular automata give good representations of the spatial dynamics of land use. He
systems.
propagation. Bryan (2000) used cellular automata to model the propagation of land
implications of the propagation process that are not considered by standard landscape
structural metrics. Sirakoulis et al. (2000) used cellular automata to study the effects of
non-linear problem.
The application of cellular automata in the simulation of spatial dynamic processes has
continued (Li Xia et al., 2001; Barredo et al', 2002; de Almeida et al', 2003; Arai et
t79
6.2.3 Prediction based on Markov Cellular automata Model within
GIS
Contemporary Geographic Information Systems (GIS) that can integrate and manage
spatial data have limited abilities to model dynamic spatial processes, and are poor at
those based on tessellation models (using raster data structures) - can be used to build
CA models. The basic form of the CA model consists of a two-dimensional array, a set
of local states, neighbourhoods and transition rules. The states of cells in the anay
undergo change according to transition rules. Transition rules are functions of the cell's
state and the state of its neighboruing cells. This basic principle is able to extend GIS
The next state of X cell is a function of the current state of X cell and the neighbour's
cells. If settlement is the state of the majority of neighbouring cells, then the X cell
will change and become settlement. The algorithm used to compute the next cell state
is referred to as the CA local rule. Usually, the same local rule applies to all cells of the
2. the width of each side of the array (w). w; is the width of theTth side of the
3. the width of the neighbornhood of the cell (d).d¡ is the width of the
neighbourhood along theTth side of the arcayi.
181
The state of the X cell, at time step (t:1), is computed according to F. F is a function of
the state ofX cell at time step (t) and the states of the cells in its neighbourhood at time
Figure 6.1. The neighbourhood ofthe (ij) cell is formed by the (ij) cell itself and the
eight adjacent cells (Sirakoulis et al., 2000).
Land use/land cover change is a dynamic spatial phenomenon. It is the result of spatial
land use/land cover type in one cell is associated with those of its neighbours. The
simple rules of CA say that the state of a cell (or cells) changes according to transition
rules, and the transition rules (a function of the cell's state and those of the
neighbouring cells) can be applied on land use/land cover change at every cell. The
future land use/land cover type (according to Markov principles) is guided by the
present land use. In other words the conditional probability of land use/land cover type
at arry given time, given all previous uses at earlier times, depends at most upon the
most recent use and not upon any earlier ones (Bell and Hinojosa, 1977). This allows
the consideration of the land use/land cover category in which a parcel or cell will be
t82
classified at time (t+1) as dependent most upon its current classification at time (t)'
Logsdon et at. (1997) conclude that Markov chains afford a powerful descriptive and
predictive model for land use change and for future land use distribution. Therefore, in
The integration of Markov based cellular automata models within GIS in order to
predict land use/land cover change is relatively new. Idrisi32 vet.2 includes a Markov-
based cellular automata model to predict land use/land cover change (Eastmar¡ 2001).
land use suitability as a constraint. In this model, the neighbouring cells that may
change are constrained by the land suitability of the neighbour cells. If neighbour cells
are suitable for settlement, for example, the cells will mostly likely change into
settlement and if neighbour cells suitable for forest it will change into forest.
6.3 Method
The Markov based cellular automata model in Idrisi32, called CA-Markov, has been
selected to model land use/land cover change in the study area. This model is based
on
the combination of cellular automata and Markov chain techniques. It adds an element
183
Three data sets are required for operation of this model. Firstly, the latest land cover
classification is used as a basis for land use/land cover prediction (see Chaptet 4).
Secondly, the Markov transition that contains the Markov transition probability matrix
of land use/land cover change. Finall¡ the distribution of land use suitability for each
land use/land cover category. The transition matrix and suitability maps will be
discussed below.
A Markovian process is one in which the state of land use/land cover at time (t+1) can
be predicted by the state of the land use/land cover at time (t) given a matrix of
transition probabilities from each cover class to every other cover class. The Markov
module in the ldrisi3.2 can be used to create such a transition probability matrix. This
matrix is a result of cross tabulation of two images adjusted by the proportional error.
As input, it takes two land use/land cover maps. It then produces the following outputs:
(Eastmar¡ 2001)
likelihood that a pixel of a given class will change to any other class (or stay
o A transition areas matrix. This expresses the total (in cells) expected to change
o A set of conditional probability images - one for each land use/land cover
class. These maps express the probability that each pixel will belong to the
184
6.3.2 Multi-criteria evaluation approach to create land suitability
Change in state of cells not only depends on the Markov transition probability, but also
on the inherent land use suitability of each pixel. The suitability map was derived
according to criteria such as proximity to roads, river, city and existing land use/land
cover. The production of these maps follows the procedure outlined in the multi-
criteria evaluation (MCE) module of Idrisi32. The criteria were developed empirically
in relation to tmderlying land use/land cover change dynamics between two different
times (in this case between 1989 and 1997). The change of settlement category for
example, was related to factors such as proximity to main roads, proximity to urban
and semi-urban centres and mostly located in area within slope less than 45Yo.
Therefore in creating the suitability for settlement, it was based on the factors of
proximity to roads and city centres as well as slope less than 45Yo and existing
settlement. Suitability for rice field was based on proximity to inigation or river, slope
less than less than 45Yo and existing rice fields. The factors that contribute to each land
185
Weighting factors have also been included in order to create the suitability map. In the
acceptable value. Using the WEIGHT module in ldrisi32, the weight value of each
Markov cellular automata combine both the concept of a cellular automata filter and
Markov model procedure. The Markov rule is that the future state can be predicted
based upon the recent state, but with the Markov model alone it provides no sense of
basis, but since there is no knowledge of the spatial distribution of occurrences within
each land use/land cover category, there is no spatial component on the modelling
outcome. A cellular automata (CA) model is used to add spatial characteristic to the
model. The CA model has the ability to change its state based upon the application of a
rule that relates the new state to its previous state and those of its neighbours. The
The contiguity filter was applied to a series of suitability maps already identifred for
each land use/land cover class. The process of the Markov Cellular automata model
was ¿Ìs follows: The first step is to create the transition probability and areas file from a
Markov Chain analysis (using the MARKOV module) of two prior land use maps,
186
establishing the quantity of expected land cover change from each existing category to
each other category in the next time period. The second step was selecting the base
image to use as the starting point for change simulation. The third step, to choose the
number of iterations according to how many years the simulation will conducted, i.e.,
if l0-year simulation, then chose 10 time iterations. Within each time step, each land
cover is considered in turn as a host category. All other land cover classes act as
claimant classes and compete for land (only from within the host class) using the
MOLA (multi-objective land allocation) procedure. The area requirements for each
claimant class within each host are equal to the total established by the transition areas
file divided by the number of iterations. The results of each MOLA operation are
overlaid to produce a new land cover map at the end of each iteration. The cellular
automaton component arises in part from the iterative process of land allocatior¡ and in
part from a filtering stage with each iteration that reduces the suitabilþ of land away
from existing areas. The filter is integral to the action of the Cellular Automata
component. Its purpose is to down-weight the suitability of pixels that are distant from
existing instances of the land cover type under consideration. The net effect is that to
be a likely choice for land cover conversior¡ the pixel must be both inherently suitable
This process is run for the number of iterations. Therefore, the result of prediction from
this model will depend on the number of iterations and filter type as well as the basic
land use/land covers data and transition probability. The filter and iteration can be
187
6.4 Selecting Parameters and Testing the Model
The Markov Cellular automata prediction model used in this study in Idrisi32 is an
experimental module. The operation of the Markov Cellular automata module to predict
land use/land cover change depends heavily on the quality of data input, such as
Markov transition probability, suitability maps and the base land use/land cover image.
The number of iterations and the contiguity filter type are other factors that have to be
considered. Therefore in order to perform suitable land use/land cover prediction these
In assessing the model in this study, the latest image that was used as a base image was
the 1989 classified image used to predict land use/land cover conditions for year 1997.
Different time intervals of the Markov probability matrix calculation were tried, as well
as dif[erent types of contiguity filters and number of iterations in order to get a suitable
result of prediction or simulation. The land suitabilþ map was also improved
especially for land suitability for settlement that very clearly associated witlì factors of
proximity to urban areas and proúmity to the main roads. The best result of prediction
is then compared with the classified 1997 image to examine the precision of the result.
The transition probability was calculated based on Markov Chain analysis (using the
MARKOV module). Many parameters are required to create the Markov transition in
this module. These options include how many years or time intervals will be used and
188
for how long the prediction or simulation will be conducted. Table 6.2 shows the
Markov transition of land use/land cover between 1989 and 1997 n the study area. It
was an 8-year difference between 1989 and 1997; therefore, the time interval selected
to produce this transition matrix w¿rs an 8-year time interval and an 8-year prediction
time. From Table 6.2 lt can be seen that the diagonal direction has the higher
probability value. For example, class I to classl (Forest), class 2 to class 2 (Plantation)
and so on. It means that the high value of transition probability occurs between the
same categories, while for a category that has no chance to be changed or low
probability to change, the value was low, such as from settlement to forest or to
plantation as well as to water/reservoir (Table 6.2). These low values make sense
unlikely.
Cells in Ex transition to
LULC ct.1 ct.2 ct.3 cr.4 cr.5 ct.6 ct.7
C1.1. Forest 57150 22300 47773 1092 116 0 26
Cl.2.Plantation 104503 234134 297490 56736 I 8865 142 0
Cl.3.Rice field 44330 287651 1985648 6101 10 271561 82092 2299
C|.4. Open/dry land 4633 80063 405784 357597 75893 2502 185
C1.5. Settlement 250 16424 68175 22583 141809 1052 75
C1.6. Fishpond 47 117 79397 4406 4476 145882 23
C1.7. Water/reservoi 1 361 581 73
Given Probabil of to
LULC c t. 1 ct.2 cr.3 ct.4 cr.5 ct.6 ct.7
C|.1. Forest o.4449 0.1 736 0.3719 0.0085 0.0009 0 0.0002
Cl.2.Plantation 0.1468 0.3289 0.4179 0.0797 0.0265 0.0002 0
189
6.4.2 Land suitability of each land use/land cover category
Suitability maps are another input parameter in Markov Cellular automata. Suitability
maps were developed in relation to the underlying land use/land cover change
dynamics in the study area between 1989-1997. Using the Multi-Criteria Evaluation
module with factors as shown in Table 6.3 and using water bodies as a constraint,
1. map
Forest suitability 4. Oper/dry land suitability map
2. map
Plantation suitability 5. Settlement suitability map
3. map
Rice field suitability 6. Fishpond suitability map
Table 6.3. The weight value of suitability factors
190
Table 6.3 above shows the weight value of factors that have been selected for every
land use/land cover category. For the category open land, for example, the existing
open land has the highest weight value, follow by existing settlement, proximity to
roads and slope. While in the settlement category, the factor proximity to road has the
highest weight value, following by proximity to city and to Jakarta metropolitan area.
These weighting processes were created based on the change characteristic associated
with land use/land cover change dynamic between years 1989-1997. For example,
settlement was mostly changed within the area close to Jakarta, urban and semi urban
Figure 6.2a and b show the suitability map for settlement category produce by multi
Lorv
+
High
Figure 6.2a. Suitable for Settlement without factor proximity to Jakarta city
t9t
Lo*'
+
High
I
::
Figure 6.2b. Suitable for Settlement with factor proximity to Jakarta city
suitability
The suitability map in Figure 6.2a and 6.2b shows the gradient of land use
for settlement. The different colours show the gradient of suitability from unsuitable
(black, red, yellow), to suitable (green) and very suitable (dark green)' It can be seen
that high suitability for settlement is located along the road and close
to cities (Figure
6.2a) and concentrated in the area close to Jakarta (on the left side image) (Figure
gradient from unsuitable
6.2b). Other category suitabitity maps have the same colour
(see APENDIX
(black, red, yellow), suitable (green) and very suitable (dark green)
VI). These suitability maps \Mere created based on examining some factors such as
study area.
192
6.4.3 The result of simulation or prediction process
Once the Markov transition area and suitability maps were constructed, the Markov
Cellular automata prediction model can be run. Land use/land cover for 1997 in the
study area was predicted by selecting the 1989 classified image as a base image,
combining it with the Markov transition and the suitability maps. In this prediction
process, different parameters such as contiguity filter type as well as suitability map
Three kernel filter types (3, 5 and 7-kernel filters) were tried in order to obtain a
suitable result. The 3x3 kernel filter was the best contiguity filter when compared with
the 5-and 7-kernel filters. The number of iteration was 8, the same as the number of
years to predict (8 years from 1989 to 1997). Different types of suitability maps have
also been used, such as before and after being adjusted for proximity to Jakarta city'
Figure 6.3a and 6.3c show the prediction result based on this parameter. Figrne 6.3a
shows the results of prediction process based on the suitability map before being
adjusted and 6.3c the prediction process based on a suitability map that has been
adjusted.
The result obtained when the suitability map was adjusted for proximity to Jakarta
(Figure 6.3c) appears better than that for the unadjusted suitability (Figure 6.a).
Settlement in areas that close to main roads and urban and semi-wban areas wers
t93
rice field as well as plantatior¡ fishpond and reservoir, look appropriate as predicted
Overall KIA of the predicted image before being adjusted was 56 %o and after being
Other land use/land cover category such as open/dry land and plantation as well as
forest have KIA value around 4l to 46 Vo, while in the image before adjustment all
categories were below 55 yo, except fishpond and water/reservoir with 7l Yo and96Vo
Plantdlwt
Rice field
0psnffly land
Fish pond
Wde/nsewoh
194
Figure 6.3. b). Original classified 1997 tmage
B ekasi
Karawang
Cikampek
Pr-rnvakarta
195
Table 6.4a .The Kappa Index of Agreement (KIA) value between
original 1997 -1997
predictedl and original 1997- 1997 predicted2'
Therefore, the KIA value between the original l9s9 and 1997 predicted images also
should have a low value - the same or close to the KAI between the 1989-197 original
t96
was better than predictedl which used the suitability map before it was adjusted. The
KIA value of predicted2 is lower than KIA value of predictedl, and has smallest
difference values with KIA value between the 1989 - 1997 original. For example, rice
value less than 0.15, and other categories such as forest, open/dry land and plantation
have difference value more than 0.31. While in predictedl all categories have
difference value more than 0.35. It means in predicted2 that the prediction result is
better than that in predictedl, especially for settlement and rice field categories. This
condition indicates that settlement and rice field was well predicted.
Modification of the Markov transition matrix improves the amount of cells that are
expected to change but did not improve the land use/land cover change on spatial
distribution. An inappropriate suitabilþ map of open/dry land category could cause the
poor result on open/dry land, forest and plantation change distribution' It was no
significant information to adjust the suitability map of this category. For example, there
are no systematic patterns of open/dry lands. Slope and proximity to tlìs roads as well
as existing open/dry land have been selected to create suiøbility maps of these
categories. A suitability map with or without existing land use category does not
improve the prediction result of open/dry land as well as forest and plantation.
r97
6.5 The potential outcome of future land use/land
influence the future rate of land use/land cover change. Therefore, to try and to
incorporate the changed circumstances, two scenarios have been developed which
will remain as idle land ("lahan tidur"), due to slow down on development
198
category as
There is no change from settlement to other land use-land
below.
The process and result ofthese simulations presented
The time frame for simulation was ten years, and the
Markov transitions for the next
1997 to 2007 is same as that between 1989 and 1997 (Markov chain
ruþ' other
t99
Table 6.5 summary of probability expected to stay as same category
lO-year scenario
From 1989-1997
lmages Scenariol Scenario2
Expected Expected Expected
unchanged unchanged unchanged
LULC
36% 60% 500Á
Forest
23o/o 50% 45o/o
Plantation
Rice field 56% 65% 60o/c
higher in
From Table 6.5 it can be seen that probability of "expected unchanged" are
10 years, land use/land cover will not change as much as between 1989-1997' For
a high probability to
to any other category from water/reservoir. Settlement also has
stay the same. Fishpond has lower value than in 1989'1997, which means it has more
use/land cover
A more detailed composition of transition probability for each land
of land use/land cover for the next 10 year from 1997 '
200
Table 6.6a. The Markov transition of land use/land cover for 10-year prediction
(1997-2007), based on 1989-1997 image'
Cells in transition to
LULC cr. 1 ct.2 ct.3 cr.4 ct.5 cr.6 cL7
70617 37349 85854 3943 277 40 59
C|.1. Forest
93561 136046 282742 19124
56372 588 0
C1.2. Plantation
53052 306456 1729823 641311 298654 88941 2497
G|.3. Rice field
8364 90688 471109 275542 90312 3571 188
C1.4. Open/dry land
C|.5. Settlement 1043 33453 142653 47324 226468 2176 181
(scenario 2007)
Table 6.6b. The Markov transition of land use/land cover for l0-year
assumed that changes willlow
Cells in transition to
LULC ct. 1 ct.2 cr.3 cl.4 ct.5 ct.6 cL7
8884 39628 29721 9907 0 0 0
1. Forest 11
201
Table 6.6c. The Markov transition of land use/land cover for 10-year (scenario 2007)
assumed that changes will very low
Cells in transition to
LULC cr. 1 ct.2 ct.3 cl.4 cr. 5 cr.6 ct.7
C|.1. Forest 99070 59442 29721 9907 0 0 0
C|.7. Water/reservoir 0 0 0 0 0 0 1
Figure 6.4 shows the predicted land use/land cover for 2007 as a result of lO-year
simulation based on the Markov transition from transition between the 1989 and 1997
images (Table 6.6a), Figure 6.5 shows land use/land cover prediction lot 2007 based
on scenariol (Table 6.6b) and Figure 6.6 shows the land use/land cover prediction
based on scenario2 (Table 6.6c). Figure 6.7 shows the original 1997 land use/land
202
Fish pond
Waìer/reseruoir
simulation based on
Figure 6.4. The result of land use/land cover change of 1O-year
Markov calculated from 1989- 1997 images (simulation2}}7)
203
Fish pond
Water/reseruolr
simulation
Figure 6.6. The result of land use/land cover change of lO-year
based on Markov scenario (Scenario2 2001)'
204
From these figures it can be seen that there are diflerences in land use/land cover
change patterns, especially for the settlement category. As it has been assumed that
land use/land cover changes for the next ten years will be low in scenariol and very
low in scenario2 compared with change between 1989-1997, the land use/land cover in
scenariol as well as in scenario2 visually was slightly different from land use/land
cover simulation for 2007 (Figure 6.4). Settlement along main roads in the simulation
does not show much difference from settlement in scenariol (Figure 6.5), but higher
Different land use/land cover changes between the original simulation and the
scenarios can also be seen by comparing the change characteristic between these
based
images (Table 6.7 and Figure 6. 8). The change characteristics were calculated
on percentage change oftotal cells for each land use/land cover category from
a cross-
From Table 6.7 and Figure 6.8 it can be seen that there are differences in the change
2007.
characteristic between simulation 2007 and scenariol 2007 as well as scenario2
40%o and increased 24% n scenario2. These change characteristics result from the
205
estates would slow down within the ten years from 1997. For example, the increased
represents a lot of idle land ('tanah tidur'). Slightly low settlement (lo/o difference) in
scenariol and very Iow (20%o difference) in scenario2 represent low activity for the
Table 6.7. The lO-year change characteristic between Simulation and Scenario.
Change characteristic Change characteristic Change characteristic
Land use- land 1997 originalto 1997 originalto 1997 originalto
cover simulation 2007 scenariol 2007 scenario2 2007
Forest 14 % increased 34olo increased 24% increased
Plantation 3% decreased 1% decreased 2% decreased
60
50
40
30
t,zo
tr
.E 10
àñ
2 4 5 7
-10 1
-20
-30 land useland covertYPe
40
E lùyear sinulaüon I lGyear sccnariol I10-¡æar scenario2
206
6.5.2 Land use/land cover Simulatio n 2017 (2O-Year)
In section 6.5.1 two scenarios were used which adjusted the transition probabilities
from the lgSg-lgg7 matrix. This was based on changing economic circumstances
following the 1998 economic crisis. The use of these adjusted probabilities probably
produced a more realistic result; however, it is not expected that these economic
conditions will remain constant.
Therefore, the simulation for both scenarios was extended to 20 years. In this case the
assumption were:
The adjusted probabilities can be seen in Table 6.8a and 6.8b as follow'
207
rabre 6 sa M*k*;åî::?üxåîr3if1:ye*iï. use/land cover for
Cells in transition to
cr. ct.2 ct.3 ct.4 ct.5 ct.6 ct.7
LULC 1
19814 9907 0 0 0
1. Forest 1 18884 49535
132397 88265 0 0 0
.2. Plantation 73554 294217
2028478 702165 234055 0 0
0 1 56037
CI.3. Rice field
394705 87955 140966 93977 0
1.4. land 28193 93977 1
0 90660 362638 0 0
5. Settlement 0 0
0 0 58848 23539 11770 141234 0
C1.6. Fishpond
1.7. Water/reservoir
Given Probab of to
cl. cl.2 cr.3 ct.4 cr.5 ct.6 ct.7
LULC 1
rable.6.sb.Markovtransitioi"ï"J31'f'::*tåîtuse/landcoverfor
Cells in transition to
cr. cL2 ct.3 cL4 ct.5 cl.6 ct.7
LULC 1
9907 0 0 0
1. Forest 148605 29721 9907
88265 0 0 0
.2. Plantation 73554 294217 132397
624147 4681 10 0 0
C|.3. Rice field 0 0 2028478
187955 281932 46989 0
4 land 0 140966 281932
0 90660 362638 0 0
5. Settlement 0 0
0 0 58848 23539 11770 141234 0
cr.6.
0 0
Water/reservoir
1.7.
Given of to
cr. ct.2 ct.3 cL4 cl.5 ct.6 cl.7
LULC 1
0.05 0.05 0.00 0.00 0.00
Cl.1. Forest 0.75 0.15
0.22 0.15 0.00 0.00 0.00
.2. Plantation 0.12 0.50
0.65 o.20 0.15 0.00 0.00
3. Rice field 0.00 0.00
0.30 0.20 0.30 0.05 0.00
C|.4. Open/dry land 0.00 0.15
0.00 0.20 0.80 0.00 0.00
.5. Settlement 0.00 0.00
0.25 0.10 0.05 0.60 0.00
.6. F 0.00 0.00
0.00 0.00 0.00 0.00 1.00
cl.7 Water/reservoir 0.00 0.00
208
When the simulation was run using these assumptions, the result showed
that land
(Simulation 2017), Figwe 6.10 shows the result of 2o-year simulation based on
scenariol and Figure 6.11 scenario2. From these figures it can be seen that land
use/land cover change was quite different, especially in scenario2 where settlement
A summary of the changes can be seen in Table 6.9 and Figure 6'L2' It can be clearly
seen that the two scenarios produce results difflerent from that of the original
simulation.
Planl¡llon
Rþê fiEld
þñd
(simulation 2017)'
Figure 6.9. The result of land use/land covef 2}-yeat simulation
209
2017)'
Figure 6.10 The result of land use/land cover 2}-year simulation (Scenariol
(Scenario2 2017)
Figure 6.1 1. The result of land use/land cover 2}-yeat simulation
2r0
Table 6.9. The 2}-year change characteristic between Simulation and Scenario
Figure 6.12 shows the comparison of land use/land cover change characteristics
scenario2, land use/land cover change was higher than in scena¡iol and the original
simulation. For example, settlement in scenario2 increased 148% and rice field
decreased 20Vo, wlfike in the original simulation and scenariol settlement increased
60Yo and 65%. This indicates that the assumption of high urbanisation in scenario2
resulted in high increased settlement and supports the assumption that within twenty
2tl
160
140
120
100
o
ctt 80
c(É
60
o 40
s
20
0
-20 1 2 4 5 6 7
40
Land useland covertYPe
The expansion of settlement along the roads was clearly recognised from the result of
simulation both in scenario2 2007 and scenario2 2017. Figure 6.13 shows an overlaid
roads network with scenario2 2007 and scenario2 2017 images. From Figure 6.13 it
can be seen that settlement increased mostly along main roads and areas that are close
to city centres, especially to Jakarta. In some places it is seen that open/dry land
212
2001
A
10 0 10 20 Kilometers
I
x Foresl
Pbnlrtlon
Rlc¡ fleld
Jrkad¡
:]I Opin¡drylmd
Scül.mênt
FlEhpond
Waler/resewok
2017
lU 0 20 Ktlomêlers
Folêsl Jrk.rlt
Plrnlallon
Rlce f¡êld
Opon dry lrnd
Såttlemêrd
Flshpond
Wdgilrasgfvolr
213
6.5.3 Analysis of Land use/land cover, 1989, 1997,2007 and 2017
It has previously been mentioned (Chapter V), that proximity to roads and proximity to
urban and semi-urban areas are the main factors influencing land use/land cover
change in the study atea, especiaþ the change of settlement and rice field categories.
Therefore, spatial analysis of the future land use/land cover in this section focuses on
land use/land cover within 10 km of main roads and 10 km of urban and semi-urban
centres. This analysis includes the 1989, 1997 classified images, 2007 and 2017
As has been mentioned previously, scenario2 2007 was based on the assumption that
land use/land cover change will slow down with very low in urbanisation due to the
economic crisis in 1998. Scenario2 2017 assumed that economic activity recovered,
and land use/land cover change increased with high rubanisation. Therefore, land
use/land cover scenario2 year 2007 and20l7 was selected, as it is assumed that these
two scenarios will represent the next 10 and 20 years of land use/land cover condition
Table 6.10 and Figure 6.14 show land use/land cover within 10 km of main roads in
214
Table 6.10. Land use/land cover within 10 km of main roads lr;r1989,1997
2007 and 2017.
rice field, and fishpond. Settlement increased. In 1989 settlement areas were less than
was more
plantations, open/dry land and rice field, but eighteen yeafs later settlement
This indicates that within l0 and 20 years from 1gg7, the agricultural areas such as
plantation along the main roads witl be less than settlement afeas'
and in 2017
of idle land ("lahan tidur"), due to slow down in development activity,
recovered'
activities increased as it was assumed that economic activity
215
1 80000
160000
140000 ! Forest
@ Plantation
120000
ag A Rice field
- 1 00000
g Open/dry land
(ú
o 80000
I Settlement
60000
n Fishpond
40000 6 Water/reservoir
20000
0
1989 ',1997 2007 2017
Year
Figure 6.14. The tend of land use/land cover change within 10 km of the main roads,
1989 to 2017 (2007 and20l7 the result of scenario2)
2t6
1989
I 9S7
20 0 20 Kllometers
Foroct
'
Fi_Ijil Jd<arta
r Plairtaüon
Rice field
openidrv land
Pror¿incid Road
StEte Road
Sêttlernent Toll W+'
Fishpond Rail Waf
Wder/reservoir
2t7
2007
2011
2û Kilomêters
20
^
0
Pldntdtion
trf,i-ì Jd(dttr
Rice tield Prw¡nc¡d Roed
Opeildry lÊnd $ete Road
Seül€mPnt TqllWåy
Fishpøtd Rail Wey
roads
Figure 6.15b. Land use/land cover within 10 km of the main
218
Settlements continued to increase and agricultr.ral areas to slightly decrease in the area
within 10 km of the main roads. It was more drastically changed in areas within 10 km
of the urban and semi-urban center (Figure 6.16 u b and c). For example, in area
within 10 km of Bekasi urban center in 1989, settlement was still less than rice field,
but after 1997 settlement was greater than agricultural areas. By 2017 rice field will be
nearly zero inthis area as it was urbanized to adjoin the Jakarta metropolitan city. The
semi-urban centers next to Bekasi are Cileungsi and Cikarang. In these two semi-urban
ateas, rice fields in years 1989 to 2007 was greater than settlement, but then decreased,
and finally between years 2007 and 2017 settlement increased and exceeds the rice
field area. In2007 Bekasi, Cikarang and Cileungsi would become one big settlement
In other urban and semi-urban centers such as Karawang, Cikampek and Purwakarta,
rice field was slightly decreased and :urrrtil 2017. Rice field was still greater than
settlement afeas, except in Karawang, where settlement areas in 2017 werc same as
219
Settlements continued to increase and agricultural areas to slightly decrease in the area
within 10 km of the main roads. It was more drasticalþ changed in areas within 10 km
of the t¡ban and semi-urban center (Figure 6.16 a, b and c). For example, in area
within 10 km of Bekasi urban center in 1989, settlement was still less than rice field,
but after 1997 settlement was greater than agricultwal areas . By 2017 rice field will be
nearly zero'tnthis area as it was urbanized to adjoin the Jakarta metropolitan city. The
semi-urban centers next to Bekasi are Cileungsi and Cikarang. In these two semi-urban
areas, rice fields in years 1989 to 2007 was greater than settlement, but then decreased,
and finally between years 2007 and 2017 settlement increased and exceeds the rice
field area. In2007 Bekasi, Cikarang and Cileungsi would become one big settlement
In other urban and semi-urban centers such as Karawang, Cikampek and Purwakarta,
rice field was slightly decreased and vrfti| 2077. Rice field was still greater than
settlement areas, except in Karawang, where settlement areas in 2017 were same as
219
25000 Bekasi
20000
rú
- 1 5000
(!
o 10000
5000
0
1989 1997 2007 2017
Year
+Forest +Plantation +Rice field +eopen / dry land
+Settlement +Fishpond +Water/ reservoir
25000
Gileungsi
20000
aú
1 s000
aú
0, 1 0000
5000
0
1 989 1997 2007 2017
Year
Figure 6.16a. Graphic of decreasing and increasing land use/land cover within 10 km
of urban and semi-urban areas.
220
25000 Gikarang
20000
$ ''uooo
E roooo
5000
0
1989 1997 2007 2017
Year
25000 Karawang
20000
-g 1 5000
G
o 1 0000
5000
0
1 989 1 997 2007 2017
Year
Figrue 6.16b. Graphic of decreasing and increasing land use/land cover within l0 km
of urban and semi-trban
221
25000 Gikampek
20000
.ú
1 5000
aú
o 1 0000
5000
0
1 989 1 997 2007 2017
Year
25000
Pururakarta
20000
f(! rsooo
o 1 0000
5000
0
1989 1997 Year 2007 2017
Figure 6.16c. Graphic of decreasing and increasing land use/land cover within 10 km
of urban and semi-urban
222
1989
€
^*o r
f
2Ð07
¡l
A
n A 20 Kíþmeters
223
-
6.6 Discussion
The Markov chain model is a powerful descrþive and predictive tool for land use
change and for future land use distribution. The Markov-Cellular automata model was
successful in predicting land use/land cover change in the study areq especially for
settlement and other categories that located close to main roads or trban centres. This
model depends heavily on the quality of data input, such as Markov transition
probability and suitability maps. The number of iterations and the contiguity filter type
Different combinations of filter types and number of iteration did not improve the
result. Selecting different filter types shows that the groups of cells/pixels become
more compact but did not improve the KIA value. The only efÊective way to improve
the prediction results was by adjusting the suitability maps. The adjusted suitability
map of settlement with added proximity factor to Jakarta and other urban areas,
improved the predicted result for the settlement category. It was shown that a modified
suitability map improves the spatial distribution of land use/land cover change.
Diffrculty in adjusting the suitability maps of open/ dry land and other categories such
as plantation and forest, resulted in poor prediction of these categories. In this case,
category in the study area needs to be considered. For example, the change
characteristics of settlements was known to relate to the proximity to main roads and
cities; therefore, factors such as proximity to roads and proximity to the cþ were
224
appropriate factors in creating a suitability map of the settlement category. Change
characteristics of open/dry land, plantation as well as forest were not clearly known.
Factors such as slope and proximity to the roads that have been used in this study are
not really dominant factors in the suitability maps for these categories. Further study is
need with regard to obtaining more suitable prediction results for categories that still
Although there were some erroneous results for some land use/land cover categories
such as forest, plantation and open/dry land, good results were obtained for settlement,
rice field and fishpond categories especially in areas close to main roads and urban
centres. Based on this outcome, further work was undertaken to predict or simulate the
future (10- and 2}-year) land use/land cover cbange in the study area. Two Markov
transition matrices for l0-year as well as 2}-year simulation were created. One was
from Markov transition based on transition land use/land cover change from 1989 and
1997 images, and the other one based on scenarios. A transition scenario was created
that was assumed to be representative of future land use/land cover change in the study
area. The economic crisis of 1998 in this country was also considered in creating the
that the future land use/land cover change was guided by the transition from the past to
the recent/present land use/land cover condition (Markov principles). The future land
use/land cover change simulation based on the Markov transition from 1989-1997 was
successful in simulating land use/land cover 2007 and20l7. The economic crisis in
225
1998 was considered as an efflect in declining some economic activities such as
order to obtain the representative transition for the next lO-and 2O-years of land
use/land cover change from 1997. Two Markov transition scenarios for 10-and 2D-year
2007 assumed that land use/land cover will be slightly lower and scenario2 assumes
very low urbanisation and high change within agricultrual area. Scenariol of the
simulation to 2017 assumed that tand use/land cover change will slightly increase and
The results show that the intensive land use category in scenariol 2007, such as
settlement, stightly increased while rice fields decreased. In scenario2 2007 there was
very low change especially for the settlement category; it is assumed in this scenario
that economic activity recovered and the change of land use/land cover was higher than
lr¡'2007.
These results indicate that Markov Cellular automata model can be used to simulate the
future land use/land cover change based on transition scenarios as far as the transition
Further analysis from overlay of land use/land cover within l0km of main roads, wban
and semi-wban centre has shown that land use/land cover, especially settlement and
rice field, changed for the years 1989 , 1997,2007 and20l7. In an area within 10 km
226
of main roads between Jakarta to Purwakarta" settlement increased and rice field
decreased, while plantation, forest and fishpond slightly decreased. In the area within
10 km of urban and semi-urban centre that are close to Jakarta such as Bekasi,
Cikarang and Cileungsi, settlement and rice field were drastically changed. In these
three areas, settlement was larger than agricultural area and n 2017 these areas will
become one a big settlement area adjoining Jakarta city. This phenomenon can easily
be understood, since proximity to main roads as well as urban and semi-urban centres
are the main factors that influence land use/land cover change, especially in settlement
The present as well as the past and the future of land use/land cover information
resulting from this study have implications for the spatial planning process in this area.
Information such as that the future settlement area will increase at the expense of
agricultgral areas such as rice field, indicates that these areas will be urbanised. This
eliminate future conflict in land use arrangement or land allocation as part of the
capacity with regard to the expansion of settled areas over agriculttral areas.
the main
Other important information from this result is that settlement increased along
roads and in areas that are close to the Jakarta metropolitan city, as the model
rules
stated. This indicates that the development of housing and new settlement
in the fringe
227
suroundings areas. This also indicates that new settlements surrounding Jakarta area
socio-economically heavily dependent on the core, and has implication on the change
Another implication is that exceed settlement areas the agricultural areas, for example ;
or a lot of open land is idle from this simulation, shows that environmental conditions
Simulation of the future land use/land cover change based on the transition Markov
transition for the future land use/land cover change, showed that the future land
that the future land use/land cover can be simulated based on some policy scenarios in
order to obtain the optimal land use planning in this area, which could have a
positive
The success of the Markov Cellular automata model to simulate the futtlre land
use/land cover in this study area also has a contribution on development of land
6.7 Gonclusion
This chapter has provided the results of an application of the Markov Cellular automata
model to predict land use/land cover n 1997 from image 1989, as well as simulation
228
and analysis of the future land use/land cover change in the study area. The result
indicates that more appropriate suitability maps especially for forest, plantation and
The prediction result shows that overall KIA value is 56Vo before the suitability map
was adjusted and 67Yo after it was adjusted, with KIA value of forest, plantation and
open/dry land below 50%, KIA values for settlement and rice field 67Yo and 56o/o, arrd
KIA values of fishpond and water being 70o/o and 90o/o. The KIA of land use/land
cover category with a value below 50% shows the prediction result of this category
was not good enough in the categories of open/dry land, forest and plantation'
Settlement and rice field categories were adequately predicted, with a KIA value over
50%. The poor prediction result on land use-land categories such as open/dry land,
plantation and forest was due to limitations in the information of the background
important in order to select factors and criteria on creating a suitability map. The
information of settlement growth along the roads and in utban and semi-urban areas,
for example, was the appropriate information in order to select factors such as
proximity to roads and city in creating a suitability map of settlement category'
This model was successful in simulating or predicting land use/land cover change
cover for
especially for settlement and rice field categories. The future land use/land
growing
years 2007 and 2077 as a result o the simulation shows that settlement was
has an
along the main roads and close to Jakarta metropolitan city. This simulation
229
implication that there \¡ras conversion from rice field into settlement especially in area
As a Markov transition principle, it was assumed that transition change from the past to
the recent time stayed the same for the future land use/land cover change, it has
weakness since the future may have significant change in economic conditions that
could affect the land use-land change behaviour. Modifred Markov transition with
some scenarios successfully created the new simulations that could represent the future
From some simulations it was found that Markov Cellular automata allow one to
simulate the future land use/land cover change, based on the Markov transition created
from information from the past to the recent time of land use/land cover change, and
also based on scenario. The appropriate background information that could have an
effect on land use/land cover change is needed in order to create a suitable transition
scenario that can represent the future land use/land cover change.
230
Chapter Seven
CONCLUSION
The main aims of this thesis were to develop a methodology effectively to detect land
use/land cover change, to understand the inter-annual dynamics of land use/land cover
changes and to analyse, as well as to predict, land use/land cover change dynamics
using remote sensing and GIS. The Post-Classification Comparison approach was
chosen because it addresses the relevant and important aspects of change: detectior¡
identification and location. A combination of PCA and NDVI transformed images was
selected instead of the original band to classifr/identifu land use/land cover in the
study area. PCA gives a strong spectral signature of settlement or built-up area, while
NDVI gives a strong spectral signature of vegetation cover. It was hypothesised that
through the two transformed images, the detail of land use/land cover category in the
Chapter 2 reviewed and defined land use/land cover change, and described some
aspects that related to land use/land cover change. It was argued that population or
land use/land cover in the study area. The demand and supply relationship of the land
use/land cover change model shows that the demand structure of land use is very
dynamic; the demand for using land for human activities increases fromtime to time
231
under the influence of population growth, community structure and economy. On the
supply side, structure is relatively static; the surface area remains constant. The
constant size of the surface area with increasing demands being placed upon it may
result in land use conflict. Land use is mostly converted from less intensive uses to
more intensive uses, such as from forest or open/dry land agriculture to rice field or to
settlement, but it is also converted from one intensive land use to other, such as from
rice field to settlement or from intensive agriculture land into non-agriculture land use,
economic growth, demands from the non-agricultr¡re sector are usually much stronger
than those fiom agriculture, due to the much higher incomes which this sector can
obtain from the same area of land. In the long tenr¡ this could result in scarcity of food
supply since insufficient land is available for agriculture activities. Land use allocation
and management supported by the appropriate techniques such as remote sensing and
Chapter 3 describes the socio-economic and biophysical conditions in the study area.
Physiographically the study area is a flat alluvial plain. Land use/land cover in the
study area has changed over time as an effect of demographic pressure and economic
growth, especially because the area is close to the Jakarta metropolitan arca. Land
use/land cover along the main roads from Jakarta to Purwakarta is very dynamic and
has become more intensive over time. S/ithin this area many rice fields have been
converted into settlement or builfup areas for housing and industrial estate
development.
232
Chapter 4 presented the analysis and results of land use/land cover change detection
using image processing and GIS. The results of the Maximum Likelihood supervised
classification of PCA and NDVI transformed image were successful, with a Kappa
index agreement (KIA) value >85yo. Land use/land cover was clearly differentiated
and able to be classiflred into land use/land cover categories namely: 1) Forest; 2)
land use/land cover change in the study area. At the annual interval, detection found
that weather or season was related to a "leaf on and leaf off' phenomena. In the dry
seasor! forest was less than plantation because of a lot of "leaf off' (leaf had fallen),
while in the wet season forest was larger than plantation. In this time leaf leaves had
regrown ("leaf oni'). Therefore the change of forest area in the study area did not
"leaf off - leaf on" phenomena. It was clearly recognized in the field that there are
present some broad leaf plantations such as rubber plantations and teak wood,
especially in the southem part of Karawang to the north of Jatiluhur reservoir, which
On the other hand, rice field (planted) and rice field (unplanted) were alternates
associated with the growing cycle. It seen at harvesting time that the land was covered
with rice plantation ready to be harvested, while at planting time there were no rice
plantation due to the fact the land is under preparation. Therefore at harvesting time
233
rice field (planted) was recognised to be larger than rice field (unplanted) but at
planting time rice fietd (unplanted) was larger than rice field (planted). This
growing cycle with changes of other land use/land cover was not recognised, except
the change in open land and dry land agriculture associated with dry and wet season.
From annual detection it was recognised that rice field (planted) and rice field
(unplanted) changed over time in association with the rice growing cycle. For the long-
term change analysis, it was necessary to reclassifr or regroup these two categories
including both the open land and dry land categories. Rice field (planted) and rice field
(unplanted) have the same function in the field as rice field area, and so were included
as rice field. Open land and dry land categories also have the same function as land
which is unused or under preparation for construction or for dry land agriculture, and
so were regrouped as open/dry land. After regrouping of these land use/land cover
categories it was found that settlement and rice field areas had permanently changed
over time. Settlement had increased while rice field had decreased. The trajectory of
land use/land cover change in this study area was recognized as follows: Rice field was
converted to open/dry land or directly to settlement, and open/dry land was converted
settlement.
Chapter 5 commences the analysis of factors influencing land use/land cover change.
The large number of possible elements were simplified into static and dynamic factors.
234
Static factors were defined as factors that might never change, such as slope, elevation,
or physiography, while dynamic factors were defined as factors that can cltange, such
transport
as population growth, population density, urban, and semi-urban areas and
routes. Overþ analysis of these two factors with land use/land cover found that
cover
roads are the dynamic factors that have a strong relationship with land use/land
in the study area. Settlement increased and rice field decreased significantly with
and semi-urban centers and proximity to roads. As the study area is mostly a
flat
(settlement
alluvial plain, land use/land cover types that represent human intervention
and rice field) were generally located in flat areas with a slope of less
than3o/o.
In urban areas such as Bekasi, Karawang and Purwakart4 decreasing rice field and
such as in
increasing settlement are clearly recognized as well as in semi-urban areas
V/ithin an 8-year time interval (1989-1997), there was a gradual change from low
population
population density to high population density. Settlement in the areas with a
density of 151-650 person/sq.km (class II) increase d 64Vo, and in areas with population
density of 651-1200 person/sqkm (class III) increased of 134%. The rice field area
decreased by -1 1% nareas with the lowest population density (class I), in class II it
decreased by -3lo/o, and in class III it decrease d by -63%. This indicates that the
rice
increasing population density is an eflect of increasing settlement and decreasing
235
field in the study area. There is no indication of a strong relationship between
increasing population density and decrease or increase in forest and plantation area.
Forest increased in areas with low population density (class I and class II), while
plantation decreased in areas with low population density (class I) and increased in
areas with high population density (class II and III). The increasing plantation areas in
areas of high population density could be due to a lot of plantations being located close
to settlement areas. This is the case with rubber and teak wood plantations that are
Another finding is that distance from the Jakarta metropolitan area has a strong
influence on increasing settlement and decreasing rice field in the study area" The
urban and semi - wban areas gradually decreases with increased distance from Jakarta.
In the Bekasi area 15 km away from Jakarta, occutred the greatest increase in
settlement density and the largest decreases in rice field, following by Cileungsi and
Cikarang, which are 30 km from Jakarta. This supports the contention that the Jakarta
metropolitan area exerts an influence on land use/land cover change on the surrounding
areas.
A positive relationship was found between total population and settlement, while the
relationship between total population and rice field, forest and plantation was negative.
It was found that for an increase in total population, the area under settlement will
increase and the area under rice field, forest as well as plantation will decease.
236
Discriminant analysis supported the strong relationship between the static and dynamic
factors that influence land use/land cover. Factors such as population density,
proximity to main roads, proximity to urban centre and slope are the factors that
and simulate future land use/land cover change in the study area. It is shown that by
selecting different filter types, the groups of cells/pixels become more compact but it
did not improve the KIA value. The only eflective way to improve the prediction
results was by adjusting suitability maps. The adjusted suitability map of settlement
which adds a proximity factor to Jakarta and other urban areas improved the predicted
The diffrculty experienced in deriving an adjusted suitability map of openl dry land and
for other categories such as plantation and forest, resulted in poor prediction of these
each land use/land cover category in the study area needs to be considered.
Markov transition probability was effective in determining the amount of area that can
be expected to change, while suitability maps are eflective in determining the pattem
of land use/land cover change distribution. Two different Markov transitions have been
used to simulate the future land use/land cover in the study area. The two Markov
transitions are a Markov transition based on information that was obtained from 1989-
1997 images (assuming that the future land use/land cover change transition stay the
237
same as that transition between 1989-1997), and a Markov transition based on a
scenario that was assumed to represents the future land use/land cover change
transition, considering the economic crisis since 1998. It was found that the future land
use-land cover change could be simulated based on the Markov transition from 1989-
1997 'mrages and the Markov transition scenario. The future land use/land cover from a
simulation based on the Markov transition scenario can be more or less changed
depending on the scenario assumed if comparing with the simulation based on Markov
rules (transition from the past land use/land cover condition to the recent condition). It
found that the trends of change which are characteristic for each land use/land cover
category in the two simulation images are diflerent. In the scenariol 2007 image as a
result of simulation based on a Markov scenario that assumed the future land use/land
cover will slightly change related to the economic crisis 1998, it was shown that the
trend of land use/land cover change was lower than that simulation based on Markov
transition 1989-1997.
a Markov scenario I and 2 was low compared with percentage change in the
amount of idle land ('tanah tidur" literally "sleeping land") within the 10 years from
2017 basedon scenariol and scenario2, the percentage change increased especially for
intensive land uses such as settlement, and rice field decreased, as it was assumed that
238
within the 20 years from 1997 Iand use/land cover change will increase due to
It was found that the future land use/land cover within 10 km of main roads and urban
increase and agricultural area slightly decreasing. Change was more drastic in areas
within l0 km of the urban and semi-urban areas. For areas within 10 km of Bekasi
urban center for example, settlement area in 1989 was less than rice field, but since
1997 settlement was more than agricultural area, and no more than rice field area in the
year 2017.In two semi-urban areas next to Bekasi, Cikarang and Cileungsi, it was seen
that rice field in 1989 was greater than settlement, but it continued to decrease and
finally 1¡-2007 settlement exceeded the rice field area. 1n2007 Bekasi, Cikarang and
Cileungsi will become one big settlement area adjoining the Jakarta metropolitan city.
which has the largest urban population in Indonesia. Construction of new housing or
industry in this area still continues, as well as the construction of the toll way between
Cikampek-Purwakarta- Padalarang.
Finally the results of future land use/land cover change simulations found that the
Markov Cellular automata model enables us to simulate the future land use/land cover
change. This simulation can run based on the Markov transition that is created from
information relating to past and recent land use/land cover change (Markov rules) and
based on a Markov transition scenario as far as the scenario could represent future land
use/land cover change. In this case, the appropriate background information and
239
significant factors that could have an eflect on the future land use/land cover change
7.2 lmplications
7 .2.1 Policy
Housing and industrial estate expansion has had an influence on the booming
development around the big cities of Indonesia, especially in Jakarta and its adjoining
areas- JABOTABEK. This pattern is being repeated around other mega city centres in
Southeast Asia such as Manila, Bangkok and Ho Chi Min city. Up until July 1997 the
Indonesian National Land Agency (BPN) had issued permits for development of about
private developers, of whict¡ developers have acquired 46000 hectares (FirmarU 1998).
This policy has had the implication that land use/land cover in this area has changed
rapidly. The recent physical development of Jakarta and its periphery is characterised
land to urban land uses such as settlement was clearly recognised in the land use/land
cover change detection in this study using satellite image for the years 1989 to 1997.
On the other hand, the future land use/land cover simulation obtained by using a
Markov Cellular automata model indicates that urbanisation in this area will continue,
from both a simulation based on Markov transition of land use/land cover from 1989-
1997 añ a Markov transition scenario which assumes that the economic crisis since
1998 will affect the rate of land use/land cover change in the study area. This
240
urbanisation will also happen around other big cities in Indonesia such as Surabaya'
l. Land use/land cover change detection using remote sensing and GIS technique
is useful to obtain the recent land use/land cover condition on the spatial base.
2. Detection of land use/land cover change by using remote sensing and GIS as
well as simulation of the future land use/land cover change are useful to support
3. The government has to consider that rapid land use change conversion from
agricultural land to urban uses that were found in this study can be used as an
input in order to arrange the land permits issued and to avoid a lot of idle land
As satellite images have different levels of spatial, as well as spectral, and temporal
resolution, further study is needed in order to compare the results from this study with
results from studies with more detailed classes of land use/land cover categories. Soft
classification using fuzzy classifiers, for example, could be used to try and obtain
better land use-land cover classification within agriculture land uses or within urban
land use whicl¡ in this study, was generalised to avoid mixed classification.
24r
Markov Cellular automata model based on Markov transition from the past to the
recent land use/land cover condition has a limitation on representing the future land
use/land cover change if the past has significantly different condition from the future,
such as declining or increasing economic conditions that could have an effect on land
use/land cover change. This limitation can be eliminated by using a modified Markov
The further study on this simulation needs to obtain the future land use/land cover
change which is simulated correctly, simulated based on the Markov transition that can
represent the future land use/land cover. For example, a Markov transition scenario
that was created based on the accurate and appropriate background information that
could affect the future land use/land cover change, such as information of relationship
between population growth or economic growth and land use/land cover change rate.
There is a pressing need to do similar studies in other mega SE Asia cities and other
smaller Indonesian cities such as Surabaya, Medan and Makasar. Further research also
land use/land cover change based on the Markov Cellular automata model.
242
REFERENCES
Adger W. Neil and Katrina Brown (1994). Land use and the causes of Global
Warming. John S/iley, England.
Alba Ma. J.B, Flaviana.H., Anita. J., and Estella. B.J. (1997). The Philippine Land
use and Landover change case study. SARCS-START-IGBP, Bangkok, Thailand.
Anderson J.R., E.E. Hardy, J.T. Roach, and R.E.Witmer (1976). A Land Use and
Land Cover Classification System for Use with Remote Sensor Data, Profesional
Paper 964,U5. Geological Survey, Reston, VA.
Angelici G.L., Brynt, N.A, and Friedman, S.Z. (1977). Techniques for Land Use
Change detection using Landsat imagery. Proceeding of the American Society of
Photo grametry, 217 -228 .
Arai T and Tetsuya A. (2003). Empirical analysis for estimating land use transition
potential functions-case in the Toþo metropolitan region. Computer, Environment
Balz¡er H., Paul'W. Braun and V/olfgang Kohler (1998). Cellular automata models
for vegetation dynamics. Ecological Modelling, 107: ll3-125.
243
Baredo J.I, Marjo.K, Nialt.Mccormick, carlo.L.(2002). Modelling dynamic
spatial processes: simulation of urban future scenarios through cellular automata.
Landscape and Urban Planning, 981 :1 -16'
Batty M and Kim K.S. (1992). Form follows function: reformulating urban
population density function. Urb an Studie s, 29 :l 043'107 0 .
Batty .M and Y.Xie (1994). From cell to Cities. Environment and Planning B:
Planning Design, 2l : 83 I -848.
Batty M, and Y. Xie (1997). Cellular Automata as the basis of integrated dynamic
regional modelting. Environment and Planning B: Planning Design,24:235'246'
Batty M. (1993). Urban evolution on the desktop: Simulation with the use of
extended cellular automata. Environment and Planning, A.30:1943-1967 .
Batty M, Yichun Xie and Zhanh Sun (1999). Modelling Urban dynamics through
GlS-based cellular aatomata.Computer, Environment and Urban System,
23(3):205-233.
Baulies X. and Szejwach, G. eds. (1997), Survey of need, gaps and priorities on
data for land use -land cover change research. LUCC Data Requirements
Workshop, LUCC Report Series No.3, Barcelona, Spain.
244
Beberoglu, C. D.Lioyd, P.M Atkinson and P.J. Curran (2000). The Integration of
Spectral and Textural Information using Neural networks for Land Cover Mapping
in the Mediterraneat. Computer and Geoscienc, 26:385'396
Bell E.J. (1977). Markov analysis of Land use change: continuous time and
Bell E.J and Hinojosa, R.C. (1977). Markov anaþsis of Land use change:
continuous time and stationary processes. Socio-Econ Planning Science. Il:.13'17.
Bicik lvan, Leos.J., Vit.Stepanek.200l Land use change and their Social
driving forces in chechia in the 19th - 20ú Centuries. Land use Policy,
18:65-73
245
Blaikie P, Cannon.T, Davis.I and ìVisner.B. (1994) At Risk: Natural Hazards,
People s vulnerability and Disaster, London.UK.
Blackard J.A and Denis J.DDD. (1999). Comparative accuracies of artificial neural
networks and discriminant analysis in predicting forest cover types from
cartographic variables. Computers and electronics in Agriculture.24:131-151.
Borak J.S., E.F.Lambin and A.H.Strahter (2000). The uses of temporal metrics for
land cover change detection at coarse spatial scales. International Journal of
Remote Sensing, 2l : l4l5 -1432.
2l:817-822.
246
Brown .D.G, B.C.Pijanowski and J.D.Duh.(2000).Modeling the relationships
between land use and land cover on private lands in the Upper Midwest, USA.
Journal of Environmental management, 59 ; 247 -263 .
Byme G.F., P.F. Crapper, and K.K. Mayo (1980). Monitoring Land Cover Change
by Principal Component Analysis of Multitemporal Landsat Data. Remote Sensing
of Environment, l0 : 17 5 -184.
Candau .J, Steen. R, Keith C.C. (2000). A couple cellular automaton model for
land use/land cover dynamics. The 4ú Intemational Conference on Integrating GIS
and Environmental Modelling (GIS/ÆM4): problem, Prospects and Research
Chapin F.S. (1972). Urban Land use planning. University of Illinois, Urbana
Charlie N, Suwit. O, Kamerun .M. (1997). The Thailand Land use and Land cover
change case study. SARCS-START-IGBP, Bangkok, Thailand.
Chou Kue-Chen and David S/. Elrod (1998). Using Discriminant Function for
Prediction of Subcellular Location of Prokaryotic Proteins, Biochemical and
B i ophy s i c al Re se ar c h C ommuni c at ion. 252 :63 -68 .
247
Cecchini A. (1996). Urban modelling by means of Cellular Automata: generalised
urban automata with the help onJine (AUGH) model. Environment and Planning
B : Planning Design, 23 :721'7 32.
Cihlar J., Qingþan Xiao, J. Chen, J. Beaubien, K. Fung and R. Latifovic (1998).
Classification by Progressive Generalization: a New automated Methodology for
Remote Sensing Mutti Channel data. International Journal of Remote Sensing, 19:
2685-2704.
clarke K.c., søcy Hoppen and Leonard J. Gaydos (1997). Methods and
Techniques for Rigorous Calibration of a cellular automaton model of urban
growth.
Cloke Paul.J. (19S9). Rural Land use planning in developed nation. Unwin
Hyman.ltd, UK
Coucelis H. (1997). From Cellular Automata to Urban models: new principles for
model development and implementation. Environment and Planning B: Planning
Design, 24:165-174.
Cohen W.B. and Maria Fiorella (1999)." Comparison of Methods for Detecting
Conifer Forest Change with Thematic Mapper Imagery". in Lunetta and
Christopher D.E.(eds). Remote Sensing Change Detection: Environmental
Monitoring Methods and Applications. Sleeping Bear, UK.
248
Coiner J.C. (1980). Using Landsat to Monitor Change in Vegetation Cover Induced
by Desertification processes. Proceedings of the l4th International Symposium on
Remote Sensing of Environment, AnrrArbor, Michigan.
Collin Clark (1977). Population growth and Land use. Macmillan. London.
Colwell J.E. and F.P. Weber (1931). Forest Change Detection. Proceedings of the
l5th International Symposium on Remote Sensing of Environment, Attrt Arbor,
Michigan.
Dawson A.H. (1934). The Land problem in the developed Economy. Croom Helm,
NSV/, Australia.
deBie CA, van Leeuwen JA and Zuidema PA.(1996). The Land Use Database. A
knowledge-based software prograrnme for structured storage and retrieval of user-
defined land use data set. Users' Reference Manual (v104, DOS). ITC, Enschede;
FAO, Rome: Agricultural University,'Wageningen.
249
Djoko S. (1996). Problems and prospects of Indonesia new towns development.
Paper presented at the International Seminar Planning and Development of
Satellite Towns in Southeast Asia. Organized by LIPI-The Japan Foundation,
Jakarta 25-27 June 1996.
Dobson, M. C., Pierce, L. E. and Ulaby, F.T. 1996, Knowledge-based land cover
classification using Ers-l/Jrs-l SAR composites, IEEE Transactions on Geoscience
and Remote Sensing, 34(l), 83 -99.
University, USA.
2s0
ERDAS. (1997). Image Manual V.8.3, 4ú edn. ERDAS, Inc. 2081 Buford
Highwas, NE, Atlanta, Georgia, USA.
Firman Tommy. 1997. Land Conversiton and Urban Development in the Northem
Region of West Java, Indone sia. Urb an Studie s, 3 4(7), I 027 -1046.
Fischer Gunther and Laixiang Sun. (2001). Model based analysis of future land-use
development in China. Agriculture, Ecosystems and Environment, 85: 183-176.
Fisher P.F. and S.Pathirana (1990). The Evaluation of Fuzzy Membership of Land
Cover Classes in the Suburban Zone. Remote Sensing of Environment, 34: l2l-132.
Fisher P.F. and S.Pathirana (1993). The Ordering of Fttzq Land - Cover
251
Foody G.M. (2001). Monitoring the Magnitude of Land - cover change around the
Southern Limits of the Satrara. Photogrametric Engineering and Remote Sensing,
67:841-847.
Foody G.M and D.S. Boyd (1999). Detection of partial land cover change
Foran B.D. (1937). Detection of Yearly Cover Change with Landsat MSS on
Pastoral Landscapes in Central Australia. Remote Sensing of Environment, 23:333'
350.
Frederick M., Iwan G., Handoko., Asep Karsidi, M. Muchlis., Andi Rahmadi.,
Hartanto Sanjaya., Laju Gandarum. (2001). Case Studies and Diagnostic Models of
the Inter-a¡nual Dynamics of Deforestation in Southeast Asia: Is the missing sink
for Carbon in Land cover Change?. Indonesian Case Study. Final report submitted
to: The United Nations OfEce of Project Services (LINOPS).
Fung T and W. Siu (2000). Environmental qualþ and its changes, an analysis
using NDVI. International Journal of Remote Sensing, 2l:I0ll'I024.
252
Fung T. and E. LeDrew (l9SS). The Determination of Optimum Threshold Levels
for Change Detection. Photogrametric Engineering and Remote Sensing, 54:.1449-
1454.
Gemmell F., Jari Varjo and Mikael S. (2001). Estimating Forest cover in a boreal
Forest test site using Thematic Mapper data from two dates. Remote Sensing of
Environment, 77 : 197 -2ll .
Gobin 4., P. Campling, J. Feyen (2002). Spatial analysis of rural land ownership.
Landscape and Urban Planning,55: 185-194.
Goldblum C., Tai-Chee Wong (2000). Growth, crisis and spatial change: a study of
haphazard urbanization in Jakarta, úrdonesia. Land Use Policy, 17:29-37 .
'Warren
Gordon D.K, R.P. and \I/illiam D.P. (1986). Fruit Tree Inventory with
Landsat Thematic Mapper Data. Photogrametric Engineering and Remote Sensing,
52: l87l-1876.
Green. K., Dick Kempka, and Lisa Lackey (1994). Using Remote Sensing to
Detect and Monitor Land Cover and Land Use Change. Photogrametric
Engineering and Remote Sensing, 60: 33I -337 .
GTZ. (2000). Final Report Private Public Parbrership (PPP) Project to Promote
Spatial Planning and Forest Mapping, Jakarta, Indonesia
253
Hashiba H., K. Kameda, S. Uesugi and S. Tanaka (2000). Land use change analysis
of Tama river basin with different spatial resolution sensor databy Landsat MSS
and TM. Advance Space Research, 26 (7):1069-1073.
Hathout S. (2002). The use of GIS for monitoring and predicting urban growth in
East and West St Paul, V/innipeg, Manitoba, Canada. Journal of Environmental
Mana ge me nt, 66 :229 -23 8.
Haub Carl (1996). Future Global Population Growth, in Ramphal, Shridath and
Steven W.S (eds). Population Growth and Environmental issues.Praeger
Publishers. London.
Helmschrot J., Wolfgang Albert Flugel Q002). Land use characterisation and
Henderson, F.M. and Xia, Z.G. 1997, SAR applications in human settlement
detection, population estimation and urban land use pattern analysis - a status
Hidding M.C and Andre, T.J.T. Q002). Beyond fragmentation: new concepts for
urban-rural development. Landscape and Urban Planning, 58:297'308.
254
Hoover E.M and Frank Gianatani (1999). An lntroduction to Regional Economics.
htþ ://www.rri.wvu.edu/WebBook/Gia¡ratani.
Howard D.M, P.J.A. Howard and D.c. Howard. (1995). A Markov Model
Projection of Soil Organic Carbon Stores Following Land Use Changes. Journal of
Environmental management, 45 : 287'302.
Howarth P.J. and G.M. Wickware (1931). Procedures for Change Detection Using
Landsat Digital Data. Inte rnational Journal of Remote Sensing, 2: 27 7 -29 I .
Hubefy C.J. (1994). Applied Discriminant Analysis. John Wiley and Sons, Inc.
New York.
Development in the ESCAP Region, Survey Manuals, New York: United Nation,
pp.l89-215.
Irwin E. G., Jacqueline Geoghegan (2001). Theory, datA and methods: developing
spatially explicit economic models of land use change. Agriculture Ecosystems and
Environment, 85:7-23.
25s
Jadhav R.N., M.M. Kimhoti and A.K.Kandya (1993). Grassland
Jensen J.R., E.W. Ramsay, H.E.Mackey, E.J. Christensen and R.R. Sharitz (1987).
Jensen J.R., K. Rutchey, M.S. Koch and S. Narumalani (1995). Inland Wetland
Change Detection in the Everglades Water Conservation Area 2A Using a Time
Series of Normalized Remotely Sensed Data. Photogrametric Engineering and
Remote Sensing, 6l : 199 -209.
Jha C.S. and N.V.M.Uni (1994). Digital Change Detection of Forest Conversion of
Dry Tropical Forest Region. Internationol Journql of Remote Sensing, 15:2543-
2552.
Joyce A.T., kvey,J.H,and Bums,G.S.(1980). The use of Landsat MSS data for
detecting Land use change in forestland. Proceedings of the l54th International
Symposium on Remote Sensing of Environment. Held in Costa Rica in 1982, pp.
979-988.
2s6
Ju Junchang, eric D.Kolackzyk, Sucharita Gopal (2003). Gaussian mixture
discriminant analysis and sub-pixel land cover charactenzation in remote sensing,
Remote Sensing of Environment. 84 :5 50 :5 60.
Kaufrnann R.K., Karen C.S. (2001). Change Detection, accufacy' and bias in a
Kwarteng A.Y and P.S. Chavez, Jr. (1998). Change Detection Study of Kuwait
City and environs Using Multi-temporal Landsat Thematic Mapper dala.
International Journal of Remote Sensing, 19:1651'1662.
257
Lachenbruch, P.A. (1973). Some Results on the multþle group discriminant
problem, in T.Cacaullus (eds), 1973, Discriminnat Analysis and Application,
Academic Press. New York.
Lambin E.F. and A.H. Strahler (1994). Change Vector Analysis in Multitemporal
Space: a Tool to Detec and Categorise Land Cover Change Processing Using High
Lindgren, David.T (1985). Land use planning and remote sensing. Martinus
Nij hoff Publishers, Dordrecht, The Netherlands.
Liu X.H, Claes A. (2003). Assessing the impact of temporal dynamics on land -use
change modeling. Computers, Environment and Urban Systems, A¡ticle in Press.
Li Xiubin and Hong Yang (2000). Cultivated land and food supply in China. Land
Use Policy,lT:73-88.
Temporal Images for the Monitoring of Rapid Urban Expansion in the Pearl River
Delta. International Journal of Remote Sensing, 19:1501-1518.
258
Li Xia, and Anthony.G.O.Yeh (2001). Calibration of cellular automata by using
neural networks for the simulation of complex urban systems. Environment and
P I annin g. A, 23 :l 445 -l 462.
Lillesand T.M. and R.W. Kiefer (19S7). Remote Sensing and Image Interpretation.
Wiley, USA.
Lillesand T.M. and R.W. Kiefer (2000). Remote Sensing and Image Interpretation.
Wiley, USA.
Logsdon.M.G, Earl.J.B and Frank. V.W. (1997). Probability Mapping of Land Use
Change: A GIS lnterface for Visualizing Transition Probabilities. Computers,
Lopez E., Gerardo Bocco, Manuel Mendoza, Emilio Duhau (2001). Predicting
land-cover and land-use change in the urban fringe, a case in Morelia city, Mexico.
Landscape and Urban Planning, 55:271-285
259
Lyon J.G., D. Yuan, R.S. Lunetta and c.D. Evidge (1998). Change Detection
Experiment Using Vegetation Indices. Photogrametric Engineering and Remote
Sensing,64:143-150.
Malila V/.4. (1930). Change vector analysis: and approach for detecting forest
changes with Landsat.Proceeding of the 6th Annual Symposium on Machine
Processing of Remotely Sensed Dafa. Purdue University,pp.326-335.
neural networks and loistic regression for predicting species distributions: a case
study with a Himalayan river bird, Ecological Modelling.l2}:331-347 .
Matsura S., Abduratrim.A., Mochtar. J. (1997). The Malaysia Land use and Land
cover change case study. SARCS-START-IGBP, Bangkok, Thailand.
Meyer V/. B. and B.L. Tumer II. (1992). Human Population Growth and Global
Land-Use / Cover Change. Annual. Rev. Ecol. System,23 :39-61.
260
Meyer W. B., and Tumer II, B. L. (eds). (1994). Change in Land Us and Land
Cover: A Globlal Perspectiv¿, Cambridge University Press, UK
McConnell W.J and Emilio F.M (eds.). (2000). Meeting in the Middle: The
challenge of meso-level integration. LUCC Report Series No.5. LUCC Focus 1
Office.
McNeil (1994. Toward a Typology and Regionalization of Land - Cover and Land
'W.
use Change: Report of Working Group B. in Meyer 8., and Tumer II, B. L.
(eds). (1994). Change in Land Us and Land Cover: A Globlal Perspective,
Cambridge University Press, UK
Morisette, J.T., Khonam, S., Maco, T. 1999, Land cover change detection
20(t4),2703 -2721.
Moser S.C. (1996). A Partial Instructional Module on Global and Regional Land
use/cover change: Assessing the Data and searching for General Relationship.
Geoj ournal, 39 (3). 241 -250.
Muchoney D.M. and B.N. Haack (1994). Change detection for Monitoring Forest
Deforestation. Photogrametric Engineering and Remote Sensing,60:1243-1251-
261
Muniyati (2000). Wetland change detection on the Kafue, Zambia, by
Classification of a multitemporal Remote Sensing Image dataset. International
Nelson R.F. (1933). Dekcting Forest Canopy Change Due to Insect Activity Using
Landsat MS S. P hoto gr ametric Engine ering and Remote Sensing, 49 :7303 -13 I 4'
Reppprot (19S9). Review of Phase I Results Java and Bali. Direktorat Bina
OECD. (1976). Land use policies and agriculture. Organization for Economic Co-
operation and Development, Paris.
Oetter D.R., Wanen B.C, Mercedes B, Thomas K.M and Robert E.K. (2001). Land
cover Mapping in an Agricultural setting using Multi-seasonal Thematic Mapper
data. Remote Sensing of Environment,76Q): I39-I55.
262
J.R (eds.). W'atershed management: balancing sustainability and environmental
change. Springer-V erlag, New York.
Pathan, S. K., 7 others, et al. 1993, Urban growth trend analysis using GIS
techniques a case study of the Bombay Metropolitan, International Journal of
Remote Sensing, 14(17), 3169 - 3179.
Pilon P.G., P.J. Howarth, R.A. Bullock and P.O. Adeniyi (1988). An Enhanced
Classification Approach to Change Detection in Semi - Arid Environment.
Photo grametri c Engineering and Remote Sensing, 54 : 77 09 - 17 16.
Pontius R.G. Jr., Joseph D. Cornell, Charles A.S. Hall (2001). Modeling the spatial
pattem of land-use change \¡rith GEOMOD2: application and validation for Costa
Rjca. Agriculture, Ecosystems and Environment, 85: 191-303.
Qi J., A.R. Huete, M.S. Moran, A. Chehbouni and R.D. Jackson (1993).
Interpretation of Vegetation Indices Derrived from Multi-Temporal SPOT Images.
Remote Sensing of Environment,44: 89-101.
Quarmby N.A. and J.L. Cushnie (1939). Monitoring Urban Land Cover Changes at
the Urban Fringe from Spot HRV Imagery in South - East England. International
Journal of Remote Sensing, l0: 953 -963.
Rais Jacub. (1997). Planning sustainable land management at regional level: the
Indonesian case, ITC Journal. 3:4.
t4(17),319t- 3200.
263
Rao, C.R.(1966). Discriminant function between composite hypotheses and related
probleq Biometrikø, 53, 339-345.
Rao K.S., Rekfia Pant (2001). Land use dynamic a Landscape change pattem in a
typical micro watershed in the mid elevation zone of central Himalay4 India.
Agricultur e, Eco sy stems and Environment, 86 :l 13 -123 .
Rhind David, Ray Hudson (1980). Land Use. Mathuen, New York.
Richard J.A. (1984). Thematic Mapping from Multitemporal Image Data Using the
Principal Components Transformation. Remote Sensing of Environment,16 3546.
Rindfuss, R.R. and Stern, P.C. 1998, Lingking remote sensing and social science:
the need and the challenges, ln People and Pixels, Eds. Diana Liverman,
Emilio,F.Moran., Ronald,R.Rendfuss. And Paul,C.Stem., National Academy Press
,Washington.
Reid R.S., Kruska R.L., Muthui N., Taye 4., Wotton S., Wilson C.J., Mulatu W.
(2000) Land use and land cover dynamics in response to changes in climatic,
biological and socio-political forces: The case of south western Ethiopia.
Landscape Ecologt, I 5 : 339-355.
Ridd M.K and Jiajun Liu (1998). A Comparison of Four Algorithms for Change
Detection in an Urban Environment. Remote Sensing of Environment,63:95-100.
264
Robinson C.J. (1979). A critical review of the change detection and urban
classification literature. Technical Memorandum CSCÆM-7916235, Computer
Sciences Corporation, Silver Spring, MD.
Robinson J, Stephen Brush, Ian Douglas, T.E. Gradael, Dean Cttactz, Winifred
Hodge, Diana Liverman, Jerry Melillo, Richard Moss, Alexei Naumov, George
Njirv, Joyce Penner, Peter Rogers, Vemon Ruffan, James Sturdevant.(l99a). Land
Use and Land Cover Projections: Report working groups C.in Changes in Land
Use and Land Cover: A Global Perspective.
Sage Colin.1994. Population and hcome in Mayer and Turner II (eds). Change in
Land use ønd land cover: A Global Perspective, Cambridge, UK: Cambridge
University Press: 263 -285 .
Indonesi4 Jakarta.
265
Air Photos and Landsat TM. International Journql of Remote Sensing, 20:1169-
1181.
SARCS (1997). The Finat Report of LULC Case Study in the four countries:
lndonesia, Malaysia, Philippine and Thailand. SARCS Secretary Bangkok.
Thailand.
Scholten H.J, Rob van de velde, Piet Rietveld and Maarten Hilferink. (1999).
Spatial Information Infrastructure for Scenario Planning: The Development of a
Land Use Planer for Holland, in Srillwell.J, Stan Geertman and Stan Openshaw
(eds.). Geo graphical Information and Planning. Springer-Verlag, B erlin.
Sembaloni .F. (1997). An Urban and regional model based on cellular automata.
Environment and Planning B: Planning Design,24:589-612.
Serneels S, Eric F.Larrbin (2001). Proximate cause of land use change in Narok
District, Kenya: as spatial statistical model. Agriculture, Ecosystems and
Environmenf. 85:65-81
266
Singh A. and A.Ha:rison (19S5). Standardized Principal Component.International
Journal of Remote Sensing, 6: 883-896.
Solecki W.D. (2001). The role of Global-to local linkages in land use/land cover
change in south Florida.Ec olo gic al Economics.3T :339 -3 5 6.
KendallÆItmt. USA.
Song C., Cutis E.S/, Karen C.S, Mary P.L and Scott A.M. (2001). Classification
and Change Detection Using Landsat TM Data: When and How to Correct
Atmo spheric Effects? Remote Sensing of Environment, 75 :230'244.
Srillwell. J, Stan Geertman and Stan Geertman and Stan Openshaw (eds.) (1999).
Geographical Information and Planning. Springer-V erlag, Berlin.
267
Schmidt, M. G., Schreier, H. E., Shalì, P. B. 1995, A GIS evaluation of land use
dynamics and forest soit fertility in a watershed in Nepal, International Journal of
Geographical Information System,9 (3), 317 -327.
Stefanov W.L., Michael S.R and Philip R.C. (2001). Monitoring Urban Land cover
Change. An expert system approach to land cover classification of semiarid to arid
urban centers. Remote Sensing of Environment,TT:173'185
Stamp, L.D. (1931). The Land utilization survey in Britain. Methuen & Co Ltd.
London.
Stow D.4., L.R Tinney and J.E.Estes (1980). Deriving Land use/Land Cover
Change Statistics From Landsat: a Study of Prime Agricultural Land. Proceedings
of the l4h International Symposium on remote Sensing of Environment, Ann Arbor
Michigan, ppl227 -1237 .
Sudradjat Adjat, Apik Karya, Atmen. Z, Dedth. F.R, Diana.H, Ekowati Ch.,
Ening.V/., Gatot.P., Giyati, Gunawan.S., Mulyono, Prijanto.P., Slamet.S','Wastra,
268
S.Y. Randha and Triyono (2001). Dasar Penerapan Kebijakan Proyek Lahan
Gambut Sejuta Hektar. Program Pasca sarjana IPB, Bogor.
91.
Tandy Clifford R.V. (1973). Industrial Land use and Dereliction, in Land use and
landscape Planning, Lovejoy Derek (eds.). Leonard Hill Books. Great Britain.
269
environment and peoples of the Kalahari, Botswana, Applied Geography, l7(l),
tt-27.
Todd W.J. (1977). Urban and Regional Land use change detected by using Landsat
datz. Journal of Research by the US. Geological Survey, 5:527'534.
Tumer B.L., Hanham, R.O. and Portaharo, A.V. (1977). Population pressure and
agricultural intensity, Annals of the Association of American Geographers, 6(3),
384-396.
Tu Qiang, Dafu Ding (2003). Detecting pathogenicity islands and anomalous gene
clusters by interactive discriminant analysi s, FE MS Micr obiolo gt Le tter s, 221 :269 -
275.
Tumer B.L, and Meyer. W.B. (1994)."Global land use and land-cover change: an
overview," in W. Meyer and B.L. Turner (eds), Changes in Land Use and Land
Cover: A Global Perspective, Cambridge, UK: Cambridge University Press, pp 3-
10.
Tumer, B.L and W.B. Meyer.199I. Land use and land cover in global environment
change: Considerations for Study, International Social Science Journal, 130:669-
679.
270
Tucker C.J. (1979). Red and Photographic Infrared Linier Combinations for
Monitoring Vegetation . Remote Sensing of Environment, I :127 -150 .
Túcker C.J, and Maxwell, E.L (1976). Sensor design for Monitoring Vegetation
Canopies. P hoto grametric Engine ering and Remote Sensing, 42: 1399 .
van Deusen P.C. (1994). Conecting Bias in Change Estimates from Thematic
Maps. Remote Sens ing of Envir onment, 50 :67'7 3 .
Veldkamp 4., E.F. Lambin (2001). Predicting Land use Change. Agriculture,
Ecosystem and Environment, 85 :1 -6.
271
Ward D., Stuart R.P. and Alan T.M. (2000). Monitoring Growth in Rapidly
Urbanizing Areas Using Remotely Sensed Data. Professional Geographer,52(3):
371-386.
rùVeismiller R.4., S.J. Kristot D.K.Scholz, P.EEE. Anuta and S.A. Momin (1977).
Change Detection in Coastal Zone Environments. Photogrametric Engineering and
Remote Sensing, 43 : 7 533 -1539.
V/eng Q. (2002). Land use change analysis in the Zhujiarg Delta of China using
satellite remote sensing, GIS and stochastic modelling. Journal of Environmental
Management, 64 :27 3 -283 .
Analysis.htþ I I encae.neUff/tamu/STEPWls.htm
White R., Engelen G. (1993). Cellular automata and fractal form: a cellular
modelling approach to the evolution of urban land use pattems. Environment and
Planning A, 25 : ll7 5 -1199.
272
White R, and G.Engele.I.U. (1997). The use of constrained cellular automata for
igh-resolution modelling of urban land-use dynamics. Environment and Planning
B : Planning Design, 24:323 -347 .
Systems,24:383400.
Development ofForest, Land, and water. The World Bank Washington D.C.
History in Jolly Carole.L and Barbara Boyle Torrey (eds) Population and Land Use
in Developing Countries. National Academy Press. Washington, D.C.
273
Wu F and C.J.Webster (1998). Simulation of Land development through the
integration of cellular automata and multi criteria evaluation. Environment and
Planning B : Planning Design, 25 : 103 -126.
V/u Fulong. (1998). Simland: a prototype to simulate land conversion through the
integrated GIS and cellular automata with AHP-derived transition rules.
Yeh A.G.O, ffid Xia Li (1993). Sustainable land development model for rapid
growth areas using GlS. Int.J.Geographical Information Science, 12 Q):169-189.
Yeh A.G.O, and Xia Li (2000). A cellular automata model to simulate development
density for urban planning Environment and Planning B: Planning Design, 29
431-450.
Yeh A.G.O, and Xia Li (2001). A constrained cellular automata model for the
simulation and planning of suitable urban form by using GIS. Environment and
Planníng B : Planning Design, 28: 733-7 53.
274
Young S.S, Benyamin,P., Jokisch,B., Ogneva,Y., and Grren,A.(I991). Global Land
use/cover: Assessment of data and some general relationships. Apendixt I of The
Report of the Ad-hoc Working group on Global Land Use Change, Committee for
Research on Global Environmental Change, Social Science Research Counsil.
NOAA GVI data for vegetation change studies over China. Photogrametric
Engineering qnd Remote Sensing, 65 :679'688.
Young S.S. and C.Y.Wang (2001). Land cover change analysis of China using
global-scale Pathfinder AVHRR Land cover (PAL) data,1982-1992. International
Journal of Remote Sensing, 2l: 1457 -1477 .
Yuan D, and Chris Elvidge (1993). NACL Land Cover Change Detection Pilot
Study: V/ashington D.C. Area Experiments Remote Sensing of Environment, 66:.
166-18.
Yue T. X., Ji Yuan Liu, Sven Erik Jorgensen, Qin Hua Ye (2003). Landscape
change detection of the newly created wetland in Yellow River Delta. Ecological
Modelling,3225: l-ll.
275
APENDIX IV
Tablel com nents transformation TM 1990
Gomp 1 Comp 2 Comp 3 Gomp 4 Gomp 5 Gomp 6
Tova¡. 89.12 7.83 2.69 0.17 0.13 0.07
Eisen value. 1895.66 166.53 57.32 3.52 2.68 1.40
Types\Slope
Forest l) 330( 1962 1701 2378 1263 62C
Types\Slope
Forest 1) 9422 2570 1811 2486 r072 418
2) 1.8i 0.51 0.36 0.49 0.21 0.08
Plantation l) 45529 3693 2078 981 156 4C
Types
Forest 42 3431 603',l t722
Plantation 985? r7135 t3369 230r6
Rice field 214006 19578 12762 s2t99
)oen/drv land 474s8 1r782 3r21 20523
Jettlement 11346 926 635 9064
20976 23 5 44
Types
Forest IT32 7216 8798 636
Plantation 9669 t6448 8722 17638
Rice field z0n09 17348 14372 50165
Open/drv land 47090 14988 3439 18923
Settlernent t9289 1425 589 1903 I
Fishpond 20791 68 32 159
Table V.3: Land use/Land cover change versus Land system
Land use types 1989 (ha) [qqz (ha) Chanee (ha) P/o change lfi.and use types I
1989 (ha) 1997 (ha
Forest sl 47s 4791 o.o9lþorest 0
land 1 0. land I I
38 0. 0 I 23 t
9 0. ish
Land use types 1989 (ha) 1997 0¡,a\ 3hanee (ha) use rqsXha) ltggt $
Forest 0 0 0 0. orest 2os5l
Plantation 0 29 29 0.01 26ß41
Rice field 2681 t3t4 -1361 -0.27lh.ice field
Open/dry land 93( 2223 1288 0.25llopen/dry land
,1
Settlement 69 62 0.0lllSettlement
ll -11 0
orest 5 0. l3
-0.1 28
LULC 1989 (ha) 1997 (ha) 1989 (ha) 1997 (ha) 1989 (ha) 1997 (ha)
Forest 770C 9961 933 1074 c c
Lrlu'
+
High
Lou'
+
High
+
High
Lon
+
Higli
+
High
Cells in transition to
LULC Class 1 Class 2 Class 3 Class 4 Class 5 Class 6 Glass 7 Total 1989
Class I 57150 22326 47773 '1092 116 0 0 12845
Class 2 104503 234134 297490 56736 19007 0 0 71187
Class 3 30538 287651 20651 13 6101 10 208186 82092 0 328369
Class 4 4633 80063 405784 360099 76079 0 0 92665
Class 5 0 1377 431 39 37605 168248 0 0 25036
Class 6 0 0 79397 4406 4453 146093 0 23434
Glass 7 0 0 251 581 0 501 59804 6113
fotal
expected '196824 625551 2938947 1070629 476088 228687 59804 559653
¡ability of :hange to
Given Probabil
LULC Class 1 Class 2 Class 3 Class 4 Class 5 Class 6 Class 7 Total
Class I 0.4449 0.1738 0.3719 0.0085 0.0009 0 0
Class 2 0.1468 0.3289 o.4179 0.0797 o.0267 0 0
Class 3 0.0093 0.0876 0.6289 0.1858 0.0634 0.025 0
Class 4 0.005 0.0864 0.437S 0.3886 0.0821 0 0
Class 5 0 0.0055 0.1723 0.1502 0.672 0 0
Class 6 0 0 0.3388 0.0188 0.019 0.6234 0
Class 7 0 0 0.0041 0.0095 0 0.0082 0.9782
Total 0.606 o.6822 2.3718 0.8411 0.8641 0.6566 0.9782
Markov transition matrix 8-year interval" 20 yeat to predict
From 1997
Cells in transition to
LULC Class 1
Class 2 Class 3 Class 4 Class 5 Class 6 Class 7 Total 1997
Glass 1 37548 30890 101903 17932 7827 1922 119 19814
LULC Class 1
Class 2 Class 3 Class 4 Class 5 Class 6 Class 7 Total 199
Class 1 148605 29721 9907 9907 0 0 0 19814
Class 2 73554 294217 132397 88265 0 0 0 58843
Class 3 0 0 2028478 624147 4681 10 0 0 3',12073