Mapping migratory bird prevalence using remote sensing data fusion

Background

Improved maps of species distributions are important for effective management of wildlife under increasing anthropogenic pressures. Recent advances in lidar and radar remote sensing have shown considerable potential for mapping forest structure and habitat characteristics across landscapes. However, their relative efficacies and integrated use in habitat mapping remain largely unexplored. We evaluated the use of lidar, radar and multispectral remote sensing data in predicting multi-year bird detections or prevalence for 8 migratory songbird species in the unfragmented temperate deciduous forests of New Hampshire, USA.

Methodology and Principal Findings

A set of 104 predictor variables describing vegetation vertical structure and variability from lidar, phenology from multispectral data and backscatter properties from radar data were derived. We tested the accuracies of these variables in predicting prevalence using Random Forests regression models. All data sets showed more than 30% predictive power with radar models having the lowest and multi-sensor synergy (“fusion”) models having highest accuracies. Fusion explained between 54% and 75% variance in prevalence for all the birds considered. Stem density from discrete return lidar and phenology from multispectral data were among the best predictors. Further analysis revealed different relationships between the remote sensing metrics and bird prevalence. Spatial maps of prevalence were consistent with known habitat preferences for the bird species.

Conclusion and Significance

Our results highlight the potential of integrating multiple remote sensing data sets using machine-learning methods to improve habitat mapping. Multi-dimensional habitat structure maps such as those generated from this study can significantly advance forest management and ecological research by facilitating fine-scale studies at both stand and landscape level.

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.

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.

The influence of vegetation height heterogeneity on forest and woodland bird species richness across the United States

Avian diversity is under increasing pressures. It is thus critical to understand the ecological variables that contribute to large scale spatial distribution of avian species diversity. Traditionally, studies have relied primarily on two-dimensional habitat structure to model broad scale species richness. Vegetation vertical structure is increasingly used at local scales. However, the spatial arrangement of vegetation height has never been taken into consideration. Our goal was to examine the efficacies of three-dimensional forest structure, particularly the spatial heterogeneity of vegetation height in improving avian richness models across forested ecoregions in the U.S. We developed novel habitat metrics to characterize the spatial arrangement of vegetation height using the National Biomass and Carbon Dataset for the year 2000 (NBCD). The height-structured metrics were compared with other habitat metrics for statistical association with richness of three forest breeding bird guilds across Breeding Bird Survey (BBS) routes: a broadly grouped woodland guild, and two forest breeding guilds with preferences for forest edge and for interior forest. Parametric and non-parametric models were built to examine the improvement of predictability. Height-structured metrics had the strongest associations with species richness, yielding improved predictive ability for the woodland guild richness models (r(2) = ∼ 0.53 for the parametric models, 0.63 the non-parametric models) and the forest edge guild models (r(2) = ∼ 0.34 for the parametric models, 0.47 the non-parametric models). All but one of the linear models incorporating height-structured metrics showed significantly higher adjusted-r2 values than their counterparts without additional metrics. The interior forest guild richness showed a consistent low association with height-structured metrics. Our results suggest that height heterogeneity, beyond canopy height alone, supplements habitat characterization and richness models of forest bird species. The metrics and models derived in this study demonstrate practical examples of utilizing three-dimensional vegetation data for improved characterization of spatial patterns in species richness.

Estimation of forest aboveground biomass and uncertainties by integration of field measurements, airborne LiDAR, and SAR and optical satellite data in Mexico

BACKGROUND

Information on the spatial distribution of aboveground biomass (AGB) over large areas is needed for understanding and managing processes involved in the carbon cycle and supporting international policies for climate change mitigation and adaption. Furthermore, these products provide important baseline data for the development of sustainable management strategies to local stakeholders. The use of remote sensing data can provide spatially explicit information of AGB from local to global scales. In this study, we mapped national Mexican forest AGB using satellite remote sensing data and a machine learning approach. We modelled AGB using two scenarios: (1) extensive national forest inventory (NFI), and (2) airborne Light Detection and Ranging (LiDAR) as reference data. Finally, we propagated uncertainties from field measurements to LiDAR-derived AGB and to the national wall-to-wall forest AGB map.

RESULTS

