Sammendrag

Key message We studied size distributions of decay-affected Norway spruce trees using cut-to-length harvester data. The harvester data comprised tree-level decay and decay severity recordings from 101 final felling stands, which enabled to analyze relationships between size distributions of all and decay-affected trees. Distribution matching technique was used to transfer the size distribution of all trees into the diameter at breast height (DBH) distribution of decay-affected trees. Context Stem decay of Norway spruce (Picea abies [L.] Karst.) results in large economic losses in timber production in the northern hemisphere. Forest management planning typically requires information on tree size distributions. However, size distributions of decay-affected trees generally remain unknown impeding decision-making in forest management planning. Aims Our aim was to analyze and model relationships between size distributions of all and decay-affected Norway spruce trees at the level of forest stands. Methods Cut-to-length harvester data of 93,456 trees were collected from 101 final felling stands in Norway. For each Norway spruce tree (94% of trees), the presence and severity of stem decay (incipient and advanced) were recorded. The stand-level size distributions (diameter at breast height, DBH; height, H) of all and decay-affected trees were described using the Weibull distribution. We proposed distribution matching (DM) models that transform either the DBH or H distribution of all trees into DBH distributions of decay-affected trees. We compared the predictive performance of DMs with a null-model that refers to a global Weibull distribution estimated based on DBHs of all harvested decay-affected trees. Results The harvester data showed that an average-sized decay-affected tree is larger and taller compared with an average-sized tree in a forest stand, while trees with advanced decay were generally shorter and thinner compared with trees having incipient decay. DBH distributions of decay-affected trees can be matched with smaller error index (EI) values using DBH (EI = 0.14) than H distributions (EI = 0.31). DM clearly outperformed the null model that resulted in an EI of 0.32. Conclusions The harvester data analysis showed a relationship between size distributions of all and decay-affected trees that can be explained by the spread biology of decay fungi and modeled using the DM technique. Keywords Root and butt rot, Heterobasidion spp., Armillaria spp., Cut-to-length harvester, Forest management and planning

Sammendrag

Accurate estimation of site productivity is essential for forest projections and scenario modelling. We present and evaluate models to predict site index (SI) and whether a site is productive (potential total stem volume production ≥ 1 m3·ha−1·year−1) in a wall-to-wall high-resolution (16 m × 16 m) SI map for Norway. We investigate whether remotely sensed data improve predictions. We also study the advantages and disadvantages of using boosted regression trees (BRT), a machine-learning algorithm, to create high-accuracy SI maps. We use climatic and topographical data, soil parent material, a land resource map, and depth to water, together with Sentinel-2 satellite images and airborne laser scanning metrics, as predictor variables. We use the SI observed at more than 10 000 National Forest Inventory (NFI) sample plots throughout Norway to fit BRT models and validate the models using 5822 independent temporary plots from the NFI. We benchmark our results against SI estimates from forest monitoring inventories. We find that the SI from BRT has root mean squared error (RMSE) ranging from 2.3 m (hardwoods) to 3.6 m (spruce) when tested against independent validation data from the NFI temporary plots. These RMSEs are similar or marginally better than an evaluation of SI estimates from operational forest management plans where SI normally stems from manual photo interpretation.

Til dokument

Sammendrag

Fine-grained information on the level of individual trees constitute key components for forest observation enabling forest management practices tackling the effects of climate change and the loss of biodiversity in forest ecosystems. Such information on individual tree crowns (ITC's) can be derived from the application of ITC segmentation approaches, which utilize remotely sensed data. However, many ITC segmentation approaches require prior knowledge about forest characteristics, which is difficult to obtain for parameterization. This can be avoided by the adoption of data-driven, automated workflows based on convolutional neural networks (CNN). To contribute to the advancements of efficient ITC segmentation approaches, we present a novel ITC segmentation approach based on the YOLOv5 CNN. We analyzed the performance of this approach on a comprehensive international unmanned aerial laser scanning (UAV-LS) dataset (ForInstance), which covers a wide range of forest types. The ForInstance dataset consists of 4192 individually annotated trees in high-density point clouds with point densities ranging from 498 to 9529 points m-2 collected across 80 sites. The original dataset was split into 70% for training and validation and 30% for model performance assessment (test data). For the best performing model, we observed a F1-score of 0.74 for ITC segmentation and a tree detection rate (DET %) of 64% in the test data. This model outperformed an ITC segmentation approach, which requires prior knowledge about forest characteristics, by 41% and 33% for F1-score and DET %, respectively. Furthermore, we tested the effects of reduced point densities (498, 50 and 10 points per m-2) on ITC segmentation performance. The YOLO model exhibited promising F1-scores of 0.69 and 0.62 even at point densities of 50 and 10 points m-2, respectively, which were between 27% and 34% better than the ITC approach that requires prior knowledge. Furthermore, the areas of ITC segments resulting from the application of the best performing YOLO model were close to the reference areas (RMSE = 3.19 m-2), suggesting that the YOLO-derived ITC segments can be used to derive information on ITC level.

Til dokument

Sammendrag

Wood resources have been essential for human welfare throughout history. Also nowadays, the volume of growing stock (GS) is considered one of the most important forest attributes monitored by National Forest Inventories (NFIs) to inform policy decisions and forest management planning. The origins of forest inventories closely relate to times of early wood shortage in Europe causing the need to explore and plan the utilisation of GS in the catchment areas of mines, saltworks and settlements. Over time, forest surveys became more detailed and their scope turned to larger areas, although they were still conceived as stand-wise inventories. In the 1920s, the first sample-based NFIs were introduced in the northern European countries. Since the earliest beginnings, GS monitoring approaches have considerably evolved. Current NFI methods differ due to country-specific conditions, inventory traditions, and information needs. Consequently, GS estimates were lacking international comparability and were therefore subject to recent harmonisation efforts to meet the increasing demand for consistent forest resource information at European level. As primary large-area monitoring programmes in most European countries, NFIs assess a multitude of variables, describing various aspects of sustainable forest management, including for example wood supply, carbon sequestration, and biodiversity. Many of these contemporary subject matters involve considerations about GS and its changes, at different geographic levels and time frames from past to future developments according to scenario simulations. Due to its historical, continued and currently increasing importance, we provide an up-to-date review focussing on large-area GS monitoring where we i) describe the origins and historical development of European NFIs, ii) address the terminology and present GS definitions of NFIs, iii) summarise the current methods of 23 European NFIs including sampling methods, tree measurements, volume models, estimators, uncertainty components, and the use of air- and space-borne data sources, iv) present the recent progress in NFI harmonisation in Europe, and v) provide an outlook under changing climate and forest-based bioeconomy objectives.

Til dokument

Sammendrag

Key message: Using satellite-based maps, Ceccherini et al. (Nature 583:72-77, 2020) report abruptly increasing harvested area estimates in several EU countries beginning in 2015. Using more than 120,000 National Forest Inventory observations to analyze the satellite-based map, we show that it is not harvested area but the map’s ability to detect harvested areas that abruptly increases after 2015 in Finland and Sweden. Keywords: Global Forest Watch, Landsat, Remote sensing, National Forest Inventory, Greenhouse Gas Inventory

Til dokument

Sammendrag

Wheel ruts, i.e. soil deformations caused by harvesting machines, are considered a negative environmental impact of forest operations and should be avoided or ameliorated. However, the mapping of wheel ruts that would be required to monitor harvesting operations and to plan amelioration measures is a tedious and time-consuming task. Here, we examined whether a combination of drone imagery and algorithms from the field of artificial intelligence can automate the mapping of wheel ruts. We used a deep-learning image-segmentation method (ResNet50 + UNet architecture) that was trained on drone imagery acquired shortly after harvests in Norway, where more than 160 km of wheel ruts were manually digitized. The cross-validation of the model based on 20 harvested sites resulted in F1 scores of 0.69–0.84 with an average of 0.77, and in total, 79 per cent of wheel ruts were correctly detected. The highest accuracy was obtained for severe wheel ruts (average user’s accuracy (UA) = 76 per cent), and the lowest accuracy was obtained for light wheel ruts (average UA = 67 per cent). Considering the nowadays ubiquitous availability of drones, the approach presented in our study has the potential to greatly increase the ability to effectively map and monitor the environmental impact of final felling operations with respect to wheel ruts. The automated mapping of wheel ruts may serve as an important input to soil impact analyses and thereby support measures to restore soil damages.

