Next Article in Journal
The 2014 Mw 6.1 Ludian Earthquake: The Application of RADARSAT-2 SAR Interferometry and GPS for this Conjugated Ruptured Event
Previous Article in Journal
Eddy Detection in HF Radar-Derived Surface Currents in the Gulf of Naples
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 12
Issue 1
10.3390/rs12010098
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Landsat Time-Series for Estimating Forest Aboveground Biomass and Its Dynamics across Space and Time: A Review

by
Trung H. Nguyen
1,2,3,4,*,
Simon Jones
1,2,
Mariela Soto-Berelov
1,2,
Andrew Haywood
5 and
Samuel Hislop
1,2
1
Remote Sensing Centre, School of Science, RMIT University, Melbourne, VIC 3000, Australia
2
Cooperative Research Centre for Spatial Information (CRCSI), Carlton, VIC 3053, Australia
3
Faculty of Resource Management, Thai Nguyen University of Agriculture and Forestry, Thai Nguyen, Vietnam
4
Department of Environment, Land, Water and Planning, Melbourne, VIC 3000, Australia
5
European Forest Institute, 08025 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(1), 98; https://doi.org/10.3390/rs12010098
Submission received: 30 November 2019 / Revised: 17 December 2019 / Accepted: 24 December 2019 / Published: 27 December 2019
(This article belongs to the Section Forest Remote Sensing)
Browse Figures
Versions Notes

Abstract

:
The free open access data policy instituted for the Landsat archive since 2008 has revolutionised the use of Landsat data for forest monitoring, especially for estimating forest aboveground biomass (AGB). This paper provides a comprehensive review of recent approaches utilising Landsat time-series (LTS) for estimating AGB and its dynamics across space and time. In particular, we focus on reviewing: (1) how LTS has been utilised to improve the estimation of AGB (for both single-date and over time) and (2) recent LTS-based approaches used for estimating AGB and its dynamics across space and time. In contrast to using single-date images, the use of LTS can benefit forest AGB estimation in two broad areas. First, using LTS allows for the filling of spatial and temporal data gaps in AGB predictions, improving the quality of AGB products and enabling the estimation of AGB across large areas and long time-periods. Second, studies have demonstrated that spectral information extracted from LTS analysis, including forest disturbance and recovery metrics, can significantly improve the accuracy of AGB models. Throughout the last decade, many innovative LTS-based approaches for estimating forest AGB dynamics across space and time have been demonstrated. A general trend is that methods have evolved as demonstrated through recent studies, becoming more advanced and robust. However, most of these methods have been developed and tested in areas that are either supported by established forest inventory programs and/or can rely on Lidar data across large forest areas. Further investigations should focus on tropical forest areas where inventory data are often not systematically available and/or out-of-date.

Graphical Abstract

1. Introduction

Landsat satellites are unique in that they have created the longest continuously acquired, consistent space-based and moderate-resolution data collection since 1972. To date, nine Landsat missions have been developed and operated by the United State Geological Survey (USGS). While the Landsat missions 1–4, launched during 1972–1982, acquired data for several years (5–10 years), Landsat 5 delivered high quality, global data across the Earth’s land surface for nearly 29 years (1984–2013). Following the failure of Landsat 6 in 1993, Landsat 7 and Landsat 8 were launched in 1999 and 2013, respectively, and continue to provide satellite data [1]. The new Landsat 9 is being developed, with a scheduled launch by December 2020 [2]. The entire historical Landsat archive has been opening for public access since 2008 [3]. As such, the Landsat archive has become one of the most valuable and cost-effective remotely sensed data sources supporting worldwide land/forest research and monitoring activities. Landsat time-series (LTS) data provide two essential elements for ecosystem monitoring: spatial and temporal dimensions. The spatial dimension is at an adequate scale that allows the capture of both natural and anthropogenic impacts across large areas. The temporal dimension enables retrospective analyses and characterisation of changes across more than 4 decades (slightly more or less in certain areas pending on data availability) captured by successive Landsat sensors [4].
Forest biomass is a key focus of environmental monitoring concerning global biogeochemical cycles and loss of biodiversity. Field methods are the most accurate approach to estimate forest biomass. However, they are unsuitable for calculating biomass across large areas and time periods due to their lack of spatial and temporal coverage, especially for large jurisdictions or remote areas [5,6]. Airborne lidar surveys (ALS) have also been used to map forest aboveground biomass (AGB) [7,8]. However, ALS generally tends to represent only a partial sample in space and time [9]. As an alternative, multispectral satellite-based Earth Observations, such as Landsat imagery, can be used for estimating forest biomass across large areas due to their wide spatial and temporal coverage. Forest AGB and Landsat spectral data are highly correlated [10,11]. AGB and its dynamics are often estimated by combining Landsat images with reference data extracted from field inventory plots. Traditionally, such studies often only used single-date Landsat images to estimate forest AGB at a single or limited number of points in time [12]. The free data policy of the Landsat archive has opened a new era for forest applications, deriving many variables, including AGB [4]. LTS has been increasingly used for large area estimation of AGB over time. This paper aims to provide the remote sensing and forest research communities with a critical review of recent approaches utilising LTS for estimating AGB and its dynamics across space and time. In particular, we address the following questions:
  • How has LTS been utilised to improve the estimation of AGB?
  • What LTS-based approaches have been demonstrated as useful for estimating AGB and its dynamics across space and time?
We first provide a brief review of preprocessing and forest change detection methods for LTS and then focus on answering the above two questions. The first question is addressed by reviewing studies using LTS for estimating AGB, for both a single-date and multiple-dates. The second question is answered by reviewing only studies using LTS for estimating and characterising AGB and its dynamics over time.

2. Advanced Preprocessing and Change Detection Methods for LTS

Recent studies have included comprehensive reviews on preprocessing and change detection methods for LTS [13,14,15]. For example, Hansen and Loveland [15] reviewed methods for processing LTS for large area land cover change monitoring. Zhu [13] reviewed frequencies, preprocessing and algorithms for change detection using LTS. In this section, we briefly summarise recent methods used for LTS preprocessing and forest change detection that facilitate estimation of AGB and its dynamics.

2.1. Robust Preprocessing Methods

Throughout the last decade, many advanced methods have been developed for preprocessing LTS data in order to provide high quality inputs for forest modelling processes. Preprocessing a LTS often includes the following steps: radiometric, geometric and atmospheric corrections, clouds/shadows detection and masking, inter-sensor harmonization, and image compositing. Landsat Level 1 terrain-corrected images have been often used to form a time-series [13]. Radiometric and geometric corrections have been performed on these products, making them suitable for time-series analysis [16]. During 2016–2018, the Landsat archive was reorganised into a formal tiered data collection structure to effectively support long-term LTS “stacking” and analysis (over 40 years) [17]. The new structure manages all Level 1 products in a consistent archive with known data quality. All Landsat acquisitions are stratified into three categories based on data quality: Tier 1, Tier 2, and Real-Time. Tier 1 includes Level 1 Precision and Terrain (L1TP) corrected images that have the highest geometric and radiometric quality (an image-to-image registration accuracy of ≤12-m Root Mean Square Error (RMSE), USGS [17]). Images in this Tier are considered the most suitable for time-series applications. Unfortunately, a large amount of Landsat 1–5 Multispectral Scanner (MSS) images cannot achieve the Tier 1 geometry specification due to less accurate orbital information, insufficient ground control, and other factors including a lack of MSS observations [18]. As a result, LTS relying on Tier 1 products may contain gaps due to the unavailability of high-quality MSS images across some areas and time periods. To address this issue, some studies improved the initial geo-registration of MSS images in Tier 2 (i.e., L1TP with RMSE > 12 m, Systematic Terrain (L1GT) and Geometric Systematic (L1GS) corrections) using some external procedures, which may include some manual steps [19].
Traditionally, atmospheric correction has been conducted using various approaches such as: relative normalization (adjusting the radiometry of other images based on a reference image) [20], the Dark-Object Subtraction (DOS) [21], and the Satellite Signal in the Solar Spectrum (6S) [22]. However, these methods often require many manual and complex steps which hinder the processing of large amount of LTS [13]. Recently, two algorithms are being implemented for automatic atmospheric correction: the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) [23] for Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+), and the Land Surface Reflectance Code (LaSRC) [24] for Landsat 8 Operational Land Imager (OLI) images. LEDAPS and LaSRC have been used to generate surface reflectance images from Landsat Level 1 products that are currently provided on the USGS archive for public access (identified as Landsat Level 2 products). These surface reflectance products have been increasingly used in recent LTS applications as they are robustly pre-processed [13]. Another benefit of using surface reflectance images is that they are provided with high quality cloud/shadow detection products. Cloud/shadow in LEDAPS and LaSRC products are detected using the CFMask algorithm, originally known as FMask developed by Zhu and Woodcock [25]. This algorithm is widely confirmed to provide a higher accuracy of cloud/shadow detection than other methods. While surface reflectance and cloud/shadow masks for TM, ETM+ and OLI images are consistently created and available on the USGS archive, additional preprocessing steps are required to derive these data from MSS images. MSS surface reflectance images are often calculated using the COST method [26], while cloud/shadow can be identified using the MSS clearview-mask (MSScvm) algorithm [27].
A challenge of LTS analysis is harmonising the different spectral, spatial, and radiometric resolutions of Landsat sensors. Most LTS studies utilise images acquired by TM and ETM+ sensors (1984–present) for their time-series analyses since they have similar spatial and spectral characteristics [19]. Utilising the full temporal record of Landsat data (including MSS and OLI data) is valuable for characterising forest dynamics, but it requires additional processing to normalise and incorporate multi-sensor images into a consistent time-series [12,28]. Several inter-sensor harmonization approaches have been recently demonstrated to expand the utility of the Landsat archive for forest dynamics applications. Vogeler, Braaten [19] presented an automated method, called LandsatLinkr, for calibrating Landsat MSS and OLI images to the spatial and spectral profiles of TM and ETM+ images. LandsatLinkr normalises MSS images to the TM images by modelling TM Tasseled Cap (TC) indices from the four MSS surface reflectance bands using multiple linear regression. Regression models are developed based on the relationship between coincident MSS and TM images simultaneously acquired by Landsat 4–5 between 1982 and 1992. They are then applied to normalise all MSS acquired by Landsat 1–3 (back to 1972). A similar approach is deployed in LandsatLinkr for harmonising OLI images to ETM+ images.
Given the limited availability of cloud-free Landsat data in many forested areas, image compositing is an important process for supporting AGB predictions, especially temporal predictions (see Section 3). This process leverages the extensive LTS data to create cloud-free and radiometrically consistent composite images that are spatially contiguous over large areas. Image compositing is also useful for reducing data volume and minimizing atmospheric influences. Since the opening of the Landsat archive in 2008, several pixel-based approaches have been developed for generating composite images from Landsat data. A robust method commonly used in LTS applications is the Best Available Pixels (BAP) compositing, developed by White, Wulder [29]. This method determines the BAP within a compositing window/interval (e.g., annual, season or month) using pixel-scoring functions including pixel-based scoring for sensor, acquisition date, distance to cloud/shadow, and atmospheric opacity [30]. Zhu, Woodcock [31] demonstrated an approach for generating daily clear-sky synthetic images from all available Landsat data by developing time-series models for each pixel and each spectral band. Other approaches such as using maximum, mean or median values of Landsat bands and indices are also used for Landsat image compositing [32,33,34].

2.2. Vegetation Change Detection Using LTS

