Spaceborne LiDAR reveals 3D structural differences between natural forests and tree plantations in China

Accurate information on the global distribution and the three-dimensional (3D) structure of Earth’s forests is needed to assess forest biomass stocks and to project the future of the terrestrial Carbon sink. In spite of its importance, the 3D structure of forests continues to be the most crucial information gap in the observational archive. The Global Ecosystem Dynamics Investigation (GEDI) Light Detection and Ranging (LiDAR) sensor is providing an unprecedented near-global sampling of tropical and temperate forest structural properties. The integration of GEDI measurements with spatially-contiguous observations from polar orbiting optical satellite data therefore provides a unique opportunity to produce wall-to-wall maps of forests’ 3D structure. Here, we utilized Visible Infrared Imaging Radiometer Suite (VIIRS) annual metrics data to extrapolate GEDI-derived forest structure attributes into 1-km resolution contiguous maps of tree height (TH), canopy fraction cover (CFC), plant area index (PAI), and foliage height diversity (FHD) for the conterminous US (CONUS). The maps were validated using an independent subset of GEDI data. Validation results for TH (r2 = 0.8; RMSE = 3.35 m), CFC (r2 = 0.79; RMSE = 0.09), PAI (r2 = 0.76; RMSE = 0.41), and FHD (r2 = 0.83; RMSE = 0.25) demonstrated the robustness of VIIRS data for extrapolating GEDI measurements across the nation or even over larger areas. The methodology developed through this study may allow multi-decadal monitoring of changes in multiple forest structural attributes using consistent satellite observations acquired by orbiting and forthcoming VIIRS instruments.

Monitoring forest structural properties is critical for a range of applications because structure is key to understanding and quantifying forest biophysical functioning, including stand dynamics, evapotranspiration, habitat, and recovery from disturbances. Monitoring of forest structural properties at desirable frequencies and cost globally is enabled by space-borne LiDAR missions such as the global ecosystem dynamics investigation (GEDI) mission. This study assessed the accuracy of GEDI estimates for canopy height, total plant area index (PAI), and vertical profile of plant area volume density (PAVD) and elevation over a gradient of canopy height and terrain slope, compared to estimates derived from airborne laser scanning (ALS) across two forest age-classes in the Central Highlands region of south-eastern Australia. ALS was used as a reference dataset for validation of GEDI (Version 2) dataset. Canopy height and total PAI analyses were carried out at the landscape level to understand the influence of beam-type, height of the canopy, and terrain slope. An assessment of GEDI’s terrain elevation accuracy was also carried out at the landscape level. The PAVD profile evaluation was carried out using footprints grouped into two forest age-classes, based on the areas of mountain ash (Eucalyptus regnans) forest burnt in the Central Highlands during the 1939 and 2009 wildfires. The results indicate that although GEDI is found to significantly under-estimate the total PAI and slightly over-estimate the canopy height, the GEDI estimates of canopy height and the vertical PAVD profile (above 25 m) show a good level of accuracy. Both beam-types had comparable accuracies, with increasing slope having a slightly detrimental effect on accuracy. The elevation accuracy of GEDI found the RMSE to be 10.58 m and bias to be 1.28 m, with an R2 of 1.00. The results showed GEDI is suitable for canopy densities and height in complex forests of south-eastern Australia.

The impact of geolocation uncertainty on GEDI tropical forest canopy height estimation and change monitoring

