Search

The Copernicus Sentinel-3 Surface Topography Mission (STM) Land Altimetry provides valuable surface elevation information over inland waters, sea ice, and land ice, thanks to its synthetic aperture radar (SAR) altimeter and its orbit that covers high-latitude polar regions. To ensure that these measurements are reliable and to maximise the return on investment, adequate validation of the geophysical retrieval methods, processing algorithms, and corrections must be performed using independent observations. The EU-ESA project St3TART (started July 2021) aims to generalise the concept of Fiducial Reference Measurements (FRMs) for the Copernicus Sentinel-3 STM. This work has gathered existing data, made new observations during field campaigns, and ensured that these observations meet the criteria of FRM standards so that they can be used to validate Sentinel-3 STM Land Altimetry products operationally. A roadmap for the operational provision of the FRM, including the definition, consolidation, and identification of the most relevant and cost-effective methods and protocols to be maintained, supported, or implemented, has been developed. The roadmap includes guidelines for SI traceability, definitions of FRM measurement procedures, processing methods, and uncertainty budget estimations.

Surface soil moisture (SSM) is a key variable for many environmental studies, including hydrology and agriculture. Synthetic aperture radar (SAR) data in the C-band are widely used nowadays to estimate SSM since the Sentinel-1 provides free-of-charge C-band SAR images at high spatial resolution with high revisit time, whereas the use of L-band is limited due to the low data availability. In this context, the main objective of this paper is to develop an operational approach for SSM estimation that mainly uses data in the C-band (Sentinel-1) with L-bands (ALOS/PALSAR) as additional data to improve SSM estimation accuracy. The approach is based on the use of the artificial neural networks (NNs) technique to firstly derive the soil roughness (Hrms) from the L-band (HH polarization) to then consider the L-band-derived Hrms and C-band SAR data (VV and VH polarizations) in the input vectors of NNs for SSM estimation. Thus, the Hrms estimated from the L-band at a given date is assumed to be constant for a given times series of C-band images. The NNs were trained and validated using synthetic and real databases. The results showed that the use of the L-band-derived Hrms in the input vector of NN in addition to C-band SAR data improved SSM estimation by decreasing the error (bias and RMSE), mainly for SSM values lower than 15 vol.% and regardless of Hrms values. Based on the synthetic database, the NNs that neglect the Hrms underestimate and overestimate the SSM (bias ranges between −8.0 and 4.0 vol.%) for Hrms values lower and higher than 1.5 cm, respectively. For Hrms 2.0 cm and SSM between 8 to 22 vol.%, the accuracy on the SSM estimation improved and the overestimation decreased by 2.2 vol.% (from 4.5 to 2.3 vol.%). From the real database, the use of Hrms estimated from the L-band brought a significant improvement of the SSM estimation accuracy. For in situ SSM less than 15 vol.%, the RMSE decreased by 1.5 and 2.2 vol.% and the bias by 1.2 and 2.6 vol.%, for Hrms values lower and higher than 1.5 cm, respectively.

Although the real timing and flow rates used for crop irrigation are controlled at the scale of individual plots by the irrigator, they are not generally known by the farm upper management. This information is nevertheless essential, not only to compute the water balance of irrigated plots and to schedule irrigation, but also for the management of water resources at regional scales. The aim of the present study was to detect irrigation timing using time series of surface soil moisture (SSM) derived from Sentinel-1 radar observations. The method consisted of assessing the direction of change of surface soil moisture (SSM) between observations and a water balance model, and to use thresholds to be calibrated. The performance of the approach was assessed on the F-score quantifying the accuracy of the irrigation event detections and ranging from 0 (none of the irrigation timing is correct) to 100 (perfect irrigation detection). The study focused on five irrigated and one rainfed plot of maize in South-West France, where the approach was tested using in situ measurements and surface soil moisture (SSM) maps derived from Sentinel-1 radar data. The use of in situ data showed that (1) irrigation timing was detected with a good accuracy (F-score in the range (80–83) for all plots) and (2) the optimal revisit time between two SSM observations was 2–4 days. The higher uncertainties of microwave SSM products, especially when the crop is well developed (normalized difference of vegetation index (NDVI) > 0.7), degraded the score (F-score = 69), but various possibilities of improvement were discussed. This paper opens perspectives for the irrigation detection at the plot scale over large areas and thus for the improvement of irrigation water management.

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.

