Remote Sensing Support for the Gain-Loss Approach for Greenhouse Gas Inventories

For tropical countries that do not have extensive ground sampling programs such as national forest inventories, the gain-loss approach for greenhouse gas inventories is often used. With the gain-loss approach, emissions and removals are estimated as the product of activity data defined as the areas of human-caused emissions and removals and emissions factors defined as the per unit area responses of carbon stocks for those activities. Remotely sensed imagery and remote sensing-based land use and land use change maps have emerged as crucial information sources for facilitating the statistically rigorous estimation of activity data. Similarly, remote sensing-based biomass maps have been used as sources of auxiliary data for enhancing estimates of emissions and removals factors and as sources of biomass data for remote and inaccessible regions. The current status of statistically rigorous methods for combining ground and remotely sensed data that comply with the good practice guidelines for greenhouse gas inventories of the Intergovernmental Panel on Climate Change is reviewed.For tropical countries that do not have extensive ground sampling programs such as national forest inventories, the gain-loss approach for greenhouse gas inventories is often used. With the gain-loss approach, emissions and removals are estimated as the product of activity data defined as the areas of human-caused emissions and removals and emissions factors defined as the per unit area responses of carbon stocks for those activities. Remotely sensed imagery and remote sensing-based land use and land use change maps have emerged as crucial information sources for facilitating the statistically rigorous estimation of activity data. Similarly, remote sensing-based biomass maps have been used as sources of auxiliary data for enhancing estimates of emissions and removals factors and as sources of biomass data for remote and inaccessible regions. The current status of statistically rigorous methods for combining ground and remotely sensed data that comply with the good practice guidelines for greenhouse gas inventories of the Intergovernmental Panel on Climate Change is reviewed.Peer reviewe

Estimating biomass and soil carbon change at the level of forest stands using repeated forest surveys assisted by airborne laser scanner data

BackgroundUnder the growing pressure to implement mitigation actions, the focus of forest management is shifting from a traditional resource centric view to incorporate more forest ecosystem services objectives such as carbon sequestration. Estimating the above-ground biomass in forests using airborne laser scanning (ALS) is now an operational practice in Northern Europe and is being adopted in many parts of the world. In the boreal forests, however, most of the carbon (85%) is stored in the soil organic (SO) matter. While this very important carbon pool is "invisible" to ALS, it is closely connected and feeds from the growing forest stocks. We propose an integrated methodology to estimate the changes in forest carbon pools at the level of forest stands by combining field measurements and ALS data.ResultsALS-based models of dominant height, mean diameter, and biomass were fitted using the field observations and were used to predict mean tree biophysical properties across the entire study area (50 km(2)) which was in turn used to estimate the biomass carbon stocks and the litter production that feeds into the soil. For the soil carbon pool estimation, we used the Yasso15 model. The methodology was based on (1) approximating the initial soil carbon stocks using simulations; (2) predicting the annual litter input based on the predicted growing stocks in each cell; (3) predicting the soil carbon dynamics of the annual litter using the Yasso15 soil carbon model. The estimated total carbon change (standard errors in parenthesis) for the entire area was 0.741 (0.14) Mg ha(-1) yr(-1). The biomass carbon change was 0.405 (0.13) Mg ha(-1) yr(-1), the litter carbon change (e.g., deadwood and leaves) was 0.346 (0.027) Mg ha(-1) yr(-1), and the change in SO carbon was - 0.01 (0.003) Mg ha(-1) yr(-1).ConclusionsOur results show that ALS data can be used indirectly through a chain of models to estimate soil carbon changes in addition to changes in biomass at the primary level of forest management, namely the forest stands. Having control of the errors contributed by each model, the stand-level uncertainty can be estimated under a model-based inferential approach

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI

A practical measure for determining if diameter (D) and height (H) should be combined into D2H in allometric biomass models

