Spatio-Temporal Distribution of Southeastern Community Forests in Togo ( West Africa )

The present study aimed to provide the basics needed to reconcile the fight against poverty with the need of a good management of community resources in the context of local and sustainable development. It focused on the identification of geographic location and spatio-temporal dynamics of these resources in the southeastern part of Togo. The methodological approach was based on the spatial analysis of the area for the years 1988, 2000 and 2018 by combining optical imagery from Landsat TM, ETM+ and Sentinel-2A MSI satellites with radar imagery from Sentinel-1A IW/GRDH. Spatial analysis showed a fragmented spatial structure undergoing major changes for the period 1988-2018. Plantations, riparian formations, agglomerations, water bodies and wetlands increased respectively by 4.61%, 2.09%, 1.07%, 0.43% and 0.35% annually, while forests, savannahs, crops and fallow lands decreased annually by 0.71%, 5.87% and 1.97%. For the “forests” class, seven community forests geographically organized in three sub-groups were identified and mapped. The analysis of their areas pointed to the fact that out of 667 ha of community forests in 1988, only 415 ha remain in 2018, which means a loss of 37.78% in forest areas over the 30 years, or an annual deforestation rate of 1.64%. The different spatial changes observed could be attributed to several unsustainable human activities. The land use maps for the years 1988, 2000 and 2018 will contribute to resource localization and protection in sensitive areas or, in other words, to the integrated and rational management of these resources. The different identified and mapped community forests could serve as management units for managers in developing their management plans. How to cite this paper: Konko, Y., Rudant, J.P., Akpamou, G.K., Noumonvi, K.D. and Kokou, K. (2018) Spatio-Temporal Distribution of Southeastern Community Forests in Togo (West Africa). Journal of Geoscience and Environment Protection, 6, 51-65. https://doi.org/10.4236/gep.2018.67004 Received: May 30, 2018 Accepted: July 8, 2018 Published: July 11, 2018 Copyright © 2018 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
Forest management planning in Togo began since the colonial era and resulted in the creation of 83 protected areas [1].Over time, these protected areas have unfortunately been managed unilaterally and repressively by the forestry administration, leading to an alienation of the rights of local communities.This situation hindered the application of regulatory texts governing the conservation of biological resources, and led to a loss in forest areas because of human activities [2] [3] [4].
In order to correct this situation and involve local populations in the conservation and management of forest resources, several policies, strategies and programs have been developed leading for instance to the introduction of the concept of Community Forests (CF) through the law No. 011 of March 13, 2007 on the decentralization and local liberty policy [5].This law proposed a favorable legal framework that allows local communities to manage, in a consensual manner with the administration, individual properties by offering them the opportunity to enjoy all the benefits they can derive from it.Individuals are defined by the Forest Code of 2008, as natural or legal persons, rural or basic communities or groups that do not fall into the category of local authorities.The CF approach is effective by empowering communities to undertake economic activities aimed at reducing poverty, improving living conditions and ensuring local development [6] [7].
The approach has already proven to be effective, as more than ten CFs have been implemented.The extension of the approach throughout the Togolese territory requires the localization of the resources, their characterization and the identification of opportunities they can offer for local development.
In southeastern Togo, and more specifically in the Lower Mono Valley (LMV), the characteristic vegetation formations of wetlands, including mangroves, semi-deciduous forests and floodable savannahs with Mitragyna inermis, are of major importance.The ecological conditions (important hydrographic network and clay soils, flooded almost every year), making the soils difficult to farm, are still assets for the durability of these community forests.However, the area is highly vulnerable since it is subject to deforestation driven by the need of fuelwood of an increasing population.The constant search for appropriate solutions to the problem of recurrent poverty is pushing local communities to exert uncontrolled pressure on local natural resources.Therefore, the adequacy between the population growth rate and the degradation of forest resources raises questions about the distribution of community forest resources in the LMV and changes they have undergone during the last three decades.
To answer these questions, the present study is based on the assumption that local communities of the LMV still have the capacity to protect a large part of the forest resources of their land, despite the problems of population growth and resource degradation that they are facing.Based on remote sensing techniques, such as the combination of optical and Radar satellite imagery, this study aims to provide the basic elements needed to reconcile the fight against poverty and the need of a good management of CFs of the LMV through the principles of sustainable forest management.Specifically, it involves studying and monitoring the spatio-temporal dynamics of the LMV community resources using satellite imagery.