Invasive species are a leading threat to biodiversity worldwide. Nypa palm (Nypa fruticans) has emerged as the predominant invasive species in the Niger Delta region of Nigeria. While endemic mangroves have high rates of carbon sequestration, stabilize coastlines and protect biodiversity, Nypa does not provide these services outside its native region of Southeast Asia. Oil exploration and urbanization in this region also exacerbate mangrove loss and Nypa spread. As Nypa is difficult to distinguish from endemic mangrove species in remotely sensed data, estimates of mangrove and ecosystem services losses in Nigeria are highly uncertain. Here, we analyse multisensor satellite data with machine learning to quantify the rapid expansion of Nypa from 2015 to 2020 in Nigeria. Using Landsat imagery and random forest classification, we quantify total potential Nypa extent in Nigeria in 2019. We then produced a Nypa extent map using iterative combinations of Sentinel‐1 SAR, Sentinel‐2 MSI and ALOS PALSAR. Random forest classifications using SAR data from ALOS and Sentinel‐1 were best suited for mapping Nypa extent with similar accuracies (78% and 75%, respectively). Based on data availability and accuracy, we focussed our change analysis on Sentinel‐1 SAR. Our results show ~28 000 ha of mangroves were converted to Nypa in Nigeria by 2020 and covered a larger extent than endemic mangroves, compounding the effect of the existing degradation and deforestation in the region. We also compared forest height and complexity estimates from Global Ecosystem Dynamics Investigation LiDAR to further distinguish between endemic mangroves and Nypa in three dimensions. Nypa structural variability, measured by top‐of‐canopy height, vegetation cover, plant area index, and foliage height diversity, was lower than that of mangroves. At current rates of Nypa expansion, the entire area of study would be invaded by Nypa by 2028, with potentially detrimental consequences to the ecosystem services provided by mangroves.

Savannas have complex, discontinuous woody vegetation structures that vary greatly in vertical and spatial arrangement and change due to climatic, ecological and management impacts. While airborne laser scanning (ALS) data have provided detailed information on vertical vegetation structure and is widely used in ecological studies, it is lacking in availability and repeat frequency. Although the Global Ecosystem Dynamics Investigation (GEDI) waveform Light Detection and Ranging (LiDAR) sensor and algorithms were optimized for measuring dense forests, it was anticipated that GEDI metrics could provide useful characterization of lower stature, sparse savannas structures. This study provided the first baseline validation of Version 2 GEDI (L2A) relative height 98 (RH98) by comparing the on-orbit GEDI-RH98orb to the simulated GEDI-RH98sim derived from ALS data across diverse savanna vegetation. It furthermore determined the influence of various factors on error, e.g. algorithm setting group (SGs), beam type, day vs. night, beam sensitivity, and vegetation phenology. After applying quality flags, 22,813 GEDI footprints were analyzed across 11 sites. SGs 4–6 that are aimed at dense forests had much larger errors than SGs 1–3. The phenological conditions at the time of GEDI data acquisition had a very large influence on the error of RH98orb. During leaf-on conditions for savanna vegetation with RH98sim < 15 m, RH98orb was very accurate with R2 = 0.61, mean bias = −0.55 m, %bias = −11.1%, RMSE = 1.64 m and %RMSE = 29.8%. In leaf-off conditions where RH98sim < 15 m, RH98orb was less accurate with R2 = 0.43, mean bias = −1.47 m, %bias = −26.5%, RMSE = 2.03 m and %RMSE = 40.9%. During leaf-off conditions, the GEDI LiDAR signal at the start of the waveform may be weaker as it interacts with denuded branches and may be truncated as noise, leading to a large negative height bias. Therefore, assessments of deciduous vegetation structures should be conducted during leaf-on periods. In leaf-on conditions, GEDI's RH98orb was very accurate between canopy heights of 3 and 7 m, with a mean bias of −0.79 m (−10%). The bias of RH98orb was not influenced by canopy cover. Due to the GEDI LiDAR pulse width of 15.6 ns, the GEDI-RH98 data product cannot reliably estimate canopy heights of shrubs below 2.34 m and will require more complex deconvolution of the waveform. GEDI's RH98 accurately estimates the canopy height of trees between 3 and 15 m allowing assessment of canopy heights over vast savanna areas.