The estimated AGB maps (NFI- and LiDAR-calibrated) showed similar goodness-of-fit statistics (R², Root Mean Square Error (RMSE)) at three different scales compared to the independent validation data set. We observed different spatial patterns of AGB in tropical dense forests, where no or limited number of NFI data were available, with higher AGB values in the LiDAR-calibrated map. We estimated much higher uncertainties in the AGB maps based on two-stage up-scaling method (i.e., from field measurements to LiDAR and from LiDAR-based estimates to satellite imagery) compared to the traditional field to satellite up-scaling. By removing LiDAR-based AGB pixels with high uncertainties, it was possible to estimate national forest AGB with similar uncertainties as calibrated with NFI data only.

CONCLUSIONS

Since LiDAR data can be acquired much faster and for much larger areas compared to field inventory data, LiDAR is attractive for repetitive large scale AGB mapping. In this study, we showed that two-stage up-scaling methods for AGB estimation over large areas need to be analyzed and validated with great care. The uncertainties in the LiDAR-estimated AGB propagate further in the wall-to-wall map and can be up to 150%. Thus, when a two-stage up-scaling method is applied, it is crucial to characterize the uncertainties at all stages in order to generate robust results. Considering the findings mentioned above LiDAR can be used as an extension to NFI for example for areas that are difficult or not possible to access.

From small-scale forest structure to Amazon-wide carbon estimates

Tropical forests play an important role in the global carbon cycle. High-resolution remote sensing techniques, e.g., spaceborne lidar, can measure complex tropical forest structures, but it remains a challenge how to interpret such information for the assessment of forest biomass and productivity. Here, we develop an approach to estimate basal area, aboveground biomass and productivity within Amazonia by matching 770,000 GLAS lidar (ICESat) profiles with forest simulations considering spatial heterogeneous environmental and ecological conditions. This allows for deriving frequency distributions of key forest attributes for the entire Amazon. This detailed interpretation of remote sensing data improves estimates of forest attributes by 20–43% as compared to (conventional) estimates using mean canopy height. The inclusion of forest modeling has a high potential to close a missing link between remote sensing measurements and the 3D structure of forests, and may thereby improve continent-wide estimates of biomass and productivity.

Improving estimates of forest biomass based on remote sensing data is important to assess global carbon cycling. Here the authors develop an approach to use forest gap models to simulate lidar waveforms and compare the outputs with ICESAT-1 GLAS profiles, showing improved estimates across the Amazon basin.

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.

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.

Integrating forest inventory and analysis data into a LIDAR-based carbon monitoring system

BACKGROUND

Forest Inventory and Analysis (FIA) data may be a valuable component of a LIDAR-based carbon monitoring system, but integration of the two observation systems is not without challenges. To explore integration methods, two wall-to-wall LIDAR-derived biomass maps were compared to FIA data at both the plot and county levels in Anne Arundel and Howard Counties in Maryland. Allometric model-related errors were also considered.

RESULTS

In areas of medium to dense biomass, the FIA data were valuable for evaluating map accuracy by comparing plot biomass to pixel values. However, at plots that were defined as "nonforest", FIA plots had limited value because tree data was not collected even though trees may be present. When the FIA data were combined with a previous inventory that included sampling of nonforest plots, 21 to 27% of the total biomass of all trees was accounted for in nonforest conditions, resulting in a more accurate benchmark for comparing to total biomass derived from the LIDAR maps. Allometric model error was relatively small, but there was as much as 31% difference in mean biomass based on local diameter-based equations compared to regional volume-based equations, suggesting that the choice of allometric model is important.

CONCLUSIONS

To be successfully integrated with LIDAR, FIA sampling would need to be enhanced to include measurements of all trees in a landscape, not just those on land defined as "forest". Improved GPS accuracy of plot locations, intensifying data collection in small areas with few FIA plots, and other enhancements are also recommended.

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.

Integrating LIDAR and forest inventories to fill the trees outside forests data gap