Study Area
The Lower Mono Valley (LMV) is located south of the Mono Transboundary Biosphere Reserve (http://www.unesco.org/) in southeastern Togo.With an average area of 107,486 ha, it lies between latitudes 6˚14'32,921" and 6˚41'24,013" North and longitudes 1˚24'6018" and 1˚48'30,631" East (Figure 1).The terrain is slightly rugged, but dominated by vast plains.The pedology is characterized by poorly developed soils, vertisols, and soils rich in iron sesquioxides [8].In terms of climate, the LMV benefits from the Guinean tropical regime characterized by an alternation of two rainy seasons and two dry seasons, well distinct and of unequal lengths.The annual rainfall varies between 1000 and 1400 mm/year.The average temperature is generally high, around 27˚C in the coastal zone [5].

Data Acquisition
The data used in this study are satellite data obtained through passive (optical) and active (Radar) remote sensing.Optical images (Landsat TM, ETM+, https://landsat.usgs.gov/;Sentinel-2A MSI, https://sentinel.esa.int/) were selected based on their availability and a cloud cover rate of less than 5%.Apart from differences in resolution between the three sensors (Table 1), there are also differences in bands.In fact, TM consists of seven bands; ETM+ has a panchromatic band in addition to those of TM; MSI has 13 bands.Whereas Landsat TM and ETM+ include a thermal band, Sentinel-2 MSI does not include a thermal band.The choice of Radar images (Sentinel-1A, https://sentinel.esa.int/) was based on seasons of the year (short and long rainy seasons, small and long dry seasons), given the fact that the Radar echo can be influenced by surface properties (Table 2).

Pre-Processing of Optical Images
The preprocessing of optical images was done under the SNAP (Sentinel Application Platform) software [9].It included contrast enhancement, resampling and mosaicking of images.Subsequently, in order to improve the discrimination of  land use units through photointerpretation, two indices were calculated from the satellite images: The Normalized Difference Vegetation Index (NDVI) [10] and the Normalized Difference Water index (NDWI) [11].These two indices were calculated as in Equations ( 1) and ( 2) where NIR is the Near Infrared channel, MIR is the mid wave infrared channel and R is the red channel.

Pre-Processing of Radar Images
Sentinel-1A Radar images were pre-processed with Sentinel 1 software.Pretreatments included geometric correction, radiometric calibration, speckles filtering using a Lee filter [12] and orthorectification.The orthorectification was performed using the Range Doppler Terrain Correction algorithm [9] and an ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) Digital Terrain Model.Afterwards, complementary images such as color composite images and resultant images (VV + VH), both filtered and at full spatial resolution, were produced and converted to decibel [13] to allow an easier discrimination of land use classes.

Visual Interpretation
The optical and radar images as well as the different indices were interpreted visually based on spectral signatures, field visits and Google Earth images.The 1956 African vegetation typology proposed by Yangambi [14] was used, with some adaptation, for naming the different land uses.Following the image interpretation, Sentinel 2A satellite sensors proved to have a better scene rendering compared to Landsat (TM/ETM +) sensors.For the filtered radar images, the color composite images allowed a better identification of land use classes compared to the VV and VH polarized images and resultant images (VV + VH).The visual interpretation of both optical and radar images made it possible to define an interpretation key for the classification of optical images.

Classification of Optical Images
The interpretation key helped define training zones by creating ROIs (Region of Interest, defining land use classes) under OrfeoToolBox (OTB).Subsequently, a supervised classification was performed using the Support Vector Machine (SVM) algorithm.This algorithm was chosen for its efficiency in classifying complex landscapes and producing good results [15].It also performs better than many other methods of classification of satellite images [16] [17].SVM is a statistical learning classification technique.SVM focuses classification decisions on the boundary between classes instead of their mean and variance values [18].

Quality Control of the Classification
The validation of the classification was carried out through 160 control points randomly selected and distinct from the ROIs, 20 points per land use class.
Ground truth was obtained through field visits and Google Earth images.By confronting the classified images with ground truth, confusion matrices, Kappa indices (K) and the global precision were computed as in [19], for the years 1988, 2000 and 2017.The kappa index is used in image classification as a metrics of the agreement between the classification results and ground truth [20].According to the scale of Landis and Koch [21], the classification is considered excellent if K is greater than 0.8, good if K is between 0.8 and 0.6, moderate when K ranges from 0.6 to 0.2 and bad when K is less than 0.2.Shapefiles and area of each landuse were obtained through the vectorization of the classified images.

Spatio-Temporal Land Use Change
The classified and validated optical images were synchronized under the SNAP software.This synchronization allowed the visualization and mapping of some of the changes that occurred between 1988 and 2018.Subsequently, a quantitative analysis was performed based on the annual evolution rate (T annual , Equation (3)) [22].S 1 and S 2 represent the area of a given land use, for the initial date and the final date respectively; d is the number of years between the two dates.

( ) ( )
For the CFs, the deforestation rate (θ) was considered instead, because it is specific to forests and takes into account cleared land.The standardized formula proposed by Puyravaud [23] was applied as in Equation (4).A 1 and A 2 represent forest areas respectively for the initial year and the final year.T 1 and T 2 correspond to the exact acquisition dates of the satellite images respectively for the initial year and the final year.

Quality Control of the Classification
The control of the classification of TM and ETM+ optical images gave an overall accuracy of 84.54% and 86.31% respectively.The Kappa index (K) calculated was 0.81 and 0.83 for the years 1988 and 2000.The control of the classification of the MSI image of 2018 showed that 147 out of the 160 randomly selected points were accurately assigned to their actual land use.This generated an overall accuracy of 91.57% and a Kappa index of 90.03% (Table 3).

Land Use Maps
Figure 2 shows maps obtained from the classification of optical images for the three dates (1988, 2000 and 2018).In total, eight land cover classes were mapped.These are forests (forest patches, sacred forests, community forests, etc.), savannahs (wooded savannahs, savannahs with few trees, grassland savannahs), plantations (palm groves, coconut groves, teak plantations), riparian vegetation, body of water, swamps, urban areas, crops and fallows.

Qualitative Analysis of Land Use Change
The comparison of the classified images of 1988 and 2018 helped perform a qualitative analysis, in order to identify changes in land use and areas where they have occurred.Figure 3 is an illustration of some of the major changes.Overall, urban areas and plantations increased in area at the expense of forests and savannahs.

Quantitative Analysis of Land Use Change
The quantitative analysis assessed areas of land use change for the study period.
In 1988, the landscape of the region was dominated by crops and fallows (36.85%), savannahs (28.08%) and plantations (21.49%).Between 1988 and 2018, there was a decrease in the area of savannahs, crops and fallows, while an increase was identified for plantations and riparian formations (Table 4).The annual growth rate during this period for plantations, riparian formations, urban areas, body of water and marshy areas was 4.61%, 2.09%, 1.07%, 0.43% and 0.35% respectively.The annual decrease rate was 5.87%, 1.97% and 0.71% for savannahs, crops and fallows and forests respectively (Table 5).In 2018, the Journal of Geoscience and Environment Protection

Spatial Distribution of Community Forests
Seven (7) semi-deciduous CFs geographically grouped into three (3) subgroups were identified and mapped.The first one is the Avélèbè community forest located on the eastern side of the town of Aveve near the Mono river.The second one includes the community forests of Amévo, Dougbanavé, Fontan and Mambui located on the western side of the city of Aveve.The third one includes the community forests of Akissa and Tètè Condji located southeast of the city of Adame near the Mono river (Figure 4).

Spatio-Temporal Dynamics of Community Forests
The change analysis of community forests (Figure 5) revealed that out of 667 ha  i.e. 37.78% of loss in forest area over 30 years.The corresponding annual deforestation rate is 1.64% (Table 6).Tètè Condji, Amévo and Manbui forests have the highest annual deforestation rates compared to the other forests.

Discussion
This study was based on multi-source satellite images.The Sentinel 2A MSI images with a spatial resolution of bands ranging from 10 to 60 m allowed the identification and characterization of the current state of the resources of South-East Togo while Landsat TM and ETM+ images of 30 m spatial resolution were used as historical data to allow spatio-temporal monitoring.Although these images came from different sensors, it was very beneficial to bring them together for the characterization and the spatio-temporal monitoring of the landscape.
However, it is important to note that the MSI sensor rendered scenes better than the TM and ETM+ sensors [18].
The results of visual interpretation of the sentinel-1A radar images helped to obtain complementary information relative to the landscape, according to the seasons of the year and the polarization (VV, VH) of the Sentinel-1A platform.
However, some difficulties were encountered in the discrimination of land use units based on the VV and VH polarized images at full spatial resolution.In fact, since these images were not filtered, the presence of speckles did not allow an easy discrimination of the land use units.On the other hand, on the filtered images, the color composite images allowed a better identification of land use classes compared to VV and VH polarized images and resultant images (VV + VH).
The coupling of Optical and Radar images was highly beneficial in this study owing to their complementarity.For instance, Radar images allowed a better identification of forest cover, whereas optical images performed well in identifying urban areas.We therefore strongly recommend the coupling of both sources of data.The quality control of the classification of optical images TM, ETM+ and MSI gave an overall accuracy of 84.54%, 86.31% and 91.57% respectively.The calculated Kappa index (K) was 0.81, 0.83 and 0.9 for the years 1988, 2000 and 2018 respectively.Therefore, based on the scale of Landis and Koch [21], the classification can be said excellent (K > 0.8) and the results are exploitable.
The spatial analysis revealed a changing spatial structure.The increase in area of water bodies, swamps and riparian formations could be explained by the fact that a considerable part of the Lower Mono Valley is often flooded by floodwaters from the Mono River during rainy seasons.In addition, the pedology of the area predominated by vertisols [8] does not favor enough infiltration of rainwater, leading to stagnation and runoff and causing flood and erosion.The expansion of urban areas reflects the construction of new housings and infrastructures.
The study of [24] on this region also showed an annual growth of 2.01% for urban areas and 1.02% for water bodies.Concerning the decrease in area of savannahs, we observed in the southeastern part of the study area that savannahs were grazed, which led to soil compaction by cattle and their transformation into

Conclusion
The present study aimed to provide the basics needed to reconcile the fight against poverty with the need of a good management of community resources in In other terms, they will be used for an integrated and rational management of resources in the Lower Mono Valley.The different community forests that were identified and mapped represent potential management units on which managers can rely for the development of management plans.

Figure 1 .
Figure 1.Geographic location of the lower mono valley.

Figure 3 .
Figure 3. Overview of some of the major changes between 1988 and 2018.

Figure 4 .
Figure 4. Geographic location of community forests of the LMV.

Figure 5 .
Figure 5. Breakdown of community forest areas in the LMV from 1988 to 2018.
marshy and/or floodable zones.In the center of the study area, savannahs are being converted into plantation (e.g.coconut and palm groves).The decrease in area of crops and fallows is explained by the rural exodus phenomenon, which led to farmland abandonment and subsequent conversion into plantations.The CFs that were identified and mapped represent majestic examples of forest in the area.They provide ecosystem services (timber and non-timber forest products collection, regulatory, cultural, and self-maintenance services) which are essential for local populations.This said, the CFs of the Lower Mono Valley constitute a development pool in economic, social, cultural and environmental terms.However, the area analysis revealed that they are being overexploited by local people, leading to a rapid decline.In fact, out of 667 ha of CF in 1988, only 415 ha remain in 2018, i.e. 37.78% of lost in forest area, which means an annual rate of deforestation of 1.64%.At this alarming rate of deforestation, there would be no more community forests in half a century if no action is taken.The proximity of the community forests to two identified agglomeration complexes also suggests that they once formed the same forest unit, which might have been fragmented and degraded because of human activities.The main factors of degradation include the uncontrolled exploitation of wood resources (timber, fuelwood, charcoal) and non-wood products (animals, medicinal plants, etc.), bush fires for the practice of agriculture and hunting.Several authors have proved that the regular return of fire does not favor the reconstitution of the original vegetation[25] [26] [27].This leads to a loss of biodiversity and consequently an important change in the structure and composition of the vegetation through the installation of herbaceous species[28] [29].The study of[30] on the sacred forests of the lower Mono Valley has also shown that conflicts in this area (population growth, uncontrolled expansion of housing, erosion of religious beliefs and current weakness of power of religious leaders) weakened resource Y. Konko et al.DOI: 10.4236/gep.2018.6700463 Journal of Geoscience and Environment Protection management systems that local people developed and conserved throughout generations.For these reasons, the resources suffer rapid and massive degradation leading to the reduction in their area or their near disappearance.
the context of local and sustainable development.It focused on the identification of geographic location and spatio-temporal dynamics of these resources in the southeastern part of Togo.Based on remote sensing techniques, this study helped to map the different land uses for the years 1988, 2000 and 2018.These maps will contribute to resource localization and preservation of sensitive areas.

Table 2 .
Characteristics of radar satellite images.

Table 3 .
Confusion matrix for the quality control of the classification.

Table 4 .
Changes in area of the different land uses.

Table 5 .
Annual rate of change (ARC) of the different land use classes.

Table 6 .
Annual rate of deforestation in community forests.