Having knowledge of the spatial and temporal patterns of forest disturbance and recovery is crucial to improving the accuracy of AGB estimates and associated dynamics (Section 3 and Section 4). Forest disturbances can be abrupt or subtle events caused by natural hazards such as wildfires, drought, and insects or human activities such as logging and mining [35]. In contrast, forest recovery is an ongoing process that can occur naturally through vegetation succession or artificially through forest management [36]. The ability of LTS to capture changes occurring in forests has been widely demonstrated, with various methods developed to characterise forest disturbance and recovery. In a recent review, Zhu [13] categorized LTS change detection methods into six broad groups based on the algorithms used for detecting changes (e.g., differencing, thresholding, trajectory classification, regression, statistical boundary and temporal segmentation). Of these, thresholding, statistical boundary and temporal segmentation have been the most commonly used. Table 1 lists some popular algorithms based on those three methods. These pixel-based algorithms have proven to be capable of detecting changes across large areas and long time-periods. In addition, they are fully automated, publicly available and well documented algorithms.
The VCT is a threshold-based algorithm. It calculates IFZ as a forest probability index for each Landsat pixel in a time-series and determines an IFZ threshold to detect forest disturbance (pixels with an IFZ value smaller than the defined threshold are considered as disturbance) [37]. VCT provides information about abrupt disturbances on an annual scale using an annual LTS. In contrast, CCDC and BFAST Monitor estimate a statistical boundary for the entire time-series and define changes as significant differences from the boundary [25,39]. Both of these algorithms utilize all available Landsat images and are designed to detect both abrupt and gradual forest change. NDVI is often used with BFAST Monitor while CCDC runs on all Landsat spectral bands. LandTrendr is a temporal segmentation algorithm utilising an annual LTS. This algorithm fits the spectral trajectory of each pixel to a series of straight-line segments to capture multiple change events as well as the overall trend of the trajectory over time. LandTrendr is capable of detecting both abrupt and gradual changes at an annual timescale. Although NBR is often used as the change index for LandTrendr, other vegetation indices (such as TC components and NDVI) and spectral bands have also been deployed [41,42]. The robustness of LandTrendr has been demonstrated across various contexts including different geographical regions, forest types and disturbance regimes [13,14]. The algorithm has been recently implemented on the Google Earth Engine platform, offering a consistent and quick approach for mapping forest disturbance and recovery from LTS data [43].
The utility of a change detection technique is dependent on its ability to detect and distinguish real changes (or changes of interest) from apparent changes caused by variations in background signal or noise [42]. Studies have demonstrated that high-impact or stand-clearing disturbances and subsequent recovery can be accurately detected from LST analysis. Whereas, changes caused by medium/low-impact disturbances (such as selective logging and droughts) are more difficult to distinguish and characterise [44]. Therefore, the selection of the most appropriate change detection algorithm for a specific application is challenging and dependent on multiple factors including target changes and the temporal density of LTS used. VTC mainly targets abrupt, high magnitude forest disturbances while CCDC and LandTrendr target a broader range of forest changes (both discrete and gradual changes, and both disturbance and recovery processes). Generally, methods requiring a higher density of LTS (i.e., CCDC and BFAST Monitor) are often computationally expensive and require long-processing times and substantial data storage. In many instances, an annual LTS should be the most suitable temporal scale for mapping forest changes as it allows for the tracking of both abrupt and gradual change without requiring as much processing effort. But regardless of the choice, a benefit of using automate change detection algorithms is that they can robustly produce a suite of spectral change metrics representing trends of forest disturbance and recovery. Change metrics are often utilised in many applications such as attributing disturbance levels and causal agents and improving forest biomass estimates [12,45,46,47].

3. How Has LTS Been Utilised to Improve the Estimation of AGB?

Since the opening of the Landsat archive in 2008, LTS data have been increasingly used for estimating forest AGB over large areas, and not only for a static date but also for multiple points in time. The contribution of LTS to improvements in forest AGB estimates has been demonstrated in a wide range of contexts. In Table 2, we summarise the benefits of using LTS for forest AGB estimation into two broad categories: (1) filling spatial and temporal data gaps in AGB predictions and (2) improving the accuracy of AGB modelling.

3.1. Filling Spatial and Temporal Data Gaps in AGB Predictions

Forest AGB is commonly estimated using high quality Landsat images (i.e., cloud-free observations) from small areas e.g., [65,66,67,68,69,70,71,72,73,74]. However, when larger areas and long time-periods are desired, it can be challenging to find cloud-free Landsat data, resulting in spatial and temporal gaps in AGB prediction maps. This prohibits creation of wall-to-wall AGB maps over large areas (i.e., across multiple Landsat scenes) and/or multiple points in time. LTS data provide an opportunity to overcome these limitations by computing image composites (Section 2.1). Image compositing methods rely on selecting representative pixels (e.g., via BAP) from a series of images rather than the best available scene [29]. Landsat image composites are increasingly used to support the development of spatially and temporally complete predictions of forest AGB (Table 2). Zald, Wulder [48] indicated that Landsat pixel-based composites are important for improving the quality of AGB predictions over large areas, regardless of the modelling technique and training data used. Matasci, Hermosilla [49] demonstrated the utility of Landsat composites for extending localised concept studies to the national scale by creating spatially explicit estimates of AGB across 552 million ha of Canadian boreal forests. The authors also extended their work to estimate AGB dynamics across three decades (1984–2016) [51]. Despite the significant increase in the number of cloud-free observations, following the launch of Landsat 8 in 2013 (and Landsat 9 expected in 2020), LTS-based image composites will remain essential for forest AGB monitoring across large areas and long-time periods, especially for regions with persistent cloud and snow cover [29,48].
Pixel-based temporal fitting is another process that can be performed on LTS to fill (or further fill) spatial and temporal data gaps [48,58,75]. This process results in a “smoothed” or “stabilized” temporal trajectory for each Landsat pixel, filling data gaps (i.e., no data observations) and removing noisy pixels (due to unscreened cloud/shadow). Pixel-based temporal fitting can be used for filling not only spatial gaps but also temporal gaps resulting from the unavailability of Landsat data at some points in time. For example, Deo, Russell [58] applied a pixel-level curve fitting approach to a LTS of 17 cloud-free images representing 26 years (1986–2011) to obtain fitted images for years where cloud-free images were not available. Pixel-base temporal fitting for a LTS can be conducted using curve fitting tools [58] and temporal regression or segmentation algorithms such as LandTrendr [9,52].

3.2. Improving the Accuracy of AGB Modelling

While a pixel-based temporal fitting process can be used for filling data gaps, a more important purpose of this analysis is to minimize temporal noise/variations (due to exogenous factors such as sun angle, phenology and atmosphere) at the pixel-level of LTS [9,51,57,58,59,76,77]. It is demonstrated that this process enables the consistent implementation of a statistical model over time and improves the accuracy of biomass predictions (Table 2). The use of calibrated predictor values is the key factor enabling the consistent implementation of a statistical model across both space and time [20]. Advanced developments in image pre-processing (Section 2.1) facilitate the generation of normalized time-series images with minimal data gaps and a consistent radiometric response. This established a foundation for estimating forest AGB measurements through time [51]. However, modelling continuous forest attributes such as AGB can be negatively impacted by residual temporal noise or variations existing in Landsat pixel-level spectral trajectories [28,48,57,58]. Thus, a pixel-based temporal fitting process is necessary to mitigate temporal fluctuations and thus improve the temporal coherence in AGB predictions. Deo, Russell [58] demonstrated that using spectral data (i.e., spectral bands and indices) extracted from a fitted LTS can improve the accuracy of AGB predictions (RMSE = 47.64 Mg·ha−1) in comparison with using observed data (RMSE = 54.18 Mg·ha−1), as the temporal noise has been removed through a fitting process.
Pixel-based temporal fitting analysis captures spectral changes and trends of forested pixels (due to forest disturbance and recovery processes) over time [37,38,40,42,78]. It is widely confirmed that LTS-derived change metrics characterising the patterns of forest disturbance and recovery can improve the accuracy of AGB estimates (Table 2). An early example of this application showing the utility of disturbance and recovery metrics derived from an annual LTS analysis (1972–2010) in predicting forest AGB is presented by Pflugmacher, Cohen [28]. Spectral change metrics such as disturbance and recovery onset years, duration, magnitude, and time since disturbance (TSD), pre- and post-change conditions were derived from a temporal segmentation analysis using the LandTrendr algorithm [40]. The study compared LTS-based models with traditional modelling approaches which rely on single-date Landsat imagery and Lidar data. They found that forest disturbance and recovery metrics substantially improved the accuracy of model predicting AGB of live trees (R2 = 0.80, RMSE = 46.9 Mg·ha−1) in comparison with single-date Landsat data (R2 = 0.58, RMSE = 65.1 Mg·ha−1). Moreover, LTS-based models were significantly better in predicting AGB of dead trees (R2 = 0.73, RMSE = 31.0 Mg·ha−1) than Lidar and single-date models (R2 = 0.21, RMSE = 43.8 Mg·ha−1; R2 = 0.00, RMSE = 47.8 Mg·ha−1, respectively). The capacity of LTS-derived change metrics for estimating forest attributes, including AGB, was further quantified by Frazier, Coops [54]. In this study, a number of disturbance and recovery metrics (both simple and complex metrics) were calculated from LandTrendr-derived spectral trajectories of TC components (TCB, TCG and TCW). Eleven RF models were then tested based on different groups of predictor variables to investigate the importance of change metrics for estimating forest biomass. The study concluded that simple change metrics (such as change onset year, duration and magnitude, pre- and post-change values) have more predictive power than complex metrics (such as longest recovery, year to year change rates, and last monotonic trends). In addition to spectral change metrics, change attributions derived from LTS such as disturbance type (or causal agent) and severity level have been recently used as predictor variables for forest AGB modelling [9,49,52]. These metrics can contribute to the improvement of prediction accuracy by distinguishing pixel-level spectral conditions. For instance, pixels can experience a similar spectral change magnitude, at the same point in time, but different casual agents (e.g., fire and logging), resulting in different forest structure and biomass [47,59]. The utility of other LTS-based change metrics such as forest age has also been demonstrated in different contexts [61,62]. Liu, Peng [62] found that including LTS-derived forest age significantly improved the accuracy of AGB model when compared to using spectral indices only (R2 value increased from 0.50 to 0.73). The use of LTS change metrics in AGB estimation has been facilitated by developments in change detection algorithms (Table 1) and innovative cloud-based methods, such as the integration of the LandTrendr algorithm on Google Earth Engine [43].
As an exception, Gómez, White [63] used “dynamic” variables extracted from LTS pixel trajectories instead of spectral change metrics. “Dynamic” variables included transformed spectral trajectories and temporal derivatives associated with those trajectories. These variables were derived as identifiers of temporal patterns in spectral trajectories using two-dimensional (2D) wavelet transforms of an ordered system. The study demonstrated that “dynamic” variables provided some improvements in modelling and explaining historical AGB variability. However, the extraction of such “dynamic” variables is dependent on the availability of re-measured AGB data.
Though most applications use an annual or near-annual LTS, some studies have demonstrated the potential of using a higher temporal density LTS to improve the estimation of forest AGB [64,79]. Zhu and Liu [64] explored the use of seasonal LTS and found that seasonal NDVI time-series can improve the accuracy of AGB estimations compared with a single-date NDVI (R2 = 0.54 vs. R2 = 0.42). The study also indicated that NDVI in the fall season (September and October) has stronger correlation with AGB (R = 0.64) than in the peak season (June, R = 0.46). Nguyen, Jung [79] used LTS to estimate seasonal AGB for mixed coniferous forests in South Korea. Results from this study indicated that using Landsat images in summer (August) can provide more accurate estimation of forest AGB than using images from other seasons, although seasonal accuracies were relatively stable. The utility of using high temporal LTS (i.e., seasonal or denser scales) should be further investigated across different contexts.

