Direct Drivers Assessment#
Note
This documentation has been produced in the framework of the global project entitled: Assessment of deforestation and forest degradation and related direct drivers using SEPAL, supported by the Central African Forest Initiative (CAFI). This project, known as the Deforestation and Degradation Drivers (DDD) project is a global methodology developed and piloted in the Central African region. The Congo Basin area is used as an example but this methodology can be applied to any other geographic area.
Background#
Deforestation and forest degradation (DD) are complex processes that have many underlying causes. A comprehensive understanding of forest conversion to other land uses is instrumental for the development of policies and actions aiming to reduce the loss of forests and its associated carbon emission. It is important to understand recurring patterns and correlations that can help countries tailor their efforts towards reducing forest loss. The causes of deforestation and forest degradation vary both regionally and temporally[1]. Different studies refer to agricultural expansion (cropland and pasture) as the largest direct cause of global deforestation[2][3][4][5][6]. Commercial agriculture is estimated to be responsible for around 70-90% of the worldwide deforestation, while in Africa, both commercial and subsistence agriculture may account for similar importance in terms of forest loss, while fuel wood collection, charcoal production, and, to a lesser extent, livestock grazing in forests are the most important drivers of degradation.
However, these studies (as well as existing national studies on the drivers of deforestation and forest degradation in the Central African region) are generally based on data acquired up to 2014 and do not consider the recent trends in observed tree cover loss. They also do not take into account the importance of the spatial fragmentation of forests and the role played by degradation induced by forest exploitation [7]. Furthermore, they simplify causes of deforestation into a single driver, which does not adequately reflect the complex local context and interacting causes, decisions that lead to anthropogenic forest disturbance at local scales [8].
These recent trends and the lack of updated studies result in little consensus on the main drivers and agents of forest change at regional scales. There is a pressing need to systematically quantify the direct causes of deforestation and degradation via a systematic and comprehensive approach. An updated assessment should provide validated evidence on the role and weight of different drivers of forest loss and support decision making to address related challenges. A spatially explicit approach should also facilitate the assessment of the efficiency of policies and actions in different contexts. Improved spatial data on deforestation and forest degradation, and an improved understanding of the drivers will support land use planning approaches in two pilot areas where the regional analysis indicated clear opportunity for supporting land use decision making and planning.
In this context, FAO has developed a robust methodology to assess the recent trends and drivers of deforestation and forest degradation via collaborative approaches in which national experts, research institutes and civil society work together and join resources and data to reach a common view on the current and future causes and drivers of anthropogenic forest disturbances. These tools generate improved information for decision making and land use planning in the relevant management context.
The FAO, in partnership with the Central African Forest Initiative (CAFI), and a coalition of donor and partner countries have developed a global, standard, large-scale methodology to assess forest dynamics using cloud-computing solutions and open-source tools. The approach maps disturbances (both deforestation and degradation) and quantifies the contribution of multiple associated direct drivers. The methodology is used to assess deforestation and forest degradation trends and their associated current and historical direct drivers in six Central Africa countries as part of an international effort to mitigate climate change, support sustainable development goals, and protect biodiversity. The project builds on a collaborative approach, in which national experts, global research institutes and civil society work together, join resources and data to provide technical evidence and reach a common view on the current and historical trends, and direct drivers of forest disturbances.
definitions#
A robust methodology uses consistent definitions. The following terms and definitions are applied throughout the workflow.
Area of interest#
The pilot study area includes the national boundaries of the six Congo Basin countries: Cameroun, Central African Republic, Equatorial Guinea, Gabon, Republic of Congo and the Democratic Republic of the Congo.
Because of consistency issues between border datasets, at national / regional / global levels, it was decided to take one global dataset, Large Scale International Boundaries (LSIB), from the U.S. Department of State.
Forest definitions#
There are a total of four different national forest definitions in the six countries of the study region. These are applied at country level using information on canopy height, tree cover and pixel area.
Scale |
Source |
Date |
MMU (ha) |
Tree height (m) |
Canopy cover (%) |
Comment |
---|---|---|---|---|---|---|
Cameroon |
REDD+ National Strategy (MINEP, 2008) |
2018 |
0.5 |
3 |
10% |
Exclusion of monospecific agro-industrial plantations with a purely economic vocation and using mainly agricultural management techniques. Are always considered as forest the areas subjected to natural disturbances which are likely to recover their past status. |
Central African Republic |
FRA |
2020 |
0.5 |
5 |
10% |
|
Gabon |
Sannier et al., 2016 |
2020 |
0.5 |
5 |
30% |
Functional definition used by national monitoring system (AGEOS) |
Democratic Republic of Congo |
FREL 2018 (Ministère de l’Environnement et Développemnt Durable, 2018) |
2018 |
0.5 |
3 |
30% |
A canopy cover criterion of around 50% for an area of 0.09 ha was used during the interpretation of the samples. |
Republic of Congo |
FREL (Coordination Nationale REDD, 2017) |
2017 |
0.5 |
5 |
30% |
Exclusion of agricultural activities, in particular palm groves in production. |
Regional Land Cover#
The baseline map for the regional forest cover was first derived from a common classification system that was validated by the project technical committee and included land cover classes reference in national system. The land cover classification has also been published in the FAO Land Cover Registry.
Note
In Central African Republic and Cameroun, shrub savannas were identified as forest, in adherence to the national forest definitions of >10% tree cover
Code |
Forest/non-Forest |
English |
French |
Spanish |
Description |
---|---|---|---|---|---|
1 |
Forest |
Dense Forest |
Forêt Dense |
Bosque denso |
Dense humid primary evergreen forest on terra firme, >60% tree cover |
2 |
Forest |
Dense Dry Forest |
Forêt Dense Sèche |
Bosque denso seco |
Dense dry forest, >60% tree cover, with dry seasons |
3 |
Forest |
Secondary Forest |
Forêt Secondaire |
Bosque secundario |
Open forest, 30-60% tree cover, degraded or secondary |
4 |
Forest |
Dry Open Forest |
Forêt Claire Sèche |
Bosque claro Seco |
Dry open forest, 30-60% tree cover, with dry seasons |
5 |
Forest |
Sub-Montane Forest |
Forêt Sub-Montagnarde |
Bosque sub-montañoso |
Forest >30% tree cover, 1100-1750m altitude |
6 |
Forest |
Montane Forest |
Forêt Montagnarde |
Bosque montañoso |
Forest >30% tree cover >1750m altitude |
7 |
Forest |
Mangrove |
Mangrove |
Manglar |
Forest >30% tree cover on saline waterlogged soils |
8 |
Forest |
Swamp Forest |
Forêt Marécageuse |
Bosque pantanoso |
Swamp mixed foret, >30% tree cover, flooded > 9 months |
9 |
Forest |
Gallery Forest |
Forêt Galerie |
Bosque en galería |
Riparian forest in valleys or along river edges |
10 |
Forest |
Mature Forest Plantation |
Plantation Forestière Mature |
Plantación forestal madura |
Tree cover >15%, cultivated or managed |
11 |
Forest |
Woodland Savanna |
Savane Arborée |
Sabana arbórea |
Woodland savanna 15-30%, tree cover > national forest definition |
12 |
Forest* |
Shrubland Savanna |
Savane Arbustive |
Sabana arbustiva |
Shrubland savanna >15% shrub cover > national forest definition |
13 |
Non-Forest |
Herbaceous Savanna |
Savane Herbacée |
Sabana herbácea |
Grassland savanna <15% tree cover |
14 |
Non-Forest |
Aquatic grassland |
Prairie Aquatique |
Pradera acuática |
Regularly flooded grassland |
15 |
Non-Forest |
Bare Land |
Sols Nus - Végétation Éparse |
Suelo desnudo-Vegetación escasa |
<15% vegetation cover |
16 |
Non-Forest |
Cultivated Areas |
Terres Cultivées |
Tierras cultivadas |
Cultivated vegetation >15% vegetation cover |
17 |
Non-Forest |
Developed Areas |
Zones Bâties |
Zonas edifiadas |
Human dominated and artificial surfaces |
18 |
Non-Forest |
Water |
Eau |
Agua |
Water > 50% |
19 |
Non-Forest |
Shrubland Savanna |
Savane Arbustive |
Sabana arbustiva |
Shrubland savanna >15% tree cover < national forest definition |
Definitions of deforestation and degradation#
In order to properly discern between deforestation and degradation, we require specific and operational definitions that can be identified from satellite image analysis.
Deforestation |
Degradation |
---|---|
Permanent reduction of forest cover below the forest definition |
A temporary or permanent reduction of forest cover that remains above the forest definition |
Conversion of forest to other land use: agriculture, pasture, mineral exploitation, development, etc… |
Includes areas where timber is exploited, or trees removes and where forest may be expected to regenerate naturally or with silvicultural methods. |
Excludes areas of planned deforestation, such as timber extraction, or in areas where the forest is expected to regenerate naturally or with silvicultural methods. |
|
Includes areas where impacts, over exploitation or environmental conditions prohibit regeneration above the forest cover definition |
Example of deforestation#
Deforestation is recognizable in images by a permanent change in forest cover. In high-resolution images, we can often see bare ground, felled trees, and sometimes the beginning of agricultural or other driving activities.

