Monitoring Land Cover Change Using Remote Sensing (RS) and Geographical Information System (GIS): A Case of Golden Pride and Geita Gold Mines, Tanzania

This study monitored land cover change in the mining sites of Golden Pride Gold Mine (GPGM) and Geita Gold Mine (GGM), Tanzania. The satellite data for land cover classification for the years 1997, 2010 and 2017 were obtained from the United States Geologic Survey Departments (USGS) online database and were analyzed using Arc GIS 10 software. Supervised classification composed of seven classes namely forest, bushland, agriculture, water, bare soil, urban area and grassland, was designed for this study, in order to classify Landsat images into thematic maps. In addition, future land cover changes for the year 2027 were simulated using a Cellular Automata (CA)Markov model after validating the model using the Land Cover for the year 2017. The results from the LULC analysis showed that forest was the most dominant land cover type in 1997 at GPGM and GGM covering 510 ha (52.1%) and 9833 ha (49.7%) respectively. In 2017, the forest area decreased and the bushland replaced forest to be the most dominant land cover type covering 219 ha (22.4%) for GPGM and 8878 ha (44.9%) for GGM. Based on the CA-Markov model, a predicted land cover map for 2027 was dominated by forest covering 340 ha (34.7%) and 8639 ha (43.7%) for GPGM and GGM respectively. An overall accuracy and kappa coefficient for GPGM were 74.7% and 70.2% respectively and for GGM were 71.4% and 66.1% respectively. Thus, land cover changes resulting from mining activities involve reduction of forest land hence endangers biodiversity. GIS and remote sensing technologies are potential to detect the trend of changes and predict future land cover. The findings are crucial as it provides basis for land use planning and intensifies monitoring programs in the mining areas of How to cite this paper: Kahangwa, C., Nahonyo, C. and Sangu, G. (2020) Monitoring Land Cover Change Using Remote Sensing (RS) and Geographical Information System (GIS): A Case of Golden Pride and Geita Gold Mines, Tanzania. Journal of Geographic Information System, 12, 387-410. https://doi.org/10.4236/jgis.2020.125024 Received: August 3, 2020 Accepted: September 7, 2020 Published: September 10, 2020 Copyright © 2020 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Mining is an extraction of valuable materials from the earth's surface usually from an ore which consequently results in loss of vegetation, loss of soil and loss of rocks [1]. Mining regions are geologically defined by economically feasible mineral resources and they often undergo abrupt and extensive changes in the land cover [2]. Large-scale mining is one of the main growing industries in Tanzania. Large-scale mining especially open-pit involves the clearing of large area of land and forest resulting in serious land cover changes [3]. During the past two decades, Golden Pride Gold Mine (GPGM) and Geita Gold Mine (GGM) among the gold mines within Lake Victoria Gold field have undergone significant land cover changes due to mining operations within the licensed mining areas. Land cover is the biotic or abiotic features that cover the earth surface such as forest, bare soil, grassland, bushland and water [3].
From environmental point of view, the dynamic process of land cover change has been described as the most remarkable anthropogenic disruption of the environment [4] [5]. It is a dynamic variable because it reflects the interaction between socio-economic activities and regional environmental changes, and for this reason, it is necessary to be updated frequently [6]. Land Cover is necessary in identifying and mapping natural resources and human activities which have a great effect on landscape and soil degradation [7] [8]. Information on land cover plays a key role in natural resources management and spatial planning [9] [10].
This study used Remote Sensing (RS) and Geographical Information System (GIS) to monitor the trend of land cover change in the mining areas and to predict the future land cover [11]. Remote sensing techniques have been used over the years and have proven to be of great value for monitoring changes at regular intervals [8]. Using satellite remote sensing data is a practical option to identify and map the land cover categories [12]. Land cover mapping using satellite remote sensing images has become widely popular in the last decades [13]. The reasons for selecting satellite image periods for this study considered the phases of mining projects from pre-mining, during active operations and post-mining phases. In addition, GIS tools are used to create the geo-database and integrate data extracted from satellite images with classes from the currently available land cover models [6]. However, developing countries still lag behind significantly compared with the developed countries in the application of remote sensing and GIS technologies. Land cover change models are very helpful for evaluating anthropogenic and natural impacts, predicting future patterns of land use changes, Journal of Geographic Information System aiding stakeholders of natural resources to achieve sustainable environments and planning and management of land transformation [11]. A GIS can be integrated to predict the future scenario of land cover using Cellular Automata (CA) Markov model. CA Markov model is widely used to predict and analyze likely future trajectory of LULC within a specified period and to characterize the dynamics of land use/cover, forest cover, coastal management, wetland landscape and so forth [14] [15] [16] [17]. The heart of a CA Markov model is the transition matrix between different temporal images to predict the future events [18] [19]. In general, various traditional modellings in geography create a spatial layout of present and past land surfaces using dynamic optimization techniques [18].
Commonly used models for estimating land cover changes are statistical models, analytical equation-based models, evolutionary models, cellular models, Markov models, hybrid models, expert system models and multi-agent models [20].
Currently, the most commonly used models in land cover change monitoring and prediction are the mixed model [20]. The CA-Markov model is the mixture of the Cellular Automata and Markov models. CA Markov model has many advantages for geographic modelling; it is capable of supporting very large parameter spaces for simulation, the model makes an implicit use of space and spatial complexity by simple local raster transformation, it is capable of combining space-time relationships, as well as supporting separate notions of space and time. Moreover, the model is capable of representing pattern, process, form and function. CA Markov model is vital in helping land use planners to make land use planning decision.
Tanzania, like many other developing countries, the large-scale gold mining activities typically decrease land covered by vegetation through the process of forest clearing for formation of mining components such as pits, Waste Rock Dumps and Tailing Storage Facilities. However, the countries are still using traditional field investigation method instead of remote sensing and GIS technologies to monitor the trend of land cover changes in the mining areas. This traditional technology is low efficient, high labor demanding and impossible to conduct continuously, especially in large scale areas. Therefore, the use of remote sensing (RS) in combination with Geographical Information System (GIS) has proved to be effective in the spatial and temporal monitoring of the land cover changes in the mining areas.  [23]. In general, studies revealed that mining activities cause extensive land use changes, which pose significant management challenges for mining companies and regulatory authorities. The effects identified in the studies include land degradation and conversion of land from forest covers to non-forest covers such as built-up areas, bare lands, Waste Rock Dumps and abandoned quarries filled with water.
In Tanzania, [24] monitored geomorphic and hydrologic change at GGM us-Journal of Geographic Information System ing Digital Elevation Model (DEM) data from the Space shuttle Radar Topography Mission (SRTM) and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). The authors assessed the effects of mining activities on farmers, village water supplies, and community forests using a hydrologic flow model. The results indicated that more than 81 million m 3 of waste has been deposited in the tailings storage facility and waste rock dumps. [1] also analyzed the impact of mining activities on land cover using Arc View GIS and ERDAS Imagine. The results indicated that there was significant conversion of forest cover into non-forest cover from 1991-2011 in Nanyala Ward. The author concluded that as the mining activities increase, the area affected by mining increases posing a serious threat to the environments. From the above background, it is evidently that, as far as land cover change due to gold mining operation is concerned, there is limited knowledge evaluating and predicting the future land covers in the gold mining environment of Tanzania. The aim of this study therefore, was to monitor land cover changes from 1997 to 2017 and then predict the future land cover for the year 2027 at GPGM and GGM mining licensed areas. The prediction of future land cover is important in the sense that, it will enhance understanding of the present land cover and will provide some foundation for effective planning of the future land use. This will help in monitoring the mining landscape for sustainable development to meet the needs of the present and future generation. We recommend effective management of degraded mining areas by intensifying progressive rehabilitation programs through planting trees as this will help to achieve positive post mine closure outcomes.