Sammendrag

Skogen i Norge har et årlig netto opptak av CO2 tilsvarende nær halvparten av de nasjonale menneskeskapte utslippene. Skogens bidrag i klimasammenheng kan økes gjennom økt opptak av CO2 i skog, men også ved økt lagring av karbon i treprodukter. Treproduktene (harvested wood products – HWP) som årlig rapporteres i Norges klimagassregnskap (National Inventory Report - NIR) for arealbrukssektoren (Land Use, Land-Use Change and Forestry - LULUCF) inkluderer trelast, trebaserte plater og papir- og kartongprodukter. Skogens årlige netto opptak av CO2 utgjorde i 2019 23,6 millioner tonn CO2 ekvivalenter. Årlig tilførsel samme år til lagring i treprodukter utgjorde 2% av dette (449 kt CO2). Totallageret av karbon i treprodukter i Norge i 2019 tilsvarer 109,1 millioner tonn CO2. Lagring av karbon i treproduker er et av virkemidlene for at Norge skal oppfylle sine klimamål under Parisavtalen. Med andre ord, en økning i årlig lagring av karbon i treprodukter vil bidra til å oppfylle Norges forpliktelser. Økt bruk av tre vil også kunne bidra til å redusere utslipp i andre sektorer gjennom at treprodukter kan erstatte materialer med høyere klimagasspåvirkning (substitusjon). Økt bruk av tre vil gjenspeiles i klimagassregnskapet for treprodukter, men den fulle effekten av substitusjonen vil ikke gjenspeiles i dette regnskapet. Målet med rapporten er å kvantifisere hvor stor andel av årlig hogst som rapporteres inn i treprodukter i Norges klimagassregnskap fra 1961 og fram til i dag. Økt forståelse av hvordan verdikjedens utnyttelse gjenspeiles i klimagassregnskapet er en nødvendig forutsetning for å bidra til økt fremtidig lagring av karbon i treprodukter. For å bedre forstå årsakene til variasjonene i rapporterte treprodukter mellom år, beskrives også den årlige materialflyten av alle typer treprodukter etter hogst (råmaterialer, halvfabrikata og bioenergi) basert på de årlige volumene (1961-2019) av: 1) produksjon (total nasjonal produksjon), 2) produksjon ekskl. eksport (nasjonalt forbruk), 3) eksport og 4) import.....

Til dokument

Sammendrag

To understand the state and trends in biodiversity beyond the scope of monitoring programs, biodiversity indicators must be comparable across inventories. Species richness (SR) is one of the most widely used biodiversity indicators. However, as SR increases with the size of the area sampled, inventories using different plot sizes are hardly comparable. This study aims at producing a methodological framework that enables SR comparisons across plot-based inventories with differing plot sizes. We used National Forest Inventory (NFI) data from Norway, Slovakia, Spain, and Switzerland to build sample-based rarefaction curves by randomly incrementally aggregating plots, representing the relationship between SR and sampled area. As aggregated plots can be far apart and subject to different environmental conditions, we estimated the amount of environmental heterogeneity (EH) introduced in the aggregation process. By correcting for this EH, we produced adjusted rarefaction curves mimicking the sampling of environmentally homogeneous forest stands, thus reducing the effect of plot size and enabling reliable SR comparisons between inventories. Models were built using the Conway–Maxell–Poisson distribution to account for the underdispersed SR data. Our method successfully corrected for the EH introduced during the aggregation process in all countries, with better performances in Norway and Switzerland. We further found that SR comparisons across countries based on the country-specific NFI plot sizes are misleading, and that our approach offers an opportunity to harmonize pan-European SR monitoring. Our method provides reliable and comparable SR estimates for inventories that use different plot sizes. Our approach can be applied to any plot-based inventory and count data other than SR, thus allowing a more comprehensive assessment of biodiversity across various scales and ecosystems.

Sammendrag

The number of people affected by snow avalanches during recreational activities has increased over the recent years. An instrument to reduce these numbers are improved terrain classification systems. One such system is the Avalanche Terrain Exposure Scale (ATES). Forests can provide some protection from avalanches, and information on forest attributes can be incorporated into avalanche hazard models such as the automated ATES model (AutoATES). The objectives of this study were to (i) map forest stem density and canopy-cover based on National Forest Inventory and remote sensing data and, (ii) use these forest attributes as input to the AutoATES model. We predicted stem density and directly calculated canopy-cover in a 20 Mha study area in Norway. The forest attributes were mapped for 16 m × 16 m pixels, which were used as input for the AutoATES model. The uncertainties of the stem number and canopy-cover maps were 30% and 31%, respectively. The overall classification accuracy of 52 ski-touring routes in Western Norway with a total length of 282 km increased from 55% in the model without forest information to 67% when utilizing canopy cover. The F1 score for the three predicted ATES classes improved by 31%, 9%, and 6%.

Sammendrag

This study aimed at estimating total forest above-ground net change (ΔAGB; Gg) over five years (2014–2019) based on model-assisted estimation utilizing freely available satellite imagery. The study was conducted for a boreal forest area (approx. 1.4 Mha) in Norway where bi-temporal national forest inventory (NFI), Sentinel-2, and Landsat data were available. Biomass change was modelled based on a direct approach. The precision of estimates using only the NFI data in a basic expansion estimator was compared to four different alternative model-assisted estimates using 1) Sentinel-2 or Landsat data, and 2) using bi- or uni-temporal remotely sensed data. We found that spaceborne optical data improved the precision of the purely field-based estimates by a factor of up to three. The most precise estimates were found for the model-assisted estimation using bi-temporal Sentinel-2 (standard error; SE = 1.7 Gg). However, the decrease in precision when using Landsat data was small (SE = 1.92 Gg). We also found that ΔAGB could be precisely estimated when remotely sensed data were available only at the end of the monitoring period. We conclude that satellite optical data can considerably improve ΔAGB estimates, when repeated and coincident field data are available. The free availability, global coverage, frequent update, and long-term time horizon make data from programs such as Sentinel-2 and Landsat a valuable data source for consistent and durable monitoring of forest carbon dynamics.

Til dokument

Sammendrag

Policy measures and management decisions aimed at enhancing the role of forests in mitigating climate change require reliable estimates of carbon (C)-stock dynamics in greenhouse gas inventories (GHGIs). The aim of this study was to assemble design-based estimators to provide estimates relevant for GHGIs using National Forest Inventory (NFI) data. We improve basic expansion (BE) estimators of living-biomass C-stock loss using only field data, by leveraging with remote sensing auxiliary data in model-assisted (MA) estimators. Our case studies from Norway, Sweden, Denmark, and Latvia covered an area of >70 Mha. Landsat-based forest cover loss (FCL) and one-time wall-to-wall airborne laser scanning (ALS) served as auxiliary data. ALS provided information on the C stock before a potential disturbance indicated by FCL. The use of FCL in MA estimators resulted in considerable efficiency gains, which in most cases were further increased by adding ALS. A doubling of efficiency was possible for national estimates and even larger efficiencies were observed at the subnational level. Average annual estimates were considerably more precise than pooled estimates of the NFI data from all years at once. The combination of remotely sensed and NFI field data yields reliable estimators, which is not necessarily the case when using remotely sensed data without reference observations.

Sammendrag

Background 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. Materials and methods 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. Results 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. Conclusions 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.

Til dokument

Sammendrag

