Evaluating the remote sensing and inventory-based estimation of biomass in the western carpathians

Understanding the potential of forest ecosystems as global carbon sinks requires a thorough knowledge of forest carbon dynamics, including both sequestration and fluxes among multiple pools. The accurate quantification of biomass is important to better understand forest productivity and carbon cycling dynamics. Stand-based inventories (SBIs) are widely used for quantifying forest characteristics and for estimating biomass, but information may quickly become outdated in dynamic forest environments. Satellite remote sensing may provide a supplement or substitute. We tested the accuracy of aboveground biomass estimates modeled from a combination of Landsat Thematic Mapper (TM) imagery and topographic data, as well as SBI-derived variables in a Picea abies forest in the Western Carpathian Mountains. We employed Random Forests for non-parametric, regression tree-based modeling. Results indicated a difference in the importance of SBI-based and remote sensing-based predictors when estimating aboveground biomass. The most accurate models for biomass prediction ranged from a correlation coefficient of 0.52 for the TM- and topography-based model, to 0.98 for the inventory-based model. While Landsat-based biomass estimates were measurably less accurate than those derived from SBI, adding tree height or stand-volume as a field-based predictor to TM and topography-based models increased performance to 0.36 and 0.86, respectively. Our results illustrate the potential of spectral data to reveal spatial details in stand structure and ecological complexity. © 2011 by the authors

United States Forest Disturbance Trends Observed Using Landsat Time Series

Disturbance events strongly affect the composition, structure, and function of forest ecosystems; however, existing U.S. land management inventories were not designed to monitor disturbance. To begin addressing this gap, the North American Forest Dynamics (NAFD) project has examined a geographic sample of 50 Landsat satellite image time series to assess trends in forest disturbance across the conterminous United States for 1985-2005. The geographic sample design used a probability-based scheme to encompass major forest types and maximize geographic dispersion. For each sample location disturbance was identified in the Landsat series using the Vegetation Change Tracker (VCT) algorithm. The NAFD analysis indicates that, on average, 2.77 Mha/yr of forests were disturbed annually, representing 1.09%/yr of US forestland. These satellite-based national disturbance rates estimates tend to be lower than those derived from land management inventories, reflecting both methodological and definitional differences. In particular the VCT approach used with a biennial time step has limited sensitivity to low-intensity disturbances. Unlike prior satellite studies, our biennial forest disturbance rates vary by nearly a factor of two between high and low years. High western US disturbance rates were associated with active fire years and insect activity, while variability in the east is more strongly related to harvest rates in managed forests. We note that generating a geographic sample based on representing forest type and variability may be problematic since the spatial pattern of disturbance does not necessarily correlate with forest type. We also find that the prevalence of diffuse, non-stand clearing disturbance in US forests makes the application of a biennial geographic sample problematic. Future satellite-based studies of disturbance at regional and national scales should focus on wall-to-wall analyses with annual time step for improved accuracy

Forest structure and aboveground biomass in the southwestern United States from MODIS and MISR

Red band bidirectional reflectance factor data from the NASA MODerate resolution Imaging Spectroradiometer (MODIS) acquired over the southwestern United States were interpreted through a simple geometric–optical (GO) canopy reflectance model to provide maps of fractional crown cover (dimensionless), mean canopy height (m), and aboveground woody biomass (Mg ha−1) on a 250 m grid. Model adjustment was performed after dynamic injection of a background contribution predicted via the kernel weights of a bidirectional reflectance distribution function (BRDF) model. Accuracy was assessed with respect to similar maps obtained with data from the NASA Multiangle Imaging Spectroradiometer (MISR) and to contemporaneous US Forest Service (USFS) maps based partly on Forest Inventory and Analysis (FIA) data. MODIS and MISR retrievals of forest fractional cover and mean height both showed compatibility with the USFS maps, with MODIS mean absolute errors (MAE) of 0.09 and 8.4 m respectively, compared with MISR MAE of 0.10 and 2.2 m, respectively. The respective MAE for aboveground woody biomass was ~10 Mg ha−1, the same as that from MISR, although the MODIS retrievals showed a much weaker correlation, noting that these statistics do not represent evaluation with respect to ground survey data. Good height retrieval accuracies with respect to averages from high resolution discrete return lidar data and matches between mean crown aspect ratio and mean crown radius maps and known vegetation type distributions both support the contention that the GO model results are not spurious when adjusted against MISR bidirectional reflectance factor data. These results highlight an alternative to empirical methods for the exploitation of moderate resolution remote sensing data in the mapping of woody plant canopies and assessment of woody biomass loss and recovery from disturbance in the southwestern United States and in parts of the world where similar environmental conditions prevail

