Incorporating supervised feature selection procedures using Landsat 8 imageries. In this regards, genetic algorithm is incorporated to choose best indices by employing kernel base one, support vector regression and linear regression methods. The proposed method revealed that there is a high degree of consistency between affected information.
Published online 2016 Jun 24. doi: 10.1186/s13021-016-0055-8
PMID: 27418944
This article has been cited by other articles in PMC.
AbstractBackground
A functional forest carbon measuring, reporting and verification (MRV) system to support climate change mitigation policies, such as REDD+, requires estimates of forest biomass carbon, as an input to estimate emissions. A combination of field inventory and remote sensing is expected to provide those data. By linking Landsat 8 and forest inventory data, we (1) developed linear mixed effects models for total living biomass (TLB) estimation as a function of spectral variables, (2) developed a 30 m resolution map of the total living carbon (TLC), and (3) estimated the total TLB stock of the study area. Inventory data consisted of tree measurements from 500 plots in 63 clusters in a 15,700 km2 study area, in miombo woodlands of Tanzania. The Landsat 8 data comprised two climate data record images covering the inventory area.
Results
We found a linear relationship between TLB and Landsat 8 derived spectral variables, and there was no clear evidence of spectral data saturation at higher biomass values. The root-mean-square error of the values predicted by the linear model linking the TLB and the normalized difference vegetation index (NDVI) is equal to 44 t/ha (49 % of the mean value). The estimated TLB for the study area was 140 Mt, with a mean TLB density of 81 t/ha, and a 95 % confidence interval of 74â88 t/ha. We mapped the distribution of TLC of the study area using the TLB model, where TLC was estimated at 47 % of TLB.
Conclusion
The low biomass in the miombo woodlands, and the absence of a spectral data saturation problem suggested that Landsat 8 derived NDVI is suitable auxiliary information for carbon monitoring in the context of REDD+, for low-biomass, open-canopy woodlands.
Keywords: Biomass, Carbon, Modeling, Miombo woodlands, REDD+, NDVI
Background
Deforestation and forest degradation in the tropics continue to be a significant source of global carbon dioxide (CO2) emissions and the largest source of CO2 in the developing world [1]. After years of negotiations, the United Nations Framework Convention for Climate Change (UNFCCC) gradually recognized a mitigation mechanism for reducing emissions from deforestation and forest degradation and forest conservation, sustainable management of forests and enhancing forest carbon stock in the tropics and subtropics (REDD+) [2]. REDD+ deploys results-based finance to incentivize emissions reduction, based on a functional forest carbon measuring, reporting and verification (MRV) system [3]. Nevertheless, technical challenges in MRV [4, 5] have, thus far, substantially contributed to the lack of progress for implementation of REDD+ on the ground.
A functional MRV to support REDD+ requires estimates of the area of forest loss and gain and the corresponding carbon stock and changes [3, 6]. These data are needed for the estimation of the actual emissions and the construction of forest reference emissions level (FREL), a benchmark against which the actual emissions are compared [7]. A combination of field inventory and remote sensing is expected to provide those data. National forest inventories (NFIs) with repeated measurements are often lacking for most of the forests and woodlands of Africa. On the other hand, there exists a wide range of remote sensing technologies, including sensors on-board aircraft and space-based platforms. Landsat is one of the most attractive remote sensors because its images are freely available and have medium resolution, large spatial and frequent temporal coverage [].
Landsat data have been widely used in forest aboveground biomass (AGB) estimation, commonly through developing empirical relationships between AGB or other forest characteristics and spectral indices such as the normalized difference vegetation index (NDVI) derived from satellite data [e.g., â17]. Large scale and coarse resolution biomass or carbon mapping using satellite imagery for tropical forests in Africa has been demonstrated [10, â20]. Fine to medium scale carbon maps offer more detailed and locally accurate auxiliary information in the context of REDD+. For instance, the majority of the REDD+ countries that have submitted a FREL proposal to the UNFCCC [21] used medium resolution maps based, largely, on Landsat to extract activity data, i.e., forest area change. Biome or forest type specific biomass models or carbon maps with higher spatial resolution can also provide estimates of carbon densities [22], that can be used as emissions factors that are more reliable than for instance the IPCC default values [23].
For Landsat based approaches, cloud and cloud shadows [24] as well as canopy reflectance rules='groups'>Variable (t/ha)MinimumMeanStd. Dev.MaximumTLB088.4964.52436.20AGB062.1048.73350.30BGB026.3916.8889.97TLC041.5930.33205.00
Values are based on plot measurements (n = 500). The minimum (0) are from plots without trees
Landsat 8 data and spectral variables for modeling biomass
Two Landsat 8 surface reflectance climate data record (CDR) images, with 30 m resolution, (Worldwide Reference System 2: path 166 and row 66 and 67) were acquired through the USGS Earth Explorer (earthexplorer.usgs.gov). Both images were from 31 July 2014, and were almost cloud free over the study area. The data acquisition date corresponded to the cool early dry season [44], which is expected to offer a good contrast between understory vegetation and the woody biomass for open woodland, as the understory vegetation begins to dry off. The two images were pre-processed by the provider, i.e., corrected for atmospheric conditions and topography, and with cloud and cloud shadow masking. The CDR data were derived from Landsat 8 Operational Land Imager (OLI) using Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) software [45], which applies moderate resolution imaging spectroradiometer (MODIS) atmospheric correction routines to Level-1 scenes [46].
For each NAFORMA plot, we extracted seven spectral reflectance values (bands 1â7) from the pixel containing the plot center, using the default extraction procedure in ArcGIS The surface reflectance specification and the center wavelengths of each Landsat 8 band are available in [46]. We also tested a second approach which extracts the weighted mean value of the four nearest pixels, weighting the contribution of each pixel by its fraction of the 707 m2 plot area. The second approach was thought to reduce sampling error in the case of spatial mis-registration of the image. However, a paired t test for the differences between spectral indices developed based on the values extracted using the two extraction methods did not show any significant difference (p > 0.05). Furthermore, the spectral variables from the two methods produced similar correlation coefficients with biomass, with no gain in using the second approach which was finally discarded.
Seven spectral indices (Table 2) were calculated using the different spectral band combinations. These particular indices were selected not only because they are related to land surface features and vegetation biomass, but also because they have some advantages for practical applications since they are readily provided from the USGS on demand [46].
Table 2
Spectral indices derived from the spectral band reflectance values
NIR near infra-red (band 5), SWIR1 short wave infra-red (band 6), SWIR2 short wave infra-red 2 (band 7), sqrt square root
Model development
In order to examine the relationship between TLB and the 14 spectral variables (seven spectral bands and seven spectral indices in Table 2), we calculated Pearsonâs correlation coefficient (Ï), and also visually analyzed the scatter plots of TLB against those spectral variables. All the 14 spectral variables correlated significantly with TLB, but they were also highly inter-correlated (p < 0.005). The latter would lead to a multicollinearity problem, if those variables were used simultaneously as predictor variables in a regression model. We, therefore, evaluated the strength of the relationship between TLB and each spectral variable using correlation and preliminary regression analysis to select the best explanatory variable for model development.
As a means to understand the form of the relationship between TLB and the spectral variables, we aggregated the plots in the respective clusters and graphically evaluated the scatter plots of TLB against the spectral values. While the aggregated cluster level data were used for graphical analysis, data from the individual plots (n = 500) were used for the final model development.
The scatter plots of the TLB against the spectral values suggested linear models may sufficiently describe the relationships. Due to the hierarchical structure of the data (plots nested in clusters) a linear mixed modeling approach was appropriate [cf., 54, 55]. We evaluated alternative mixed models with the random intercept, random slopes or both, on the basis of improvements in Akaikeâs Information Criterion (AIC) and the Bayesian Information Criterion (BIC). The mixed models with random intercept and slope (that varies with cluster) (Eq. 4) produced the smallest AIC and BIC.
4
where Y denotes the dependent variables (TLB), X denotes the independent spectral variables and i indexes plots within clusters (1, 2,â¦, 8), and j indexes clusters (1, 2,â¦,63). Coefficients α and β represent the intercept and the slope of the fixed part of the model, respectively, and aj and bj are the random intercept and slope, respectively (deviations around the value α and β). The aj and bj are assumed to be distributed normally with mean zero, and G-varianceâcovariance matrix given as:
In the G varianceâcovariance matrix, the Ïint2 and Ïslope2 represent variance components for the random intercept and slope, respectively; the (Ïint,slope2) are covariance components representing the correlation between the random intercept and the random slope.
The eij are the residuals of ith plot in the jth cluster and are assumed to be normally distributed with mean 0 and variance Ï2, i.e., eij ~ N (0, Ï2). We fitted the model using the Restricted Maximum Likelihood (REML) with unstructured varianceâcovariance, using the mixed procedure in SAS [55]. We applied the empirical option for the mixed effects procedure in SAS, also known as the âsandwichâ estimator which computes a robust varianceâcovariance matrix of the fixed effects parameters using an asymptotically consistent estimator [56â59].
We evaluated model performance based on root mean square error (RMSE) and RMSE as a percent of the mean of the field inventory TLB (%RMSE), and absolute bias (Eqs. 5, 6, 7) and residual diagnostics. To evaluate the accuracy of the final TLB model, we used a leave-one-out cross validation (LOOCV) in which we left out a cluster at a time.
5
6
7
where n is number of plots, y^i is predicted TLB for the ith plot using model (4) fitted without the plots in the same cluster as the ith plot, and (yi) is the observed TLB for the ith plot.
Mapping and estimation of total living carbon
We applied the best TLB model to the NDVI composite of the two Landsat 8 CDR images and estimated the mean TLB density and the TLB stock of the entire study area. To estimate a 95 % confidence interval of the mean TLB density and total biomass stock, we used a model-based variance estimator for the population (the study area) as described in McRoberts [60], i.e.,
8
where, Var^Y¯^ is the variance estimator of the population mean estimate Y¯^. Y¯^ was estimated as the mean of the pixel-wise model predictions across the entire study area, X¯ is the vector for the population mean of the explanatory variables including the intercept, T stands for transpose of a matrix, and Σ^ is the varianceâcovariance matrix of the fixed parameter estimates (intercept and slope of the model, Eq. 4). We compared the TLB density estimate of the study area with a value estimated from three recent pan-tropical biomass maps of Baccini et al. [20], Saatchi et al. [] and Avitabile et al. []. The later was developed by combining the data sets of the earlier two [, 20]. Since these biomass maps are for AGB, the mean ratio of BGB to AGB of 0.43 from our field data (n = 500) was used to calculate the corresponding TLB.
We applied the TLB model to the two NDVI composite images covering the study area to construct a 30 m resolution TLC map, where TLC is estimated at 47 % of TLB. To further demonstrate the utility of the model for local carbon mapping and stock estimation, we extracted the carbon map of 13 wards, using the ward boundaries as a mask. We compared the estimated TLC values with values based on field data from previous works in the region. In the absence of other maps for that particular year and season, we visually compared extracts of the carbon distribution pattern from our map with image extracts from Google Earth [62] of the same year and season. The 13 wards have a combined area of 4780 km2 and included Liwale town, the major town in Liwale district. We selected this area, particularly because it consisted of settlements as well as a number of human activities, such as land clearing for agriculture, fuel wood and charcoal, as well as fire, which are likely to influence the forest cover and thus carbon stock and its distribution.
ResultsRelationships and regression models for TLB
All the 14 spectral variables were significantly correlated to TLB and inter-correlated amongst each other (p < 0.005). The Pearson correlation (Ï) between the spectral variables and TLB ranged from |0.18| to |0.50| (Table 3). The indices NDVI, EVI, SAVI, and MSAVI and the spectral bands-Blue and Red correlated most strongly with TLB (Ï â¥ |0.45|). It appeared, therefore, that any one of these variables, particularly the EVI, would be a useful predictor in the model. Nevertheless, the correlation (Ï) of EVI, although slightly better, was not convincingly larger (Table 3) to forgo NDVI, which is popular, intuitive and easier to interpret and to compare with other works. NDVI was thus selected as a single predictor for the final model. The large variability in the preliminary scatter plots of TLB against spectral data of individual plots limited our ability to detect absence or presence of a data saturation effect. It is likely that the small size of the field plots contributed to the observed large variability. However, the scatter plots of the NDVI against TLB, particularly those based on mean values at cluster level data showed no sign of saturation at higher biomass (Fig. 2) and suggested that the observed relationships were linear.
Table 3
Pearson correlation coefficients between spectral variables and TLB (t/ha) of plots (n = 500)
Prob > |Ï| under Ho: Ï = 0, and Ï is significant (p < 0.005) for all the spectral variables. The level of confidence α = 0.05
Scatter plots of the total living biomass (TLB) (t/ha) against the NDVI [plot values (left), n = 500, and cluster mean values (right), n = 63]
Parameter estimates and the associated statistics of the linear mixed model are presented in Table 4, along with the fit statistics. Plots of the residuals against the fitted TLB values, except for few data points of NDVI values less than 0.2, showed that the residuals had a constant variance over the range of TLB. The estimated covariance matrix used for the variance estimators of the mean and the total TLB for the entire study is given in Table 5.
Table 4
Parameter estimates of the fixed part of the TLB model (p < 0.001) with NDVI as a predictor
Statistics presented for the model construction and cross validation (LOOCV)
Table 5
Empirical covariance matrix for the fixed effect parameters (Σ^) where values in the diagonal are the variances of the intercept and the slope, respectively, and the off diagonal elements are the covariance between the two
Biomass and carbon mapping and estimation
The estimated mean TLB density for the study area was 81 t/ha with a 95 % confidence interval of ±7 t/ha. The total TLB stock was estimated at 140 Mt, with a 95 % confidence interval of ±14 Mt. We used the TLC map, developed based on the TLB model, to estimate carbon stock by ward and visualize the carbon distribution pattern among 13 wards for the year 2014 (Fig. 3). The mean TLC varied considerably among the 13 wards, ranging from 24 to 51 t/ha. For example, the four most eastern wards in the map had carbon densities of 25â29 t/ha, which was less than the average of the entire study area. In contrast, the four most western wards had carbon densities of 42â51 t/ha, which was higher than the average for the study area. The spatial distribution of predicted TLC corresponded well to the extracts from Google Earth images (Fig. 4). As can be expected, settlement areas and farmlands had low predicted TLC, and the dark green areas representing forests depict higher values.
Total living carbon (TLC) density map of 13 wards in Liwale district. Numbers in each ward are mean TLC (t/ha). The reference, Liwale town falls within the boundaries of four wards, Nangando, Likongowele, Liwale Mjini, and Liwale B
Visual comparison between the extracts of the 30 m resolution total living carbon (TLC) map in Liwale district left panels (based on this study)âand extracts of Google Earth imagesâright panels, showing settlement areas (Liwale town)âtop two panels, clear cut areas in Liwale Mjini wardâmiddle panels, and a dense forest in Makata wardâbottom panels. The TLC density map and google earth image extracts are from the year 2014. In the extracts of the map (left panels), color ranges are from low (deep red) to high (green) TLC values
DiscussionRelationship between biomass and NDVI
The linear relationship between the TLB and NDVI suggested that spectral data saturation is not a problem when using Landsat for biomass and carbon monitoring in miombo woodlands. Data saturation related to optical sensors is reported as a major limitation or source of uncertainty in forest biomass estimation in complex stand structures and closed canopy forests [10, 22, 24, 63]. A similar study in miombo woodlands in Mozambique [64] found a significant linear correlation (Ï = 0.30) between NDVI and the AGB. In dry land open canopy vegetation in Argentina as well, NDVI data saturation was, not observed [65]. Pixel-level NDVI is influenced by the degree of absorption of the red wavelengths by chlorophyll, which is proportional to the amount of chlorophyll in the leaf and the reflectance of near infrared (NIR) radiation, which is proportional to green leaf density [47]. The relatively low aboveground biomass in the open canopy woodlands might have resulted in the low reflectance values in the NIR band, which overcomes the saturation problem of the NDVI. The NIR, although statistically significant, produced the weakest correlation with TLB (Ï = 0.18), suggesting less green leaf density and thus low reflectance in open canopy woodlands. For boreal forests with evergreen conifers and relatively closed canopy forests in Canada [66, 67] and mountain birch forests in Finland [68], NIR was considered a strong predictor of stand attributes such as height, crown closure or AGB.
Geographic factors, such as elevation, slope and aspect as well as shadows are also known to affect spectral reflectance values, and thus compromise biomass estimation performances of models based on spectral variables [24]. However, NDVI, as a ratio of two bands, is expected to correct potential errors due to topography and shade and compensate the variation in illumination resulting from terrain features [28]. Lyon et al. [27] compared seven vegetation indices to detect land-cover change in a Chiapas, Mexico and reported that the NDVI was least affected by local geographic factors such as slope and aspect.
The absence of data saturation facilitated development of a linear model between TLB and NDVI. Model fit statistics, i.e., RMSE (%) and less bias in the cross validation (Table 4), along with the constant residual errors (not shown) suggested the model has a comparable accuracy with previous works based on similar approach. These results are comparable to many other models for vegetation properties based on Landsat data, where model %RMSE values were in the range of 21â50 % [13, 28, 67, 69]. Other studies that have reported higher model accuracies [e.g., 68, 70] were based on data from single species stands or plantations. Furthermore, many of these studies attempt to model AGB, rather than TLB, mostly because of the relative ease of estimating AGB or availability of allometric models for AGB only. The model presented here enables direct estimation of TLB, including the BGB. In the absence of BGB data, the default ratio of 0.28 between the BGB and AGB [36] is recommended for tropical dry forests for estimating TLB. This default ratio is much lower than the mean BGB to AGB ratio of 0.43 in our data. The BGB, although significantly related to the AGB, in this study was independently calculated for each plot based or the below ground tree allometric equations. The TLB model can, therefore, be considered an important contribution for miombo woodlands, where data availability for the BGB remains a challenge.
The performance of the model showed that NDVI is a good predictor of TLB in open woodlands. Yet, there remains a sizeable amount of unexplained variation in the data, which can be attributed to a number of sources. Miombo woodlands are subject to distinct seasonal processes, and the beginning of the dry season marks a change in phenological cycles and leaf density might decrease, compared to the rainy season when ecosystem productivity is higher. NDVI, as a greenness indicator, cannot clearly distinguish variations in biomass related to phenology, the complexity of forest species and canopy structure. Shrub and grass layers were not measured in the field for the BGB and AGB calculations, but their contribution to the total biomass may be small. Nevertheless, in the spectral variable extraction, the effects of the shrub and grass layer on canopy reflectance are unknown, but could be large, depending on the intensity of greenness or the grass layer phenology during image acquisition. Other sources of uncertainty include a potential mismatch between the field coordinates and the points of spectral variable extraction. The field plots are circular with a size of 707 m2 while the spectral variables are extracted from square pixels of size 900 m2. Furthermore, tree level BGB and AGB were derived from allometric models with sizeable errors that are likely to propagate through to plots, pixels and thus to the TLB model [e.g. 71]. Therefore, model sensitivity to those factors, including phenology, stand structure, and their application in areas of low biomass need to be investigated. Incorporating auxiliary biophysical predictors, such as land-use, land-cover or vegetation structure information may improve model accuracy, by reducing the spectral variability within those biophysical predictors [e.g. see 10, 67, 72].
Other remote sensing data sources with higher resolution may produce more accurate models and estimates of better precision. In a recent study, conducted in parts of our study area, Næsset et al. [73] estimated the relative efficiency (RE), i.e., variance of the field-based estimate relative to the variance of that of the different remotely sensed data to quantify contribution of different remote sensing data sources to improve precision of AGB estimates. Their results showed that RE values were 3.6 for ALS, 3.3 for RapidEye, and 2.8 for InSAR; while for Landsat and PALSAR the RE values were only 1.3â1.4. The latter two contributed only marginally to improve precision. The limited spatial and temporal coverage as well as the costs of acquiring and processing those types of data, however, limit their extensive applicability to larger areas. In contrast, the use of Landsat 8 data for modeling purposes can be justified as the data cover large areas and are freely available. Furthermore, Landsat 8 came with considerable improvements in sensor signal-to-noise performance and associated improvements in radiometric resolution [35]. Dube and Mutanga [70] reported an improved accuracy for estimation of AGB, for instance for E. dunii plantations, by using Landsat 8 (RMSE = 26.54 % of the mean) compared to Landsat 7 ETM+ (RMSE = 35.30 % of the mean). Furthermore, the availability of higher order data (Landsat 8 CDR) products, for free, reduces the efforts and costs of data preprocessing, which otherwise require complex procedures.
Carbon mapping and estimation
The carbon mapping and estimation was based on the TLB model and estimation. The mean TLB density of 81 t/ha for the study area is higher than the estimate based on Avitabile et al. [], 48 t/ha, and lower than that of Baccini et al. [20], 138 t/ha; but closer to that of Saatchi et al. [], 73 t/ha. Since these tropical biomass maps and our maps were based on different data acquisition systems and methods, calibrated with different field data sets and have different spatial accuracies, the large differences among them in the biomass estimates for smaller localities are not unexpected. Previous studies showed that maps of Baccini [20] and Saatchi [] often show very strong local differences, but at national and global scales their estimates tend to converge [, ], mostly due to compensation of contrasting estimates when averaging over large areas []. Moreover, data for these biomass maps were acquired on average around 2007 [], while that of our study are from 2014. The dynamics in the miombo woodlands such as clearing and fire as well as reforestation means that, maps based on data from the same year or season are needed for a more confident comparison.
The developed medium spatial resolution (30 m) map for prediction of TLC, along with the model-based variance estimates for the predicted values of the TLC, demonstrated the potential utility of Landsat 8 for forest carbon monitoring. One of the major applications of carbon maps includes use of the map as auxiliary data in estimation of carbon density and stocks, and the corresponding land area. Forest carbon accounting principles in the context of REDD+ draws on the basic formula where the emissions from REDD+ activities are calculated as estimates of areas (activity data) multiplied by estimates of change in carbon density in those areas (emissions factors) [36]. Both activity data and emissions factors are required to be as accurate as practicable [36], since the total uncertainty of emissions is calculated by combining the uncertainties of the two.
The lack of ecosystem specific knowledge on carbon stock (the basis for emissions factors), and methods to determine the activity data (areas) are considered the major limitations in the progress for REDD+, particularly in Africa [4, 5]. The approach presented here can help determining local carbon density and stock for a defined area, in combination with available land use or land cover maps as a mask. But for a given reporting region or sub-region such as REDD+ project areas, carbon maps based on Landsat 8 calibrated with field observations can be used to estimate carbon density and total carbon stock within the boundary, so as to maximize potential financial benefits from payments or credits for carbon.
The estimated mean carbon density of 38 t/ha for the entire study area is very close to the IPCC Tier 1 default value of 36 t/ha for tropical dry forests in Africa [23]. But the wide range of TLC (25â51 t/ha) among the 13 wards indicates a large spatial variation and suggests the utility of spatial maps to visualize the carbon distributions and estimate carbon density or stock for smaller administrative units or localities where carbon reporting might be required. The estimated mean and range for the 13 wards are close to the living carbon for miombo woodlands in Mozambique, (average 29.8 t/ha ± 13.07) [] and 28â36 t/ha [76]. Tanzanian studies report only the AGB estimates. In an attempt to compare with our results, we added the BGB to their AGB estimates, using the BGB to AGB ratio of 0.43 in our data. Our estimates of TLC density range correspond with the 39â54 t/ha estimated by Kashindye et al. [77] and 15â47 t/ha estimated by Shirima et al. [78]. While the two studies are from different areas of the miombo woodlands in Tanzania, another study from three villages located within the boundary of our study area presented a range of 24â28 t/ha [79], which is at the lower range of our estimate. It is, however, important to note that the other studies considered only plots with trees while our estimates came from all land use types, including agricultural lands.
The other potential application of the carbon map includes detection of spatial patterns of carbon density in connection with human induced disturbances such as settlements, fire or forest clearing. As expected, lower carbon densities are exhibited near and inside settlements such as Liwale town (Fig. 4, top panels) and in wards located east of Liwale town (Fig. 3). Furthermore, Google Earth images from June 2014 and our carbon map (Fig. 4, middle panels) show areas of low carbon corresponding to patches of clear cut areas, perhaps cleared for agriculture or wood extraction and roads. Clearing for agriculture was estimated to result in an annual forest loss of 10â25 ha in Liwale district [80], which would amount to 380â950 tons of carbon, given a carbon density estimate of 38 t/ha in this study. On the other hand, higher biomass vegetation can be shown for instance, in dense forests (Fig. 4, the bottom two panels). Fire is also an important component of the miombo woodland ecosystem, affecting up to 30 % of the woodlands in Liwale district [80]. The low carbon densities in eastern wards such as Kiangara and Mirui (Fig. 3) correspond with the presence of large forest fires in these areas during the year 2013 [62]. Furthermore, there was an indication of fire scar detected using visual analysis of the fire-sensitive short-wave infrared (band 7) complemented by the Near Infra-Red (band 5) and the visible (band 2) [cf., 81] of the same Landsat 8 data used for this study.
Our approach may also offer a potential for detection of carbon change due to deforestation, fire or land use change. This can be accomplished by image-differencing or calculating the differences in TLB values between two points in time. Nevertheless, several factors influence the reflectance values from different satellite images, including vegetation phenology and stand structure, atmospheric conditions when images were taken, and differences in sensors and image pre-processing procedures. This requires recalibration of the model coefficients to apply to other NDVI composite images of different areas and dates of acquisition. These factors are not addressed here, but understanding and identifying these sources of uncertainty will help to refine the model and improve its applicability to carbon mapping and estimation and applications to other seasons or regions.
The regression model slightly under predicts carbon for pixels with low NDVI, but will still be reliable for miombo woodlands, because biomes typically comprised of open woodlands and savannah exhibit NDVI values greater than 0.4 [82], and 95 % of the field plots in our study also exhibited and NDVI larger than 0.4. Although plots without trees and thus lower NDVI are few in the data, we had a true probability sample and the sampling frame covered other land use types including agriculture, and burnt areas which are common in the miombo woodlands.
Conclusion
Landsat 8 CDR provides suitable data for monitoring of forest biomass and carbon in miombo woodlands. The open canopy and low biomass of the miombo woodlands means that there was low or no data saturation problem, which otherwise is a well-known challenge in using Landsat sensors. This property facilitated the development of a simple linear model which provides the basis for mapping forest carbon, estimating carbon stock and detecting its spatial distribution. The developed TLB model and carbon map also includes the below ground biomass, which otherwise are often estimated using a default root to shoot ratio from AGB. The Approaches presented here i.e., modeling, carbon mapping, and estimation can assist to estimate carbon and detect areas of low or high forest carbon, which are relevant to REDD+ activities, particularly, deforestation and enhancing carbon stock.
Authorsâ contributions
All the authors have made substantial contribution towards the successful completion of this manuscript. BG designed the study, carried out the Landsat data processing, the statistical analyses and modeling, and did the majority of the writing. SS and EN have contributed to the design of the study, statistical analyses, writing and reviewing the manuscript. OMB contributed to the population inference and estimation and commented the manuscript. EN, TG, OMB, EZ and EWM contributed to the field survey design and inventory data collection, and commented the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The Norwegian Institute of Bioeconomy Research (NIBIO) supported data analysis and writing up for the lead author (Project Number: 341110-25). The Royal Norwegian Embassy in Tanzania is acknowledged for funding the field work of this study through the project âEnhancing the measuring, reporting and verification (MRV) of forests in Tanzania through the application of advanced remote sensing techniquesâ. The project is also a contribution to the Government of Norwayâs International Climate and Forest Initiative. Comments and suggestions of two anonymous reviewers have improved the manuscript.
Competing interests
The authors declare that they have no competing interests.
Abbreviations
Contributor Information
Belachew Gizachew, Email: on.oibin@geb, Email: moc.oohay@zwehcazig.
Svein Solberg, Email: [email protected].
Erik Næsset, Email: [email protected].
Terje Gobakken, Email: [email protected].
Ole Martin Bollandsås, Email: [email protected].
Johannes Breidenbach, Email: [email protected].
Eliakimu Zahabu, Email: moc.oohay@eubahaz.
Ernest William Mauya, Email: moc.oohay@ayuamtsenre.
References
1. IPCC. Climate change 2014: mitigation of climate change: IPCC working group III contribution to the Fifth Assessment Report of the Intergovernmental Panel on climate change. IPCC WGIII; 2014.
2. UNFCCC. Decision 2/CP. 13: reducing emissions from deforestation in developing countries: approaches to stimulate action. United Nations Framework Convention on Climate Chgange, Bonn. 2007. http://www.unfccc.int/resource/docs/2007/cop13/eng/06a01.pdf. Accessed 19 Sept 2014.
3. UNFCCC. Decision 4/CP.15. Methodological guidance for activities relating to reducing emissions from deforestation and forest degradation and the role of conservation, sustainable management of forests and enhancement of forest carbon stocks in developing countries. United Nations Framework Convention on Climate Change, Bonn. 2009. http://www.unfccc.int/resource/docs/2009/cop15/eng/11a01.pdf.
4. Herold M. An assessment of national forest monitoring capabilities in tropical non-Annex I countries: recommendations for capacity building. GOFC-GOLD Land Cover Project Office. Final Report; 2009.
5. Romijn E, et al. Assessing capacities of non-Annex I countries for national forest monitoring in the context of REDD+ Environ Sci Policy. 2012;19â20:33â48. doi: 10.1016/j.envsci.2012.01.005. [CrossRef] [Google Scholar]
6. UNFCCC. Report of the conference of the parties on its nineteenth session, held in Warsaw from 11 to 23 November 2013 Bonn; 2014.
7. UNFCCC. United Nations Framework Convention on Climate Change (UNFCCC) 2011 Decision 12/CP.17 on guidance on systems for providing information on how safeguards are addressed and respected and modalities relating to forest reference emission levels and forest reference levels as referred to in decision 1/CP.16: appendix I COP 17 decisions. 2011. http://www.unfccc.int/files/meetings/durban_nov_2011/decisions/application/pdf/cop17_safeguards.pdf.
8. Hansen MC, et al. High-resolution Global maps of 21st-century forest cover change. Science. 2013;342(6160):850â853. doi: 10.1126/science.1244693. [PubMed] [CrossRef] [Google Scholar]
9. Anderson GL, Hanson JD, Haas RH. Evaluating Landsat thematic mapper derived vegetation indexes for estimating aboveground biomass on semiarid rangelands. Remote Sens Environ. 1993;45(2):165â175. doi: 10.1016/0034-4257(93)90040-5. [CrossRef] [Google Scholar]
10. Avitabile V, et al. Capabilities and limitations of Landsat and land cover data for aboveground woody biomass estimation of Uganda. Remote Sens Environ. 2012;117:366â380. doi: 10.1016/j.rse.2011.10.012. [CrossRef] [Google Scholar]
11. Gomez C, et al. Historical forest biomass dynamics modelled with Landsat spectral trajectories. ISPRS J Photogramm Remote Sens. 2014;93:14â28. doi: 10.1016/j.isprsjprs.2014.03.008. [CrossRef] [Google Scholar]
12. Helmer EH, Lefsky MA, Roberts DA. Biomass accumulation rates of Amazonian secondary forest and biomass of old-growth forests from Landsat time series and the Geoscience Laser Altimeter System. J Appl Remote Sens. 2009;3:033505. doi: 10.1117/1.3082116. [CrossRef] [Google Scholar]
13. Ji L, et al. Estimating aboveground biomass in interior Alaska with Landsat data and field measurements. Int J Appl Earth Obs Geoinf. 2012;18:451â461. doi: 10.1016/j.jag.2012.03.019. [CrossRef] [Google Scholar]
14. Kelsey KC, Neff JC. Estimates of aboveground biomass from texture analysis of Landsat imagery. Remote Sens. 2014;6(7):6407â6422. doi: 10.3390/rs6076407. [CrossRef] [Google Scholar]
15. Propastin P. Large-scale mapping of aboveground biomass of tropical rainforest in Sulawesi, Indonesia, using Landsat ETM plus and MODIS data. Gisci Remote Sens. 2013;50(6):633â651.[Google Scholar]
16. Zheng DL, et al. Estimating aboveground biomass using Landsat 7 ETM+ data across a managed landscape in northern Wisconsin, USA. Remote Sens Environ. 2004;93(3):402â411. doi: 10.1016/j.rse.2004.08.008. [CrossRef] [Google Scholar]
17. Zhu XL, Liu DS. Improving forest aboveground biomass estimation using seasonal Landsat NDVI time-series. ISPRS J Photogramm Remote Sens. 2015;102:222â231. doi: 10.1016/j.isprsjprs.2014.08.014. [CrossRef] [Google Scholar]
18. Saatchi SS, et al. Benchmark map of forest carbon stocks in tropical regions across three continents. Proc Natl Acad Sci USA. 2011;108(24):9899â9904. doi: 10.1073/pnas.1019576108.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
19. Baccini A, et al. A first map of tropical Africaâs above-ground biomass derived from satellite imagery. Environ Res Lett. 2008;3(4):045011. doi: 10.1088/1748-9326/3/4/045011. [CrossRef] [Google Scholar]
20. Baccini A, et al. Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nat Clim Chang. 2012;2(3):182â185. doi: 10.1038/nclimate1354. [CrossRef] [Google Scholar]
21. UN-REDD. Emerging approaches to forest reference emission levels and/or forest reference levels for REDD+. 2014. http://www.un-redd.org/FRELPublication/tabid/794487/Default.aspx. Accessed 13 Nov 2015.
22. Lu D, et al. Aboveground forest biomass estimation with Landsat and LiDAR data and uncertainty analysis of the estimates. Int J For Res. 2012; 2012.
23. Aaron R, Gibbs HK. New IPCC Tier-1 Global biomass carbon map for the year 2000. Carbon dioxide information analysis center. Oak Ridge: Oak Ridge National Laboratory. 2008. http://www.cdiac.ornl.gov.
24. Lu D. Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon. Int J Remote Sens. 2005;26(12):2509â2525. doi: 10.1080/01431160500142145. [CrossRef] [Google Scholar]
25. Powell SL, et al. Quantification of live aboveground forest biomass dynamics with Landsat time-series and field inventory data: a comparison of empirical modeling approaches. Remote Sens Environ. 2010;114(5):1053â1068. doi: 10.1016/j.rse.2009.12.018. [CrossRef] [Google Scholar]
26. Steininger MK. Satellite estimation of tropical secondary forest above-ground biomass: data from Brazil and Bolivia. Int J Remote Sens. 2000;21(6â7):1139â1157. doi: 10.1080/014311600210119. [CrossRef] [Google Scholar]
27. Lyon JG, et al. A change detection experiment using vegetation indices. Photogramm Eng Remote Sens. 1998;64(2):143â150.[Google Scholar]
28. Meng Q, et al. A linear mixed-effects model of biomass and volume of trees using Landsat ETM+ images. For Ecol Manage. 2007;244(1â3):93â101. doi: 10.1016/j.foreco.2007.03.056. [CrossRef] [Google Scholar]
29. Salvador R, Pons X. On the applicability of Landsat TM images to Mediterranean forest inventories. For Ecol Manage. 1998;104(1â3):193â208. doi: 10.1016/S0378-1127(97)00264-8. [CrossRef] [Google Scholar]
30. Carreiras JMB, Pereira JMC, Pereira JS. Estimation of tree canopy cover in evergreen oak woodlands using remote sensing. For Ecol Manage. 2006;223(1â3):45â53. doi: 10.1016/j.foreco.2005.10.056. [CrossRef] [Google Scholar]
31. Tomppo E, et al. A sampling design for a large area forest inventory: case Tanzania. Can J For Res. 2014;44(8):931â948. doi: 10.1139/cjfr-2013-0490. [CrossRef] [Google Scholar]
32. Campbell BM. The miombo in transition: woodlands and welfare in Africa. Bogor: Center for International Forestry Research; 1996. [Google Scholar]
33. Lupala ZJ, Lusambo LP, Ngaga YM. Management, Growth, and carbon storage in miombo woodlands of Tanzania. Int J For Res. 2014;2014:11.[Google Scholar]
34. Irons JR, Dwyer JL, Barsi JA. The next Landsat satellite: the Landsat data continuity mission. Remote Sens Environ. 2012;122:11â21. doi: 10.1016/j.rse.2011.08.026. [CrossRef] [Google Scholar]
35. Roy DP, et al. Landsat-8: science and product vision for terrestrial global change research. Remote Sens Environ. 2014;145:154â172. doi: 10.1016/j.rse.2014.02.001. [CrossRef] [Google Scholar]
36. IPCC. Guidelines for National greenhouse gas inventoriesâvolume 4âagriculture, forestry and other land use. Institute for Global Environmental Strategies; 2006.
37. Mauya EW, et al. Modelling aboveground forest biomass using airborne laser scanner data in the miombo woodlands of Tanzania. Carbon Balance Manag. 2015;10(28):16.[PMC free article] [PubMed] [Google Scholar]
38. Dondeyne S, et al. No short cuts to sound forest management: experiences from a participatory survey in Angai Forest, Tanzania. Land. 1998;2(3):181â190.[Google Scholar]
39. Mukama K, Mustalahti I, Zahabu E. Participatory forest carbon assessment and REDD+: learning from Tanzania. Int J For Res; 2012.
40. Tarimo B, et al. Spatial distribution of temporal dynamics in anthropogenic fires in miombo savanna woodlands of Tanzania. Carbon Balance Manag. 2015;10(1):1â15. doi: 10.1186/s13021-015-0029-2.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
41. Cochran WG. Sampling techniques. New York: Wiley; 1977. [Google Scholar]
42. Mugasha WA, Bollandsas OM, Eid T. Relationships between diameter and height of trees in natural tropical forest in Tanzania. Southern Forests. 2014;76(1):65. doi: 10.2989/20702620.2013.869876. [CrossRef] [Google Scholar]
43. Mugasha WA, et al. Allometric models for prediction of above- and belowground biomass of trees in the miombo woodlands of Tanzania. For Ecol Manage. 2013;310:87â101. doi: 10.1016/j.foreco.2013.08.003. [CrossRef] [Google Scholar]
44. Pereira JMC, et al. A simulation analysis of the detectability of understory burns in miombo woodlands. Remote Sens Environ. 2004;93(3):296â310. doi: 10.1016/j.rse.2004.01.009. [CrossRef] [Google Scholar]
45. Masek JG, et al. LEDAPS calibration, reflectance, atmospheric correction preprocessing code, Version 2. Model product. Oak Ridge: Oak Ridge National Laboratory Distributed Active Archive Center. 2013. doi: 10.3334/ORNLDAAC/1146. http://www.daac.ornl.gov.
46. USGS. U. S Geological survey product guide. Provisional Landsat 8 surface reflectance product. Version 1.4. 2015. http://www.landsat.usgs.gov/documents/provisional_l8sr_product_guide.pdf.
47. Rouse JW. Monitoring the vernal advancement of retrogradation of natural vegetation. NASA/GSFC, Type III, Final Report, Greenbelt. 1974. p. 371.
48. Huete AR, et al. A comparison of vegetation indices global set of TM images for EOS-MODIS. Remote Sens Environ. 1997;59(3):440â451. doi: 10.1016/S0034-4257(96)00112-5. [CrossRef] [Google Scholar]
49. Huete AR. A soil-adjusted vegetation index (Savi) Remote Sens Environ. 1988;25(3):295â309. doi: 10.1016/0034-4257(88)90106-X. [CrossRef] [Google Scholar]
50. Qi J, et al. A modified soil adjusted vegetation index. Remote Sens Environ. 1994;48(2):119â126. doi: 10.1016/0034-4257(94)90134-1. [CrossRef] [Google Scholar]
51. Rondeaux G, Steven M, Baret F. Optimization of soil-adjusted vegetation indices. Remote Sens Environ. 1996;55(2):95â107. doi: 10.1016/0034-4257(95)00186-7. [CrossRef] [Google Scholar]
52. Gao BC. NDWIâa normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens Environ. 1996;58(3):257â266. doi: 10.1016/S0034-4257(96)00067-3. [CrossRef] [Google Scholar]
53. Miller JD, Thode AE. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR) Remote Sens Environ. 2007;109(1):66â80. doi: 10.1016/j.rse.2006.12.006. [CrossRef] [Google Scholar]
54. Davidian M, Giltinan DM. Nonlinear models for repeated measurement data. Monographs on statistics and applied probability. Boca Raton: Chapman & Hall/CRC; 1998. [Google Scholar]
55. Littell RC, et al. SAS for mixed models. 2. Cary: SAS Institute, Inc; 2006. [Google Scholar]
56. Diggle PJ, Liang K-Y, Zeger SL. Analysis of longitudinal data. London: Chapman and Hall; 1994. [Google Scholar]
57. Huber PJ. The behavior of maximum likelihood estimates under nonstandard conditions. Proceedings of the fifth Berkeley symposium on mathematical statistics and probability. vol I; 1967. p. 221â33.
58. Liang K-Y, Zeger SL. Longitudnal data analysis using generalized linear models. Biometrika 73; 1986. p. 13â22.
59. White H. A heteroskedasticity-consistent Covariance matrix estimator and a direct test for heteroskedasticty. Econometrica. 1980;48(4):817â838. doi: 10.2307/1912934. [CrossRef] [Google Scholar]
60. McRoberts RE. A model-based approach to estimating forest area. Remote Sens Environ. 2006;103(1):56â66. doi: 10.1016/j.rse.2006.03.005. [CrossRef] [Google Scholar]
61. Avitabile V, et al. An integrated panâtropical biomass map using multiple reference datasets. Global Change Biol; 2016. [PubMed]
62. Digital Globe, Google Earth images 2015. Image Accessed 10 Sept 2015. 2015.
63. Labrecque S, et al. A comparison of four methods to map biomass from Landsat-TM and inventory data in western Newfoundland. For Ecol Manage. 2006;226(1â3):129â144. doi: 10.1016/j.foreco.2006.01.030. [CrossRef] [Google Scholar]
64. Ribeiro NS, et al. Aboveground biomass and Leaf Area Index (LAI) mapping for Niassa Reserve, northern Mozambique. J Geophys Res Biogeosci. 2008;113(G2).
65. Gasparri NI, et al. Assessing multi-temporal Landsat 7 ETM + images for estimating above-ground biomass in subtropical dry forests of Argentina. J Arid Environ. 2010;74(10):1262â1270. doi: 10.1016/j.jaridenv.2010.04.007. [CrossRef] [Google Scholar]
66. Gerylo GR, et al. Empirical relations between Landsat TM spectral response and forest stands near fort simpson, Northwest Territories, Canada. Can J Remote Sens. 2002;28(1):68â79. doi: 10.5589/m02-007. [CrossRef] [Google Scholar]
67. Hall RJ, et al. Modeling forest stand structure attributes using Landsat ETM+ data: application to mapping of aboveground biomass and stand volume. For Ecol Manage. 2006;225(1â3):378â390. doi: 10.1016/j.foreco.2006.01.014. [CrossRef] [Google Scholar]
68. Heiskanen J. Estimating aboveground tree biomass and leaf area index in a mountain birch forest using ASTER satellite data. Int J Remote Sens. 2006;27(5â6):1135â1158. doi: 10.1080/01431160500353858. [CrossRef] [Google Scholar]
69. Myeong S, Nowak DJ, Duggin MJ. A temporal analysis of urban forest carbon storage using remote sensing. Remote Sens Environ. 2006;101(2):277â282. doi: 10.1016/j.rse.2005.12.001. [CrossRef] [Google Scholar]
70. Dube T, Mutanga O. Evaluating the utility of the medium-spatial resolution Landsat 8 multispectral sensor in quantifying aboveground biomass in uMgeni catchment, South Africa. ISPRS J Photogramm Remote Sens. 2015;101:36â46. doi: 10.1016/j.isprsjprs.2014.11.001. [CrossRef] [Google Scholar]
71. Chen Q, Laurin GV, Valentini R. Uncertainty of remotely sensed aboveground biomass over an African tropical forest: propagating errors from trees to plots to pixels. Remote Sens Environ. 2015;160:134â143. doi: 10.1016/j.rse.2015.01.009. [CrossRef] [Google Scholar]
72. Chen XX, et al. Estimating aboveground forest biomass carbon and fire consumption in the U.S. Utah High Plateaus using data from the Forest Inventory and Analysis Program, Landsat, and LANDFIRE. Ecol Ind. 2011;11(1):140â148. doi: 10.1016/j.ecolind.2009.03.013. [CrossRef] [Google Scholar]
73. Næsset E, et al. Mapping and estimating forest area and aboveground biomass in miombo woodlands in Tanzania using data from airborne laser scanning, TanDEM-X, RapidEye, and global forest maps: a comparison of estimated precision. Remote Sens Environ. 2016;175:282â300. doi: 10.1016/j.rse.2016.01.006. [CrossRef] [Google Scholar]
74. Mitchard ET, et al. Uncertainty in the spatial distribution of tropical forest biomass: a comparison of pan-tropical maps. Carbon Balance Manag. 2013;8(1):10. doi: 10.1186/1750-0680-8-10.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
75. Ribeiro N, et al. Monitoring vegetation dynamics and carbon stock density in miombo woodlands. Carbon Balance Manag. 2013;8:11. doi: 10.1186/1750-0680-8-11.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
76. Ryan CM, Williams M, Grace J. Above- and Below-ground carbon stocks in a miombo woodland landscape of Mozambique. Biotropica. 2011;43(4):423â432. doi: 10.1111/j.1744-7429.2010.00713.x. [CrossRef] [Google Scholar]
77. Kashindye A, et al. Multi-temporal assessment of forest cover, stocking parameters and above-ground tree biomass dynamics in miombo woodlands of Tanzania. Afr J Environ Sci Technol. 2013;7(7):611â623.[Google Scholar]
78. Shirima DD, et al. Carbon storage, structure and composition of miombo woodlands in Tanzaniaâs Eastern Arc Mountains. Afr J Ecol. 2011;49(3):332â342. doi: 10.1111/j.1365-2028.2011.01269.x. [CrossRef] [Google Scholar]
79. Kusaga MM. Participatory forest carbon assessment in Angai village land forest reserve in Liwale district, Lindi region, Tanzania. S.U.O. Agriculture, editor. Masters Thesis. Morogoro; 2010. p. 1â167.
80. The Katoomba Group. Characterisation of high ranking districts for the MJUMITA/TFCG REDD project: making REDD work for communities and forest conservation in Tanzania. 2010. http://www.tfcg.org/pdf/. Accessed 30 Jul 2015.
81. Schroeder W, et al. Active fire detection using Landsat-8/OLI data. Remote Sens Environ; 2015.
82. Hill MJ, Aspinall R. Spatial information for land use management. Amsterdam: Gordon & Breach Science Publishers; 2000. [Google Scholar]
Articles from Carbon Balance and Management are provided here courtesy of BioMed Central
Comments are closed.
|
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |