Monitoring Forest Cover Change and Fragmentation Using Remote Sensing and Landscape Metrics in Nyungwe-Kibira Park

The objective of this study was to evaluate forest cover change and forest degradation in Nyungwe-Kibira Park, a natural reserve straddling Rwanda and Burundi from 1986 to 2015. Landsat TM, ETM+ and 8OLI images of 30 m spatial resolution were used as primary datasets. Geographic Information System (GIS) techniques were used for forest cover mapping and landscape metrics were calculated by using FRAGSTATS software. Classification and change analysis of forest cover type and landscape patterns analysis were carried out. In addition, to analyze the correlated external disturbances, the buffer zone of 5 Km was delineated outside the boundary of Nyungwe-Kibira Park. The results revealed that in among 5 land cover classes considered within the Park, the dominant one was dense forest class covering over 70% of the entire Park area while in the buffer zone cultivated and open land dominated at over 90% between the years 1986 and 2015. Change detection highlighted that within Nyungwe-Kibira forest, approximately 0.27% (4.97 Km) of forest cover was cleared while 0.07% (1.22 Km) was regenerated annually. In the buffer zone, the annual cleared forest cover was about 0.76% (13.02 Km). The five landscape indices chosen at class level indicated a considerable fragmentation of forest inside the Park and the highest fragmentation in the buffer zone. Indeed, these results shed a bleak image over the future of the Nyungwe-Kibira forest that should be helpful for the policy-makers and managers of these natural parks to establish adequate policies to mitigate the forest loss and degradation by implementing quick and effective solutions. How to cite this paper: Kayiranga, A., Kurban, A., Ndayisaba, F., Nahayo, L., Karamage, F., Ablekim, A., Li, H.W. and Ilniyaz, O. (2016) Monitoring Forest Cover Change and Fragmentation Using Remote Sensing and Landscape Metrics in Nyungwe-Kibira Park. Journal of Geoscience and Environment Protection, 4, 13-33. http://dx.doi.org/10.4236/gep.2016.411003 Received: September 30, 2016 Accepted: November 1, 2016 Published: November 4, 2016 Copyright © 2016 by authors 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
Tropical deforestation has been claimed as the primary contributor to global environmental change [1] which substantially threatens biodiversity, climate and livelihoods [2].Deforestation in African tropics has recently accounted for over 23% of forest loss worldwide [3].Protected Areas (PAs) are found as a cornerstone of forest conservation policy in the developing countries [4].The establishment of PAs is the main means to conserve forests and the associated ecological benefits including biodiversity habitats, hydrological sources and carbon sequestrations [5] [6].In principle, PAs are established prior to stem forests clearing and degradation by restricting land-use changes and extractive activities within their borders [5].However, in some countries, PAs are often established after forests have been cleared and degraded or yet the restrictions may not have been enforced because of political issues, financial insufficiency of their surrounding community, uncertainty of land tenure and conflicts within local communities.These factors have been inventoried among the reasons of large deforestation and forest degradation in Africa [7].
In central Africa, especially in Albertine Rift zone, deforestation is reported to have increased in several parks like in Uganda where, in the last-half century, the forest around Kibale National Park was cleared in attribution of charcoal production and the formerly impenetrable forest around Bwindi in South-West of Uganda vanished due to agricultural expansion [8].In the Democratic Republic of Congo (DRC), thousands of hectares of forest have been cleared for mining activities, charcoal production and agricultural activities [9].Rwanda and Burundi experienced similar cases of forest degradation, where the PAs have decreased in size as a result of agricultural expansion and ambiguously nudged forest boundaries [10].As highlighted by Haino, M et al., [11], rapid population growth is correlated with the forest losses over the whole Albertine Rift region [12].In the southern Albertine Rift landscape, the large forest of Nyungwe-Kibira Park, which is rich in biodiversity, was reported to be protection sensitive [13].However, as the forest is enclosed by public and private landholdings, it has been reported that agricultural expansion, legal and illegal timber harvesting, unclear land tenure system are the main constraints to this forest conservation [14].Despite a number of studies which have been carried out focusing on forest loss in the Nyungwe-Kibira corridor [14] [15], little to zero ranging emphasis has been put on highlighting the spatial and temporal patterns of forest cover changes within and outside the boundaries, especially within the last three decades.Therefore, historical forest cover changes and landscape scale status could be as-sessed to provide critical information for future planning [16].Hence, spatial and temporal pattern research is needed to assess the net effects on the protected areas on both forest cover change and the analysis of its landscape scale generally [17].The main issue here is to identify the forest cover change specifically the rate of deforestation and degradation and rate of forest regeneration.Analysis of landscape scale is also necessary to distinguish the state of patches of forest and understand the spatial temporal complexity of forest cover fragments.Satellite-based tools have been proven as the appropriate way to study forest cover changes due to their capacity to cover a large-scale area and to transcend country borders when the study is carried out on a transboundary scale.Additionally, it provides the time-series information from the start to end points of the period [18].In that framework, the generation of forest cover change cartography from the satellite image processing is one of the most substantial and widespread remote sensing applications [19].From this cartography, the possibility to obtain and analyze the spatial patterns of forest cover in a region is higher.Satellite based remotely sensed data have been widely proven to be most useful in landscape ecology studies and in understanding the ecological processes with strong spatial temporal components [20].Landsat data have been mostly used for determining forest cover and measuring forest cover changes and the rate of change [21] [22].Therefore, the objectives of this study are 1) to describe the land cover changes in Nyungwe-Kibira forest from 1986 to 2015 2) to monitor the forest cover changes and 3) to analyze landscape dynamics within and outside Nyungwe-Kibira Park boundary.