Assessing North American Forest Disturbance from the Landsat Archive

Forest disturbances are thought to play a major role in controlling land-atmosphere fluxes of carbon. Under the auspices of the North American Carbon Program, the LEDAPS (Landsat Ecosystem Disturbance Adaptive Processing System) and NACP-FIA projects have been analyzing the Landsat satellite record to assess rates of forest disturbance across North America. In the LEDAPS project, wall-to-wall Landsat imagery for the period 1975-2000 has been converted to surface reflectance and analyzed for decadal losses (disturbance) or gains (regrowth) in biomass using a spectral "disturbance index". The NACP-FIA project relies on a geographic sample of dense Landsat image time series, allowing both disturbance rates and recovery trends to be characterized. Preliminary results for the 1990's indicate high rates of harvest within the southeastern US, Eastern Canada, and the Pacific Northwest, with spatially averaged (approx.50x50 km) turnover periods as low as 25-40 years. Lower rates of disturbance are found in the Rockies and Northeastern US

Assessing small area estimates via artificial populations from KBAABB: a kNN-based approximation to ABB

Comparing and evaluating small area estimation (SAE) models for a given application is inherently difficult. Typically, we do not have enough data in many areas to check unit-level modeling assumptions or to assess unit-level predictions empirically; and there is no ground truth available for checking area-level estimates. Design-based simulation from artificial populations can help with each of these issues, but only if the artificial populations (a) realistically represent the application at hand and (b) are not built using assumptions that could inherently favor one SAE model over another. In this paper, we borrow ideas from random hot deck, approximate Bayesian bootstrap (ABB), and k nearest neighbor (kNN) imputation methods, which are often used for multiple imputation of missing data. We propose a kNN-based approximation to ABB (KBAABB) for a different purpose: generating an artificial population when rich unit-level auxiliary data is available. We introduce diagnostic checks on the process of building the artificial population itself, and we demonstrate how to use such an artificial population for design-based simulation studies to compare and evaluate SAE models, using real data from the Forest Inventory and Analysis (FIA) program of the US Forest Service. We illustrate how such simulation studies may be disseminated and explored interactively through an online R Shiny application

Improving estimates of forest disturbance by combining observations from Landsat time series with U.S. Forest Service Forest Inventory and Analysis data