Overview of the Study Area
GPGM and GGM are located in the Sukuma land Greenstone belt within Lake Victoria Goldfield in Northern Tanzania. Lake Victoria Goldfield is among the three principle goldfields which are the main producers of gold in Tanzania.
Other goldfields are Mpanda and Lupa located in the Western and Southern Tanzania respectively [25]. Geographically, the GPGM is located at latitude: 040 23'31"S and longitude 0320 53'55"E in Nzega District, Tabora Region, Tanzania.
It is approximately 18 km north of Nzega Township, and 200 km south of Mwanza regional center. Average annual rainfall is between 700 -800 mm per year and average temperatures range from 22˚C to 27˚C [26]. Whereas, GGM is located at latitude: 020 52'03"S and longitude: 0320 11'11"E it is approximately 4 km west of Geita town and 90 km south-west of Mwanza City in north-western Tanzania ( Figure 1). There is very little annual variation in daytime temperatures throughout the year with an annual average daily temperature of 25˚C as recorded at GGM. Evaporation rates throughout the year are related to the monthly rainfall and are much greater during the dry season, yet are overall relatively low, ranging between 70 and 140 mm per month.

Research Materials
This paper covers Land Cover classification and Accuracy Assessment. In order to determine the effects of large-scale mining on the land cover of the study area, land cover satellite data sets for the years 1997 and 2010 were delivered from Landsat Thematic Mapper and for the year 2017 were delivered from Landsat 8 Operational Land Imager (OLI). The satellite data sets were obtained from the United States Geologic Survey Department (USGS) online database through http//earthexplorer.gov. Various digital image processing techniques were employed to prepare land cover maps of the study areas.

