Acquiring high spatial and temporal resolution imagery from small unmanned aerial systems (sUAS) provides new opportunities for inventorying forests at small scales. Only a few studies have investigated the use of UASs in forest inventories, and the results are inconsistent and incomplete. The present study used three-dimensional (3D) variables derived from UAS imagery in combination with ground reference data to fit linear models for Lorey’s mean height (h_L), dominant height (h_dom), stem number (N), basal area (G), and stem volume (V). Plot-level cross validation revealed adjusted R² values of 0.71, 0.97, 0.60, 0.60, and 0.85 for h_L, h_dom, N, G, and V, respectively, with corresponding RMSE values of 1.4 m, 0.7 m, 538.2 ha⁻¹, 4.5 m²∙ha⁻¹, and 38.3 m³∙ha⁻¹. The respective relative RMSE values were 13.3%, 3.5%, 39.2%, 15.4%, and 14.5% of the mean ground reference values. The mean predicted values did not differ significantly from the reference values. The results revealed that the use of UAS imagery can provide relatively accurate and timely forest inventory information at a local scale. In addition, the present study highlights the practical advantages of UAS-assisted forest inventories, including adaptive planning, high project customization, and rapid implementation, even under challenging weather conditions.

Despite the popularity of random forests (RF) as a prediction algorithm, methods for constructing confidence intervals for population means using this technique are still only sparsely reported. For two regional study areas (Spain and Norway) RF was used to predict forest volume or aboveground biomass using remotely sensed auxiliary data obtained from multiple sensors. Additionally, the changes per unit area of these forest attributes were estimated using indirect and direct methods. Multiple inferential frameworks have attracted increased recent attention for estimating the variances required for confidence intervals. For this study, three different statistical frameworks, design-based expansion, model-assisted and model-based estimators, were used for estimating population parameters and their variances. Pairs and wild bootstrapping approaches at different levels were compared for estimating the variances of the model-based estimates of the population means, as well as for mapping the uncertainty of the change predictions. The RF models accurately represented the relationship between the response and remotely sensed predictor variables, resulting in increased precision for estimates of the population means relative to design-based expansion estimates. Standard errors based on pairs bootstrapping within or internal to RF were considerably larger than standard errors based on both pairs and wild external bootstrapping of the entire RF algorithm. Pairs and wild external bootstrapping produced similar standard errors, but wild bootstrapping better mimicked the original structure of the sample data and better preserved the ranges of the predictor variables.

Recent developments in remote sensing (RS) technology have made several sources of auxiliary data available to support forest inventories. Thus, a pertinent question is how different sources of RS data should be combined with field data to make inventories cost-efficient. Hierarchical model-based estimation has been proposed as a promising way of combining: (i) wall-to-wall optical data that are only weakly correlated with forest structure; (ii) a discontinuous sample of active RS data that are more strongly correlated with structure; and (iii) a sparse sample of field data. Model predictions based on the strongly correlated RS data source are used for estimating a model linking the target quantity with weakly correlated wall-to-wall RS data. Basing the inference on the latter model, uncertainties due to both modeling steps must be accounted for to obtain reliable variance estimates of estimated population parameters, such as totals or means. Here, we generalize previously existing estimators for hierarchical model-based estimation to cases with non-homogeneous error variance and cases with correlated errors, for example due to clustered sample data. This is an important generalization to take into account data from practical surveys. We apply the new estimation framework to case studies that mimic the data that will be available from the Global Ecosystem Dynamics Investigation (GEDI) mission and compare the proposed estimation framework with alternative methods. Aboveground biomass was the variable of interest, Landsat data were available wall-to-wall, and sample RS data were obtained from an airborne LiDAR campaign that produced simulated GEDI waveforms. The results show that generalized hierarchical model-based estimation has potential to yield more precise estimates than approaches utilizing only one source of RS data, such as conventional model-based and hybrid inferential approaches.

∙ Key message: The study presents novel model-based estimators for growing stock volume and its uncertainty estimation, combining a sparse sample of field plots, a sample of laser data, and wall-to-wall Landsat data. On the basis of our detailed simulation, we show that when the uncertainty of estimating mean growing stock volume on the basis of an intermediate ALS model is not accounted for, the estimated variance of the estimator can be biased by as much as a factor of three or more, depending on the sample size at the various stages of the design. ∙ Context: This study concerns model-based inference for estimating growing stock volume in large-area forest inventories, combining wall-to-wall Landsat data, a sample of laser data, and a sparse subsample of field data. ∙ Aims: We develop and evaluate novel estimators and variance estimators for the population mean volume, taking into account the uncertainty in two model steps. ∙ Methods: Estimators and variance estimators were derived for two main methodological approaches and evaluated through Monte Carlo simulation. The first approach is known as two-stage least squares regression, where Landsat data were used to predict laser predictor variables, thus emulating the use of wall-to-wall laser data. In the second approach laser data were used to predict field-recorded volumes, which were subsequently used as response variables in modeling the relationship between Landsat and field data. Results: ∙ The estimators and variance estimators are shown to be at least approximately unbiased. Under certain assumptions the two methods provide identical results with regard to estimators and similar results with regard to estimated variances. ∙ Conclusion: We show that ignoring the uncertainty due to one of the models leads to substantial underestimation of the variance, when two models are involved in the estimation procedure.

A functional forest carbon measuring, reporting and verification (MRV) system to support climate change mitigation policies, such as REDD+, requires estimates of forest biomass carbon, as an input to estimate emissions. A combination of field inventory and remote sensing is expected to provide those data. By linking Landsat 8 and forest inventory data, we (1) developed linear mixed effects models for total living biomass (TLB) estimation as a function of spectral variables, (2) developed a 30 m resolution map of the total living carbon (TLC), and (3) estimated the total TLB stock of the study area. Inventory data consisted of tree measurements from 500 plots in 63 clusters in a 15,700 km2 study area, in miombo woodlands of Tanzania. The Landsat 8 data comprised two climate data record images covering the inventory area. We found a linear relationship between TLB and Landsat 8 derived spectral variables, and there was no clear evidence of spectral data saturation at higher biomass values. The root-mean-square error of the values predicted by the linear model linking the TLB and the normalized difference vegetation index (NDVI) is equal to 44 t/ha (49 % of the mean value). The estimated TLB for the study area was 140 Mt, with a mean TLB density of 81 t/ha, and a 95 % confidence interval of 74–88 t/ha. We mapped the distribution of TLC of the study area using the TLB model, where TLC was estimated at 47 % of TLB. The low biomass in the miombo woodlands, and the absence of a spectral data saturation problem suggested that Landsat 8 derived NDVI is suitable auxiliary information for carbon monitoring in the context of REDD+, for low-biomass, open-canopy woodlands.