4. What LTS-Based Approaches Have Been Demonstrated for Estimating AGB and Its Dynamics across Space and Time?

In this section, we review recent studies using LTS data as inputs for estimating forest AGB and its dynamics across space and time (Table 3). Of note, all studies investigated boreal and temperate forests in North America (US and Canada) and Europe, none of them was conducted in tropical forests. Study areas ranged from regional scales (several Landsat scenes) to continental scales such as over 650 million ha of Canadian forests [51]. Commonly investigated periods ranged from 20 to 40 years. Most of the studies used Landsat TM/ETM+ data to form a LTS stack (starting in the 1980s) while few of them utilized MSS data to extend the time-series further back to 1972 [12,62]. Studies often developed LTS stacks at an annual timescale using cloud-free observations or image composites. However, multi-year epochal or seasonal-scale approaches were also used in some applications. In addition, all studies estimated historical AGB, except Ma, Domke [80] who made AGB future projections.
As shown in Table 3, many LTS-based approaches for estimating and characterising forest AGB dynamics have been carried out in different contexts and they have become more advanced and robust over time. A conceptual diagram for estimating AGB dynamics from LTS is shown in Figure 1. LTS data are processed to extract explanatory variables (e.g., spectral indices and change metrics). These are then used in a modelling framework as predictor variables along with reference data (derived from forest inventory or Lidar-based plots). The implementation of this step results in predictions of AGB across space and time. Patterns of AGB dynamics are often characterized by integrating AGB changes with forest disturbance and recovery history, which are normally derived from LTS. Accuracy of predictions is often assessed using independent field plots or Lidar-based AGB data. A few studies, however, followed other approaches (e.g., Ma, Domke [80] and Zhang, Lu [84]) to directly predict AGB dynamics from LTS (Table 3).

4.1. Modelling Approaches

4.1.1. Explanatory Data

Explanatory data used for multi-temporal AGB modelling generally include three groups of variables: (1) Landsat image-based spectral variables (i.e., Landsat surface reflectance bands, vegetation indices, band ratios, and texture metrics), (2) LTS-derived spectral change/trend metrics and (3) other ancillary data such as topographic and climatic layers (Table 3).
Image-based spectral variables are key predictors for estimating forest AGB from Landsat data. These tend to incorporate the near-infrared (NIR) and shortwave infrared (SWIR) wavelengths such as TM/ETM+ band 5, OLI band 6, NBR, NDVI and EVI, given they often have strong relationships with observed AGB values [10,85]. The TC transformations, including TCB (overall reflectance), TCG (contrast between NIR and visible reflectance), and TCW (contrast between NIR and SWIR) [86,87,88,89], are also used due to their abilities in revealing key forest attributes, including AGB [28,54,85,90]. In addition to indices directly derived from Landsat bands, the use of LTS has facilitated the introduction of some pseudo/integrated spectral indices that are useful for AGB modelling, such as TCA, TCD, DI and IFZ (see Appendix A for descriptions). These indices are often used with an assumption that using an integrated index may provide more useful information than using a group of individual indices. Healey, Cohen [91] showed that using DI alone provided better results in a land cover change analysis than using a group of TC components or TM reflectance bands. Since first introduced by Powell, Cohen [81], TCA, a combination of TCB and TCG, has been utilised as an important predictor for estimating forest AGB dynamics in a wide range of research contexts (Table 3). Though IFZ (representing the likelihood of a pixel being forested) was initially used for forest change detection [37], it has been used for AGB modelling in some recent studies [57,58,80]. According to Deo, Russell [58], this metric provides a normalized predictor that can significantly reduce the spatial and temporal variability of spectral signatures caused by atmospheric conditions and sensor issues, improving the accuracy of AGB predictions. However, the utility of this index in estimating forest AGB has only been demonstrated in deciduous and coniferous forests (North America) and thus needs to be further investigated in other forested areas.
In the context of estimating forest AGB over time, the selection of spectral indices is not only dependent on their relationship with observed AGB values but also their ability to capture the forest’s history of disturbance and recovery. Among reviewed studies (Table 3), spectral indices such as NBR, NDVI, TCW and TCA are the most popular given their sensitivity to changes occurring in forests. These indices are often analysed through a temporal fitting process to extract forest disturbance and recovery metrics that contribute to the improvement of AGB predictions (Section 3). Another consideration when selecting spectral indices for estimating forest AGB over time is which sensors are being used. Studies utilising MSS data often relied on simple band ratios [62] or modelled TC indices [12]. Some spectral indices that are useful for AGB estimation and forest change detection such as NBR cannot be calculated from MSS images due to the lack of SWIR bands. This suggests that the inclusion of MSS data in LTS should be specifically considered for each case study. Spectral indices that are the most important for AGB estimation, and also the most sensitive to forest change, should be determined before processing MSS data. If selected indices are not available from MSS data, then the inclusion of MSS images can negatively impact temporal predictions of forest AGB.
Forest AGB estimates can be enhanced by incorporating spectral change metrics representing disturbance and recovery trends prior to the prediction date (Table 2, Section 3). However, the use of disturbance and recovery metrics for multiple-date AGB predictions is different from a single-date estimation. A single-date estimation can utilise temporal trends extracted from the full length of LTS. In contrast, when estimating AGB over multiple dates, the length of LTS diminishes as the prediction date moves further back in time and consequently, change metrics contain less information on disturbance and recovery which may negatively impact prediction accuracy [12]. Pflugmacher, Cohen [12] found that model accuracy gradually reduced when decreasing the length of LTS. Furthermore, they noticed that the appropriate length of LTS needed to derive meaningful change metrics is dependent on the disturbance frequency and intensity of a given study area (in their case between 10–20 years). Gómez, White [63] demonstrated that 15–25 years is sufficient for identifying significant temporal patterns in spectral trajectories to support AGB estimates in pine forests, Spain. This suggests that the balance between the prediction accuracy and the required temporal density will depend on the specific application. If only periodic dates are required for making predictions (e.g, more than 10 years), spectral disturbance and recovery metrics should be included to optimize the accuracy of biomass predictions. In contrast, if predictions are required at higher temporal densities (e.g., annual), single-date predictor variables should be used to maintain the temporal accuracy. Recent studies have incorporated some change attributions, such as disturbance causal agents and TSD, in annual modelling of forest biomass [9,51,52,62]. Though the predictive power of these metrics is normally not significant, they can benefit AGB modelling by distinguishing pixels with similar spectral information [48].
In addition to LTS-based variables, ancillary data such as topographic and climatic information are often included in AGB estimating across large areas (Table 3). The importance of these variables for predicting AGB has been widely highlighted given topographic and climatic conditions are often diverse over large areas, resulting in different forest ecosystems.

4.1.2. Training Data

For predicting AGB across space and time, LTS data need to be combined with plot training data. Studies often use two main sources of AGB training data: forest inventory and Lidar-based plots. Observed AGB estimates are often directly calculated from forest inventory plots using species-based allometric relationships between tree metrics and their biomass. On the other hand, Lidar-based AGB plots are predicted through a modelling process that combines Lidar with forest inventory data.
Forest inventory is the most common source of AGB reference data. A significant benefit of using forest inventory plots is that they provide the most accurate observations of AGB. Unfortunately, the number of available field plots is often limited across large and remote areas. When developing a LTS-based model, inventory plot data collected from multiple points in time and data sources can be incorporated, increasing the sample density for AGB modelling. For example, Kennedy, Ohmann [9] integrated all available plots from several inventories surveyed during 1991–2011 into a database of 4318 plots with 8454 measurements. The authors then developed a monolithic model to estimate annual AGB by combining the plot database with LTS data (extracting spectral values from image dates coinciding with plot measurement dates). In another study, Boisvenue, Smiley [50] used 1381 field plots across Saskatchewan, Canada (originating from 1949) and LTS data to estimate AGB dynamics from 1984 to 2012. It is important to note that uncertainties arising from ingesting different data sources need to be carefully considered given the inconsistency in inventory methods including plot sizes and survey techniques [9,50]. The use of forest inventory data is also dependent on modelling methods used for estimating AGB dynamics. Some approaches rely on data from remeasured plots which are often not available in many forested areas [50,80,84].
Airborne Lidar data provide the most accurate remote sensing-based predictions of forest AGB, but they are often not available wall-to-wall and over time given high acquisition costs and computational requirements. Alternatively, recent studies often employ airborne Lidar data as a sampling tool to improve the efficiency of forest inventory data. Sample-based Lidar plots can significantly increase the number of accurate samples of AGB estimates in areas without sufficient field inventory data [92], reducing the variability in plot-level AGB values. Lidar-based plots can be used as an alternative for training models based on LTS data, which are temporally and spatially complete [48,49].
The utility of Lidar-based plots in estimating forest AGB dynamics has been demonstrated in several studies (Table 3). Pflugmacher, Cohen [12] used airborne Lidar to sample a wide range of forest disturbance histories, as an alternative of 51 field plots within their study area. Lidar-based AGB plots were then combined with LTS data for large-area modelling of annual forest AGB over 38 years (1972–2009). Matasci, Hermosilla [51] extracted nearly 85,000 Lidar plots from 34 survey transects across Canadian forest ecosystems and integrated these plots with LTS to produce annual AGB estimates and related dynamics from 1984 to 2016 at the national scale. However, Lidar-based AGB predictions contain uncertainties that should be considered when describing the uncertainties in Landsat-based predictions [12,48]. Although Lidar data often show a strong relationship with observed AGB values, the consequences of introducing additional uncertainty into Landsat-based AGB maps need to be further investigated. Furthermore, Lidar-based plots are often spatially constrained by Lidar transects, which may not cover the full range of the forest population (e.g., transects may not be systematically distributed across different forest types). This will add more uncertainty into AGB predictions, especially in the case of a kNN imputation model, with k = 1, being used (i.e., predicted values cannot exceed the range of sample dataset) [48]. Another limitation is that Lidar-based plots cannot provide species composition information due to the lack of tree measurements; therefore, they cannot be used when deriving species specific AGB estimates.

4.1.3. Statistical Modelling Technique