Remote Sensing Data Collection
To prepare the base maps, the Landsat satellite images (1997, 2010 and 2017) were collected from the official website of US Geological Survey (USGS). Landsat Path 171 row 62 covered GGM and Landsat Path 170 row 63 covered GPGM.
The pixel size of the images was 30 × 30 m. Accordingly, the land cover was classified into seven classes namely forest, bushland, agriculture, water, bare soil, urban area and grassland, for the three sampling periods i.e. 1997, 2010 and 2017. Journal of Geographic Information System

Data Acquisition and Processing
Time analyses of land cover change require a proper selection and preparation to ensure the compatibility of the Landsat images. In this study, the images were selected from the same season (July-September) and with minimal cloud cover (<10%). The images were from different sensors with similar spectral resolution i.e. 30 metres.

Land Cover Change Detection
The land cover change detection was done using ArcGIS10.3 software. The two classified Land cover layers i.e. Land cover 1997 and 2010 and Land cover 2010 and 2017 were used. Spatial analysis was done using the zonal tabulate area function and was used to generate land cover change matrix. The function calculates cross-tabulated areas between two datasets and provides output of a table that displays a record for each unique value of the zone dataset and a field for each unique value of the class dataset. Also, the land cover change analyses were performed and cross tabulation was done using Spatial analysis function in Arc-GIS 10 software [5].

Land Cover Change Modelling
Coupled Cellular Automata (CA)-Markov model was employed to conduct Land cover change modelling. The CA-Markov model combination represents advancement in spatio-temporal dynamic modelling and forecasting, and enable achieving a better simulation of land cover changes both in quantity and space based on studies of land use changes in the past [14] [20]. The algorithms in the IDRIS Andes package integrate the functions of the CA filter and Markov process, using conversion tables and conditional probabilities from the conversion map applied to simulate and forecast the states of Land cover change [15]. Therefore, to simulate future land cover changes for this study sites using a CA-Markov model, the following specific processes were followed: 1) After converting the vector data to raster, the classified land cover maps for the years 1997, 2010 and 2017 were used to obtain the transition matrices for the land cover categories between 1997 and 2010 as well as 2010 and 2017 based on the first-order Markov model. 2) The transition suitability maps, which were used to predict the land cover in 2017 and to simulate the distribution in 2027 were generated based on the main transitions that occurred among the land cover categories from 2010 to 2017. In addition, to determine CA filters, the regular 5 × 5 contiguity filter was used as the neighbourhood definition. 3) Based on the CA-Markov model approach, the land cover for the year 2017 was modelled using the transition probabilities from 1997 to 2010 with the land cover base map from the year 2010. Kappa statistics were used to assess the accuracy of the forecasted 2017 land cover map to evaluate its agreement with the actual 2017 land cover map. 4) Finally, following the same process, the land cover for the year 2027 was projected with the CA-Markov model in IDRISI using the transition probabilities from 2010 to 2017 and the land cover base map from the year 2017. Journal of Geographic Information System

Cross Tabulation
Idrisi software was used to perform a cross-tabulation analysis. It compared images containing categorical variables of two types. Cross Tabulation provides the information on the frequencies with which each land cover classes remained either unchanged or has changed to one of the other classes, using two thematic maps of different dates [27]. Three cross-tabulation tables were generated from thematic maps 1997-2010, 2010-2017 and 2017-2027.

Accuracy Assessment
The process of determining accuracy for a classified image is one of the most important steps undertaken in post-classification [1]. To acquire a better accuracy, ground truthing, or physical appearance in the study site is necessary [28].
However, in these cases, training points located in remote areas could not be visited in the field due to time and cost constraints. Kappa statistics was used to assess the accuracy of the forecasted land cover map to evaluate its agreement with the actual land cover map and the analysis was done using validate tool in Idris Selva software. The classified land cover for year 2017 was compared using projected land cover for year 2017.Kappa coefficient estimates the agreement between map and reality and it ranges from 0 to 1, where 0 represent total disagreement and 1 total agreement [6]. [29] also provides the interpretation of Kappa statistics as presented in Table 1.

