Global biomass maps can increase the precision of (sub)national aboveground biomass estimates: A comparison across tropical countries.

Countries within the tropics face ongoing challenges in completing or updating their national forest inventories (NFIs), critical for estimating aboveground biomass (AGB) and for forest-related greenhouse gas (GHG) accounting. While previous studies have explored the integration of map information with local reference data to fill in data gaps, limited attention has been given to the specific challenges presented by the clustered plot designs frequently employed by NFIs when combined with remote sensing-based biomass map units. This research addresses these complexities by conducting four country case-studies, encompassing a variety of NFI characteristics within a range of AGB densities. Examining four country case-studies (Peru, Guyana, Tanzania, Mozambique), we assess the potential of European Space Agency's Climate Change Initiative (CCI) global biomass maps to increase precision in (sub)national AGB estimates. We compare a baseline approach using NFI field-based data with a model-assisted scenario incorporating a locally calibrated CCI biomass map as auxiliary information. The original CCI biomass maps systematically underestimate AGB in three of the four countries at both the country and stratum level, with particularly weak agreement at finer map resolution. However, after calibration with country-specific NFI data, stratum and country-level AGB estimates from the model-assisted scenario align well with those obtained solely from field-based data and official country reports. Introducing maps as a source of auxiliary information fairly increased the precision of stratum and country-wise AGB estimates, offering greater confidence in estimating AGB for GHG reporting purposes. Considering the challenges tropical countries face with implementing their NFIs, it is sensible to explore the potential benefits of biomass maps for climate change reporting mechanisms across biomes. While country-specific NFI design assumptions guided our model-assisted inference strategies, this study also uncovers transferable insights from the application of global biomass maps with NFI data, providing valuable lessons for climate research and policy communities.

Spatial predictions and uncertainties of forest carbon fluxes for carbon accounting.

Countries have pledged to different national and international environmental agreements, most prominently the climate change mitigation targets of the Paris Agreement. Accounting for carbon stocks and flows (fluxes) is essential for countries that have recently adopted the United Nations System of Environmental-Economic Accounting - ecosystem accounting framework (UNSEEA) as a global statistical standard. In this paper, we analyze how spatial carbon fluxes can be used in support of the UNSEEA carbon accounts in five case countries with available in-situ data. Using global multi-date biomass map products and other remotely sensed data, we mapped the 2010-2018 carbon fluxes in Brazil, the Netherlands, the Philippines, Sweden and the USA using National Forest Inventory (NFI) and local biomass maps from airborne LiDAR as reference data. We identified areas that are unsupported by the reference data within environmental feature space (6-47% of vegetated country area); cross-validated an ensemble machine learning (RMSE=9-39 Mg C [Formula: see text] and [Formula: see text]=0.16-0.71) used to map carbon fluxes with prediction intervals; and assessed spatially correlated residuals (<5 km) before aggregating carbon fluxes from 1-ha pixels to UNSEEA forest classes. The resulting carbon accounting tables revealed the net carbon sequestration in natural broadleaved forests. Both in plantations and in other woody vegetation ecosystems, emissions exceeded sequestration. Overall, our estimates align with FAO-Forest Resource Assessment and national studies with the largest deviations in Brazil and USA. These two countries used highly clustered reference data, where clustering caused uncertainty given the need to extrapolate to under-sampled areas. We finally provide recommendations to mitigate the effect of under-sampling and to better account for the uncertainties once carbon stocks and flows need to be aggregated in relatively smaller countries. These actions are timely given the global initiatives that aim to upscale UNSEEA carbon accounting.

Biomes of the world under climate change scenarios: increasing aridity and higher temperatures lead to significant shifts in natural vegetation.

The global potential distribution of biomes (natural vegetation) was modelled using 8,959 training points from the BIOME 6000 dataset and a stack of 72 environmental covariates representing terrain and the current climatic conditions based on historical long term averages (1979-2013). An ensemble machine learning model based on stacked regularization was used, with multinomial logistic regression as the meta-learner and spatial blocking (100 km) to deal with spatial autocorrelation of the training points. Results of spatial cross-validation for the BIOME 6000 classes show an overall accuracy of 0.67 and R2logloss of 0.61, with "tropical evergreen broadleaf forest" being the class with highest gain in predictive performances (R2logloss = 0.74) and "prostrate dwarf shrub tundra" the class with the lowest (R2logloss = -0.09) compared to the baseline. Temperature-related covariates were the most important predictors, with the mean diurnal range (BIO2) being shared by all the base-learners (i.e.,random forest, gradient boosted trees and generalized linear models). The model was next used to predict the distribution of future biomes for the periods 2040-2060 and 2061-2080 under three climate change scenarios (RCP 2.6, 4.5 and 8.5). Comparisons of predictions for the three epochs (present, 2040-2060 and 2061-2080) show that increasing aridity and higher temperatures will likely result in significant shifts in natural vegetation in the tropical area (shifts from tropical forests to savannas up to 1.7 ×105 km2 by 2080) and around the Arctic Circle (shifts from tundra to boreal forests up to 2.4 ×105 km2 by 2080). Projected global maps at 1 km spatial resolution are provided as probability and hard classes maps for BIOME 6000 classes and as hard classes maps for the IUCN classes (six aggregated classes). Uncertainty maps (prediction error) are also provided and should be used for careful interpretation of the future projections.

Forest tree species distribution for Europe 2000-2020: mapping potential and realized distributions using spatiotemporal machine learning.

This article describes a data-driven framework based on spatiotemporal machine learning to produce distribution maps for 16 tree species (Abies alba Mill., Castanea sativa Mill., Corylus avellana L., Fagus sylvatica L., Olea europaea L., Picea abies L. H. Karst., Pinus halepensis Mill., Pinus nigra J. F. Arnold, Pinus pinea L., Pinus sylvestris L., Prunus avium L., Quercus cerris L., Quercus ilex L., Quercus robur L., Quercus suber L. and Salix caprea L.) at high spatial resolution (30 m). Tree occurrence data for a total of three million of points was used to train different algorithms: random forest, gradient-boosted trees, generalized linear models, k-nearest neighbors, CART and an artificial neural network. A stack of 305 coarse and high resolution covariates representing spectral reflectance, different biophysical conditions and biotic competition was used as predictors for realized distributions, while potential distribution was modelled with environmental predictors only. Logloss and computing time were used to select the three best algorithms to tune and train an ensemble model based on stacking with a logistic regressor as a meta-learner. An ensemble model was trained for each species: probability and model uncertainty maps of realized distribution were produced for each species using a time window of 4 years for a total of six distribution maps per species, while for potential distributions only one map per species was produced. Results of spatial cross validation show that the ensemble model consistently outperformed or performed as good as the best individual model in both potential and realized distribution tasks, with potential distribution models achieving higher predictive performances (TSS = 0.898, R2 logloss = 0.857) than realized distribution ones on average (TSS = 0.874, R2 logloss = 0.839). Ensemble models for Q. suber achieved the best performances in both potential (TSS = 0.968, R2 logloss = 0.952) and realized (TSS = 0.959, R2 logloss = 0.949) distribution, while P. sylvestris (TSS = 0.731, 0.785, R2 logloss = 0.585, 0.670, respectively, for potential and realized distribution) and P. nigra (TSS = 0.658, 0.686, R2 logloss = 0.623, 0.664) achieved the worst. Importance of predictor variables differed across species and models, with the green band for summer and the Normalized Difference Vegetation Index (NDVI) for fall for realized distribution and the diffuse irradiation and precipitation of the driest quarter (BIO17) being the most frequent and important for potential distribution. On average, fine-resolution models outperformed coarse resolution models (250 m) for realized distribution (TSS = +6.5%, R2 logloss = +7.5%). The framework shows how combining continuous and consistent Earth Observation time series data with state of the art machine learning can be used to derive dynamic distribution maps. The produced predictions can be used to quantify temporal trends of potential forest degradation and species composition change.

Exploring characteristics of national forest inventories for integration with global space-based forest biomass data.

