Monitoring the biophysical parameters of agricultural surfaces is a key issue for food security in the context of climate change. Since 2016, agricultural surfaces can be monitored from space at high spatial resolution (similar to 10/30 m) in the microwave and optical domains owing to radiometer and SAR sensors onboard Sentinel-1, -2 and Landsat-8 satellites. This paper draws on multi-temporal acquisitions over a six-year period to analyze satellite time series for the main winter and summer crops (corn, sunflower, soybean, sorghum, rapeseed, wheat) grown in southwestern France and more widely cultivated around the world. From January 2016 to December 2021, satellite signals extracted at the field spatial scale offer a unique opportunity to monitor agricultural surfaces with a high temporal resolution (every 1 or 2 days) never achieved before thanks to the combination of multi-sensor and multi-orbit data. Analyses on the impact of the topography and satellites' viewing angles showed that the NDVI values derived from Sentinel-2 and Landsat-8 are very close (r > 0.92) and can be merged to construct multi-annual time series. Angular sensitivity is much more pronounced for radar images; while it demonstrates a weaker cross-polarization and polarization ratio, it is greater for co-polarization. Optical and radar time series are modulated in time and amplitude depending on yearly climatic conditions and agricultural practices. The combined use of the ascending and descending orbits of the two Sentinel-1 satellites makes it possible to detect specific periods (harvest, flowering) for certain crops (wheat and sunflower). The long-term approach has enabled the modeling of satellite time series using double logistic functions with good performance (r > 0.92 on average), allowing the identification of interannual variations of crop development driven by climatic conditions and agricultural practices.
Climate change greatly impacts forest carbon and water dynamics in the boreal biome. To understand the vulnerability of the North American boreal biome to climate change disturbances, including droughts and temperature extremes, a biome-wide characterization of near-surface soil moisture (SM) is essential. SM is an important control on photosynthetic carbon dioxide uptake and the associated evaporative loss of water through transpiration. Two missions with passive microwave radiometers operating at L-band (1-2 GHz) were launched to estimate SM at a global scale: soil moisture and ocean salinity (SMOS) launched in 2009 and soil moisture active passive (SMAP) launched in 2015. In this study, we analyzed the relationships between various microwave observations [brightness temperature (TB), soil reflectivity (r), and soil moisture (SM)] derived from SMOS and SMAP with in situ SM measurements (2016-2021) over 16 sites across the boreal biome of North America. Our results show good agreement between SMAP and SMOS in terms of relationships between TB and in situ SM measurements (with a coefficient of determination R (2) > 0.7). Analysis comparing microwave variables (TB, r, and SM) and in situ SM data demonstrated considerable differences across sites and years. We observed moderate to strong correlation coefficients (R_TB > 0.5) between the TB and in situ SM data, generally when a seasonal soil drying trend emerges over the summer. However, soils with high water content and frequent precipitation input showed weaker correlations between TB and in situ SM. By analyzing the correlations between r and in situ SM data, we demonstrated that the soil effective temperature partly affects performances in terms of R_TB. Comparisons with other SMAP SM retrievals (multitemporal dual-channel algorithm and multichannel collaborative algorithm) revealed similar results in correlations and showed high spatiotemporal differences across all SM retrievals and study sites. Our findings highlight the complexity of the interactions between microwave observations and SM dynamics across the boreal biome, thus challenging the L-band passive microwave SM retrievals.
This study analyses a long-time series (from January 2016 to December 2021) of optical and SAR signatures. The data were processed into two contexts: 1) over a station-monitored field called Aurade in southwestern France (part of the ICOS network) to examine the impact of surface states on satellite signals and, 2) over 939 fields of sunflower to investigate the impact of climatic conditions on crop development. Results show that backscatters can be used with moderate confidence to replace or gap fill NDVI (r(6-years)(Aurade) = 0.73 over the Aurade field and rr(6-years)(sunflower)=0.77 over the 939 fields of sunflower). Nevertheless, additional analyses highlighted the effectiveness of combining gamma(VH)-V-0, gamma(VH)-V-0/VV with NDVI, GAI (Green Area Index), and Fcover indices for complementary monitoring of crop phenology. Superpositions of satellite signals observed on sunflowers suggest that it is possible to detect development anomalies using optical or radar signals based on yearly climatic conditions. Results also show the importance of considering both sunflower architecture (orientation, inter-row spacing, F-cover) and radar acquisition geometry (orbit pass, viewing angles), in contrast to denser or more covering crops. The joint use of the 2 SAR orbits enabled the detection of phenological stages (from flowering to harvest) not as well detectable with a single orbit. This long-time series analysis provided insights under multiple climatic conditions, including the hottest year in France since the beginning of meteorological records.
The present study aims to investigate the potential of multi-configuration Sentinel-1 (S-1) synthetic aperture radar (SAR) images for characterizing four wheat parameters: total fresh mass (TFM), total dry mass (TDM), plant heights (He), and water content (WC). Because they are almost independent on the weather conditions, we have chosen to use only SAR. Samples of wheat parameters were collected over seven fields (three irrigated and four rainfed fields) in Southwestern France. We first analyzed the temporal behaviors of wheat parameters (He, TDM, TFM and WC) between February and June 2016. Then, the temporal profiles of the S-1 backscattering coefficients (VV, VH), the difference (VH − VV), the sum of the polarizations (VH + VV) and their cumulative values are analyzed for two orbits (30 and 132) during the wheat-growing season (from January to July 2016). After that, S-1 signals were statistically compared with all crop parameters considering the impact of pass orbit, irrigation and two vegetative periods in order to identify the best S-1 configuration for estimating crop parameters. Interesting S-1 backscattering behaviors were observed with the various wheat parameters after separating irrigation impacts and vegetative periods. For the orbit 30 (mean incidence angle of 33.6°); results show that the best S-1 configurations (with coefficient of determination (R2) > 0.7) were obtained using the VV and VH + VV as a function of the He, TDM and WC, over irrigated fields and during the second vegetative period. For the orbit 132 (mean incidence angle of 43.4°), the highest dynamic sensitivities (R2 > 0.8) were observed for the VV and VH + VV configurations with He, TDM and TFM over irrigated fields during the first vegetative period. Overall, the sensitivity of S-1 data to wheat variables depended on the radar configuration (orbits and polarizations), the vegetative periods and was often better over irrigated fields in comparison with rainfed ones. Significant improvements of the determination coefficients were obtained when the cumulative (VH + VV) index was considered for He (R² > 0.9), TDM (R² > 0.9) and TFM (R² > 0.75) for irrigated fields, all along the crop cycle. The estimate of WC was more limited (R² > 0.6) and remained limited to the second period of the vegetation cycle (from flowering onwards). Whatever parameters were considered, the relative errors never exceeded 23%. This study has shown the importance of considering the agricultural practices (irrigation) and vegetative periods to effectively monitor some wheat parameters with S-1 data.
Soil texture is a key parameter in agricultural processes and an important measure for agricultural prediction, water cycle, filtering of pollutants and carbon storage. Besides, its estimation is essential for agronomists, hydrologists, geologists and environmentalists and for modeling in these application areas. Several studies have been based on understanding and modeling the biological, physical and chemical processes in the soil. Regarding the texture of the soil, few researches propose soil texture spatialization, and are generally based on ground measurements. Among other things, field observations or laboratory analyzes are very expensive and are not very representative. Indeed, the soil texture presents a strong heterogeneity even at the scale of a field. It is then necessary to use precise and spatialized information on soils. These methods are generally based on remote sensing data and particularly optical data to restore soil component. However, these techniques are strongly affected by atmospheric conditions. This constraint is not valid for Radar sensors (Radio Detection And Ranging). Radar data are mainly sensitive to soil moisture and soil roughness, and has also been evaluated for its ability to perform texture measurements. The aim of this study is evaluate the potential of these techniques based on optical and radar data for soil texture estimation. By its composition, its structure, its texture and its porosity, soil moisture is strongly influenced by the soil nature. With the arrival of Sentinel-1 (S-1) and Sentinel-2 (S-2) ESA spatial missions, data are acquired with high spatial and temporal resolution between July and early December 2017, on a semi-arid area in central Tunisia. This study is therefore conducted using S-2 SWIR (Short-Wave Infrared) bands (B11 and B12, most sensitive to clay) and soil moisture products derived from radar data. And algorithms based on the support vector machine (SVM) and random forest (RF) methods are proposed for the classification and mapping of clay content. In order to evaluate the approach and determine the adequate data (between optical and radar data) allowing to precisely characterize the clay content, a cross-validation was used. The SWIR bands lead to less satisfactory outcomes compared to soil moisture. With an overall accuracy of approximately 65%, soil moisture achieved the best performance for estimating soil texture. The results also showed that RF and SVM are robust classifiers for texture estimation despite the small number of training data. However, RF displays greater accuracy and speed of simulation compared to SVM.
Soil texture estimation is important in several applications. However, field sampling or laboratory analyzes are very expensive and not very representative. Besides, existing soil maps are neither exhaustive nor sufficiently precise for modeling and meeting needs at field scale. This study aimed to explore the potential of Sentinel satellites to predict topsoil texture, and more precisely to produce clay content map at fine spatial resolution. With its components and its porosity, soil texture is directly linked to soil moisture. In this context and with the arrival of Sentinel constellation, data are acquired with high spatial and temporal resolution. And soil moisture is retrieved from a synergetic use of Sentinel-l (S-1) and Sentinel-2 (S-2) data between July and early December 2017, over a semi-arid area in central Tunisia. Relationship between soil moisture and clay content is studied and used to produce texture map. Classification algorithm based on random forest (RF) is used for the mapping of clay content classes. The results showed the potential of S-1 and S-2 products to predict soil texture.
This paper discusses the combined use of remotely sensed optical and radar data for the estimation and mapping of soil texture. The study is based on Sentinel-1 (S-1) and Sentinel-2 (S-2) data acquired between July and early December 2017, on a semi-arid area about 3000 km2 in central Tunisia. In addition to satellite acquisitions, texture measurement samples were taken in several agricultural fields, characterized by a large range of clay contents (between 13% and 60%). For the period between July and August, various optical indicators of clay content Short-Wave Infrared (SWIR) bands and soil indices) were tested over bare soils. Satellite moisture products, derived from combined S-1 and S-2 data, were also tested as an indicator of soil texture. Algorithms based on the support vector machine (SVM) and random forest (RF) methods are proposed for the classification and mapping of clay content and a three-fold cross-validation is used to evaluate both approaches. The classifications with the best performance are achieved using the soil moisture indicator derived from combined S-1 and S-2 data, with overall accuracy (OA) of 63% and 65% for the SVM and RF classifications, respectively.
The aim of this paper was to demonstrate the potential of multi-temporal X-band Synthetic Aperture Radar (SAR) moisture products to be used for the calibration of a model reproducing soil moisture (SM) variations. We propose the MHYSAN model (Modele de bilan HYdrique des Sols Agricoles Nus) for simulating soil water balance of bare soils. This model was used to simulate surface evaporation fluxes and soil moisture content (SMC) at daily time scale over a semi-arid, bare agricultural site in Tunisia (North Africa). The MHYSAN model was calibrated using seven very high-resolution SAR (TerraSAR-X) SM outputs ranging over only two months. Then, the simulated SM was validated using continuous thetaprobe measurements during 15 months. The high performances observed could be explained by the fact that although images were acquired during a short time range, the important soil moisture variation captured allowed a good calibration of the soil parameters. These results highlight the potential of Sentinel-1 images for daily soil moisture monitoring using simple models.
The aim of this paper is to estimate geometric, water and physical surface soil parameters from typical semi-arid regions made over bare study area (North Africa) using multi-temporal X-band SAR images (TerraSAR-X). For spatial and temporal surface roughness estimation, empirical relationships between radar and soil roughness parameters (rms height "Hrms", and Zg parameter) were proposed. Two roughness classes are identified through radar signal inversion (smooth and ploughed soils). For the retrieval of surface soil moisture at a high spatial resolution, an algorithm combing TerraSAR-X images with continuous thetaprobe measurements was proposed. Two assumptions were studied: (1) roughness variations during the radar acquisition campaigns were not accounted for; (2) a simple correction for temporal variations in roughness was included. Finally, an empirical relationship was established between the mean moisture values retrieved from the SAR images and the percentage of clay over several test fields. Results showed that highly accurate clay estimations can be achieved.
The present study highlights the potential of multi-temporal X-band Synthetic Aperture Radar (SAR) moisture products to be used for the calibration of a model reproducing soil moisture (SM) variations. We propose the MHYSAN model (Modèle de bilan HYdrique des Sols Agricoles Nus) for simulating soil water balance of bare soils. This model was used to simulate surface evaporation fluxes and SM content at daily time scale over a semi-arid, bare agricultural site in Tunisia (North Africa). Two main approaches are considered in this study. Firstly, the MHYSAN model was successfully calibrated for seven sites using continuous thetaprobe measurements at two depths. Then the possibility to extrapolate local SM simulations at distant sites, based on soil texture similarity only, was tested. This extrapolation was assessed using SAR estimates and manual thetaprobe measurements of SM recorded at these distant sites. The results reveal a bias of approximately 0.63% and 3.04%, and an RMSE equal to 6.11% and 4.5%, for the SAR volumetric SM and manual thetaprobe measurements, respectively. In a second approach, the MHYSAN model was calibrated using seven very high-resolution SAR (TerraSAR-X) SM outputs ranging over only two months. The simulated SM were validated using continuous thetaprobe measurements during 15 months. Although the SM was measured on only seven different dates for the purposes of calibration, satisfactory results were obtained as a result of the wide range of SM values recorded in these seven images. This led to good overall calibration of the soil parameters, thus demonstrating the considerable potential of Sentinel-1 images for daily soil moisture monitoring using simple models.
The aim of the present study is to analyze the sensitivity of X-band SAR (TerraSAR-X) signals as a function of different physical bare soil parameters (soil moisture, soil roughness), and to evaluate the accuracy of change detection approach proposed for soil moisture estimation. Firstly, we presented a brief description of our ground and satellite database. Secondly, we considered the main results of our statistical analysis of the relationships between radar and soil parameters: soil moisture and different roughness parameters (the rms height, Zs parameter, and a new roughness parameter Zg. Finally, we proposed an algorithm combing multi-temporal X-band SAR images (TerraSAR-X) with different continuous thetaprobe measurements for the retrieval of surface soil moisture at a high spatial resolution.
The goal of this paper is to analyze the sensitivity of X-band SAR (TerraSAR-X) signals as a function of different physical bare soil parameters (soil moisture, soil roughness), and to demonstrate that it is possible to estimate of both soil moisture and texture from the same experimental campaign, using a single radar signal configuration (one incidence angle, one polarization). Firstly, we analyzed statistically the relationships between X-band SAR (TerraSAR-X) backscattering signals function of soil moisture and different roughness parameters (the root mean square height Hrms, the Zs parameter and the Zg parameter) at HH polarization and for an incidence angle about 36°, over a semi-arid site in Tunisia (North Africa). Results have shown a high sensitivity of real radar data to the two soil parameters: roughness and moisture. A linear relationship is obtained between volumetric soil moisture and radar signal. A logarithmic correlation is observed between backscattering coefficient and all roughness parameters. The highest dynamic sensitivity is obtained with Zg parameter. Then, we proposed to retrieve of both soil moisture and texture using these multi-temporal X-band SAR images. Our approach is based on the change detection method and combines the seven radar images with different continuous thetaprobe measurements. To estimate soil moisture from X-band SAR data, we analyzed statistically the sensitivity between radar measurements and ground soil moisture derived from permanent thetaprobe stations. Our approaches are applied over bare soil class identified from an optical image SPOT / HRV acquired in the same period of measurements. Results have shown linear relationship for the radar signals as a function of volumetric soil moisture with high sensitivity about 0.21 dB/vol%. For estimation of change in soil moisture, we considered two options: (1) roughness variations during the three-month radar acquisition campaigns were not accounted for; (2) a simple correction for temporal variations in roughness was included. The results reveal a small improvement in the estimation of soil moisture when a correction for temporal variations in roughness is introduced. Finally, by considering the estimated temporal dynamics of soil moisture, a methodology is proposed for the retrieval of clay and sand content (expressed as percentages) in soil. Two empirical relationships were established between the mean moisture values retrieved from the seven acquired radar images and the two soil texture components over 36 test fields. Validation of the proposed approach was carried out over a second set of 34 fields, showing that highly accurate clay estimations can be achieved.
The aim of this paper is to propose a methodology combing multi-temporal X-band SAR images (TerraSAR-X) with continuous ground thetaprobe measurements, for the retrieval of surface soil moisture and texture at a high spatial resolution. Our analysis is based on seven radar images acquired at a 36° incidence angle in the HH polarization, over a semi-arid site in Tunisia (North Africa). The soil moisture estimations are based on an empirical change detection approach using TerraSAR-X data and ground auxiliary thetaprobe network measurements. Two assumptions were tested: (1) roughness variations during the three-month radar acquisition campaigns were not accounted for; (2) a simple correction for temporal variations in roughness was included. The results reveal a small improvement in the estimation of soil moisture when a correction for temporal variations in roughness is introduced. By considering the estimated temporal dynamics of soil moisture, a methodology is proposed for the retrieval of clay and sand content (expressed as percentages) in soil. Two empirical relationships were established between the mean moisture values retrieved from the seven acquired radar images and the two soil texture components over 36 test fields. Validation of the proposed approach was carried out over a second set of 34 fields, showing that highly accurate clay estimations can be achieved. Maps of soil moisture, clay and sand percentages at the studied site are derived.
The aim of this paper is to analyze the potential of X-band SAR measurements (COSMO-SkyMed and TerraSAR-X) made over bare soils for the estimation of soil moisture and surface geometry parameters at a semi-arid site in Tunisia (North Africa). Radar signals acquired with different configurations (HH and VV polarizations, incidence angles of 26° and 36°) are statistically compared with ground measurements (soil moisture and roughness parameters). The radar measurements are found to be highly sensitive to the various soil parameters of interest. A linear relationship is determined for the radar signals as a function of volumetric soil moisture, and a logarithmic correlation is observed between the radar signals and three surface roughness parameters: the root mean square height (Hrms), the parameter Zs = Hrms2/l (where l is the correlation length) and the parameter Zg = Hrms × (Hrms/l)α (where α is the power of the surface height correlation function). The highest dynamic sensitivity is observed for Zg at high incidence angles. Finally, the performance of different physical and semi-empirical backscattering models (IEM, Baghdadi-calibrated IEM and Dubois models) is compared with SAR measurements. The results provide an indication of the limits of validity of the IEM and Dubois models, for various radar configurations and roughness conditions. Considerable improvements in the IEM model performance are observed using the Baghdadi-calibrated version of this model.
The goal of this paper is to analyze the potential of COSMO-SkyMed and TerraSAR-X SAR measurements over bare soils in order to estimate correctly soil parameters. We analyzed statistically the relationships between X-SAR backscattering signals function of soil moisture and different roughness parameters (the root mean square height Hrms, the Zs parameter and the Zg parameter) at HH polarization and for an incidence angle about 35.5°. Results have shown a high sensitivity of real radar data to the two soil parameters: roughness and moisture. A linear relationship is obtained between volumetric soil moisture and radar signal with the strongest correlation observed with gravimetric moisture measurements. A logarithmic correlation is observed between backscattering coefficient and all roughness parameters. The highest dynamic sensitivity is obtained with Zg parameter.
The goal of this study is to discuss the effect of multi-frequency radar configurations on the relationship between surface soil moisture and the nature of radar backscatter over bare soils. In an attempt to answer this question, the Advanced Integral Equation Model (AIEM) was used to simulate backscatter from soil surfaces with various moisture vertical profiles, for three frequency bands: L, C and X. In these computations, we investigated the influence of the vertical heterogeneity of soil moisture on the characteristics of the backscattered signals. The effect of radar frequency is distinctly demonstrated. A database produced from Envisat ASAR and TerraSAR-X data, acquired over bare soils with in situ measurements of moisture content and ground surface roughness, was used to validate the usefulness of taking the soil moisture heterogeneity into account in the backscattering model.