Rheumatoid arthritis (RA) is a complex disease with a wide range of underlying susceptibility factors. Recently, dysregulation of microRNAs (miRNAs) in RA have been reported in several immune cell types from blood. However, B cells have not been studied in detail yet. Given the autoimmune nature of RA with the presence of autoantibodies, CD19+ B cells are a key cell type in RA pathogenesis and alterations in CD19+ B cell subpopulations have been observed in patient blood. Therefore, we aimed to reveal the global miRNA repertoire and to analyze miRNA expression profile differences in homogenous RA patient phenotypes in blood-derived CD19+ B cells. Small RNA sequencing was performed on CD19+ B cells of newly diagnosed untreated RA patients (n=10), successfully methotrexate (MTX) treated RA patients in remission (MTX treated RA patients, n=18) and healthy controls (n=9). The majority of miRNAs was detected across all phenotypes. However, significant expression differences between MTX treated RA patients and controls were observed for 27 miRNAs, while no significant differences were seen between the newly diagnosed patients and controls. Several of the differentially expressed miRNAs were previously found to be dysregulated in RA including miR-223-3p, miR-486-3p and miR-23a-3p. MiRNA target enrichment analysis, using the differentially expressed miRNAs and miRNA-target interactions from miRTarBase as input, revealed enriched target genes known to play important roles in B cell activation, differentiation and B cell receptor signaling, such as STAT3, PRDM1 and PTEN. Interestingly, many of those genes showed a high degree of correlated expression in CD19+ B cells in contrast to other immune cell types. Our results suggest important regulatory functions of miRNAs in blood-derived CD19+ B cells of MTX treated RA patients and motivate for future studies investigating the interactive mechanisms between miRNA and gene targets, as well as the possible predictive power of miRNAs for RA treatment response.

Sammendrag

Diameter at breast height (DBH) distributions offer valuable information for operational and strategic forest management decisions. We predicted DBH distributions using Norwegian national forest inventory and airborne laser scanning data and compared the predictive performances of linear mixed-effects (PPM), generalized linear-mixed (GLM), and k nearest-neighbor (NN) models. While GLM resulted in smaller prediction errors than PPM, both were clearly outperformed by NN. We therefore studied the ability of the NN model to improve the precision of stem frequency estimates by DBH classes in the 8.7 Mha study area using a model-assisted (MA) estimator suitable for systematic sampling. MA estimates yielded greater than or approximately equal efficiencies as direct estimates using field data only. The relative efficiencies (REs) associated with the MA estimates ranged between 0.95–1.47 and 0.96–1.67 for 2 and 6 cm DBH class widths, respectively, when dominant tree species were assumed to be known. The use of a predicted tree species map, instead of the observed information, decreased the REs by up to 10%.

Sammendrag

Butt rot (BR) damage of a tree results from a decay caused by a pathogenic fungus. BR damages associated with Norway spruce (Picea abies [L.] Karst.) account for considerable economic losses in timber production across the northern hemisphere. While information on BR damages is critical for optimal decision-making in forest management, maps of BR damages are typically lacking in forest information systems. Timber volume damaged by BR was predicted at the stand-level in Norway using harvester information of 186,026 stems (clear-cuts), remotely sensed, and environmental data (e.g. climate and terrain characteristics). This study utilized Random Forests models with two sets of predictor variables: (1) predictor variables available after harvest (theoretical case) and (2) predictor variables available prior to harvest (mapping case). Our findings showed that forest attributes characterizing the maturity of forest, such as remote sensing-based height, harvested timber volume and quadratic mean diameter at breast height, were among the most important predictor variables. Remotely sensed predictor variables obtained from airborne laser scanning data and Sentinel-2 imagery were more important than the environmental variables. The theoretical case with a leave-stand-out cross-validation resulted in an RMSE of 11.4 m3 · ha−1 (pseudo-R2: 0.66) whereas the mapping case resulted in a pseudo-R2 of 0.60. When spatially distinct clusters of harvested forest stands were used as units in the cross-validation, the RMSE value and pseudo-R2 associated with the mapping case were 15.6 m3 · ha−1 and 0.37, respectively. The findings associated with the different cross-validation schemes indicated that the knowledge about the BR status of spatially close stands is of high importance for obtaining satisfactory error rates in the mapping of BR damages.

Sammendrag

Key message Large-scale forest resource maps based on national forest inventory (NFI) data and airborne laser scanning may facilitate synergies between NFIs and forest management inventories (FMIs). A comparison of models used in such a NFI-based map and a FMI indicate that NFI-based maps can directly be used in FMIs to estimate timber volume of mature spruce forests. Context Traditionally, FMIs and NFIs have been separate activities. The increasing availability of detailed NFI-based forest resource maps provides the possibility to eliminate or reduce the need of field sample plot measurements in FMIs if their accuracy is similar. Aims We aim to (1) compare a timber volume model used in a NFI-based map and models used in a FMI, and (2) evaluate utilizing additional local sample plots in the model of the NFI-based map. Methods Accuracies of timber volume estimates using models from an existing NFI-based map and a FMI were compared at plot and stand level. Results Estimates from the NFI-based map were similar to or more accurate than the FMI. The addition of local plots to the modeling data did not clearly improve the model of the NFI-based map. Conclusion The comparison indicates that NFI-based maps can directly be used in FMIs for timber volume estimation in mature spruce stands, leading to potentially large cost savings.

Sammendrag

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 m2 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. Landuse 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.

Til dokument

Sammendrag

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Til dokument

Sammendrag

Background: Large area forest inventories often use regular grids (with a single random start) of sample locations to ensure a uniform sampling intensity across the space of the surveyed populations. A design-unbiased estimator of variance does not exist for this design. Oftentimes, a quasi-default estimator applicable to simple random sampling (SRS) is used, even if it carries with it the likely risk of overestimating the variance by a practically important margin. To better exploit the precision of systematic sampling we assess the performance of five estimators of variance, including the quasi default. In this study, simulated systematic sampling was applied to artificial populations with contrasting covariance structures and with or without linear trends. We compared the results obtained with the SRS, Matérn’s, successive difference replication, Ripley’s, and D’Orazio’s variance estimators. Results: The variances obtained with the four alternatives to the SRS estimator of variance were strongly correlated, and in all study settings consistently closer to the target design variance than the estimator for SRS. The latter always produced the greatest overestimation. In populations with a near zero spatial autocorrelation, all estimators, performed equally, and delivered estimates close to the actual design variance. Conclusion: Without a linear trend, the SDR and DOR estimators were best with variance estimates more narrowly distributed around the benchmark; yet in terms of the least average absolute deviation, Matérn’s estimator held a narrow lead. With a strong or moderate linear trend, Matérn’s estimator is choice. In large populations, and a low sampling intensity, the performance of the investigated estimators becomes more similar. Keywords: Spatial autocorrelation, Linear trend, Model based, Design biased, Matérn variance, Successive difference replication variance, Geary contiguity coefficient, Random site effects

Sammendrag

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−2 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 ( RMSD% ) 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.

Sammendrag

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.

Til dokument

Sammendrag

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t ha−1 ). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

Sammendrag

Nation-wide Sentinel-2 mosaics were used with National Forest Inventory (NFI) plot data for modelling and subsequent mapping of spruce-, pine-, and deciduous-dominated forest in Norway at a 16 m × 16 m resolution. The accuracies of the best model ranged between 74% for spruce and 87% for deciduous forest. An overall accuracy of 90% was found on stand level using independent data from more than 42 000 stands. Errors mostly resulting from a forest mask reduced the model accuracies by ∼10%. The produced map was subsequently used to generate model-assisted (MA) and poststratified (PS) estimates of species-specific forest area. At the national level, efficiencies of the estimates increased by 20% to 50% for MA and up to 90% for PS. Greater minimum numbers of observations constrained the use of PS. For MA estimates of municipalities, efficiencies improved by up to a factor of 8 but were sometimes also less than 1. PS estimates were always equally as or more precise than direct and MA estimates but were applicable in fewer municipalities. The tree species prediction map is part of the Norwegian forest resource map and is used, among others, to improve maps of other variables of interest such as timber volume and biomass.

Til dokument

Sammendrag

Forest inventories provide predictions of stand means on a routine basis from models with auxiliary variables from remote sensing as predictors and response variables from field data. Many forest inventory sampling designs do not afford a direct estimation of the among-stand variance. As consequence, the confidence interval for a model-based prediction of a stand mean is typically too narrow. We propose a new method to compute (from empirical regression residuals) an among-stand variance under sample designs that stratify sample selections by an auxiliary variable, but otherwise do not allow a direct estimation of this variance. We test the method in simulated sampling from a complex artificial population with an age class structure. Two sampling designs are used (one-per-stratum, and quasi systematic), neither recognize stands. Among-stand estimates of variance obtained with the proposed method underestimated the actual variance by 30-50%, yet 95% confidence intervals for a stand mean achieved a coverage that was either slightly better or at par with the coverage achieved with empirical linear best unbiased estimates obtained under less efficient two-stage designs.

