Spatial heterogeneity of global forest aboveground carbon stocks and fluxes constrained by spaceborne lidar data and mechanistic modeling.

Forest carbon is a large and uncertain component of the global carbon cycle. An important source of complexity is the spatial heterogeneity of vegetation vertical structure and extent, which results from variations in climate, soils, and disturbances and influences both contemporary carbon stocks and fluxes. Recent advances in remote sensing and ecosystem modeling have the potential to significantly improve the characterization of vegetation structure and its resulting influence on carbon. Here, we used novel remote sensing observations of tree canopy height collected by two NASA spaceborne lidar missions, Global Ecosystem Dynamics Investigation and ICE, Cloud, and Land Elevation Satellite 2, together with a newly developed global Ecosystem Demography model (v3.0) to characterize the spatial heterogeneity of global forest structure and quantify the corresponding implications for forest carbon stocks and fluxes. Multiple-scale evaluations suggested favorable results relative to other estimates including field inventory, remote sensing-based products, and national statistics. However, this approach utilized several orders of magnitude more data (3.77 billion lidar samples) on vegetation structure than used previously and enabled a qualitative increase in the spatial resolution of model estimates achievable (0.25° to 0.01°). At this resolution, process-based models are now able to capture detailed spatial patterns of forest structure previously unattainable, including patterns of natural and anthropogenic disturbance and recovery. Through the novel integration of new remote sensing data and ecosystem modeling, this study bridges the gap between existing empirically based remote sensing approaches and process-based modeling approaches. This study more generally demonstrates the promising value of spaceborne lidar observations for advancing carbon modeling at a global scale.

Toward a forest biomass reference measurement system for remote sensing applications.

Forests contribute to climate change mitigation through carbon storage and uptake, but the extent to which this carbon pool varies in space and time is still poorly known. Several Earth Observation missions have been specifically designed to address this issue, for example, NASA's GEDI, NASA-ISRO's NISAR and ESA's BIOMASS. Yet, all these missions' products require independent and consistent validation. A permanent, global, in situ, site-based forest biomass reference measurement system relying on ground data of the highest possible quality is therefore needed. Here, we have assembled a list of almost 200 high-quality sites through an in-depth review of the literature and expert knowledge. In this study, we explore how representative these sites are in terms of their coverage of environmental conditions, geographical space and biomass-related forest structure, compared to those experienced by forests worldwide. This work also aims at identifying which sites are the most representative, and where to invest to improve the representativeness of the proposed system. We show that the environmental coverage of the system does not seem to improve after at least the 175 most representative sites are included, but geographical and structural coverages continue to improve as more sites are added. We highlight the areas of poor environmental, geographical, or structural coverage, including, but not limited to, Canada, the western half of the USA, Mexico, Patagonia, Angola, Zambia, eastern Russia, and tropical and subtropical highlands (e.g. in Colombia, the Himalayas, Borneo, Papua). For the proposed system to succeed, we stress that (1) data must be collected and processed applying the same standards across all countries and continents; (2) system establishment and management must be inclusive and equitable, with careful consideration of working conditions; and (3) training and site partner involvement in downstream activities should be mandatory.

The use of sun elevation angle for stereogrammetric boreal forest height in open canopies.