Continuous large-scale mapping of forest canopy height is crucial for estimating and reporting forest carbon content, analyzing forest degradation and restoration, or to model ecosystem variables such as aboveground biomass. Over the last years, the spaceborne Light Detection and Ranging (LiDAR) sensor specifically designed to acquire forest structure information, Global Ecosystem Dynamics Investigation (GEDI), has been used to extract forest canopy height information over large areas. Yet, GEDI has no spatial coverage for most forested areas in Canada and other high latitude regions. On the other hand, the spaceborne LiDAR called Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) provides a global coverage but was not specially developed to study forested ecosystems. Nonetheless, both spaceborne LiDAR sensors obtain point-based information, making spatially continuous forest canopy height estimation very challenging. This study compared the performance of both spaceborne LiDAR, GEDI and ICESat-2, combined with ALOS-2/PALSAR-2 and Sentinel-1 and -2 data to produce continuous canopy height maps in Canada for the year 2020. A set-aside dataset and airborne LiDAR (ALS) from a national LiDAR campaign were used for accuracy assessment. Both maps overestimated canopy height in relation to ALS data, but GEDI had a better performance than ICESat-2 with a mean difference (MD) of 0.9 m and 2.9 m, and a root mean square error (RMSE) of 4.2 m and 5.2 m, respectively. However, as both GEDI and ALS have no coverage in most of the hemi-boreal forests, ICESat-2 captures the tall canopy heights expected for these forests better than GEDI. PALSAR-2 HV polarization was the most important covariate to predict canopy height, showing the great potential of L-band in comparison to C-band from Sentinel-1 or optical data from Sentinel-2. The approach proposed here can be used operationally to produce annual canopy height maps for areas that lack GEDI and ICESat-2 coverage.

Identification of forest fuels is a key step for forest fire prevention since they provide valuable information of fire behavior. This study assesses NASA’s Global Ecosystem Dynamics Investigation (GEDI) system to classify fuel types in Mediterranean environments according to the Prometheus model in a forested area of NE Spain. We used 59,554 GEDI footprints and extracted variables related to height metrics, canopy profile metrics, and aboveground biomass density estimates from products L2A, L2B, and L4A, respectively. Four quality filters were applied to discard high uncertainty data, reducing the initial footprints to 9,703. Spectral indices from Landsat-8 OLI scenes were created to test the effect of their integration with GEDI variables on fuel types estimation. Ground-truth data were comprised of Prometheus fuel types estimated in two previous studies. Only the types that matched in each GEDI footprint in both studies were used, resulting in a final sample of 1,112 footprints. Spearman’s correlation coefficient, Kruskal-Wallis and Dunn’s tests determined the variables to be included in the classification models: the relative height at the 85th percentile, the Plant Area Index, and the Aboveground Biomass Density from GEDI, and the brightness from Landsat-8 OLI. Best performances were achieved with Random Forest (RF) and Support Vector Machine with radial kernel (SVM-R), which were lower including only GEDI variables (accuracies: RF and SVM-R = 61.54 %) than integrating the brightness from Landsat-8 OLI (accuracies: RF = 83.71 %, SVM-R = 81.90 %). These results allow validating GEDI for fuel type classification of Prometheus model, constituting a promising information for forest management over large areas.

Tracking savanna vegetation structure in South Africa by extension of GEDI canopy metrics with Landsat, Sentinel-2, and PALSAR

