Land Use-Land Cover Change and Drivers of Deforestation in the Patako Protected Area (Center-West of Senegal)
Simon Sambou1, Anne Mette Lykke2, Hyacinthe Sambou1, 3, Idrissa Guiro1, Bienvenu Sambou1, Cheikh Mbow1, 3, 4
1Institute of Environmental Sciences, Faculty of Sciences and Technology, Cheikh Anta Diop University of Dakar. Dakar Fann, Senegal
2Department of Biological sciences, Aarhus University, Vejlsovej, Silkeborg, Denmark
3Laboratory of Teaching and Research in Geomatics, Polytechnic High School, Cheikh Anta Diop University of Dakar, Dakar-Fann, Sénégal
4World Agroforestry Centre (ICRAF), Nairobi, Kenya
To cite this article:
Simon Sambou, Anne Mette Lykke, Hyacinthe Sambou, Idrissa Guiro, Bienvenu Sambou, Cheikh Mbow. Land Use-Land Cover Change and Drivers of Deforestation in the Patako Protected Area (Center-West of Senegal). American Journal of Environmental Protection. Vol. 4, No. 6, 2015, pp. 306-317. doi: 10.11648/j.ajep.20150406.17
Abstract: This study focus on land use-land cover changes relation to the multidimensional characteristics of spatial and temporal patterns of land dynamics and social practices. The objectives are to characterize the dynamics of land-use and land-cover changes from 1992 to 2015, and to identify the drivers of deforestation and land degradation. Nine (09) types of classes were identified: humid vegetation (HV), savanna/wooded savanna (WS), savanna/shrub savanna (SSS), shrub savanna (SS), mangrove (MA), grass savanna (GS), cropland (CP), water (WA), and bare soil (BS). The process of change was analysed within a 5614.6 ha which includes the Patako protected area. The results have shown significant changes in land use-land cover for more than two decades (1992 to 2015). A high variability of land cover transfer was recorded between periods and class categories. At least, 55% of the HV changed to WS that was manifested by a loss of 16% during the first decade. The most important decline in area was recorded annually within the HV (4.8%), and the WS (2.6%), while the SSS experienced a substantial annual increase of 0.5%. The appearance of the SS (10.3%), during the last fifteen years, is a witness to degradation process as a result of anthropogenic pressure mainly for subsistence raisons. This pressure resulted to a conversion of small forested land to agricultural land on the boundary of the protected area. Tree cover represented 99.4% of the entire forest area in 1992, 97.4% in 2010, and slightly more (+1.2%) in 2015. According to this evolution, the annual rate of deforestation was estimated to 0.09% with an average of 0.11±0.08% for the three defined periods. The drivers leading to changes in land use-land cover stresses the complexity that is related to sustainable management of protected areas. Urgent action is necessary to reduce loss of biodiversity due to deforestation and land degradation. For this purpose, a particular attention must be paid to the implementation of the Patako management plan under a Reducing Emissions from Deforestation and Forest Degradation (REDD+) project.
Keywords: Land Use Change, Savanna, Saloum Delta, West Africa
Livelihoods and productive systems in African countries are highly related to vegetation changes [1,2,3]. Gathering landscape scale information to assess land use types, their evolution and deforestation rate requires earth observation methods that allow for large spatial coverage and historical tracking of spatial dynamics [4,5]. The use of remote sensing enables better design of management plans and thereby improved management of forest resources . Various applied remote sensing methods exist, but despite real progress in applied remote sensing, big challenges are faced in connecting maps with drivers of land use change. The link between image outputs and casual factors of change is sometimes fuzzy, nonlinear or incidental and requires steady disaggregated information on local human practices. Additionally, there is many intertwined connections between drivers with cascade effects that make it difficult to single-out individual factors. In many reports, the most described drivers are forest fires [7,1,8,9,10], pasture, farming land extension, fuelwood and timber extraction.
In this context, reversal of land degradation in semi-arid zones requires a clear understanding of the synergetic effect of natural and socio-economic factors . Several studies described these complex dynamics of land use-land cover change in savanna ecosystems of Senegal, e.g. in the groundnut cultivation zone . The rapid development of groundnut production and the subsequent economic crisis and migration are significant causes of change . Natural vegetation has been replaced by farm lands, where trees have been eliminated because of farming techniques and agronomist approaches that aim at modern intensification.
The conversion of forests to agricultural land and/or pasture induce in short and medium terms significant decrease in carbon stock, thereby promoting an increase in CO2 emissions in the atmosphere. Under UNDESERT (understanding and combating desertification to mitigate its impact on ecosystem services) program, a particular emphasis was given to the restoration of degraded lands, restoration of biodiversity and increasing carbon sequestration. A buffer zone has been created to mitigate population impact on Patako protected forest and several activities such as forestry and agroforestry have been implemented. Therefore understanding of drivers that lead to deforestation at local scale is essential for designing these interventions and develop best policies and measures. Addressing the drivers is also important as strengthening the conservation of forests and its resources are an important part of the development objectives of Senegal. This study focus on land use-land cover changes relation to the multidimensional characteristics of spatial and temporal patterns of land dynamics and social practices. The objectives are to characterize the dynamics of land-use and land-cover changes from 1992 to 2015 in a protected area, and to identify the drivers of deforestation and land degradation.
2. Material and Methods
2.1. Description of the Study Area
The study was conducted in Patako protected area and its surrounding agricultural areas in the center-west of Senegal (13º40’59.59’’N and 16º15’59.88’’W). Patako covers about 5614.6 hectares and was classified as a protected area in 1934. Patako was originally a tropical dry forest, but has changed into savanna mainly because of human impact such as fire and cutting. The vegetation has changed within the last 80 years.
The precipitation (600-800 mm) is unimodal with a short rainy season from June to October (4-5 months) and a dry season from October to May (8-9 months) . The interannual variability is high ranging from 562 mm to 813 mm within the last 10 years. Mean annual temperature is 26.58°C.
The area is dominated by ferruginous tropical soils, poorly evolved soils, hydromorphic soils on alluvial plains and gley-salted soils along valleys. Croplands with groundnuts, millet and maize occupy vast spaces, with mosaics of other land use-land cover types, mainly savannas, plantations of Anacardium occidentale, humid vegetation and vegetable gardens. There is a sharp border between the agrosystem with open cultivated land and the protected area. The agricultural lands have an average density of less than 2 trees/ha (range from 0 to 4 trees/ ha) and a deforestation rate of 3% per year .
2.2. Satellite Data Used
The whole area is covered by the Landsat scene (Path 205, Rows 50) for each of relevant dates. The selected images are acquired in the same period of the year during the dry season to reduce problems related to phrenology change in vegetation (exaggeration of the herbaceous plant cover) and soil humidity difference. The advantage of this period is to have a low cloud cover. Thus, the images used are: an old cover TM of 5 March 1992, ETM + of January 9, 2001 and December 28, 2010 and a recent Landsat 8 OLI (Operational Land Imager) of January 9, 2015. All images were taken from the database of the United State Geological Survey (USGS).
2.2.1. Images Processing and Classification
The first step of image processing was the visual analysis. It corresponds to a visual interpretation, ie the application of the conventional method of photo-interpretation to satellite image. This method involves identifying the different homogenous units . The interpretation was often facilitated by knowledge of the study area and field surveys. The images used for this study were already georeferenced (UTM WGS 84), with radiometric corrections and did not need geometric corrections. In order to avoid mismatches due to difference in the Landsat sensors and optimize the superposition of images for diachronic analysis, the data from 1992, 2001 and 2010 were orthorectified in reference to 2015 on the basis of seamarks points. The images were mosaicked and a smoothing of the dynamic performed. In order to harmonize and facilitate the superposition of images for diachronic analysis, the image of 1992 and 2001 have been readjusted compared to that of 2015. The wedging error (RMS = 0.124) is less than a Landsat pixel.
After wedging images, we created colored compositions by combining channels 5 for infrared [0.75-0.90 µm], 4 for red [0.63 to 0.69 µm] and 3 for green [0.52 to 0.60 µm]. In the colorimetric order red, green, blue were chosen for the images 1992, 2001 and 2010. As for the image 2015, the channels 6 for infrared [1.560 to 1.660 µm], 5 for red [0.845 -0.885 µm] and 4 for green [0.630 to 0.680 µm] were chosen. On this basis the study area was extracted in order to quantify the land use and its evolution over time.
The supervised classification method by maximum likelihood was used for all images. The choice of this method can be explained by the heterogeneity of land use in the study area and the probability of confusion for some components, namely bare soil and tans. The hierarchical classification starts with the determination of the number of classes and the choice of training plots to achieve a cartographic synthesis through successive selections of the best classified elements. The determination of the number of classes is based on the clustering performed upstream and information gathered in the field.
Furthermore, to eliminate isolated pixels or small unclassified areas and to homogenize the results of the thematic classification a filter of 3x3 pixels was applied from majority analysis. To check the validity and performance of the classifications, witness areas as well as reference data were collected in the field. The selected classes were also confronted with exogenous documents that provided information on the type of land use. These documents are thematic and synthesis maps. The analysis of the colored compositions allowed to identify at least nine (09) classes interpreted using the following abbreviation codes: humid vegetation (HV), savanna/wooded savanna (WS), savanna/shrub savanna (SSS), shrub savanna (SS), mangrove (MA), grass savanna (GS), cropland (CP), water (WA), and bare soil (BS). The type of vegetation was defined in the basis of the nomenclature of Yangambi  modified by .
2.2.2. Land Use-Land Cover Change Analysis
The result obtained from the supervised classification was exported to vector format in ARC GIS software. This transfer of map files helped to complete the processing and analysis of images and produce land use land cover maps. Changes detection were perform based on these maps. The crossing of the 1992-2001, 2001-2010 and 2010-2015 land use-land cover maps allowed to establish the maps of change. This cross was made based on classes’ codification and a code comparison between two maps. Statistics of detected changes were computed by comparing image class shifts between the four images 1992, 2001, 2010 and 2015.
The trends of changes were calculated based on typical rate of land cover transition for each category. Transition intensities were estimated relatively to the size of a land cover class in the initial year to its size in the subsequent year. These changes can be positive (increase or perspective of gains) or negative (decrease or perspective of losses). Three situations have been observed and taken into account in the analysis (e.g. land cover stability, conversion and modification). Land cover stability refers to all classes that have not been affected by changes during the studied period. The changes was expressed as conversion and/or modification. By conversion, we mean the transition from one category to another (e.g. from woodland to cropland), while modification represents a change within one land use or land cover category [18,19] due to changes in its physical or functional attributes (e.g. from savana to shrub savana to degraded savana). Matrices of land use-land cover transfer was generated using changes from one layer to another.
2.2.3. Annual Rate of Deforestation
Deforestation is explained as a conversion of forested areas to nonforest land use. This type of degradation notify a reduction of forest quality and a decrease of its overall cover due to native forest land clearing for different purposes . The annual rate of deforestation was calculated by comparing the area under forest cover at two different times. The period covered in estimating the rate of deforestation is from 1992 to 2015. Different methods have been used to estimate the annual rate of change. But in order to standardize the methods and avoid confusion,  devised a new formula which is used in the context of this study:
where S1 is the class area at the initial time (t1) and S2 is the class area at the ﬁnal time (t2) according to the land use-land cover maps.
2.3. Drivers of Deforestation
Characterizing the drivers of deforestation first involve their identification and grouping in different categories. We identified two main drivers (direct and indirect) which were divided into two categories. The direct drivers concerned natural effects that includes climate variability and change, flood and parasitic attack. The direct anthropogenic drivers are actions or factors that directly impact forest cover. They are considered as subsistence drivers. The indirect or underlying drivers of deforestation are complex interactions of demographic, economic, cultural and institutional that in combination intensify direct human causes. Semi-structured interviews were conducted with resource people for three target groups of local stakeholders (local population, local authorities and local technical staff of forest service). The interviewees were selected from 6 villages located in a 0–3 km radius around the protected area to help better understand the complexity of the drivers. Furthermore, on the basis of field knowledge observations were made into and outside Patako protected area for triangulating the obtained information.
3.1. Land Use-Land Cover from 1992 to 2015
The land use-land cover maps show that three characteristic classes were dominant during the studied period (1992-2015). These classes related to the CP in the periphery, the SSS and the WS in the Patako protected area. The land cover is highly dynamic (Table 1). In 1992, the periphery of the protected area was covered by a mixture of CP and SSS. Within the protected area the SSS represented (56.6%) while the WS and the HV occupied respectively 28.1% and 14.5%. A significant difference was observed between the latter and the rest of classes (WA, GS, CP and BS) as totaling less than 1%. These classes have increased slightly (1.2%) nine years after (in 2001). As for the WS a substantial increase of its size was observed (51.1%) at the expense of the SSS (44.3%) and the HV (3.4%). In opposition to the protected area the periphery was highly deforested and substituted to agricultural lands.
|WS||Woodland to wooded Savanna||1578.3||28.1||2869.5||51.1||1091.9||19.4||844.1||15.0|
|SSS||Savanna to shrub savanna||3177.9||56.6||2489.5||44.3||3688.4||65.7||3603.1||64.2|
In 2010 there is almost an elimination of the SSS in the periphery. Within the protected area the latter class cover remained the most significant type (65.7%). Therefore the WS and the HV covered respectively 19.4% and 3.5%. Marginal classes covered a total of 2.9% of the area of the protected area, an increase of more than 41% compared to 2001. Between 2001 and 2010, we noted the emergence of a new class, the SS representing (8.4%). It results from the fragmentation of the SSS and generally located in the center and north-eastern parts (Figure 2c).
Concerning 2015 there are some increase in land cover for many classes compared to 2001 and 2010 (Figure 2d). The SSS remained the most significant class (64.2%) despite the fragmentation it undergoes and which causes the increase of the SS (14.1%). This relative extension came at the expense of the WS (15.3%). The MA was characterized by a mean increase of 0.3% ranging from 0.2% to 0.4% between 1992 and 2015. This is also the case of the CP which substantially increased for 5 to 95 ha between 1992 in 2010 before decreasing to 48 ha in 2015. Other land cover (WA, GS, etc.) were poorly represented inside the protected area.
3.2. Land Use-Land Cover Change Analysis
Results from LULCC showed significant changes in land use and land cover. Between 1992 and 2001, many land use were converted to other types in the agrosystem. These changes affected mostly the SSS and GS (Figure 3a, b) and processes of change continued to occur during the period 2001-2010. Between 2010 and 2015 (Figure 3c) the agrosystem became more stable due to the almost total use of farmland resulting in a scarcity of fallows. As the changes in the agricultural lands were unidirectional we focused the analysis within the protected area.
Regarding the land use-land cover transfer between class categories it was found that generally modification was more important in term of frequency within the Patako protected area. This type of change affected mostly the HV from 1992 to 2001. At least 83% of its area moved to WS. The latter moved for 31.6% to the SSS and inversely for 36.9%. For this, a stability of 64.6% and 62.3% was respectively observed for these classes. Concerning the conversion the most important land use-land cover transfer was found in the CP and GS. More than 56.6% moved from the CP to GS and in the opposite direction 91.4%.
From 2001 to 2010 the HV became more stable (67.3%) compared to the first period. An average modification (19.3%) was observed towards the WS. This latter was the class that decreased in favor of the SSS (58.5%). This transfer caused a progressive shrinking of the class cover. As for the conversion it was characterized mainly by the GS which stability had become quasi null. 25.8% of the CP was abandoned and replaced by the SSS. Changes from vegetation cover to small classes (WA, GS, GS and CP) was generally weak (less than 10 %) but the transfer in the opposite direction was more significant despite the good stability (60.4%) of the CP.
From 2010 to 2015 the stability of some land use-land cover classes is rather good (> 50%) compared to others like the GS (3.4%). The vegetation cover modification affected more the WS and the SS which areas were respectively transfered towards the SSS for more than 36% and the 54%. The results from the matrix of land use-land cover transfer reveal that conversions were mainly observed in the GS and the CP. Despite its small area 54.5% of the GS was converted to SS while 18%, 12% change respectively towards the CP and the SSS. These changes explain the low stability of that class. As for the cultuvated land more than 35% were converted to the SSS and 15% to the SS.
When we consider the matrix of land use-land cover transfer from the overall period of study (1992 to 2015) we note that the HV and the WS experienced a low stability. For the first class at least 55% was modified to WS and 18% to SSS. For the second one the major modification took place towards the SSS as representing 72.2%. The transfer of area from the SSS towards the SS was done gradually going from 12.6% in 2010 to 14.9% in 2015 with an average of 20.2% between 1992 and 2015. It is also important to note that the intrusion of the arable lands generally took place in the border between the protected area and the agricultural zone especially where the farmers do not have a precise reference mark of delimitation.
Several changes occurred within the protected area (Figure 4). The most important decline in vegetation cover were held from 1992 to 2001 in the HV (16.0%). We also noted a substantial decline in the same period for the SSS (2.7%) while the WS experienced a gain of 6.6%. By cons, it recorded a loss of 10% of its area between 2001 and 2010 and 4.8% between 2010 and 2015. During the last period, substantial losses (0.4%) have been also demonstrated in the SSS. This decline was offset by a gain in area in the SS (10.3%). Globally the rate of change corresponding to these types showed a deline of 4.8%.yr-1in HV and 2.6%.yr-1in the WS, while the SSS recorded an extension of its area during the study period (1992-2015). Regarding the dynamics, MA stand out with significative changes. However, one of the most disturbing change was the influence of the agricultural zone of the protected area. Between 1992 and 2010, the CP experienced an increase of 16.3% before declining to 13.8% for the last five years (2010 and 2015).
3.3. Rate of Deforestation
Deforestation is the process by which forested land is replaced by another type of non-forested land, for example agricultural land. The total area of the Patako protected forest is estimated to 5614.6 ha. The conversion of forested land in this area is mainly connected with small-scale shifting agriculture or land clearing by sedentary landholders who want to expand their crop.
Taking into account the vegetation cover an annual rate of deforestation was estimated. The forested area decreased slightly during the study period. It covered 5578.4 ha (99.4%) in 1992, 5564.5 ha (99.1%) in 2001 and 5468.5 ha (97.4%) in 2010. In 2015 a slight increase (1.2%) was recorded compared to the coverage of 2010.
The results showed a differential annual rates of deforestation between the defined intervals of the studied period which are 0.03% from 1992 to 2001. During the following 9 years (2001-2010) the rate has increased to 0.19% before decreasing to 0.13% from 2010 to 2015. The average deforestation rate of the three periods is estimated to 0.11±0.08% /year, while the overall annual rate of the studied period (1992-2015) is 0.09%, which are globally a slow rate.
3.4. Drivers of Deforestation
Deforestation implies a long-term or permanent loss of forest cover and its transformation into another land use. There are a multitude of drivers that can be classified into two categories: direct and indirect. At the local level the drivers of deforestation are subjected to complex interactions with short or long-term impacts. Figure 5 summarizes the main drivers and their probable impacts in the study area.
The localities around the protected area did not have the same impact on forest degradation. Villages of both sides within a radius of 3 km were grouped into four locations. Those located in the western, northern and southern parts of Patako protected area had more facilities for access than the eastern part. This has resulted by an increase of the deforestation front line. Djikoye river is a natural barrier that reduced the flow of forest operators from the eastern side. The nearest protected area of Baria played an important role in meeting fuelwood and pasture for local communities.
4.1. Land Use-Land Cover from 1992 to 2015
The supervised classification approach resulted in good classification. The HV category presented both dense trees, thickets and aquatic invasive species (Typha domingensis and ferns). Pixels falling into this mixed category of vegetation had high similarity in spectral signature and were difficult to discriminate. It have been reported that degraded forest may be difficult to distinguish from certain types of grassland [22,23,24]. Such difficulties arise particularly in wetlands and riparian areas [25,26]. The 30 m resolution of images used may be a limit in spatial analysis and can influence the accuracy of the classification of small areas. Therefore field visits are invaluable to avoid lousy visual interpretation of images.
The LULC maps showed a progressive increase in land under cultivation over the years of study. The results have shown a rapid evolution of WS between 1992 and 2010. In 2001, the observed extend of this category might be explained by different phonological process of some species. In the other hand, this may be as a result of increased precipitation in 1999 and 2000 in the sudanian zone compared to earlier years. An increase of grass cover was also observed at the expense of the SSS. This category included areas currently under crop and land under preparation. In 2010 and 2015, it was highly difficult to discriminate GS in the agrosystem area because of the agricultural practices. Famers used to apply one time of weed control to avoid competition with grass until a certain height of millet crops. As for peanut cultivation, they usually weed two times and export straw after harvest (food for livestock). Such practices show a character of bare soil during the dry season.
4.2. Dynamics of Land Use-Land Cover Change
Land use and land cover changes (LULCC) are complex processes involving multiple driving forces that are location specific and context dependent. This complexity becomes even more important that the interactions between direct, indirect and underlying drivers are dominated by anthropogenic activities or influences.
Significant changes in the LULC have been observed in the agrosystem and within the protected area during the period 1992-2015. Some categories have been systematically transformed in varied proportions. These transformations concerned mainly land conversion. This type of change, from SSS or GS to CP, was quite intense and common in the studied landscape. Indeed, since the 1970s, many protected areas in the sudanian zone of Senegal have known an important spatial and temporal dynamics . However, despite the human pressure agricultural incursion within the protected area was negligible. In other hand land modification affected HV, WS, SSS and SS, but these changes can be subtle . Some authors supported that a detailed analysis of spatial dynamics allows a better knowledge of human activities impact on habitats [27,28]. Changes resulting from land use are generally complex and contributes to alter natural habitats.  described this alteration as leading to the establishment of a mosaic of savanna in Patako protected area.
Changes affected mostly the periphery of the protected area. It has been shown that almost all lands in the agrosystem were cultivated in 2010. An important increase of land converted to agricultural land was observed, which supports earlier studies from 1972 to 2002 . These changes resulted from an intensification of the use of animals for cultivation during many years ago in the Groundnut Basin. We argued that such increasing trend of agricultural production was to fulfill growing food demand and commercial opportunities. The cultivation techniques used promote land clearing. Likewise, canopy fragmentation had an impact on diameter structure its opening led to tree density decreasing and habitat fragmentation [31,32]. Combining field observations of tree cover, height and diameter structure helped to characterize subtle changes in land cover.
4.3. Rate of Land Cover Change
In the study area, the combination of drivers led to a fast rate of land cover change. These drivers constituted a wide variety of human activities including climate variability and change. The most affected vegetation types in Patako protected area were HV and WS. This could be related to the types of species exploited in these categories, mainly large-diameter trees for timber (Khaya senegalensis, Daniellia oliveri, Pterocarpus erinaceus and Cordyla pinnata). However, the SSS have been increasing after two decades despite of its serious modification at the middle and the northeast. Knowing how and why forest area changes over time is important for managing forests sustainably because such changes may result in long-term forest conversion to agriculture. From 1992 to 2015 the current study found an annual rate of deforestation which is globally very low compared to that estimated in the whole of Senegal by , approximately for the same period (0.47% annually). Compared to the protected area the agrosystem has a slightly higher rate of deforestation as estimated by . Nowadays, this area is dominated by a low tree density of Cordyla pinnata with aging population. Indeed, such aged stands of trees increased the vulnerability to fire and an absence of natural regeneration recruitment.
4.4. Drivers of Deforestation
This section reviews some of the most important drivers of deforestation that are either direct or indirect. Some direct drivers are natural. It is the case of climate variability and change and its repercussions on wetlands. Since the 70s, as other Sahelian countries in West Africa, Senegal has experienced a long period of rainfall deficits that affected natural ecosystems. Mean rainfall decreased by 20–49% between the periods 1931 – 1960 and 1968 –1997 . The declining rainfall has modified the hydrological regime of the rivers around Patako protected area. This change has resulted in the development of vegetation in the riverbeds. However, since the early 2000s there has been an increase in rainfall despite a sharp intermonth and interannual variability.
Floods became a cause of tree mortality in the valleys. It was perceived that the number of flood events has been increased on average from less than 2 per year before 1990 to more than 10 during the 2000s . The deposits of organic matter and silt drained by the water favored the establishment of an invasive perennial herbaceous plant (Typha domingensis). This invasive species has an important asphyxiating power that inhibits the development of other species. Such ecological conditions increased tree mortality along the river banks and created small openings on the canopy.
Parasite attacks, although weakly expressed are not negligible for describing natural drivers of deforestation. The genus Tapinanthus is the most common and often lead to tree mortality in the study area. As observed by , the susceptibility for trees to be infested was aided by senescence, fissure of bark and decrease of ability for resistance. This phenomenon is common in the local area where the Cordyla pinnata parkland consists of aging population.
An important part of direct drivers of deforestation was determined by human activities such as wood extraction, pasture, NTFP exploitation and the use of fire. Agricultural expansion was the main driving force of deforestation in the periphery. Contrary to the agrosystem, this problem is mostly dominated, in the core protected area, by wood harvesting which is driven by illegal logging.
Wood extraction has been demonstrated as one of the strong drivers of deforestation in the study area. The surrounding communities use wood as a primary source of energy for cooking. Unlike fuelwood, charcoal usually comes from trunks or large limbs and requires cutting trees. Even if charcoal may not always be the primary driver of deforestation, it can contribute significantly to degradation and destruction of the vegetation strata that have already been disturbed. As for the timber woods, there is a selective logging both among species and individuals into a cohort of one species. It was supported that this practice constituted a major agent of landscape modification due to the intensity use of machines for illegal timber extraction . The most threatened species in the savanna formations in relation to wood timber extraction are Pterocarpus erinaceus, Cordyla pinnata, Bombax costatum, etc. By cons, in the HV, Khaya senegalensis is abusively exploited due to the economic value of its wood. Even if proximity to Patako protected area is determinant in considering the problem of deforestation, it has been perceived that people engaged in illegal activities are mostly non-residents. Indeed, beyond the satisfaction of the local market in furniture, the attraction is also due to supply some fisheries areas for making boats.
These types of forest exploitation has been documented as a problem in many rural areas in Africa [37,38,39,40].  considered wood extraction for fuelwood and timber as contributing to forest degradation even without significant changes especially when the initial area was well bordered. As stated by  continuous and/or widespread human impacts can lead to vegetative changes. The damages can be an irreversible loss of species and carbon stock. Indeed, repeated actions of tree cutting and overgrazing resulted in open savanna types (SSS and DSS) with a significant frequency of stumps and cut tree stems. These elements have been considered as a key indicator of forest degradation [43,44,45]. The presence of several tracks on any direction in Patako protected area illustrated the pressure level on natural vegetation, as it has been shown in other protected areas in the region .
Human pressure was also manifested by the incursion of the Fulani transhumant pastoralists during the dry season. This incursion was considered a driving force behind degradation because requiring a significant animal load. The environmental problems caused by overgrazing in Patako area are the destruction of vegetation and soil erosion.
The predominant motivation of actors in NTFP harvesting is the quest for income generation. Exudates from plants mainly such as resins are processed to give local communities added value in the marketplace. As an example, human pressure on the perfumes products of Daniellia oliveri was manifested by unsustainable practices of harvesting. Other activities such as the use of plants (leaves, barks, roots) in traditional medicines can be a source of deforestation or degradation because they become more vulnerable to fire.
Unsustainable practices of forest exploitation are often a key variable of annual fire regimes in Patako area. Here, fires are caused accidentally or simply used by forest operators to burn grasses in order to facilitate access during the dry season. Even if its impact was considered negligible on adult trees , fire is one of the significant driving force in the dynamics of savannas [48,49].
Indirect drivers contributed mostly to increase deforestation in Patako area. Four categories of indirect drivers were described: demographic, economic, cultural and institutional.
Rapid growth in population and increase of socio-economic needs implied a cumulative pressure on forest resources. The study area recorded also an important flow of migration due to space availability and social context. This context explained the spatial distribution of land use . National agricultural policies in the 1960s and 1970s encouraged farmers to clear their fields of most trees to facilitate crop planting and harvesting. Then, trees density and diversity has declined as well because of the use of animal traction technologies for cultivation.
The economic value of some species used for timber and fuelwood has dramatically increased in recent years because of the local and external strong demand. In fact, the creation of weekly markets in areas around a radius of 8 km from the border of the protected area opened increasingly a way to fraudulent and opportunistic exploitation of forest resources. Some villages have a major advantage due to the presence of fertile valleys with fresh water, a favourable framework for market gardening. Other ones, where this opportunity did not exist people developed logging activities to supply economic needs in dry season. In this context,  considered the existence of the communities in many surrounded villages as closely linked to the use of forest resources.
Cultural drivers can be catalysts that influence flows entering into the protected areas. A significant part of the informants considered the natural resources in the protected area as belonging to local people. This representation has led to attitudes hardly compatible with a sustainable management of forest resources .
The institutional drivers of deforestation are multiform and complex. Two technical agents were mandated by the Forest Service for managing the Patako protected area. However, they are frequently faced with a problem of logistics and the fact their service is 22 km far from the site of intervention, justify the irregularities on forest survey. The lack of close monitoring paved the way for fraudulent exploitation. The consequence is that savannas support two weights (anthropogenic pressure and devastation caused by fires) causing a constant spatial dynamics. Underlying drivers must be treated with cautious because they have a cumulative effect on forest resource.
In this study, the combination of direct and indirect anthropogenic drivers has demonstrated the complexity of land cover change analysis. The results have shown significant changes in land use-land cover for more than two decades (1992 to 2015). These changes are characterized by two main situations which are either an increase or a decrease in the area of the specified classes. Within the Patako protected area, the rate of change varied from one type of land use to another. Globally, the area of HV and WS declined due probably to the selective logging of some species exploited for their wood timber. These species are subjected to a covetousness and constitute one of the main causes of the loss in vegetable cover. The SSS has experienced a small extension of its area despite some fragmentations in stands. The persistence of this fragmentation has contributed to the emergence of the SS which is a key indicator of human pressure. Some limitations have been shown using Landsat images of medium spatial resolution for discriminating riparian tree vegetation to invasive aquatic grasses. Detecting all stages of degradation was difficult because the scale at which the degradation takes place is often at sub-pixel resolution. This means that the nature of the degradation affects areas smaller than the detection capability of the remotely sensed pixel. So mapping land cover change can introduce substantial errors in data processing and interpretation. To avoid inappropriate interpretation a generic name HV (humid vegetation) was given to this land cover instead of gallery forest. In the agrosystem the rate of land cover change was high because of the agricultural practices. To solve the problem of deforestation, a reforestation project has been implemented in a participatory way with rural communities. The main goal was to create a buffer zone by restoring plant biodiversity in order to reduce the pressure on the protected area for which in spite of a low rate of deforestation experienced many modifications.
The authors are grateful to UNDESERT Program (www.undesert.neri.dk) for the financial support. UNDESERT (EU FP7 243906), "Understanding and combating desertification to mitigate its impact on ecosystem services" is funded by the European Commission, Directorate General for Research and Innovation, Environment Program. We acknowledge the Department of Bioscience of Aarhus University (Denmark) for the collaboration.