Study Area
Nyungwe-Kibira forest (Figure 1) was the first gazetted forest reserve by the Belgian colonial administration in 1933 combining Nyungwe forest and the Kibira forest.It lies within the border at the South-West of Rwanda and stretches to the Kibira National Park in the North-West of Burundi and is counted among the located tropical mountain rainforests of East-Africa [23].Between 1962 and1988 the massif of Nyungwe-Kibira forest was classified as the straddle of Nile and Congo River Basin and the principal Rwanda's water supply [24].In 1996 Kibira (North western Burundi) was officially declared a National Park while its conjoint Nyungwe in Rwanda was officially recognized as a National Park in 2005 [24] [25].Nyungwe-Kibira forest covers an approximate area of 1600 Km 2 and is characterized by a tropical climate with daily temperatures varying between 14.4˚C to 19.9˚C and a mean annual rainfall of 1744 mm [26] [27] with two seasons, dry and rainy [28].This forest is rich in both flora and fauna biodiversity; about 644 plant species, 98 species of mammals and 43 bird families, along more than 200 species, were identified within this forest as reported by the United Nations Food and Agriculture Organization (FAO) in its document repository 2010 [29].The study site is situated on the altitude ranging from 1012 to 2944 meters above sea level (m.a.s.l) (Figure 1). Figure 1.Location map of Nyungwe-Kibira forest.The extension area of the park was 1686.5 Km 2 whereas that of the buffer zone was 1700 Km 2 calculated for stratified 5 Km around the park.

