Harmonised statistics and maps of forest biomass and increment in Europe.

Forest biomass is an essential resource in relation to the green transition and its assessment is key for the sustainable management of forest resources. Here, we present a forest biomass dataset for Europe based on the best available inventory and satellite data, with a higher level of harmonisation and spatial resolution than other existing data. This database provides statistics and maps of the forest area, biomass stock and their share available for wood supply in the year 2020, and statistics on gross and net volume increment in 2010-2020, for 38 European countries. The statistics of most countries are available at a sub-national scale and are derived from National Forest Inventory data, harmonised using common reference definitions and estimation methodology, and updated to a common year using a modelling approach. For those counties without harmonised statistics, data were derived from the State of Europe's Forest 2020 Report at the national scale. The maps are coherent with the statistics and depict the spatial distribution of the forest variables at 100 m resolution.

A novel cubic-exp evaluation algorithm considering non-symmetrical axial response signals of confocal microscopes.

The depth discrimination in confocal microscopy is based on the digital analysis of depth response signals obtained by each camera pixel during measurement. Various signal-processing algorithms are used for this purpose. The accuracy of these algorithms is inter alia restricted by the axial symmetry of the signals. However, in practice response signals are rather asymmetrical especially in case of measurement objects with critical surface structures such as edges or steep flanks. We present a novel signal-processing algorithm based on an exponential function with a cubic argument to handle asymmetrical and also symmetrical depth response signals. Results obtained by this algorithm are compared to those of commonly used signal processing algorithms. It turns out that the novel algorithm is more robust, more accurate and exhibits a repeatability of a similar order compared to other algorithms. RESEARCH HIGHLIGHTS: A novel, more robust algorithm with improved accuracy in peak extraction especially for asymmetrical response signals in confocal microscopy is introduced and validated. Improved accuracy is demonstrated for height and layer thickness measurements.

Plot size matters: Toward comparable species richness estimates across plot-based inventories.

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.

MicroRNA Expression Differences in Blood-Derived CD19+ B Cells of Methotrexate Treated Rheumatoid Arthritis Patients.

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.

A century of National Forest Inventory in Norway - informing past, present, and future decisions.

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. Land-use changes are monitored and trees outside forests are recorded. FUTURE: Considerable research efforts towards the integration of remote sensing technologies enable the publication of the Norwegian Forest Resource Map since 2015, which is also used for small area estimation at the municipality level. On the analysis side, capacity and software for long term growth and yield prognosis are being developed. Furthermore, we foresee the inclusion of further variables for monitoring ecosystem services, and an increasing demand for mapped information. The relatively simple NFI design has proven to be a robust choice for satisfying steadily increasing information needs and concurrently providing consistent time series.

An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits.

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.

Systematic assessment of commercially available low-input miRNA library preparation kits.

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.

Mapping and estimating the total living biomass and carbon in low-biomass woodlands using Landsat 8 CDR data.

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.

Modeling Soil Carbon Dynamics in Northern Forests: Effects of Spatial and Temporal Aggregation of Climatic Input Data.

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.

Empirical coverage of model-based variance estimators for remote sensing assisted estimation of stand-level timber volume.

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