Past: In the early twentieth century, forestry was one of the most important sectors in Norway and an agitated discussion about the perceived decline of forest resources due to over-exploitation was ongoing. To base the discussion on facts, the young state of Norway established Landsskogtakseringen – the world’s first National Forest Inventory (NFI). Field work started in 1919 and was carried out by county. Trees were recorded on 10 m wide strips with 1–5 km interspaces. Site quality and land cover categories were recorded along each strip. Results for the first county were published in 1920, and by 1930 most forests below the coniferous tree line were inventoried. The 2nd to 5th inventories followed in the years 1937–1986. As of 1954, temporary sample plot clusters on a 3 km × 3 km grid were used as sampling units. Present: The current NFI grid was implemented in the 6th NFI from 1986 to 1993, when permanent plots on a 3 km × 3 km grid were established below the coniferous tree line. As of the 7th inventory in 1994, the NFI is continuous, and 1/5 of the plots are measured annually. All trees with a diameter ≥ 5 cm are recorded on circular, 250 m² plots. The NFI grid was expanded in 2005 to cover alpine regions with 3 km × 9 km and 9 km × 9 km grids. In 2012, the NFI grid within forest reserves was doubled along the cardinal directions. Clustered temporary plots are used periodically to facilitate county-level estimates. As of today, more than 120 variables are recorded in the NFI including bilberry cover, drainage status, deadwood, and forest health. Land-use changes are monitored and trees outside forests are recorded. Future: Considerable research efforts towards the integration of remote sensing technologies enable the publication of the Norwegian Forest Resource Map since 2015, which is also used for small area estimation at the municipality level. On the analysis side, capacity and software for long term growth and yield prognosis are being developed. Furthermore, we foresee the inclusion of further variables for monitoring ecosystem services, and an increasing demand for mapped information. The relatively simple NFI design has proven to be a robust choice for satisfying steadily increasing information needs and concurrently providing consistent time series.

Laser scanning data from unmanned aerial vehicles (UAV-LS) offer new opportunities to estimate forest growing stock volume ( $V$ ) exclusively based on the UAV-LS data. We propose a method to measure tree attributes and using these measurements to estimate $V$ without the use of field data for calibration. The method consists of five steps: i) Using UAV-LS data, tree crowns are automatically identified and segmented wall-to-wall. ii) From all detected tree crowns, a sample is taken where diameter at breast height (DBH) can be recorded reliably as determined by visual assessment in the UAV-LS data. iii) Another sample of crowns is taken where tree species were identifiable from UAV image data. iv) DBH and tree species models are fit using the samples and applied to all detected tree crowns. v) Single tree volumes are predicted with existing allometric models using predicted species and DBH, and height directly obtained from UAV-LS. The method was applied to a Riegl-VUX data set with an average density of 1130 points m⁻² and 3 cm orthomosaic acquired over an 8.8 ha managed boreal forest. The volumes of the identified trees were aggregated to estimate plot-, stand-, and forest-level volumes which were validated using 58 independently measured field plots. The root-mean-square deviance ( $RMS{D}_{\%}$ ) decreased when increasing the spatial scale from the plot (32.2%) to stand (27.1%) and forest level (3.5%). The accuracy of the UAV-LS estimates varied given forest structure and was highest in open pine stands and lowest in dense birch or spruce stands. On the forest level, the estimates based on UAV-LS data were well within the 95% confidence interval of the intense field survey estimate, and both estimates had a similar precision. While the results are encouraging for further use of UAV-LS in the context of fully airborne forest inventories, future studies should confirm our findings in a variety of forest types and conditions.

