A framework for estimating forest disturbance intensity from successive remotely sensed biomass maps: moving beyond average biomass loss estimates
Carbon Balance and Management volume 10, Article number: 27 (2015)
The success of satellites in mapping deforestation has been invaluable for improving our understanding of the impacts and nature of land cover change and carbon balance. However, current satellite approaches struggle to quantify the intensity of forest disturbance, i.e. whether the average rate of biomass loss for a region arises from heavy disturbance focused in a few locations, or the less severe disturbance of a wider area. The ability to distinguish between these, very different, disturbance regimes remains critical for forest managers and ecologists.
We put forward a framework for describing all intensities of forest disturbance, from deforestation, to widespread low intensity disturbance. By grouping satellite observations into ensembles with a common disturbance regime, the framework is able to mitigate the impacts of poor signal-to-noise ratio that limits current satellite observations. Using an observation system simulation experiment we demonstrate that the framework can be applied to provide estimates of the mean biomass loss rate, as well as distinguish the intensity of the disturbance. The approach is robust despite the large random and systematic errors typical of biomass maps derived from radar. The best accuracies are achieved with ensembles of ≥1600 pixels (≥1 km2 with 25 by 25 m pixels).
The framework we describe provides a novel way to describe and quantify the intensity of forest disturbance, which could help to provide information on the causes of both natural and anthropogenic forest loss—such information is vital for effective forest and climate policy formulation.
Tropical deforestation has been estimated to occur at a rate of 13 million ha per year , with an associated net loss of forest biomass of 1.3 ± 0.7 Pg C year−1 . Remote sensing has been successful at mapping global deforestation [3, 4]. However, deforestation presents a simplified view of forest disturbance that ignores the many graduations of lower intensity, but often widespread, forest disturbance and degradation . Forest disturbance refers to the mechanisms which limit biomass by causing its destruction . The impact of forest disturbance is highly variable, leading to the total or partial loss of biomass through a diverse range of natural (e.g. disease, droughts, fires, herbivory and windstorms) and/or anthropogenic processes (e.g. urbanisation, agriculture, selective logging, and fires). Remotely sensed information on the spatial extent and intensity of forest disturbance would be extremely useful for managers and ecologists in attempts to develop a mechanistic understanding of forest degradation and the processes of forest disturbance [7, 8].
Unfortunately, whilst satellites have the coverage to provide global information on forest disturbance, very few studies have attempted to do so . The focus of remote sensing has remained on mapping deforestation (e.g. ) and arises from a number of technological limitations which include, amongst other factors, observation precision and the relative “directness” of the observation [10–12]. Ideally remote sensing measurements would have a direct dependence on the quantity being estimated; however current optical satellite derived estimates of forest disturbance are not able to achieve this ideal. Instead, optical satellite measurements of forest disturbance rely on observing physical properties that are expected to correlate with disturbance (e.g. changes in leaf area that can be detected by satellites due to a change in the absorption in the chlorophyll spectral bands). When these correlations change, indirect measures cannot be expected to provide a robust estimate of disturbance.
As Synthetic aperture radar (SAR) is sensitive to the structural properties of forests it can be thought of as a more direct measure of forest biomass than the passive optical alternatives [8, 13, 14]. SAR provides one of the only viable data sources currently available for global monitoring of forest disturbance at moderate spatial resolutions. However, the low signal to noise ratio of SAR leads to poor precision and large uncertainty in the estimated biomass of individual pixels  and this equates to high levels of uncertainty when detecting biomass change at the pixel level. Filtering can be used to improve the signal to noise ratio of SAR, but this improvement comes at the expense of a reduction in the effective spatial resolution of the biomass change estimates. Therefore remote sensing can provide regional estimates of biomass change or fine scale maps of deforestation, but it cannot yet be said to truly determine the intensity of forest degradation across scales .
In this study we set out a novel framework for quantifying the intensity of forest disturbance from successive remotely sensed biomass maps. This approach can describe both focused high intensity forest loss and low intensity, widespread degradation. We provide an example methodology to exploit this framework using SAR data with realistic observation errors. We explore the strengths and limitations of this approach using an observation system simulation experiment.
Results and discussion
A framework for quantifying the intensity of forest disturbance
The framework that we propose describes forest disturbance for an area in which all satellite observation pixels are assumed to experience the same disturbance regime. That is, whilst each pixel might or might not be disturbed, each pixel within each area has the same probability of being disturbed and the same relative loss of biomass when disturbed. We consider an ensemble of n remote sensing pixels, which are observed on two successive dates (t = 1 and t = 2). This approach builds on an earlier framework set out in Williams  that used the biomass distribution at a single point in time. The first advantage of this new approach is that it allows ensemble statistics to be calculated, negating the limitations of poor signal to noise of individual biomass pixels by estimating ensemble’s mean fractional loss of biomass per year (E M). The second advantage of the framework is that it permits E M to be split into two factors for the ensemble: the probability of disturbance per year for each ensemble pixel member (E P) and the fractional loss of biomass per disturbance for each disturbed ensemble pixel member (E I), Eq. 1.
Both E M and its factors E P and E I can take values ranging from 0–1. Where, for example, an E P = 0.05 would imply a 1 in 20 chance of each pixel being disturbed in a year. An E I = 0.2 implies that, if disturbed, a pixel will lose 20 % of its biomass. Combining these example factors would lead to the expectation of E M = 0.01, or a 1 % reduction in mean biomass for the ensemble. The ensemble size can be picked to balance the competing requirements of high precision on estimates of E M, E P and E I and meeting the assumption of a continuous disturbance regime.
The inclusion of the factors E P and E I allows a flexible description of biomass loss, without the arbitrary distinction between deforestation and lower intensity forest disturbance. Forest disturbance (other than deforestation) can be represented by the parameter space 0 < E M < 1, 0 < E I < 1 and 0 < E P ≤ 1 (Fig. 1). There are also several special cases: the total deforestation of the ensemble area (E M = 1, E I = 1 and E P = 1); partial deforestation within the area (0 < E M < 1, E I = 1 and 0 < E P < 1); and no disturbance (E M = 0).
Describing forest disturbance using this framework has a number of advantages over traditional descriptors as it avoids the need for an arbitrary threshold for a forest cover loss used in other studies, e.g. . In turn, this allows for a more nuanced description of forest disturbance than is possible with categorical land cover classes, or measures of forest cover.
Using the framework to estimate biomass loss and disturbance intensity
Using the disturbance framework, our observation system simulation experiments (OSSE) show it is possible to robustly estimate the mean biomass loss for the ensemble (E M) and also the disturbance regime, as described by E P and E I (Fig. 2; Table 1). E M is better constrained than either E I or E P, reflecting the direct impact of E M on the mean of the biomass distribution versus the more variable impacts of the E I or E P on the second and third order moments: standard deviation and skew. With larger mean disturbances the constraints on E I or E P improve, presumably due to the larger number of pixels effected and/or greater impact on the observed biomass.
How the framework is affected by observation bias
It is highly likely that biomass maps used with the framework will not only have random noise, but also systematic observation errors , therefore we test the framework’s sensitivity to bias. The inclusion of a realistic observation bias of ±160 gC m−2 (±1.6 tC ha−1)  had the largest impacts on the estimates of E M, which showed a bias of 0.004 (0.4 %). However E I and E P were less affected (Fig. 3; Table 1) and it was still possible to distinguish high and low intensity disturbance regimes. The bias in the estimate for E M is consistent with change in the mean biomass that is implied by a bias of 160 gC m−2, and would therefore apply to any other approach to estimating biomass loss.
What is the optimal ensemble size?
The accuracy of predictions using the framework are best when considering an ensemble size of n ≥ 1600, whilst for n ≤ 400, the precision of our estimates drops rapidly (Fig. 4; Table 1). Therefore, for the 25 m pixels typical of current SAR biomass estimates, the recommended ensemble size covers an area of at least 1 km2 or 100 ha. The implication is that it must be reasonable to assume that the disturbance regime is constant across areas of at least 100 ha. We expect this minimum area to be robust, at least for ALOS PALSAR data, but the sensitivity to ensemble size should checked in each new study. It is worth noting that there is no requirement for these areas to be rectangular, or even contiguous. The grainy texture evident in the estimates from the analyses is due to the simulation of stochastic disturbances (Figs. 2, 3, 4). This graining is reduced with increasing ensemble size, but can be further reduced, at significant computational expense, by averaging repeat runs of the maximum likelihood estimation.
Alternative formulations of the framework
We formulated the framework in terms of a fractional change in biomass for both E M and E I. This is not the only possible formulation, nor is it necessarily the optimal for all situations; however, it is a mathematically simple approach to maintaining a positive (i.e. plausible) biomass for each pixel. The implicit assumptions of the formulation are: 1) that all pixels are equally likely to be disturbed, and 2) that when disturbed a fixed fraction is lost, irrespective of the starting biomass. It is possible to imagine scenarios that would not be well described by our scheme where (say) selective logging only targets the largest trees (i.e. E P is high for pixels above a threshold biomass and zero for all others). These scenarios do not contradict our first assumption of common disturbance regime for the ensemble and it should be possible to restate the parameters used in the framework to accommodate a particular set of assumptions about the disturbance regime. However the inclusion of more complicated mathematical representations and, specifically, more parameters, can be expected to increase the challenge of estimating the parameters of any new framework formulation and decrease its general applicability.
Limitations to the framework
Finally there are two notable caveats: Firstly, the success of the approach is strongly tied to the ability to characterise the random error of biomass estimates. The design of the analysis mitigates some of the errors as biomass errors that are consistent between the two biomass maps will be removed by this differencing. Whilst Ryan et al. (2012) report normal errors in the biomass domain, other studies assume errors will actually be normal in the dB (i.e. log10) domain  which would result in log-normal, asymmetric errors on biomass estimates. We therefore reran the synthetic experiments with 0.5 dB error and achieved similar results. The second caveat is that the biomass differencing approach assumes the carbon model (A-DALEC) is unbiased. This assumption is, by definition, valid in a observation system simulation experiment (OSSE). However in practice the productivity of the ecosystem is not known perfectly, and so estimates of E M can be expected to be biased. However the assumption is not as crude as it at first might seem: ignoring the production term (as is implicitly done by most deforestation algorithms) makes the assumption that the forest otherwise in steady state. This issue is likely to be more severe in field sites where less information is available on which to provide independent estimate of the rates of aboveground biomass accumulation, e.g. .
From a management and policy perspective it is important to be able to distinguish between the different intensities of biomass loss as they may be associated with different disturbances mechanisms (e.g. low intensity disturbances are likely driven by a need for timber or fuel, and high intensity disturbances driven by a need for agricultural land) . However, current estimates of deforestation and biomass loss are not adequate for estimating lower intensity forest disturbance. Theoretically, high resolution biomass loss estimates could provide fine-scale estimates of forest degradation, but the current precision is not adequate and there is no immediate prospect of this changing, partly due to the speckle and other noise in SAR imagery . The framework that we have described is a pragmatic representation of forest degradation that uses ensemble statistics to mitigate the poor precision of current SAR biomass estimates [13, 17, 20]. SAR is expected to remain a key technique for global mapping of forest biomass, and the framework we propose is compatible with the upcoming BIOMASS and new L-band satellites . Using an OSSE we have shown that is possible to robustly estimate the parameters of the framework to describe forest disturbance, provided that two successive biomass maps, separated by at least one year, are available. We suggest that using a similar framework will allow remote sensing studies to provide more relevant constraints on estimates of land-use change.
Estimating the framework parameters
A number of approaches could be taken to estimating E M, E I and E P, we chose to use maximum likelihood estimation. We identify the combinations of E I and E P that allow a simulations to most closely match the observed changes in the biomass distribution from time t = 1 to time t = 2, for an ensemble of n pixels.
The observed biomass (O t i ) of the ith pixel at time t is sum of the actual biomass (B t i ) and the measurement noise (N t i ), Eq. 2. The observed biomass at time t + 1, follows the same logic (Eq. 3). The actual biomasses are related via the growth G t i and disturbance D t i of each pixel in the ensemble (Eq. 4).
Similarly we are able to simulate the observed biomass at time t (Os t i ) based on the simulation biomass (Bs t i ) added to the simulated noise (Ns t i ), Eq. 5. The observed biomass at time t (Os t+1 i ) is calculated from the simulation biomass (Bs t i ) through the sum of the simulated growth (Gs t i ), the simulated disturbance (Ds t i ), and the simulated noise (Ns t+1 i ), Eq. 6.
where Ds t i is related to the disturbance parameters in Eq. 7.
We model Gs t i using the A-DALEC model, a biogeochemical model of carbon cycling in forests, which runs at an annual time step (see Additional file 1 A for full details) [15, 21]. Ns t i and Ns t+1 i are based on observed pixel uncertainty in SAR biomass estimates . We cannot base Bs t i directly on the observed biomass O t i , as the observation noise N t i cannot be determined (and thus removed) for each pixel, and so the A-DALEC model is used to spin-up an estimate of Bs t i , see Additional file 1 section B for full details. The equations above assume that t = 1 and t = 2 are separated by 1 year. When the separation is more than 1 year, as in the case of our OSSE, the base Bs t i should have the annual growth and disturbance applied once for each year of separation.
95 % confidence intervals for the factors E P and E I are found using a maximum likelihood estimation routine, see Additional file 1 C for full details. This routine finds the most likely combinations of E P and E I based on minimising the differences between the observed and simulated ensemble distributions; ΔO = O t+1 − O t and ΔOs = Os t+1 − Os t.
Testing the framework parameter estimates
To allow the technique to be assessed against a known ‘truth’ we use a set of observation system simulation experiment (OSSE). In the OSSE we simulate two successive (noisy) biomass observations which have observation noise and a definable disturbance regime applied. The OSSE observations are derived by spinning up the A-DALEC model to match ALOS observations within the Gorongosa and Nhamatanda districts of Sofala province in central Mozambique. The area is dominated by dry miombo woodland, the dominant woodland type in Southern Africa, see Additional file 1 B. The first OSSE observation, t = 1, is taken at the end of the 1st year, the second observation, t = 2, is taken at the end of the 5th year, i.e. separated by 4 years. Separate spin-ups are performed for the ΔO and ΔS simulations. Three tests are performed; the first to assess the ability of the approach to identify E M, E I and E P, the second to assess the impact of observation bias and the third to determine impact of ensemble size.
Determining the accuracy of the framework parameter estimates: test 1
A set of nine analyses were used to test the ability of our approach to identify E M, E I and E P. Each analysis was performed for the same size ensemble of n = 1600 pixels, which equates to an area of 1000 by 1000 m given the 25 by 25 m pixel size of the ALOS biomass maps . We perform OSSEs with four levels of mean disturbance E M = 0, 0.018, 0.045, 0.090 and 0.180. These disturbance fractions equate to annual biomass loss rates of 0, 1.8, 4.5, 9 and 18 %. For each nonzero E M we simulate a high intensity (E I ≫ E P) and low intensity (E I ≪ E P) disturbance scenario. In the high intensity cases only a few pixels are disturbed, but when disturbed 90 % of the biomass is lost (i.e. E I = 0.9). In the low intensity cases, 90 % of the pixels are disturbed (i.e. E P = 0.9), but each disturbance results in a small loss of biomass at the pixel level. The various combinations of factors are shown in Table 1.
Determining the impact of observation bias on the framework parameter estimates: test 2
To simulate the impact of observation bias we perturb the synthetic experiments using the bias observed in the actual biomass maps . Only the first biomass observation (i.e. at time t = 1) is biased by −160, 0, and 160 gC m−2. The impact of the bias is assessed using both a low intensity case (i.e. E I = 0.05 and E P = 0.9), and a high intensity case (i.e. E I = 0.9 and E P = 0.05). In both cases the OSSE has an E M = 0.045 and an ensemble size of n = 1600 pixels.
Determining the impact of ensemble size on the framework parameter estimates: test 3
The impact of varying the ensemble size is tested for n = 6400, 1600, 400, and 100. These ensembles equate to areas of 2000 by 2000 m (6400 pixels), 1000 by 1000 m (1600 pixels), 500 by 500 m (400 pixels), and 250 by 250 m (100 pixels). The areas are centred on the same location. We analyses these areas for a single low intensity disturbance regime E M = 0.045, E I = 0.05 and E P = 0.9.
FAO. Global forest resources assessment 2010: main report. FAO forestry paper. Rome: Food and agriculture organization of the United Nations; 2010. p. 2010.
Pan Y, Birdsey RA, Fang J, Houghton R, Kauppi PE, Kurz WA, et al. A large and persistent carbon sink in the world’s forests. Science. 2011;333(6045):988–93. doi:10.1126/science.1201609.
Achard F, Beuchle R, Mayaux P, Stibig H-J, Bodart C, Brink A, et al. Determination of tropical deforestation rates and related carbon losses from 1990 to 2010. Glob Chang Biol. 2014;20(8):2540–54. doi:10.1111/gcb.12605.
Hansen MC, Stehman SV, Potapov PV. Quantification of global gross forest cover loss. Proc Natl Acad Sci USA. 2010;107(19):8650–5. doi:10.1073/pnas.0912668107.
Mertz O, Müller D, Sikor T, Hett C, Heinimann A, Castella J-C, et al. The forgotten D: challenges of addressing forest degradation in complex mosaic landscapes under REDD+. Geografisk Tidsskrift-Danish J Geogr. 2012;112(1):63–76. doi:10.1080/00167223.2012.709678.
Grime JP. Evidence for existence of three primary strategies in plants and its relevance to ecological and evolutionary theory. Am Nat. 1977;111(982):1169–94. doi:10.1086/283244.
Agarwal DK, Silander JA Jr, Gelfand AE, Dewar RE, Mickelson JG Jr. Tropical deforestation in Madagascar: analysis using hierarchical, spatially explicit, Bayesian regression models. Ecol Model. 2005;185(1):105–31. doi:10.1016/j.ecolmodel.2004.11.023.
Ryan CM, Berry NJ, Joshi N. Quantifying the causes of deforestation and degradation and creating transparent REDD + baselines: a method and case study from central Mozambique. Appl Geogr. 2014;53:45–54. doi:10.1016/j.apgeog.2014.05.014.
Herold M, Roman-Cuesta RM, Mollicone D, Hirata Y, Van Laake P, Asner GP, et al. Options for monitoring and estimating historical carbon emissions from forest degradation in the context of REDD+. Carbon Balance Manag. 2011;6(1):13. doi:10.1186/1750-0680-6-13.
Hill TC, Williams M, Bloom AA, Mitchard ETA, Ryan CM. Are inventory based and remotely sensed above-ground biomass estimates consistent? PLoS One. 2013;8(9):e74170. doi:10.1371/journal.pone.0074170.
Woodhouse IH, Mitchard ETA, Brolly M, Maniatis D, Ryan CM. Radar backscatter is not a ‘direct measure’ of forest biomass. Nature Clim Chang. 2012;2(8):556–7.
Saatchi S, Ulander L, Williams M, Quegan S, LeToan T, Shugart H, et al. Forest biomass and the science of inventory from space. Nature Clim Chang. 2012;2(12):826–7.
Ryan CM, Hill TC, Woollen E, Ghee C, Mitchard E, Cassells G, et al. Quantifying small-scale deforestation and forest degradation in African woodlands using radar imagery. Glob Chang Biol. 2012;18:243–57. doi:10.1111/j.1365-2486.2011.02551.x.
Le Toan T, Quegan S, Davidson MWJ, Balzter H, Paillou P, Papathanassiou K, et al. The BIOMASS mission: mapping global forest biomass to better understand the terrestrial carbon cycle. Remote Sens Environ. 2011;115(11):2850–60. doi:10.1016/j.rse.2011.03.020.
Williams M, Hill TC, Ryan CM. Using biomass distributions to determine probability and intensity of tropical forest disturbance. Plant Ecol Divers. 2013;6(1):1–13. doi: 10.1080/17550874.2012.692404.
Harris NL, Brown S, Hagen SC, Saatchi SS, Petrova S, Salas W, et al. Baseline map of carbon emissions from deforestation in tropical regions. Science. 2012;336(6088):1573–6. doi:10.1126/science.1217962.
Mitchard ETA, Saatchi SS, Woodhouse IH, Nangendo G, Ribeiro NS, Williams M et al. Using satellite radar backscatter to predict above-ground woody biomass: a consistent relationship across four different African landscapes. Geophys Res Lett. 2009;36. doi: 10.1029/2009GL040692.
Williams M, Ryan CM, Rees RM, Sarnbane E, Femando J, Grace J. Carbon sequestration and biodiversity of re-growing miombo woodlands in Mozambique. For Ecol Manage. 2008;254(2):145–55. doi:10.1016/j.foreco.2007.07.033.
Frolking S, Palace MW, Clark DB, Chambers JQ, Shugart HH, Hurtt GC. Forest disturbance and recovery: a general review in the context of spaceborne remote sensing of impacts on aboveground biomass and canopy structure. J Geophys Res Biogeosci. 2009;114. doi: 10.1029/2008jg000911.
Williams M, Hill TC, Ryan CM, Peylin P. Final report BIOMASS level 2 product to flux study. ESA contract No 4000102042/10/NL/CT. Nordwijk: The European Space Research and Technology Centre; 2011.
Williams M, Schwarz PA, Law BE, Irvine J, Kurpius MR. An improved analysis of forest carbon dynamics using data assimilation. Glob Chang Biol. 2005;11(1):89–105.
The study was designed by all authors. TH performed the analysis and produced the figures and tables and wrote the manuscript. All co-authors contributed substantially to the interpretation of results and revision of the manuscript. All authors read and approved the final manuscript.
This work was funded by ESA (Contract No 4000102042/10/NL/CT), the NERC National Centre for Earth Observation, the NERC CarbonFusion project (NE/D000874/1), the Mpingo Conservation and Development Initiative and the EU Framework 7 I-REDD + project. ESA and JAXA provided the ALOS PALSAR imagery through C1P.7493 and the 4th Research Agreement for the Advanced Land Observing Satellite-2. The authors would like to thank Klaus Scipal and Shaun Quegan for advice during the study. Finally we would like to thank the three anonymous reviewers for their helpful comments.
The authors declare that they have no competing interests.
About this article
Cite this article
Hill, T.C., Ryan, C.M. & Williams, M. A framework for estimating forest disturbance intensity from successive remotely sensed biomass maps: moving beyond average biomass loss estimates. Carbon Balance Manage 10, 27 (2015). https://doi.org/10.1186/s13021-015-0039-0
- Remote sensing