Land Cover Classification
The Landsat images for the study area were classified in order to identify the changes in the cover between the four periods, i.e. 1997 Table 2 and Table 3 for GPGM and GGM respectively.
Transition matrices are presented on Tables 4-6 for GPGM and Tables 7-9 for GGM. The land cover categories delineated in the study areas reveal the changes in land cover (in hectares and percentages) that have taken place during the period between 1997-2017 as being due to the expansion of mining activities. Classified Accuracy Assessment report and Kappa statistics for GPGM are provided in Table 10 and Table 11 respectively and for GGM are presented on Table 12 and Table 13 respectively. Journal of Geographic Information System

Land Cover Map at GPGM
The examination of the land cover map in 1997, showed the forest to be the dominant land cover in the study area. During that period there were no mining activities in the area. The project commenced operations in 1998. Thus, by 2010 part of the forest has been replaced by other land covers and bushland is seen as the dominant land cover, due to the fact that, re-vegetation of the degraded mining sites was done concurrently with the mining operations. The satellite image of 2017 has shown that the area which was initially covered by forest has not returned into its original state and bushland is seen to dominate the area ( Figure 2). During this period the project had already ceased and the areas which were seen covered by bushland are the rehabilitated areas.

Change Detection Analysis at GPGM
Gold Mining in GPGM resulted in widespread land cover change between 1997-2017. Data in Table 2    probability is forest and the lowest unchanged transition probability is water. Tables 7-9   ha and water bodies 86.6 ha ( Table 9). The results of the transition matrices for the years 1997-2010, 2010-2017 and 2017-2027 indicate clearly that, the highest unchanged transition probability is forest and the lowest unchanged transition probability is water.

Analysis of Future Trend of Land Cover Classification
In this study, the 1997, 2010 and 2017 land cover maps were used to classify the future land cover of 2027. Future land covers can be observed in Table 3, Figure   4 for GPGM and Table 4, Figure 5 for GGM. It is seen in the projected land cover map of 2027 of GPGM that, the large part of the area is covered by forest followed by bushland then agricultural land with scattered settlements ( Figure   2). These projected results may be useful references to the regulatory agencies and the mining companies in dealing with management challenges and planning for future land use. The future land cover map within the specified period was analyzed using CA-Marcov model.
In the projected land cover map of 2027 of GGM, the large part of the area is seen to be covered by forest followed by bushland then grassland (Figure 3).
These projected results may be useful references to the regulatory agencies and the mining companies in dealing with management challenges and planning for future land use. Marcov chain model was also used to analyze the future land cover map within the specified period.

Classification of Accuracy Assessment for GPGM
The final stage of classification is to determine accuracy assessment. Each of the land cover map was compared to the reference data to assess the accuracy of the classified images. Accuracy Assessment classified for GPGM presented in the Classified accuracy assessment report (Table 10) indicated that the overall accuracy is 74.7% and the accuracy assessment of the forecasted land cover map to evaluate its agreement with the actual land cover using Kappa statistics was 70.2% (Table 11). These assessments indicate that there is substantial agreement between the classified land cover and actual land cover.

Classification Accuracy Assessment for GGM
Accuracy Assessment classified for GGM presented in the Classified accuracy assessment report (Table 12) indicate that the overall accuracy is 71.4% which is Journal of Geographic Information System acceptable and the accuracy assessment of the forecasted land cover map to evaluate its agreement with the actual land cover using Kappa statistics was 66.1% (Table 13) indicating almost perfect agreement.

Land Cover Change Detection
The land cover changes observed in GPGM and GGM have superficial impacts on plant species in natural habitats in the past 20 years. The multi-temporal satellite data generated from the mining areas provide evidence on the detected changes in land cover types based on the development stage of the mining activities. The impacts of mining operations in the study area within the period under review indirectly influenced adjacent land uses [2]. The results from the Land Cover analysis showed that, for the year 1997 at GPGM before the project started operations in 1998, forest was found to be the dominant land cover type, covering over 50% of the total study area, followed by bushland, then followed by agriculture with scattered settlements, grassland, built up area, bare soil and the water were not detected in this period. processing plant [24], infrastructure (access and energy) and construction of camps and town [30]. This is supported in the study conducted by Perpetual [16] who found that, between 1990-2010 high density forest and sparse forest in the Prestea mining area in Ghana were lost to pave way for mining operations.
For the case of decreased agricultural land, this study is in line with the study conducted by [31], on the effects of loss of agricultural land due to large-scale gold mining on agriculture in Ghana. The authors confirmed that, a total of 4935.3 ha of agricultural land, representing 25.5% of the three concessions, were lost due to large-scale gold mining operations in the study period. The results are also supported in the studies conducted by [23] who found out that, there was substantial loss of forest resulting from mining activities to pave way for mining operations such as pits development. Moreover, [32] pointed out that, mining activities have disturbed biodiversity at the mining sites of Nigeria. The major focus in this study is the natural forest cover which is a major habitat for fauna.
The decrease in the area of forest and agriculture in the active mining phase results in the increase of the areas for bushland, grassland, built up area, water and bare soil. For example, the built-up area increased because of the construction of staging area that would house project personnel, equipment and the processing plant. This is in line with the study conducted by [33], who evaluated land use/cover change with temporal satellite data and information systems and Journal of Geographic Information System