The Norwegian forest resource map (SR16) maps forest attributes by combining national forest inventory (NFI), airborne laser scanning (ALS) and other remotely sensed data. While the ALS data were acquired over a time interval of 10 years using various sensors and settings, the NFI data are continuously collected. Aims of this study were to analyze the effects of stratification on models linking remotely sensed and field data, and assess the accuracy overall and at the ALS project level. The model dataset consisted of 9203 NFI field plots and data from 367 ALS projects, covering 17 Mha and 2/3 of the productive forest in Norway. Mixed-effects regression models were used to account for differences among ALS projects. Two types of stratification were used to fit models: 1) stratification by the three main tree species groups spruce, pine and deciduous resulted in species-specific models that can utilize a satellite-based species map for improving predictions, and 2) stratification by species and maturity class resulted in stratum-specific models that can be used in forest management inventories where each stand regularly is visually stratified accordingly. Stratified models were compared to general models that were fit without stratifying the data. The species-specific models had relative root-mean-squared errors (RMSEs) of 35%, 34%, 31%, and 12% for volume, aboveground biomass, basal area, and Lorey’s height, respectively. These RMSEs were 2–7 percentage points (pp) smaller than those of general models. When validating using predicted species, RMSEs were 0–4 pp. smaller than those of general models. Models stratified by main species and maturity class further improved RMSEs compared to species-specific models by up to 1.8 pp. Using mixed-effects models over ordinary least squares models resulted in a decrease of RMSE for timber volume of 1.0–3.9 pp., depending on the main tree species. RMSEs for timber volume ranged between 19%–59% among individual ALS projects. The stratification by tree species considerably improved models of forest structural variables. A further stratification by maturity class improved these models only moderately. The accuracy of the models utilized in SR16 were within the range reported from other ALS-based forest inventories, but local variations are apparent.

Due to the availability of good and reasonably priced auxiliary data, the use of model-based regression-synthetic estimators for small area estimation is popular in operational settings. Examples are forest management inventories, where a linking model is used in combination with airborne laser scanning data to estimate stand-level forest parameters where no or too few observations are collected within the stand. This paper focuses on different approaches to estimating the variances of those estimates. We compared a variance estimator which is based on the estimation of superpopulation parameters with variance estimators which are based on predictions of finite population values. One of the latter variance estimators considered the spatial autocorrelation of the residuals whereas the other one did not. The estimators were applied using timber volume on stand level as the variable of interest and photogrammetric image matching data as auxiliary information. Norwegian National Forest Inventory (NFI) data were used for model calibration and independent data clustered within stands were used for validation. The empirical coverage proportion (ECP) of confidence intervals (CIs) of the variance estimators which are based on predictions of finite population values was considerably higher than the ECP of the CI of the variance estimator which is based on the estimation of superpopulation parameters. The ECP further increased when considering the spatial autocorrelation of the residuals. The study also explores the link between confidence intervals that are based on variance estimates as well as the well-known confidence and prediction intervals of regression models.

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.

Background: The age of forest stands is critical information for forest management and conservation, for example for growth modelling, timing of management activities and harvesting, or decisions about protection areas. However, area-wide information about forest stand age often does not exist. In this study, we developed regression models for large-scale area-wide prediction of age in Norwegian forests. For model development we used more than 4800 plots of the Norwegian National Forest Inventory (NFI) distributed over Norway between latitudes 58° and 65° N in an 18.2 Mha study area. Predictor variables were based on airborne laser scanning (ALS), Sentinel-2, and existing public map data. We performed model validation on an independent data set consisting of 63 spruce stands with known age. Results: The best modelling strategy was to fit independent linear regression models to each observed site index (SI) level and using a SI prediction map in the application of the models. The most important predictor variable was an upper percentile of the ALS heights, and root mean squared errors (RMSEs) ranged between 3 and 31 years (6% to 26%) for SI-specific models, and 21 years (25%) on average. Mean deviance (MD) ranged between − 1 and 3 years. The models improved with increasing SI and the RMSEs were largest for low SI stands older than 100 years. Using a mapped SI, which is required for practical applications, RMSE and MD on plot level ranged from 19 to 56 years (29% to 53%), and 5 to 37 years (5% to 31%), respectively. For the validation stands, the RMSE and MD were 12 (22%) and 2 years (3%), respectively. Conclusions: Tree height estimated from airborne laser scanning and predicted site index were the most important variables in the models describing age. Overall, we obtained good results, especially for stands with high SI. The models could be considered for practical applications, although we see considerable potential for improvements if better SI maps were available.