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

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