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.

Driver

example

characteristics

Artisanal agriculture

Small-scale agriculture is composed of small, informal, unstructured and irregular agricultural plots covering an area of less than 5ha. The presence of fires (slash-and-burn agriculture) can be observed, and the land is often soil cover in various stages of cultivation.

Industrial agriculture

Industrial agriculture is characterized by agricultural areas larger than 5 ha that tend to be homogeneous and often consist of a single crop. In some cases, agriculture may be more varied and consist of many fields closely packed together. Therefore, large areas consisting of many small fields cultivated at the same time are also considered industrial agriculture under the definition.

Infrastructure

Roads are visible in the images with linear features and are identified as motorized when they are wide enough (5m) to carry vehicle traffic. Small irregular paths through vegetation are not included. Roads can be large highways, or logging trails, and are most often found with other engines such as villages, mining facilities.

Settlements

Villages and settlements can be hard or soft roofed, they can be buildings or huts, and they are often accompanied by roads and other drivers such as small-scale agriculture. This engine can be an urban area (left image), or a small isolated village in a forest stand (right image).

Artisanal forestry

Small-scale or artisanal logging is characterized by the selective extraction of trees in an irregular manner, leaving a tree cover. These are areas that are not visibly cultivated. These areas are often found in places accessible by small roads or villages.

Industrial forestry

Large-scale or industrial forestry is recognizable by the presence of logging roads, along which selective logging degradation occurs. These roads may be recent or old, and the canopy can quickly cover them, so all years of imagery acquired over the entire study period are evaluated.

Artisanal mine

Small-scale mining is characterized by muddy clearings, and usually ponds or water catchments and may feature turbid water. Artisanal in nature, the clearings are generally small, isolated, and often located along streams.

Industrial mine

Large-scale mining is characterized by large ponds, open pits and clearings, as well as extensive infrastructure and roads.

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:

  1. 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

  2. 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.

  3. 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.

  4. 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.