Tree diameter at breast height (D) and tree height (H) are often used as predictors of individual tree biomass. Because D and H are correlated, the combined variable D2H is frequently used in regression models instead of two separate independent variables, to avoid collinearity related issues. The justification for D2H is that aboveground biomass is proportional to the volume of a cylinder of diameter, D, and height, H. However, the D2H predictor constrains the model to produce parameter estimates for D and H that have a fixed ratio, in this case, 2.0. In this paper we investigate the degree to which the D2H predictor reduces prediction accuracy relative to D and H separately and propose a practical measure, Q-ratio, to guide the decision as to whether D and H should or should not be combined into D2H. Using five training biomass datasets and two fitting approaches, weighted nonlinear regression and linear regression following logarithmic transformations, we showed that the D2H predictor becomes less efficient in predicting aboveground biomass as the Q-ratio deviates from 2.0. Because of the model constraint, the D2H-based model performed less well than the separate variable model by as much as 12 per cent with regard to mean absolute percentage residual and as much as 18 per cent with regard to sum of squares of log accuracy ratios. For the analysed datasets, we observed a wide variation in Q-ratios, ranging from 2.5 to 5.1, and a large decrease in efficiency for the combined variable model. Therefore, we recommend using the Q-ratio as a measure to guide the decision as to whether D and H may be combined further into D2H without the adverse effects of loss in biomass prediction accuracy

A Comparative Study between Two Regression Methods on LiDAR Data: A Case Study

Airborne LiDAR (Light Detection and Ranging) has become an excellent tool for accurately assessing vegetation characteristics in forest environments. Previous studies showed empirical relationships between LiDAR and field-measured biophysical variables. Multiple linear regression (MLR) with stepwise feature selection is the most common method for building estimation models. Although this technique has provided very interesting results, many other data mining techniques may be applied. The overall goal of this study is to compare different methodologies for assessing biomass fractions at stand level using airborne Li- DAR data in forest settings. In order to choose the best methodology, a comparison between two different feature selection techniques (stepwise selection vs. genetic-based selection) is presented. In addition, classical MLR is also compared with regression trees (M5P). The results when each methodology is applied to estimate stand biomass fractions from an area of northern Spain show that genetically-selected M5P obtains the best results

Taking Stock of Circumboreal Forest Carbon With Ground Measurements, Airborne and Spaceborne LiDAR

The boreal forest accounts for one-third of global forests, but remains largely inaccessible to ground-based measurements and monitoring. It contains large quantities of carbon in its vegetation and soils, and research suggests that it will be subject to increasingly severe climate-driven disturbance. We employ a suite of ground-, airborne- and space-based measurement techniques to derive the first satellite LiDAR-based estimates of aboveground carbon for the entire circumboreal forest biome. Incorporating these inventory techniques with uncertainty analysis, we estimate total aboveground carbon of 38 +/- 3.1 Pg. This boreal forest carbon is mostly concentrated from 50 to 55degN in eastern Canada and from 55 to 60degN in eastern Eurasia. Both of these regions are expected to warm >3 C by 2100, and monitoring the effects of warming on these stocks is important to understanding its future carbon balance. Our maps establish a baseline for future quantification of circumboreal carbon and the described technique should provide a robust method for future monitoring of the spatial and temporal changes of the aboveground carbon content

Estimation of Aboveground Biomass in Alpine Forests: A Semi-Empirical Approach Considering Canopy Transparency Derived from Airborne LiDAR Data

