Updating national forest inventory estimates of growing stock volume using hybrid inference

Abstract International organizations increasingly require estimates of forest parameters to monitor the state of and changes in forest resources, the sustainability of forest practices and the role of forests in the carbon cycle. Most countries rely on data from their national forest inventories (NFI) to produce these estimates. However, because NFI survey years may not match the required reporting years, techniques for updating NFI-based estimates are necessary. The main aim was to develop an unbiased method to update NFI estimates of mean growing stock volume (m 3 /ha) using models to predict annual plot-level volume change, and to estimate the associated uncertainties. Because the final large area volume estimates were based on plot-level model predictions rather than field observations, hybrid inference was necessary to accommodate both model prediction uncertainty and sampling variation. Specific objectives were to compare modelling approaches, to assess the utility of Landsat data for increasing model prediction accuracy, to select the most accurate method, and to compare model-based and design-based uncertainty components. For four monospecific forest types, data from the 2nd and 3rd Spanish NFI surveys together with site variables and Landsat images were used to construct models to predict NFI information for the year of the 4th NFI survey. Data from the 3rd and 4th surveys were used to assess the accuracy of the model predictions at both plot-level and large area spatial scales. The most accurate method used a set of three models: one to predict the probability of volume removals, one to predict the amount of removed volume, and one to predict gross annual volume. Incorporation of Landsat-based variables in the models increased prediction accuracy. Differences between large area estimates based on plot-level field observations for the 4th NFI survey and estimates based on the model predictions were minimal for all four forest types. Further, the standard errors of the estimates based on the model predictions were only slightly greater than standard errors based on the field observations. Thus, model predictions of plot-level growing stock volume based on field and satellite image data as auxiliary information can be used to update large area NFI estimates for reporting years for which spectral data are available but field observations are not. Finally, variances of means are under-estimated unless hybrid inferential methods are used to incorporate both model prediction uncertainty and sampling variation. For the two forest types for which the two sources of uncertainty were of the same order of magnitude, the under-estimation was non-negligible.

[1]  R. Hall,et al.  Modeling forest stand structure attributes using Landsat ETM+ data: Application to mapping of aboveground biomass and stand volume , 2006 .

[2]  Göran Ståhl,et al.  Hybrid estimators for mean aboveground carbon per unit area , 2016 .

[3]  Oscar García,et al.  Evaluating forest Growth Models , 1997 .

[4]  Sean P. Healey,et al.  Application of two regression-based methods to estimate the effects of partial harvest on forest structure using Landsat data , 2006 .

[5]  A. Lanz,et al.  The role of European National Forest Inventories for international forestry reporting , 2016, Annals of Forest Science.

[6]  Daniel E. Irwin,et al.  Estimating proportional change in forest cover as a continuous variable from multi-year MODIS data , 2008 .

[7]  Ronald E. McRoberts,et al.  Probability- and model-based approaches to inference for proportion forest using satellite imagery as ancillary data , 2010 .

[8]  J. Beauchamp,et al.  Corrections for Bias in Regression Estimates After Logarithmic Transformation , 1973 .

[9]  J. Ellis,et al.  Identifying and addressing gaps in the UNFCCC reporting framework , 2015 .

[10]  N. Coops Improvement in Predicting Stand Growth of Pinus radiafa (D. Don) across Landscapes Using NOAA AVHRR and Landsat MSS Imagery Combined with a Forest Growth Process Model (3-PGS) , 1999 .

[11]  Sarah C. Goslee,et al.  Analyzing Remote Sensing Data in R: The landsat Package , 2011 .

[12]  Warren B. Cohen,et al.  Spatial, spectral and temporal patterns of tropical forest cover change as observed with multiple scales of optical satellite data , 2007 .

[13]  S. Merino-de-Miguel,et al.  Forest biomass estimation through NDVI composites. The role of remotely sensed data to assess Spanish forests as carbon sinks , 2006 .

[14]  Gert-Jan Nabuurs,et al.  Future harvesting pressure on European forests , 2007, European Journal of Forest Research.

[15]  P. Teillet,et al.  On the Slope-Aspect Correction of Multispectral Scanner Data , 1982 .

[16]  David Riaño,et al.  Assessment of different topographic corrections in Landsat-TM data for mapping vegetation types (2003) , 2003, IEEE Trans. Geosci. Remote. Sens..

[17]  R. McRoberts,et al.  A general method for assessing the effects of uncertainty in individual-tree volume model predictions on large-area volume estimates with a subtropical forest illustration , 2015 .

[18]  Göran Ståhl,et al.  Use of models in large-area forest surveys: comparing model-assisted, model-based and hybrid estimation , 2016, Forest Ecosystems.

[19]  M. Schelhaas,et al.  Overview of methods and tools for evaluating future woody biomass availability in European countries , 2016, Annals of Forest Science.

[20]  Changhui Peng,et al.  Growth and yield models for uneven-aged stands: past, present and future , 2000 .

[21]  R. Nelson,et al.  Model-based inference for biomass estimation in a LiDAR sample survey in Hedmark County, Norway , 2011 .

[22]  F. Houllier,et al.  Forecasting wood resources on the basis of national forest inventory data. Application to Pinus pinaster Ait. in southwestern France , 2001 .

[23]  Piermaria Corona,et al.  Estimation of standing wood volume in forest compartments by exploiting airborne laser scanning information: model-based, design-based, and hybrid perspectives , 2014 .

[24]  T. Gschwantner,et al.  Annual volume increment of the European forests—description and evaluation of the national methods used , 2016, Annals of Forest Science.

[25]  W. Cohen,et al.  Landsat's Role in Ecological Applications of Remote Sensing , 2004 .

[26]  Joanne C. White,et al.  Integration of Landsat time series and field plots for forest productivity estimates in decision support models , 2016 .

[27]  S. T. Buckland,et al.  An Introduction to the Bootstrap. , 1994 .

[28]  Dengsheng Lu,et al.  Relationships between forest stand parameters and Landsat TM spectral responses in the Brazilian Amazon Basin , 2004 .

[29]  N. Nagelkerke,et al.  A note on a general definition of the coefficient of determination , 1991 .

[30]  C. Tucker Red and photographic infrared linear combinations for monitoring vegetation , 1979 .

[31]  Göran Ståhl,et al.  Adapting National Forest Inventories to changing requirements - the case of the Swedish National Forest Inventory at the turn of the 20th century , 2014 .

[32]  Göran Ståhl,et al.  Preparing emission reporting from forests: use of National Forest Inventories in European countries , 2008 .

[33]  David A. Ratkowsky,et al.  Nonlinear regression modeling : a unified practical approach , 1984 .

[34]  C. Woodcock,et al.  Classification and Change Detection Using Landsat TM Data: When and How to Correct Atmospheric Effects? , 2001 .

[35]  Terje Gobakken,et al.  Inference for lidar-assisted estimation of forest growing stock volume , 2013 .

[36]  Göran Ståhl,et al.  Estimating biomass in Hedmark County, Norway using national forest inventory field plots and airborne laser scanning , 2012 .

[37]  D. Leckie,et al.  Prediction Models for Estimating the Area, Volume, and Age of Rubber (Hevea brasiliensis) Plantations in Malaysia using Landsat TM Data , 2004 .

[38]  E. Tomppo National Forest Inventories : pathways for common reporting , 2010 .