Data Acquisition
Landsat Thematic Mapper (TM), Landsat Enhanced Thematic Mapper Plus (ETM+) and Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) remotely sensed data were used to map land cover change in Nyungwe-Kibira forest and its landscape.Thus, Landsat satellite images for 1986, 1990, 1995, 2000, 2010 and 2015 were acquired from the U.S Geological Survey (USGS) to evaluate the change in forest cover from 1986 to 2015 in 5 years time-steps.Good quality scenes from Landsat 172 path and 61 row were primarily targeted during this period.However, due to the obscure of clouds over the forests cover and high mountains on the surface of the study area, which may induce changes in surface features' spectral characteristics, some images were deemed unusable.Thus, images with the lowest cloud cover percentage were screened.Images either full of haze and clouds or affected by the Scan Line Corrector (SLC-off) defections were disqualified; hence, 2005 images and scenes prior to 1986 were precluded from the usable dataset in this particular study.Furthermore, in order to improve the discrimination during classification, only dry season images were selected since rainy season images in tropical areas are more complex to distinguish vegetation cover classes because of higher greenness [30].

Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) version
4 was used to derive the slopes and the aspects.The two tiles of DEM of 30 m spatial resolution have been downloaded from the CGIAR Consortium for Spatial Information (CGIAR-CSI)'s website, (http://srtm.csi.cgiar.org/selection/listImages.asp); one at Center point latitude of 2˚50' South and Longitude of 27˚50' East and the other at Latitude of 2˚50' South and Longitude of 32˚50' East.This SRTM data is available as 1 arc second (approx.30 m resolution) DEMs to the equator and 1 × 1 degrees tiles in WGS84 Geoid reference datum.This dataset has been chosen due to its finer spatial resolution that closely matches with the geo-referenced Landsat images and at comparable accuracies [30].Nyungwe-Kibira Park boundary was delineated based on the "Atlas change of our environment Program-UNEP" [31].And from this boundary, we created a 5 Km buffer zone to describe the surrounding landscape.The buffer zone has been established to evaluate the amount of pressures (socio-economic effects) coming from the forest's vicinity [32].This method has been applied in other studies and the 5 Km distance has been hypothesized as enough to analyze the external effects to the forest cover [32] [33] [34].

Image Processing
Image processing was performed using ENVI5.2 and ArcGis10.2into three phases: 1) pre-processing; 2) image classification; 3) post-classification and change detection.In the first phase radiometric and atmospheric corrections have been performed to correct atmospheric conditions from sensors' scanning errors as well as distortions from solar angle and sensors' angle in the Landsat images.DEM, Slope and Aspect were used for topographical correction to compensate the high topography and to eliminate classification errors.All images from each study year were clipped to match with the study area and buffer zone boundaries.
In the second phase, image classification and land cover change detection were performed.Supervised classification method was applied to six Landsat images (1986, 1990, 1995, 2000, 2010 and 2015) using the Maximum likelihood Classification method (MLC) a well-known parametric classifier for supervised classification [35].Training samples were independently selected from each image in order to produce related spectral signatures that depicted five classes in the final classification (dense forest, dispersed forest and shrub, water, perennial croplands (tea plantation) and other agricultural land including cultivated and open lands).The third phase contained the post classification process in which kernel size of 7 × 7 majority algorithm was applied and used to smooth and eliminate noise in classes.Moreover, change detection was performed to assess the changes through the entire land cover during the study period for both the Park area and buffer zone.Land cover maps were derived from each year of study.To address the concern of detecting forest cover change in Nyungwe-Kibira Park from 1986 to 2015, land cover types of interest were forested areas (dense forest and dispersed forest).In this case, the forest classes were combined into one class called forest cover and was assigned to the values of "1" and the combination of all other classes was named "Non-forest" and assigned to the values of "0" to create a binary image in ArcGIS 10.2 where the processes were applied to all classified images .Image differencing was then undertaken and recent scenes were subtracted from the older ones to detect spatial and temporal patterns of forest changes in time using Raster calculator tool developed in ArcGIS 10.2 to separate and statistically evaluate the forested area by non-forest areas and estimate forest losses (deforestation) and forest gains (regeneration).

Accuracy Assessment
The accuracy assessment was performed to determine the quality of information obtained from data analysis and classification processes of each images [36].The overall, kappa, producer and user accuracies of land cover and land use map were developed from Landsat data by comparing the land cover results with Google Earth (GE) based on Very High Resolution Satellite Images (VHRSI).Systematic initial reference points were geo-referenced on GE precisely according to the field knowledge of the study area.
Similar studies applied this particular approach in the recent past and it was proven successful [37].A total number of 320 points were stratified from the available Very high resolution images on GE (for 2010 and 2015) where at least 75 points for the classes covering big areas and 30 points for those covering small areas were selected.

