A comprehensive framework for assessing the accuracy and uncertainty of global above-ground biomass maps

Over the past decade, several global maps of above-ground biomass (AGB) have been produced, but they exhibit significant differences that reduce their value for climate and carbon cycle modelling, and also for national estimates of forest carbon stocks and their changes. The number of such maps is anticipated to increase because of new satellite missions dedicated to measuring AGB. Objective and consistent methods to estimate the accuracy and uncertainty of AGB maps are therefore urgently needed. This paper develops and demonstrates a framework aimed at achieving this. The framework provides a means to compare AGB maps with AGB estimates from a global collection of National Forest Inventories and research plots that accounts for the uncertainty of plot AGB errors. This uncertainty depends strongly on plot size, and is dominated by the combined errors from tree measurements and allometric models (inter-quartile range of their standard deviation (SD) = 30–151 Mg ha−1). Estimates of sampling errors are also important, especially in the most common case where plots are smaller than map pixels (SD = 16–44 Mg ha−1). Plot uncertainty estimates are used to calculate the minimum-variance linear unbiased estimates of the mean forest AGB when averaged to 0.1∘. These are used to assess four AGB maps: Baccini (2000), GEOCARBON (2008), GlobBiomass (2010) and CCI Biomass (2017). Map bias, estimated using the differences between the plot and 0.1∘ map averages, is modelled using random forest regression driven by variables shown to affect the map estimates. The bias model is particularly sensitive to the map estimate of AGB and tree cover, and exhibits strong regional biases. Variograms indicate that AGB map errors have map-specific spatial correlation up to a range of 50–104 km, which increases the variance of spatially aggregated AGB map estimates compared to when pixel errors are independent. After bias adjustment, total pantropical AGB and its associated SD are derived for the four map epochs. This total becomes closer to the value estimated by the Forest Resources Assessment after every epoch and shows a similar decrease. The framework is applicable to both local and global-scale analysis, and is available at https://github.com/arnanaraza/PlotToMap. Our study therefore constitutes a major step towards improved AGB map validation and improvement

A comprehensive framework for assessing the accuracy and uncertainty of global above-ground biomass maps

International audienceOver the past decade, several global maps of above-ground biomass (AGB) have been produced, but they exhibit significant differences that reduce their value for climate and carbon cycle modelling, and also for national estimates of forest carbon stocks and their changes. The number of such maps is anticipated to increase because of new satellite missions dedicated to measuring AGB. Objective and consistent methods to estimate the accuracy and uncertainty of AGB maps are therefore urgently needed. This paper develops and demonstrates a framework aimed at achieving this. The framework provides a means to compare AGB maps with AGB estimates from a global collection of National Forest Inventories and research plots that accounts for the uncertainty of plot AGB errors. This uncertainty depends strongly on plot size, and is dominated by the combined errors from tree measurements and allometric models (inter-quartile range of their standard deviation (SD) = 30–151 Mg ha−1). Estimates of sampling errors are also important, especially in the most common case where plots are smaller than map pixels (SD = 16–44 Mg ha−1). Plot uncertainty estimates are used to calculate the minimum-variance linear unbiased estimates of the mean forest AGB when averaged to 0.1∘. These are used to assess four AGB maps: Baccini (2000), GEOCARBON (2008), GlobBiomass (2010) and CCI Biomass (2017). Map bias, estimated using the differences between the plot and 0.1∘ map averages, is modelled using random forest regression driven by variables shown to affect the map estimates. The bias model is particularly sensitive to the map estimate of AGB and tree cover, and exhibits strong regional biases. Variograms indicate that AGB map errors have map-specific spatial correlation up to a range of 50–104 km, which increases the variance of spatially aggregated AGB map estimates compared to when pixel errors are independent. After bias adjustment, total pantropical AGB and its associated SD are derived for the four map epochs. This total becomes closer to the value estimated by the Forest Resources Assessment after every epoch and shows a similar decrease. The framework is applicable to both local and global-scale analysis, and is available at https://github.com/arnanaraza/PlotToMap. Our study therefore constitutes a major step towards improved AGB map validation and improvement

Aboveground forest biomass varies across continents, ecological zones and successional stages: refined IPCC default values for tropical and subtropical forests