Leaf or plant area index (LAI, PAI) information is frequently used to describe vegetation structure in environmental science. While field measurements are time-consuming and do not scale to landscapes, model-based air- or space-borne remote-sensing methods have been used for many years for area-wide monitoring. As of 2019, NASA’s Global Ecosystem Dynamics Investigation (GEDI) mission delivers a point-based LAI product with 25 m footprints and periodical repetition. This opens up new possibilities in integrating GEDI as frequently generated training samples with high resolution (spectral) sensors. However, the foreseeable duration of the system installed on the ISS is limited. In this study we want to test the potential of GEDI for regional comprehensive LAI estimations throughout the year with a focus on its usability beyond the lifespan of the GEDI mission. We study the landscape of Hesse, Germany, with its pronounced seasonal changes. Assuming a relationship between GEDI’s PAI and Sentinel-1 and -2 data, we used a Random Forest approach together with spatial variable selection to make predictions for new Sentinel scenes. The model was trained with two years of GEDI PAI data and validated against a third year to provide a robust and temporally independent model validation. This ensures the applicability of the validation for years outside the training period, reaching a total RMSE of 1.12. Predictions for the test year showed the expected seasonal and spatial patterns indicated by RMSE values ranging between 0.75 and 1.44, depending on the land cover class. The overall prediction performance shows good agreement with the test data set of the independent year which supports our assumption that the usage of GEDI’s PAI beyond the mission lifespan is feasible for regional studies.

Forests are widely recognized as critical to combating climate change due to their ability to sequester and store carbon in the form of biomass. In recent years, the combined use of data from ground-based forest inventories and remotely sensed data from light detection and ranging (lidar) has proven useful for large-scale assessment of forest biomass, but airborne lidar is expensive and data acquisition is infeasible for many countries. By contrast, the spaceborne Global Ecosystem Dynamics Investigation (GEDI) lidar instrument has collected freely available data for most of the world’s temperate and tropical forests since 2019. GEDI’s biomass products rely on models calibrated with a global network of field plots paired with GEDI waveforms simulated from airborne lidar to predict biomass. While this calibration strategy minimizes spatial and temporal offsets between field measurements and corresponding lidar returns, calibration data are sparse in many regions. Paraguay’s forests are known to be poorly represented in GEDI’s current calibration dataset, and here we demonstrate that local models calibrated opportunistically with on-orbit GEDI data and field surveys from Paraguay’s national forest inventory can be used with GEDI’s statistical estimators of aboveground biomass density (AGBD). We specify a protocol for opportunistically matching GEDI observations with field plots to calibrate a field-to-GEDI biomass model for use in GEDI’s hybrid statistical framework. Country-specific calibration using on-orbit data resulted in relatively accurate and unbiased predictions of footprint-level biomass, and importantly, supported the assumption underlying model-based inference that the model must ‘apply’ to the area of interest. Using a locally calibrated biomass model, we estimate that the mean AGBD in Paraguay is 65.55 Mg ha−1, which coincides well with the design-based approach employed by the national forest inventory. The GEDI estimates for individual forest strata range from 52.34 Mg ha−1 to 103.88 Mg ha−1. On average, the standard errors are 47% lower for estimates based on GEDI than the forest inventory, representing a significant gain in precision. Our research demonstrates that GEDI can be used by national forest inventories in countries that seek reliable estimates of AGBD, and that local calibration using existing field plots may be more appropriate in some applications than using GEDI global models, especially in regions where those models are sparsely calibrated.

The Global Ecosystem Dynamics Investigation (GEDI), a new spaceborne LiDAR system of the National Aeronautics and Space Administration (NASA), has the potential to revolutionize global measurements of vertical vegetation structure. However, GEDI performance among different forest types and factors influencing GEDI performance needs to be evaluated against similar measurements from existing airborne LiDAR platforms. Ideally, comparisons across diverse forest types will inform future work quantifying biomass or mapping species habitats. Thus, we compared the second version of GEDI L2A product (GEDI V2) with Airborne Observation Platform (AOP) leaf-on LiDAR data across 33 National Ecological Observation Network (NEON) sites. Comparisons were made for ground elevation and relative height (RH) of GEDI with simulated airborne laser scanning (ALS) waveforms from discrete point cloud LiDAR. Results indicated that GEDI V2 obtained high accuracy on ground elevation and RH100 estimations (3σ) with RMSEs of 1.38 m and 2.62 m, respectively. GEDI produced forest height estimations (RH100) for all 12 forest types with a %RMSE below 25%. GEDI RHs were sensitive to ground finding accuracy, and GEDI performance of RH estimation varied from forest profiles of different forest types. For factors influencing GEDI performance, greater than 21% of GEDI RH95 and 33% of ground elevation variations can be explained by land surface attributes, observing sensor system characteristics, and the collection time differences between GEDI and NEON LiDAR. Furthermore, geolocation error remains an essential factor affecting GEDI performance, which varies among forest and land cover types, especially for canopy height estimation. The findings reported here can provide insights to guide and enhance future GEDI-based global forest structure mapping and applications.