Til dokument

Sammendrag

Global Forest Watch (GFW) provides a global map of annual forest cover loss (FCL) produced from Landsat imagery, offering a potentially powerful tool for monitoring changes in forest cover. In managed forests, FCL primarily provides information on commercial harvesting. A semi-autonomous method for providing data on the location and attributes of harvested sites at a landscape level was developed which could significantly improve the basis for catchment management, including risk mitigation. FCL in combination with aerial images was used for detecting and characterising harvested sites in a 1607 km2 mountainous boreal forest catchment in south-central Norway. Firstly, the forest cover loss map was enhanced (FCLE) by removing small isolated forest cover loss patches that had a high probability of representing commission errors. The FCLE map was then used to locate and assess sites representing annual harvesting activity over a 17-year period. Despite an overall accuracy of >98%, a kappa of 0.66 suggested only a moderate quality for detecting harvested sites. While errors of commission were negligible, errors of omission were more considerable and at least partially attributed to the presence of residual seed trees on the site after harvesting. The systematic analysis of harvested sites against aerial images showed a detection rate of 94%, but the area of the individual harvested site was underestimated by 29% on average. None of the site attributes tested, including slope, area, altitude, or site shape index, had any effect on the accuracy of the area estimate. The annual harvest estimate was 0.6% (standard error 12%) of the productive forest area. On average, 96% of the harvest was carried out on flat to moderately steep terrain (<40% slope), 3% on steep terrain (40% to 60% slope), and 1% on very steep terrain (>60% slope). The mean area of FCLE within each slope category was 1.7 ha, 0.9 ha, and 0.5 ha, respectively. The mean FCLE area increased from 1.0 ha to 3.2 ha on flat to moderate terrain over the studied period, while the frequency of harvesting increased from 249 to 495 sites per year. On the steep terrain, 35% of the harvesting was done with cable yarding, and 62% with harvester-forwarder systems. On the very steep terrain (>60% slope), 88% of the area was harvested using cable yarding technology while harvesters and forwarders were used on 12% of the area. Overall, FCL proved to be a useful dataset for the purpose of assessing harvesting activity under the given conditions.

Sammendrag

This paper describes the development and utility of the Norwegian forest resources map (SR16). SR16 is developed using photogrammetric point cloud data with ground plots from the Norwegian National Forest Inventory (NFI). First, an existing forest mask was updated with object-based image analysis methods. Evaluation against the NFI forest definitions showed Cohen's kappa of 0.80 and accuracy of 0.91 in the lowlands and a kappa of 0.73 and an accuracy of 0.96 in the mountains. Within the updated forest mask, a 16×16 m raster map was developed with Lorey's height, volume, biomass, and tree species as attributes (SR16-raster). All attributes were predicted with generalized linear models that explained about 70% of the observed variation and had relative RMSEs of about 50%. SR16-raster was segmented into stand-like polygons that are relatively homogenous in respect to tree species, volume, site index, and Lorey's height (SR16-vector). When SR16 was utilized in a combination with the NFI plots and a model-assisted estimator, the precision was on average 2–3 times higher than estimates based on field data only. In conclusion, SR16 is useful for improved estimates from the Norwegian NFI at various scales. The mapped products may be useful as additional information in Forest Management Inventories.

Sammendrag

I denne rapporten presenteres framskrivninger for opptak og utslipp fra arealbrukssektoren (eng. Land Use, Land-Use Change and Forestry; LULUCF) frem til 2100. Framskrivninger av opptak og utslipp av CO2 og andre klimagasser fra arealbrukssektoren er utført i tråd med metodikken brukt i klimagassregnskapet for Norge i 2019 (Miljødirektoratet mfl. 2019), og basert på data rapportert for 2010 – 2017 som referanseperiode. Framskrivningen for opptak og utslipp i skog er basert på tilsvarende metodikk som i referansebanen for forvaltede skogarealer (eng. Forest Reference Level, FRL), som publisert i National Forest Accounting Plan (Klima- og miljødepartementet 2019), men basert på nyeste tilgjengelige data og med implementert politikk. Framskrivningene er utført basert på rapporteringen under FNs klimakonvensjon og Kyotoprotokollen, samt EUs LULUCF-forordning.

Til dokument

Sammendrag

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

Sammendrag

Monitoring changes in forest height, biomass and carbon stock is important for understanding the drivers of forest change, clarifying the geography and magnitude of the fluxes of the global carbon budget and for providing input data to REDD+. The objective of this study was to investigate the feasibility of covering these monitoring needs using InSAR DEM changes over time and associated estimates of forest biomass change and corresponding net CO2 emissions. A wall-to-wall map of net forest change for Uganda with its tropical forests was derived from two Digital Elevation Model (DEM) datasets, namely the SRTM acquired in 2000 and TanDEM-X acquired around 2012 based on Interferometric SAR (InSAR) and based on the height of the phase center. Errors in the form of bias, as well as parallel lines and belts having a certain height shift in the SRTM DEM were removed, and the penetration difference between X- and C-band SAR into the forest canopy was corrected. On average, we estimated X-band InSAR height to decrease by 7 cm during the period 2000–2012, corresponding to an estimated annual CO2 emission of 5 Mt for the entirety of Uganda. The uncertainty of this estimate given as a 95% confidence interval was 2.9–7.1 Mt. The presented method has a number of issues that require further research, including the particular SRTM biases and artifact errors; the penetration difference between the X- and C-band; the final height adjustment; and the validity of a linear conversion from InSAR height change to AGB change. However, the results corresponded well to other datasets on forest change and AGB stocks, concerning both their geographical variation and their aggregated values.

Sammendrag

Background: Generalized height-diameter curves based on a re-parameterized version of the Korf function for Norway spruce (Picea abies (L.) Karst.), Scots pine (Pinus sylvestris L.) and silver birch (Betula pendula Roth) in Norway are presented. The Norwegian National Forest Inventory (NFI) is used as data base for estimating the model parameters. The derived models are developed to enable spatially explicit and site sensitive tree height imputation in forest inventories as well as future tree height predictions in growth and yield scenario simulations. Methods: Generalized additive mixed models (gamm) are employed to detect and quantify potentially non-linear effects of predictor variables. In doing so the quadratic mean diameter serves as longitudinal covariate since stand age, as measured in the NFI, shows only a weak correlation with a stands developmental status in Norwegian forests. Additionally the models can be locally calibrated by predicting random effects if measured height-diameter pairs are available. Based on the model selection of non-constraint models, shape constraint additive models (scam) were fit to incorporate expert knowledge and intrinsic relationships by enforcing certain effect patterns like monotonicity. Results: Model comparisons demonstrate that the shape constraints lead to only marginal differences in statistical characteristics but ensure reasonable model predictions. Under constant constraints the developed models predict increasing tree heights with decreasing altitude, increasing soil depth and increasing competition pressure of a tree. A two-dimensional spatially structured effect of UTM-coordinates accounts for the potential effects of large scale spatially correlated covariates, which were not at our disposal. The main result of modelling the spatially structured effect is lower tree height prediction for coastal sites and with increasing latitude. The quadratic mean diameter affects both the level and the slope of the height-diameter curve and both effects are positive. Conclusions: In this investigation it is assumed that model effects in additive modelling of height-diameter curves which are unfeasible and too wiggly from an expert point of view are a result of quantitatively or qualitatively limited data bases. However, this problem can be regarded not to be specific to our investigation but more general since growth and yield data that are balanced over the whole data range with respect to all combinations of predictor variables are exceptional cases. Hence, scam may provide methodological improvements in several applications by combining the flexibility of additive models with expert knowledge.

Til dokument

Sammendrag