Stereogrammetry applied to globally available high resolution spaceborne imagery (HRSI; < 5 m spatial resolution) yields fine-scaled digital surface models (DSMs) of elevation. These DSMs may represent elevations that range from the ground to the vegetation canopy surface, are produced from stereoscopic image pairs (stereopairs) that have a variety of acquisition characteristics, and have been coupled with lidar data of forest structure and ground surface elevation to examine forest height. This work explores surface elevations from HRSI DSMs derived from two types of acquisitions in open canopy forests. We (1) apply an automated mass-production stereogrammetry workflow to along-track HRSI stereopairs, (2) identify multiple spatially coincident DSMs whose stereopairs were acquired under different solar geometry, (3) vertically co-register these DSMs using coincident spaceborne lidar footprints (from ICESat-GLAS) as reference, and (4) examine differences in surface elevations between the reference lidar and the co-registered HRSI DSMs associated with two general types of acquisitions (DSM types) from different sun elevation angles. We find that these DSM types, distinguished by sun elevation angle at the time of stereopair acquisition, are associated with different surface elevations estimated from automated stereogrammetry in open canopy forests. For DSM values with corresponding reference ground surface elevation from spaceborne lidar footprints in open canopy northern Siberian Larix forests with slopes < 10°, our results show that HRSI DSMs acquired with sun elevation angles > 35° and < 25° (during snow-free conditions) produced characteristic and consistently distinct distributions of elevation differences from reference lidar. The former include DSMs of near-ground surfaces with root mean square errors < 0.68 m relative to lidar. The latter, particularly those with angles < 10°, show distributions with larger differences from lidar that are associated with open canopy forests whose vegetation surface elevations are captured. Terrain aspect did not have a strong effect on the distribution of vegetation surfaces. Using the two DSM types together, the distribution of DSM-differenced heights in forests (µ = 6.0 m, σ = 1.4 m) was consistent with the distribution of plot-level mean tree heights (µ = 6.5 m, σ = 1.2 m). We conclude that the variation in sun elevation angle at time of stereopair acquisition can create illumination conditions conducive for capturing elevations of surfaces either near the ground or associated with vegetation canopy. Knowledge of HRSI acquisition solar geometry and snow cover can be used to understand and combine stereogrammetric surface elevation estimates to co-register and difference overlapping DSMs, providing a means to map forest height at fine scales, resolving the vertical structure of groups of trees from spaceborne platforms in open canopy forests.

Intergovernmental Panel on Climate Change (IPCC) Tier 1 forest biomass estimates from Earth Observation.

Aboveground biomass density (AGBD) estimates from Earth Observation (EO) can be presented with the consistency standards mandated by United Nations Framework Convention on Climate Change (UNFCCC). This article delivers AGBD estimates, in the format of Intergovernmental Panel on Climate Change (IPCC) Tier 1 values for natural forests, sourced from National Aeronautics and Space Administration's (NASA's) Global Ecosystem Dynamics Investigation (GEDI) and Ice, Cloud and land Elevation Satellite (ICESat-2), and European Space Agency's (ESA's) Climate Change Initiative (CCI). It also provides the underlying classification used by the IPCC as geospatial layers, delineating global forests by ecozones, continents and status (primary, young (≤20 years) and old secondary (>20 years)). The approaches leverage complementary strengths of various EO-derived datasets that are compiled in an open-science framework through the Multi-mission Algorithm and Analysis Platform (MAAP). This transparency and flexibility enables the adoption of any new incoming datasets in the framework in the future. The EO-based AGBD estimates are expected to be an independent contribution to the IPCC Emission Factors Database in support of UNFCCC processes, and the forest classification expected to support the generation of other policy-relevant datasets while reflecting ongoing shifts in global forests with climate change.

Toward a roadmap for space-based observations of the land sector for the UNFCCC global stocktake.

Space-based remote sensing can make an important contribution toward monitoring greenhouse gas emissions and removals from the agriculture, forestry, and other land use (AFOLU) sector, and to understanding and addressing human-caused climate change through the UNFCCC Paris Agreement. Space agencies have begun to coordinate their efforts to identify needs, collect and harmonize available data and efforts, and plan and maintain a long-term roadmap for observations. International cooperation is crucial in developing and realizing the roadmap, and the Committee on Earth Observation Satellites (CEOS) is a key coordinating driver of this effort. Here, we first identify the data and information that will be useful to support the global stocktake (GST) of the Paris Agreement. Then, the paper explains how existing and planned space-based capabilities and products can be used and combined, particularly in the land use sector, and provides a workflow for their harmonization and contribution to greenhouse gas inventories and assessments at the national and global level.

New 3D measurements of large redwood trees for biomass and structure.

Large trees are disproportionately important in terms of their above ground biomass (AGB) and carbon storage, as well as their wider impact on ecosystem structure. They are also very hard to measure and so tend to be underrepresented in measurements and models of AGB. We show the first detailed 3D terrestrial laser scanning (TLS) estimates of the volume and AGB of large coastal redwood Sequoia sempervirens trees from three sites in Northern California, representing some of the highest biomass ecosystems on Earth. Our TLS estimates agree to within 2% AGB with a species-specific model based on detailed manual crown mapping of 3D tree structure. However TLS-derived AGB was more than 30% higher compared to widely-used general (non species-specific) allometries. We derive an allometry from TLS that spans a much greater range of tree size than previous models and so is potentially better-suited for use with new Earth Observation data for these exceptionally high biomass areas. We suggest that where possible, TLS and crown mapping should be used to provide complementary, independent 3D structure measurements of these very large trees.