National forest inventories (NFIs) are a reliable source for national forest measurements. However, they are usually not developed for linking with remotely sensed (RS) biomass information. There are increasing needs and opportunities to facilitate this link towards better global and national biomass estimation. Thus, it is important to study and understand NFI characteristics relating to their integration with space-based products; in particular for the tropics where NFIs are quite recent, less frequent, and partially incomplete in several countries. Here, we (1) assessed NFIs in terms of their availability, temporal distribution, and extent in 236 countries from FAO's Global Forest Resources Assessment (FRA) 2020; (2) compared national forest biomass estimates in 2018 from FRA and global space-based Climate Change Initiative (CCI) product in 182 countries considering NFI availability and temporality; and (3) analyzed the latest NFI design characteristics in 46 tropical countries relating to their integration with space-based biomass datasets. We observed significant NFI availability globally and multiple NFIs were mostly found in temperate and boreal countries while most of the single NFI countries (94 %) were in the tropics. The latest NFIs were more recent in the tropics and many countries (35) implemented NFIs from 2016 onwards. The increasing availability and update of NFIs create new opportunities for integration with space-based data at the national level. This is supported by the agreement we found between country biomass estimates for 2018 from FRA and CCI product, with a significantly higher correlation in countries with recent NFIs. We observed that NFI designs varied greatly in tropical countries. For example, the size of the plots ranged from 0.01 to 1 ha and more than three-quarters of the countries had smaller plots of ≤0.25 ha. The existing NFI designs could pose specific challenges for statistical integration with RS data in the tropics. Future NFI and space-based efforts should aim towards a more integrated approach taking advantage of both data streams to improve national estimates and help future data harmonization efforts. Regular NFI efforts can be expanded with the inclusion of some super-site plots to enhance data integration with currently available space-based applications. Issues related to cost implications versus improvements in the accuracy, timeliness, and sustainability of national forest biomass estimation should be further explored.

Retrieval of Hyperspectral Information from Multispectral Data for Perennial Ryegrass Biomass Estimation.

The use of spectral data is seen as a fast and non-destructive method capable of monitoring pasture biomass. Although there is great potential in this technique, both end users and sensor manufacturers are uncertain about the necessary sensor specifications and achievable accuracies in an operational scenario. This study presents a straightforward parametric method able to accurately retrieve the hyperspectral signature of perennial ryegrass (Lolium perenne) canopies from multispectral data collected within a two-year period in Australia and the Netherlands. The retrieved hyperspectral data were employed to generate optimal indices and continuum-removed spectral features available in the scientific literature. For performance comparison, both these simulated features and a set of currently employed vegetation indices, derived from the original band values, were used as inputs in a random forest algorithm and accuracies of both methods were compared. Our results have shown that both sets of features present similar accuracies (root mean square error (RMSE) ≈490 and 620 kg DM/ha) when assessed in cross-validation and spatial cross-validation, respectively. These results suggest that for pasture biomass retrieval solely from top-of-canopy reflectance (ranging from 550 to 790 nm), better performing methods do not rely on the use of hyperspectral or, yet, in a larger number of bands than those already available in current sensors.

Data synthesis for crop variety evaluation. A review.

Crop varieties should fulfill multiple requirements, including agronomic performance and product quality. Variety evaluations depend on data generated from field trials and sensory analyses, performed with different levels of participation from farmers and consumers. Such multi-faceted variety evaluation is expensive and time-consuming; hence, any use of these data should be optimized. Data synthesis can help to take advantage of existing and new data, combining data from different sources and combining it with expert knowledge to produce new information and understanding that supports decision-making. Data synthesis for crop variety evaluation can partly build on extant experiences and methods, but it also requires methodological innovation. We review the elements required to achieve data synthesis for crop variety evaluation, including (1) data types required for crop variety evaluation, (2) main challenges in data management and integration, (3) main global initiatives aiming to solve those challenges, (4) current statistical approaches to combine data for crop variety evaluation and (5) existing data synthesis methods used in evaluation of varieties to combine different datasets from multiple data sources. We conclude that currently available methods have the potential to overcome existing barriers to data synthesis and could set in motion a virtuous cycle that will encourage researchers to share data and collaborate on data-driven research.

Optimization of rain gauge sampling density for river discharge prediction using Bayesian calibration.