Various techniques have been deployed for modelling AGB across space and time. These can be grouped into three broad model types: parametric regression, nonparametric regression and nonparametric imputation. Parametric regression is a traditional method for estimating AGB based on remote sensing, which finds the linear relationship between observed AGB and predictor variables. Few studies opt for this method (Table 3), given it is known to overestimate and underestimate low and high AGB values, respectively. Furthermore, Powell, Cohen [81] and Zhang, Lu [84] compared different empirical approaches for modelling AGB dynamics from LTS and both found that regression models, such as Linear Regression and Reduced Major Axis, provide less accurate results than non-parametric methods such as Random Forest (RF) and Gradient Boost Regression Tree.
Non-parametric approaches are increasingly used to derive forest structural attributes from remote sensing data as they do not have restrictive assumptions on the data analysed. Two nonparametric regression approaches commonly used for modelling AGB dynamics are regression trees and kNN [93] (Table 3). Regression tree methods are widely used for forest applications due to their tremendous analytical and operational flexibility. RF is a regression tree algorithm widely used for forest modelling. It is an ensemble learning algorithm that develops a number of small regression trees that vote on predictions. Each tree is constructed from a random subsample from the original dataset and nodes on trees are split using a random subset of predictor variables [94]. RF is capable of accounting for complex and non-linear relationships among variables, and robust in preventing over-fitting. An advantage of using RF is that it can handle a large number of variables/metrics extracted from LTS data. In addition to RF, other regression tree algorithms such as Decision tree (known as CART) and Gradient Boost Regression Tree have also been used for modelling AGB dynamics [63,84]. kNN regression predicts the AGB value for a target sample (pixel) as the mean of the k most similar reference samples, defined using a distance metric. Although it is a commonly used modelling method for forest applications [95], few studies have used kNN regression for estimating AGB over time based on LTS [79]. Nevertheless, Wilson, Knight [55] demonstrated that kNN regression and RF produced comparable accuracies when modelling forest AGB from LTS data.
Non-parametric imputation approaches have recently received considerable attention for forest inventory mapping applications. These multivariate modelling approaches can simultaneously map large assemblages of inventory attributes, therefore retaining the variance structure to that of the observations [96]. For predicting forest AGB dynamics from LTS, the kNN imputation has been frequently employed (Table 3). This method imputes (or shares) observed values of one or multiple response variables to target samples/pixels [8]. Kennedy, Ohmann [9] demonstrated that the kNN imputation approach performed well to overcome the saturation of Landsat signals known to occur in forests with high biomass. Unlike regression, where a predicted value is often different from any observed values, an imputed value will be the most similar observed value when the single nearest neighbour is used (k = 1). When k > 1, an imputed value will be the mean of observed values associated to the k nearest neighbours. While using a k > 1 can improve the imputation accuracy and produce more smoothed predictions, the use of k = 1 is preferable to maintain the variance structure of predictions that are similar to the observations [76]. Most studies deploying kNN for AGB dynamic mapping used the single nearest neighbour as they preferred to utilize the full range of AGB values in training data [9,51,52]. Deo, Russell [58] suggested that a k value of 3 to 5 is the most suitable in order to improve model accuracy while retaining the variance structure of predictions closer to observations. Nguyen, Jones [52] achieved this using k = 1 when processing the kNN model though a bootstrapping analysis with multiple repetitions.
In kNN imputation, the nearest neighbours are found based on a distance metric derived using different techniques such as Gradient Nearest Neighbour (GNN) [68] and RF [97]. GNN is a weighted Euclidean distance-based technique which is a combination of the canonical correspondence analysis and k = 1. In contrast, RF is a machine learning-based technique in which the distance metric is calculated based on a proximity matrix [98]. Several studies have shown the utility of GNN in modelling AGB dynamics [9,81]. In a recent study, Kennedy, Ohmann [9] integrated LandTrendr and GNN-based imputation in a robust empirical forest biomass monitoring system. RF-based kNN imputation has been increasingly used in various contexts for estimating AGB dynamics based on LTS data. For example, Matasci, Hermosilla [51] and Nguyen, Jones [52] used this method for large area predictions of AGB dynamics through 30-year periods in Canada and southeast Australia, respectively. Studies comparing different kNN distance techniques used for imputing forest attributes often indicated the outperformance of RF [8,56,96].
When using kNN imputation, AGB can be predicted using either a direct or indirect imputation approach. In the direct approach, AGB is imputed by directly linking observed AGB (as a response variable) with Landsat-derived predictors through a kNN model [58,81]. However, recent studies often used an indirect imputation approach for estimating AGB. In this method, a kNN model is trained by forest structural variables extracted from inventory plots [9,52] or Lidar metrics extracted from Lidar plots [51]. AGB is not included in the model and is attached as an ancillary variable to the imputation process for target samples. In other words, the nearest neighbours are found according to the relationship between structural variables and LTS-derived predictors and then indirectly transferred to AGB. Response variables included in the model should appropriately capture forest conditions and can be explained by Landsat spectral data. The accuracy of AGB predictions can be enhanced as it is imputed based on another variable (or group of variables) that has higher correlation with predictors. Nguyen, Jones [56] found that AGB can be better predicted using an indirect imputation approach combining a group of basal area and stem density variables with LTS. In addition, other forest attributes, along with AGB, can be mapped directly without new model developments. Kennedy, Ohmann [9] combined field plot measurements, including basal area by species and tree size classes, with LTS-derived predictors through a GNN-based imputation model. This resulted in yearly inventory-like maps (each Landsat pixel was assigned tree measurements from the nearest inventory plot), from which forest AGB and other attributes were calculated. Matasci, Hermosilla [51] developed a RF-based imputation model trained by six Lidar measurements of return height and percentage. Using this model, forest structural variables, including AGB, were then indirectly imputed across space and time.
Ma, Domke [80] introduced a novel approach using matrix models for estimating AGB dynamics from LTS. The study first developed a matrix model based on remeasured inventory data to predict AGB from forest structural variables including stand basal area. The basal area in the matrix was then replaced by Landsat spectral variables to project large-scale AGB dynamics short (5 years) and long-term (30 years). This modelling method incorporates various combinations of remote sensing and inventory data, avoiding a loss in accuracy from only using inventory variables. However, this approach is dependent on the availability of cyclical inventory data.

4.2. Accuracy Assessment

Depending on the availability of validation data, different approaches can be used to assess the performance of AGB models (Table 3). Where the sample size is large enough, studies often split the reference data into two sets (training and testing data), so the model can be cross-validated using the testing data. However, this approach often increases the variance in model training data, causing negative impacts on model accuracy. Therefore, where the sample size of reference data is small, studies often use internal assessment methods such as RF out-of-bag (OOB) errors and n-fold cross-validation (also known as leave-one-out cross-validation when n = 1). Some studies compare AGB perdition maps with independent AGB products such as national biomass maps [58,82] or Lidar-based AGB data [9,49,52]. Lidar-based AGB data are often used as they are widely confirmed as the most accurate remote sensing-based predictions.
Perhaps more important than model accuracy, when extending single model through time, is the need to evaluate (1) the temporal transferability of the model and (2) the robustness of the model in capturing AGB changes resulting from forest disturbance and recovery. Unfortunately, this is often dependent on the availability of validation data in multi-points in time (i.e., cycle inventory or Lidar data), which are normally not available in many forested regions. Few studies that estimated historical AGB have examined the temporal transferability of the model extension (Table 3). Powell, Cohen [81] performed a Landsat scene-level assessment by comparing scene-level predicted AGB trajectories with temporal trends of observed AGB derived from remeasured inventory plots. Matasci, Hermosilla [51] conducted a multi-temporal assessment of annual AGB predictions using Lidar data acquired from 2006 to 2012. In a recent study, Nguyen, Jones [52] conducted a time-series assessment of 30-year annual AGB prediction maps using un-changed, forested Lidar-based plots. Interestingly, these studies all confirmed that a single-date LTS-based model is transferable to estimate AGB over time.
Validating of AGB change is the most challenging as it requires repeated AGB observations. Where remeasured inventory plots are available, studies often extracted and compared observed AGB change (i.e., the difference of AGB between measurement rounds) with predicted AGB change [12,50,63,80,81,84]. Powell, Cohen [81] indicated that the accuracy of AGB change prediction improved when applying a pixel-based temporal fitting to AGB trajectories. Boisvenue, Smiley [50] fit a linear mixed-effect model to both observed and predicted AGB change datasets to estimate the long-term mean values of AGB change for each forest type ( Δ AGB   =   f ( a g e ) ), which were then compared together. Where remeasured inventory plots are unavailable, multi-temporal Lidar data can be used as an alternative. Nguyen, Jones [52] did this when evaluating the ability of a LTS-based model in capturing AGB change resulting from forest disturbance and recovery processes. If repeated AGB observations are not available for validation, comparing with results from other studies could be an alternative. For example, Powell, Cohen [82] compared LTS-based AGB loss from disturbance (at a stratum-level) against results from several studies conducted within the same study area.

4.3. Characterising Spatial and Temporal Patterns of Forest AGB Dynamics

Forest AGB maps are normally created at the spatial resolution of 30 m, aligning with Landsat pixels, and at an annual or near-annual temporal density (few studies produced seasonal or multi-year epochal maps, Table 3). AGB change analysis across space and time is a crucial component of a forest biomass monitoring system. Along with developing robust modelling methods for estimating historical AGB, studies have proposed various approaches for quantifying spatial and temporal patterns of AGB dynamics to effectively support forest management (Table 3).
Methods for characterising patterns of AGB dynamics can be stratified into two broad groups: difference-based and trajectory-based. Difference is a traditional method for calculating AGB change by subtracting AGB values between two prediction dates. Though this is a simple and straightforward method, it often provides only generic information on AGB change. This method has been normally used when the temporal density of AGB predictions is low (often multiple-year epochal) and the intervals between prediction dates are wide enough to capture actual AGB changes in forests (5–10 years) [58,63,83,84]. As an exception, Boisvenue, Smiley [50] used the difference-based method to compute AGB changes from 29 annual AGB maps, but the authors then developed a mixed-effect model to summarize long-term trends of AGB at a stratum-level.
In contrast, most studies predicting AGB at an annual scale often used a trajectory-based method for characterizing AGB dynamics (Table 3). This method allows consistently and comprehensively tracking and understanding AGB changes across space and time according to forest disturbance and recovery processes. Some studies applied a temporal fitting process to annual AGB trajectories to create forest AGB disturbance and recovery maps [72,81]. Others integrated fitted AGB trajectories with LTS-derived forest disturbance maps to quantify AGB loss from disturbance [82] and according to different disturbance causal agents [9]. Matasci, Hermosilla [51] combined annual AGB predictions with an LTS-derived disturbance map to characterise AGB dynamics of Canadian forests across three decades. Temporal AGB trends for each ecozone were summarised for three scenarios: undisturbed forests, forests impacted by wildfires and harvesting. Nguyen, Jones [52] characterised spatial and temporal patterns of AGB dynamics by computing change metrics, including AGB loss and gain from disturbance and recovery, Recovery Indicator (RI), and Years to Recovery (Y2R). These allowed authors to robustly track AGB change according to different disturbance scenarios (i.e., disturbance severity and causal agents).

5. Conclusions and Future Opportunities

The opening to the public of the Landsat archive since 2008 has facilitated the use of LTS for forest applications. In this review, we provide an overview of recent studies utilizing LTS for estimating forest AGB and its dynamics. In summary, the use of LTS has the potential to improve the estimates of forest AGB in comparison with using single-date images. It provides a unique opportunity to create spatially and temporally complete products of forest AGB across large areas and long time-series. Furthermore, time-series processing and variables extracted from LTS can help improve the consistency and accuracy of AGB modelling. Many innovative approaches for estimating forest AGB across space and time have been demonstrated throughout the last decade, with methods becoming more advanced and robust in more recent studies. Thus, LTS-based products of AGB dynamics have become more reliable to support forest managers and researchers.
Estimating forest AGB dynamics across space and time is an ongoing topic of interest. LTS is increasingly used for this purpose as it provides the longest (and free) collection of imagery. Most current studies focused on boreal and temperate forests, highlighting the need for exploring the utility of LTS for estimating AGB dynamics in tropical forests. In contrast to boreal and temperate forests, using LTS to calculate AGB dynamics is all the more challenging in tropical forest regions due to persistent cloud cover. In addition, tropical forests are mainly found in developing countries, which often lack a systematic inventory network and Lidar data. Implementing a single model to predict AGB over time can produce temporal uncertainties (e.g., the balance of the model over time), which has mostly been ignored in current studies. More research investigating both spatial and temporal accuracies of LTS-derived AGB change are anticipated over coming years. Along with Landsat, other satellites such as Sentinel 2 and SPOT have been providing imagery as a time-series, with higher spatial and temporal resolutions. Incorporating such remote sensing data with LTS for improving the estimates of forest AGB dynamics will be a main area of research in the near future. Some studies are also investigating opportunities to combine LTS with active remote sensing data such as synthetic aperture radar (SAR) and the Global Ecosystem Dynamics Investigation (GEDI) for forest biomass estimation [99]. This is especially promising for tropical areas that suffer from extensive cloud/shadow issues. Last but not least, some scientists are attempting to project future forest AGB dynamics using LTS [80]. This is a challenging task for the forest remote sensing community, expected to be one of the main foci in the next few years. Trends of AGB over the last 40 years derived from LTS will be fundamental for predicting future AGB dynamics.