With earth's surface temperature and human population both on the rise a new emphasis has been placed on monitoring changes to forested ecosystems the world over. In the United States the U.S. Forest Service Forest Inventory and Analysis (FIA) program monitors the forested land base with field data collected over a permanent network of sample plots. Although these plots are visited repeatedly through time there are large temporal gaps (e.g. 5–10 years) between remeasurements such that many forest canopy disturbances go undetected. In this paper we demonstrate how Landsat time series (LTS) can help improve FIA's capacity to estimate disturbance by 1.) incorporating a new, downward looking response variable which is more sensitive to picking up change and 2.) providing historical disturbance maps which can reduce the variance of design-based estimates via post-stratification. To develop the LTS response variable a trained analyst was used to manually interpret 449 forested FIA plots located in the Uinta Mountains of northern Utah, USA. This involved recording cause and timing of disturbances based on evidence gathered from a 26-year annual stack of Landsat images and an 18-year, periodically spaced set of high resolution (~1 m) aerial photographs (e.g. National Aerial Image Program, NAIP and Google Earth). In general, the Landsat data captured major disturbances (e.g. harvests, fires) while the air photos allowed more detailed estimates of the number of trees impacted by recent insect outbreaks. Comparing the LTS and FIA field observations, we found that overall agreement was 73%, although when only disturbed plots were considered agreement dropped to 40%. Using the non-parametric Mann–Whitney test, we compared distributions of live and disturbed tree size (height and DBH) and found that when LTS and FIA both found non-stand clearing disturbance the median disturbed tree size was significantly larger than undisturbed trees, whereas no significant difference was found on plots where only FIA detected disturbance. This suggests that LTS interpretation and FIA field crews both detect upper canopy disturbances while FIA crews alone add disturbances occurring at or below canopy level. The analysis also showed that plots with only LTS disturbance had a significantly greater median number of years since last FIA measurement (6 years) than plots with both FIA and LTS disturbances (2.5 years), indicating that LTS improved detection on plots which had not been field sampled for several years. Next, to gauge the impact of incorporating LTS disturbances into the FIA estimation process we calculated design-based estimates of disturbance (for the period 1995–2011) using three response populations 1.) LTS observations, 2.) FIA field observations, and 3.) Combination of FIA and LTS observations. The results showed that combining the FIA and LTS observations led to the largest and most precise (i.e. smallest percent standard error) estimates of disturbance. In fact, the estimate based on the combined observations (486,458 ha, +/−47,101) was approximately 65% more than the estimate derived solely with FIA data (294,295 ha, +/−44,242). Lastly, a Landsat forest disturbance map was developed and tested for its ability to post-stratify the design-based estimates. Based on relative efficiency (RE), we found that stratification mostly improved the estimates derived with the LTS response data. Aside from insects (RE = 1.26), the estimates of area affected by individual agents saw minimal gain, whereas the LTS and combined FIA + LTS estimates of total disturbance saw modest improvement, with REs of 1.43 and 1.50 respectively. Overall, our results successfully demonstrate two ways LTS can improve the completeness and precision of disturbance estimates derived from FIA inventory data.Keywords: Disturbance mapping, Landsat time series, Design-based estimation, Forest disturbance, Post-stratification, Forest Inventory and Analysis (FIA)Keywords: Disturbance mapping, Landsat time series, Design-based estimation, Forest disturbance, Post-stratification, Forest Inventory and Analysis (FIA

Changes in timber haul emissions in the context of shifting forest management and infrastructure

<p>Abstract</p> <p>Background</p> <p>Although significant amounts of carbon may be stored in harvested wood products, the extraction of that carbon from the forest generally entails combustion of fossil fuels. The transport of timber from the forest to primary milling facilities may in particular create emissions that reduce the net sequestration value of product carbon storage. However, attempts to quantify the effects of transport on the net effects of forest management typically use relatively sparse survey data to determine transportation emission factors. We developed an approach for systematically determining transport emissions using: 1) -remotely sensed maps to estimate the spatial distribution of harvests, and 2) - industry data to determine landscape-level harvest volumes as well as the location and processing totals of individual mills. These data support spatial network analysis that can produce estimates of fossil carbon released in timber transport.</p> <p>Results</p> <p>Transport-related emissions, evaluated as a fraction of transported wood carbon at 4 points in time on a landscape in western Montana (USA), rose from 0.5% in 1988 to 1.7% in 2004 as local mills closed and spatial patterns of harvest shifted due to decreased logging on federal lands.</p> <p>Conclusion</p> <p>The apparent sensitivity of transport emissions to harvest and infrastructure patterns suggests that timber haul is a dynamic component of forest carbon management that bears further study both across regions and over time. The monitoring approach used here, which draws only from widely available monitoring data, could readily be adapted to provide current and historical estimates of transport emissions in a consistent way across large areas.</p

Forest canopy height from the Multiangle Imaging SpectroRadiometer (MISR) assessed with high resolution discrete return lidar

In this study retrievals of forest canopy height were obtained through adjustment of a simple geometricoptical (GO) model against red band surface bidirectional reflectance estimates from NASA\u27s Multiangle Imaging SpectroRadiometer (MISR), mapped to a 250 m grid. The soil-understory background contribution was partly isolated prior to inversion using regression relationships with the isotropic, geometric, and volume scattering kernel weights of a Li-Ross kernel-driven bidirectional reflectance distribution function (BRDF) model. The height retrievals were assessed using discrete return lidar data acquired over sites in Colorado as part of the Cold Land Processes Experiment (CLPX) and used with fractional crown cover retrievals to obtain aboveground woody biomass estimates. For all model runs with reasonable backgrounds and initial b/r (vertical to horizontal crown radii) values \u3c2.0, root mean square error (RMSE) distributions were centered between 2.5 and 3.7 m while R2 distributions were centered between 0.4 and 0.7. The MISR/ GO aboveground biomass estimates predicted via regression on fractional cover and mean canopy height for the CLPX sites showed good agreement with U.S. Forest Service Interior West map data (adjusted R2=0.84). The implication is that multiangle sensors such as MISR can provide spatially contiguous retrievals of forest canopy height, cover, and aboveground woody biomass that are potentially useful in mapping distributions of aboveground carbon stocks, tracking disturbance, and in initializing, constraining, and validating ecosystem models. This is important because the MISR record is spatially comprehensive and extends back to the year 2000 and the launch of the NASA Earth Observing System (EOS) Terra satellite; it might thus provide a ~10-year baseline record that would enhance exploitation of data from the NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission, as well as furthering realization of synergies with active instruments