Forest inventories are commonly used to estimate total tree biomass of forest land even though they are not traditionally designed to measure biomass of trees outside forests (TOF). The consequence may be an inaccurate representation of all of the aboveground biomass, which propagates error to the outputs of spatial and process models that rely on the inventory data. An ideal approach to fill this data gap would be to integrate TOF measurements within a traditional forest inventory for a parsimonious estimate of total tree biomass. In this study, Light Detection and Ranging (LIDAR) data were used to predict biomass of TOF in all "nonforest" Forest Inventory and Analysis (FIA) plots in the state of Maryland. To validate the LIDAR-based biomass predictions, a field crew was sent to measure TOF on nonforest plots in three Maryland counties, revealing close agreement at both the plot and county scales between the two estimates. Total tree biomass in Maryland increased by 25.5 Tg, or 15.6%, when biomass of TOF were included. In two counties (Carroll and Howard), there was a 47% increase. In contrast, counties located further away from the interstate highway corridor showed only a modest increase in biomass when TOF were added because nonforest conditions were less common in those areas. The advantage of this approach for estimating biomass of TOF is that it is compatible with, and explicitly separates TOF biomass from, forest biomass already measured by FIA crews. By predicting biomass of TOF at actual FIA plots, this approach is directly compatible with traditionally reported FIA forest biomass, providing a framework for other states to follow, and should improve carbon reporting and modeling activities in Maryland.

The effectiveness of global protected areas for climate change mitigation

Forests play a critical role in stabilizing Earth's climate. Establishing protected areas (PAs) represents one approach to forest conservation, but PAs were rarely created to mitigate climate change. The global impact of PAs on the carbon cycle has not previously been quantified due to a lack of accurate global-scale carbon stock maps. Here we used ~412 million lidar samples from NASA's GEDI mission to estimate a total PA aboveground carbon (C) stock of 61.43 Gt (+/- 0.31), 26% of all mapped terrestrial woody C. Of this total, 9.65 + /- 0.88 Gt of additional carbon was attributed to PA status. These higher C stocks are primarily from avoided emissions from deforestation and degradation in PAs compared to unprotected forests. This total is roughly equivalent to one year of annual global fossil fuel emissions. These results underscore the importance of conservation of high biomass forests for avoiding carbon emissions and preserving future sequestration.

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.

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.

Animal soundscapes reveal key markers of Amazon forest degradation from fire and logging

Safeguarding tropical forest biodiversity requires solutions for monitoring ecosystem structure over time. In the Amazon, logging and fire reduce forest carbon stocks and alter habitat, but the long-term consequences for wildlife remain unclear, especially for lesser-known taxa. Here, we combined multiday acoustic surveys, airborne lidar, and satellite time series covering logged and burned forests (n = 39) in the southern Brazilian Amazon to identify acoustic markers of forest degradation. Our findings contradict expectations from the Acoustic Niche Hypothesis that animal communities in more degraded habitats occupy fewer “acoustic niches” defined by time and frequency. Instead, we found that aboveground biomass was not a consistent proxy for acoustic biodiversity due to the divergent patterns of “acoustic space occupancy” between logged and burned forests. Ecosystem soundscapes highlighted a stark, and sustained reorganization in acoustic community assembly after multiple fires; animal communication networks were quieter, more homogenous, and less acoustically integrated in forests burned multiple times than in logged or once-burned forests. These findings demonstrate strong biodiversity cobenefits from protecting burned Amazon forests from recurrent fire. By contrast, soundscape changes after logging were subtle and more consistent with acoustic community recovery than reassembly. In both logged and burned forests, insects were the dominant acoustic markers of degradation, particularly during midday and nighttime hours, which are not typically sampled by traditional biodiversity field surveys. The acoustic fingerprints of degradation history were conserved across replicate recording locations, indicating that soundscapes may offer a robust, taxonomically inclusive solution for digitally tracking changes in acoustic community composition over time.

Characterizing the structural complexity of the Earth's forests with spaceborne lidar

Forest structural complexity is a key element of ecosystem functioning, impacting light environments, nutrient cycling, biodiversity, and habitat quality. Addressing the need for a comprehensive global assessment of actual forest structural complexity, we derive a near-global map of 3D canopy complexity using data from the GEDI spaceborne lidar mission. These data show that tropical forests harbor most of the high complexity observations, while less than 20% of temperate forests reached median levels of tropical complexity. Structural complexity in tropical forests is more strongly related to canopy attributes from lower and middle waveform layers, whereas in temperate forests upper and middle layers are more influential. Globally, forests exhibit robust scaling relationships between complexity and canopy height, but these vary geographically and by biome. Our results offer insights into the spatial distribution of forest structural complexity and emphasize the importance of considering biome-specific and fine-scale variations for ecological research and management applications. The GEDI Waveform Structural Complexity Index data product, derived from our analyses, provides researchers and conservationists with a single, easily interpretable metric by combining various aspects of canopy structure.