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

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

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

Evaluation of the potential index model to predict habitat suitability of forest species: the potential distribution of mountain pine (Pinus uncinata) in the Iberian peninsula

Characterization of the suitability or potentiality of a territory for forest tree species is an important source of information for forest planning and managing. In this study, we compared a relatively simple methodology to generate potential habitat distribution areas that has been traditionally used in Spain (the potential index model) with a statistical modelling approach (generalized linear model). We modelled the potential distribution of mountain pine (Pinus uncinata) in the Iberian peninsula as a working example. The potential index model generated a map of habitat suitability according to the values of an index of potentiality, whose distribution has usually divided into four categories based on quartiles (from optimum to low suitability). Considering all values of the index of potentiality as presences of mountain pine resulted in a low to moderate degree of agreement between the potential index model and the generalized linear model according to the kappa coefficient. Using the cut-off value of the index of potentiality that maximized the degree of agreement between both modelling approaches resulted in a substantial similarity between the maps of the predicted distribution of mountain pine. This cut-off value did lie in the upper-third quartile of the potential index distribution (high suitability category), and roughly coincided with the upper 30th percentile. The use of statistical techniques, which have proved to be powerful and versatile for species distribution modelling, is recommended. However, the potential index model, together with the adjustments proposed here, could be a reasonably simple methodology to predict the potential distribution of forest tree species that forest managers should take into account when evaluating forestation and afforestation projects

Temporal, spatial, and structural patterns of adult trembling aspen and white spruce mortality in Quebec's boreal forest

Temporal, spatial, and structural patterns of adult trembling aspen (Populus tremuloides Michx.) and white spruce (Picea glauca (Moench) Voss) mortality were studied in intact 150-year-old stands in the southwestern boreal forest of Quebec. For both species, mortality decreases (number of dead trees/total number of trees) with distance from the lake edge until 100-150 m, from which point it slightly increases. Strong peaks in mortality were found for 40- to 60-year-old aspen mainly between 1974 and 1992. Such mortality in relatively young aspen is likely related to competition for light from the dominant canopy trees. Also, the recruitment of this young aspen cohort is presumably the result of a stand breakup that occurred when the initial aspen-dominated stand was between 90 and 110 years old. For spruce, strong peaks in mortality were found in 110- to 150-year-old trees and they occurred mainly after 1980. No clear explanation could be found for these peaks, but we suggest that they may be related to senescence or weakening of the trees following the last spruce budworm outbreak. Suppressed and codominant aspen had a much higher mortality ratio than spruce in the same height class, while more surprisingly, no difference in mortality rate was found between dominant trees of the two species. Most spruce trees were found as standing dead, which leads us to reject the hypothesis that windthrow is an important cause of mortality for spruce in our forests

Crown Plasticity and Competition for Canopy Space: A New Spatially Implicit Model Parameterized for 250 North American Tree Species

BACKGROUND: Canopy structure, which can be defined as the sum of the sizes, shapes and relative placements of the tree crowns in a forest stand, is central to all aspects of forest ecology. But there is no accepted method for deriving canopy structure from the sizes, species and biomechanical properties of the individual trees in a stand. Any such method must capture the fact that trees are highly plastic in their growth, forming tessellating crown shapes that fill all or most of the canopy space. METHODOLOGY/PRINCIPAL FINDINGS: We introduce a new, simple and rapidly-implemented model--the Ideal Tree Distribution, ITD--with tree form (height allometry and crown shape), growth plasticity, and space-filling, at its core. The ITD predicts the canopy status (in or out of canopy), crown depth, and total and exposed crown area of the trees in a stand, given their species, sizes and potential crown shapes. We use maximum likelihood methods, in conjunction with data from over 100,000 trees taken from forests across the coterminous US, to estimate ITD model parameters for 250 North American tree species. With only two free parameters per species--one aggregate parameter to describe crown shape, and one parameter to set the so-called depth bias--the model captures between-species patterns in average canopy status, crown radius, and crown depth, and within-species means of these metrics vs stem diameter. The model also predicts much of the variation in these metrics for a tree of a given species and size, resulting solely from deterministic responses to variation in stand structure. CONCLUSIONS/SIGNIFICANCE: This new model, with parameters for US tree species, opens up new possibilities for understanding and modeling forest dynamics at local and regional scales, and may provide a new way to interpret remote sensing data of forest canopies, including LIDAR and aerial photography

Aboveground biomass density models for NASA's Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (similar to 25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available