For monitoring and reporting forest carbon stocks and fluxes, many countries in the tropics and subtropics rely on default values of forest aboveground biomass (AGB) from the Intergovernmental Panel on Climate Change (IPCC) guidelines for National Greenhouse Gas (GHG) Inventories. Default IPCC forest AGB values originated from 2006, and are relatively crude estimates of average values per continent and ecological zone. The 2006 default values were based on limited plot data available at the time, methods for their derivation were not fully clear, and no distinction between successional stages was made. As part of the 2019 Refinement to the 2006 IPCC Guidelines for GHG Inventories, we updated the default AGB values for tropical and subtropical forests based on AGB data from >25 000 plots in natural forests and a global AGB map where no plot data were available. We calculated refined AGB default values per continent, ecological zone, and successional stage, and provided a measure of uncertainty. AGB in tropical and subtropical forests varies by an order of magnitude across continents, ecological zones, and successional stage. Our refined default values generally reflect the climatic gradients in the tropics, with more AGB in wetter areas. AGB is generally higher in old-growth than in secondary forests, and higher in older secondary (regrowth >20 years old and degraded/logged forests) than in young secondary forests (20 years old). While refined default values for tropical old-growth forest are largely similar to the previous 2006 default values, the new default values are 4.0-7.7-fold lower for young secondary forests. Thus, the refined values will strongly alter estimated carbon stocks and fluxes, and emphasize the critical importance of old-growth forest conservation. We provide a reproducible approach to facilitate future refinements and encourage targeted efforts to establish permanent plots in areas with data gaps

Combined fit to the spectrum and composition data measured by the Pierre Auger Observatory including magnetic horizon effects

The measurements by the Pierre Auger Observatory of the energy spectrum and mass composition of cosmic rays can be interpreted assuming the presence of two extragalactic source populations, one dominating the flux at energies above a few EeV and the other below. To fit the data ignoring magnetic field effects, the high-energy population needs to accelerate a mixture of nuclei with very hard spectra, at odds with the approximate E$^{-2}$ shape expected from diffusive shock acceleration. The presence of turbulent extragalactic magnetic fields in the region between the closest sources and the Earth can significantly modify the observed CR spectrum with respect to that emitted by the sources, reducing the flux of low-rigidity particles that reach the Earth. We here take into account this magnetic horizon effect in the combined fit of the spectrum and shower depth distributions, exploring the possibility that a spectrum for the high-energy population sources with a shape closer to E$^{-2}$ be able to explain the observations

Studies of the mass composition of cosmic rays and proton-proton interaction cross-sections at ultra-high energies with the Pierre Auger Observatory

In this work, we present an estimate of the cosmic-ray mass composition from the distributions of the depth of the shower maximum (Xmax) measured by the fluorescence detector of the Pierre Auger Observatory. We discuss the sensitivity of the mass composition measurements to the uncertainties in the properties of the hadronic interactions, particularly in the predictions of the particle interaction cross-sections. For this purpose, we adjust the fractions of cosmic-ray mass groups to fit the data with Xmax distributions from air shower simulations. We modify the proton-proton cross-sections at ultra-high energies, and the corresponding air shower simulations with rescaled nucleus-air cross-sections are obtained via Glauber theory. We compare the energy-dependent composition of ultra-high-energy cosmic rays obtained for the different extrapolations of the proton-proton cross-sections from low-energy accelerator data

Study of downward Terrestrial Gamma-ray Flashes with the surface detector of the Pierre Auger Observatory

The surface detector (SD) of the Pierre Auger Observatory, consisting of 1660 water-Cherenkov detectors (WCDs), covers 3000 km2 in the Argentinian pampa. Thanks to the high efficiency of WCDs in detecting gamma rays, it represents a unique instrument for studying downward Terrestrial Gamma-ray Flashes (TGFs) over a large area. Peculiar events, likely related to downward TGFs, were detected at the Auger Observatory. Their experimental signature and time evolution are very different from those of a shower produced by an ultrahigh-energy cosmic ray. They happen in coincidence with low thunderclouds and lightning, and their large deposited energy at the ground is compatible with that of a standard downward TGF with the source a few kilometers above the ground. A new trigger algorithm to increase the TGF-like event statistics was installed in the whole array. The study of the performance of the new trigger system during the lightning season is ongoing and will provide a handle to develop improved algorithms to implement in the Auger upgraded electronic boards. The available data sample, even if small, can give important clues about the TGF production models, in particular, the shape of WCD signals. Moreover, the SD allows us to observe more than one point in the TGF beam, providing information on the emission angle

The dynamic range of the upgraded surface-detector stations of AugerPrime

The detection of ultra-high-energy cosmic rays by means of giant detector arrays is often limited by the saturation of the recorded signals near the impact point of the shower core at the ground, where the particle density dramatically increases. The saturation affects in particular the highest energy events, worsening the systematic uncertainties in the reconstruction of the shower characteristics. The upgrade of the Pierre Auger Observatory, called AugerPrime, includes the installation of an 1-inch Small PhotoMultiplier Tube (SPMT) inside each water-Cherenkov station (WCD) of the surface detector array. The SPMT allows an unambiguous measurement of signals down to about 250m from the shower core, thus reducing the number of events featuring a saturated station to a negligible level. In addition, a 3.8m2 plastic scintillator (Scintillator Surface Detector, SSD) is installed on top of each WCD. The SSD is designed to match the WCD (with SPMT) dynamic range, providing a complementary measurement of the shower components up to the highest energies. In this work, the design and performances of the upgraded AugerPrime surface-detector stations in the extended dynamic range are described, highlighting the accuracy of the measurements. A first analysis employing the unsaturated signals in the event reconstruction is also presented