The objective of this study was to analyze the sensitivity of radar signals in the X-band in irrigated grassland conditions. The backscattered radar signals were analyzed according to soil moisture and vegetation parameters using linear regression models. A time series of radar (TerraSAR-X and COSMO-SkyMed) and optical (SPOT and LANDSAT) images was acquired at a high temporal frequency in 2013 over a small agricultural region in southeastern France. Ground measurements were conducted simultaneously with the satellite data acquisitions during several grassland growing cycles to monitor the evolution of the soil and vegetation characteristics. The comparison between the Normalized Difference Vegetation Index (NDVI) computed from optical images and the in situ Leaf Area Index (LAI) showed a logarithmic relationship with a greater scattering for the dates corresponding to vegetation well developed before the harvest. The correlation between the NDVI and the vegetation parameters (LAI, vegetation height, biomass, and vegetation water content) was high at the beginning of the growth cycle. This correlation became insensitive at a certain threshold corresponding to high vegetation (LAI ~2.5 m2/m2). Results showed that the radar signal depends on variations in soil moisture, with a higher sensitivity to soil moisture for biomass lower than 1 kg/m². HH and HV polarizations had approximately similar sensitivities to soil moisture. The penetration depth of the radar wave in the X-band was high, even for dense and high vegetation; flooded areas were visible in the images with higher detection potential in HH polarization than in HV polarization, even for vegetation heights reaching 1 m. Lower sensitivity was observed at the X-band between the radar signal and the vegetation parameters with very limited potential of the X-band to monitor grassland growth. These results showed that it is possible to track gravity irrigation and soil moisture variations from SAR X-band images acquired at high spatial resolution (an incidence angle near 30°).

The ability of Synthetic Aperture Radar (SAR) Sentinel-1 data to detect the main wheat phenological phases was investigated in the Bekaa plain of Lebanon. Accordingly, the temporal variation of Sentinel-1 (S1) signal was analyzed as a function of the phenological phases’ dates observed in situ (germination; heading and soft dough), and harvesting. Results showed that S1 data, unlike the Normalized Difference Vegetation Index (NDVI) data, were able to estimate the dates of theses phenological phases due to significant variations in S1 temporal series at the dates of germination, heading, soft dough, and harvesting. Particularly, the ratio VV/VH at low incidence angle (32−34°) was able to detect the germination and harvesting dates. VV polarization at low incidence angle (32−34°) was able to detect the heading phase, while VH polarization at high incidence angle (43−45°) was better than that at low incidence angle (32−34°), in detecting the soft dough phase. An automated approach for main wheat phenological phases’ determination was then developed on the western part of the Bekaa plain. This approach modelled the S1 SAR temporal series by smoothing and fitting the temporal series with Gaussian functions (up to three Gaussians) allowing thus to automatically detect the main wheat phenological phases from the sum of these Gaussians. To test its robustness, the automated method was applied on the northern part of the Bekaa plain, in which winter wheat is harvested usually earlier because of the different weather conditions. The Root Mean Square Error (RMSE) of the estimation of the phenological phases’ dates was 2.9 days for germination, 5.5 days for heading, 5.1 days soft dough, 3.0 days for West Bekaa’s harvesting, and 4.5 days for North Bekaa’s harvesting. In addition, a slight underestimation was observed for germination and heading of West Bekaa (−0.2 and −1.1 days, respectively) while an overestimation was observed for soft dough of West Bekaa and harvesting for both West and North Bekaa (3.1, 0.6, and 3.6 days, respectively). These results are encouraging, and thus prove that S1 data are powerful as a tool for crop monitoring, to serve enhanced crop management and production handling.