Author Contributions

Conceptualization, T.H.N., S.J., M.S.-B. and A.H.; formal analysis, T.H.N. and S.H.; supervision, S.J., A.H. and M.S.-B.; writing—original draft, T.H.N.; writing—review and editing, S.J., M.S.-B., S.H. and A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This research was partly funded by the Australian Award Scholarship (AAS) and the Cooperative Research Centre for Spatial Information (CRCSI) under Project 4.104 (A Monitoring and Forecasting Framework for the Sustainable Management of southeast Australian Forests at the Large Area Scale). CRCSI activities are funded by the Australian Commonwealth’s Cooperative Research Centres Programme.

Conflicts of Interest

The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; or the writing of the manuscript.

Appendix A

Table
Table A1. Landsat spectral indices commonly used for forest AGB estimates.
Table A1. Landsat spectral indices commonly used for forest AGB estimates.
Landsat Spectral IndexCalculation
Normalized Difference Vegetation Index (NDVI) [100]NDVI = (NIR − R)/(NIR + R)
Normalized Burn Ratio (NBR) [101]NBR = (NIR − SWIR)/(NIR + SWIR)
Normalized Difference Moisture Index (NDMI)NDMI = (NIR − SWIR)/(NIR + SWIR)
Enhanced Vegetation Index (EVI) [102]EVI = G * ((NIR − R)/(NIR + C1 * R – C2 * B + L))
L = value to adjust for canopy background, C = coefficients for atmospheric resistance, B = the blue band
Soil Adjusted Vegetation Index (SAVI) [103]SAVI = ((NIR − R)/(NIR + R + L)) * (1 + L)
Chlorophyll Vegetation Index (CVI) [104]CVI = (NIR x R)/G
G = the green band
Difference Vegetation Index (DVI) [105]DVI = NIR − R
Linear transform of multiple bands [84]VIS123 = B + G + R
MID57 = TM band 5 + TM band 7 (SWIR)
Integrated Forest Z-score (IFZ) [37]z-score measure of a pixel likelihood of being forested, using TM bands 3, 5 and 7
Tasseled Cap (TC) transformations:
TC brightness (TCB); TC greenness (TCG); TC wetness (TCW) [86,87,88,89]		TCW, TCB, and TCG are calculated by multiplying Landsat band pixel values with TC coefficients. See the coefficients in references.
TC angle (TCA) [81]TCA = arctan(TCG/TCB)
TC distance (TCD) [106]   TCD   = TCB 2 + TCG 2    
TC Disturbance Index (DI) [91]DI = TCBr – (TCGr + TCWr)
r = denotes rescaled TC indices based upon the mean and standard deviation of the scene’s forest values