The GEDI Simulator: A Large-Footprint Waveform Lidar Simulator for Calibration and Validation of Spaceborne Missions.

NASA's Global Ecosystem Dynamics Investigation (GEDI) is a spaceborne lidar mission which will produce near global (51.6°S to 51.6°N) maps of forest structure and above-ground biomass density during its 2-year mission. GEDI uses a waveform simulator for calibration of algorithms and assessing mission accuracy. This paper implements a waveform simulator, using the method proposed in Blair and Hofton (1999; https://doi.org/10.1029/1999GL010484), and builds upon that work by adding instrument noise and by validating simulated waveforms across a range of forest types, airborne laser scanning (ALS) instruments, and survey configurations. The simulator was validated by comparing waveform metrics derived from simulated waveforms against those derived from observed large-footprint, full-waveform lidar data from NASA's airborne Land, Vegetation, and Ice Sensor (LVIS). The simulator was found to produce waveform metrics with a mean bias of less than 0.22 m and a root-mean-square error of less than 5.7 m, as long as the ALS data had sufficient pulse density. The minimum pulse density required depended upon the instrument. Measurement errors due to instrument noise predicted by the simulator were within 1.5 m of those from observed waveforms and 70-85% of variance in measurement error was explained. Changing the ALS survey configuration had no significant impact on simulated metrics, suggesting that the ALS pulse density is a sufficient metric of simulator accuracy across the range of conditions and instruments tested. These results give confidence in the use of the simulator for the pre-launch calibration and performance assessment of the GEDI mission.

New Opportunities for Forest Remote Sensing Through Ultra-High-Density Drone Lidar.

Current and planned space missions will produce aboveground biomass density data products at varying spatial resolution. Calibration and validation of these data products is critically dependent on the existence of field estimates of aboveground biomass and coincident remote sensing data from airborne or terrestrial lidar. There are few places that meet these requirements, and they are mostly in the northern hemisphere and temperate zone. Here we summarize the potential for low-altitude drones to produce new observations in support of mission science. We describe technical requirements for producing high-quality measurements from autonomous platforms and highlight differences among commercially available laser scanners and drone aircraft. We then describe a case study using a heavy-lift autonomous helicopter in a temperate mountain forest in the southern Czech Republic in support of calibration and validation activities for the NASA Global Ecosystem Dynamics Investigation. Low-altitude flight using drones enables the collection of ultra-high-density point clouds using wider laser scan angles than have been possible from traditional airborne platforms. These measurements can be precise and accurate and can achieve measurement densities of thousands of points · m-2. Analysis of surface elevation measurements on a heterogeneous target observed 51 days apart indicates that the realized range accuracy is 2.4 cm. The single-date precision is 2.1-4.5 cm. These estimates are net of all processing artifacts and geolocation errors under fully autonomous flight. The 3D model produced by these data can clearly resolve branch and stem structure that is comparable to terrestrial laser scans and can be acquired rapidly over large landscapes at a fraction of the cost of traditional airborne laser scanning.

Monitoring individual tree-based change with airborne lidar.

Understanding the carbon flux of forests is critical for constraining the global carbon cycle and managing forests to mitigate climate change. Monitoring forest growth and mortality rates is critical to this effort, but has been limited in the past, with estimates relying primarily on field surveys. Advances in remote sensing enable the potential to monitor tree growth and mortality across landscapes. This work presents an approach to measure tree growth and loss using multidate lidar campaigns in a high-biomass forest in California, USA. Individual tree crowns were delineated in 2008 and again in 2013 using a 3D crown segmentation algorithm, with derived heights and crown radii extracted and used to estimate individual tree aboveground biomass. Tree growth, loss, and aboveground biomass were analyzed with respect to tree height and crown radius. Both tree growth and loss rates decrease with increasing tree height, following the expectation that trees slow in growth rate as they age. Additionally, our aboveground biomass analysis suggests that, while the system is a net source of aboveground carbon, these carbon dynamics are governed by size class with the largest sources coming from the loss of a relatively small number of large individuals. This study demonstrates that monitoring individual tree-based growth and loss can be conducted with multidate airborne lidar, but these methods remain relatively immature. Disparities between lidar acquisitions were particularly difficult to overcome and decreased the sample of trees analyzed for growth rate in this study to 21% of the full number of delineated crowns. However, this study illuminates the potential of airborne remote sensing for ecologically meaningful forest monitoring at an individual tree level. As methods continue to improve, airborne multidate lidar will enable a richer understanding of the drivers of tree growth, loss, and aboveground carbon flux.

