Remote Sensing of Forest Biomass Using GNSS Reflectometry

In this study, the capability of Global Navigation Satellite System Reflectometry in evaluating forest biomass from space has been investigated by using data coming from the TechDemoSat-1 (TDS-1) mission of Surrey Satellite Technology Ltd. and from the Cyclone Satellite System (CyGNSS) mission of NASA. The analysis has been first conducted using TDS-1 data on a local scale, by selecting five test areas located in different parts of the Earth's surface. The areas were chosen as examples of various forest coverages, including equatorial and boreal forests. Then, the analysis has been extended by using CyGNSS to a global scale, including any type of forest coverage. The peak of the Delay Doppler Map calibrated to retrieve an "equivalent" reflectivity has been exploited for this investigation and its sensitivity to forest parameters has been evaluated by a direct comparison with vegetation optical depth (VOD) derived from the Soil Moisture Active Passive L-band radiometer, with a pantropical aboveground biomass (AGB) map and then with a tree height (H) global map derived from the Geoscience Laser Altimeter System installed on-board the ICEsat satellite. The sensitivity analysis confirmed the decreasing trend of the observed equivalent reflectivity for increasing biomass, with correlation coefficients 0.31 ≤ R ≤ 0.54 depending on the target parameter (VOD, AGB, or H) and on the considered dataset (local or global). These correlations were not sufficient to retrieve the target parameters by simple inversion of the direct relationships. The retrieval has been therefore based on Artificial Neural Networks making it possible to add other inputs (e.g., the incidence angle, the signal to noise ratio, and the lat/lon information in case of global maps) to the algorithm. Although not directly correlated to the biomass, these inputs helped in improving the retrieval accuracy. The algorithm was tested on both the selected areas and globally, showing a promising ability to retrieve the target parameter, either AGB or H, with correlation coefficients R ≃ 0.8