The Nordic countries have long traditions in forest inventory and remote sensing (RS). In sample-based national forest inventories (NFIs), utilization of aerial photographs started during the 1960s, satellite images during the 1980s, laser scanning during the 2000s, and photogrammetric point clouds during the 2010s. In forest management inventories (FMI), utilization of aerial photos started during the 1940s and laser scanning during the 2000s. However, so far, RS has mostly been used for map production and research rather than for estimation of regional parameters or inference on their accuracy. In recent years, the RS technology has been developing very fast. At the same time, the needs for information are constantly increasing. New technologies have created possibilities for cost-efficient production of accurate, large area forest data sets, which also will change the way forest inventories are done in the future. In this study, we analyse the state-of-the-art both in the NFIs and FMIs in the Nordic countries. We identify the benefits and drawbacks of different RS materials and data acquisition approaches with different user perspectives. Based on the analysis, we identify the needs for further development and emerging research questions. We also discuss alternatives for ownership of the data and cost-sharing between different actors in the field.

Sammendrag

Denne rapporten sammenstiller informasjon om tilstand og utviklingstendenser for foryngelse i skog, med utgangspunkt i data fra Resultatkartleggingen, Landsskogtakseringen og Økonomisystem for Skogordningene (ØKS). Oppdraget har videre omfattet å vurdere styrker og svakheter ved dagens systemer for innhenting av tilstandsdata om foryngelsesaktiviteten i Norge, og foreslå metoder og eventuelle verktøy som vil gi bedre oversikt over foryngelsestilstanden på årlig basis...

Sammendrag

In many applications, estimates are required for small sub-populations with so few (or no) sample plots that direct estimators that do not utilize auxiliary variables (e.g. remotely sensed data) are not applicable or result in low precision. This problem is overcome in small area estimation (SAE) by linking the variable of interest to auxiliary variables using a model. Two types of models can be distinguished based on the scale on which they operate: i) Unit-level models are applied in the well-known area-based approach (ABA) and are commonly used in forest inventories supported by fine-resolution 3D remote sensing data such as airborne laser scanning (ALS) or digital aerial photogrammetry (AP); ii) Area-level models, where the response is a direct estimate based on a sample within the domain and the explanatory variables are aggregated auxiliary variables, are less frequently applied. Estimators associated with these two model types can make use of sample plots within domains if available and reduce to so-called synthetic estimators in domains where no sample plots are available. We used both model types and their associated model-based estimators in the same study area with AP data as auxiliary variables. Heteroscedasticity, i.e. for continuous dependent variables typically an increasing dispersion of re- siduals with increasing predictions, is often observed in models linking field- and remotely sensed data. This violates the model assumption that the distribution of the residual errors is constant. Complying with model assumptions is required for model-based methods to result in reliable estimates. Addressing heteroscedasticity in models had considerable impacts on standard errors. When complying with model assumptions, the precision of estimates based on unit-level models was, on average, considerably greater (29%–31% smaller standard errors) than those based on area-level models. Area-level models may nonetheless be attractive because they allow the use of sampling designs that do not easily link to remotely sensed data, such as variable radius plots.

Sammendrag

I henhold til det Norske klimagassregnskap leder avskoging til en betydelig del av de nasjonale klimagassutslipp. Målet med denne rapporten er en kartlegging av størrelse og årsaker til avskoging som kan forbedre forståelsen av avskogingsprosesser, og på sikt kan være et første steg for å redusere utslippene fra avskogingen. I Kyotoprotokollen er avskoging betegnet som menneskeskapte endringer fra skog til en annen arealkategori siden 1990. I Norge har avskoging siden 1990 vært på om lag 58 km2 per år. På grunn av påskoging (på aktivt forvaltede arealer) og skogutvidelse (naturlig etablering på ikke forvaltede arealer) har skogarealet ikke forandret seg nevneverdig. Men den teoretiske produksjonsevnen, altså skogens evne til å produsere biomasse og dermed også til å ta opp karbon fra atmosfæren i et gitt tidsrom, av det samlede arealet av påskoging og skogutvidelse er mindre enn produksjonsevnen av avskogingsarealet. Hovedgrunnen til avskoging var utbygging (68 % av avskogingsarealet), men også omlegging til beite (18 %) eller nydyrking (13 %) bidro. I denne rapporten er alle areal og utslippsestimater basert på Landsskogtakseringen som er en landsdekkende utvalgsundersøkelse. Grunnet det lille totale areal av avskoging i Norge er arealestimatene assosiert med relativ stor usikkerhet relatert til antall prøvefelter i utvalgskartleggingen. Blant utbyggingskategoriene var vei og bebyggelse de viktigste grunnene til avskoging. Traktor- og skogsbilveier var de største enkeltkategoriene blant veikategoriene og til sammen står de for om lag 13 % av avskogingsarealet. Bolig og fritidsbolig var de største enkeltkategoriene blant bebyggelseskategoriene og til sammen står de for om lag 13 % av avskogingsarealet...

Til dokument

Sammendrag

Small-area estimation is a subject area of growing importance in forest inventories. Modelling the link between a study variable Y and auxiliary variables X— in pursuit of an improved accuracy in estimators—is typically done at the level of a sampling unit. However, for various reasons, it may only be possible to formulate a linking model at the level of an area of interest (AOI). Area-level models and their potential have rarely been explored in forestry. This study demonstrates, with data (Y = stem volume per ha) from four actual inventories aided by aerial laser scanner data (3 cases) or photogrammetric point clouds (1 case), application of three distinct models representing the currency of area-level modelling. The studied AOIs varied in size from forest management units to forest districts, and municipalities. The variance explained by X declined sharply with the average size of an AOI. In comparison with a direct estimate mean of Y in an AOI, all three models achieved practically important reduction in the relative root-mean-squared error of an AOI mean. In terms of the reduction in mean-squared errors, a model with a spatial location effect was overall most attractive. We recommend the pursuit of a spatial model component in area-level modelling as promising within the context of a forest inventory.

Til dokument

Sammendrag

Many parties to the United Nation's Framework Convention on Climate Change (UNFCCC) base their reporting of change in Land Use, Land-Use Change and Forestry (LULUCF) sector carbon pools on national forest inventories. A strong feature of sample-based inventories is that very detailed measurements can be made at the level of plots. Uncertainty regarding the results stems primarily from the fact that only a sample, and not the entire population, is measured. However, tree biomass on sample plots is not directly measured but rather estimated using regression models based on allometric features such as tree diameter and height. Estimators of model parameters are random variables that exhibit different values depending on which sample is used for estimating model parameters. Although sampling error is strongly influenced by the sample size when the model is applied, modeling error is strongly influenced by the sample size when the model is under development. Thus, there is a trade-off between which sample sizes to use when applying and developing models. This trade-off has not been studied before and is of specific interest for countries developing new national forest inventories and biomass models in the REDD+ context. This study considers a specific sample design and population. This fact should be considered when extrapolating results to other locations and populations.

Sammendrag

The use of digital aerial photogrammetry (DAP) for forest inventory purposes has been widely studied and can produce comparable accuracy compared with airborne laser scanning (ALS) in small, homogeneous areas. However, the accuracy of DAP for large scale applications with heterogeneous terrain and forest vegetation has not yet been reported. In this study we examined the accuracy of timber volume, biomass and basal area prediction models based on DAP and national forest inventory (NFI) data on a large area in central Norway. Two separate point clouds were derived from aerial image acquisitions of 2010 and 2013. Vegetation heights were extracted by subtracting terrain elevation derived from ALS. A large number of NFI sample plots (483) measured between 2010 and 2014 were used as reference data to fit linear models for timber volume, biomass and basal area with height metrics derived from the DAP data as explanatory variables. Variables describing the heterogeneous environmental and image acquisition conditions were calculated and their influence on the model accuracy was tested. The results showed that forest parameter prediction using DAP works well when applied to a large area. The model fits of the timber volume, biomass and basal area models were good with R2 of 0.80, 0.81, 0.81 and RMSEs of 41.43 m3 ha−1 (55% of the mean observed value), 32.49 t ha−1 (47%), 5.19 m2 ha−1 (41%), respectively. Only a small proportion of the variation could be attributed to the heterogeneous conditions. The inclusion of the relative sun inclination led to an improvement of the model RMSEs by 2% of the mean observed values. The relatively low cost and stability across large areas make DAP an attractive source of auxiliary information for large scale forest inventories.