Implications of allometric model selection for county-level biomass mapping.

BACKGROUND: Carbon accounting in forests remains a large area of uncertainty in the global carbon cycle. Forest aboveground biomass is therefore an attribute of great interest for the forest management community, but the accuracy of aboveground biomass maps depends on the accuracy of the underlying field estimates used to calibrate models. These field estimates depend on the application of allometric models, which often have unknown and unreported uncertainties outside of the size class or environment in which they were developed. RESULTS: Here, we test three popular allometric approaches to field biomass estimation, and explore the implications of allometric model selection for county-level biomass mapping in Sonoma County, California. We test three allometric models: Jenkins et al. (For Sci 49(1): 12-35, 2003), Chojnacky et al. (Forestry 87(1): 129-151, 2014) and the US Forest Service's Component Ratio Method (CRM). We found that Jenkins and Chojnacky models perform comparably, but that at both a field plot level and a total county level there was a ~ 20% difference between these estimates and the CRM estimates. Further, we show that discrepancies are greater in high biomass areas with high canopy covers and relatively moderate heights (25-45 m). The CRM models, although on average ~ 20% lower than Jenkins and Chojnacky, produce higher estimates in the tallest forests samples (> 60 m), while Jenkins generally produces higher estimates of biomass in forests < 50 m tall. Discrepancies do not continually increase with increasing forest height, suggesting that inclusion of height in allometric models is not primarily driving discrepancies. Models developed using all three allometric models underestimate high biomass and overestimate low biomass, as expected with random forest biomass modeling. However, these deviations were generally larger using the Jenkins and Chojnacky allometries, suggesting that the CRM approach may be more appropriate for biomass mapping with lidar. CONCLUSIONS: These results confirm that allometric model selection considerably impacts biomass maps and estimates, and that allometric model errors remain poorly understood. Our findings that allometric model discrepancies are not explained by lidar heights suggests that allometric model form does not drive these discrepancies. A better understanding of the sources of allometric model errors, particularly in high biomass systems, is essential for improved forest biomass mapping.

Author Correction: ISS observations offer insights into plant function.

In the version of this Comment previously published, in Box 1, the spacing of the GEDI footprints should have read 60 m along the track, not 25 m. Also the second affiliation for Susan Ustin was incorrect, she is only associated with the University of California, Davis. These errors have now been corrected.

Local discrepancies in continental scale biomass maps: a case study over forested and non-forested landscapes in Maryland, USA.

BACKGROUND: Continental-scale aboveground biomass maps are increasingly available, but their estimates vary widely, particularly at high resolution. A comprehensive understanding of map discrepancies is required to improve their effectiveness in carbon accounting and local decision-making. To this end, we compare four continental-scale maps with a recent high-resolution lidar-derived biomass map over Maryland, USA. We conduct detailed comparisons at pixel-, county-, and state-level. RESULTS: Spatial patterns of biomass are broadly consistent in all maps, but there are large differences at fine scales (RMSD 48.5-92.7 Mg ha-1). Discrepancies reduce with aggregation and the agreement among products improves at the county level. However, continental scale maps exhibit residual negative biases in mean (33.0-54.6 Mg ha-1) and total biomass (3.5-5.8 Tg) when compared to the high-resolution lidar biomass map. Three of the four continental scale maps reach near-perfect agreement at ~4 km and onward but do not converge with the high-resolution biomass map even at county scale. At the State level, these maps underestimate biomass by 30-80 Tg in forested and 40-50 Tg in non-forested areas. CONCLUSIONS: Local discrepancies in continental scale biomass maps are caused by factors including data inputs, modeling approaches, forest/non-forest definitions and time lags. There is a net underestimation over high biomass forests and non-forested areas that could impact carbon accounting at all levels. Local, high-resolution lidar-derived biomass maps provide a valuable bottom-up reference to improve the analysis and interpretation of large-scale maps produced in carbon monitoring systems.