River discharges are often predicted based on a calibrated rainfall-runoff model. The major sources of uncertainty, namely input, parameter and model structural uncertainty must all be taken into account to obtain realistic estimates of the accuracy of discharge predictions. Over the past years, Bayesian calibration has emerged as a suitable method for quantifying uncertainty in model parameters and model structure, where the latter is usually modelled by an additive or multiplicative stochastic term. Recently, much work has also been done to include input uncertainty in the Bayesian framework. However, the use of geostatistical methods for characterizing the prior distribution of the catchment rainfall is underexplored, particularly in combination with assessments of the influence of increasing or decreasing rain gauge network density on discharge prediction accuracy. In this article we integrate geostatistics and Bayesian calibration to analyze the effect of rain gauge density on river discharge prediction accuracy. We calibrated the HBV hydrological model while accounting for input, initial state, model parameter and model structural uncertainty, and also taking uncertainties in the discharge measurements into account. Results for the Thur basin in Switzerland showed that model parameter uncertainty was the main contributor to the joint posterior uncertainty. We also showed that a low rain gauge density is enough for the Bayesian calibration, and that increasing the number of rain gauges improved model prediction until reaching a density of one gauge per 340 km2. While the optimal rain gauge density is case-study specific, we make recommendations on how to handle input uncertainty in Bayesian calibration for river discharge prediction and present the methodology that may be used to carry out such experiments.

Using household survey data to identify large-scale food security patterns across Uganda.

To target food security interventions for smallholder households, decision makers need large-scale information, such as maps on poverty, food security and key livelihood activities. Such information is often based on expert knowledge or aggregated data, despite the fact that food security and poverty are driven largely by processes at the household level. At present, it is unclear if and how household level information can contribute to the spatial prediction of such welfare indicators or to what extent local variability is ignored by current mapping efforts. A combination of geo-referenced household level information with spatially continuous information is an underused approach to quantify local and large-scale variation, while it can provide a direct estimate of the variability of welfare indicators at the most relevant scale. We applied a stepwise regression kriging procedure to translate point information to spatially explicit patterns and create country-wide predictions with associated uncertainty estimates for indicators on food availability and related livelihood activities using household survey data from Uganda. With few exceptions, predictions of the indicators were weak, highlighting the difficulty in capturing variability at larger scale. Household explanatory variables identified little additional variation compared to environmental explanatory variables alone. Spatial predictability was strongest for indicators whose distribution was determined by environmental gradients. In contrast, indicators of crops that were more ubiquitously present across agroecological zones showed large local variation, which often overruled large-scale patterns. Our procedure adds to existing approaches that often only show large-scale patterns by revealing that local variation in welfare is large. Interventions that aim to target the poor must recognise that diversity in livelihood activities for income generation within any given area often overrides the variability of livelihood activities between distant regions in the country.

Design and Implementation of an Interactive Web-Based Near Real-Time Forest Monitoring System.

This paper describes an interactive web-based near real-time (NRT) forest monitoring system using four levels of geographic information services: 1) the acquisition of continuous data streams from satellite and community-based monitoring using mobile devices, 2) NRT forest disturbance detection based on satellite time-series, 3) presentation of forest disturbance data through a web-based application and social media and 4) interaction of the satellite based disturbance alerts with the end-user communities to enhance the collection of ground data. The system is developed using open source technologies and has been implemented together with local experts in the UNESCO Kafa Biosphere Reserve, Ethiopia. The results show that the system is able to provide easy access to information on forest change and considerably improves the collection and storage of ground observation by local experts. Social media leads to higher levels of user interaction and noticeably improves communication among stakeholders. Finally, an evaluation of the system confirms the usability of the system in Ethiopia. The implemented system can provide a foundation for an operational forest monitoring system at the national level for REDD+ MRV applications.

Spatial relationship between climatologies and changes in global vegetation activity.