Til dokument

Sammendrag

This study presents an approach for predicting stand-level forest attributes utilizing mobile laser scanning data collected as a nonprobability sample. Firstly, recordings of stem density were made at point locations every 10th metre along a subjectively chosen mobile laser scanning track in a forest stand. Secondly, kriging was applied to predict stem density values for the centre point of all grid cells ina5m×5m lattice across the stand. Thirdly, due to nondetectability issues, a correction term was computed based on distance sampling theory. Lastly, the mean stem density at stand level was predicted as the mean of the point-level predictions multiplied with the correction factor, and the corresponding variance was estimated. Many factors contribute to the uncertainty of the stand-level prediction; in the variance estimator, we accounted for the uncertainties due to kriging prediction and due to estimating a detectability model from the laser scanning data. The results from our new approach were found to correspond fairly well to estimates obtained using field measurements from an independent set of 54 circular sample plots. The predicted number of stems in the stand based on the proposed methodology was 1366 with a 12.9% relative standard error. The corresponding estimate based on the field plots was 1677 with a 7.5% relative standard error.

Til dokument

Sammendrag

Forest stands are important units of management. A stand-by-stand estimation of the mean and variance of an attribute of interest (Y) remains a priority in forest enterprise inventories. The advent of powerful and cost effective remotely sensed auxiliary variables (X) correlated with Y means that a census of X in the forest enterprise is increasingly available. In combination with a probability sample of Y, the census affords a modeldependent stand-level inference. It is important, however, that the sampling design affords an estimation of possible stand-effects in the model linking X to Y.We demonstrate, with simulated data, that failing to quantify non-zero stand-effects in the intercept of a linear population-level model can lead to a serious underestimation of the uncertainty in a model-dependent estimate of a stand mean, and by extension a confidence interval with poor coverage.We also provide an approximation to the variance of stand-effects in an intercept for the case when a sampling design does not afford estimation. Furthermore, we propose a method to correct a potential negative bias in an estimate of the variance of stand-effects when a sampling design prescribes few stands with small within-stand sample sizes.

Sammendrag

Skog er en viktig del av den globale karbonsyklusen, både som lager og som opptaker av karbon fra atmosfæren. Norge rapporterer årlig utslipp og opptak av klimagasser i skog til FNs klimakonvensjon, samt til Kyotoprotokollen. Skog rapporteres under landsektoren (Land Use, Land-Use Change and Forestry; LULUCF). I 2015 var netto-opptaket i skog 29,0 millioner tonn CO2-ekvivalenter, mens det totale utslippet av klimagasser i Norge i de øvrige sektorene var 53,9 millioner tonn. Netto opptak i skog tilsvarer dermed 54 prosent av klimagassutslippene i de øvrige sektorene.

Til dokument

Sammendrag

Estimates of stand averages are needed by forest management for planning purposes. In forest enterprise inventories supported by remotely sensed auxiliary data, these estimates are typically derived exclusively from a model that does not consider stand effects in the study variable. Variance estimators for these means may seriously underestimate uncertainty, and confidence intervals may be too narrow when a model used for computing a stand mean omits a nontrivial stand effect in one or more of the model parameters, a nontrivial spatial distance dependent autocorrelation in the model residuals, or both. In simulated sampling from 36 populations with stands of different sizes and differing with respect to (i) the correlation between a study variable (Y) and two auxiliary variables (X), (ii) the magnitude of stand effects in the intercept of a linear population model linking X to Y, and (iii) a first-order autoregression in Y and X, we learned that none of the tested designs provided reliable estimates of the within-stand autocorrelation among model residuals. More-reliable estimates were possible from stand-wide predictions of Y. The anticipated bias in an estimated autoregression parameter had a modest influence on estimates of variance and coverage of nominal 95% confidence intervals for a synthetic stand mean.

Sammendrag

App og webside basert på R Shiny for å gi en alminnelig og fleksibel tilgang til Landsskogtakseringens estimater. Språk: norsk og engelsk. Flere oppdateringer i 2016, 2017 og 2018. Online: https://landsskog.nibio.no/

Til dokument

Sammendrag

Background: 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. Results: 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. Conclusion: 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.

Til dokument

Sammendrag

Boreal forests contain 30% of the global forest carbon with the majority residing in soils. While challenging to quantify, soil carbon changes comprise a significant, and potentially increasing, part of the terrestrial carbon cycle. Thus, their estimation is important when designing forest-based climate change mitigation strategies and soil carbon change estimates are required for the reporting of greenhouse gas emissions. Organic matter decomposition varies with climate in complex nonlinear ways, rendering data aggregation nontrivial. Here, we explored the effects of temporal and spatial aggregation of climatic and litter input data on regional estimates of soil organic carbon stocks and changes for upland forests. We used the soil carbon and decomposition model Yasso07 with input from the Norwegian National Forest Inventory (11275 plots, 1960–2012). Estimates were produced at three spatial and three temporal scales. Results showed that a national level average soil carbon stock estimate varied by 10% depending on the applied spatial and temporal scale of aggregation. Higher stocks were found when applying plot-level input compared to country-level input and when long-term climate was used as compared to annual or 5-year mean values. A national level estimate for soil carbon change was similar across spatial scales, but was considerably (60–70%) lower when applying annual or 5-year mean climate compared to long-term mean climate reflecting the recent climatic changes in Norway. This was particularly evident for the forest-dominated districts in the southeastern and central parts of Norway and in the far north. We concluded that the sensitivity of model estimates to spatial aggregation will depend on the region of interest. Further, that using long-term climate averages during periods with strong climatic trends results in large differences in soil carbon estimates. The largest differences in this study were observed in central and northern regions with strongly increasing temperatures.

Til dokument

Sammendrag

This paper focuses on the use of models for increasing the precision of estimators in large-area forest surveys. It is motivated by the increasing availability of remotely sensed data, which facilitates the development of models predicting the variables of interest in forest surveys. We present, review and compare three different estimation frameworks where models play a core role: model-assisted, model-based, and hybrid estimation. The first two are well known, whereas the third has only recently been introduced in forest surveys. Hybrid inference mixes design-based and model-based inference, since it relies on a probability sample of auxiliary data and a model predicting the target variable from the auxiliary data..We review studies on large-area forest surveys based on model-assisted, model-based, and hybrid estimation, and discuss advantages and disadvantages of the approaches. We conclude that no general recommendations can be made about whether model-assisted, model-based, or hybrid estimation should be preferred. The choice depends on the objective of the survey and the possibilities to acquire appropriate field and remotely sensed data. We also conclude that modelling approaches can only be successfully applied for estimating target variables such as growing stock volume or biomass, which are adequately related to commonly available remotely sensed data, and thus purely field based surveys remain important for several important forest parameters.

Til dokument

Sammendrag

Modelling stem taper and volume is crucial in many forest management and planning systems. Taper models are used for diameter prediction at any location along the stem of a sample tree. Furthermore, taper models are flexible means to provide information on the stem volume and assortment structure of a forest stand or other management units. Usually, taper functions are mean functions of multiple linear or nonlinear regression models with diameter at breast height and tree height as predictor variables. In large-scale inventories, an upper diameter is often considered as an additional predictor variable to improve the reliability of taper and volume predictions. Most studies on stem taper focus on accurately modelling the mean function; the error structure of the regression model is neglected or treated as secondary. We present a semi-parametric linear mixed model where the population mean diameter at an arbitrary stem location is a smooth function of relative height. Observed tree-individual diameter deviations from the population mean are assumed to be realizations of a smooth Gaussian process with the covariance depending on the sampled diameter locations. In addition to the smooth random deviation from the population average, we consider independent zero mean residual errors in order to describe the deviations of the observed diameter measurements from the tree-individual smooth stem taper. The smooth model components are approximated by cubic spline functions with a B-spline basis and a small number of knots. The B-spline coefficients of the population mean function are treated as fixed effects, whereas coefficients of the smooth tree-individual deviation are modelled as random effects with zero mean and a symmetric positive definite covariance matrix. The taper of a tree is predicted using an arbitrary number of diameter and corresponding height measurements at arbitrary positions along the stem to calibrate the tree-individual random deviation from the population mean estimated by the fixed effects. This allows a flexible application of the method in practice. Volume predictions are calculated as the integral over cross-sectional areas estimated from the calibrated taper curve. Approximate estimators for the mean squared errors of volume estimates are provided. If the tree height is estimated or measured with error, we use the “law of total expectation and variance” to derive approximate diameter and volume predictions with associated confidence and prediction intervals. All methods presented in this study are implemented in the R-package TapeR.