Landscape Metrics Measurement and the Evaluation of Landscape Patterns at Class Level
To quantitatively analyze Nyungwe-Kibira forest and its surrounding landscape's forest fragmentation, the landscape metrics at landscape and class levels were applied.The five indices were calculated from FRAGSTAT4.2 for both the Park and buffer zone to further explore the degree of landscape cover types and spatial dispersion.In this study, we only focused on forest cover type which was the main focus of interest for this research.The indices were evaluated in each year based on the set of landscape range values as shown in Table 1.
The process of quantifying the habitat fragmentation involves both loss and change in habitat pattern.The fragmentation of landscape by habitat types would go from 100% to 0% [39].We found that the indices indicated in Table 1 were suitable to the landscape analysis at class level in conformity with the objectives of this study, consisting of the assessment of forest degradation inside Nyungwe-Kibira Park and the surrounding buffer zone.Mean-Proximity Index-PROX_MN (m) PROX_MN ≥ 0, PROX_MN = 0 if a patch has no neighbors of the same type at specific radius.Lower number of patches at a specified radius distance signifies the state of fragmentation.

Spatial Patterns of Land Cover Change
Land cover classification results showed that the dense forest (intact forest) over 70%

Accuracy and Validation
The accuracy assessment was systematically performed for each image; the results revealed that the overall accuracy and kappa coefficients represented for each classified image were greater than 75% for all images from 1986 to 2015.Error matrix was used to determine producer's accuracy and user's accuracy among the Land Cover (LC) classes.
The results from producer and user's accuracy were represented for variation among LC classes and periods, however, the confusion was caused by the increase of omission and commission errors only in Cultivated and open land which was at some instances confused with dispersed forest and shrub class.The classification accuracies are summarized in Table 3. respectively.However, within the buffer zone the increase in forest cover was only registered during the period between 2000 and 2010 by approximately 3.73% (62.73 Km 2 ) while the other time-series represented negative values as illustrated in Figure 3 and

Annual Rate of Forest Cover Change within the Park
The changes in forest cover within the Park were interchangeably represented by losses (deforestation) and gains (forest regeneration).The results from the analysis showed that between 1986 and 2015 the deforested and degraded area was approximately 10.68% and the regenerated forest was approximately 2.14% over the entire forest cover size inside Nyungwe-Kibira Park.This study has used Pyravaud's formula [40] that provides a standardized method to calculate the net annual rate of change from forest cover at different time periods while reducing confusion and misinterpretation.The annual rate of change abbreviated as "r" is extrapolated for each time period as follows: where A1 and A2 represent the area of forest cover at time one (t1) and time two (t2) respectively.
It has been revealed that an annual rate of 0.4% of initial forest cover was deforested/ degraded while an average of 0.07% of forest was regenerated.These statistical trends imply that 8% of the forest cover was cleared within Nyungwe-Kibira Park where about 0.27% of forest which is equal to 4.97 Km 2 was lost annually between 1986 and 2015.
Despite the greater loss found compared to the gain of forest cover inside the Park during the entire study period, this situation was balanced by the middle periods that represented higher gains than losses; for instance, we found that between 1990 and 1995 forest cover loss was 3.63% while forest gain was 7.82% and during the period between 2000 and 2010, the loss was 1.72% against 7.17% of gain, this expressed that the rate of deforestation was less than forest regeneration only in these mentioned periods as shown in Figure 3 and Figure 4.

Annual Rate of Forest Cover Change in Buffer Zone
Forest cover quantification within 5 km outside the park depicted the highest loss of approximately 22.11% over 29 years.This means that 0.76% of forest cover was lost annually, notwithstanding a moderate increase of 5.45% (63.72 Km 2 ) registered during the period from 2000 to 2010.The major losses incurred during the other two intervening periods (1995-2000 and 2010-2015) magnified the loss in total forest cover throughout the entire study period.Hence, the overall decline of 377.60 Km 2 in total forest cover originated from the mixed temporal patterns (see Figure 3 and Figure 4).