Mapping irrigated plots is essential for better water resource management. Today, the free and open access Sentinel-1 (S1) and Sentinel-2 (S2) data with high revisit time offers a powerful tool for irrigation mapping at plot scale. Up to date, few studies have used S1 and S2 data to provide approaches for mapping irrigated plots. This study proposes a method to map irrigated plots using S1 SAR (synthetic aperture radar) time series. First, a dense temporal series of S1 backscattering coefficients were obtained at plot scale in VV (Vertical-Vertical) and VH (Vertical-Horizontal) polarizations over a study site located in Catalonia, Spain. In order to remove the ambiguity between rainfall and irrigation events, the S1 signal obtained at plot scale was used conjointly to S1 signal obtained at a grid scale (10 km × 10 km). Later, two mathematical transformations, including the principal component analysis (PCA) and the wavelet transformation (WT), were applied to the several SAR temporal series obtained in both VV and VH polarization. Irrigated areas were then classified using the principal component (PC) dimensions and the WT coefficients in two different random forest (RF) classifiers. Another classification approach using one dimensional convolutional neural network (CNN) was also performed on the obtained S1 temporal series. The results derived from the RF classifiers with S1 data show high overall accuracy using the PC values (90.7%) and the WT coefficients (89.1%). By applying the CNN approach on SAR data, a significant overall accuracy of 94.1% was obtained. The potential of optical images to map irrigated areas by the mean of a normalized differential vegetation index (NDVI) temporal series was also tested in this study in both the RF and the CNN approaches. The overall accuracy obtained using the NDVI in RF classifier reached 89.5% while that in the CNN reached 91.6%. The combined use of optical and radar data slightly enhanced the classification in the RF classifier but did not significantly change the accuracy obtained in the CNN approach using S1 data.

This study proposes an effective method to map rice crops using the Sentinel-1 SAR (Synthetic Aperture Radar) time series over the Camargue region, Southern France. First, the temporal behavior of the SAR backscattering coefficient over 832 plots containing different crop types was analyzed. Through this analysis, the rice cultivation was identified using metrics derived from the Gaussian profile of the VV/VH time series (3 metrics), the variance of the VV/VH time series (one metric), and the slope of the linear regression of the VH time series (one metric). Using the derived metrics, rice plots were mapped through two different approaches: decision tree and Random Forest (RF). To validate the accuracy of each approach, the classified rice map was compared to the available national data. Similar high overall accuracy was obtained using both approaches. The overall accuracy obtained using a simple decision tree reached 96.3%, whereas an overall accuracy of 96.6% was obtained using the RF classifier. The approach, therefore, provides a simple yet precise and powerful tool to map paddy rice areas.

This paper presents a technique for the mapping of soil moisture and irrigation, at the scale of agricultural fields, based on the synergistic interpretation of multi-temporal optical and Synthetic Aperture Radar (SAR) data (Sentinel-2 and Sentinel-1). The Kairouan plain, a semi-arid region in central Tunisia (North Africa), was selected as a test area for this study. Firstly, an algorithm for the direct inversion of the Water Cloud Model (WCM) was developed for the spatialization of the soil water content between 2015 and 2017. The soil moisture retrieved from these observations was first validated using ground measurements, recorded over 20 reference fields of cereal crops. A second method, based on the use of neural networks, was also used to confirm the initial validation. The results reported here show that the soil moisture products retrieved from remotely sensed data are accurate, with a Root Mean Square Error (RMSE) of less than 5% between the two moisture products. In addition, the analysis of soil moisture and Normalized Difference Vegetation Index (NDVI) products over cultivated fields, as a function of time, led to the classification of irrigated and rainfed areas on the Kairouan plain, and to the production of irrigation maps at the scale of individual fields. This classification is based on a decision tree approach, using a combination of various statistical indices of soil moisture and NDVI time series. The resulting irrigation maps were validated using reference fields within the study site. The best results were obtained with classifications based on soil moisture indices only, with an accuracy of 77%.