Til dokument

Sammendrag

Bark beetles cause widespread damages in the coniferous-dominated forests of central Europe and North America. In the future, areas affected by bark beetles may further increase due to climate change. However, the early detection of the bark beetle green attack can guide management decisions to prevent larger damages. For this reason, a field-based bark beetle monitoring program is currently implemented in Germany. The combination of remote sensing and field data may help minimizing the reaction time and reducing costs of monitoring programs covering large forested areas. In this case study, RapidEye and TerraSAR-X data were analyzed separately and in combination to detect bark beetle green attack. The remote sensing data were acquired in May 2009 for a study site in south-west Germany. In order to distinguish healthy areas and areas affected by bark beetle green attack, three statistical approaches were compared: generalized linear models (GLM), maximum entropy (ME) and random forest (RF). The spatial scale (minimum mapping unit) was 78.5 m2. TerraSAR-X data resulted in fair classification accuracy with a cross-validated Cohen’s Kappa Coefficient (kappa) of 0.23. RapidEye data resulted in moderate classification accuracy with a kappa of 0.51. The highest classification accuracy was obtained by combining the TerraSAR-X and RapidEye data, resulting in a kappa of 0.74. The accuracy of ME models was considerably higher than the accuracy of GLM and RF models.

Sammendrag

There is a need for monitoring methods for forest volume, biomass and carbon based on satellite remote sensing. In the present study we tested interferometric X-band SAR (InSAR) from the Tandem-X mission. The aim of the study was to describe how accurate volume and biomass could be estimated from InSAR height and test whether the relationships were curvilinear or not. The study area was a spruce dominated forest in southeast Norway. We selected 28 stands in which we established 192 circular sample plots of 250 m2, accurately positioned by a Differential Global Positioning System (dGPS). Plot level data on stem volume and aboveground biomass were derived from field inventory. Stem volume ranged fromzero to 596 m3/ha, and aboveground biomass up to 338 t/ha.We generated 2 Digital Surface Models (DSMs) fromInSAR processing of two co-registered, HH-polarized TanDEM-X image pairs – one ascending and one descending pair.We used a Digital TerrainModel (DTM) from airborne laser scanning (ALS) as a reference and derived a 10 m × 10 m Canopy Height Model (CHM), or InSAR height model. We assigned each plot to the nearest 10 m × 10 m InSAR height pixel. We applied a nonlinear, mixed model for the volume and biomass modeling, and from a full model we removed effects with a backward stepwise approach. InSAR heightwas proportional to volume and aboveground biomass, where a 1 m increase in InSAR height corresponded to a volume increase of 23 m3/ha and a biomass increase of 14 t/ha. Root Mean Square Error (RMSE) values were 43–44% at the plot level and 19–20% at the stand level.

Til dokument

Sammendrag

National Forest Inventories (NFIs) provide estimates of forest parameters for national and regional scales. Many key variables of interest, such as biomass and timber volume, cannot be measured directly in the field. Instead, models are used to predict those variables from measurements of other field variables. Therefore, the uncertainty or variability of NFI estimates results not only from selecting a sample of the population but also from uncertainties in the models used to predict the variables of interest. The aim of this study was to quantify the model-related variability of Norway spruce (Picea abies [L.] Karst) biomass stock and change estimates for the Norwegian NFI. The model-related variability of the estimates stems from uncertainty in parameter estimates of biomass models as well as residual variability and was quantified using a Monte Carlo simulation technique. Uncertainties in model parameter estimates, which are often not available for published biomass models, had considerable influence on the model-related variability of biomass stock and change estimates. The assumption that the residual variability is larger than documented for the models and the correlation of within-plot model residuals influenced the model-related variability of biomass stock change estimates much more than estimates of the biomass stock. The larger influence on the stock change resulted from the large influence of harvests on the stock change, although harvests were observed rarely on the NFI sample plots in the 5-year period that was considered. In addition, the temporal correlation between model residuals due to changes in the allometry had considerable influence on the model-related variability of the biomass stock change estimate. The allometry may, however, be assumed to be rather stable over a 5-year period. Because the effects of model-related variability of the biomass stock and change estimates were much smaller than those of the sampling-related variability, efforts to increase the precision of estimates should focus on reducing the sampling variability. If the model-related variability is to be decreased, the focus should be on the tree fractions of living branches as well as stump and roots.

Sammendrag

Landsskogtakseringen gir i dag estimater av skoglige parametre på nasjonalt og regionalt nivå ved hjelp av et systematisk nettverk av prøveflater. Den genererte informasjonen blir brukt til en rekke formål, herunder utforming av nasjonal og regional skogpolitikk, rapportering til internasjonale organer og avtaler slik som Kyoto-protokollen, strategisk planlegging for skognæringen, samt overvåking av viktige biologiske indikatorer.

Sammendrag

K-nearest neighbor (kNN) approaches are popular statistical methods for predicting forest attributes in airborne laser scanning (ALS) based inventories. Their main upsides are the simplicity to predict multivariate response variables and their freeness of distributional assumptions on the conditional response.One of their largest draw-backs is that predictions outside the range of the reference data inherently result in an under- or overestimation. This property of kNN approaches is known as extrapolation bias and aggravates with an increasing number of neighbors (k) used for the prediction.This study presents one possibility to reduce extrapolation biases of predictions based on the area-based approach (ABA) by using individual tree crown (ITC) approaches within those specific areas of a low density ALS acquisition where the point density might be sufficiently high for using ITC methods.In the proposed strategy, additional (or artificial) reference plots augmented field measured plots. Artificial plots were created by applying ITC segmentation to a canopy height model derived from high density ALS data. The response variable biomass per hectare was predicted for every segment following a semi-ITC approach.The segment predictions were aggregated on the artificial plot level. The artificial plots were then treated in the same way as the original reference data to make predictions in areas with low density ALS data based on the ABA. It was hereby assumed that the predicted plot level response on the artificial plots is equivalent with the observed plot level response on the original reference data.The data consisted of 110 reference plots with a smaller data range than the 201 independent validation plots. Considerable extrapolation bias was visible if only the reference plots were used for the prediction. Almost no extrapolation bias was found if the prediction was based on reference plots augmented by artificial plots. The root mean squared error (RMSE) of the biomass predictions based on the reference plots was 39.1%. The RMSE reduced to 29.8% if the reference plots were augmented by artificial plots.

Sammendrag

The Norwegian National Forest Inventory (NNFI) provides estimates of forest parameters on national and regional scales by means of a systematic network of permanent sample plots. One of the biggest challenges for the NNFI is the interest in forest attribute information for small sub-populations such as municipalities or protected areas. Frequently, too few sampled observations are available for such small areas to allow estimates with acceptable precision. However, if an auxiliary variable exists that is correlated with the variable of interest, small area estimation (SAE) techniques may provide means to improve the precision of estimates. The study aimed at estimating the mean above-ground forest biomass for small areas with high precision and accuracy, using SAE techniques. For this purpose, the simple random sampling (SRS) estimator, the generalized regression (GREG) estimator, and the unit-level empirical best linear unbiased prediction (EBLUP) estimator were compared. Mean canopy height obtained from a photogrammetric canopy height model (CHM) was the auxiliary variable available for every population element. The small areas were 14 municipalities within a 2,184 km2 study area for which an estimate of the mean forest biomass was sought. The municipalities were between 31 and 527 km2 and contained 1–35 NNFI sample plots located within forest. The mean canopy height obtained from the CHM was found to have a strong linear correlation with forest biomass. Both the SRS estimator and the GREG estimator result in unstable estimates if they are based on too few observations. Although this is not the case for the EBLUP estimator, the estimators were only compared for municipalities with more than five sample plots. The SRS resulted in the highest standard errors in all municipalities. Whereas the GREG and EBLUP standard errors were similar for small areas with many sample plots, the EBLUP standard error was usually smaller than the GREG standard error. The difference between the EBLUP and GREG standard error increased with a decreasing number of sample plots within the small area. The EBLUP estimates of mean forest biomass within the municipalities ranged between 95.01 and 153.76 Mg ha−1, with standard errors between 8.20 and 12.84 Mg ha−1.

