Forest aboveground biomass mapping and estimation across multiple spatial scales using model-based inference

Remotely sensed data have been widely used in recent years for mapping and estimating biomass. However, the characterization of the uncertainty of mapped or estimated biomass in previous studies was either based on ad-hoc approaches (e.g., using model fitting statistics such root mean square errors derived from purposive samples) or mostly limited to the analysis of mean biomass for the whole study area. This study proposed a novel uncertainty analysis method that can characterize biomass uncertainty across multiple spatial scales and multiple spatial resolutions. The uncertainty analysis method built on model-based inference and can propagate errors from trees to field plots, individual pixels, and small areas or large regions that consist of multiple pixels (up to all pixels within a study area). We developed and tested this method over northern Minnesota forest areas of approximately 69,508 km2 via a unique combination of several datasets for biomass mapping and estimation: wall-to-wall airborne lidar data, national forest inventory (NFI) plots, and destructive measurements of tree aboveground biomass (AGB). We found that the pixel-level AGB prediction error is dominated by lidar-based AGB model residual errors when the spatial resolution is near 380 m or finer and by model parameter estimate errors when the spatial resolution is coarser. We also found that the relative error of AGB predicted from lidar can be reduced to approximately 11% (or mean 5.1 Mg/ha; max 43.6 Mg/ha) at one-hectare scale (or at 100 m spatial resolution) over our study area. Because our uncertainty analysis method uses model-based inference and does not require probability samples of field plots, our methodology has potential applications worldwide, especially over tropics and developing countries where NFI systems are not well-established.