This paper assesses the potential of Synthetic Aperture Radar (SAR) in the C and L bands to penetrate into the canopy cover of wheat, maize and grasslands. For wheat and grasslands, the sensitivity of the C and L bands to in situ surface soil moisture (SSM) was first studied according to three levels of the Normalized Difference Vegetation Index (NDVI < 0.4, 0.4 < NDVI < 0.7, and NDVI > 0.7). Next, the temporal evolution of the SAR signal in the C and L bands was analyzed according to SSM and the NDVI. For wheat and grasslands, the results showed that the L-band in HH polarization penetrates the canopy even when the canopy is well-developed (NDVI > 0.7), whereas the penetration of the C-band into the canopy is limited for an NDVI < 0.7. For an NDVI less than 0.7, the sensitivity of the radar signal to SSM is approximately 0.27 dB/vol.% for the L-band in HH polarization and approximately 0.12 dB/vol.% for the C-band (in both VV and VH polarizations). For highly developed wheat and grassland cover (NDVI > 0.7), the sensitivity of the L-band in HH polarization to SSM is approximately 0.19 dB/vol.%, whereas as the C-band is insensitive to SSM. For maize, only the temporal evolution of the C-band according to SSM and the NDVI was studied because the swath of SAR images in the L-band did not cover the maize plots. The results showed that the C-band in VV polarization is able to penetrate the maize canopy even when the canopy is well developed (NDVI > 0.7) due to high-order scattering along the soil-vegetation pathway that contains a soil contribution. According to results obtained in this paper, the L-band would penetrate a well-developed maize cover since the penetration depth of the L-band is greater than that of the C-band.

The objective of this paper is to evaluate the potential of Sentinel-1 Synthetic Aperture Radar “SAR” data (C-band) for monitoring agricultural frozen soils. First, investigations were conducted from simulated radar signal data using a SAR backscattering model combined with a dielectric mixing model. Then, Sentinel-1 images acquired at a study site near Paris, France were analyzed using temperature data to investigate the potential of the new Sentinel-1 SAR sensor for frozen soil mapping. The results show that the SAR backscattering coefficient decreases when the soil temperature drops below 0 °C. This decrease in signal is the most important for temperatures that ranges between 0 and −5 °C. A difference of at least 2 dB is observed between unfrozen soils and frozen soils. This difference increases under freezing condition when the temperature at the image acquisition date decreases. In addition, results show that the potential of the C-band radar signal for the discrimination of frozen soils slightly decreases when the soil moisture decreases (simulated data were used with soil moisture contents of 20 and 30 vol%). The difference between the backscattering coefficient of unfrozen soil and the backscattering coefficient of frozen soil decreases by approximately 1 dB when the soil moisture decreases from 30 to 20 vol%). Finally, the results show that both VV and VH allow a good detection of frozen soils but the sensitivity of VH is higher by approximately 1.5 dB. In conclusion, this study shows that the difference between a reference image acquired without freezing and an image acquired under freezing conditions is a good tool for detecting frozen soils.

This study evaluates the accuracy of several recent remote sensing Surface Soil Moisture (SSM) products at sites in southwestern France. The products used are Soil Moisture Active Passive “SMAP” (level 3: 36 km × 36 km, level 3 enhanced: 9 km × 9 km, and Level 2 SMAP/Sentinel-1: 1 km × 1km), Advanced Scatterometer “ASCAT” (level 2 with three spatial resolution 25 km × 25 km, 12.5 km × 12.5 km, and 1 km × 1 km), Soil Moisture and Ocean Salinity “SMOS” (SMOS INRA-CESBIO “SMOS-IC”, SMOS Near-Real-Time “SMOS-NRT”, SMOS Centre Aval de Traitement des Données SMOS level 3 “SMOS-CATDS”, 25 km × 25 km) and Sentinel-1(S1) (25 km × 25 km, 9 km × 9 km, and 1 km × 1 km). The accuracy of SSM products was computed using in situ measurements of SSM observed at a depth of 5 cm. In situ measurements were obtained from the SMOSMANIA ThetaProbe (Time Domaine reflectometry) network (7 stations between 1 January 2016 and 30 June 2017) and additional field campaigns (near Montpellier city in France, between 1 January 2017 and 31 May 2017) in southwestern France. For our study sites, results showed that (i) the accuracy of the Level 2 SMAP/Sentinel-1 was lower than that of SMAP-36 km and SMAP-9 km; (ii) the SMAP-36 km and SMAP-9 km products provide more precise SSM estimates than SMOS products (SMOS-IC, SMOS-NRT, and SMOS-CATDS), mainly due to higher sensitivity of SMOS to RFI (Radio Frequency Interference) noise; and (iii) the accuracy of SMAP-36 km and SMAP-9 km products was similar to that of ASCAT (ASCAT-25 km, ASCAT-12.5 km and ASCAT-1 km) and S1 (S1-25 km, S1-9 km, and S1-1 km) products. The accuracy of SMAP, Sentinel-1 and ASCAT SSM products calculated using the average of statistics obtained on each site is defined by a bias of about −3.2 vol. %, RMSD (Root Mean Square Difference) about 7.6 vol. %, ubRMSD (unbiased Root Mean Square Difference)about 5.6 vol. %, and R coefficient about 0.57. For SMOS products, the station average bias, RMSD, ubRMSD, and R coefficient were about −10.6 vol. %, 12.7 vol. %, 5.9 vol. %, and 0.49, respectively.