Example of degradation#
Degradation is more difficult to determine because changes are more subtle, sometimes a few trees removed, and the tree cover remains above the national definition. It is therefore necessary to look at the whole time series and make sure that the changes are not deforestation. Degradation is also not the same everywhere and will differ by forest type and environmental and human context.

Date convention#
The time period for this pilot study is 2015-2022, with an assessment of changes encompassing 31/12/2015 to 31/12/2022. The year 2015 was used as the baseline, with the detection of annual changes in deforestation and degradation starting in 2016 through 2022. This fits with the availability of Sentinel satellite imagery in 2015 (although incomplete for that year), as well as new monthly high-resolution mosaics available for the tropics from Planet, which are available from 2015 and are used for additional validation.
The following date convention was adopted:
A product for the year YYYY is considered as of 31/12/YYYY.
This convention allows a consistent approach to assessing change products. A change map from year1 to year2 will be consistent with both year1 and year2 maps. The status of the year takes into account any changes that occurred during the year.
Direct Driver definitions#
A total of eight direct drivers were defined by their specific characteristics identifiable in high resolution satellite imagery from Planet.
To address the overlap of drivers in the same location and interpret local context, our approach identifies archetypes, or common driver combinations which represent realities and processes on the ground. The most common archetype consists of four drivers, which include artisanal agriculture, artisanal forestry, roads and settlements, which is representative of the agricultural mosaic, or so-called “rural complex” commonly observed in the region[9].
The observed combinations of drivers are grouped into thematic classes or archetypes.
Deforestation |
Degradation |
---|---|
Rural complex |
Artisanal agriculture with roads and settlements, with or without artisanal forestry, and no industrial drivers |
Artisanal forestry |
Artisanal forestry with or without “other” driver, or with settlements or roads without any artisanal agriculture |
Industrial Agriculture |
Industrial agriculture and other non-industrial drivers |
Industrial forestry |
Industrial forestry and other non-industrial drivers |
Industrial Forestry and Agriculture |
Industrial Forestry and Agriculture identified together |
Industrial mining |
Presence of industrial mining without other industrial drivers |
Artisanal mining |
No more than 2 drivers, including artisanal mining, no industrial drivers present |
Human infrastructure |
Roads, settlements observed alone or together, no other drivers present |
Infrastructure related agriculture |
Infrastructure and artisanal agriculture observed together |
Methodology#
The major components of this this methodology include the generation of wall-to-wall geospatial data on forest cover types, changes, and discerning areas of deforestation from degradation for the entire Central African region. Next, these products are validated via visual interpretation and the presence of various direct drivers are identified to evaluate the direct causes of disturbance, and interpreted in the context of strategic investments for climate change mitigation and support for national efforts for emissions reduction.
The methodology uses FAO’s OpenForis suite of tools including the SEPAL platform for satellite data analysis, Collect Earth Online and Google Earth Engine. The approach analyses dense satellite time-series to generate geospatial data on forest changes which are then validated and interpreted for direct drivers in 5 major steps:
Creating cloud-free mosaics: processing of optical (Landsat 4/5/7/8) and radar (Sentinel 1/ALOS PALSAR) satellite images to create mosaics for the classification of wall-to-wall maps of vegetation types, recoded to a binary forest mask (following national forest definitions), and forest fragmentation assessment for the baseline year 2015
Time-series analysis: processing of optical satellite (Landsat 4/5/7/8) time series data covering 2012-2020 (2012-2015 is the historical time period, monitoring is from 2016-2020), using seasonal models and break detection algorithms to produce a forest change map for 2015-2020 at regional scale identifying areas of both deforestation and degradation.
Sample Stratification: Stratified random sampling is conducted on the change map from step 2. Systematic validation for all points identified as change, plus a sample of stable points is conducted in Collect Earth Online, evaluating land cover types, changes and dates of change and the identification of the presence of direct drivers.
Quantify Direct Drivers: The quantification of direct drivers by forest types, fragmentation class

