The Global Ecosystem Dynamics Investigation LiDAR (GEDI) is a new full waveform (FW) based LiDAR system that presents a new opportunity for the observation of forest structures globally. The backscattered GEDI signals, as all FW systems, are distorted by topographic conditions within their footprint, leading to uncertainties on the measured forest variables. In this study, we explore how well several approaches based on waveform metrics and ancillary digital elevation model (DEM) data perform on the estimation of stand dominant heights ($H_{dom}$) and wood volume (V) across different sites of Eucalyptus plantations with varying terrain slopes. In total, five models were assessed on their ability to estimate $H_{dom}$ and four models for V. Results showed that the models using the GEDI metrics, such as the height at different energy quantiles with terrain data from the shuttle radar topography mission’s (SRTM) digital elevation model (DEM) were still dependent on the topographic slope. For $H_{dom}$, an RMSE increase of 14% was observed for data acquired over slopes higher than 20% in comparison to slopes between 10 and 20%. For V, a 74% increase in RMSE was reported between GEDI data acquired over slopes between 0–10% and those acquired over slopes higher than 10%. Next, a model relying on the height at different energy quantiles of the entire waveform ($H{T}_n$) and the height at different energy quartiles of the bare ground waveform ($H{G}_n$) was assessed. Two sets of the $H{G}_n$ metrics were generated, the first one was obtained using a simulated waveform representing the echo from a bare ground, while the second one relied on the actual ground return from the waveform by means of Gaussian fitting. Results showed that both the simulated and fitted models provide the most accurate estimates of $H_{dom}$ and V for all slope ranges. The simulation-based model showed an RMSE that ranged between 1.39 and 1.66 m (between 26.76 and 39.26 m³·ha⁻¹ for V) while the fitting-based method showed an RMSE that ranged between 1.26 and 1.34 m (between 26.78 and 36.29 m³·ha⁻¹ for V). Moreover, the dependency of the GEDI metrics on slopes was greatly reduced using the two sets of metrics. As a conclusion, the effect of slopes on the 25-m GEDI footprints is rather low as the estimation on canopy heights from uncorrected waveforms degraded by a maximum of 1 m for slopes between 20 and 45%. Concerning the wood volume estimation, the effect of slopes was more pronounced, and a degradation on the accuracy (increased RMSE) of a maximum of 20 m³·ha⁻¹ was observed for slopes between 20 and 45%.

Radar altimetry is now commonly used to provide long-term monitoring of inland water levels in complement to or for replacing disappearing in situ networks of gauge stations. Recent improvements in tracking and acquisition modes improved the quality the water retrievals. The newly implemented Open Loop mode is likely to increase the number of monitored water bodies owing to the use of an a priori elevation, especially in hilly and mountainous areas. The novelty of this study is to provide a comprehensive evaluation of the performances of the past and current radar altimetry missions according to their acquisition (Low Resolution Mode or Synthetic Aperture Radar) and tracking (close or open loop) modes, and acquisition frequency (Ku or Ka) in a mountainous area where tracking losses of the signal are likely to occur, as well as of the recently launched ICESat-2 and GEDI lidar missions. To do so, we evaluate the quality of water level retrievals from most radar altimetry missions launched after 1995 over eight lakes in Switzerland, using the recently developed ALtimetry Time Series software, to compare the performances of the new tracking and acquisition modes and also the impact of the frequency used. The combination of the Open Loop tracking mode with the Synthetic Aperture Radar acquisition mode on SENTINEL-3A and B missions outperforms the classical Low Resolution Mode of the other missions with a lake observability greater than 95%, an almost constant bias of (−0.17 ± 0.04) m, a RMSE generally lower than 0.07 m and a R most of the times higher than 0.85 when compared to in situ gauge records. To increase the number of lakes that can be monitored and the temporal sampling of the water level retrievals, data acquired by lidar altimetry missions were also considered. Very accurate results were also obtained with ICESat-2 data with RMSE lower than 0.06 and R higher than 0.95 when compared to in situ water levels. An almost constant bias (0.42 ± 0.03) m was also observed. More contrasted results were obtained using GEDI. As these data were available on a shorter time period, more analyses are necessary to determine their potential for retrieving water levels.

Estimating forest canopy height from large-footprint satellite LiDAR waveforms is challenging given the complex interaction between LiDAR waveforms, terrain, and vegetation, especially in dense tropical and equatorial forests. In this study, canopy height in French Guiana was estimated using multiple linear regression models and the Random Forest technique (RF). This analysis was either based on LiDAR waveform metrics extracted from the GLAS (Geoscience Laser Altimeter System) spaceborne LiDAR data and terrain information derived from the SRTM (Shuttle Radar Topography Mission) DEM (Digital Elevation Model) or on Principal Component Analysis (PCA) of GLAS waveforms. Results show that the best statistical model for estimating forest height based on waveform metrics and digital elevation data is a linear regression of waveform extent, trailing edge extent, and terrain index (RMSE of 3.7 m). For the PCA based models, better canopy height estimation results were observed using a regression model that incorporated both the first 13 principal components (PCs) and the waveform extent (RMSE = 3.8 m). Random Forest regressions revealed that the best configuration for canopy height estimation used all the following metrics: waveform extent, leading edge, trailing edge, and terrain index (RMSE = 3.4 m). Waveform extent was the variable that best explained canopy height, with an importance factor almost three times higher than those for the other three metrics (leading edge, trailing edge, and terrain index). Furthermore, the Random Forest regression incorporating the first 13 PCs and the waveform extent had a slightly-improved canopy height estimation in comparison to the linear model, with an RMSE of 3.6 m. In conclusion, multiple linear regressions and RF regressions provided canopy height estimations with similar precision using either LiDAR metrics or PCs. However, a regression model (linear regression or RF) based on the PCA of waveform samples with waveform extent information is an interesting alternative for canopy height estimation as it does not require several metrics that are difficult to derive from GLAS waveforms in dense forests, such as those in French Guiana.