Transition Matrices
GPGM and GGM are among the large-scale gold mines in Tanzania, where the forest covers have been significantly changed to non-forest covers leading to superficial loss of vegetation. Before the commencement of the projects local people from the surrounding environments used the areas for anthropogenic activities such as agriculture, cattle grazing and timber harvesting. These activities were among the drivers of land cover changes in the areas [26]. After the Journal of Geographic Information System projects took off these activities ceased and mining and its associated activities have been among key drivers of land cover changes experienced in the licensed areas.
In the 1997-2010 transition matrix at GPGM, large part of forest has been transformed into bushland and built up area. The main driver of forest change was clearing of vegetation for construction of roads, construction of buildings for project personnel and the development of other infrastructure. This indicates that, the large part of the forest with tree species of more than 10% crown cover and the minimum height of 3 m could be removed hence changed to bushland.
Also, the forest cover was cleared off to pave way for construction of project infrastructures and accommodation for project personnel. At GGM during 1997-2010 and 2010-2017 periods, the main driver of forest land cover change was also bushland. This indicates that, these changes are associated with intense deforestation resulting from transforming forest to bushland.

Land Cover Projection (Future Land Cover) for the Year 2027
The findings have further revealed that, in both studied sites, it was projected that in the year 2027, forest will remain to be the most dominant land cover type while water bodies will be the least land cover type. Although, the forest is projected to be the most dominant land cover type in 2027, it will not reach the total area it covered in 1997. For example, at GPGM in 1997 the area covered by forest was 52.1% and, in the projection of 2027, it will be 34.7% and at GGM the area covered by forest in 1997 was 49.7% and in the projection, it will be 43.7%.
In addition, uncontrolled gold mining operations led to loss of biodiversity and ecosystem services. Notwithstanding the increased mining industry due to the government effort on promoting investment to achieve economic development of Tanzania, there should be strong enforcement of the laws to ensure that mining companies rehabilitate the degraded mining sites to the level close to its original condition. Also, it is important for the government and mining operators to enhance monitoring measures in the revegetated areas to avoid trespassing. Thus, there is a need of involving local people in the management of the rehabilitated areas.

Accuracy Assessment for Both GPGM and GGM
In this study, the accuracy assessment classification was done with reference to the raw satellite images. The aaccuracy assessment classified for GPGM indicated that the overall accuracy is 74.7% and the assessment of land cover model using Kappa statistics is 70.2%. Whereas, the accuracy assessment classified for GGM indicated that the overall accuracy is 71.4% and Kappa statistics is 66.1%.
The overall accuracy values for GPGM and GGM falls within standards suggested by [36] and the Kappa statistics were interpreted as being in substantial agreement than by chance alone [8] [29] [37]. This is confirmed in the study conducted by [38] who provided land use and land cover accuracy of 82.00% Journal of Geographic Information System and Kappa (K) statistics of 77.02% which is acceptable in both accuracy total and Kappa statistics. According to [29] the overall accuracy is interpreted as almost perfect agreement and Kappa Statistics accuracy is interpreted as substantial agreement. In addition, the study conducted by [39] achieved overall classification accuracies of 95.32% and 95.13% and overall kappa statistics of 0.9237 and 0.9070 respectively for the classification of 1992 and 2012 images. Moreover, the change detection approach conducted by [23] yielded a reliable change map with an overall accuracy of 83.3% and a kappa coefficient of 0.79 which is acceptable.

Conclusion
The study aimed at monitoring land cover changes of GPGM and GGM toward 1997-2017 and prediction of 2027 by using remote sensing and GIS technologies.

Area for Further Research
Remote sensing and GIS technologies used in this study serve as a reference point for further studies; the following area needs further research: Determining and quantifying the health of vegetation in the gold mining areas using Normalized Difference Vegetation Index (NDVI) by measuring the difference between near-infrared and red light.