Assessing and Predicting Changes in the Status of Gambari Forest Reserve, Nigeria Using Remote Sensing and GIS Techniques ()
1. Introduction
Tropical forests are the richest and the most diverse ecosystem on earth but are unsustainably over exploited despite legislation to control their exploitation [1] [2] . FAO [3] revealed that ninety percent of the original moist forest in West Africa is gone and the leftover is heavily fragmented and degraded, while in Nigeria, Federal Environmental Protection Agency [4] showed that over 43% of the total geographic area of 923,768 km2 has been lost in 12 years to human activities [5] . The threat to the rainforest compelled the establishment of forest reserves with the introduction of exotic tree species in the 1960s. Nigeria has 1160 constituted forest reserves covering a total land area of 10,752,702 hectares and this represents about 10% of the total land area [6] . Exotic tree species such as teak (Tectona grandis) and Gmelina were established in the forest reserves by the Federal and State governments in conjunction with the World Bank to provide raw materials and products in the form of poles, timber, veneer, wood particles, pit props, pulp, fuelwood and for research [7] [8] . However, as far back as thirty-one years ago, observation of tree population in a permanent sample plots in the forest reserves, for example the Gambari Forest Reserve located close to Ibadan, southwest Nigeria revealed that there was consistent decline in the number of trees per hectare [9] . In fact, FAO report on the rate of deforestation in Nigeria that between 1990 and 2010, Nigeria lost an average of 409,650 ha or 2.38% per year [10] . Presently, satellite imageries have shown serious deforestation and degradation of the Gambari Forest Reserve like many other reserves in the country due to the expansion of agriculture, fire, grazing, natural mortality, bush burning and more importantly indiscriminate exploitation of timber by illegal fellers [11] [12] . Adedeji [13] observed that the dereservation of forest reserves and poor forest management practices in the Gambari Forest Reserve among other were responsible for some of the indiscriminate exploitation of tree species at alarming rates. Forests are multifunctional ecosystems providing material and intangible services as carbon sinks, enhanced regional biodiversity, environmental services, forest-based economies [15] , and functions to living organisms, hence their losses are taken seriously [11] [14] . Loss of forest resources poses a significant threat to the achievement of Millennium Development Goals due to decline in ecosystem health [16] [17] . In view of the steady and rapid decline of the natural and plantation forest and the attendant environmental and health problems, it is essential that this land- change dynamic be appropriately monitored and analyzed.
Despite the importance of vegetation mapping as one of the basic components of ecological database of any region, Adejuwon [18] lamented that vegetation has continued to be inaccurately mapped in Nigeria. Furthermore, most decisions about how much land to use were and for what purposes (and the related consequences) are still poorly understood [19] . Therefore, there is an urgent need for proper geo-management of land and the concomitant availability of a detailed, accurate and up-to-date geo-information [20] . Understanding the dynamics of forest loss is critical for the planning and protection of biodiversity [21] [22] . Insights into patterns, rates, and trends of landscape changes are necessary to realize forest dynamics, enable preservation, and evaluate the effectiveness of management approaches [22] [23] . Land cover changes may result from human-induced land use changes or natural processes, i.e. climatic variability and natural disturbances [24] . Recent advances have demonstrated the feasibility of using coarse spatial resolution remote sensing imagery to assess tropical forest over large areas [25] - [27] . Remotely sensed data have in the past four decades steadily become invaluable information source for ecological characterization and survey [28] . This is due to the effectiveness of air-borne and space-borne remote sensing platforms and sensors that facilitate observation of biophysical attributes over extensive areas at multiple spatial, spectral and temporal scales [29] . Studying the dynamics of land cover requires maps that reflect the status of land cover at different times [30] . Land cover and land use changes are recognized as main agents of global change via their interactions with climate, ecosystem processes, biogeochemical cycles, biodiversity and human activities [31] . Change detection is the concept of identifying contrasts or discrepancies in the state of an object or phenomenon by observing it at different times [32] . Change detection studies are centred on identifying the biotic and abiotic components of the spectral and temporal changes that are occurring within ecosystems [33] . It involves the ability to quantify temporal alteration and transformation using multi-temporal data sets. Detecting, describing and understanding such changes are of considerable interest, not only to ecologists or environmental conservationists, but also to resource managers [34] . In fact, understanding of land cover change especially in relation to the forest is a pressing need for sustainable development because these changes can lead to land use conflicts due to the need for resources and space and the capacity of the land to absorb and support these needs [35] . Land-cover mapping determines the current composition and distribution of landscape attributes, and this is subsequently used as the basis for assessing future change [36] - [39] . Satellite remote sensing is an important source of data source for monitoring, detection, quantification, and mapping of forest cover patterns and changes, because of its repetitive data acquisition, digital format appropriate for computer processing, and accurate georeferencing approaches [40] - [44] . The advent of geographic information system (GIS) has made it possible to integrate multisource and multidate data for the generation of land use and land cover changes involving such information as the trend, rate, nature, location and magnitude of the changes [45] . However, there is a general lack of accurate and current LULC maps in Nigeria, especially concerning forest and other natural resources [46] . The aim of this study is to assess and predict in changes in status of Gambari forest reserve, Nigeria using remote sensing and GIS techniques. It also model future change in the area using Markov Chain Modeller in Idrisi [47] .
2. Methodology
2.1. Study Area
The Gambari Forest Reserve (Figure 1) is classified as IUCN category IV reserve. It was declared to be Ibadan
Figure 1. Gambari forest reserve, Nigeria.
forest reserve by a resolution of the then Ibadan council on 14th September 1899 [48] and take its name from a nearby village called Onigambari. The reserve is located on latitude 7˚25' and 7˚55'N and longitude 3˚53' and 3˚9'E. The forest reserve, which used to be part of the lowland rainforest, is now in the Guinea savanna and Derived savanna zones [11] . It covers a total land area of 13932.18 hectares consisting of five zones: Mamu, Onigambari, Busogboro, Odo-Ona/Onipe and Olonde [49] . The Mamu and Onigambari zones were predominantly made up of Gmelina plantations. The topography is generally undulating, lying at altitude between 120-150 m above sea level and the area is bounded by two Rivers (Rivers Ona and Awon). Mean annual rainfall ranges between 1200 mm to 1300 mm with two distinct wet seasons, May-July and September-November and a short dry period in August. The major dry season occurs between December and March. Temperatures are high throughout the year with a mean annual value of about 27˚C and an annual range of 30˚C [50] . Monthly relative humidity at 0900 hr GMT is more than 80 per cent, and at 1500 hr, the value fluctuates between 60 and 80 per cent except during January and February when it drops below 60 per cent [51] . The soil is generally brown sandy loam and described as having a depth of 45 cm with overlaying clay and iron stone gravel-soil that is about 1.21 m deep but varies from place to place and grouped as an Alfisols [52] . It is derived from the basement complex of metamorphic and granite rocks with intrusions of quartzite, schists, gneisses and iron concretions [53] . The clay mineral is kaolinite. Soil pH is about 6.5 - 7.0 which tends to rise with accumulation of litter. A large part of the original natural rainforest consisting of indigenous species such as Terminalia spp, Triplochiton scleroxylon, Irvingia garbonensis, and Treculia africana were cleared and replanted with exotic tree species of Gmelina arborea and Tectona grandis [54] . The natural and plantation forests presently lay side-by-side in the area.
2.2. Data Acquisition and Processing of Images
The images used in this study were downloaded from the official website of Global Land Cover Facility (GLCF) ―(http://www.glcf.umiacs.umd.edu). Detailed characteristics of the sources of data used for the study are shown in Table 1.
Landsat TM for 1984 and 2000 and Operational Land Imager/Thermal Infrared Sensor (OLI/TIRS) 2014 images were employed in the study (Table 1). The three scenes corresponded with path 191 and row 055 of the WRS-2 (Worldwide Reference System). The scenes were cloud-free and had good radiometric quality as reported by the data provider. The scenes were the best available satellite images for the study period. The imageries were processed in the Geographic Information System (GIS) environment. Image pre-processing including geometric, atmospheric and topographic corrections were carried out to ensure spatial and temporal comparability of the datasets [55] . The three sets of images were first geo-corrected and geo-referenced because effective image processing is critical to successful urban land use/land cover mapping and change detection [56] . Image pre-processing allow for conformity between multi-temporal imagery necessary for quantification and spatial comparisons [57] . The 2000 image was used as the base image for geometrical correction due to its better visual quality. Geometric correction of the 2000 scene was based on Ground Control Points (GCPs) identified on the topographical maps of the area. The 1984 and 2014 scenes were co-registered to the base image using additional GCPs into UTM projection with geometric errors of less than one pixel, so that all the images have the same coordinate system. Image enhancement method was carried out using the histogram equalization and the differences in the pixels are compensated for to produce uniformly distributed pixels along the output axis [58] . The preparation and cropping of the images to cover the specified area required was achieved using subset prepared using Erdas Imagine 10.0 software. Ground truth information were collected and combined with the images to assess the accuracy of image classification. All data were developed into Universal Transverse Mercator (UTM) coordinate system, zone 31, with World Geocoded System (UTM WGS 84) projection parameters. Figure 2
Table 1. Characteristics of the Landsat images used in the study.
Figure 2. Flowchart of methodology in the study.
shows the Flowchart of methodology of the study.
To map land use/land cover types as well as to analyse land use changes in the Gambari Forest Reserve overtime, false colour composite (FCC) images of Landsat-TM (1984, 2000) and Landsat-OLI/TIRS data (2014), were made and processed using ArcMap 10.0 software. On the images, water bodies appeared dark as it would absorb all the wavelengths, bare soil and bare rock surfaces appeared green, vegetation (forests and cultivated land) showed up in red and settlements and roads showed up in shades of cyan. A false colour composite from Landsat-TM uses an RGB combination of band 4 representing the NIR infrared, band 3 belongs to red and band 2 to green while in Landsat-OLI/TIRS uses an RGB combination of band 5 representing the NIR infrared, band 4 belongs to red and band 3 to green (Figure 3).
2.3. Image Classification
Image Classification is a process of categorising an image into a fewer number of individual classes, based on the reflectance values [59] . The goal of image classification is to extract landscape information from the satellite images [60] . Jensen [61] have discussed in detail the fundamental elements of image interpretation including grayscale tone, colour , height and depth, size, shape, texture, pattern shadow, site, association and arrangement. Image classification refers to the extraction of differentiated land cover and land use categories classes from raw remotely sensed digital satellite data [62] . The imageries were classified based on sample set created according to training samples. Training samples are representative of the desired land use classes [63] and were determined based on ground truthing, researcher’s personal experience and physiographical knowledge of the study area [61] .
Figure 3. False Colour Composite (FCC) of satellite imageries of Gambari Forest Reserve for 1984, 2000, and 2014.
Twenty (20) training samples were collected for each Landsat imagery and the Maximum Likelihood Classifier [64] was used in performing supervised classification in the GIS environment using ArcGIS 10.0 software. This procedure has proven to be a robust and consistent classifier for multi-date classifications [65] . The classified land use images of Gambari forest reserve for 1984, 2000 and 2014 were vectorised to show discrete land use types as units on the land use/land cover maps that depict predominant land cover type within relatively homogeneous areas in the study area. The land cover classes are important to make the classification easier [66] . Land cover classes adopted in the study were based on the classes proposed by the Coordination of Information on the Environment (CORINE). CORINE land cover is a map of the European environmental landscape based on interpretation of satellite images. It gives useful geo-referenced information for disaggregation [67] . The Level 1 classification system was followed and the data was classified into four land use/land cover classes spread over the study area. The classes are shown in Table 2.
A post-classification processing was applied using ancillary data to discriminate between forests and timber plantations as well as other land cover classes. Change trajectory of post classification comparison was used to map the patterns and extents of land use and land cover in the forest reserve as well as determine the magnitude of changes between the years of interest, 1984, 2000 and 2014, respectively. The classified images of different periods were then superimposed in GIS environment with the objective of filling a matrix showing transitions between the classified land use/cover categories. The evolving land use/land cover types were quantified using cross tabulation statistics to carry out land use/land cover change.
2.4. Accuracy Assessment of the Image
It is important to evaluate the accuracy of classification results and for this confusion/error matrices were used. This is an important method for evaluating per-pixel classification [62] . Accuracy assessments of the classified
Table 2. Description of major Land cover types defined in this study.
images were assessed using the three multi-temporal images. Overall accuracy of 94.62%) was achieved for the Landsat TM image of 1984, while overall accuracies of the Landsat TM of 2000 and Landsat OLI/TIRS of 2014 were 90.05% and 87.42% respectively. These indicated excellent overall classifications.
2.5. Change Detection Analysis
Change detection is defined as a process that is used to quantify the changes that are associated with land use and land cover in the landscape using geo-referenced multi-temporal remote sensing information [68] . To calculate the extent of each land cover class, we analysed classified maps using ArcGIS 10.0. Change detection analysis was carried out on Landsat images of different years (i.e. 1984, 2000 and 2014) to analyze the pattern and trend of change analysis in the study area using Land Change Modeller (LCM) for ecological sustainability and Markov chain model embedded in Idrisi [69] . Land Change Modeller is innovative land planning and decision support software, which allows rapid analysis of land cover change and simulate future land change scenarios [70] - [73] . The LCM relies on separate multi-temporal image classification and subsequent image comparison [63] - [64] . It is suited for analysis and prediction of LULC types and evaluation of implications of the changes on the entire ecosystem [69] [74] . LCM can be used to know detail spatial increase and loss, net change, net change drivers, tendencies of change and landscape prediction [75] . Furthermore, Change Analysis and Map Transition Option in LCM is a mapping tool to visualize the change that occurred from one land class all the other land classes [76] and predicting the course of change into the future. Using LCM requires mainly two time categorical maps and so the classified maps of 1984 (time-1) and 2014 (time-2) were used as inputs for the Change Analysis Tab of IDRISI. This enabled us to understand the gains and losses and the transition of areas among the land use/land cover classes; and to quantify the changes occurred from time-1 to time-2 [70] . Finally, pairs of LULC maps were overlaid to produce the LULC change map.
2.6. Markov Chain Modelling
A Markovian process is one in which the future state of a system can be modelled purely based on the immediately preceding state. Markovian chain analysis will describe land use change from one period to another and use this as the basis to project future changes. Modelling help to provide a better understanding of the factors that drive forest changes and generate future forest cover scenarios to support the design of policy responses [77] . Markov Chain helps to determines exactly how much land would be expected to transit from the later date to the predicted date based on a projection of the transition potentials into the future and creates a transition probabilities file [70] - [71] [78] . This is achieved by developing a transition probability matrix of land use change from time one to time two, which shows the nature of change while still serving as the basis for projecting to a later time period [70] [79] . The Markov chain projection model was implemented on two classified LULC maps of different times [80] . The first image (2000) as the earlier image and the second image (2014) as the latter image was inputted. The model according to Iacono et al. [80] is based on the assumption that a future state (t2) can be determined by its current state (t1). In LULC modelling, the process determines the t1 to t2 LULC distribution using a transition matrix. It is expressed as:
(1)
where vt1 is the LULC proportion vector input, vt2 is the LULC proportional vector output and M is the m × m transition matrix for the time difference Δt = t2 − t1. The transition probabilities file is a matrix that records the probability that each land cover category will change to every other category.
2.7. Statistical Analysis
Land cover maps were employed to quantity land cover changes between 2000 and 2014 according to Melendez- Pastor et al. [81] . Land cover change statistics were computed in two different ways: 1) as absolute percentage increments of the whole study area calculated by subtracting percentage areas among 2014 and 2000; and 2) as relative percentage increments from 2000. The negative symbol in the statistics indicated a loss of surface. Both approaches provided complementary information.
3. Results
3.1. Land Change Process
Land cover change detected from comparison of CORINE Land Cover data for 1984, 2000 and 2014 showed that over the whole study period, forest cover in the Gambari reserve has decreased dramatically. The most extensive land cover class in 1984 was the thick forest, which occupied more than a half of the area, followed by the plantation of exotic trees such as teak and Gmelina. Within the sixteen tears from 1984 and 2000 land cover change dynamic was dominated by a considerable increase plantation and agricultural land use (Table 3). However, there was reduction in the extent of the plantation due to dereservation policies and illegal felling of trees in the reserve. The thick forest have drastically reduced from 7514.4 ha in 1984 to 2274.39 ha by 2014, while there was increase in bare ground from mere 139.42 ha in 1984 to 1673.1 ha as of 2014. This an indication of rapid deforestation and degradation of land in the forest reserve overtime.
Agriculture land use also declined from 3177.03 ha in 1984 to 1431.18 ha by 2014, as many farm were abandon after the decline in soil fertility due to long time of tillage. Increase in the bareground is also due to increasing quarrying activities in parts of the reserve, which despoil the lands.
3.2. Land Use/Land Cover Matrix
A cross-tabulation procedure between the classifications was processed with IDRISI Selva [82] ; area change, gains, losses and persistence were calculated as proposed by Pontius et al. [83] . Systematic transitions among classes were calculated and examined through the off-diagonal entries of the cross-tabulation matrix. Table 4 shows the summary of the major land use /land cover conversion in the Gambari Forest Reserve that took place between1984 and 2000. The diagonal of the table shows the land use/land cover that remained unchanged during the period, which is a total area of 5203.5 ha. This represented about 37.34% of the study area.
Thick forest had the highest conversion of 2906.74 ha to cultivated land, representing 20.79% of the area. About 22.33 ha of the area formerly occupied by plantation representing 0.16% of the total area were laid bare due to removal of trees without replanting. The thick forest area observed further disturbances as 45.65 ha (0.33%) and 372.57 (2.66%) were converted to light forest and agricultural lands respectively during the two periods.
The area covered by plantations of exotic trees of teak and gmelina declined from 2566.53 ha in 1984 to
Table 3. Size and proportion of land cover classes from 1984-2014.
Table 4. Land use/Land cover matrix of Gambari Forest Reserve between 1984 and 2000 (Hectares).
Grey hade means permanence, underline indicate significant transition.
2291.63 ha in 2000. Further land use conversions were witnessed in the period between 2000 and 2014 as the thick forest correspondingly decreased from 3762.9 ha in 1984 to 2274.39 has in 2014 (Table 5). Within the 30-year period, the forest was lost due to illegal logging and other land uses including land clearance for agriculture. Furthermore, 514.17 hectares (5.05%) of the land used for cultivation as of 2000 became bare ground fourteen years later in 2014 due to inappropriate land management such as bushing burning, shifting cultivation, and land abandonment. The area degraded by various anthropogenic activities including farming and quarrying within the study area increased from 111.73 ha in 984 to 751.86 ha in 2000 and further to 1673.1 ha by 2014. The structure of the vegetation in the Gambari Forest Reserve is rapidly changing because of human activities in all the zones making up the reserve. These activities including illegal logging, unsustainable agricultural practices (e.g. bush burning), and soil erosion are causing the loss of forest lands. Depletion of unreserved forests has resulted in the concentration of logging activities within the forest reserves by both local and foreign timber merchants in conjunction with corrupt forestry officers. In fact, Kio [84] indicated that for the former state of western Nigeria the areas outside forest reserves accounted for 52% of the total timber volume output in 1960 and that this percentage was reduced to 16% five years later. Total area occupied by natural forest type in the Gambari forest reserve was about 38 percent [85] . Land cover maps were obtained from satellite imagery for 1984, 2000 and 2014 (Figure 2).
4. Discussions and Interpretation
Land use/land cover changes are critical issue that degrade biodiversity and create impact on human life [86] . It should be noted that great environmental consequences could result with the continuous decline in the forest resources with gradual reduction of biodiversity and landscape complexity and a higher vulnerability to natural hazards [87] . Major driver of land use change in the study area include increasing demand for timber products, proximity of the forest reserve to big cities such as Lagos and Ibadan where there is demand for timber for both local use and export. Accessibility to forest is another driver because the forest reserve is well linked with good and accessible roads.
4.1. Gain and Losses between 1984 and 2000, 2000 and 2014
Figure 5 and Figure 6 show the gain and losses in land use/ land cover that occurred in the forest reserve from 1984 to 2000 and 2000 to 2014. The green bars (right side) indicate the gain per class in hectares, while the left side define the loss of each class. In terms of absolute gain, the light forest gained much from the felling of trees in the thick forest which latter replanted or grow on its own.
The light forest had the highest amount of gains between 1984 and 2000 (5542 hectares) while about 4248 hectares of the thick forest were lost within the same period (Figure 4). Between 2000 and 2014 more of the thick forest were lost with increase in the amount of bare ground (1520 hectares).
Figure 7 shows the net changes in land use/land cover from 2000 to 2014, while Figure 8 shows the contributions to net changes of each land use categories between 1984 to 2000 and 2000 to 2014 .
The study also calculated the contributions of other categories to their net change (Figure 8). This showed which classes in the 1984 map were identified as thick forest, light forest, plantation, agricultural land use or
Figure 4. Classified satellite imageries of Gambari Forest Reserve for 1984, 2000, and 2014.
Table 5. Land use / Land cover matrix of Gambari Forest Reserve between 2000 and 2014 (Hectares).
Grey shade means permanence, underline indicate significant transition.
bare ground in the 2000 map and what proportion of the total change for a class and proportion of the study area they explain. For instance light forest explain the majority of the total decrease in thick forest (2854 hectares), while losses in plantation, agricultural land use and thick forest explain the increase in the bare ground. Figure 9 show the major contributors to net changes in land use/land cover from 2000 to 2014. Forest recovery through succession was detected by the increase in the light forest after the thick forest was removed.
4.2. Trends and Magnitude of Land Use Categories
Table 6 provides information on the overall land cover changes between 2000 and 2014 based on the land cover
Table 6. Land use/land covers change: Trend and Magnitude 2000 and 2014.
Figure 5. Gain and losses (Ha) in land use/land cover from 1984 to 2014.
Figure 6. Gain and losses (Ha) in land use/land cover from 2000 to 2014.
Figure 7. Net changes (Ha) in land use/land cover at Gambari Forest Reserve from 1984 to 2014.
maps. Both the absolute percentage increments between the two years and relative percentage increments from 2000 were calculated. Negative symbols in the statistics indicated a loss of surface [81] .
The plantation, thick forest and agricultural land use suffered losses in both absolute and relative terms between 2000 and 2014. In absolute term, 10.69% of the thick forest and 5.63% of the plantation were lost during the fourteen-year period. These losses translated into increase in the amount of bare ground and light forest. Figure 10 shows the transition map between 2000 and 2014.
Figure 8. Contributions (Ha) to net changes (1984-2000).
Figure 9. Contributions (Ha) to net changes (2000-2014).
Figure 10. Land use transition between 2000 and 2014.
4.3. Land Cover Change Prediction with Land Change Modeller
The analysis of past changes or LULC distributions with regard to spatial explanatory variables enables to assess the degree to which locations might likely change in the future. Markov analysis was performed for the multi-temporal land cover images of 1984-2014 including predicting for 2044 as shown in Figure 11. Results of the Markov models are focused predominantly on providing the knowledge of how much, where, what type of land use changes has occurred from 2014 to 2044. The prediction of the future scenarios in land use change in the study area suggested a continued degradation of the plantation and thick forest resulting in more degraded lands. It should be noted that vegetation loss and degradation could reduce precipitation infiltration and runoff regulation in the forest reserve and promotes soil erosion.
5. Conclusion
This study made use of Landsat TM imageries of 1984 and 2000, and OLI/TIRS of 2014 to identify, classify, assess and interpret changes in the Gambari Forest Reserve for the year 1984, 2000 and 2014 respectively. The land cover categories and their changes for these years were generated and analysed in the GIS environment using ArcGIS 10.0 and the land use/land cover scenarios for 30 years from 2014 to 2044 was performed using Markov analysis in Idrisi. The result shows that the thick forest and plantation of exotic tree species i.e. teak and gmelina were fast declining due to several anthropogenic activities such as illegal felling of trees for timber, fuelwood and farming. The degradation can be attributed to the pressure of urbanization, poor management of the forest reserve and increasing human activities in the area. Future prediction showed an almost complete absence of thick forest in the area with increase decline in the plantation, which can be very dangerous for biodiversity, and sustainability of the forest resources. The study have shown the versatility of remote sensing and
GIS in mapping and monitoring the dynamics of land use/land cover of the study area and provide a useful background for identifying the hot spot and prepare appropriate conservation and management programmes for ecological sustainability.
Acknowledgements
The authors wish to thank several contributors for the supply of data, reading the manuscript and making useful contributions.
NOTES
*Corresponding author.