LAND COVER CHANGE AND HYDROLOGICAL REGIMES IN THE SHIRE RIVER CATCHMENT, MALAWI Lobina Getrude Chozenga Palamuleni Student Number 920418079 A thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy in the Department of Geography, Environmental Management and Energy Studies, University of Johannesburg. Supervisor: Professor Harold John Annegarn Johannesburg, 20 March 2009
191
Embed
LAND COVER CHANGE AND HYDROLOGICAL REGIMES IN THE …
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
LAND COVER CHANGE AND HYDROLOGICAL
REGIMES IN THE SHIRE RIVER CATCHMENT,
MALAWI
Lobina Getrude Chozenga Palamuleni
Student Number 920418079
A thesis submitted in fulfilment of the requirements for the degree of Doctor of
Philosophy in the Department of Geography, Environmental Management and
Energy Studies, University of Johannesburg.
Supervisor: Professor Harold John Annegarn
Johannesburg, 20 March 2009
i
Declaration
I declare that the work contained in this thesis is my own original writing. Sources
referred to in the creation of this work have been appropriately acknowledged by
explicit references or footnotes. Other assistance received has been acknowledged. I
have not knowingly copied or used the words or ideas of others without such
acknowledgement.
Signed ___________________ Date __________________
Sections of this work have been presented at conferences and have been submitted for
Journal publication:
Palamuleni, L. G., T. Landmann and H. J. Annegarn, (2006), Land cover
mapping for the Shire River catchment in Malawi using Landsat
satellite data, Palamuleni, L. G., T. Landmann and H. J. Annegarn,
(2008), Awarded the Best Paper at the 6th African Association of
Remote Sensing of the Environment (AARSE) Conference, 30 October
-2 November 2006, Cairo, Egypt.
Palamuleni, L. G., T. Landmann and H. J. Annegarn, (2007), Mapping rural
savanna woodlands, a comparison of maximum likelihood and fuzzy
classifiers, International Geoscience and Remote Sensing Symposium
(IGARSS’07), 23-27 July 2007, Barcelona, Spain
Palamuleni, L. G., T. Landmann and H. J. Annegarn, (2008), An assessment
of land cover change using multi-temporal Landsat imagery for the
Shire River catchment, Malawi, 7th African Association of Remote
Sensing of the Environment (AARSE) Conference, 27 - 30 October
2008, Accra, Ghana.
Palamuleni, L. G., T. Landmann and H. J. Annegarn, (2008), Processing
changes in land cover using Landsat imagery for the Shire River
catchment, Malawi. Paper submitted to Journal of Applied Earth
Observation and Geoinformation.
Palamuleni, L. G., H. J. Annegarn, T. Landmann amd P. M. Ndomba, (2008),
Application of the AVSWATX Tool for the Shire River Sub-Catchment in
Malawi. Paper submitted to Hydrological Sciences Journal.
ii
Dedication
To my husband, Dr Martin E. Palamuleni, for showing me the path to greatness
and
to my sister, Rhoda, and my lovely daughters Tadala and Tamanda
for encouraging me against all odds
iii
Abstract
Land cover changes associated with growing human populations and expected
changes in climatic conditions are likely to accelerate alterations in hydrological
phenomena and processes on various scales. Subsequently, these changes could
significantly influence the quantity and quality of water resources for both nature and
human society. Documenting the distribution of land cover types within the Shire
River catchment is the foundation for applications in this study of the hydrology of the
Shire catchment.
The aim of this study is to investigate the relationships between the measured land
cover changes and hydrological regimes in the Shire River Catchment in Malawi.
Maps depicting land cover dynamics for 1989 and 2002 were derived from multi-
spectral and multi-temporal Landsat 5 (1989) and Landsat 7 ETM+ (2002) satellite
remote sensing data for this catchment. Other spectral-independent data sets included
the 90-m resolution Shuttle Radar Topographic Mission (SRTM) digital elevation
model (DEM), Geographical Information System (GIS) layers of soils, geology and
archived land cover. Core image-derived data sets such as individual Landsat bands,
Normalized Difference Vegetation Index (NDVI), Principal Components Analysis and
Tasseled Cap transformations were computed. From generated composite images,
land cover classes were identified using a maximum likelihood algorithm. Eight land
cover classes were mapped.
A hierarchical multispectral shape classifier with an object conditional approach
determined by the Food and Agriculture Organisation (FAO) Land Cover
Classification System (LCCS) legend structure was used to map land cover variables.
LCCS was used as a basis for classification to achieve legend harmonization within
Africa and on a global scale. Flexibility of the hierarchical system allowed
incorporation of digital elevation objects, soil and underlying geological features as
well as other available geographical data sets. This approach improved classification
accuracy and can be adopted to discriminate land cover features at several scales,
which are internally relatively homogeneous. In addition to compatibility with the
FAO/LCCS classification system, the derived land cover maps have provided recent
iv
and improved classification accuracy, and added thematic detail compared to the
existing 1992 land cover maps.
Fieldwork was conducted to validate the land cover classes identified during
classification. Accuracy assessment was based on the correlation between ground
reference samples collected during field exercise and the satellite image classification.
The overall mapping accuracy was 87%, with individual classes being mapped at
accuracies of above 77% for both user and producer accuracy. The combination of
Landsat images, vector data and detailed ground truthing information was used
successfully to classify land cover of the Shire River catchment for years 1989 and
2002.
Quantitative changes in the areas of various land cover categories and the direction of
change were determined. Land cover change detection was carried out by Multi-date
visual compositing, followed by Post-classification analysis. For the first step,
degradation of vegetation was chosen as the main indicator of change, while post
classification statistical analysis was employed to determine the specific nature of
changes in each land cover type. Multi-date visual composites were found to detect
areas of change and of no change better than post-classification. Using the post-
classification procedure, areal statistics and direction of change in each land cover
class were derived using a combination of both methods. This activity highlights areas
where there are major changes of land cover (i.e. "hot spots"), both in temporal and
spatial aspects. The study revealed significant changes in magnitude and direction that
have occurred in the catchment between 1989 and 2002, mainly in areas of human
habitation. Trends in land cover change in the upper Shire River catchment depict
land cover transition from woodlands to mostly cultivated/grazing and built-up areas.
Twelve per cent of the total land surface of the study area had been converted to
cultivation/grazing over a 13-year interval.
Positive changes (referring to reforestation of degraded areas) in woody closed areas
especially within the former refugee areas close to the Mozambican border, provides
some evidence of the ecological sustainability of the resource. However, the reversal
of the decreasing trend in woody open and savanna shrubs has raised some questions
regarding the possible continuation of the observed trends in future. As subsistence
farming continues to play a dominant role in land cover conversion, degradation, from
v
evergreen Brachystegia woodlands to more open, dry vegetation, and to grassland
formations, will continue.
Considering the present scale of temporal and spatial change of the land cover in the
area, more continuous and comprehensive land cover change monitoring is required
with multi-spectral and multi-temporal satellite data merging. This study has provided
insights into the kind of landscape transformations that have taken place over 13
years, and will serve as input for the monitoring and proper utilisation of the Shire
River catchment for sustainable socio-economic development and water resources
management.
The land cover mapping derived from satellite images served as input for hydrological
modelling within the Shire River catchment. A GIS interface for SWAT, the ArcView
Soil and Water Assessment Tool eXtendable (AVSWATX) tool was used to model the
hydrology of the Shire catchment. Input variables for AVSWATX included digital
elevation data, soil and land cover grids, and weather data (daily rainfall, temperature,
relative humidity and wind speed). Available catchment streamflow data from 1977 to
1981 (5 years) were used for model calibration, while data from 1984 to 1985 were
used for model validation. The calibration was done at daily time-steps, for which
observed and modelled outputs were compared at Liwonde gauging station, the outlet
point of the catchment. Statistical evaluation of simulated catchment streamflows for
the calibration yielded Nash and Sutcliffe efficiencies (ENS) of 86% and 42% for the
monthly and daily predictions respectively. The Nash and Sutcliffe efficiency for the
validation period were considered acceptable, since the model was capable of
capturing 64% of the variance on monthly, and 42% on daily, observed records.
This validated simulation for 2002 land cover was used as a baseline for scenario
development of three scenarios: (i) continued land cover change at current trends
(business as usual); (ii) accelerated land cover degradation, associated with extensive
deforestation; and (iii) land cover restoration, reflection large-scale land restoration
and reforestation. Average annual and monthly modeled outputs from the alternative
scenarios were compared to the business as usual values to compute percent change in
annual values of surface flow, baseflow, and total channel discharge.
The condition of water resources in the Shire River catchment, Malawi, has been
affected adversely by rapid changes in land cover over the last two decades. A record
vi
of land cover changes that have taken place in the Shire River catchment has been
produced. The study has quantified the relationships that exist between land cover
changes and long-term changes in streamflow yield. A cost-effective set of techniques
has been demonstrated, combining satellite remote sensing for land cover mapping
and hydrological monitoring, which can be used in the formulation of policies for
sustainable land and water resources management in Malawi, and similar
environments elsewhere in Africa.
vii
Acknowledgements
I would like to express my genuine gratitude and appreciation to my thesis supervisor,
Professor Harold Annegarn of University of Johannesburg. His willingness to take me in
despite his busy schedule and not only coming in to scrutinize finished products but also
working through each chapter with me taught me a lot about thinking scientifically. Thanks to
Dr. Tobias Landmann from the University of Würzburg, Remote Sensing Chair, Germany,
who guided me in shaping my proposal into a workable project and for his continued support
during the PhD period. In spite of their tight schedule, they were able to provide guidance,
advice and help in this research project. I also really appreciate the confidence they showed in
me for launching into a research area that was completely new to me. I am indebted to your
understanding and kindness for helping me throughout this study, from beginning to
submission of this thesis.
I would also like to extend my gratitude to the University of Malawi, Chancellor College, for
allowing me study leave to do this PhD degree. In addition, my heartfelt gratitude goes to my
colleagues in the Department of Geography and Earth Sciences for their support.
My special thanks go to Mrs Melanie Kneen, Research Assistant at the University of
Johannesburg, for her major contribution to this research study. She was a great help
throughout, always providing good suggestions, advice, guidance and showing great patience.
I appreciate the journey she took me through TNTmips software applications and for
proofreading my drafts. I am also grateful to Dr P. M. Ndomba from the University of Dar es
Salaam, Tanzania for the very useful help concerning the use of AVSWATX.
I would like to thank Mr David Stevens, from United Nations Office for Outer Space Affairs
(UNOOSA), and US Geological Survey, for providing access to the Landsat TM and ETM+
images used for this study.
I would also like extend my sincere gratitude to Fatima Ferraz, a colleague who first
appreciated my work in remote sensing and water resources management. Through her
inspiration, I was able to participate in the European Space Agency “Tiger Africa Initiative”.
This initiative brought a lot of experience and inspiration in my experience as a beginner in
remote sensing and hydrology.
viii
Fieldwork would have been impossible without the support of Mr Jonathan Gwaligwali, GIS
Technician at the University of Malawi, Chancellor College, for gracefully enduring the
July/August heat associated with overhead sun and bushfires in Malawi.
I shared many good moments with fellow PhD students in the Department of Geography,
Environmental Management and Energy Studies who made my stay at the University of
Johannesburg a memorable period of my life. Thanks to the many colleagues I shared offices
with for the support and friendship: Dr. Patience Gwaze, Julião Cumbane, Matthew Ojelede,
Olusola Ololade, Joseph Kanyanga, Charles Ntui, Micky Josipovic, Philip Goyns and Charles
Paradzayi.
I wish to acknowledge the silent presence of my dad, Donald A. Chozenga for his love and
guidance in life. I wish he were still alive to share this memorable academic achievement. My
deepest gratitude and appreciation go to my dearest mum, Rennie Chozenga. I would like to
thank my lovely sisters Rhoda, Rose, Grace and brothers Dalitso and Justice for always caring
for me.
My deepest thanks go to my husband, Dr Martin Palamuleni, for his great comprehension,
love, and support during this study. To my daughters, Tadala and Tamanda, thank you so
much for enduring my absence.
Special thanks to the Malawian community studying in South Africa, Johannesburg who has
been a special pillar of strength. Thanks for all the encouragement and the prayers. My friends
who have been there for me, with even a number of you physically coming over to check on
me, thank you for your love.
I would also like to appreciate the help from Kate Pendlebury for the detailed editing and
proofreading on the draft thesis.
During my PhD study I received support from a number of people and organizations. I thank
Deutscher Akademischer Austausch Dienst (DAAD) through the African Network of
Scientific and Technological Institutions (ANSTI) for the PhD fellowship and the University
of Johannesburg for support during my stay at the University. Support was also given for
operating expenses and conference travel from the National Research Foundation through a
Focus Area Grant: FA2005040600018 “Sustainability Studies Using GIS and Remote
Sensing” to Prof H Annegarn. I am also grateful to the African Association of Remote
Sensing of the Environment (AARSE) and ITC, Netherlands, for a conference/workshop
fellowship to attend the 6th AARSE conference and the refresher course on Innovative
Applications of Remote Sensing and Geoinformation Sciences for female professionals in
Earth Sciences, in Cairo, Egypt.
Above all, I am grateful to the Lord God Almighty who is my source of strength and
wisdom. “I can do all things through Christ who strengthens me”.
ix
Contents
Declaration........................................................................................................... i
Dedication ........................................................................................................... ii
Abstract.............................................................................................................. iii
Acknowledgements ........................................................................................... vii
Contents ............................................................................................................. ix
List of Figures .................................................................................................... xi
Table of Abbreviations and Acronyms...............................................................xiv
1 Introduction 1
1.1 Water resources and sustainable development in Africa...............................1
1.2 Challenges within the Shire River catchment...............................................3
1.3 Aim and objectives......................................................................................6
1.4 Concepts and definitions..............................................................................9
1.5 Structure of thesis......................................................................................10
2 Land Cover Dynamics in the upper Shire River Catchment 11
2.1 Shire River catchment ...............................................................................11
2.2 Overview of land cover mapping...............................................................17
2.2.1 Classification system ...............................................................................19 2.2.2 Earlier land cover mapping in Malawi.....................................................21
2.3.1 Selection of satellite images.....................................................................24 2.3.2 Image processing .....................................................................................27 2.3.3 Image classification .................................................................................33 2.3.4 Land Cover Classification System...........................................................36
2.4 Results and discussions .............................................................................36
2.4.1 Transformation results .............................................................................37 2.4.2 Land cover maps ......................................................................................40 2.4.3 Description of land cover classes.............................................................43 2.4.4 Distribution of land cover categories .......................................................49 2.4.5 Thematic accuracy assessment.................................................................53
3.3.1 Input for change detection .......................................................................65 3.3.2 Approaches ..............................................................................................65
3.4 Results and discussion...............................................................................67
3.4.1 Image overlay...........................................................................................67 3.4.2 Post classification and land cover change areas.......................................70
4.1.1 Land cover and hydrological processes ...................................................86 4.1.2 Hydrological modelling approaches ........................................................87 4.1.3 Overview of the AVSWATX model ..........................................................90 4.1.4 Application of the AVSWATX model to the Shire River
4.2.1 Data..........................................................................................................97 4.2.2 Model setup............................................................................................114 4.2.3 Modelling the Shire River catchment ....................................................118 4.2.4 Testing effects of land cover change......................................................120 4.2.5 Scenario generation................................................................................120
4.3 Results and discussion.............................................................................124
4.3.1 Hydrological characterization of Shire River catchment using
hydrological variables ............................................................................124 4.3.2 Modelling of the Shire River catchment ................................................126 4.3.3 Scenario outcomes .................................................................................142
Figure 1: Shire River hydrograph at Liwonde, 1948 to 2002 - water year beginning November each year ..................................................................................................5
Figure 2: Total annual rainfall from five rainfall gauging stations within the Shire
River catchment.........................................................................................................6
Figure 3: Location map of Shire River catchment, Malawi....................................................11
Figure 4: Population growth in three districts located within the study area..........................12
Figure 5: Location of the sample sites for primary data collection.........................................30
Figure 6: False colour images - 1989 and 2002 ......................................................................37
Figure 7: Principal component images - 2002 and 1989 ........................................................39
Figure 8: NDVI images - 1989 and 2002................................................................................40
Figure 9: Land cover maps - 1989 and 2002 ..........................................................................41
Figure 30: Spatial distribution of soils within the Shire River catchment ..............................100
Figure 31: Rainfall variability in the Shire River catchment – first pattern (Data from
Department of Meteorology, Malawi)...................................................................103
xii
Figure 32: Rainfall variability in the Shire River catchment – second pattern (Data from the Department of Meteorology, Malawi) ....................................................104
Figure 33: Rainfall variability in the Shire River catchment – third pattern (Data from
the Department of Meteorology, Malawi).............................................................105
Figure 34: Weather stations and river gauging stations in the Shire River catchment ...........106
Figure 35: Time series streamflow for the Shire River Mangochi (inflow) and
Liwonde (outflow) for the period 1976 - 1981: data as received ..........................108
Figure 36: Streamflow data for 1977 - 1981, data as received ...............................................109
Figure 37: Smoothed daily streamflow data from 1977 - 1981 ..............................................111
Figure 45: Comparison of daily catchment streamflows for calibration period 1977 - 1981......................................................................................................................131
Figure 46: Comparison of daily catchment streamflow for validation period: 1984 -
Figure 56: Monthly mean, standard deviations and maxima of daily simulated
catchment streamflows for 2002 land cover scenario 7 and 8 simulations ...........150
Figure 57: Monthly mean, standard deviations and maxima of daily simulated
catchment streamflows for 2002 land cover, scenario 9 and 10 simulations ........151
Figure 58: Monthly mean, standard deviations and maxima of daily simulated
catchment streamflows for 2002 land cover, scenario 11 and 12 simulations ............................................................................................................152
xiii
List of Tables
Table 1: Refugees statistics by December 1992 ....................................................................13
Table 2: Description of Landsat 5 TM ..................................................................................26
Table 3: Description of Landsat ETM+.................................................................................26
Table 5: Correlation matrix for Landsat 5 reflective bands...................................................38
Table 6: Percentage of variance and correlation mapped to each principal
components in study area ........................................................................................38
Table 7: Land cover classes and their definitions..................................................................42
Table 8: Spatial distribution of land cover classes – 1989 and 2002.....................................48
Table 9: Error matrix for land cover classes..........................................................................55
Table 10: Assignment of Change Classification Codes of land cover for the Shire
River Catchment for 1989 and 2002 .......................................................................66
Table 11: Land cover changes of the Shire River catchment during 1989 to 2002.................72
Table 12: Areas changed into cultivated or grazing areas between 1989 and 2002................73
Table 13: Areas changed into grassland areas between 1989 and 2002..................................76
Table 14: Areas changed into savanna shrubs areas between 1989 and 2002.........................77
Table 15: Areas changed into built-up areas between 1989 and 2002 ....................................77
Table 16: Areas changed into woody open areas between 1989 and 2002 .............................81
Table 17: Areas changed into woody closed areas between 1989 and 2002...........................82
Table 18: Examples of large-scale hydrologic model applications .........................................88
Table 19: Data sets and sources for input into the AVSWATX model......................................97
Table 20: Spatial distribution of land cover classes and SWAT land cover class
codes for 1989 and 2002 .........................................................................................99
Table 21: Major soil types of the Shire River catchment and percent area covered .............101
Table 22: Soil parameters required by AVSWATX ................................................................101
Table 23: Weather stations and available data ......................................................................106
Table 24: Daily river flow data..............................................................................................107
Table 25: Characteristics of 2002 land cover data and deforestation scenarios ....................123
Table 26: Characteristics of simulated land cover forestation scenarios...............................124
Table 27: Relative sensitivity values of the optimised parameters........................................126
Table 28: Parameter values calibrated in SWAT using the auto-calibration tool..................127
Table 29: Average annual volumes obtained from calibration for 1977-1981......................128
Table 30: Parameters obtained from annual simulations for 1989 and 2002 land cover ......................................................................................................................134
Table 31: Parameters obtained from daily simulations for 1989 and 2002 land cover .........138
Table 32: Surface run-off simulated from 1989 and 2002 land cover...................................139
Table 33: Baseflow simulated from 1989 and 2002 land cover ............................................140
Table 34: Simulation results from bounding cases scenarios ................................................143
Table 35: Simulation results from land degradation scenarios..............................................144
Table 36: Simulation results obtained from land conservation scenarios .............................146
xiv
Table of Abbreviations and Acronyms AVHRR Advanced Very High Resolution Radiometer
AVSWATX ArcView Soil and Water Assessment Tool eXtendable
CGIAR Consultative Group on International Agricultural Research
DEM Digital Elevation Model
ENSO El Niño Southern Oscillation
FAO/LCCS Food and Agricultural Organisation/Land Cover Classification System
FAO/UNESCO Food and Agricultural Organisation/United Nations Education
Scientific and Cultural Organisation
GIS Geographical Information System
GLCN Global Land Cover Network
GLCF Global Land Cover Facility
GPS Global Positioning System
HRUs Hydrological Response Units
IDA International Development Association
IDP Integrated Development Planning
IGBP International Geosphere Biosphere Programme
IWRM Integrated Water Resource Management
Landsat TM Land satellite Thematic Mapper
Landsat ETM+ Land satellite Enhanced Thematic Mapper Plus
LH-OAT Latin Hypercube One-Factor-At-a-Time
MDGs Millennium Development Goals
MSG METEOSAT Second Generation
MODIS Moderate Resolution Imaging Spectroradiometer
NDBI Normalised Difference Built-up Index
NDVI Normalized Difference Vegetation Index
NOAA National Oceanic and Atmospheric Administration
PR Precipitation Radar
SADC Southern African Development Community
SAR Synthetic Aperture Radar
SCE-UA Shuffled Complex Evolution - University of Arizona
SEVIRI Spinning Enhanced Visible and Infrared Imager
SPOT Satellite Pour l'Observation de la Terre
SWAT Soil and Water Assessment Tool
SWRRB Simulator for Water Resources in Rural Basins
TRMM Tropical Rainfall Measuring Mission
USDA-ARS United States Department of Agriculture - Agricultural Research Service
UNEP United Nations Environmental Programme
UN/FAO United Nations, Food and Agricultural Organization
USDA-SCS United States Department of Agriculture - Soil Conservation Service
USGS United States Geological Survey
UTM Universal Transverse Mercator
WSSD World Summit of Sustainable Development
1
Chapter 1
1 INTRODUCTION
Chapter 1 presents the framework for integrated water resources management while
conceptualising the challenges of water resources in Africa in an era of growing
population and climate change. This research addresses these challenging needs within
the upper Shire River catchment in Malawi. Fundamental issues relating to land cover
change and land surface hydrological response have been summarised. Within this
context, the research hypothesis and objectives are articulated. The chapter ends with
an outline of the thesis structure.
1.1 Water resources and sustainable development in Africa
Water resources are inextricably linked with climate change, population growth and
rainfall variability, so the prospect of global climate change has serious implications for
water resources and regional development in Africa [IPCC, 2001]. Efforts to provide
adequate water resources for Africa will confront several challenges, including population
pressure; problems associated with land use, such as erosion/siltation; and possible
ecological consequences of land-use change on the hydrological cycle [Riebsame et al.,
1994]. Assessment of management issues relating to the distribution and use of water
resources requires an integrated approach. Integrated water resources management is not
only required for analyzing consequences of the adverse natural conditions such as floods
or drought but also to assess possible strategies to make the area less vulnerable to
environmental constraints and changing climate [Tolba, 1982]. In addition, the notion of
water resources management requires the matching of water availability and water use in a
river basin [Terpstra and van Mazijk, 2001]. River basins are the preferred land surface
units for water-related regional scale studies because their drainage areas represent natural
spatial integrators or accumulators of water and associated material transports and thus
allow for the investigation of cumulative effects of human activities on the environment
[Lahmer et al., 2001]. This is the system of integrated river basin management endorsed in
Agenda 21 of Rio, 1992 and echoed at the World Summit of Sustainable Development
(WSSD), Johannesburg, South Africa in 2002.
Water availability is generally a consequence of precipitation and catchment run-off. An
aspect often forgotten in this respect is the impact of new patterns of land use and land
cover on the hydrological availability of water. There are many connections between land
surface characteristics and the water cycle. Firstly, land cover can affect both the degree of
2
infiltration and run-off following precipitation events. Secondly, the degree of vegetation
cover and the albedo (degree of absorption/reflection of sun's rays) of the surface can affect
rates of evaporation, humidity levels and cloud formation [Newson, 1992]. Any change in
land use and land cover will have correlated effects in the hydrological regimes, and
possible impacts on the habitat and ecological communities [Calder, 1992; Lorup et al.,
1998]. In essence, the degree and type of land cover influences the initiation of surface
run-off, the rate of infiltration and consequently the rate of ground water recharge [Calder,
1992].
Recent studies in the Southern African Development Community [SADC, 1995] region
have revealed that climatic and land cover changes threaten to undermine the integrity of
riverine habitats, the availability and quality of water, and agricultural productivity
[Headstreams Project, 2004]. Moreover, there are several indicators of water stress and
scarcity in the SADC region, including the amount of water available per person and the
volume ratio of water withdrawn and potentially available [IPCC, 2001]. This situation has
been attributed to increasing population, which translates into increased demand for water
supply (for domestic, agricultural and industrial use), as well as to climate change. Global
warming would induce changes in precipitation and wind patterns, changes in the
frequency and intensity of storms, ecosystem stress and species loss, reduced availability
of fresh water, and a rising global mean sea level [Ominde and Juma, 1991]. Although the
impacts may not be easily predicted, changes in weather patterns may lead to the
prevalence of severe drought conditions or extreme flood events in the SADC region. The
existence of prolonged drought periods vis-à-vis water scarcity will seriously affect
agricultural production and the socio-economic activities in the region. Malawi is one of
the southern African countries likely to experience absolute water scarcity by 2025 [SADC,
1995], which is a challenge for water resources management to sustain economic
development in the country.
To balance supply and demand for water resources and to reduce negative or undesired
effects for the environment and society, changes of actual land cover have to be studied at
all spatial scales. The land surface provides a critical role in the water cycle as it is the
level at which precipitation is redistributed into evaporation, run-off or soil moisture
storage [Verburg et al., 1999]. Thus, land use and land cover studies should be viewed as
responding to the complex interactions and feedbacks linking social and biophysical
3
processes that occur on the land [Dolman and Verhagen, 2003; Maidment, 1993]. With
increasing human activities vis-à-vis water conflicts, it is important to understand the
interactions between hydrological regimes and associated land use and land cover changes
in catchments [Rockström et al., 2002]. Such an understanding can be achieved by
integrating land use planning and water resources management. Land cover and land cover
change data represents a key variable in the management and understanding of the
environment, as well as driving many environmental models such as hydrological models
within large river basins or even for particular smaller catchments. Therefore, there is a
need to develop proper planning and management approaches within the context of
Integrated Water Resource Management (IWRM). IWRM as defined by the Global Water
Partnership [Global Water Partnership, 2005] is a process that considers the co-ordination
of development and management of water, land and related resources to enhance economic
and social welfare without jeopardising the sustainability of the ecosystem. Thus,
sustainable development of water resources is a key to the maintenance of the natural
ecosystem that supports the well-being of human populations.
1.2 Challenges within the Shire River catchment
The Shire River system, the only outlet of Lake Malawi, is probably the most important
water resource for Malawi. Hydro-electric power plants of about 200 MW generation
output, based on a firm flow of 170 m³ s-1
, have been developed on Shire River providing
98% of electricity produced and used in Malawi [Malawi Government, 2001]. This
electricity is the primary source driving the economic and industrial infrastructure and
services in the country. An estimated 20-25 m³ s-1
of water is abstracted for irrigation in the
Lower Shire valley and government sponsored smallholder schemes. Blantyre City
abstracts 1 m³ s-1
of water for both domestic and industrial use. The Shire River has also
led to the development of fisheries, water-transport and tourism industries. This translates
into an increased demand for water for diverse needs and values. When its supply is
limited in quantity or quality or its distribution is uneven, water can be a source of both
cooperation and contestation among its different users [Mulwafu et al., 2003].
Over the last three decades, the Shire River catchment has undergone considerable changes
in the structure and composition of land use and land cover [Malawi Government, 1998b].
The major driving forces are related to human population increases and rainfall variability.
The national population growth rate has been increasing from 2.7% during the 1977 census
4
to 3.2% in 1998 and is likely to double in the next twenty years [National Statistical Office,
2000]. The population density in Malawi is high – the national average density was 87
people km-² and 171 people km
-² of arable land during the 1977 and 1987 censuses
respectively [National Statistical Office, 1991]. Population density in settlements within
the Shire River catchment was recorded at over 275 people km-² during the 1998 census
[National Statistical Office, 2000].
The high population growth has translated into rapidly increasing demands from land in
terms of food, shelter, energy (fuelwood) and construction materials. Some of the
woodlands are now replaced by agricultural crops, while the grass-covered dambos have
been either overgrazed or cultivated and are left bare. Swamp vegetation has been drained
and cultivated. Studies done in 1967 estimated the woodland cover to be 74% while in
1990 the cover was estimated at 61% [Green and Nanthambwe, 1992]. These observations
suggest that woodland cover has declined by 13% between 1967 and 1990 and between
1981 and 1992 Mwanza district alone experienced 1.8% average annual deforestation rate
[Hudak and Wessman, 2000]. Much of the deforestation has been linked to conversion of
communally owned miombo woodlands into agricultural land [Desanker et al., 1997;
Place and Otsuka, 2001], while high wood demands for energy has exacerbated the
situation. Although land for agricultural production is limited to only 37% of the land area
under rain-fed cultivation at traditional management level, as much as 48% of the land was
found to be under cultivation by 1989/90 growing season [Green and Nanthambwe, 1992].
Inappropriate agricultural practices including overgrazing, mono-cropping, cultivation on
steep slopes and river banks and other marginal areas have degraded land through soil
erosion, reduced water retention and the loss of soil nutrients. Aggravating this situation is
the subsequent decrease in land holding sizes, estimated at 0.5 ha per household [Malawi
Government, 1998b].
Consequently, processes of the land hydrology such as run-off, infiltration,
evapotranspiration and interception have been modified. In most cases, this has resulted in
increased run-off, accelerated soil loss with sedimentation problems leading to reduction of
baseflows and increased incidences of flood disasters during heavy storms [Malawi
Government, 1998b]. In addition, hydropower supplies are threatened by low water flows
and sedimentation, hence power disruptions occur frequently especially in dry years
[Kaluwa et al., 1997]. Aggravating the situation is an increase in demand for water, by
5
different groups with diverse needs and values. Furthermore, it is important to note that the
flow from Lake Malawi into the Shire discontinued for a period of 22 years from 1915 to
1937 [Kidd, 1983] and almost dried up in 1997 [Malawi Government, 2001]. On the one
hand, it is hypothesised that the lack of outflow was due to the vegetation growth and
piling of sediments from the small tributaries near the source, while on the other hand, low
rainfall in the catchment area during the period prior to 1937 is said to be responsible for
the lowering of the lake levels [Sheila, 1995]. However, it is unlikely that sedimentation
would have affected the 1915-1937 occurrences due to low population and associated
agricultural activities and deforestation. A significant decline in the flow of Shire River has
been observed since 1992 (Figure 1) with mean flows as low as 130 m³ s-1
in 1997
compared to 825 m³ s-1
in 1980 and 634 m³ s-1
in 1990 [Malawi Government, 2001].
0
100
200
300
400
500
600
700
800
900
1948/1949
1951/1952
1954/1955
1957/1958
1960/1961
1963/1964
1966/1967
1969/1970
1972/1973
1975/1976
1978/1979
1981/1982
1984/1985
1987/1988
1990/1991
1993/1994
1996/1997
1999/2000
Streamflow (m3 s-1)
Figure 1: Shire River hydrograph at Liwonde, 1948 to 2002 - water year beginning November each year
In addition, recent droughts in southern Africa have been associated with the drop in river
flows [SADC, 1995]. In the last three decades, Malawi has experienced variability and
unpredictability in seasonal rainfall. There have been three significant droughts (in
6
1978/79, 1981/82, and most severe was in the 1991/92 season), frequent and increasingly
long dry spells, and an erratic onset and cessation of rainfall [Malawi Government, 2001].
Further discussion on rainfall variability within the Shire River catchment is in section
4.2.1 of this thesis. Rainfall data collected from gauges within the catchment (from 1977 to
1981) is plotted in Figure 2. This is the period that has been utilised for the hydrological
modelling in this study.
1977 1978 1979 1980 1981
500
1000
1500
2000
2500
Total annual rainfall (mm)
Ntaja
Salima
Mangochi
Balaka
Chancellor college
Figure 2: Total annual rainfall from five rainfall gauging stations within the Shire River catchment
Unprecedented rainfall variability means unexpected droughts or flooding which may in
turn produce changes in land use and land cover [Meyer and Turner, 1994]. Precipitation
variability, water scarcity and changes to the water regimes through land cover change will
seriously affect agricultural production and the socio-economic activities in the country. In
addition, such variations influence the temporal phenology and chlorophyll characteristics
of the vegetation in the area, which is a challenge in remote sensing studies.
1.3 Research question, Hypothesis, Aim and Objectives
Research question
This study is intended to address the following question: What are the effects of significant
land cover changes over the past two decades on river flow characteristics that are
important for water resources, environmental functioning and hydrological processes
7
within the upper Shire River catchment? Land cover changes associated with growing
human populations and expected changes in climatic conditions are likely to accelerate
alterations in hydrological phenomena and processes on various scales. Subsequently,
these changes could significantly influence the quantity and quality of water resources for
both nature and human society. This aim will be pursued in the context of developing
integrated land use planning and water resources management in Malawi.
Hypothesis
The following general hypothesis is proposed for the Shire River catchment in Malawi:
Unsustainable changes in land cover due to human activities are significantly
altering aggregate catchment conditions, giving rise to long-term, potentially
irreversible changes in river flow characteristics.
Aim and objectives
This research hypothesis will be tested through a structured sequence of land cover change
analyses and hydrological model simulations. Accordingly, the objectives are set out as:
• To map land cover within the upper Shire River catchment for 1989 and
2002, using Landsat TM and ETM satellite imagery;
• To quantify land cover changes in the catchment between 1989 and 2002;
• To model the hydrological regimes in the upper Shire River catchment and
its sub-basins;
• To challenge the thesis hypothesis by using the hydrological model to
evaluate effects of derived quantitative land cover changes on hydrological
processes;
• To simulate likely changes to hydrological processes in response to
continued land cover changes; and
• To discuss the implications of land use management on stabilising water
regimes of the Shire River catchment.
1.4 Research Design
Land cover mapping and change detection will be based on analyses of two Landsat
images captured 13 years apart. Supplementary digital mapping data sets were obtained
from the Department of Surveys in Malawi. To map land cover dynamics, pixel based
8
classification was undertaken using Maximum Likelihood algorithm. Accuracy assessment
was carried out using producer and user accuracies for each class along with overall
accuracies [Congalton and Green, 1999]. The UN Food and Agricultural Organisation
Land Cover Classification System [Food and Agriculture Organisation, 2005] was used to
label land cover variables to achieve legend harmonisation within Africa and on a global
scale. The new classification is also internally consistent, allowing for scalability and
flexibility that can be used at different scales and different levels of detail to distinguish
land cover features.
Two approaches were utilised to detect and compare changes in the upper Shire River
catchment, namely multi-date visual composite and post-classification analysis. In the
multi-date approach, vegetation reduction was chosen as the main indicator of land cover
change. This technique is not meant to be quantitative, but rather was used to identify and
explore areas of change. Using a post-classification approach, Landsat TM and Landsat
ETM+ images were classified and labelled individually. Later, classification results were
compared on a pixel-by-pixel basis using a change detection matrix where areas of change
were extracted. Quantitative statistics were compiled to determine specific changes
between the two images i.e. magnitude and direction of change in each land cover type
[Calder, 2002].
Hydrological responses were tackled using an existing physically based hydrological
model, the Soil and Water Assessment Tool (SWAT) [Arnold et al., 1994]. This model
incorporates key features of catchment properties, including links between land cover
hydrologic responses. A Geographical Information System (GIS) interface for SWAT, the
ArcView Soil and Water Assessment Tool eXtendable (AVSWATX) tool was used to prepare
parameter input values for the Shire catchment. Input variables for AVSWATX included
digital elevation data, soil and land cover grids and weather data (daily rainfall,
temperature, relative humidity and wind speed). Five rainfall gauges were used to provide
input daily rainfall data to AVSWATX. Available catchment streamflow data from 1977 to
1981 (5 years) was used for model calibration, while data from 1984 to 1985 was used for
model validation. The calibration was done at a daily time-step where observed and
measured outputs were compared at the same outlet point on the catchment, Liwonde
gauging station. Model runs were validated using the parameterised 2002 land cover data.
This validated simulation for land cover in 2002 was used as a baseline for scenario
9
development of three scenarios: (i) continued land cover change at current trends; (ii)
accelerated land cover change associated with extensive deforestation; (iii) reduced land
cover change due to management and reforestation. Average annual outputs from three
alternative futures were then differenced from the baseline values to compute percent
change in annual values of surface flow, baseflow, and total channel discharge.
1.5 Concepts and definitions
This section describes some of the basic terminologies used in land use and land cover
research. The definitions are fundamental to fully understand and apply research results to
a broader readership. The definitions are based mainly from FAO [2005].
Land is any delineable area of the Earths’ terrestrial surface involving all attributes of the
biosphere immediately above and below this surface. It encompasses the near-surface
climate, soils and terrain, surface hydrology and human settlements patterns and physical
results of human activities. Land can be considered in two domains: (i) land in its natural
condition, (ii) land that has been modified by human beings to suit a particular use or a
range of uses.
Land use is the manner in which human beings utilise the land and its resources. Examples
of land use include agriculture, urban development, grazing, logging, and mining.
Land cover describes the physical state of the land surface. Land cover categories include
cropland, forests, wetlands, pasture, roads, and urban areas. Land cover is taken to mean a
physical description of space, of the observed (bio)-physical cover of the Earths’ surface. It
indicates what covers the land such as forest, bushes, uncultivated areas and water bodies.
Land cover classification is the process of defining land cover and land use classes based
on well-defined diagnostic criteria. A classification describes the systematic framework
with the names of the classes and the criteria used to distinguish them, and the relationship
between classes. Such information is taken from ground surveys or through remote
sensing.
Land cover change can be categorised into two types: modification and conversion. Land
cover modifications entail the changes that affect the character of the land without
changing its overall classification and can either be human induced, for example, tree
removal for logging; or have natural origins resulting from, for example, flooding, drought
10
and disease epidemics. Land cover conversion is the complete replacement of one cover
type by another such as deforestation to create cropland or pasture.
1.6 Structure of thesis
This thesis is divided into five chapters. Chapter 1 comprise of the general introduction,
which also outlines the research hypothesis and the objectives. Chapter 2 is a discussion of
land cover dynamics within the upper Shire River catchment based on a supervised
Maximum Likelihood classification of Landsat 5 TM (1989) and Landsat 7 ETM+ (2002).
The variability in spatial land cover extents for each classified land cover class between the
two periods has been examined. The results are used in chapter 3 and 4. Chapter 3 is an
examination of land cover changes within the upper Shire River catchment. The adopted
change detection methods quantitatively reveal the major changes that have occurred in the
catchment between 1989 and 2002. Chapter 4 integrates the work of the previous chapters
(2 and 3) into the preparation of parameters for the Soil Water Assessment Tool (SWAT)
hydrological model. This chapter presents the calibration, validation, and application of
SWAT model for predicting the hydrological response from land cover activities within the
catchment. Critical land cover change simulations demonstrate the capability of the model
in guiding spatially distributed land cover change and precipitation events. The final
Chapter involves the critical discussion of the main research findings and recommends
future investigations to advance the field of physically based hydrological modelling for
the management of water resources.
11
Chapter 2
2 LAND COVER DYNAMICS IN THE UPPER SHIRE RIVER
CATCHMENT
Chapter 2 provides a description of the study area, including topographical, climatic
and hydrological characteristics. An overview of land cover mapping and concepts
with regard to application of satellite data are discussed. This is followed by analyses
of land cover dynamics within the upper Shire River catchment based on a supervised
Maximum Likelihood classification of two images captured 13 years apart:
Landsat 5 TM (1989) and Landsat 7 ETM+ (2002). Differences in spatial land cover extents for each classified land cover class between the two times are examined.
2.1 Shire River catchment
Location
The Shire River catchment lies in the southern part of the Great East African Rift Valley
system and is the outlet of Lake Malawi. The river flows approximately 400 km from
Mangochi on the southern extremity of Lake Malawi, to Ziu Ziu in Mozambique at the
confluence with the Zambezi River (Figure 3). The catchment area of the basin is
18,000 km2 and is divided into upper, middle and lower sections.
Salima
Chancellor College
Balaka
Mangochi
Ntaja
Liwonde
Lake
Malawi
Lake
Malombe
Lilongwe
Zomba
Mangochi
Liwonde
Catchment boundary
Country boundary
Lilongwe
Zomba
Mangochi
Liwonde
Lilongwe
Zomba
Mangochi
Liwonde
Catchment boundary
Country boundary
Catchment boundary
Country boundary
Salima
Chancellor College
Balaka
Mangochi
Ntaja
Liwonde
Lake
Malawi
Lake
Malombe
Salima
Chancellor College
Balaka
Mangochi
Ntaja
Liwonde
Lake
Malawi
Lake
Malombe
Salima
Chancellor College
Balaka
Mangochi
Ntaja
Liwonde
Salima
Chancellor College
Balaka
Mangochi
Ntaja
Liwonde
Lake
Malawi
Lake
Malombe
Lilongwe
Zomba
Mangochi
Liwonde
Catchment boundary
Country boundary
Lilongwe
Zomba
Mangochi
Liwonde
Lilongwe
Zomba
Mangochi
Liwonde
Catchment boundary
Country boundary
Catchment boundary
Country boundary
Figure 3: Location map of Shire River catchment, Malawi
12
The upper Shire River catchment is between Mangochi and Matope, with a total channel
bed drop of about 15 m over a distance of 130 km. The focus of this study is the uppermost
reach from Mangochi to Liwonde, which is almost flat at 465 – 600 m above mean sea
level over a distance of 87 km. It forms a catchment area of 4,500 km2, located between
latitudes 14° 20' S; 15° 12' S and longitudes 34° 59' E; 35° 30' E. The river flows through
Lake Malombe, which is 1.8 m below Lake Malawi.
Population
According to administrative boundaries, the study area is located within three districts
namely: Mangochi, Machinga and Balaka. The population within these three districts has
increased from 644 177 in 1977 to 1 011 843 in 1987 and 1 218 177 in 1998 [National
Statistical Office, 2000]. The southern region of Malawi, which forms the catchment area
of the Shire River, has the highest population density ranging between 53 and 275
people km-² of arable land, varying from district to district [National Statistical Office,
2000]. Population increased by 46% from 1977 to 1998 as depicted in Figure 4. Given that
a high proportion of the population is in subsistence agriculture, an increase of population
has serious implication for the degradation of the environment.
0
100
200
300
400
500
600
700
1977
1987
1998
Population ('000)
Mangochi Machinga Balaka
Figure 4: Population growth in three districts located within the study area
13
The increase in population between the 1977 and 1987 census is not only due to natural
population growth, but also due to the influx of refugees from Mozambique. Between 1976
and 1992, Mozambique was ravaged by civil wars such that many of its neighbouring
countries, including Malawi, became home to refugees. The 1992 refugee statistics for the
district within the study area are shown in Table 1.
Table 1: Refugees statistics by December 1992
District Malawians Mozambican Refugees
Mangochi 551 190 46 973
Machinga 577 860 33 300
Balaka 400 057 126 869
Source: Office of the President and Cabinet December 1992.
Rainfall
The dominant climate in Malawi is tropical savanna with distinct dry and wet seasons.
Rainfall is governed by the movement of the Inter-Tropical Convergence Zone (ITCZ) and
other belts of distribution. The onset of rain is not usually predictable but falls between
October/November and ends in April/May of the following year. Almost 90% of rainfall
occurs between December and March and most of the country receives 800 –
1,200 mm a-1
, with some exceptions [Hutcheson, 1998]. Rainfall statistics from stations
within the Shire River catchment receives an average rainfall of 996 mm a-1
, but there are
variations in the amount of rain, its onset, duration and intensity during the wet season.
Further discussion on rainfall variability within the Shire River catchment may be found
section 4.2.1 of this thesis.
Temperature
Temperatures vary with altitude. During the cold season from May to August, the
highlands have mean temperatures of 15 – 18°C and the low-lying rift valley of 20 – 24°C.
The margins of the river have long hot seasons and high humidity, with mean daily
temperatures ranging from 26°C in January to 21°C in July. During the hot season from
September to January, the highest temperatures are recorded in the Shire Valley and along
the lakeshore with the daily average maximum reaching approximately 32°C in October.
The lowest temperatures (26°C) are recorded over high altitude areas particularly the Shire
Highlands.
14
2.1.1 Vegetation
The natural vegetation in Malawi is part of the extensive dry forest Miombo woodland eco-
region, covering most of the southern and eastern parts of Africa [Abbot et al., 1995;
Desanker et al., 1997]. The Miombo woodland classification is characterised by mixed
deciduous woodlands, with dominant species from the family Caesalpinacea –
Brachystegia, Julbernardia and Isoberlinia. The vegetation types vary depending on
altitude, rainfall pattern, soil types and locations, and range from lowland rain forests,
mopane and sub-mopane woodlands, dry evergreen forests, wooded grasslands, wooded
farmlands, and swamp forests [Moyo et al., 1993].
The dominant vegetation in the study area is mopane woodland, which varies in density
from tall open woodland to dense scrub. Mopane woodlands may form pure stands
excluding other species, but are generally associated with several other prominent trees and
shrubs such as Kirkia acuminata, Dalbergia melanoxylon, Adansonia digitata, Combretum
apiculatum, C. imberbe, Acacia nigrescens, Cissus cornifolia, and Commiphora spp. The
herbaceous component of mopane communities differs according to soil conditions and
vegetation structure: dense swards are found beneath gaps in the mopane canopy on
favourable soils, while grasses are almost completely absent in shrubby mopane
communities on mopanosols [Low and Rebelo, 1996; Smith, 1998; White, 1983]. Within
the Shire highlands and the slopes along the catchment, lowland forest and Brachystegia
woodland are found.
Grasses include large tussocks of Festuca costata and Maximuella davyi interspersed with
cushions of Eragrostis volkensii and the Alloeochaete oreogena. Tall grasses are associated
with low altitude woodland, including Hyparrhenia gazensis, Hyparrhenia variabilis,
Hyparrhenia dichroa, Andropogon gayanus, Setaria palustris and Panicum maximum. In
densely settled and cultivated locations, tall reedy grasses are replaced by Urochloa
pullulans and Urochloa mosambicensis. Woodlands are characterised by Sterculia
africana, Colophospermum mopane, Acacia tortilis and Faidherbia albida according to
locality. Acacia woodland provides valuable grazing from pods to supplement grasses in
the dry season. Mature trees may stand within a dense understory, which includes
Commiphora spp., Bauhinia tomentosa, and Popowia obovata. The understory is likely to
be man-induced since lone-standing mature trees are found in open areas of cultivated
land, and in some cases trees are selectively retained by farmers to maturity (e.g.
Faidherbia albida). Base rich soils support Euphorbia ingens and Commiphora thicket,
15
whilst Hyphaene ventricosa, Hyphaene crinita and Borassus aethiopium palms occur
where the water table is high [Low and Rebelo, 1996; Smith, 1998; White, 1983].
Miombo woodlands form an integral part of the livelihood and farming systems of southern
Africa including the Shire catchment [Frost, 1996]. For most rural communities, the
woodlands are a primary source of energy in the form of fuelwood and charcoal and a
NDVI is calculated from a scene by taking the ratio of the difference of the near infrared
and red reflection and the sum of these two bands, as shown in Equation 9 [De Jong,
1994]:
NDVI=( )( )34
34
TMTM
TMTM
+−
(9)
35
where NDVI is the Normalised Difference Vegetation Index; TM 4 is the TM spectral
band 4; and TM 3 is the TM spectral band 3.
Normalised Difference Built-up Index
In monitoring the spatial distribution of built-up areas, the Normalised Difference Built-up
Index (NDBI) [Jensen, 2005] was used to map the possible areas under settlements
(Equation 10).
45
45
TMTM
TMTM
NIRMidIR
NIRMidIRNDBI
+−
= (10)
where NDBI is the Normalised Difference Built-up Index, MidIR is the Mid-Infrared,
NIR is the Near Infrared, TM is the Thermal Band.
After computing the Normalised Difference Built-up Index (Equation 11), a Boolean
image is created by subtracting the Normalised Difference Built-up Index from Normalised
Difference Vegetation Index i.e through differencing two real number images with a value
(0-1).
NDVINDBIupBuilt area −=
(11)
where NDBI is the Normalised Difference Built-up Index and NDVI is the Normalised
Difference Vegetation Index.
Maximum Likelihood Classification
Various methods are used for classifying land cover such as supervised and unsupervised
classifications. Supervised classification uses area statistics based on sample training to
classify an image, whereas unsupervised classifications involve algorithms that examine a
large number of unknown pixels and divide them into a number of classes based on natural
groupings [Lillesand et al., 2004]. In this study, the Maximum Likelihood algorithm was
used for the classification. Maximum Likelihood is a supervised classifier, i.e. the analyst
supervises the classification by identifying representative areas, called training areas.
These areas are then described numerically and presented to the computer algorithm,
which classifies the pixels of the entire scene into the respective spectral class that appears
to be most alike. In a Maximum Likelihood classification, the distribution of the response
pattern of each class is assumed to be normal (Gaussian). The data should include all
spectral variation within each class. In theory, a statistically based algorithm requires a
minimum of n+1 pixels for training in each class, where n is the number of wavelength
36
bands. However, in practice, the use of a minimum of 10 n to 100 n is advised by Lillesand
et al. [2004].
The Maximum Likelihood algorithm consists of two steps:
Estimation of model parameters, which determines the a posteriori probability that a given
pixel belongs to class i, given that the pixel has feature f. This probability is calculated
using Bayes Rule of conditional probability (Equation 12):
( ) ( ) ( )( ) ( )∑
=
j
ipifp
ipifpfip (12)
where p(f|i) is the probability of a pixel having feature f, given that it belongs to class i,
and p(i) is the probability that class i occurs on the image of interest (also referred to as the
a priori probability).
Discriminant functions, which establish the decision rule to classify a pixel is usually set to
be equal to the a posteriori probability for optimal results. It is stated as follows (Equation
13): if a pixel satisfies the equation, then it is assigned to class i.
( ) ( )fDfD ki ≤= (13)
2.3.4 Land Cover Classification System
To achieve land cover harmonisation within Africa and on a global scale, the FAO/LCCS
legend structure Land Cover Classification System (LCCS) [Food and Agriculture
Organisation, 2005] was used as has been discussed in section 2.2.1. All modifications in
the categories were done keeping in view the area under investigation and application of
the derived land cover maps i.e. hydrological modelling.
2.4 Results and discussions
This section describes results obtained from false colour composite and principal
components analysis. Visual interpretation of the images during land cover classification
was improved by means of these transformations. Land cover classification results are
presented with a display of land cover maps generated from the 1989 and 2002 Landsat
images for the Shire River catchment. This is followed by a description of the
characteristics of each classified land cover class based on the LCCS classification system.
Related information portraying spatial extents and distribution of each classified land cover
37
class has been highlighted. The section ends with an interpretation of the accuracy
assessment results.
2.4.1 Transformation results
To aid visual interpretation, visual appearances of the objects in the image were improved
by image enhancement techniques. The goal of image enhancement was to improve visual
interpretability of an image by increasing the apparent distinction between features.
Training areas were derived from the enhanced images in combination with an
interpretation of the indices.
False Colour Composite
A false colour composite image of bands {7, 4, 2} (RGB) was selected as providing the
best visualization, as discussed in Section 2.3.3. Results of this transformation applied to
the 1989 and 2002 Landsat images are presented in Figure 6.
1989
2002 1989
2002
Figure 6: False colour images - 1989 and 2002
This combination distinctively distinguished different types of land cover features:
cultivated/grazing areas are shown in magenta; fresh water bodies emerge in bluish or
black; built-up areas are revealed as shades of pink; and vegetated areas appear in different
shades of green.
38
Principal components
The six standardised principal components (PC1-6) were calculated from the
transformation of the original TM and ETM+ reflective bands 1 – 5, and 7 and the
correlations are shown in (Table 4 and Table 5). The percentage variance that was mapped
to each component after PCA shows that PC1, PC2 and PC3 had higher variance
percentages. These first three components accounted for 98.4% for TM and 98.3% for
ETM of the variability in the data as shown in Table 6. Thus, higher order components (3 –
6) were dropped from further analysis as they constitute noise in the data set.
Table 4: Correlation matrix for Landsat 7 ETM+ reflective bands
b1 b2 b3 b4 b5 b7
b1 1.0 0.9 0.8 0.3 0.7 0.7
b2 1.0 0.9 0.5 0.8 0.8
b3 1.0 0.5 0.9 0.9
b4 1.0 0.9 0.6
b5 1.0 0.9
b7 1.0
Table 5: Correlation matrix for Landsat 5 reflective bands
b1 b2 b3 b4 b5 b7
B1 1.0 0.9 0.8 0.4 0.7 0.8
B2 1.0 0.9 0.6 0.9 0.9
B3 1.0 0.6 0.9 0.9
B4 1.0 0.8 0.6
B5 1.0 0.9
B7 1.0
Table 6: Percentage of variance and correlation mapped to each principal components
in study area
Percentage of variance for ETM+
Percentage of variance for TM
PC1 83.1 86.7
PC2 15.2 11.7
PC3 0.83 0.79
PC4 0.67 0.61
PC5 0.13 0.13
PC6 0.07 0.06
39
Images compiled from the retained principal components are displayed in Figure 7, for
1989 and 2002 respectively. It is apparent that the principal component transformation
displays the land cover classes better than the false colour composite. The PCA images
show more clearly differences between bare soil, green vegetated areas and water bodies.
Water bodies appear as dark blue to black; recently burned grasslands appear in bright
blue; cultivated areas appear as purple; and built-up areas appear in pinkish red colour.
Varying shades of yellow and light green allow discrimination of stages of phenology of
the different types of vegetation.
1989 2002
1989 2002
Figure 7: Principal component images - 2002 and 1989
NDVI
NDVI maps were produced as a measure of biomass distribution over the landscape. The
maps provide an insight of vegetation distribution in the area as shown in Figure 8.
Vegetation indices provide values that are indicative of the spectral reflectance of the
vegetation at a given place. Vegetated areas yielded high values of NDVI because of their
relatively high near infrared reflectance and low visible reflectance. In contrast, water,
clouds and iron roofed built-up areas have larger visible reflectance than near-infrared
reflectance. Thus, these features yielded negative index values. Rock and bare soil areas
have similar reflectance in the two bands (3 and 4) and result in vegetation indices near
40
zero shown as pale brown to brown. Varying shades of green depict different types of
vegetation. For purposes of this study, our focus was on categorical vegetation mapping, as
such NDVI values derived from these images were not used for quantitative analysis.
Results from NDVI was an initial procedure in the classification process enabling land
cover features to be grouped into two broad categories: vegetated and non-vegetated areas.
1989 20021989 2002
Figure 8: NDVI images - 1989 and 2002
2.4.2 Land cover maps
The upper Shire River catchment land cover classification presented here is a result of the
Maximum Likelihood classifier. Eight land cover classes were identified in the upper Shire
River catchment as displayed in Figure 9.
41
19892002
19892002
Figure 9: Land cover maps - 1989 and 2002
The classified land cover categories in the upper Shire River catchment were woody
closed, woody open, savanna shrubs, grasslands, marshes, cultivated or grazing areas
built-up areas and fresh water. Land cover classes and their definitions developed from
FAO/LCCS classification system are shown in Table 7.
42
Table 7: Land cover classes and their definitions
Land Cover LCCS Code LCCS user label
LCCS Own description LCCS Label attributes
Built-up Areas 5003-14 Built-up areas
Includes village settlements (both clustered and scattered) with grass thatched roofing and mud walls
Medium Density Urban Area(s)
Natural Waterbodies
8011-1 Fresh water Lake and river water Deep To Medium Deep Perennial Natural Waterbodies (Flowing)
Cultivated and managed
terrestrial areas
11291-12771-S0305S0403S0
503
Cultivated/ grazing areas
Cultivated areas during the wet season, communal land for grazing (cattle and goats) after harvesting before the next growing season. Land holding sizes >0.5 ha per
household, on average.
Small Sized Field(s) of Rain fed Graminoid Crop(s) (Two Additional Crops) (Two Herbaceous Terrestrial Crops both with Simultaneous Period) Dominant Crop: Cereals - Maize (Zea mays L.) 2nd Crop: Roots & Tubers -Sweet potato (Ipomoea batatas (L.) Lam) 3rd Crop: Pulses & Vegetables - Beans (Vigna spp.)
Forest 20599-13225-L23L8N2N4P1
0
Woody closed
Gazetted as protected areas but prone to encroachment
Broadleaved Evergreen High Trees with Dwarf Shrubs Major Landclass: Sloping Land, Medium-Gradient Escarpm. Zone, Slope Class: Hilly Soils: Soil Surface, Stony (5 - 40%) Altitude: 1000-1500 m
Grasslands 40397-4732 Marshes Mostly drained and cultivated during dry season in the shallow areas.
Closed Medium Tall Grassland On Permanently Flooded Land
Woodland 20808-97744-L15L5N2N1109
P8
Woody open Remnants of the vegetation mostly in cultivated areas available as fruits trees and agroforestry trees (msangu and masau and mango
trees)
Semi-Deciduous (40 - (20-10)%) Woodland with Open Medium to Tall Herbaceous Layer and Sparse Medium High Shrubs Major Landclass: Level Land, Valley Floor, Slope Class: Flat to almost flat Soils: Soil Surface, Subsurface: Ferralsols Altitude: 300-600m
Grasslands 20441-12289-L11L5N2N4P7
Grasslands Prone to annual dry season burning Medium Tall Grassland with Medium High Trees Major Landclass: Level Land, Plain, Slope Class: Flat to Almost Flat Soils: Soil Surface, Stony (5 - 40%), Altitude: 100-300 m
Shrubland 21104-40064-L11L5N2N1109
P8
Savanna shrubs
Dominated by scattered baobab trees (woodlands) and man induced
thicket common in cultivated parklands
Broadleaved Deciduous Medium High Shrubland with Closed Medium to Tall Herbaceous and Medium High Emergents Major Landclass: Level Land, Plain, Slope Class: Flat To Almost Flat Soils: Soil Surface, Subsurface: Ferralsols Altitude: 300-600 m
43
2.4.3 Description of land cover classes
Detailed descriptions of each of the land cover class characteristics follow:
This woodland category covers those communities referred to as Miombo in which
Brachystegia-Jubernalia species are dominant (Figure 10). It is predominantly composed
of mixed evergreen or deciduous forest; and dry evergreen broadleaf types Combretum,
Acacia and Piliostigma tree species. The few emergent shrubs and grass layers are
depressed by the relatively light-crowned trees. From a phenological point of view, this
type of deciduous forest looks evergreen, as the diverse deciduous species do not shed their
leaves at the same time. Woody closed areas are common in the hilly escarpments
dominated by stony lithosols.
Figure 10: Woody closed
Woody open (LCCS code: 20808-97744-L15L5N2N1109P8)
Woody open consists of plant formations comprising a continuous woody stratum and an
herbaceous stratum (Figure 11). The woody trees dominate the herbaceous layer, with a
total woody cover of between 35 and 60 percent. This category is associated with short
grass and deciduous broad-leaved woodlands, common on areas of ferralsols soils.
Selective felling has given rise to various communities with single species dominance of
trees such as Zizyphus jujube (masau), Faidherbia albida (msangu) and Mangifera indica
(mango) common in cultivated parklands. These trees are left intentionally - farmers
collect the fruits for sale or domestic use (masau and mango). Faidherbia albida is rich in
nitrogen and mostly left in the farmlands for nitrogen fixing purposes. This vegetation
exemplifies the regrowth after exploitation by over-logging, clear-cutting or cultivation.
44
woody open with mature trees of F.albida woody open with shrubs and regrowthswoody open with mature trees of F.albida woody open with shrubs and regrowths
are in leaf during the dry season. The plausible way to separate them was using the DEM
where marshes are located in very low areas close to water bodies while woody trees are
mostly located in upper areas in association with cultivated areas.
Shrublands and grasslands in the catchment form a complex matrix that is patchy in some
places and homogeneously mixed in others, and spectral separation was difficult. To
separate these types, transformations and contextual knowledge of the area were
incorporated. Digital elevation and soil data helped to identify boundaries between woody
54
open and woody closed areas. The closed woodlands are dominant in hills and plateau
areas while the open woodlands are common in the lower escarpment. From a remote
sensing point of view, discriminating open mixed deciduous with an area cleared under
shifting cultivation or logging activities proved difficult. Because of the dry non-growing
season, the field supposed to belong to cultivation category appeared mostly as shrubland.
The most accurately classified land cover was the fresh water category due to its
contiguous nature that exhibited a unique spectral signature. Fresh water had 100%
accuracy for both producer and user accuracy. The user accuracy of 100% for grassland
could be explained by the fact that the land cover class was broad and large homogeneous
polygons were purposively selected for the exercise.
2.5 Conclusion
The basis of this research comprised multi-temporal classification of Landsat satellite
imagery to provide a recent perspective of land cover types within the upper Shire River
catchment. Results from this study indicate successfully the synergy between Landsat data,
vector data and detailed ground information in mapping land cover features. It was found
that by integrating contextual information and ancillary data, discrimination of
heterogeneous land cover classes dominant in savanna woodlands was improved. Eight
land cover classes were mapped for the Shire River catchment, with an overall accuracy of
87%. Thus, the classification procedure not only separated land cover into different
classes, but also illustrated the patterns that exist across the landscape.
In addition, this project has produced land cover maps at 1:50 000 scale for the Shire River
catchment, covering 4,500 km2. These maps were prepared using Landsat satellite data,
acquired in 1989 and 2002 as the main data source and thus represent the land cover
existing at that time. The land cover maps have compatible digital formats hence they can
easily be applied to a variety of future GIS applications. Additional themes can be
incorporated as more resource information becomes available, or as new management
needs are identified.
55
Table 9: Error matrix for land cover classes
Ground truth data
Classification Fresh water
Built-up areas
Cultivated /grazing
Marshes Grass-lands
Savanna shrubs
Woody open
Woody closed
TR Accuracy
Fresh water 11 0 0 0 0 0 0 0 11 100.0%
Built-up areas 0 20 5 0 0 0 1 0 26 76.9%
Cultivated /grazing 0 2 40 0 1 2 5 1 51 78.4%
Marshes 0 2 0 23 0 0 0 0 25 92.0%
Grasslands 0 0 0 0 29 0 0 0 100 100.0%
Savanna shrubs 0 0 0 0 0 20 2 0 22 90.9%
Woody open 0 0 1 6 0 0 31 1 39 79.5%
Woody closed 0 0 0 1 0 0 0 24 25 96.0%
TC 11 24 46 30 30 22 39 26 228
Total 100% 83.3% 87% 76.7% 96.7% 90.9% 79.5% 92.3%
where TD = sum of major diagonal, TC = column totals, TR = row totals. The overall classification accuracy was TD/TR (198/228) = 87%. Errors were considered consistent with limits of the available technology and ancillary data.
56
The present study has adopted the hierarchical legend structure determined by the Food
and Agriculture Organisation (FAO) Land Cover Classification System (LCCS) to label
land cover variables [Di Gregorio and Jansen, 2005]. LCCS was used as a basis for the
classification to achieve legend harmonisation within Africa and on a global scale.
Currently, the system serves as an internationally agreed reference base for land cover
mapping. LCCS is a relatively new classification system and has been applied in the
Africover project [Di Gregorio and Jansen, 2005]. Through the Africover Programme and
the GLCN initiative, a number of countries were mapped using the LCCS classification
system including: Burundi, Democratic Republic of Congo, Egypt, Eritrea, Kenya,
Rwanda, Somalia, Sudan, Tanzania and Uganda. No land cover mapping study to date has
been conducted in Malawi using LCCS as a coding classification system. The adoption of
LCCS carried out in this research may thus be an important step towards a rigorous update
and translation of the existing land cover maps for Malawi.
The new classification system is also internally consistent, allowing scalability and
mappability that can be used at different scales and levels of detail to discriminate land
cover features. The methodology is applicable at any scale and is comprehensive in the
sense that any land cover identified anywhere in the world can be readily accommodated.
In addition to compatibility with the FAO/LCCS, the derived land cover maps have
provided recent and improved classification accuracy, added thematic detail compared to
the existing 1992 land cover maps.
LCCS presents an advantage for mapping heterogeneous landscapes such as those present
within the Shire River catchment in Malawi. The distributions of savanna woodlands, rural
residential areas, both clustered and scattered (with grass thatched roofing, mud walls and
occasionally iron sheet roofing) and cultivated or grazing areas, represent classes which
have similar spectral signatures (especially during the dry season). They occur in similar
environments and are often in adjacent or mixed stands. Flexibility of LCCS allows
incorporation of geographical data sets to facilitate definition of class boundaries that are
explicit and clear. In this study, digital elevation objects, soil and underlying geological
features were used to identify boundaries between spectrally similar land cover types, for
example, open woodlands and closed woodlands. The closed woodlands are dominant in
hills and plateau areas associated with lithosols while open woodlands are common in the
lower escarpment. Overall, the hierarchical classification technique presents an advantage
57
for mapping heterogeneous landscapes such as those present within the Shire River
catchment.
In terms of classification methodology, the overall approach developed and applied in the
study comprises a variety of techniques ranging from integration of contextual information
such as terrain height and soil characteristics to Maximum Likelihood classifier. The
combination of analytical tools grants considerable insights in land cover mapping of
heterogeneous savanna woodland environments. The approach has the potential of being
adapted and replicated in other similar regions of the world, hence enhancing knowledge of
land cover mapping and dynamic processes at the policy-relevant meso-scale.
The derived land cover maps have added thematic detail enabling spatial and temporal
analysis of patterns and trends in land cover in the upper Shire River catchment. Mapping
of built-up areas was significant as it provided information for establishing the association
between the spreading pattern of the built-up areas with the intensity of cultivation system
and other changes in land cover classes. As such, the characteristically localised nature of
catchment changes can easily be accounted for, making it possible for interpretation and
prediction of the effects of land cover change. These changes are further discussed in
Chapter 3. Results of land cover dynamics reported here will be used as input data for
hydrological modelling in Chapter 4 to establish relationships between fluctuations in the
hydrological regimes and each land cover class.
58
Chapter 3
3 LAND COVER CHANGE ASSESSMENT 1989-2002
Chapter 3 examines of land cover changes in the upper Shire River catchment. Using
Landsat classifications, a discussion of the changes in landscape diversity and
fragmentation between 1989 and 2002 is presented. Areal statistics and the direction
of change in each land cover class were derived. Combination of both image overlay
and post-classification change detection methods reveals significant changes that have
occurred in the Shire River catchment between 1989 and 2002. Changes in land cover
have potential effects on the hydrological processes of the catchment.
3.1 Land cover change
Rapid changes in land cover have a significant impact on conditions of catchment
ecosystems. Accurate information on the status of and trends in land cover changes is
needed to develop strategies for sustainable development and to improve the livelihood of
communities. The ability to monitor catchment land cover changes is highly desirable by
both local communities and policy decision makers. With increased availability and
improved quality of multi-spatial and multi-temporal remote sensing data as well as new
analytical techniques, it is now possible to monitor catchment land cover changes and their
potential hydrological responses in a timely and cost-effective way.
In Malawi, as in many other countries, the landscape is continually changing under the
influence of several factors (including population growth, agricultural expansion and
urbanisation), and as a result, land cover maps rapidly become out of date. In the area
under consideration - the Shire River catchment - two factors may be related to land cover
changes in the past two decades: population growth and the expansion of subsistence
agriculture. Population growth causes changes in land cover through the expansion of
human settlements. The conversion of land to cropland for food security is a cause of
major concern. Forest clearance has resulted from an increased demand for forest products
such as fuelwood, commercial logging and construction materials. Many current land
cover practices in the catchment have the potential to adversely effect and degrade the
environment with respect to forests, soil, water, and biodiversity resources. Previous
studies show increasing land cover change without quantifying the degree and direction of
change in terms of land cover categories [Green and Nanthambwe, 1992]. The impact of
such land cover changes on land surface hydrological processes is of major interest in this
thesis.
59
3.1.1 Land cover change and hydrological response
Land cover plays a critical role in the hydrological cycle, as water availability is generally
a consequence of the distribution of precipitation into evaporation, run-off and soil
moisture storage [Dolman and Verhagen, 2003]. There are many connections between land
surface characteristics and the water cycle. Firstly, land cover can affect both the degree of
infiltration and run-off following precipitation events. Secondly, the degree of vegetation
cover and the albedo of the surface can affect rates of evaporation, humidity levels and
cloud formation. Any change in land cover will have correlated effects in the hydrological
regimes, and possible impacts on the habitat and ecological communities [Calder, 1992;
Lorup et al., 1998].
Agricultural developments have resulted in the widespread deterioration of soil structure, a
process which favours soil sealing and crusting, and related modifications in the rates of
infiltration and storage [Boardman and Favis-Mortlock, 1993; De Roo et al., 2001;
Robinson, 1990]. Specifically enhanced grazing pressure and intensive cultivation
practices have led to soil compaction, reduced infiltration and ground water recharge and
excessive run-off [Boardman and Favis-Mortlock, 1993; Evans, 1990; Forher et al., 2001;
Moussa et al., 2002].
In catchments which were traditionally cultivated or grazing areas that have been changed
to woodland, increases in evaporative loss and decreases in discharge to the outlet have
been observed [Dagnachew et al., 2003; De Roo et al., 2001; Robinson, 1990]. The effect
of land drainage on river discharge can enhance peak flows by increasing the density of
rivers or inhibiting water flow (infiltration) and storage within the soil matrix [Maidment,
1993]. Consequently, seasonal variations in river discharge have been associated with
decreased infiltration during the wet season, which can influence the streamflow during the
dry season by inhibiting ground water recharge.
The consequences of forestation include an increase in infiltration and a reduction of the
incidence of surface run-off influencing the seasonal regime of rivers. However,
infiltration generally depends upon the combined effect of soil properties, climatic
conditions as well as land cover and the developmental stage of the vegetation.
Demonstrations from catchment studies showed that the pine forestation of former
grassland not only reduces annual streamflow by 440 mm but also reduces the dry season
flow by 15 mm [Bosch, 1979]. Thus, benefits gained by the forestation of degraded or
eroded catchments will be dependent on the situation and the management methods
employed [Calder, 1998].
60
Various microclimate changes resulting from land cover changes have been documented.
Deforestation alters the disposition of radiant energy by increasing surface albedo and
daytime long-wave emission by the land surface resulting in lower net radiation [Bastable
et al., 1993; Culf et al., 1995; Larkin, 2002]. The characteristics and variable state of the
land surface control the proportions of net radiation extended upon latent
evapotranspiration and sensible energy flux. A drastic change in vegetation cover, such as
clear cutting in the Pacific north-west, can produce 90% more run-off than in catchments
unaltered by human practices [Franklin, 1992]. Replacing forests with other land cover
types changes the energy partitioning because of differences in leaf area, aerodynamic
roughness, root depth, and stomata behaviour. Field studies have verified that the
replacement of forest by pasture or non-irrigated crops reduces evapotranspiration which
generally increases streamflow, possibly increasing flood hazards [Giambelluca et al.,
2000; Jipp et al., 1998; Wright, 1992]. However, as Bruinjzeel [2001] emphasises, the
effects of deforestation on evapotranspiration and streamflow are not uniform, and depend
on the original forest type, characteristics of replacement land cover, climate, exposure and
soil depth.
Canopy characteristics such as canopy depth, crop height and spacing, leaf area and shape
will have an influence on the hydrological performance of vegetation cover through its
effects on interception. Studies conducted in wet conditions [Bouten et al., 1991] indicate
that interception losses will be higher from forests than from shorter crops resulting in the
reduction of run-off from forested areas compared with those under shorter vegetation. In
addition, the amount of vegetation controls the partitioning of incoming solar energy into
sensible and latent heat fluxes thereby affecting evapotranspiration rates. Yet in deforested
areas, interception and evaporation rates will be lower while increasing surface run-off and
sediment production which consequently lower flow discharge to rivers [Calder, 1992;
Newson, 1992; Shaw, 1990]. However, the amount of precipitation intercepted and
evaporated will vary depending on the vegetation type and species [Hall and Calder,
1993].
Land cover conversions such as the change from forests or vegetated areas into urban land
uses have many hydrologic effects with significant ecological and sociological
ramifications [Booth, 1991; Hollis, 1975; Leopold, 1968]. These changes reduce
interception, infiltration, subsurface flow, evapotranspiration, storm water storage on hill
61
slopes, and the time required for storm water to travel over and through a hill slope to a
stream [Burges et al., 1998; Dinicola, 1990]. The percentage of the catchment surface that
is impermeable due to urban and road surfaces influences the volume of water that runs
and increases the amount of sediment that can be moved [Arnold and Gibbons, 1996].
Urbanisation is known to transform permeable areas such as forest and farmland into
impermeable areas, thereby reducing the amount of water infiltrated into the soil and the
amount of water on the ground surface [Cheng and Wang, 2002; DeFries R. and
Eshleman, 2004]. As run-off volumes in urban channels increase, the duration of high flow
decreases because groundwater is no longer contributing to the flow. In addition, urban
development causes a decrease in lag time between rainfall and run-off by increasing the
hydraulic efficiency of the drainage system (water can reach the channel more swiftly
when it travels over smooth, hard surfaces). In addition, in urban areas, which are largely
covered by roads and buildings, water evaporation from the ground surface and plants into
the atmosphere is generally lower than in agricultural areas or forests. Recent results
regarding the impact of land cover change on hydrological regimes in urban areas have
been discussed elsewhere [Acreman et al., 2000; De Roo et al., 2001; Jennings and
Jarnagin, 2002; Lahmer et al., 2001]. Theoretically, urban built-up areas - especially those
with impervious surfaces - could affect the flow behaviour within a catchment.
The size of a catchment will also influence the hydrological response of precipitation
where small catchments show obvious responses to specific land uses. In large catchments,
meanwhile, hydrological responses are affected by the complex water storage and release
mechanisms [Calder, 1992; Forher et al., 2001]. In small catchments, land use changes to
urbanisation will reduce evapotranspiration and infiltration and increase surface run-off
[De Roo et al., 2001]. Land use change to forest will increase evaporation which may
result in low soil moisture and possibly increase infiltration during rainfall resulting in a
decreased run-off [Calder, 1992; Dagnachew et al., 2003; De Roo et al., 2001]. However,
management activities associated with forestry such as cultivation, drainage, road
construction and soil compaction during logging are more likely to influence hydrological
responses than the presence or absence of the forests themselves.
The magnitude of land cover effects will be dictated by the nature and scale of the land
cover change. This is also linked to physiographic, climatic and management differences.
Land cover change studies should be viewed as responding to the complex interactions and
62
feedbacks linking social, or indirect and direct biophysical processes that occur on the land
[Dolman and Verhagen, 2003; Maidment, 1993]. To balance the supply and demand for
water resources and to reduce negative or undesired effects for the environment and
society, changes in actual land cover have to be studied at all spatial scales due to the
heterogeneous patterns of the distribution of terrestrial vegetation and soils over land
surfaces [Verburg et al., 1999]. Therefore, establishing clear linkages between specific
land cover changes and hydrological responses at the catchment scale poses a challenging
problem in hydrological studies.
3.2 Land cover change detection
The use of multi-date satellite remote sensing data to detect land cover change began in the
early 1970s [Singh, 1989]. Change detection is the process of identifying differences in the
state of an object or phenomenon by observing it at different times [Singh, 1989]. The
process encompasses the quantification of multi-date imagery to derive changes over two
time periods [Coppin et al., 2004]. It is an important process in monitoring and managing
natural resources and urban development because it provides a quantitative analysis of the
spatial distribution of the population of interest. Change detection is useful in such diverse
applications as land cover change analysis, the monitoring of shifting cultivation, the
assessment of deforestation, the study of changes in vegetation phenology, seasonal
changes in pasture production, damage assessment, crop stress detection, disaster
monitoring, day/night analysis of thermal characteristics as well as other environmental
changes [Singh, 1989].
Techniques to perform change detection with satellite imagery have become numerous as a
result of the increasing versatility of digital data and increasing computing power. A wide
range of approaches to change detection analysis have been reported [Fung, 1990; Green et
al., 1994; Howarth and Wickware, 1981; Jensen, 1986; Jensen et al., 1993; Lillesand et al.,
2004; Milner, 1988; Mouat et al., 1993; Singh, 1989]. Approaches include transparency
compositing [Crapper and Hynson, 1983], image differencing [Dale et al., 1996;
Muchoney and Haack, 1994; Price et al., 1992], image ratio-ing, classification
comparisons [Dale et al., 1996; Jensen, 1995; Muchoney and Haack, 1994], image
enhancement techniques, such as principal component analysis [Fung and LeDrew, 1987],
Normalized Difference Vegetation Index [Mikkola, 1996], image algebra change detection
[Green et al., 1994], and post-classification comparison [Rutchy and Vilchek, 1999].
63
Although these methods have been successful in monitoring change for a myriad of
applications, there is no consensus as to the ‘best’ change detection approach. The type of
change detection method employed is largely dependent on the purpose of the change
investigation, data availability, the geographic area of study, time and computing
constraints, and the type of application [Coppin et al., 2004].
Rather than attempt to review all published change detection techniques, only techniques
that have been selected for use in this study are reviewed, namely image overlay and post-
classification analysis. The purpose of change detection in this exercise is focused on
factors that may influence the behaviour of the hydrological regimes in the catchment.
Image overlay
This change detection algorithm uses a digital enhancement technique for on-screen
change delineation. It provides a simple mechanism to display changes between two dates
of imagery quickly and efficiently [Jensen, 2005]. The image is prepared by making a
photographic two-colour composite showing the two dates in separate colour overlays. The
colours of the resulting image indicate changes in reflectance values between the two
dates. For instance, features which are bright (high reflectance) on date one, but dark (low
reflectance) on date two, will appear in the colour of the first photographic overlay and
vice versa. Features which are unchanged between the two dates will be equally bright in
both overlays and hence will appear as the colour sum of the two overlays.
To represent the bi-temporal variations within one single image product, Alwashe and
Bokhari merged TM bands 2, 4, and 5 of two different acquisition dates via an intensity-
hue-saturation (IHS) transformation. On such an image, vegetation differences showed up
in distinctly different colours [Alwashe and Bokhari, 1993]. Similarly, Sunar [1998]
performed change detection between Landsat TM data of 1984 and 1992 and observed that
areas of no change were represented by values of 127 (mid-grey), while areas that were
darker in 1992 than they were in 1984 had values between 128 and 255. However, this
technique provides the analyst with little information regarding the nature of the change
[Jensen, 2005].
Post-classification analysis
Post classification analysis is the most commonly used method of land cover change
detection. It is a comparative analysis of spectral classifications for times t1 and t2 (or
more) produced independently [Singh, 1989]. Classification results are compared on a
64
pixel-by-pixel basis using a change detection matrix, and the areas of change extracted
[Jensen, 1996; Jensen, 2005; Singh, 1989; Yuan and Elvidge, 1998]. The classified images
are combined to create a new change image classification. The new image indicates the
changes “from” and “to” that took place and the kind of landscape transformations that
have occurred are calculated and mapped. The advantage of this method includes the
detailed “from-to” information that can be extracted. Individual classification of two image
dates minimises the problem of normalizing for atmospheric and sensor differences
between two dates [Jensen, 2005; Singh, 1989] although accurate geo-referencing and
accurate classifications are crucial to ensure precise change-detection results [Augenstein et
al., 1991; Foody, 2001; Rutchy and Vilchek, 1999].
A useful example of the use of the post-classification approach is the work of Munyati
[2000] in which Landsat images acquired in 1984, 1988, 1991 and 1994 were used to
assess change on a section of the Zambian Kafue Flats floodplain wetland system. Similar
supervised Maximum Likelihood classification procedures were employed on all images.
The classified images produced were analysed for change in each land cover category by
overlaying them in a GIS framework, and transition rates between the classes were
calculated. The change detection results provided a reliable indication of the long-term
change that the Kafue Flats wetland area has undergone.
Land cover maps derived from Landsat 5 and Landsat 7 ETM+, discussed in Chapter 2,
form the input data for change detection. Two approaches were used to detect changes in
the upper Shire River catchment, namely image overlay and post-classification analysis
and their results compared. Using image overlay, vegetation characteristics were chosen as
the main indicator of land cover change. Post classification statistical analysis was
employed to determine the specific nature of changes between dates of imagery (1989 and
2002) in each land cover class. Image overlay was found to detect areas of change and no
change within the catchment. By using post classification procedures, areal statistics and
the direction of change in each classified land cover class were derived. A combination of
both methods reveals the significant changes that have occurred in the Shire River
catchment between 1989 and 2002. These activities also highlighted areas where there are
major changes in land cover (i.e. "hot spots"), in both temporal and spatial aspects.
65
3.3 Methodology
3.3.1 Input for change detection
Land cover maps derived and discussed in Chapter 2 form the basis of data for the land
cover change detection.
3.3.2 Approaches
Image overlay
In this approach, vegetation reduction was chosen as the main indicator of land cover
change. Therefore, band 4 of Landsat TM and ETM+ (near-infrared) respectively were
selected. For the 1989 image, band 4 was displayed through the R and B (blue) channels of
the computer monitor. Band 4 of the 2002 image was displayed through the G (green)
channel of the computer monitor. This technique is not intended to be quantitative but
rather was used to identify qualitatively and explore the areas of change in the Shire River
catchment.
Post-classification
The starting point for the post-classification approach is a pair of images for which land
cover classification and labelling have been carried out, in this case the Landsat TM and
Landsat ETM+ images discussed in Chapter 2. The next step is to compare the resultant
classification images on a pixel-by-pixel basis and to extract areas of change using a
change detection matrix. A change detection matrix is an {N x N} table of change
detection classification codes produced by assigning a unique integer to each possible
change from one land cover category to every other category, where N is the number of
thematic land cover classes in the study.
For this study, eight land classes were identified. Each of the eight identified land classes is
assigned a base value in the sequence {20, 40 …160}. A 64-cell table (8 x 8) of change
detection codes, was created for each possible pair of “changed from” and “changed to”
combinations by incrementing a base value by unity, e.g. from the base value 20 assigned
to built-up areas, the vector {21, 22 …28} is created. The full set of assigned codes is
presented in Table 10. A change from savanna shrubs to built-up areas, for example,
would be represented by the code 26 shown in column 1. Diagonal elements of the matrix
are assigned to pixels that have not undergone change. Certain changes are physically
improbable, e.g. open water to woody closed, but values are assigned to complete the
matrix.
66
By applying these change detection codes to the comparison of the two classified images, a
new layer is produced, the post classification change layer, with each pair of pixels
represented by a change detection code. From the change detection layer, change detection
maps may be produced to visualise areas of change by assigning colours to single or
grouped categories of change. For example, any pixel in the 1989 map that changed to
built-up area by 2002 is assigned the value red; any pixel that changed into cultivated or
grazing land by 2002 is assigned the value pink. The assigned colours to be used in the
results section are superimposed on the matrix shown in Table 10.
Table 10: Assignment of Change Classification Codes of land cover for the Shire River
Catchment for 1989 and 2002
Built-up areas
Cultivated/grazin
g
Fresh Water
Grass-lands
Marshes Savanna shrubs
Woody open
Woody closed
To 2002
From 1989 20 40 60 80 100 120 140 160
Built-up areas 21 41 61 81 101 121 141 161
Cultivated/grazing 22 42 62 82 102 122 142 162
Fresh water 23 43 63 83 103 123 143 163
Grasslands 24 44 64 84 104 124 144 164
Marshes 25 45 65 85 105 125 145 165
Savanna shrubs 26 46 66 86 106 126 146 166
Woody open 27 47 67 87 107 127 147 167
Woody closed 28 48 68 88 108 128 148 168
� no change in land cover between dates � new built-up areas � land cover change to cultivated or grazing areas � water related changes not selected for display � water � change in land cover to vegetated areas � vegetation regeneration from previously built-up areas � land cover change to grassland areas � change in land cover to vegetated areas
This table is interpretable by presenting tabular results quantifying change between classes,
such as the change from woody open to cultivated or grazing land. The table shows the
degree of change, the pixel prior classification, and the pixel post classification. This
matrix can be used to distinguish the entire change analysis in one view and is thus a
powerful tool for examining landscape change.
67
Land cover change areal statistics were also summarised to express the net land cover lost
and gained by each class. This was accomplished by subtracting the total amount of land
cover gain from the total amount of loss. Statistics were compiled in the Microsoft Excel
programme to determine the specific nature of changes between dates of imagery, i.e. the
size of the differences between the dates and the direction (positive and negative) in each
land cover type.
3.4 Results and discussion
3.4.1 Image overlay
The visual overlay composite image (Figure 21) indicates visual changes in land cover
change that took place in the upper Shire River catchment between 1989 and 2002. This
image highlights the changes in vegetation and qualitatively assesses the extent of the
change. This technique showed clearly areas of no change, and areas where changes have
taken place, in the upper Shire River catchment. The resulting image is grey in areas where
no change has occurred, such as to the west of the Shire River in the Liwonde National
Park. Where vegetation has decreased, the image is magenta/purple and where vegetation
has increased, it is green.
In Figure 21, the magenta regions represent vegetation that has been cleared for new
cultivation and settlements. Major vegetation reductions are apparent within Mangochi
Township near the outlet of the Shire River (Figure 22).
Mangochi is one of the fastest growing towns in Malawi due to the development of tourism
around the Shire River and Lake Malawi [National Statistical Office, 2000]. Its population
had increased from 496 578 in 1987 to 599 935 in 1998 as captured by population censuses
[National Statistical Office, 2000]. The high population growth has translated into rapidly
increasing demands from land in terms of food, shelter, energy (in particular, fuelwood)
and construction materials. Disappearing vegetation cover lessens the landscape's ability to
intercept, retain and transport precipitation. Instead of trapping precipitation, which then
percolates to groundwater systems, deforested areas become sources of surface water
run-off, which moves much faster than subsurface flows.
68
Figure 21: Image overlay for the Shire River catchment: 1989 — 2002
69
Mangochi TownshipMangochi Township
Figure 22: Expansion of Mangochi Township into previously vegetated areas
Another area of interest is the riverline below Lake Malombe, which shows a decreasing
trend in marshes, in contrast to an increase of marshland above the lake. Marshes are
critical because they recharge groundwater supplies and moderate streamflow. This is an
especially important function during periods of drought. The presence of marshes in a
catchment helps to reduce damage caused by floods by slowing and storing flood water. It
was verified during the field exercise (Section 2.3.2) that the grass-covered marshes have
either been drained and cultivated or overgrazed. These agricultural activities are part of an
initiative by Save the Children (UK), a non-governmental organisation, as a means of
increasing food security in the area.
In some parts of the south-eastern and north-western side of the catchment, some degree of
vegetation recovery is depicted. These areas were previously occupied by Mozambican
refugees during the civil war between 1976 and 1992. After the civil war ended, most
refugees returned to Mozambique. The harvesting of natural forest products by the
occupants of these refugee camps resulted in significant loss of forest area and increase in
fragmentation as captured by the 1989 image. By 2002, vegetation re-growth can be
observed in most of what were previously refugee camps. The green regions represent
70
higher reflectance in the near-infrared wavelength signifying vegetation re-growth through
natural processes and government intervention through reforestation programmes. At the
initial stage of the reforestation programmes, exotic trees such as bluegum and gmelina
were planted. However, water resources were being compromised due to these fast
growing species. The Malawi Government and NGOs working in this area, now
rehabilitate the camps with indigenous trees [Kafakoma, 1996]. As a result, the presence of
vegetation can change the quantity of water on the surface, in the soil or groundwater
recharge. This in turn affects run-off rates and the availability of water for either ecosystem
functions or human services.
3.4.2 Post classification and land cover change areas
Results from post-classification analysis are presented using a series of change maps for
visualisation, and statistical tables to provide quantitative measures of change. The change
maps use the colour coded from-to indices of change defined in Table 10.
Change Map
A map of overall land cover changes that have occurred within the Shire River catchment
is shown in Figure 23. The change information is overlaid onto the 2002 Landsat natural
colour image {3, 2, 1} for orientation purposes. Shades of bright green depict areas that did
not change between 1989 and 2002. Areas that have changed to cultivated or grazing lands
are shown in purple. Changes from other types of land cover to vegetated areas are
displayed in yellow and green. This map confirms the qualitative changes indicated in the
change overlay image of Figure 21.
71
no land cover change
� new built up areas
� change cultivated or grazing areas
�water related changesnot selected for display
??� water
� change to vegetated areas
vegetation regeneration from previously built up areas
� change to grassland areas
� change to vegetated areas
no land cover change
� new built up areas
� change cultivated or grazing areas
�water related changesnot selected for display
??� water
� change to vegetated areas
vegetation regeneration from previously built up areas
� change to grassland areas
� change to vegetated areas
no land cover change
� new built up areas
� change cultivated or grazing areas
�water related changesnot selected for display
??� water
� change to vegetated areas
vegetation regeneration from previously built up areas
� change to grassland areas
� change to vegetated areas
Figure 23: Post classification change map of the Shire River catchment between 1989 and 2002
Overall change statistics
Changes in each land cover category over the thirteen-year period 1989 to 2002 are
summarised in tables (Table 11 to Table 17) and discussed in this section. Table 11 shows
the net gains and net losses for the various land cover categories in hectares and
percentages. There are significant increases in the spatial extent of grasslands, cultivated
or grazing, and of woody closed classes, while savanna shrubs, woody open and built-up
areas decreased in extent. Unsurprisingly, the analysis shows insignificant changes in the
extent of fresh water and marshes.
72
Table 11: Land cover changes of the Shire River catchment during 1989 to 2002
Total area in land
class at 1989 ( ha )
Total area changed to per class
(ha)
Total area changed from per land class
(ha)
Net gain (loss) (ha)
Net change (%)
Built-up areas 39 813 456 2 005 -1 549 -3.9
Cultivated or grazing
95 428 4 622 3 158 1 464 1.5
Fresh water 38 353 0 108 -108 -0.3
Grasslands 15 127 3 616 469 3 147 20.8
Marshes 6 444 108 180 -72 -1.1
Savanna shrubs 152 791 2 584 4 370 -1 786 -1.2
Woody open 70 612 1 312 3 017 -1 705 -2.4
Woody closed 37 930 1 652 1 043 608 1.6
Total area 456 498 14 350 14 350 0
Changed areas for each land cover type are calculated by summing the respective numbers
of from-to change pixels occurring in the change detection matrix for each change
classification code (Table 10). Reference is made to the changed land area of the
catchment (excluding water bodies) amounting to 14 350 ha. The entire catchment
excluding open water bodies is 418 145 ha, thus the total changed area represents 3.43% of
the land area of the catchment. The following sections discuss changes within each land
cover class.
Cultivated or grazing areas
The total land area within the upper Shire catchment that was converted from various
classes to cultivated or grazing land amounted to 4 622 ha (Table 12). The increase in
cultivated or grazing land occurred mainly at the expense of savanna shrubs, woody open
areas and built-up areas: cultivated or grazing areas expanded by 1 671 ha (1.0%) from
previously savanna shrubs region; and 1 359 ha (1.9%) from previously woody open areas.
Surprisingly, there is a significant change from built-up areas (1 156 ha) to cultivated or
grazing land.
73
Table 12: Areas changed into cultivated or grazing areas between 1989 and 2002
From Total area in land class at 1989
( ha )
Area changed to cultivated ( ha )
% change
per class♣♣♣♣
Built-up areas 39 813 1 156 2.9
Grasslands 15 127 153 1.0
Marsh 6 444 24 0.4
Savanna shrubs 152 791 1 672 1.0
Woody open 70 612 1 359 1.9
Woody closed 37 930 259 0.7
Total 322 717 4 622 1.4
♣ Percentages calculated as fraction of the original area of the class.
The expansion of cultivated or grazing land indicates increased subsistence agriculture.
Savanna shrubs have experienced degradation and exploitation due to cultivation and
demand for wood resources. The rich calcimorphic alluvial soils render the savanna more
vulnerable to agricultural expansion. Most of the soils in the Shire rift valley are of alluvial
origin, rich in nutrients and ideal for agricultural production. The 2002 land cover map
shows that most of the remaining savanna shrubs are found in the Liwonde National Park
(protected area) and south east of Lake Malombe, which is sparsely populated (Figure 24).
A study of the two classification images shows that cultivated or grazing areas are also
expanding from the lowlands into the higher areas predominantly occupied by woodlands.
These processes are closely related to increases in the demand for food production.
Inevitably, cultivation and encroachment for exploitation of wood for construction and
fuelwood has spread onto such areas as can be noted from the 1989 compared to the 2002
land cover map. This has given rise to a transition from closed savanna woodlands to open
and sparse savanna, and it is mostly economic shrubs that are left in the fields (Figure 24).
74
1989
2002
1989
2002
Figure 24: Expansion of cultivated or grazing areas into predominantly savanna areas
An interesting result is the loss from built-up areas amounting to 1 156 ha which was
converted to cultivated or grazing land. Overall, this represents a net loss of -1 549 ha
(-3.9%) into other land cover categories (Table 11). Built-up areas (especially grass-
thatched) and cultivated or grazing areas occur in similar environments, often in adjacent
or mixed stands and have similar spectral signatures. The decrease could be an aspect of
misclassification due to similarity between roofing materials and the surrounding pixels
dominated by shallow soils and rock outcrops.
Additionally, the ravages of civil war in Mozambique had left large scars on the landscape,
particularly in the forest areas of the eastern mountains, which form the international
boundary between Malawi and Mozambique (Figure 25). Some bare patches where once
there were refugees camps have not yet been reforested, as can be seen in the 2002 image.
These bare patches from dismantled refugee camps are indistinguishable from fallow
cultivated fields and have probably been classified as change to cultivated or grazing areas.
75
Declining forest cover poses serious threats to water supplies within the area. A cause of
concern for this trend in land cover conversion is the significant role the forest reserves
play in catchment management. The forest reserves, Namizimu forest and Machinga hills
(Figure 20) serve some of the big rivers flowing into the Shire River. Overexploited
forestlands are threatened by accelerated run-off and soil erosion, and consequently
degrade the agricultural productivity of the land. Unless monitored, this trend may impact
negatively on the livelihoods of people and the survival of other natural resources in this
area, including water.
Machinga Forest
Namizimu Forest
Forest fragments Machinga Forest
Namizimu Forest
Forest fragments
Figure 25: Forest fragmentation around forest reserve areas
The predominant subsistence agricultural practices signify a gradual but chronic
degradation of the landscape along the valleys and adjoining slopes. These findings concur
with the observations by Kalipeni [1996], which indicates that pressures on land are
increasing, resulting in cultivation of marginal land, more fragile soils and steeper slopes.
Grasslands
From the post-classification land cover change results, 3 616 ha of land was lost from other
classes and converted to grasslands (Table 13). Taking into account losses of grasslands to
other types, the results in a net gain of 3 147 ha (+20.8%) in the extent of grasslands
76
(Table 11). High proportions of change originated from savanna shrubs (1 368 ha) and
cultivated or grazing lands (1 073 ha). Other changes include transition from woody open
land (553 ha), built-up areas (439 ha) and woody closed (158 ha). This is a considerable
change, in view of the short period over which the land cover changes have occurred.
Table 13: Areas changed into grassland areas between 1989 and 2002
From Total area in land class at 1989
( ha )
Area changed to grasslands
( ha )
% change
per class♣♣♣♣
Built-up areas 39 813 439 1.1
Cultivated or grazing 95 428 1 073 1.1
Marshes 6 444 25 0.4
Savanna shrubs 152 791 1 368 0.9
Woody open 70 612 553 0.8
Woody closed 37 930 158 0.4
Total 403 018 3 616 0.9
♣ Percentages calculated as fraction of the original area of the class.
The increase in grasslands from previously woody closed, woody open and savanna shrub
areas is a result of grass thriving best without canopy cover. Increasing demand for wood
resources places substantial strain on landscape integrity. However, the loss of built-up
areas to grassland could be a result of misclassification where the rural building materials
are from dry grass yielding similar spectral response to the natural grassland. In addition,
the surrounding pixels are made of bare soil, which increases the uniformity of spectral
appearance. As a result, spectral classification would assign all pixels including
grasslands, built-up areas and bare soil as a single class.
Changes in fire frequency and timing are known to transform vegetation structure and
composition. Frequent late dry season fires eventually transform woodland into open
grasslands with only isolated, fire-tolerant canopy trees, scattered understory and shrubs
[Desanker et al., 1997]. As the landscape becomes altered, run-off velocity and flow rates
in river systems are also affected.
Savanna shrubs
By analysing the changes that occurred to the landscape, a total of 2 586 ha (0.9%) was
converted to savanna shrub (Table 14). This increase does not compensate for the loss of
77
savanna shrubs to cultivated or grazing areas as discussed above. However, it could be a
result of misclassification as savanna shrubs and cultivated or grazing areas occurring in
similar environments.
Table 14: Areas changed into savanna shrubs areas between 1989 and 2002
From Total area in land class at 1989 ( ha )
Area changed to savanna shrubs ( ha )
% change
per class♣♣♣♣
Built-up areas 39 813 217 0.5
Cultivated or grazing 95 428 1 373 1.4
Grasslands 15 127 267 1.8
Marshes 6 444 32 0.5
Woody open 70 612 377 0.5
Woody closed 37 930 320 0.9
Total 265 354 2 586 0.97
♣ Percentages calculated as fraction of the original area of the class.
Built-up areas
Spatial extents in built-up areas show an overall subtle increase in settlements. Changes in
this land cover category are shown in Table 15 depicting higher land cover changes from
woody open areas (127 ha), cultivated or grazing areas (187 ha) and savanna shrubs
(87 ha). A detailed examination of built-up areas category across the catchment revealed
four important subclasses of change, which are discussed below.
Table 15: Areas changed into built-up areas between 1989 and 2002
From Total area in land class at 1989 ( ha )
Area changed to built-up areas ( ha )
% change
per class♣♣♣♣
Cultivated or grazing 95 428 187 0.47
Grasslands 15 127 28 0.18
Marshes 6 444 7 0.11
Savanna shrubs 152 791 87 0.06
Woody open 70 612 123 0.17
Woody closed 37 930 25 0.06
Total 378 322 457 0.12
♣ Percentages calculated as fraction of the original area of the class.
The spatial distribution of change in built-up areas is of particular interest in this thesis.
Most of the emerging towns, for example Liwonde Township, Mangochi Township and
Balaka Township portray an increase in built-up area signifying an aspect of rural-urban
78
migration (Figure 26). With the national annual urbanisation growth rate estimated at 6.7%
in 1998, many rural dwellers are migrating to urban and expanding district centres in
anticipation of wage employment [National Statistical Office, 2000]. There is a sparse
distribution of spontaneous rural settlements on the lower escarpments, across valleys and
flat lands. With the increasing population, slight expansions into the higher escarpment are
visible on the change image.
Mangochi Township
Built up areas
Mangochi Township
Built up areasBuilt up areas
Figure 26: Built-up areas expanding around Mangochi Township
A decrease in built-up areas is apparent in areas especially formally occupied by
Mozambican refugees. Most of the refugee camps were located to the north western part of
the study area around Mangochi district and Namwera Township. This area is close to the
border with Mozambique such that in 1989 during the civil war most refugee camps had
occupied this area. However, by 2002 the civil war was over and upon the repatriation of
Mozambican refugees, the Malawi Government embarked on land rehabilitation at the
79
former refugee camps, which has resulted in the slight increase of woody closed areas
(Figure 27).
Woody closed areasWoody closed areas
Figure 27: Increase in woody closed areas
Another note of significance is the effects of the electric fence constructed around Liwonde
National Park, a government reserved area, between 2002 and 2004. During the time of the
satellite observation in 1989, there were built-up areas (settlements) close to Liwonde
National Park. Due to prowling animals from the park, many settlements were abandoned,
leading to the regeneration of savanna shrubs as seen in the 2002 image (Figure 9). Field
observations and consultation with members of the surrounding community (in 2006)
indicated that since the containment of animals by the game fence, people have started
moving back to areas close to the park boundary [Palamuleni et al., 2008]. The area is
endowed with calcimorphic alluvial soils, which are very attractive for subsistence farming
(Figure 28).
The boundaries of the reserved area are assumed not to have changed, though subtle
increases in built-up areas do occur within the reserve. These negligible changes are due to
80
the expanding tourism industry within Liwonde National Park, which are visible in Figure
28. New built-up areas along the rivers are primarily composed of tourist resorts and camp
sites.
Liwonde National
Park
Tourist resorts and camp sites
Built up areas expanding outside Liwonde National
park
Tourist resorts and camp sites
Built up areas
Liwonde national park boundary
Liwonde National
Park
Liwonde National
Park
Tourist resorts and camp sites
Built up areas expanding outside Liwonde National
park
Tourist resorts and camp sites
Built up areas
Tourist resorts and camp sites
Built up areas expanding outside Liwonde National
park
Tourist resorts and camp sites
Tourist resorts and camp sites
Built up areas expanding outside Liwonde National
park
Tourist resorts and camp sites
Built up areas
Liwonde national park boundary
Figure 28: Increase in built-up areas and tourism expansion around Liwonde National Park
According to the image classification, the size of built-up areas has decreased between
1989 and 2002 with a net loss of -1 549 ha (-3.9%) into other land cover categories (Table
11). The decrease could be a result of misclassification because the population increased
from 1 011 843 in 1987 to 1 218 177 in 1998 as recorded by population censuses [National
Statistical Office, 2000] and population increases result in the expansion in built-up and
cultivated areas. The misclassification may be a result of local building materials, which
include grass (thatching) and sand (hand made clay bricks) such that built-up areas and
cultivated or grazing areas where there is little green grass and significant exposed soils
are spectrally identical. During summer, most of the areas are non-vegetated with exposed
soils due to ploughing by farmers in preparation for the next crop season. This condition
yields similar spectral values to those of cultivated areas; grass thatched houses and
81
exposed bare areas. Another plausible explaination could be increase in cluster and
nucleated settlement pattern as opposed to scattered settlements.
Percentage of the catchment surface that is impermeable due to settlement and road
surfaces influences the volume of water that runs and increases the amount of sediment
that can be moved [Arnold and Gibbons, 1996]. Generally, such trends in the landscape
could significantly alter the hydrology of the system.
Woody open
A total land area of 1 312 ha (0.38%) was converted from other land cover categories to
woody open land (Table 16). The transitions to woody open areas originated from land
areas previously occupied by marshes (455 ha), cultivated or grazing areas (377 ha),
woody closed areas (281 ha) and built-up areas (166 ha). Although this was a positive
trend, the extent of the transition from for example cultivated or grazing areas (377 ha) to
woody open areas was negligible in comparison to the 1 359 ha loss of woody open areas
to cultivated or grazing areas.
Table 16: Areas changed into woody open areas between 1989 and 2002
From Total area in land class at 1989
( ha )
Area changed to woody open
( ha )
% change
per class♣♣♣♣
Built-up areas 39 813 166 0.42
Cultivated or grazing 95 428 377 0.39
Grasslands 15 127 12 0.08
Marshes 6 444 455 0.30
Savanna shrubs 152 791 21 0.33
Woody closed 37 930 281 0.74
Total 347 533 1 312 0.38
♣ Percentages calculated as fraction of the original area of the class.
Classification errors could be attributed to this transition, which may not be a true
indication of the actual trends in this land cover category. Marshes, woody closed and
woody open areas have similar spectral characteristics as they both have high chlorophyll
during the dry season. The woody trees: Zizyphus jujube (masau), Mangifera indica
(mango) and Faidherbia albida (msangu) are in leaf during the dry season whereas
marshes are mostly evergreen due to presence of water throughout the year.
82
Woody closed
There was a slight increase of 1 652 ha (0.48%) in the area of the woody closed areas
(Table 17). A unique aspect observed in this land cover class is the reclamation of areas
previously occupied by Mozambican refugees and reforestation of school compounds. This
process is illustrated by regeneration from fallow fields especially those areas formerly
occupied by the Mozambican refugees in 1989 (Figure 27). Additionally, annual tree
planting exercises involve the planting of both indigenous and exotic trees such as
Eucalyptus globulus, Gmelina aborea and Acacia nigrescens in woodlots around schools.
As a result, there is a net gain of 608 ha (+1.6%) in woody closed areas across the
catchment area (Table 11). The woodlots had a similar spectral reflectance to the woody
closed areas as verified during field data collection and site context (Section 2.4.5).
Table 17: Areas changed into woody closed areas between 1989 and 2002
From Total area in land class at 1989
( ha )
Area changed to woody closed
( ha )
% change
per class♣♣♣♣
Built-up areas 39 813 28 0.1
Cultivated or grazing 95 428 149 0.2
Grasslands 15 127 10 0.1
Marshes 6 444 70 1.1
Savanna shrubs 152 791 789 0.5
Woody open 70 612 606 0.9
Total 347 533 1 652 0.5
♣ Percentages calculated as fraction of the original area of the class.
The proportion of rainfall that directly reaches the ground surface is affected by the
coverage of vegetation cover. Tree root systems hold the soil together and thus slow the
rate of run-off and reduce erosion. Trees also absorb water during the rainy season.
Because of the good ground cover, surface run-off is retarded and water can infiltrate the
topsoil, leading to high levels of infiltration and recharging of the aquifer. The positive
change in the extent of woody closed areas signifies ecological sustainability for water
resources management.
Although the increase in woody closed areas is a positive development in relation to
hydrological processes, this amount may be somewhat misleading because of the presence
of shadows in some parts of the images, which created difficulties in image classification.
83
It is possible to misclassify topographic shadows and woodlands as marshes. The
absorption effects of water in the marshes are spectrally similar to the lack of light
collected from the shadowed areas, while the strong reflection from the surrounding semi-
deciduous Miombo woodlands and the evergreen Eucalypts from woodlots exhibit spectral
similarity with the marshes. Most likely, misclassifications were expected. Nonetheless,
these areas have been identified in the DEM data for the region and they do not constitute a
large area. This potential misclassification is a result of the inherent limitations of the
technology and ancillary data at hand since during pre-processing no shadows were
removed. For instance, the 70 ha of marshland converted to woody closed area may not be
an accurate indication of trends in land cover change. The increase in the area classified as
woodlands during the period 1989-2002, however, is a positive development.
Marshes
The size of the marshy areas appears to be increasing, although it is possible non-marsh
areas were misclassified. However, water bodies do change in response to climatic
conditions and records of rainfall distribution in the Shire River catchment depict
continuous decreasing variations [Malawi Government, 1999]. As such, the water levels
along the shores of the river and the lake have been decreasing giving rise to waterlogged
areas, hence the increase of 108 ha.
Because marshes retain water, they support the vigorous growth of grass and provide good
dry season grazing areas when other forms of grazing are in short supply. Marsh margins
are also used for gardens (during the dry season) providing a more reliable crop output to
supplement rainfed harvests. However, marshy areas act as hydrological stores, holding
water and releasing it as base flow to the headwater streams during the dry season.
Continuous water extraction for dry season cultivation may be subject to uncertainty in the
variability of seasonal river recharge.
Fresh water
The changes in the extents of fresh water and marshes are insignificant, possibly due to
rainfall variations between the two time periods. The 1989 image was preceded by a La
Nina episode while the 2002 it was strong La Nina episode [SADC: Drought Monitoring
Centre]. In the 1989 image some very small dams were mapped which appear to have
dried up or have significantly reduced in surface area during the peak of the dry season in
2002.
84
3.5 Conclusion
Scientific approaches and holistic decisions regarding land management are required for
the sustainable development of the Shire River catchment. Change detection techniques
using temporal remote sensing data provide detailed information for detecting and
assessing land cover dynamics. Landsat classifications were used to produce accurate
landscape change maps and statistics. Two change detection techniques (multi-date image
overlay and post classification analysis) were applied to monitor land cover changes based
on data from the two dates (1989 and 2002). The result of this study documents a
significant combination of deforestation and fragmentation of forests in the study area.
General patterns and trends of land cover change in the Shire River catchment depict
transition towards degradation of woodlands and an increase in patches within the forest.
The change analysis showed an increase of 4 622 ha (7.9%) in cultivated or grazing lands
emanating from declining natural vegetated areas. The main driving factor is subsistence
agricultural expansion and demand for wood resources as evidenced by current agricultural
practices and population growth. As agriculture continues to play a dominant role in land
cover conversion and degradation from Brachystegia woodlands to more open and dry
vegetation formations will continue to evolve. To this end, it is apparent that the rapid
increase in cultivated areas within the catchment will not only decrease the amount of
forests and vegetated areas but also increase run-off potential thereby diminishing the
overall quality and quantity of water resources. Distinguishing and quantifying where
potentially risky changes occur is critical to the initiation of regular monitoring of
resources and the environment in general.
The patterns of land cover change observed in this study do not provide evidence of
impacts of industrialisation and commercialisation in the catchment over the last decade.
The strongest signals of landscape change detected thus far appear to correspond to an
apparent expansion of subsistence agricultural practices and changes in the amount and
extent of forest resources. Hence, the rapid growth in population, coupled with an absence
of industrialisation, will continue to put severe strains on the sustainability of
environmental resources.
Change data can be used to generate information tied to geographic coordinates and thus to
update of maps. The selected algorithms used in this study allow for the efficient thematic
updating of land cover change maps vis-à-vis land cover information for the Shire River
85
catchment and other areas in Malawi. The maps can also be used to estimate the rates of
land cover change. Recent information about land cover variables and the nature of the
transformation of land cover can provide a valuable guide for formulating appropriate
policies and for the effective implementation of programs for natural resource allocation,
land husbandry, conservation, management, sustainable use and combating deforestation.
Results from this study have demonstrated that satellite remote sensing approaches provide
a cost-effective alternative when more information is needed, but financial resources are
limited. Land cover change information generated from this study could serve the
requirements of the Department of Environmental Affairs in Malawi in their goal of
characterising the existing conditions of the Shire River catchment, both as a baseline for
later research and as the starting point for the development of future scenarios. Overall,
considering the present scale of temporal and spatial land cover change in the Shire River
catchment, a more continuous and comprehensive land cover change monitoring system is
required.
Land cover change information derived from remotely sensed data discussed in this
chapter will be used to examine the relationship between land cover variables and
hydrological regimes in the upper Shire River catchment. Chapter 4 provides an analysis of
the hydrological response from land cover activities.
86
Chapter 4
4 HYDROLOGICAL MODELLING BASED ON THE LAND COVER
ANALYSIS
Chapter 4 integrates the findings of the previous chapters into a conceptually based
spatially distributed ArcView Soil and Water Assessment Tool eXtendable version
(AVSWATX) hydrological model. This chapter presents the calibration, validation, and
application of the AVSWATX model for predicting the catchment hydrological
responses to changes in land cover in the Shire River catchment. Land cover change
simulations and scenarios demonstrate the model’s ability to integrate spatially distributed land cover change and precipitation events into output responses – in the
form of run-off from the catchment and total water yield at the outlet. Model
predictions are compared to observed hydrological records. The AVSWATX model provides a new approach for evaluating relative land cover changes across catchment
landscapes.
4.1 Introduction
4.1.1 Land cover and hydrological processes
Land surface heterogeneity is characteristic of many regions of the world [Lahmer et al.,
2001]. It plays an important role in partitioning incoming radiation at the land surface into
latent and sensible heat, and in partitioning precipitation into percolation, run-off and
evaporation [Calder, 2002]. At the larger scale typified by orographic features such as
mountains and valleys, topographic variation directly affects precipitation and surface
temperature. Precipitation and surface temperature are the main drivers of surface
hydrology through their effects on local and regional atmospheric circulation, the vertical
distribution of atmospheric moisture and temperature, and condensation [Maidment, 1993].
At a smaller scale, topographic variation can modify surface and subsurface run-off
through down-slope redistribution of soil water.
Precipitation and land hydrological processes maintain the water balance in a river basin.
Land surface performs a role in the hydrological cycle, as water availability is generally a
consequence of precipitation redistributed into evaporation, run-off and soil moisture
storage [Dolman and Verhagen, 2003]. The majority of precipitation must pass over the
land surface or drain through the soil and bedrock to translate into river flows. The spatial
heterogeneity associated with land cover, soil properties and localized precipitation
influences soil moisture and surface fluxes. Land cover change and the effects of land
management on the hydrological response of a catchment are most likely where the change
87
alters the surface characteristics of a basin. The degree and type of land cover influences
surface run-off and the rate of infiltration, and consequently the rate of ground water
recharge [Calder, 1992; Shaw, 1990]. Changes in these hydrological variables may have
implications for water resources.
Studies of the relationships between changes in land cover, environmental change, the
amount of run-off and percolation at the landscape scale can be used to compare
catchments, identify at-risk communities, and aid management attempts to limit undesired