Quantification of the above ground carbon stock (AGC) is important in sustainable forest management and policy advice on climate change mitigation. Traditional ground vegetation survey methods have been used to provide data for estimation of AGC stock but constrained by inadequate time and often too costly. Remote sensing when combined with few ground collected data has the potential of improving forest resource assessment even though this opportunity has not well been utilised. In this study, we mapped AGC through combination of ground survey data collected from 51 permanent sapling plots with Normalized Difference Vegetation Index (NDVI) derived from Landsat 5 Thematic Mapper image. Linkage of the two data sources was made during a training stage of supervised classification. The overall classification accuracy was 98%, suggesting that reliable estimate of AGC for a large area can be made through combination of medium resolution satellite imagery and few samples from the ground.
Forests play an important role in the global carbon cycle. They can either be a source of atmospheric carbon in the case of biomass combustion, or a sink in the case of carbon sequestration from growth. The global forest pool has been estimated to contain about 80% of the aboveground, and 40% of the belowground carbon stored in terrestrial ecosystems (Dixon et al., 1994) . About 43% of global forests are found in the tropics, of which 42% are located in arid and semi-arid areas (dry forests, woodlands and rangelands) (Brown et al., 2005) . In general, dry forests have lower biomass stocks than wetter forests. However, the more widespread coverage of dry forests make them a considerable terrestrial carbon store (Glenday, 2008) .
The estimation of carbon in forest woody biomass has importance in global climate mitigation policy and processes as well as for sustainable planning of resource conservation, utilization and management. The conventional way of ground vegetation survey involves an intensive field work data collection. It is an expensive procedure requiring more resources such as manpower and equipment (forest inventory instruments and vehicles). Field work data collection is time consuming hence exacerbating monetary expenditure for the whole process in the form of subsistence and overtime allowances (Hall et al., 2002) . Remote sensing overcomes most of challenges encountered in ground field surveys. When combined with few ground collected data, remote sensing has the potential of improving forest resource management, even though this opportunity has not well been utilised (Wulder & Franklin, 2003) . Advances in remote sensing techniques allow track changes in real time and pinpoint hotspots and specific areas of concern, providing critical knowledge that was not possible prior to this technological development. Furthermore, advantages of remote sensing become more pronounced as the region to be inventoried gets larger or more remote or both.
This study link ground data collected through ground vegetation survey and satellite imagery in order to map AGC of Kasane Forest Reserve in 2011. KFR is among the most important biodiversity in Botswana, but it’s constantly under threat of being depleted due to human activities such as township expansion, infrastructure development, fuelwood collection, and recurrent wildfires.
The study area is Kasane Forest Reserve (KFR) located within the Chobe District of Botswana and its total area is 75,040 ha (
The estimation of forest carbon in the Chobe Forest Reserves in Botswana has
in the past been neglected with only a few small scale project level forest inventories conducted. These inventories relied on using a small number of field plots mainly due to high unsustainable costs of field based inventories thus not capturing a true representation of the state of forest carbon in the forest reserve. The lack of carbon stocks information has made it difficult to monitor emissions changes therefore, presented a challenge to Botswana in meeting international reporting obligations and making sound decision for sustainable forest management.
Vegetation data was collected from 51 permanent plots by the Department of Forestry and Range Resources inventory in 2011 (
Global Positioning System (GPS) facilitated orienting direction to the next plots. Concentric plot (National Forest Resources Monitoring and Assessment of Tanzania-NAFORMA, 2010) was adopted and used in this study. To ease the counting process, each sample plot was sub-divided into two sub-plots (concentric plots) of radius 8 m and 30 m located 300 m apart. In the 8m circle, all trees of diameter at breast height (dbh 1.3 m) of 5 to 19.9 cm were measured while for 30m circle trees of 20 cm and above were measured for dbh. Tree diameters were measured using veneer caliper and/or diameter tape while tree heights were measured using Suunto hypsometer.
Vegetation data was analysed for biomass through allometric equations for Miombo woodlands with coefficient of determination (R2) of 0.95 (Malimbwi et al., 2016) . The equation is as follows:
B = 0. 1 0 27 × dbh 2 . 4798 (1)
where:
B = biomass (kg).
dbh = diameter at breast height (cm).
Above ground carbon was converted to carbon using a factor of 0.49 (National Forest Resource Monitoring and Assessment of Tanzania-NAFORMA, 2010) .
A 30 m resolution Landsat 5 Thematic Mapper (TM) satellite image of 18th April, 2011 was downloaded from USGS Earth explorer, Path 173 and row 73. The April image was acquired because it has minimum cloud cover and coincides well with the duration when ground vegetation survey was conducted. Normalized Difference Vegetation Index (NDVI) was derived from the Landsat Image using Equation (2) below. NDVI was selected because is the most commonly used metrics for predicting and estimating different biophysical parameters such as forests canopy, biomass, volume and carbon (Vafaei et al., 2018) .
NDVI = ( NIR − R ) ( NIR + R ) (2)
where:
NIR is near infra-red band value.
R is Red band value.
In this study, the NDVI map values ranged from −0.00083317 to 0.773484 (
Linkage between ground survey and satellite imagery was made through a training stage of supervised classification. Firstly, the average carbon stock and NDVI for each permanent plot was obtained. In the process, 50% of AGC was matched with NDVI and training samples obtained (
PSP ID | Biomass t/ha | Carbon t/ha | NDVI Values |
---|---|---|---|
1 | 7.05 | 3.45 | 0.26 |
4 | 45.83 | 22.46 | 0.33 |
5 | 18.97 | 9.29 | 0.22 |
6 | 12.16 | 5.96 | 0.23 |
7 | 20.51 | 10.05 | 0.31 |
8 | 33.15 | 16.25 | 0.25 |
9 | 0.53 | 0.26 | 0.20 |
11 | 2.25 | 1.10 | 0.19 |
14 | 49.07 | 24.04 | 0.36 |
17 | 46.63 | 22.85 | 0.31 |
18 | 24.71 | 12.11 | 0.30 |
19 | 18.20 | 8.92 | 0.21 |
20 | 8.26 | 4.05 | 0.23 |
22 | 36.58 | 17.92 | 0.32 |
23 | 27.21 | 13.33 | 0.23 |
24 | 21.35 | 10.46 | 0.22 |
26 | 3.24 | 1.59 | 0.21 |
27 | 26.66 | 13.07 | 0.22 |
30 | 33.24 | 16.29 | 0.33 |
32 | 9.68 | 4.74 | 0.22 |
35 | 30.88 | 15.13 | 0.23 |
36 | 14.37 | 7.04 | 0.25 |
37 | 14.03 | 6.87 | 0.20 |
Where: NDVI: Normalized Difference Vegetation Index, ABC: Above Biomass Carbon.
Supervised classification was performed using Support Vector Machine algorithms in ENVI 5.1 to produce AGC map of the entire forest (
because of shrubs and grasslands while low carbon is caused by mainly bare ground that is degraded land among other factors.
Carbon maps derived from classification of images usually contain some sort of errors due to several factors that range from classification techniques to methods of satellite data capture. Hence, evaluation of classification results is an important process in the classification procedure Yesserie (2009) . The classification accuracy was performed in ENVI 5.3 using 50% of AGC data from permanent plot (that were not used in the classification (
The findings of the present study show that field data acquired through ground survey can be successfully combined with freely available satellite imagery to estimate carbon for a larger area and provide spatially explicit maps useful for management, planning, and reporting purposes. The approach presented in the study can be applied to map carbon stock of forest with similar biophysical properties.
Class | 0 - 20 | 21 - 40 | 41 - 60 | Above-61 | Total |
---|---|---|---|---|---|
Unclassified | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
0 - 20 | 99.31 | 0.00 | 0.00 | 0.45 | 57.29 |
21 - 40 | 0.00 | 96.77 | 2.43 | 8.04 | 18.40 |
41 - 60 | 0.00 | 0.75 | 97.57 | 0.00 | 14.37 |
Above-61 | 0.69 | 2.49 | 0.00 | 91.52 | 9.93 |
Total | 100.00 | 100.00 | 100.00 | 100.00 | 100.00 |
We would like to thank Keletso Mike Seabo and Ephraim Kagiso, who provided shapefiles for the research. Appreciation also goes to Boiki Ramatho, who helped in rereading this write-up. We also thank the Department of Forestry and Range Resources of Botswana for provision of the secondary field data for this project.
The authors declare that there is no conflict of interests regarding the publication of this paper.
Basalumi, L., Kilawe, C.J. and Mauya, E.W. (2018). Linking Ground Forest Inventory and NDVI in Mapping above Ground Carbon Stock in Kasane Forest Reserve, Botswana. Open Journal of Forestry, 8, 429-438. https://doi.org/10.4236/ojf.2018.83027