In this study, a semi-empirical model that was originally developed for stem volume estimation is used for aboveground biomass (AGB) estimation of a spruce dominated alpine forest. The reference AGB of the available sample plots is calculated from forest inventory data by means of biomass expansion factors. Furthermore, the semi-empirical model is extended by three different canopy transparency parameters derived from airborne LiDAR data. These parameters have not been considered for stem volume estimation until now and are introduced in order to investigate the behavior of the model concerning AGB estimation. The developed additional input parameters are based on the assumption that transparency of vegetation can bemeasured by determining the penetration of the laser beams through the canopy. These parameters are calculated for every single point within the 3D point cloud in order to consider the varying properties of the vegetation in an appropriate way. Exploratory Data Analysis (EDA) is performed to evaluate the influence of the additional LiDAR derived canopy transparency parameters for AGB estimation. The study is carried out in a 560 km2 alpine area in Austria, where reference forest inventory data and LiDAR data are available. The investigations show that the introduction of the canopy transparency parameters does not change the results significantly according to R2 (R2 = 0.70 to R2 = 0.71) in comparison to the results derived from, the semi-empirical model, which was originally developed for stem volume estimation

Geovisualization of forest simulation modelling results: A case study of carbon sequestration and biodiversity

Sustainable forest management requires new tools to analyze spatial and temporal forest dynamics and to examine those forest parameters that are related to sustainability. We built a prototype system for data analysis and decision-making at forest enterprise level by integrating a forest ecosystem model EFIMOD-PRO (long-term prediction of forest growth and soil development) with an interactive visualization system CommonGIS for analysis of spatially and temporally related data. Using the prototype, a case study in Central European Russia simulated four silvicultural regimes over 200 years: natural development, selective forestry, legal forestry according to the Russian forestry legislation, and illegal forest practice. Exploratory analysis of the simulation results demonstrated that (1) natural stand development is the best alternative for carbon sequestration; (2) legal forest management is the best regime for timber production; (3) selective forestry combines the advantages of two previous strategies, and can be the best strategy for implementing sustainable forest management; and (4) illegal forest practices lead to a fast decrease in forest productivity and decreasing biodiversity. Interactive and dynamic visualizations with maps and statistical graphics played a crucial role in data cleaning, model validation, and analysis of simulation results. The case study demonstrated the potential of integrating forest ecosystem models with exploratory data visualization for the analysis and expert evaluation at the local level. The prototype can be used to present ecological and silvicultural consequences of various management practices to stakeholders and differing social groups, thus stimulating effective decision-making for sustainable forestry

Accounting for density reduction and structural loss in standing dead trees: Implications for forest biomass and carbon stock estimates in the United States

<p>Abstract</p> <p>Background</p> <p>Standing dead trees are one component of forest ecosystem dead wood carbon (C) pools, whose national stock is estimated by the U.S. as required by the United Nations Framework Convention on Climate Change. Historically, standing dead tree C has been estimated as a function of live tree growing stock volume in the U.S.'s National Greenhouse Gas Inventory. Initiated in 1998, the USDA Forest Service's Forest Inventory and Analysis program (responsible for compiling the Nation's forest C estimates) began consistent nationwide sampling of standing dead trees, which may now supplant previous purely model-based approaches to standing dead biomass and C stock estimation. A substantial hurdle to estimating standing dead tree biomass and C attributes is that traditional estimation procedures are based on merchantability paradigms that may not reflect density reductions or structural loss due to decomposition common in standing dead trees. The goal of this study was to incorporate standing dead tree adjustments into the current estimation procedures and assess how biomass and C stocks change at multiple spatial scales.</p> <p>Results</p> <p>Accounting for decay and structural loss in standing dead trees significantly decreased tree- and plot-level C stock estimates (and subsequent C stocks) by decay class and tree component. At a regional scale, incorporating adjustment factors decreased standing dead quaking aspen biomass estimates by almost 50 percent in the Lake States and Douglas-fir estimates by more than 36 percent in the Pacific Northwest.</p> <p>Conclusions</p> <p>Substantial overestimates of standing dead tree biomass and C stocks occur when one does not account for density reductions or structural loss. Forest inventory estimation procedures that are descended from merchantability standards may need to be revised toward a more holistic approach to determining standing dead tree biomass and C attributes (i.e., attributes of tree biomass outside of sawlog portions). Incorporating density reductions and structural loss adjustments reduces uncertainty associated with standing dead tree biomass and C while improving consistency with field methods and documentation.</p