Modelling forest volume with small area estimation of forest inventory using GEDI footprints as auxiliary information

Estimating biomass consumption and carbon emissions by integrating GEDI, Sentinel-1, and Sentinel-2 data

A comparative discussion of the advantages and disadvantages of natural stands and plantations, including in terms of their productivity and stability, began from the moment of the first forest plantings and continues to this day. In the context of the progressive replacement of natural forests by plantations due to deforestation, the question of how will change the carbon storage capacity of forest cover when replacing natural forests with artificial ones in a changing climate becomes extremely relevant. This article presents the first attempt to answer this question at the transcontinental level on a special case for two-needled pine trees (subgenus Pinus L.). The research was carried out using the database compiled by the authors on the single-tree biomass structure of forest-forming species of Eurasia, in particular, data of 1880 and 1967 of natural and plantation trees, respectively. Multi-factor regression models are calculated after combining the matrix of initial data on the structure of tree biomass with the mean January temperature and mean annual precipitation, and their adequacy indices allow us to consider them reproducible. It is found that the aboveground and stem biomass of equal-sized and equal-aged natural and plantation trees increases as the January temperature and precipitation rise. This pattern is only partially valid for the branches biomass, and it has a specific character for the foliage one. The biomass of all components of plantation trees is higher than that of natural trees, but the percent of this excess varies among different components and depends on the level of January temperatures, but does not depend at all on the level of annual precipitation. A number of uncertainties that arose during the modeling process, as well as the preliminary nature of the obtained regularities, are noted.

The replacement of native forests by high-density tree plantations affects the richness and composition of animal communities through the modification of the resource availability, hydrological regimes, nutrient cycles and soil structure. Previous studies have been conducted mainly on birds and mammals, whereas few have explored the response of anurans, a taxon considered at risk and in steep decline. With this study we aimed to evaluate changes in anuran diversity associated with native forest replacement by pine plantations (Pinus taeda) in the southern Atlantic forest of Argentina, a highly diverse and threatened ecosystem. Additionally, we intended to explore the role of the micro-habitat (terrestrial, aquatic) in the response. We characterized vegetation and habitat structure and sampled anurans in terrestrial and aquatic micro-habitats in the native forest and pine plantations monthly over two consecutive years (2012–2014), using three complementary techniques (pitfall traps, audio strip transects and larval sampling). A total of 964 individuals of 21 species were captured: Physalaemus cuvieri, Odontophrynus americanus and Elachistocleis bicolor were the most abundant species. Replacement of the structurally complex native forest by extremely simplified monoculture tree plantations influenced the patterns of anuran alpha and beta diversity. However, these changes were micro-habitat-dependent: changes in anuran diversity in the terrestrial habitat were explained by species loss, while those in the aquatic habitat were explained by both species loss and turnover. The arboreal species of the family Hylidae (eight) were found absent from tree plantations. Both the hydroperiod instability of the water bodies and the simplified vegetation structure of the tree plantations are probably limiting the suitability for both reproduction and larval development of some specialist species. The native vegetation surrounding the water bodies in the tree plantations constitutes a central element to maintain anuran diversity, through increasing the hydroperiod and providing a habitat for species reproduction in the southern Atlantic forest and to facilitate the movement of individuals among native forest remnants.