References

  1. U.S. Geological Survey. Landsat—Earth Observation Satellites. In Fact Sheet; U.S. Geological Survey: Reston, VA, USA, 2015; p. 4. [Google Scholar]
  2. U.S. Geological Survey. Landsat 9. In Fact Sheet; U.S. Geological Survey: Reston, VA, USA, 2019; p. 2. [Google Scholar]
  3. Woodcock, C.E.; Allen, R.; Anderson, M.; Belward, A.; Bindschadler, R.; Cohen, W.; Gao, F.; Goward, S.N.; Helder, D.; Helmer, E. Free access to Landsat imagery. Science 2008, 320, 1011. [Google Scholar] [CrossRef]
  4. Wulder, M.A.; Masek, J.G.; Cohen, W.B.; Loveland, T.R.; Woodcock, C.E. Opening the archive: How free data has enabled the science and monitoring promise of Landsat. Remote Sens. Environ. 2012, 122, 2–10. [Google Scholar] [CrossRef]
  5. Soenen, S.A.; Peddle, D.R.; Hall, R.J.; Coburn, C.A.; Hall, F.G. Estimating aboveground forest biomass from canopy reflectance model inversion in mountainous terrain. Remote Sens. Environ. 2010, 114, 1325–1337. [Google Scholar] [CrossRef]
  6. Wulder, M.; Skakun, R.; Kurz, W.; White, J. Estimating time since forest harvest using segmented Landsat ETM+ imagery. Remote Sens. Environ. 2004, 93, 179–187. [Google Scholar] [CrossRef]
  7. Hudak, A.T.; Strand, E.K.; Vierling, L.A.; Byrne, J.C.; Eitel, J.U.; Martinuzzi, S.; Falkowski, M.J. Quantifying aboveground forest carbon pools and fluxes from repeat LiDAR surveys. Remote Sens. Environ. 2012, 123, 25–40. [Google Scholar] [CrossRef] [Green Version]
  8. Hudak, A.T.; Crookston, N.L.; Evans, J.S.; Hall, D.E.; Falkowski, M.J. Nearest neighbor imputation of species-level, plot-scale forest structure attributes from LiDAR data. Remote Sens. Environ. 2008, 112, 2232–2245. [Google Scholar] [CrossRef] [Green Version]
  9. Kennedy, R.E.; Ohmann, J.; Gregory, M.; Roberts, H.; Yang, Z.; Bell, D.M.; Kane, V.; Hughes, M.J.; Cohen, W.B.; Powell, S.; et al. An empirical, integrated forest biomass monitoring system. Environ. Res. Lett. 2018, 13, 025004. [Google Scholar] [CrossRef]
  10. Lu, D. The potential and challenge of remote sensing—Based biomass estimation. Int. J. Remote Sens. 2006, 27, 1297–1328. [Google Scholar] [CrossRef]
  11. Franklin, J. Thematic Mapper analysis of coniferous forest structure and composition. Int. J. Remote Sens. 1986, 7, 1287–1301. [Google Scholar] [CrossRef]
  12. Pflugmacher, D.; Cohen, W.B.; Kennedy, R.E.; Yang, Z. Using Landsat-derived disturbance and recovery history and lidar to map forest biomass dynamics. Remote Sens. Environ. 2014, 151, 124–137. [Google Scholar] [CrossRef]
  13. Zhu, Z. Change detection using landsat time series: A review of frequencies, preprocessing, algorithms, and applications. ISPRS J. Photogramm. Remote Sens. 2017, 130, 370–384. [Google Scholar] [CrossRef]
  14. Banskota, A.; Kayastha, N.; Falkowski, M.J.; Wulder, M.A.; Froese, R.E.; White, J.C. Forest Monitoring Using Landsat Time Series Data: A Review. Can. J. Remote Sens. 2014, 40, 362–384. [Google Scholar] [CrossRef]
  15. Hansen, M.C.; Loveland, T.R. A review of large area monitoring of land cover change using Landsat data. Remote Sens. Environ. 2012, 122, 66–74. [Google Scholar] [CrossRef]
  16. Zhu, Z.; Woodcock, C.E. Continuous change detection and classification of land cover using all available Landsat data. Remote Sens. Environ. 2014, 144, 152–171. [Google Scholar] [CrossRef] [Green Version]
  17. USGS. Landsat Collection1 Level 1 Product Definition; United State Geological Survey: Reston, VA, USA, 2019.
  18. Gutman, G.; Huang, C.; Chander, G.; Noojipady, P.; Masek, J.G. Assessment of the NASA–USGS global land survey (GLS) datasets. Remote Sens. Environ. 2013, 134, 249–265. [Google Scholar] [CrossRef]
  19. Vogeler, J.C.; Braaten, J.D.; Slesak, R.A.; Falkowski, M.J. Extracting the full value of the Landsat archive: Inter-sensor harmonization for the mapping of Minnesota forest canopy cover (1973–2015). Remote Sens. Environ. 2018, 209, 363–374. [Google Scholar] [CrossRef]
  20. Song, C.; Woodcock, C.E.; Seto, K.C.; Lenney, M.P.; Macomber, S.A. Classification and change detection using Landsat TM data: When and how to correct atmospheric effects? Remote Sens. Environ. 2001, 75, 230–244. [Google Scholar] [CrossRef]
  21. Song, C.; Woodcock, C.E. Monitoring forest succession with multitemporal Landsat images: Factors of uncertainty. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2557–2567. [Google Scholar] [CrossRef]
  22. Vermote, E.F.; Tanre, D.; Deuze, J.L.; Herman, M.; Morcette, J. Second Simulation of the Satellite Signal in the Solar Spectrum, 6S: An overview. IEEE Trans. Geosci. Remote Sens. 1997, 35, 675–686. [Google Scholar] [CrossRef] [Green Version]
  23. Masek, J.G.; Vermote, E.F.; Saleous, N.E.; Wolfe, R.; Hall, F.G.; Huemmrich, K.F.; Gao, F.; Kutler, J.; Lim, T.-K. A Landsat surface reflectance dataset for North America, 1990–2000. IEEE Geosci. Remote Sens. Lett. 2006, 3, 68–72. [Google Scholar] [CrossRef]
  24. Vermote, E.; Justice, C.; Claverie, M.; Franch, B. Preliminary analysis of the performance of the Landsat 8/OLI land surface reflectance product. Remote Sens. Environ. 2016, 185, 46–56. [Google Scholar] [CrossRef]
  25. Zhu, Z.; Woodcock, C.E. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sens. Environ. 2012, 118, 83–94. [Google Scholar] [CrossRef]
  26. Chavez, P.S., Jr. An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data. Remote Sens. Environ. 1988, 24, 459–479. [Google Scholar] [CrossRef]
  27. Braaten, J.D.; Cohen, W.B.; Yang, Z. Automated cloud and cloud shadow identification in Landsat MSS imagery for temperate ecosystems. Remote Sens. Environ. 2015, 169, 128–138. [Google Scholar] [CrossRef] [Green Version]
  28. Pflugmacher, D.; Cohen, W.B.; Kennedy, R.E. Using Landsat-derived disturbance history (1972–2010) to predict current forest structure. Remote Sens. Environ. 2012, 122, 146–165. [Google Scholar] [CrossRef]
  29. White, J.; Wulder, M.; Hobart, G.; Luther, J.; Hermosilla, T.; Griffiths, P.; Coops, N.; Hall, R.; Hostert, P.; Dyk, A. Pixel-based image compositing for large-area dense time series applications and science. Can. J. Remote Sens. 2014, 40, 192–212. [Google Scholar] [CrossRef] [Green Version]
  30. Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W. Regional detection, characterization, and attribution of annual forest change from 1984 to 2012 using Landsat-derived time-series metrics. Remote Sens. Environ. 2015, 170, 121–132. [Google Scholar] [CrossRef]
  31. Zhu, Z.; Woodcock, C.E.; Holden, C.; Yang, Z. Generating synthetic Landsat images based on all available Landsat data: Predicting Landsat surface reflectance at any given time. Remote Sens. Environ. 2015, 162, 67–83. [Google Scholar] [CrossRef]
  32. Roy, D.P.; Ju, J.; Kline, K.; Scaramuzza, P.L.; Kovalskyy, V.; Hansen, M.; Loveland, T.R.; Vermote, E.; Zhang, C. Web-enabled Landsat Data (WELD): Landsat ETM+ composited mosaics of the conterminous United States. Remote Sens. Environ. 2010, 114, 35–49. [Google Scholar] [CrossRef]
  33. Potapov, P.; Turubanova, S.; Hansen, M.C. Regional-scale boreal forest cover and change mapping using Landsat data composites for European Russia. Remote Sens. Environ. 2011, 115, 548–561. [Google Scholar] [CrossRef]
  34. Flood, N. Seasonal Composite Landsat TM/ETM+ Images Using the Medoid (a Multi-Dimensional Median). Remote Sens. 2013, 5, 6481–6500. [Google Scholar] [CrossRef] [Green Version]
  35. Pickell, P.D.; Hermosilla, T.; Frazier, R.J.; Coops, N.C.; Wulder, M.A. Forest recovery trends derived from Landsat time series for North American boreal forests. Int. J. Remote Sens. 2015, 37, 138–149. [Google Scholar] [CrossRef]
  36. Bartels, S.F.; Chen, H.Y.H.; Wulder, M.A.; White, J.C. Trends in post-disturbance recovery rates of Canada’s forests following wildfire and harvest. For. Ecol. Manag. 2016, 361, 194–207. [Google Scholar] [CrossRef] [Green Version]
  37. Huang, C.; Goward, S.N.; Masek, J.G.; Thomas, N.; Zhu, Z.; Vogelmann, J.E. An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks. Remote Sens. Environ. 2010, 114, 183–198. [Google Scholar] [CrossRef]
  38. Zhu, Z.; Woodcock, C.E.; Olofsson, P. Continuous monitoring of forest disturbance using all available Landsat imagery. Remote Sens. Environ. 2012, 122, 75–91. [Google Scholar] [CrossRef]
  39. Devries, B.; Verbesselt, J.; Kooistra, L.; Herold, M. Robust monitoring of small-scale forest disturbances in a tropical montane forest using Landsat time series. Remote Sens. Environ. 2015, 161, 107–121. [Google Scholar] [CrossRef]
  40. Kennedy, R.E.; Yang, Z.; Cohen, W.B. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr—Temporal segmentation algorithms. Remote Sens. Environ. 2010, 114, 2897–2910. [Google Scholar] [CrossRef]
  41. Cohen, W.B.; Yang, Z.; Healey, S.P.; Kennedy, R.E.; Gorelick, N. A LandTrendr multispectral ensemble for forest disturbance detection. Remote Sens. Environ. 2018, 205, 131–140. [Google Scholar] [CrossRef]
  42. Kennedy, R.E.; Yang, Z.; Cohen, W.B.; Pfaff, E.; Braaten, J.; Nelson, P. Spatial and temporal patterns of forest disturbance and regrowth within the area of the Northwest Forest Plan. Remote Sens. Environ. 2012, 122, 117–133. [Google Scholar] [CrossRef]
  43. Kennedy, R.; Yang, Z.; Gorelick, N.; Braaten, J.; Cavalcante, L.; Cohen, W.; Healey, S. Implementation of the LandTrendr Algorithm on Google Earth Engine. Remote Sens. 2018, 10, 691. [Google Scholar] [CrossRef] [Green Version]
  44. Cohen, W.; Healey, S.; Yang, Z.; Stehman, S.; Brewer, C.; Brooks, E.; Gorelick, N.; Huang, C.; Hughes, M.; Kennedy, R.; et al. How Similar Are Forest Disturbance Maps Derived from Different Landsat Time Series Algorithms? Forests 2017, 8, 98. [Google Scholar] [CrossRef]
  45. Schroeder, T.A.; Schleeweis, K.G.; Moisen, G.G.; Toney, C.; Cohen, W.B.; Freeman, E.A.; Yang, Z.; Huang, C. Testing a Landsat-based approach for mapping disturbance causality in U.S. forests. Remote Sens. Environ. 2017, 195, 230–243. [Google Scholar] [CrossRef] [Green Version]
  46. Kennedy, R.E.; Yang, Z.; Braaten, J.; Copass, C.; Antonova, N.; Jordan, C.; Nelson, P. Attribution of disturbance change agent from Landsat time-series in support of habitat monitoring in the Puget Sound region, USA. Remote Sens. Environ. 2015, 166, 271–285. [Google Scholar] [CrossRef]
  47. Nguyen, T.H.; Jones, S.D.; Soto-Berelov, M.; Haywood, A.; Hislop, S. A spatial and temporal analysis of forest dynamics using Landsat time-series. Remote Sens. Environ. 2018, 217, 461–475. [Google Scholar] [CrossRef]
  48. Zald, H.S.J.; Wulder, M.A.; White, J.C.; Hilker, T.; Hermosilla, T.; Hobart, G.W.; Coops, N.C. Integrating Landsat pixel composites and change metrics with lidar plots to predictively map forest structure and aboveground biomass in Saskatchewan, Canada. Remote Sens. Environ. 2016, 176, 188–201. [Google Scholar] [CrossRef] [Green Version]
  49. Matasci, G.; Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W.; Zald, H.S.J. Large-area mapping of Canadian boreal forest cover, height, biomass and other structural attributes using Landsat composites and lidar plots. Remote Sens. Environ. 2018, 209, 90–106. [Google Scholar] [CrossRef]
  50. Boisvenue, C.; Smiley, B.P.; White, J.C.; Kurz, W.A.; Wulder, M.A. Integration of Landsat time series and field plots for forest productivity estimates in decision support models. For. Ecol. Manag. 2016, 376, 284–297. [Google Scholar] [CrossRef]
  51. Matasci, G.; Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W.; Bolton, D.K.; Tompalski, P.; Bater, C.W. Three decades of forest structural dynamics over Canada’s forested ecosystems using Landsat time-series and lidar plots. Remote Sens. Environ. 2018, 216, 697–714. [Google Scholar] [CrossRef]
  52. Nguyen, T.H.; Jones, S.D.; Soto-Berelov, M.; Haywood, A.; Hislop, S. Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data. Int. J. Appl. Earth Obs. Geoinf. 2020, 84. [Google Scholar] [CrossRef]
  53. Morel, A.C.; Fisher, J.B.; Malhi, Y. Evaluating the potential to monitor aboveground biomass in forest and oil palm in Sabah, Malaysia, for 2000–2008 with Landsat ETM+ and ALOS-PALSAR. Int. J. Remote Sens. 2012, 33, 3614–3639. [Google Scholar] [CrossRef]
  54. Frazier, R.J.; Coops, N.C.; Wulder, M.A.; Kennedy, R. Characterization of aboveground biomass in an unmanaged boreal forest using Landsat temporal segmentation metrics. ISPRS J. Photogramm. Remote Sens. 2014, 92, 137–146. [Google Scholar] [CrossRef]
  55. Wilson, B.T.; Knight, J.F.; McRoberts, R.E. Harmonic regression of Landsat time series for modeling attributes from national forest inventory data. ISPRS J. Photogramm. Remote Sens. 2018, 137, 29–46. [Google Scholar] [CrossRef]
  56. Nguyen, T.; Jones, S.; Soto-Berelov, M.; Haywood, A.; Hislop, S. A Comparison of Imputation Approaches for Estimating Forest Biomass Using Landsat Time-Series and Inventory Data. Remote Sens. 2018, 10, 1825. [Google Scholar] [CrossRef] [Green Version]
  57. Deo, R.; Russell, M.; Domke, G.; Andersen, H.-E.; Cohen, W.; Woodall, C. Evaluating Site-Specific and Generic Spatial Models of Aboveground Forest Biomass Based on Landsat Time-Series and LiDAR Strip Samples in the Eastern USA. Remote Sens. 2017, 9, 598. [Google Scholar] [CrossRef] [Green Version]
  58. Deo, R.K.; Russell, M.B.; Domke, G.M.; Woodall, C.W.; Falkowski, M.J.; Cohen, W.B. Using Landsat Time-Series and LiDAR to Inform Aboveground Forest Biomass Baselines in Northern Minnesota, USA. Can. J. Remote Sens. 2017, 43, 28–47. [Google Scholar] [CrossRef]
  59. Zald, H.S.J.; Ohmann, J.L.; Roberts, H.M.; Gregory, M.J.; Henderson, E.B.; McGaughey, R.J.; Braaten, J. Influence of lidar, Landsat imagery, disturbance history, plot location accuracy, and plot size on accuracy of imputation maps of forest composition and structure. Remote Sens. Environ. 2014, 143, 26–38. [Google Scholar] [CrossRef]
  60. Bolton, D.K.; White, J.C.; Wulder, M.A.; Coops, N.C.; Hermosilla, T.; Yuan, X. Updating stand-level forest inventories using airborne laser scanning and Landsat time series data. Int. J. Appl. Earth Obs. Geoinf. 2018, 66, 174–183. [Google Scholar] [CrossRef]
  61. Lefsky, M.A. 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. [Google Scholar] [CrossRef]
  62. Liu, L.; Peng, D.; Wang, Z.; Hu, Y. Improving artificial forest biomass estimates using afforestation age information from time series Landsat stacks. Environ. Monit. Assess. 2014, 186, 7293–7306. [Google Scholar] [CrossRef]
  63. Gómez, C.; White, J.C.; Wulder, M.A.; Alejandro, P. Historical forest biomass dynamics modelled with Landsat spectral trajectories. ISPRS J. Photogramm. Remote Sens. 2014, 93, 14–28. [Google Scholar] [CrossRef] [Green Version]
  64. Zhu, X.; Liu, D. Improving forest aboveground biomass estimation using seasonal Landsat NDVI time-series. ISPRS J. Photogramm. Remote Sens. 2015, 102, 222–231. [Google Scholar] [CrossRef]
  65. Hall, R.J.; Skakun, R.S.; Arsenault, E.J.; Case, B.S. Modeling forest stand structure attributes using Landsat ETM+ data: Application to mapping of aboveground biomass and stand volume. For. Ecol. Manag. 2006, 225, 378–390. [Google Scholar] [CrossRef]
  66. Zheng, D.; Rademacher, J.; Chen, J.; Crow, T.; Bresee, M.; Le Moine, J.; Ryu, S.-R. Estimating aboveground biomass using Landsat 7 ETM+ data across a managed landscape in northern Wisconsin, USA. Remote Sens. Environ. 2004, 93, 402–411. [Google Scholar] [CrossRef]
  67. Dong, J.R.; Kaufmann, R.K.; Myneni, R.B.; Tucker, C.J.; Kauppi, P.E.; Liski, J.; Buermann, W.; Alexeyev, V.; Hughes, M.K. Remote sensing estimates of boreal and temperate forest woody biomass: Carbon pools, sources, and sinks. Remote Sens. Environ. 2003, 84, 393–410. [Google Scholar] [CrossRef] [Green Version]
  68. Ohmann, J.L.; Gregory, M.J. Predictive mapping of forest composition and structure with direct gradient analysis and nearest- neighbor imputation in coastal Oregon, U.S.A. Can. J. For. Res. 2002, 32, 725–741. [Google Scholar] [CrossRef]
  69. Brian, T.; Kim, W.; Jürgen, B.; Geoff, C.; Cris, B.; Peter, K. Woody Biomass Methods for Estimating Changes; The Australian National University: Canberra, Australia, 1999. [Google Scholar]
  70. Wulder, M.; White, J.; Fournier, R.; Luther, J.; Magnussen, S. Spatially Explicit Large Area Biomass Estimation: Three Approaches Using Forest Inventory and Remotely Sensed Imagery in a GIS. Sensors 2008, 8, 529–560. [Google Scholar] [CrossRef] [Green Version]
  71. Zheng, G.; Chen, J.M.; Tian, Q.J.; Ju, W.M.; Xia, X.Q. Combining remote sensing imagery and forest age inventory for biomass mapping. J. Environ. Manag. 2007, 85, 616–623. [Google Scholar] [CrossRef]
  72. Main-Knorn, M.; Cohen, W.B.; Kennedy, R.E.; Grodzki, W.; Pflugmacher, D.; Griffiths, P.; Hostert, P. Monitoring coniferous forest biomass change using a Landsat trajectory-based approach. Remote Sens. Environ. 2013, 139, 277–290. [Google Scholar] [CrossRef]
  73. Main-Knorn, M.; Moisen, G.G.; Healey, S.P.; Keeton, W.S.; Freeman, E.A.; Hostert, P. Evaluating the Remote Sensing and Inventory-Based Estimation of Biomass in the Western Carpathians. Remote Sens. 2011, 3, 1427–1446. [Google Scholar] [CrossRef] [Green Version]
  74. Jiménez, E.; Vega, J.A.; Fernández-Alonso, J.M.; Vega-Nieva, D.; Ortiz, L.; López-Serrano, P.M.; López-Sánchez, C.A. Estimation of aboveground forest biomass in Galicia (NW Spain) by the combined use of LiDAR, LANDSAT ETM+ and National Forest Inventory data. iForest-Biogeosci. For. 2017, 10, 590–596. [Google Scholar] [CrossRef]
  75. Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W. An integrated Landsat time series protocol for change detection and generation of annual gap-free surface reflectance composites. Remote Sens. Environ. 2015, 158, 220–234. [Google Scholar] [CrossRef]
  76. Ohmann, J.L.; Gregory, M.J.; Roberts, H.M. Scale considerations for integrating forest inventory plot data and satellite image data for regional forest mapping. Remote Sens. Environ. 2014, 151, 3–15. [Google Scholar] [CrossRef]
  77. Bell, D.M.; Gregory, M.J.; Ohmann, J.L. Imputed forest structure uncertainty varies across elevational and longitudinal gradients in the western Cascade Mountains, Oregon, USA. For. Ecol. Manag. 2015, 358, 154–164. [Google Scholar] [CrossRef]
  78. DeVries, B.; Decuyper, M.; Verbesselt, J.; Zeileis, A.; Herold, M.; Joseph, S. Tracking disturbance-regrowth dynamics in tropical forests using structural change detection and Landsat time series. Remote Sens. Environ. 2015, 169, 320–334. [Google Scholar] [CrossRef]
  79. Nguyen, H.C.; Jung, J.; Lee, J.; Choi, S.U.; Hong, S.Y.; Heo, J. Optimal Atmospheric Correction for Above-Ground Forest Biomass Estimation with the ETM+ Remote Sensor. Sensors (Basel) 2015, 15, 18865–18886. [Google Scholar] [CrossRef] [Green Version]
  80. Ma, W.; Domke, G.M.; D’Amato, A.W.; Woodall, C.W.; Walters, B.F.; Deo, R.K. Using matrix models to estimate aboveground forest biomass dynamics in the eastern USA through various combinations of LiDAR, Landsat, and forest inventory data. Environ. Res. Lett. 2018, 13, 125004. [Google Scholar] [CrossRef]
  81. Powell, S.L.; Cohen, W.B.; Healey, S.P.; Kennedy, R.E.; Moisen, G.G.; Pierce, K.B.; Ohmann, J.L. 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, 1053–1068. [Google Scholar] [CrossRef]
  82. Powell, S.L.; Cohen, W.B.; Kennedy, R.E.; Healey, S.P.; Huang, C. Observation of Trends in Biomass Loss as a Result of Disturbance in the Conterminous U.S.: 1986–2004. Ecosystems 2013, 17, 142–157. [Google Scholar] [CrossRef]
  83. Badreldin, N.; Sanchez-Azofeifa, A. Estimating Forest Biomass Dynamics by Integrating Multi-Temporal Landsat Satellite Images with Ground and Airborne LiDAR Data in the Coal Valley Mine, Alberta, Canada. Remote Sens. 2015, 7, 2832–2849. [Google Scholar] [CrossRef] [Green Version]
  84. Zhang, J.; Lu, C.; Xu, H.; Wang, G. Estimating aboveground biomass of Pinus densata-dominated forests using Landsat time series and permanent sample plot data. J. For. Res. 2018, 30, 1689–1706. [Google Scholar] [CrossRef]
  85. Cohen, W.B.; Goward, S.N. Landsat’s role in ecological applications of remote sensing. Bioscience 2004, 54, 535–545. [Google Scholar] [CrossRef]
  86. Crist, E.P. A TM tasseled cap equivalent transformation for reflectance factor data. Remote Sens. Environ. 1985, 17, 301–306. [Google Scholar] [CrossRef]
  87. Baig, M.H.A.; Zhang, L.; Shuai, T.; Tong, Q. Derivation of a tasselled cap transformation based on Landsat 8 at-satellite reflectance. Remote Sens. Lett. 2014, 5, 423–431. [Google Scholar] [CrossRef]
  88. Huang, C.; Wylie, B.; Yang, L.; Homer, C.; Zylstra, G. Derivation of a tasselled cap transformation based on Landsat 7 at-satellite reflectance. Int. J. Remote Sens. 2002, 23, 1741–1748. [Google Scholar] [CrossRef]
  89. Kauth, R.J.; Thomas, G. The tasselled cap—A graphic description of the spectral-temporal development of agricultural crops as seen by Landsat. In Proceedings of the LARS Symposia, West Lafayette, IN, USA, 1 July 1976; p. 159. [Google Scholar]
  90. Labrecque, S.; Fournier, R.; Luther, J.; Piercey, D. A comparison of four methods to map biomass from Landsat-TM and inventory data in western Newfoundland. For. Ecol. Manag. 2006, 226, 129–144. [Google Scholar] [CrossRef]
  91. Healey, S.P.; Cohen, W.B.; Zhiqiang, Y.; Krankina, O.N. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection. Remote Sens. Environ. 2005, 97, 301–310. [Google Scholar] [CrossRef]
  92. Wulder, M.A.; White, J.C.; Bater, C.W.; Coops, N.C.; Hopkinson, C.; Chen, G. Lidar plots—A new large-area data collection option: Context, concepts, and case study. Can. J. Remote Sens. 2014, 38, 600–618. [Google Scholar] [CrossRef]
  93. Altman, N.S. An Introduction to Kernel and Nearest-Neighbor Nonparametric Regression. Am. Stat. 1992, 46, 175–185. [Google Scholar] [CrossRef] [Green Version]
  94. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  95. Chirici, G.; Mura, M.; McInerney, D.; Py, N.; Tomppo, E.O.; Waser, L.T.; Travaglini, D.; McRoberts, R.E. A meta-analysis and review of the literature on the k-Nearest Neighbors technique for forestry applications that use remotely sensed data. Remote Sens. Environ. 2016, 176, 282–294. [Google Scholar] [CrossRef]
  96. Eskelson, B.N.I.; Temesgen, H.; Lemay, V.; Barrett, T.M.; Crookston, N.L.; Hudak, A.T. The roles of nearest neighbor methods in imputing missing data in forest inventory and monitoring databases. Scand. J. For. Res. 2009, 24, 235–246. [Google Scholar] [CrossRef] [Green Version]
  97. Liaw, A.; Wiener, M. Classification and regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  98. Crookston, N.L.; Finley, A.O. yaImpute: An R package for kNN imputation. J. Stat. Softw. 2008, 23, 16. [Google Scholar] [CrossRef] [Green Version]
  99. Ho Tong Minh, D.; Ndikumana, E.; Vieilledent, G.; McKey, D.; Baghdadi, N. Potential value of combining ALOS PALSAR and Landsat-derived tree cover data for forest biomass retrieval in Madagascar. Remote Sens. Environ. 2018, 213, 206–214. [Google Scholar] [CrossRef]
  100. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  101. Key, C.; Benson, N. Landscape Assessment: Remote Sensing of Severity, the Normalized Burn Ratio and Ground Measure of Severity, the Composite Burn Index; FIREMON: Fire Effects Monitoring and Inventory System; USDA Forest Service, Rocky Mountain Res. Station: Ogden, UT, USA, 2005. [Google Scholar]
  102. Liu, H.Q.; Huete, A. A feedback based modification of the NDVI to minimize canopy background and atmospheric noise. IEEE Trans. Geosci. Remote Sens. 1995, 33, 457–465. [Google Scholar] [CrossRef]
  103. Huete, A.R. A soil-adjusted vegetation index (SAVI). Remote Sens. Environ. 1988, 25, 295–309. [Google Scholar] [CrossRef]
  104. Hunt, E.R.; Daughtry, C.; Eitel, J.U.; Long, D.S. Remote sensing leaf chlorophyll content using a visible band index. Agron. J. 2011, 103, 1090–1099. [Google Scholar] [CrossRef] [Green Version]
  105. Phua, M.-H.; Saito, H. Estimation of biomass of a mountainous tropical forest using Landsat TM data. Can. J. Remote Sens. 2003, 29, 429–440. [Google Scholar] [CrossRef]
  106. Duane, M.V.; Cohen, W.B.; Campbell, J.L.; Hudiburg, T.; Turner, D.P.; Weyermann, D.L. Implications of alternative field-sampling designs on Landsat-based mapping of stand age and carbon stocks in Oregon forests. For. Sci. 2010, 56, 405–416. [Google Scholar]