Creating cloud-free mosaics#
To accurately determine disturbances within forest ecosystems and distinguish from other dynamics occurring in non-forest areas, a baseline forest mask is required. This is achieved by classifying cloud-free image mosaics, which are created using the optical and radar mosaic recipes.
As you can see in this online animation, clouds are persistent in the Congo Basin region. For this reason we will produce mosaics of optical cloud-free imagery, and radar (cloud independent) composites for the best observations of the study region.
Optical cloud-free composite#
Multi-temporal image mosaics are compiled from data collected over several months or years. Cloud-free pixels from multiple images are integrated into an image with fewer clouds, haze and shadows using the pixel quality band provided with image metadata.
We evaluated the availability of Landsat 4,5,7 and 8 images for the creation of optical mosaics for the baseline year or 2015. As you can see from the figure below, only certain sensors are available for certain time periods – from 2003 onwards the Landsat 7 sensor experienced a malfunction which results in data gaps in strips. This sensor should be only included when necessary, i.e. when not enough imagery is available. Luckily in SEPAL the selection of sensors is automatic based on the selected date and only provides the available options.

The coverage of Landsat over time is shown below. The western part of the study region along the coast, results in cloudy or data gaps in Gabon, Equatorial Guinea and Cameroon.

To create our optical mosaic, we will use the SEPAL optical mosaic recipe. To learn more about the different available parameters and how to use the recipe, please see Optical mosaics.
In this example we will use a custom asset from GEE for the AOI parameter: projects/cafi_fao_congo/aoi/cafi_countries_buffer_simple
. It includes an ISO column to select Congo Basin countries according to their three digit code. See Area of Interest Selection for more AOI selection methods.
In the DAT section you can select the dates of interest.
For later years (after 2018), the sensor coverage is good so you can safely select all images of a single year.
For earlier eras, e.g. 2015 use the advanced option to add images from prior years from a targeted season (in this case the full year). This will help to fill the gaps in cloudy areas.

