Accurate information about the irrigated surface is essential to help assess the impact of irrigation on water consumption, the hydrological cycle and regional climate. In this study, we compare ...recently developed operational and spatially transferrable classification models proposed for irrigation mapping. The first model suggests the use of spatio-temporal soil moisture indices derived from the Sentinel-1/2 soil moisture product (S2MP) at plot scale to map irrigated areas using the unsupervised K-means clustering algorithm (Dari model). The second model called the Sentinel-1/2 Irrigation mapping (S2IM) is a classification model based on the use the Sentinel-1 (S1) and Sentinel-2 (S2) time series data. Five study cases were examined including four studied years in a semi-oceanic area in north-central France (between 2017 and 2020) and one year (2020) in a Mediterranean context in south France. Main results showed that the soil-moisture based model using K-means clustering (Dari model) performs well for irrigation mapping but remains less accurate than the S2IM model. The overall accuracy of the Dari model ranged between 72.1% and 78.4% across the five study cases. The Dari model was found to be limited over humid conditions as it fails to correctly distinguish rain-fed plots from irrigated plots with an accuracy of the rain-fed class reaching 24.2% only. The S2IM showed the best accuracy in the five study cases with an overall accuracy ranging between 72.8% and 93.0%. However, for humid climatic conditions, the S2IM had an accuracy of the rain-fed class reaching 62.0%. The S2IM is thus superior in terms of accuracy but with higher complexity for application than the Dari model that remains simple yet effective for irrigation mapping.
The Global Ecosystem Dynamics Investigation (GEDI) Light Detection And Ranging (LiDAR) altimetry mission was recently launched to the International Space Station with a capability of providing ...billions of high-quality measurements of vertical structures globally. This study assesses the accuracy of the GEDI LiDAR altimetry estimation of lake water levels. The difference between GEDI’s elevation estimates to in-situ hydrological gauge water levels was determined for eight natural lakes in Switzerland. The elevation accuracy of GEDI was assessed as a function of each lake, acquisition date, and the laser used for acquisition (beam). The GEDI elevation estimates exhibit an overall good agreement with in-situ water levels with a mean elevation bias of 0.61 cm and a standard deviation (std) of 22.3 cm and could be lowered to 8.5 cm when accounting for instrumental and environmental factors. Over the eight studied lakes, the bias between GEDI elevations and in-situ data ranged from −13.8 cm to +9.8 cm with a standard deviation of the mean difference ranging from 14.5 to 31.6 cm. Results also show that the acquisition date affects the precision of the GEDI elevation estimates. GEDI data acquired in the mornings or late at night had lower bias in comparison to acquisitions during daytime or over weekends. Even though GEDI is equipped with three identical laser units, a systematic bias was found based on the laser units used in the acquisitions. Considering the eight studied lakes, the beams with the highest elevation differences compared to in-situ data were beams 1 and 6 (standard deviations of −10.2 and +18.1 cm, respectively). In contrast, the beams with the smallest mean elevation difference to in-situ data were beams 5 and 7 (−1.7 and −2.5 cm, respectively). The remaining beams (2, 3, 4, and 8) showed a mean difference between −7.4 and +4.4 cm. The standard deviation of the mean difference, however, was similar across all beams and ranged from 17.2 and 22.9 cm. This study highlights the importance of GEDI data for estimating water levels in lakes with good accuracy and has potentials in advancing our understanding of the hydrological significance of lakes especially in data scarce regions of the world.
The difficulty of calculating the daily water budget of irrigated fields is often due to the uncertainty surrounding irrigation amounts and timing. The automated detection of irrigation events has ...the potential to greatly simplify this process, and the combination of high-resolution SAR (Sentinel-1) and optical satellite observations (Sentinel-2) makes the detection of irrigation events feasible through the use of a surface soil moisture (SSM) product. The motivation behind this study is to utilize a large irrigation dataset (collected during the ESA Irrigation + project over five sites in three countries over three years) to analyze the performance of an established algorithm and to test potential improvements. The study’s main findings are (1) the scores decrease with SSM observation frequency; (2) scores decrease as irrigation frequency increases, which was supported by better scores in France (more sprinkler irrigation) than in Germany (more localized irrigation); (3) replacing the original SSM model with the force-restore model resulted in an improvement of about 6% in the F-score and narrowed the error on cumulative seasonal irrigation; (4) the Sentinel-1 configuration (incidence angle, trajectory) did not show a significant impact on the retrieval of irrigation, which supposes that the SSM is not affected by these changes. Other aspects did not allow a definitive conclusion on the irrigation retrieval algorithm: (1) the lower scores obtained with small NDVI compared to large NDVI were counter-intuitive but may have been due to the larger number of irrigation events during high vegetation periods; (2) merging different runs and interpolating all SSM data for one run produced comparable F-scores, but the estimated cumulative sum of irrigation was around −20% lower compared to the reference dataset in the best cases.
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.
In semi-arid areas characterized by frequent drought events, there is often a strong need for an operational grain yield forecasting system, to help decision-makers with the planning of annual ...imports. However, monitoring the crop canopy and production capacity of plants, especially for cereals, can be challenging. In this paper, a new approach to yield estimation by combining data from the Simple Algorithm for Yield estimation (SAFY) agro-meteorological model with optical SPOT/ High Visible Resolution (HRV) satellite data is proposed. Grain yields are then statistically estimated as a function of Leaf Area Index (LAI) during the maximum growth period between 25 March and 5 April. The LAI is retrieved from the SAFY model, and calibrated using SPOT/HRV data. This study is based on the analysis of a rich database, which was acquired over a period of two years (2010⁻2011, 2012⁻2013) at the Merguellil site in central Tunisia (North Africa) from more than 60 test fields and 20 optical satellite SPOT/HRV images. The validation and calibration of this methodology is presented, on the basis of two subsets of observations derived from the experimental database. Finally, an inversion technique is applied to estimate the overall yield of the entire studied site.
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.
A new approach based on the change detection technique is proposed for the estimation of surface soil moisture (SSM) from a time series of radar measurements. A new index of reflectivity (IR) is ...defined that uses radar signals and Fresnel coefficients. This index is equal to 0 in the case of the smallest value of the Fresnel coefficient, corresponding to the driest conditions and the weakest radar signal, and is equal to 1 for the highest value of the Fresnel coefficient, corresponding to the wettest soil conditions and the strongest radar signal. The integrated equation model is used to simulate the behavior of radar signals as a function of soil moisture and roughness. This approach validates the greater usefulness of the IR compared with that of the commonly used index of SSM (ISSM), which assumes that the SSM varies linearly as a function of radar signal strength. The IR-based approach was tested using Sentinel-1 radar data recorded over three regions: Banizombou (Niger), Merguellil (Tunisia), and Occitania (France). The IR approach was found to perform better for the estimation of SSM than the ISSM approach based on comparisons with ground measurements over bare soils.
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 (Hdom) 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 Hdom 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 Hdom, 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 (HTn) and the height at different energy quartiles of the bare ground waveform (HGn) was assessed. Two sets of the HGn 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 Hdom 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 m3·ha−1 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 m3·ha−1 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 m3·ha−1 was observed for slopes between 20 and 45%.
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 <1.5 cm and most SSM values higher than 10 vol.%, the use of Hrms as an input in the NNs decreases the underestimation of the SSM (bias ranges from −4.5 to 0 vol.%) and provides a more accurate estimation of the SSM with a decrease in the RMSE by approximately 2 vol.%. Moreover, for Hrms values between 1.5 and 2.0 cm, the overestimation of SSM slightly decreases (bias decreased by around 1.0 vol.%) without a significant improvement of the RMSE. In addition, 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.
This article aims to estimate the water footprint (WF) of cereals—specifically, wheat and barley—in the Kairouan plain, located in central Tunisia. To achieve this objective, two components must be ...determined: actual evapotranspiration (ETa) and crop yield. The study covers three growing seasons from 2010 to 2013. The ETa estimation employed the S-SEBI (simplified surface energy balance index) model, utilizing Landsat 7 and 8 optical and thermal infrared spectral bands. For yield estimation, an empirical model based on the normalized difference vegetation index (NDVI) was applied. Results indicate the effectiveness of the S-SEBI model in estimating ETa, demonstrating an R2 of 0.82 and an RMSE of 0.45 mm/day. Concurrently, yields mapped over the area range between 6 and 77 qx/ha. Globally, cereals’ average WF varied from 1.08 m3/kg to 1.22 m3/kg over the three study years, with the majority below 1 m3/kg. Notably in dry years, the importance of the blue WF is emphasized compared to years with average rainfall (WFb-2013 = 1.04 m3/kg, WFb-2012 = 0.61 m3/kg, WFb-2011 = 0.41 m3/kg). Moreover, based on an in-depth agronomic analysis combining yields and WF, four classes were defined, ranging from the most water efficient to the least, revealing that over 30% of cultivated areas during the study years (approximately 40% in 2011 and 2012 and 29% in 2013) exhibited low water efficiency, characterized by low yields and high WF. A unique index, the WFI, is proposed to assess the spatial variability of green and blue water. Spatial analysis using the WFI highlighted that in 2012, 40% of cereal plots with low yields but high water consumption were irrigated (81% blue water compared to 6% in 2011).