LiDAR technology has been widely used to characterize structural parameters of forest ecosystems, which in turn are valuable information for forest monitoring. GEDI is a spaceborne LiDAR system specifically designed to measure vegetation’s vertical structure, and it has been acquiring waveforms on a global scale since April 2019. In particular, canopy height is an important descriptor of forest ecosystems, as it allows for quantifying biomass and other inventory information. This paper analyzes the accuracy of canopy height estimates from GEDI data over tropical forests in French Guiana and Gabon. The influence of various signal acquisition and processing parameters is assessed to highlight how they impact the estimation of canopy heights. Canopy height models derived from airborne LiDAR data are used as reference heights. Several linear and non-linear approaches are tested given the richness of the available GEDI information. The results show that the use of regression models built on multiple GEDI metrics allows for reaching improved accuracies compared to a direct estimation from a single GEDI height metric. In a notable way, random forest improves the canopy height estimation accuracy by almost 80% (in terms of RMSE) compared to the use of rh_95 as a direct proxy of canopy height. Additionally, convolutional neural networks calibrated on GEDI waveforms exhibit similar results to the ones of other regression models. Beam type as well as beam sensitivity, which are related to laser penetration, appear as parameters of major influence on the data derived from GEDI waveforms and used as input for canopy height estimation. Therefore, we recommend the use of only power and high-sensitivity beams when sufficient data are available. Finally, we note that regression models trained on reference data can be transferred across study sites that share identical environmental conditions.

In the context of monitoring and assessment of water consumption in the agricultural sector, the objective of this study is to build an operational approach capable of detecting irrigation events at plot scale in a near real-time scenario using Sentinel-1 (S1) data. The proposed approach is a decision tree-based method relying on the change detection in the S1 backscattering coefficients at plot scale. First, the behavior of the S1 backscattering coefficients following irrigation events has been analyzed at plot scale over three study sites located in Montpellier (southeast France), Tarbes (southwest France), and Catalonia (northeast Spain). To eliminate the uncertainty between rainfall and irrigation, the S1 synthetic aperture radar (SAR) signal and the soil moisture estimations at grid scale (10 km × 10 km) have been used. Then, a tree-like approach has been constructed to detect irrigation events at each S1 date considering additional filters to reduce ambiguities due to vegetation development linked to the growth cycle of different crops types as well as the soil surface roughness. To enhance the detection of irrigation events, a filter using the normalized differential vegetation index (NDVI) obtained from Sentinel-2 optical images has been proposed. Over the three study sites, the proposed method was applied on all possible S1 acquisitions in ascending and descending modes. The results show that 84.8% of the irrigation events occurring over agricultural plots in Montpellier have been correctly detected using the proposed method. Over the Catalonian site, the use of the ascending and descending SAR acquisition modes shows that 90.2% of the non-irrigated plots encountered no detected irrigation events whereas 72.4% of the irrigated plots had one and more detected irrigation events. Results over Catalonia also show that the proposed method allows the discrimination between irrigated and non-irrigated plots with an overall accuracy of 85.9%. In Tarbes, the analysis shows that irrigation events could still be detected even in the presence of abundant rainfall events during the summer season where two and more irrigation events have been detected for 90% of the irrigated plots. The novelty of the proposed method resides in building an effective unsupervised tool for near real-time detection of irrigation events at plot scale independent of the studied geographical context.

LiDAR data has been successfully used to estimate forest parameters such as canopy heights and biomass. Major limitation of LiDAR systems (airborne and spaceborne) arises from their limited spatial coverage. In this study, we present a technique for canopy height mapping using airborne and spaceborne LiDAR data (from the Geoscience Laser Altimeter System (GLAS)). First, canopy heights extracted from both airborne and spaceborne LiDAR were extrapolated from available environmental data. The estimated canopy height maps using Random Forest (RF) regression from airborne or GLAS calibration datasets showed similar precisions (~6 m). To improve the precision of canopy height estimates, regression-kriging was used. Results indicated an improvement in terms of root mean square error (RMSE, from 6.5 to 4.2 m) using the GLAS dataset, and from 5.8 to 1.8 m using the airborne LiDAR dataset. Finally, in order to investigate the impact of the spatial sampling of future LiDAR missions on canopy height estimates precision, six subsets were derived from the initial airborne LiDAR dataset. Results indicated that using the regression-kriging approach a precision of 1.8 m on the canopy height map was achievable with a flight line spacing of 5 km. This precision decreased to 4.8 m for flight line spacing of 50 km.