For data sources, more is generally better. Select all Landsat options for a consistent mosaic. If you like, Sentinel-2 can be added for more data, but as the tiling system of the 2 sensors are different you will be forced to use all available images - the option to select images will not be available.
If you are only working with Landsat (or only with Sentinel), you can manually select scenes if you want to tailor your mosaic and you have a lot of time to devote to your mosaic. USE ALL SCENES is the quicker and simpler approach and recommended for large areas.
For the composite options we recommend SR and BRDF, you can exclude pixels with low NDVI (particularly if you have a long time period) and select the following options:
You can then retrieve the mosaic as a Google asset at 30m resolution. We select the original bands as all other indices can be recalculated later: BLUE, GREEN, RED, NIR, SWIR1, SWIR2, THERMAL
Once the export is finished, you can view the asset in Google Earth Engine or SEPAL. Here is the 2015 mosaic of the Congo Basin using the above parameters:

ALOS Palsar mosaics#
Radar imagery has the added benefit of being cloud-free by design as active sensors are not influenced by clouds.
Alos Palsar is a L-band radar that gives good results for monitoring forest ecosystems. Data is provided by the Kyoto & Carbon Initiative from the Japanese Space Agency (JAXA) for the year 2015 onward. SEPAL provides an application to select, process and download them to your user space pr Google Earth Engine Account.
For more information about the parameters, Please see ALOS Kyoto & Carbon Mosaics by JAXA.
Sentinel-1 mosaics#
You can use the Sentinel-1 recipe to create a mosaic from ESA Copernicus radar data.
The aoi selection is the same as for the optical mosaic. For the dates you can enter a year, a date range, or a single date. When you add a year or date range, SEPAL will provide a “time-scan” composite which includes bands which are statistical metrics of the range of data including phase and amplitude which assess the phenology and variations within the time period.
For the best results in the Congo Basin the following parameters are proposed:
Both ascending and descending orbits will ensure complete coverage of the AOI
The terrain correction will mask any errors due to topography, or terrain “shadows”
We don’t need to apply a speckle filter
moderate outlier removal will provide the most consistent results
Select which bands to export in the retrieve window, you may select all of them depending on the space available in your GEE repository or SEPAL workspace. Resolution can also be selected accordingly - you can choose 30 to be at the same scale as the optical mosaic, which will be classified in the next step.
Time-series analysis#
Warning
This part of the documentation is still under construction.
Sample Stratification#
Warning
This part of the documentation is still under construction.
Quantify Direct Drivers#
Warning
This part of the documentation is still under construction.