Til dokument

Sammendrag

Climate change is a factor that largely contributes to the increase of forest areas affected by natural damages. Therefore, the development of methodologies for forest monitoring and rapid assessment of affected areas is required. Space-borne synthetic aperture radar (SAR) imagery with high resolution is now available for large-scale forest mapping and forest monitoring applications. However, a correct interpretation of SAR images requires an adequate preprocessing of the data consisting of orthorectification and radiometric calibration. The resolution and quality of the digital elevation model (DEM) used as reference is crucial for this purpose. Therefore, the primary aim of this study was to analyze the influence of the DEM quality used in the preprocessing of the SAR data on the mapping accuracy of forest types. In order to examine TerraSAR-X images to map forest dominated by deciduous and coniferous trees, High Resolution SpotLight images were acquired for two study sites in southern Germany. The SAR images were preprocessed with a Shuttle Radar Topography Mission (SRTM) DEM (resolution approximately 90 m), an airborne laser scanning (ALS) digital terrain model (DTM) (5 m resolution), and an ALS digital surface model (DSM) (5 m resolution). The orthorectification of the SAR images using high resolution ALS DEMs was found to be important for the reduction of errors in pixel location and to increase the classification accuracy of forest types. SAR images preprocessed with ALS DTMs resulted in the highest classification accuracies, with kappa coefficients of 0.49 and 0.41, respectively. SAR images preprocessed with ALS DTMs resulted in greater accuracy than those preprocessed with ALS DSMs in most cases. The classification accuracy of forest types using SAR images preprocessed with the SRTM DEM was fair, with kappa coefficients of 0.23 and 0.32, respectively.Analysis of the radar backscatter indicated that sample plots dominated by coniferous trees tended to have lower scattering coefficients than plots dominated by deciduous trees. Leaf-off images were only slightly better suited for the classification than leaf-on images. The combination of leaf-off and leaf-on improved the classification accuracy considerably since the backscatter changed between seasons, especially in deciduous-dominated forest.

Sammendrag

Digital aerial images over Vestfold county were acquired by TerraTec in summer 2007 with a Vexcel UltraCamX sensor. The flying height above-ground was approximately 2800-3000 m which resulted in images of approximately 1880x2880 m size. The images were acquired in north-south oriented flight strips with a 20% side and 60% within-strip overlap. Panchromatic image data were acquired in 20 cm ground sampling distance (GSD). Near infrared, red, green and blue image bands were acquired in 60 cm GSD but were pansharpened to a 20 cm pixel size by the data vendor. The original radiometric resolution of the images (12 bit) was resampled to 8 bit for archival storage. The plane location and orientation during image acquisition were logged using a GPS and an inertial navigation system (INS). To increase the accuracy of the external orientation, an aerial triangulation was performed based on 34 ground control points using the software Match-AT.

Sammendrag

Vegetation height information is one of the most important variables for predicting forest attributes such as timber volume and biomass. Although airborne laser scanning (ALS) data are operationally used in forest planning inventories in Norway, a regularly repeated acquisition of ALS data for large regions has yet to be realized. Therefore, several research groups analyze the use of other data sources to retrieve vegetation height information. One very promising approach is the photogrammetric derivation of vegetation heights from overlapping digital aerial images. Aerial images are acquired over almost all European countries on a regular basis making image data readily available. The Norwegian Forest and Landscape Institute (Skog og Landskap) invited researchers and practitioners that produce and utilize photogrammetric data to share their experiences. More than 30 participants followed the invitation and contributed to a successful event with interesting presentations and discussions. We wish to thank the speakers for their contributions and hope that all participants found the seminar useful. These short proceedings of the seminar include summaries of the talks. The presentations, which provide more information, can be found at the end of this document.

Sammendrag

Top dieback and mortality of Norway spruce is a particular forest damage that has severe occurrences in scattered forest stands in southeast Norway. As a part of a project to study the extent and causes of the damage we are working on an algorithm for automatic detection dead and declining spruce trees for an entire county, - Vestfold. The data set is aerial imagery. The county was covered in 2007. Preliminary tests showed a considerable confusion between dead trees and bare ground. In order to avoid this confusion we have had the imagery automatically processed into a photogrammetric digital surface model (DSM) and true orthophotos. The data set derived from this processing was a 5 layer file, containing blue, green, red, and near-infrared, as well as the height above ground of the canopy height model (a DSM normalized by the terrain height, nDSM)

Sammendrag

The aim in the analysis of sample surveys is frequently to derive estimates of sub-population characteristics. This task is denoted small area estimation (SAE). Often, the sample available for the sub-population is, however, too small to allow a reliable estimate. Frequently, auxiliary variables exist that are correlated with the variable of interest. Several estimators can make use of auxiliary information which may reduce the variance of the estimate.

Sammendrag

The Norwegian National Forest Inventory (NNFI) provides estimates of forest parameters on national and regional scales by means of a systematic network of permanent sample plots. One of the biggest challenges for the NNFI is the interest in forest attribute information for small subpopulations such as municipalities or protected areas. Frequently, too few sampled observations are available for those small areas to allow an estimate with acceptable precision. However, if an auxiliary variable exists that is correlated with the variable of interest, small area estimation (SAE) techniques may provide means to improve the precision of estimates.

Til dokument

Sammendrag

While forest inventories based on airborne laser scanning data (ALS) using the area based approach (ABA) have reached operational status, methods using the individual tree crown approach (ITC) have basically remained a research issue. One of the main obstacles for operational applications of ITC is biased results often experienced due to segmentation errors. In this article, we propose a new method, called "semi-ITC" that overcomes the main problems related to ITC by imputing ground truth data within crown segments from the nearest neighboring segment. This may be none, one, or several trees. The distances between segments were derived based on a set of explanatory variables using two nonparametric methods, i.e., most similar neighbor inference (MSN) and random forest (RF). RF favored the imputation of common observations in the data set which resulted in significant biases. Main conclusions are therefore based on MSN. The explanatory variables were calculated by means of small footprint ALS and multispectral data. When testing with empirical data the new method compared favorably to the well-known ABA. Another advantage of the new method over the ABA is that it allowed for the modeling of rare tree species. The results of predicting timber volume with the semi-ITC method were unbiased and the root mean squared error (RMSE) on plot level was smaller than the standard deviation of the observed response variables. The relative RMSEs after cross validation using semi-ITC for total volume and volume of the individual species pine, spruce, birch, and aspen on plot level were 17, 38, 40, 101, and 222%, respectively. Due to the unbiasedness of the estimation, this study is a showcase for how to use crown segments resulting from ITC algorithms in a forest inventory context. (C) 2009 Elsevier Inc. All rights reserved.

Sammendrag

The semi-individual tree crown approach (semi-ITC) was used to predict crown base heights (CBH) on the level of single crown segments based on airborne laser scanning (ALS) derived metrics. The root-mean-squared-differences (RMSD) on the segment level were smallest for spruce. However, they were larger than the standard deviation of the measured CBH for pine and birch. The RMSD values were also larger compared to other studies. This can in part be explained by the fact that the semi-ITC approach incorporates errors of the segmentation algorithm. As a consequence, all instead of only correctly identified trees were considered in modeling which results in more realistic RMSD values. After aggregating the individual segment predictions to the plot level, the RMSD values were smaller than the standard deviations of the field measurements and comparable to other studies. The relative RMSD values for birch, spruce, pine and all species were 51.61, 35.22, 49.28, and 13.89%, respectively.