Remotesensing 12 00098 g001 550
Figure 1. A common concept for estimating AGB dynamics using LTS data.
Figure 1. A common concept for estimating AGB dynamics using LTS data.
Remotesensing 12 00098 g001
Table
Table 1. Common publicly available and automated change detection algorithms using LTS. See Table A1 for a description of spectral indices.
Table 1. Common publicly available and automated change detection algorithms using LTS. See Table A1 for a description of spectral indices.
AlgorithmMethodTemporal FrequencySpatial ScaleChange IndexChange Type
Vegetation Change Tracker (VTC, [37])ThresholdingAnnualPixelIntegrated Forest Z-score (IFZ)Abrupt
Continuous Change Detection and Classification (CCDC, [38])Statistical boundaryAll available imagesPixelLandsat spectral bandsAbrupt, gradual
Breaks for additive Season and Trend Monitor (BFAST Monitor, [39])Statistical boundaryAll available imagesPixelNormalized difference vegetation index (NDVI)Abrupt, gradual
Landsat-based detection of Trends in Disturbance and Recovery (LandTrendr, [40])Temporal segmentationAnnualPixelNormalized burn ratio (NBR), NDVI, TC components (Brightness (TCB), Greenness (TCG), Wetness (TCW)), spectral bandsAbrupt, gradual
Table
Table 2. Benefits of using LTS to improve the estimation of forest AGB (both single-date and over time).
Table 2. Benefits of using LTS to improve the estimation of forest AGB (both single-date and over time).
BenefitSpecific Improvement/FindingReference
Filling spatial and temporal data gaps in AGB predictionsImage composites derived from LTS are necessary for spatially complete estimations of AGB over large areas, regardless modelling approach adoptedZald, Wulder [48]; Matasci, Hermosilla [49]
Demonstrating the utility of LTS-derived image composites in consistently monitoring AGB over long time-periods (20–40 years)Boisvenue, Smiley [50]; Matasci, Hermosilla [51]; Kennedy, Ohmann [9]; Nguyen, Jones [52]
Utilising LTS-derived image composites for AGB predictions in different contextsMorel, Fisher [53]; Frazier, Coops [54]; Wilson, Knight [55]; Nguyen, Jones [56]
Filling spatial and temporal data gaps in LTS using a pixel-based temporal fitting processDeo, Russell [57]; Deo, Russell [58]
Improving the accuracy of AGB modellingUsing disturbance and recovery metrics derived from LTS can improve the accuracy of AGB predictions in comparison with using single-date imagesPflugmacher, Cohen [28]; Pflugmacher, Cohen [12]
Quantifying the capacity of LTS-derived change metrics for estimating forest AGBFrazier, Coops [54]
Demonstrating the utility of LTS-derived change attributions (e.g., disturbance severity and agent) in modelling forest structure and AGBZald, Ohmann [59]; Zald, Wulder [48]; Nguyen, Jones [56]; Bolton, White [60]
Forest age derived from LTS can improve forest AGB estimatesLefsky [61]; Liu, Peng [62]
Identifiers of temporal patterns in spectral trajectories provide improvements in modelling and explaining historical AGB dynamicsGómez, White [63]
Seasonal NDVI time series improving AGB estimations compared with a single-date NDVIZhu and Liu [64]
Fitted LTS data improving the accuracy of AGB predictions compared with raw/observed data.Deo, Russell [58]
A pixel-based temporal fitting process improves the temporal consistency of AGB predictions.Matasci, Hermosilla [51], Deo, Russell [58]; Kennedy, Ohmann [9]
Table
Table 3. A summary of studies that used LTS-based approaches for estimating forest AGB dynamics. See Table A1 for a description of spectral indices.
Table 3. A summary of studies that used LTS-based approaches for estimating forest AGB dynamics. See Table A1 for a description of spectral indices.
Study (Year)Ecosystem, LocationLTS StackModelling ApproachesPredicting AGB and Its DynamicsAccuracy Assessment
Explanatory DataReference DataModelling TechniqueTemporal EstimatesCharacterising AGB Dynamics
Powell, Cohen [81]
(2010)		Coniferous and mixed forests in northern Arizona and Minnesota, USAnnual cloud-free images (1985–2006)-Six spectral bands;
-Spectral indices: NDVI, TCB, TCW, TCG, TC angle (TCA), TC distance (TCD), disturbance index (DI);
-Topographic and climatic variables.		Field inventory (FI) plotsReduced Major Axis regression, Gradient Nearest Neighbour (GNN) imputation (k = 1), Random Forest (RF) regressionAnnual-Smoothing AGB trajectories using linear segmentation algorithms to detect changes
-Creating forest AGB disturbance maps		-Model cross-validation using withheld FI plots
-Validating AGB change using re-measured field plots
-Comparing scene-level AGB trajectories
Powell, Cohen [82]
(2013)		Multiple forest types across 50 Landsat scenes, USAnnual least cloudy images (1986–2004)-Six spectral bands;
-Spectral indices: NDVI, TCB, TCW, TCG, TCA, TCD;
-Topographic and climatic variables.		FI plotsRF regressionAnnual-Smoothing AGB trajectories using LandTrendr
-Integrating AGB loss estimates with VCT-derived disturbance maps
-Estimating national and stratum levels of AGB loss from disturbance.
-Evaluating trends in AGB loss on the context of large forest carbon sink in the US.		-Model assessment using out-of-bag (OOB) errors
-Independent validations of AGB predictions using FI plots and other AGB products
-Comparing AGB loss estimates with results from other studies
Main-Knorn, Cohen [72]
(2013)		Coniferous forests, Western Carpathian Mountains, EuropeAnnual cloud-free images (1985–2010)-Six spectral bands
-Topographic variables		FI plotsRF regressionAnnual-Fitting AGB trajectories using LandTrendr to detect changes
-Mapping AGB loss and gain according to disturbance and recovery trends		-Model assessment using a 10-fold cross-validation
-Validating AGB disturbance maps using TimeSync
Gómez, White [63]
(2014)		Pine forests, Spain
Annual cloud-free images (1984–2009)-Spectral indices: TCA, TCD, NDVI
-Change metrics: dynamic variables derived by transformations of spectral trajectories		Remeasured FI plotsDecision trees (CART)2 maps of 1990 and 2000-Calculating AGB change over the period 1990–2000 as the difference between AGB values of the 2 years with uncertainty estimates.-Model cross-validation using withheld FI plots
-Validating AGB change using re-measured FI plots
Liu, Peng [62]
(2014)		Artificial mixed forests, north-western ChinaNear-annual cloud-free images (1974–2013)-Spectral indices: simple ratio index (SR = NIR/Red)
-Forest afforestation age derived using the VTC algorithm.		FI plotsLinear regressionEight years during 1974–2013-Calculating mean AGB by different afforestation ages over the 40 years
-Temporal statistics of biomass at a regional level		-Leave-one-out cross-validation model
Pflugmacher, Cohen [12]
(2014)		Mixed-conifer forests, Oregon, USAnnual cloud-free images (1972–2009)-Spectral indices: TCA, TCD
-Change metrics: LandTrendr-derived relative change magnitude, pre- and post-change conditions; and TSD		Airborne Lidar-based plotsRF regressionAnnual-Estimating AGB change by subtracting predictions between two points of time-Model assessment using RF OOB errors
-Cross-validation of AGB and its changes using re-measured field plots
Badreldin and Sanchez-Azofeifa [83]
(2015)		Coniferous forests, Coal Valley Mine, Alberta, CanadaAnnual cloud-free composites (1988–2011)-Spectral indices: NDVI, enhanced vegetation index (EVI and EVI2), chlorophyll vegetation index (CVI), TCB, TCG, TCWGround Lidar plotsMultiple linear regressionFive years during 1988–2011-Estimating AGB change as differences between two points of time
-Trend analysis using the Mann-Kendall test		-Model validation using ground Lidar plots
Nguyen, Jung [79]
(2015)		Mixed coniferous forests, South KoreaFive seasonal images (2010–2011)-Six spectral bandsFI plotskNN regression
(k = 1–20)		Seasonal-None-Accuracy assessment using 10-fold cross-validation
Boisvenue, Smiley [50]
(2016)		Boreal forests, Saskatchewan, CanadaAnnual cloud-free composites (1984–2012)-Six fitted spectral bands
-Topographic variables		Remeasured FI plots (1381 plots since 1949)RF regressionAnnual-Computing annual AGB changes from annual AGB maps
-Developing a mixed-effect model to summarize long-term trends of AGB at a stratum-level		-Model assessment using RF OOB accuracy
-Validating predictions of AGB change using remeasured FI plots
Deo, Russell [58]
(2017)		Mixed forests, Minnesota, the USNear-annual cloud-free images (1990–2011)-Fitted spectral band 5
-Fitted spectral indices: NDVI, NBR, TCA, DI, IFZ		Remeasured FI plots (2000–2011)RF-based k-Nearest Neighbour (kNN) imputation
(k = 1,3,5)		Six years during 1990–2011-Estimating AGB change as the differences between temporal predicted maps.-Model cross-validation using FI plots (2000,2010)
-Comparing AGB predictions with other products at the stand-level
Zhang, Lu [84]
(2018)		Pinus densata-dominated forests, southwest ChinaFive-year intervals least cloudy TM images
(1987–2007)		-Six spectral bands
-Spectral indices: SRs, normalized ratios (NDs), difference vegetation index (DVI), NDVI, TCW, TCB, TCG, linear transform of multiple bands (VIS123, MID57)
-Texture measures of 5×5 and 9×9 pixel-windows
-Topographic variables		Remeasured FI plots
(5-year cycles during 1987–2007)		Multiple Linear Regression, Partial Least Square Regression, Gradient Boost Regression Tree, and RF regressionFive years during 1987–2007-Directly modelling AGB dynamics by linking changes of spectral variables with changes in AGB measurements
-Combining AGB change maps with land cover classification maps and estimating the total AGB change for each of 5-year intervals		-Validating AGB predictions and changes using remeasured FI plots.
Ma, Domke [80]
(2018)		Mixed deciduous and coniferous forests in the eastern USAnnual cloud-free images (1984–2015)-Spectral band 4 (SWIR)
-Spectral indices: DI, EVI, IFZ, NDVI, NBR, TCW, TCA, TCB, TCG, normalized difference moisture index (NDMI), soil adjusted vegetation index (SAVI)		Remeasured FI plotsMatrix models -Using matrix models to project future AGB dynamics-Evaluating the accuracy of matrix models using FI plots
-Using fuzzy sets representing uncertainty
Kennedy, Ohmann [9]
(2018)		Coniferous forests, Western USAnnual cloud-free composites (1990–2012)-LandTrendr-derived data:
-Fitted spectral indices: NBR, TCB, TCG, TCW
-Spectral change metrics, disturbance causal agents
-Topographic and climatic variables.		Multi-temporal remeasured FI plots (1991–2011)GNN imputation
(k = 1,2,3)		Annual-Consistently tracking AGB changes by forest disturbance and recovery
-Estimating AGB loss according to disturbance agents with uncertainties.		-Estimating uncertainty based on different scenarios of LandTrendr and GNN
-Evaluating AGB outcomes using 6018 FI plots across the study area
-Comparing AGB predictions with airborne Lidar AGB data
Matasci, Hermosilla [51]
(2018)		Canadian forest ecosystems (650 million ha)Annual cloud-free composites (1984–2016)-Fitted spectral indices: TCB, TCG, TCA, NBR
-Change metrics: years since the greatest change
-Topographic and climatic variables; pixel locations		Airborne Lidar-based plotsRF-based kNN imputation
(k = 1)		Annual-Integrating AGB prediction maps with LTS-derived disturbance maps
-Sample-based summarising AGB dynamics for three scenarios: undisturbed forests, forests impacted by wildfires, and harvesting		-Model cross-validation using Lidar-based plots
-Multi-temporal, stand-level assessments using airborne Lidar data
Nguyen, Jones [52]
(2020)		Mixed temperate forests, Victoria, AustraliaAnnual cloud-free composites (1988–2017)-Fitted spectral indices: NBR, TCB, TCG, TCA
-Change metrics: disturbance levels and causal agents, TSD
-Topographic and climatic variables; pixel locations (X and Y)		FI plotsRF-based kNN imputation
(k = 1)		Annual-Fitting temporal AGB trajectories according to fitted NBR trajectories.
-Integrating AGB trajectories with forest disturbance and recovery history
-Characterising spatial and temporal patterns of AGB dynamics using spatial change metrics, including: AGB loss and gain from disturbance and recovery, Recovery Indicator (RI = AGB gain/AGB loss), Years to Recovery (Y2R)
-Determining how AGB responds to different disturbance scenarios (severities and causal agents)		-Internal assessment using bootstrapping sample
-Time-series validations of AGB predictions using airborne multi-temporal Lidar data
-Validating AGB changes according to forest disturbance and recovery history using Lidar-based AGB data

Share and Cite

MDPI and ACS Style

H. Nguyen, T.; Jones, S.; Soto-Berelov, M.; Haywood, A.; Hislop, S. Landsat Time-Series for Estimating Forest Aboveground Biomass and Its Dynamics across Space and Time: A Review. Remote Sens. 2020, 12, 98. https://doi.org/10.3390/rs12010098

AMA Style

H. Nguyen T, Jones S, Soto-Berelov M, Haywood A, Hislop S. Landsat Time-Series for Estimating Forest Aboveground Biomass and Its Dynamics across Space and Time: A Review. Remote Sensing. 2020; 12(1):98. https://doi.org/10.3390/rs12010098

Chicago/Turabian Style

H. Nguyen, Trung, Simon Jones, Mariela Soto-Berelov, Andrew Haywood, and Samuel Hislop. 2020. "Landsat Time-Series for Estimating Forest Aboveground Biomass and Its Dynamics across Space and Time: A Review" Remote Sensing 12, no. 1: 98. https://doi.org/10.3390/rs12010098

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop