Abstract:Understanding change in forest carbon (C) is important for devising strategies to reduce emissions of greenhouse gases. National forest inventories (NFIs) are important to meet international accounting goals, but data are often incomplete going back in time, and the amount of time between remeasurements can make attribution of C flux to specific events difficult. The long time series of Landsat imagery provides spatially comprehensive, consistent information that can be used to fill the gaps in ground measurements with predictive models. To evaluate such models, we relate Landsat spectral changes and disturbance interpretations directly to C flux measured on NFI plots and compare the performance of models with and without ground-measured predictor variables. The study was conducted in the forests of southwest Oregon State, USA, a region of diverse forest types, disturbances, and landowner management objectives. Plot data consisted of 676 NFI plots with remeasured individual tree data over a mean interval (time 1 to time 2) of 10.0 years. We calculated change in live aboveground woody carbon (AWC), including separate components of growth, mortality, and harvest. We interpreted radiometrically corrected annual Landsat images with the TimeSync (TS) tool for a 90 m × 90 m area over each plot. Spectral time series were divided into segments of similar trajectories and classified as disturbance, recovery, or stability segments, with type of disturbance identified. We calculated a variety of values and segment changes from tasseled cap angle and distance (TCA and TCD) as potential predictor variables of C flux. Multiple linear regression was used to model AWC and net change in AWC from the TS change metrics. The TS attribution of disturbance matched the plot measurements 89% of the time regarding whether fire or harvest had occurred or not. The primary disagreement was due to plots that had been partially cut, mostly in vigorous stands where the net change in AWC over the measurement was positive in spite of cutting. The plot-measured AWC at time 2 was 86.0 ± 78.7 Mg C ha−1 (mean and standard deviation), and the change in AWC across all plots was 3.5 ± 33 Mg C ha−1 year−1. The best model for AWC based solely on TS and other mapped variables had an R2 = 0.52 (RMSE = 54.6 Mg C ha−1); applying this model at two time periods to estimate net change in AWC resulted in an R2 = 0.25 (RMSE = 28.3 Mg ha−1) and a mean error of −5.4 Mg ha−1. The best model for AWC at time 2 using plot measurements at time 1 and TS variables had an R2 = 0.95 (RSME = 17.0 Mg ha−1). The model for net change in AWC using the same data was identical except that, because the variable being estimated was smaller in magnitude, the R2 = 0.73. All models performed better at estimating net change in AWC on TS-disturbed plots than on TS-undisturbed plots. The TS discrimination of disturbance between fire and harvest was an important variable in the models because the magnitude of spectral change from fire was greater for a given change in AWC. Regional models without plot-level predictors produced erroneous predictions of net change in AWC for some of the forest types. Our study suggests that, in spite of the simplicity of applying a single carbon model to multiple image dates, the approach can produce inaccurate estimates of C flux. Although models built with plot-level predictors are necessarily constrained to making predictions at plot locations, they show promise for providing accurate updates or back-calculations of C flux assessments.

A new approach for Spatializing the CAnadian National Forest Inventory (SCANFI) using Landsat dense time series

Improved k-NN Mapping of Forest Attributes in Northern Canada Using Spaceborne L-Band SAR, Multispectral and LiDAR Data

Monitoring seedling stands using national forest inventory and multispectral airborne laser scanning data

Updating Forest Stand Inventories: Integration of Photo-Interpreted and Airborne Laser Scanning Forest Attributes Using Generic Region Merging Segmentation and kNN Imputation

Spatial distribution of carbon sources and sinks in Canada's forests

Mapping vegetation height and identifying the northern forest limit across Canada using ICESat-2, Landsat time series and topographic data

Fine-Scale Mapping of Natural Ecological Communities Using Machine Learning Approaches

Harmonizing GEDI and LVIS Data for Accurate and Large-Scale Mapping of Foliage Height Diversity

Automated forest inventory: Analysis of high-density airborne LiDAR point clouds with 3D deep learning

Integrating TimeSync Disturbance Detection and Repeat Forest Inventory to Predict Carbon Flux

Combining visual interpretation and image segmentation to derive canopy cover index from high resolution satellite imagery in functionally diverse coniferous forests

Fine-resolution Forest Tree Height Estimation Across the Sierra Nevada Through the Integration of Spaceborne LiDAR, Airborne LiDAR, and Optical Imagery.

Ecoregion-wise fractional mapping of tree functional composition in temperate mixed forests with sentinel data: Integrating time-series spectral and radar data

Characterizing post-fire northern boreal forest height dynamics

Forest Change Mapping using Multi-Source Satellite SAR, Optical, and LiDAR Remote Sensing Data

A Deep Learning Approach to Estimate Canopy Height and Uncertainty by Integrating Seasonal Optical, SAR and Limited GEDI LiDAR Data over Northern Forests

Annual maps of forest and evergreen forest in the contiguous United States during 2015–2017 from analyses of PALSAR-2 and Landsat images

Use of Remote Sensing Data to Improve the Efficiency of National Forest Inventories: A Case Study from the United States National Forest Inventory

Land cover classification in an era of big and open data: Optimizing localized implementation and training data selection to improve mapping outcomes

A Vegetation Mapping Strategy for Conifer Forests by Combining Airborne LiDAR Data and Aerial Imagery

Development of a Forest Health Index Using Multispectral Airborne Digital Camera Imagery