Soil moisture mapping at a high spatial resolution is very important for several applications in hydrology, agriculture and risk assessment. With the arrival of the free Sentinel data at high spatial and temporal resolutions, the development of soil moisture products that can better meet the needs of users is now possible. In this context, the main objective of the present paper is to develop an operational approach for soil moisture mapping in agricultural areas at a high spatial resolution over bare soils, as well as soils with vegetation cover. The developed approach is based on the synergic use of radar and optical data. A neural network technique was used to develop an operational method for soil moisture estimates. Three inversion SAR (Synthetic Aperture Radar) configurations were tested: (1) VV polarization; (2) VH polarization; and (3) both VV and VH polarization, all in addition to the NDVI information extracted from optical images. Neural networks were developed and validated using synthetic and real databases. The results showed that the use of a priori information on the soil moisture condition increases the precision of the soil moisture estimates. The results showed that VV alone provides better accuracy on the soil moisture estimates than VH alone. In addition, the use of both VV and VH provides similar results, compared to VV alone. In conclusion, the soil moisture could be estimated in agricultural areas with an accuracy of approximately 5 vol % (volumetric unit expressed in percent). Better results were obtained for soil with a moderate surface roughness (for root mean surface height between 1 and 3 cm). The developed approach could be applied for agricultural plots with an NDVI lower than 0.75.

In a perspective to develop an inversion approach for estimating surface soil moisture of crop fields from Sentinel-1/2 data (radar and optical sensors), the Water Cloud Model (WCM) was calibrated from C-band Synthetic Aperture Radar (SAR) data and Normalized Difference Vegetation Index (NDVI) values collected over crops fields and grasslands. The soil contribution that depends on soil moisture and surface roughness (in addition to SAR instrumental parameters) was simulated using the physical backscattering model IEM (Integral Equation Model). The vegetation descriptor used in the WCM is the NDVI because it can be directly calculated from optical images. A large dataset consisting of radar backscattered signal in Vertical transmit and Vertical receive (VV) and Vertical transmit and Horizontal receive (VH) polarizations with wide range of incidence angle, soil moisture, surface roughness, and NDVI-values was used. It was collected over two agricultural study sites. Results show that the soil contribution to the total radar backscattered signal is lower in VH than in VV because VH is more sensitive to vegetation cover. Thus, the use of VH alone or in addition to VV for retrieving the soil moisture is not advantageous in presence of well-developed vegetation cover.

Mapping forest AGB (Above Ground Biomass) is of crucial importance to estimate the carbon emissions associated with tropical deforestation. This study proposes a method to overcome the saturation at high AGB values of existing AGB map (Vieilledent’s AGB map) by using a map of correction factors generated from GLAS (Geoscience Laser Altimeter System) spaceborne LiDAR data. The Vieilledent’s AGB map of Madagascar was established using optical images, with parameters calculated from the SRTM Digital Elevation Model, climatic variables, and field inventories. In the present study, first, GLAS LiDAR data were used to obtain a spatially distributed (GLAS footprints geolocation) estimation of AGB (GLAS AGB) covering Madagascar forested areas, with a density of 0.52 footprint/km2. Second, the difference between the AGB from the Vieilledent’s AGB map and GLAS AGB at each GLAS footprint location was calculated, and additional spatially distributed correction factors were obtained. Third, an ordinary kriging interpolation was thus performed by taking into account the spatial structure of these additional correction factors to provide a continuous correction factor map. Finally, the existing and the correction factor maps were summed to improve the Vieilledent’s AGB map. The results showed that the integration of GLAS data improves the precision of Vieilledent’s AGB map by approximately 7 t/ha. By integrating GLAS data, the RMSE on AGB estimates decreases from 81 t/ha (R2 = 0.62) to 74.1 t/ha (R2 = 0.71). Most importantly, we showed that this approach using LiDAR data avoids underestimating high biomass values (new maximum AGB of 650 t/ha compared to 550 t/ha with the first approach).