Vegetation forms a main component of the terrestrial biosphere and plays a crucial role in land-cover and climate-related studies. Activity of vegetation systems is commonly quantified using remotely sensed vegetation indices (VI). Extensive reports on temporal trends over the past decades in time series of such indices can be found in literature. However, little remains known about the processes underlying these changes at large spatial scales. In this study, we aimed at quantifying the spatial relationship between changes in potential climatic growth constraints (i.e. temperature, precipitation and incident solar radiation) and changes in vegetation activity (1982-2008). We demonstrate an additive spatial model with 0.5° resolution, consisting of a regression component representing climate-associated effects and a spatially correlated field representing the combined influence of other factors, including land-use change. Little over 50% of the spatial variance could be attributed to changes in climatologies; conspicuously, many greening trends and browning hotspots in Argentina and Australia. The nonassociated model component may contain large-scale human interventions, feedback mechanisms or natural effects, which were not captured by the climatologies. Browning hotspots in this component were especially found in subequatorial Africa. On the scale of land-cover types, strongest relationships between climatologies and vegetation activity were found in forests, including indications for browning under warming conditions (analogous to the divergence issue discussed in dendroclimatology).

Where and when should sensors move? Sampling using the expected value of information.

In case of an environmental accident, initially available data are often insufficient for properly managing the situation. In this paper, new sensor observations are iteratively added to an initial sample by maximising the global expected value of information of the points for decision making. This is equivalent to minimizing the aggregated expected misclassification costs over the study area. The method considers measurement error and different costs for class omissions and false class commissions. Constraints imposed by a mobile sensor web are accounted for using cost distances to decide which sensor should move to the next sample location. The method is demonstrated using synthetic examples of static and dynamic phenomena. This allowed computation of the true misclassification costs and comparison with other sampling approaches. The probability of local contamination levels being above a given critical threshold were computed by indicator kriging. In the case of multiple sensors being relocated simultaneously, a genetic algorithm was used to find sets of suitable new measurement locations. Otherwise, all grid nodes were searched exhaustively, which is computationally demanding. In terms of true misclassification costs, the method outperformed random sampling and sampling based on minimisation of the kriging variance.

Mobile devices for community-based REDD+ monitoring: a case study for Central Vietnam.

Monitoring tropical deforestation and forest degradation is one of the central elements for the Reduced Emissions from Deforestation and Forest Degradation in developing countries (REDD+) scheme. Current arrangements for monitoring are based on remote sensing and field measurements. Since monitoring is the periodic process of assessing forest stands properties with respect to reference data, adopting the current REDD+ requirements for implementing monitoring at national levels is a challenging task. Recently, the advancement in Information and Communications Technologies (ICT) and mobile devices has enabled local communities to monitor their forest in a basic resource setting such as no or slow internet connection link, limited power supply, etc. Despite the potential, the use of mobile device system for community based monitoring (CBM) is still exceptional and faces implementation challenges. This paper presents an integrated data collection system based on mobile devices that streamlines the community-based forest monitoring data collection, transmission and visualization process. This paper also assesses the accuracy and reliability of CBM data and proposes a way to fit them into national REDD+ Monitoring, Reporting and Verification (MRV) scheme. The system performance is evaluated at Tra Bui commune, Quang Nam province, Central Vietnam, where forest carbon and change activities were tracked. The results show that the local community is able to provide data with accuracy comparable to expert measurements (index of agreement greater than 0.88), but against lower costs. Furthermore, the results confirm that communities are more effective to monitor small scale forest degradation due to subsistence fuel wood collection and selective logging, than high resolution remote sensing SPOT imagery.

Development of a Dynamic Web Mapping Service for Vegetation Productivity Using Earth Observation and in situ Sensors in a Sensor Web Based Approach.

This paper describes the development of a sensor web based approach which combines earth observation and in situ sensor data to derive typical information offered by a dynamic web mapping service (WMS). A prototype has been developed which provides daily maps of vegetation productivity for the Netherlands with a spatial resolution of 250 m. Daily available MODIS surface reflectance products and meteorological parameters obtained through a Sensor Observation Service (SOS) were used as input for a vegetation productivity model. This paper presents the vegetation productivity model, the sensor data sources and the implementation of the automated processing facility. Finally, an evaluation is made of the opportunities and limitations of sensor web based approaches for the development of web services which combine both satellite and in situ sensor sources.