Estimating Degradation of Forest in Nyungwe-Kibira Park and the Surroundings, Including Landscape Analysis
The analysis per class level was done to evaluate the ecological landscape measurement  1 were calculated and the results were reported in Table 5 and Table 6 shown below.The results of forest cover type analyzed within and outside Nyungwe-Kibira Park represented a high variation along the metric indices which were used.Significant discrepancies were depicted between the forest covers inside and outside the park.NP (the number of patches) highly decreased in the period from 1986 to 2015 and these changes were proportional to the changes in Patch density and size in both inside and outside the park.MPS values found were greater within Nyungwe-kibira Park while found in the buffer zone were lower; signaling out the severe fragmentation outside and the slight fragmentation inside the park.These levels of forest degradation outside the park can also be emphasized by the patch density analysis.The results from this analysis indicated that patch density (PD) in forest within Park was lower at PD < 1; meaning that forest patch was not subdivided while in buffer zone due to high deforestation and degradation, PD is higher at PD > 1.This expresses that forest outside the Park boundary was highly subdivided into many patches which represent a high forest fragmentation as illustrated in Table 5 and Table 6 above.

Discussion
Analysis of forest cover changes and landscape dynamics in Nyungwe-Kibira Park through time series as well as Land Cover mapping revealed that the interchanges among the cover types were highly rendered by changes of forest cover from 1986 to 2015.The results revealed a moderate deforestation inside the Park where approximately 8% (144.42Km 2 ) was cleared within the Park.This expresses an annual loss of 4.97 Km 2 (0.27%) from 1986 to 2015 (Figure 3 and Figure 4).This loss in forest cover implying to an approximate annual loss of 0.76% equivalent to 17.06 Km 2 (Figure 3 and Figure 4).This high deforestation rate in the area surrounding Nyungwe-Kibira was associated to the intensive agriculture that was approximately 71.3% in 1986, before climbing to 93.1% in 2015.This analysis highlighted that conversions among the land cover types basically triggered the reduction of forest cover size both within and outside the park.Given that the major transitions mainly occurred among three classes during the entire study period, we analyzed the land cover transitions inside and outside the park; the results indicated that inside the park, dense forest as the preponde-  .The analysis of MPS indicated that forest inside the Park was slightly fragmented as reported by previous studies [18] [41] [42].Similar investigations highlighted that the greater MPS the lower degradation; conforming to the findings of this study especially the analysis of number of patches and edge density (Table 5).Normally due to the lowest change found in forest cover, without external circumstances, the inner forest inside the Park could have remained in its intact state.However NP revealed that forest was affected and slightly patched.
Despite the lowest level of forest fragmentation found inside the Park, in the buffer zone MPS were represented by the lower values which indicated the high level of forest fragmentation; hence, in buffer zone, the forest was highly fragmented; following the assumption that the greater the number of patches the greater subdivision of forest cover may have incurred (Table 6).Association between forest and other land cover types using Interspersion Juxtaposition Index (IJI) as suggested by Zhou, Q. et al. [43] as the most useful way to describe spatial association between land cover classes was applied in this study to assume the associability between forest patches, and the results indicated a high correlation between IJI and forest cover losses in both inside the park and in buffer zone.The correlation of forest change and IJI was high with R 2 = 0.756 and P-value = 0.0244 inside the park and R 2 = 0.6278 and P-value = 0.0602 in buffer zone.The IJI has changed apparently each year, which implied the spatial forest patches adjacency changes where the calculated IJI showed medium variations both inside the park and in buffer zone (Figure 5).
The aggregation of forest patches also was analyzed both inside the park and in buffer zone to check for the number of forest patches that may have clumped within a search radius and the threshold distance of 100 meters.PROX-MN is the most commonly used metric to evaluate the inter-patches distance and the level of isolation as reported by Bortoleto et al., [44].This method helped to verify the degradation of forest  within the park, precisely in Kamiranzovu wetland as detailed by Masozera, M.K, et al. [47]; the national tarred roads creating edges such as the 52 Km road linking the western and the southern provinces of Rwanda as well as the road linking Rwanda and Burundi, cutting through the forest [48].In addition, it is worthwhile to note the stretching of tea croplands located in the region along the parks' outskirts which cause enormous pressure to the forest reserve where Rwegura, Teza, Rwabidege, Gisakura, Gisovu and Gabegi in Mudasomwa tea croplands showed a significant penetration into the park boundary (Figure 2).As reported by NICOLE, D. et al. [49], the impacts of agricultural dominance and the encroachment exerted on forest reserves were clearly illustrated by the statistics where approximately 90% of the households in Nyungwe buffer zone depend on agriculture for livelihoods and per capita wood demand for fuel amounted to at least 14,881 kg day −1 which in turn affect the forest cover.Moreover, several analyses have been conducted on the sensitivity of forests with regards to the conflicts inflicted by the surrounding communities [50].
This study only dealt with the spatial and temporal pattern analysis and depicted the proximate forest changes within and outside the Park in a predetermined distance (buffer zone).The results discussed herein, land cover maps and statistical forest cover changes provided were not yet assessed and presented so far.Therefore, there is a great task for forest and Park management on both Rwandan and Burundian sides to effectively halt the degradation processes through joint collaboration and timely responses while promoting sustained research and mapping.

Conclusion
The paper.Therefore, this expresses the need for conservation efforts to be oriented towards this fact by promoting favorable measures such as the reduction of community dependencies and reforestation in the buffer zone.Finally, it is noted that the quantitative assessments of deforestation and degradation using spatial and temporal patterns constitute a valuable tool in terms of forest sustainability and management policies.
remained stable from 1986 to 2010 while the high forest cover reduction was recorded in 2015 falling to 63.73% inside Nyungwe-Kibira Park.Moreover, the dispersed forest cover became larger in 2015; expanding to cover about 18.73% of Nyungwe-Kibira Park total area compared to its former size (13%) in 1986.This degradation of forest cover increased relatively with the increase of cultivated and open lands which largely expanded from 1.5% in 1986 to 16.40% in 2015.Tea plantation, which was a predominant perennial crop in Nyungwe-Kibira corridor relatively decreased with the decrease of dense forest.The statistics representing changes within the Park are totally different from the results found within 5 kilometers in the landscape surrounding the Park (Buffer zone).Intensive agriculture land use was highly dominating the land cover (over 70% of occupancy from 1986 to 1995, growing to over 80% from 2000 to 2010 and reaching over 90% in 2015).Perennial crop was also characterized by an increasing cover area, expanding from 0.6% in 1986 to over 1% of the total area in 2015 as shown in Table 2 and Figure 2.

Figure 2 .
Figure 2. Thematic map showing land cover within Nyungwe-Kibira Park and within 5 Km of buffer zone in landscape outside the Park boundary from 1986 to 2015.

Figure 3 .
Figure 3. Forest cover change trends within Nyungwe-Kibira Park and in buffer zone (negative values = decrease, positive values = increase (values expressed in square Km).

Figure 4 .
Figure 4. Forest cover change results; change in forest cover was quantified within Park and in the buffer zone by differentiating the adjacent periods, the latest period minus the older in the entire study period.Areas of Forest lost (deforestation) were depicted in tan color.Forest extent gained (reforestation) was depicted in very light-green color.

within
Nyungwe-Kibira forest and in buffer zone of 5 Km outside the Park, using five indices to qualitatively and quantitatively analyze the spatial temporal patterns of change in forest cover type for 1986, 1990, 1995, 2000, 2010 and 2015.The indices illustrated in Table inside the Park has been converted into different classes among the Land Cover types found in the Nyungwe-Kibira Park after classification especially "open lands".The situation inside the Park was totally different from the outer parts (5Km buffer zone delineated outside the Park) which were characterized by extensive deforestation/degradation. Geo-statistical analysis indicated that approximately 22.11% (377.60 Km 2 ) of forest cover type in the buffer zone was lost (deforested and degraded) from 1986 to 2015;

Figure 5 .
Figure 5. Variation of the IJI representing the spatial adjacency of forest cover type patches in-side the park and in buffer zone.

Figure 6 .
Figure 6.Mean Proximity index (PROX_MN) values showing the aggregation or clumping level of forest patches within radius distance of 100 meters (m).(a) PROX-MN forest aggregation changes within Park and (b) showing the aggregation change in buffer zone.
forest cover changes and landscape dynamics within and outside Nyungwe-Kibira Park boundary were quantified and analyzed by means of remotely sensed data.This study provided deforestation/degradation and regeneration statistics at Nyungwe-Kibira Park and its buffer zone, in a period of29 years (1986-2015).The results indicated that 180.18 Km 2 of forest cover within the Park were lost while the gained or regenerated forest extended over 35.43 Km 2 .Also, a total of 412.08 Km 2 of forest cover were lost compared to 36.62 Km 2 of forest regenerated within 5 Km buffer zone delineated outside the Park.The analysis of forest degradation using fragmentation measurement methods indicated a slight degradation inside the Park and the highest degradation outside the buffer zone.The high deforestation and the way that forest cover was greatly fragmented in the buffer zone directly reflect the dominance of agriculture around the park.This study has only focused on the spatial and temporal changes of forest cover without thoroughly delving into the driving forces of forest degradation.Nevertheless, previous researches as discussed above, clearly established the direct links between anthropogenic activities and forest degradation and how agriculture takes the lead when land conversions are targeted to maintain livelihoods as demonstrated in this

Table 1 .
Landscape metrics at class level used to quantify forest cover in Nyungwe-Kibira corridor and the surrounding landscape (5 Km buffer zone).

Table 2 .
Spatial distribution of Land Cover change results for Nyungwe-Kibira Park between 1986 and 2015 within and Outside the Park (Buffer zone within 5 Kilometers).1:Dense Forest, 2: Dispersed and shrub, 3: Cultivated and Open land, 4: Perennial crop and 5: Water.

Table 3 .
Summary of accuracy assessment results.
1: Dense Forest, 2: Dispersed and shrub, 3: Cultivated and Open land, 4: Perennial crop and 5: Water, P.A: Producer's Accuracy in percentage (%) and U.A: User's Accuracy in percentage (%), O.A: Overall Accuracy and Kappa: Kappa coefficient.3.3.Land Cover Change Detection inNyungwe-Kibira Park and in Surrounding Landscape-Scale Within the Park and the outside landscape, the land cover classification highlighted different scenarios according to which changes were high in the dominant Land cover classes (Forest and cultivated/open land).From 1986 to 2015, dense forest reduced by approximately 14.32% within the Park while a decline of 5.04% was observed outside the Park (in buffer zone).The declines of forest was boasted by the great increase of cultivated or open lands where the open land and agriculture practices occurring within the Park's boundaries greatly increased up to 8.5% inside and highly increased up to 21.36% outside the Park.The statistical trends of the resulting land cover changes are illustrated in Table4.

Table 4 .
Land cover change detection within Nyungwe-Kibira Park and its surrounding buffer zone, (negative = decrease, positive values = increase).

Forest Covers Change Detection at Nyungwe-Kibira Corridor (Within and Outside the Park)
Approximately 8% (144.42Km 2 ) of forest was lost within Nyungwe and Kibira National Parks between 1986 and 2015 and a highest decrease of forest, about 22.12%, was recorded in the outer areas (buffer zone).The periods from 1990 to 1995 and 2000 to 2010 were characterized by the expansion of forest area inside the park.The magnitude of expansion was approximately evaluated at 4.29% and 5.48% (72.3 Km 2 and 92.5 Km 2 )

Table 6 .
Forest degradation estimates within the buffer zone outside Nyungwe-Kibira Park.
Due to landscape analysis at class level, the results revealed that the values of MPS increased proportionally with decrease of NP, lowest MPS values were calculated in the buffer zone opposite to the MPS found inside the Park during the whole study period