The objective of this paper is to propose a new semi-empirical radar backscattering model for bare soil surfaces based on the Dubois model. A wide dataset of backscattering coefficients extracted from synthetic aperture radar (SAR) images and in situ soil surface parameter measurements (moisture content and roughness) is used. The retrieval of soil parameters from SAR images remains challenging because the available backscattering models have limited performances. Existing models, physical, semi-empirical, or empirical, do not allow for a reliable estimate of soil surface geophysical parameters for all surface conditions. The proposed model, developed in HH, HV, and VV polarizations, uses a formulation of radar signals based on physical principles that are validated in numerous studies. Never before has a backscattering model been built and validated on such an important dataset as the one proposed in this study. It contains a wide range of incidence angles (18°–57°) and radar wavelengths (L, C, X), well distributed, geographically, for regions with different climate conditions (humid, semi-arid, and arid sites), and involving many SAR sensors. The results show that the new model shows a very good performance for different radar wavelengths (L, C, X), incidence angles, and polarizations (RMSE of about 2 dB). This model is easy to invert and could provide a way to improve the retrieval of soil parameters.

Land monitoring using temporal series of Synthetic Aperture Radar (SAR) images requires radiometrically well calibrated sensors. In this paper, the radiometric stability of the new SAR Sentinel-1A “S-1A” sensor was first assessed by analyzing temporal variations of the backscattering coefficient (σ°) returned from invariant targets. Second, the radiometric level of invariant targets was compared from S-1A and Radarsat-2 “RS-2” data. The results show three stable sub-time series of S-1A data. The first (between 1 October 2014 and 19 March 2015) and third (between 25 November 2015 and 1 February 2016) sub-time series have almost the same mean σ°-values (a difference lower than 0.3 dB). The mean σ°-value of the second sub-time series (between 19 March 2015 and 25 November 2015) is higher than that of the first and the third sub-time series by roughly 0.9 dB. Moreover, our results show that the stability of each sub-time series is better than 0.48 dB. In addition, the results show that S-1A images of the first and third sub-time series appear to be well calibrated in comparison to RS-2 data, with a difference between S-1A and RS-2 lower than 0.3 dB. However, the S-1A images of the second sub-time series have σ°-values that are higher than those from RS-2 by roughly 1 dB.

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.

The objective of this paper is to propose a new semi-empirical radar backscattering model for bare soil surfaces based on the Dubois model. A wide dataset of backscattering coefficients extracted from synthetic aperture radar (SAR) images and in situ soil surface parameter measurements (moisture content and roughness) is used. The retrieval of soil parameters from SAR images remains challenging because the available backscattering models have limited performances. Existing models, physical, semi-empirical, or empirical, do not allow for a reliable estimate of soil surface geophysical parameters for all surface conditions. The proposed model, developed in HH, HV, and VV polarizations, uses a formulation of radar signals based on physical principles that are validated in numerous studies. Never before has a backscattering model been built and validated on such an important dataset as the one proposed in this study. It contains a wide range of incidence angles (18°-57°) and radar wavelengths (L, C, X), well distributed, geographically, for regions with different climate conditions (humid, semi-arid, and arid sites), and involving many SAR sensors. The results show that the new model shows a very good performance for different radar wavelengths (L, C, X), incidence angles, and polarizations (RMSE of about 2 dB). This model is easy to invert and could provide a way to improve the retrieval of soil parameters. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]