Empirical fitting of forward backscattering models for multitemporal retrieval of soil moisture from radar data at L-band

A multitemporal algorithm, originally conceived for the C-band radar aboard the Sentinel-1 satellite, has been updated to retrieve soil moisture from L-band radar data, such as those provided by the National Aeronautics and Space Administration Soil Moisture Active/Passive (SMAP) mission. This type of algorithm may deliver more accurate soil moisture maps that mitigate the effect of roughness and vegetation changes. Within the multitemporal inversion scheme based on the Bayesian maximum a posteriori probability (MAP) criterion, a dense time series of radar measurements is integrated to invert a forward backscattering model. The model calibration and validation tasks have been accomplished using the data collected during the SMAP validation experiment 12 spanning several soil conditions (pasture, wheat, corn, and soybean). The data have been used to update the forward model for bare soil scattering at L-band and to tune a simple vegetation scattering model considering two different classes of vegetation: those producing mainly single scattering effects and those characterized by a significant multiple scattering involving terrain surface and vegetation elements interaction. The algorithm retrievals showed a root mean square difference (RMSD) around 5% over bare soil, soybean, and cornfields. As for wheat, a bias was observed; when removed, the RMSD went down from 7.7% to 5%

Quantifying Irrigated Winter Wheat LAI in Argentina Using Multiple Sentinel-1 Incidence Angles

Synthetic aperture radar (SAR) data provides an appealing opportunity for all-weather day or night Earth surface monitoring. The European constellation Sentinel-1 (S1) consisting of S1-A and S1-B satellites offers a suitable revisit time and spatial resolution for the observation of croplands from space. The C-band radar backscatter is sensitive to vegetation structure changes and phenology as well as soil moisture and roughness. It also varies depending on the local incidence angle (LIA) of the SAR acquisition’s geometry. The LIA backscatter dependency could therefore be exploited to improve the retrieval of the crop biophysical variables. The availability of S1 radar time-series data at distinct observation angles holds the feasibility to retrieve leaf area index (LAI) evolution considering spatiotemporal coverage of intensively cultivated areas. Accordingly, this research presents a workflow merging multi-date S1 smoothed data acquired at distinct LIA with a Gaussian processes regression (GPR) and a cross-validation (CV) strategy to estimate cropland LAI of irrigated winter wheat. The GPR-S1-LAI model was tested against in situ data of the 2020 winter wheat campaign in the irrigated valley of Colorador river, South of Buenos Aires Province, Argentina. We achieved adequate validation results for LAI with R2CV = 0.67 and RMSECV = 0.88 m2 m−2. The trained model was further applied to a series of S1 stacked images, generating temporal LAI maps that well reflect the crop growth cycle. The robustness of the retrieval workflow is supported by the associated uncertainties along with the obtained maps. We found that processing S1 smoothed imagery with distinct acquisition geometries permits accurate radar-based LAI modeling throughout large irrigated areas and in consequence can support agricultural management practices in cloud-prone agri-environments.EEA Hilario AscasubiFil: Caballero, Gabriel. Technological University of Uruguay (UTEC). Agri-Environmental Engineering; UruguayFil: Caballero, Gabriel. University of Valencia. Image Processing Laboratory (IPL); EspañaFil: Pezzola, Alejandro. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Winschel, Cristina Ines. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Casella, Alejandra. Instituto Nacional de Tecnología Agropecuaria (INTA). Instituto de Clima y Agua; ArgentinaFil: Sanchez Angonova, Paolo Andres. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Orden, Luciano. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Orden, Luciano. Universidad Miguel Hernández. Centro de Investigación e Innovación Agroalimentaria y Agroambiental. GIAAMA Reseach Group; EspañaFil: Berger, Katja. University of Valencia. Image Processing Laboratory (IPL); EspañaFil: Berger, Katja. Mantle Labs GmbH; AustriaFil: Verrelst, Jochem. University of Valencia. Image Processing Laboratory (IPL); EspañaFil: Delegido, Jesús. Universidad de Valencia. Image Processing Laboratory (IPL); Españ

Parameterizing the modified water cloud model to improve soil moisture data retrieval using vegetation models

L

Crop biophysical parameter retrieval from Sentinel-1 SAR data with a multi-target inversion of Water Cloud Model

Estimation of bio-and geophysical parameters from Earth observation (EO) data is essential for developing applications on crop growth monitoring. High spatio-temporal resolution and wide spatial coverage provided by EO satellite data are key inputs for operational crop monitoring. In Synthetic Aperture Radar (SAR) applications, a semi-empirical model (viz., Water Cloud Model (WCM)) is often used to estimate vegetation descriptors individually. However, a simultaneous estimation of these vegetation descriptors would be logical given their inherent correlation, which is seldom preserved in the estimation of individual descriptors by separate inversion models. This functional relationship between biophysical parameters is essential for crop yield models, given that their variations often follow different distribution throughout crop development stages. However, estimating individual parameters with independent inversion models presume a simple relationship (potentially linear) between the biophysical parameters. Alternatively, a multi-target inversion approach would be more effective for this aspect of model inversion compared to an individual estimation approach. In the present research, the multi-output support vector regression (MSVR) technique is used for inversion of the WCM from C-band dual-pol Sentinel-1 SAR data. Plant Area Index (PAI, m2 m−2) and wet biomass (W, kg m−2) are used as the vegetation descriptors in the WCM. The performance of the inversion approach is evaluated with in-situ measurements collected over the test site in Manitoba (Canada), which is a super-site in the Joint Experiment for Crop Assessment and Monitoring (JECAM) SAR inter-comparison experiment network. The validation results indicate a good correlation with acceptable error estimates (normalized root mean square error–nRMSE and mean absolute error–MAE) for both PAI and wet biomass for the MSVR approach and a better estimation with MSVR than single-target models (support vector regression–SVR). Furthermore, the correlation between PAI and wet biomass is assessed using the MSVR and SVR model. Contrary to the single output SVR, the correlation between biophysical parameters is adequately taken into account in MSVR based simultaneous inversion technique. Finally, the spatio-temporal maps for PAI and W at different growth stages indicate their variability with crop development over the test site.This research was supported in part by Shastri Indo-Candian Institute, New Delhi, India and the Spanish Ministry of Economy, Industry and Competitiveness, in part by the State Agency of Research (AEI), in part by the European Funds for Regional Development under project TEC2017-85244-C2-1-P

Microwave Indices from Active and Passive Sensors for Remote Sensing Applications

Past research has comprehensively assessed the capabilities of satellite sensors operating at microwave frequencies, both active (SAR, scatterometers) and passive (radiometers), for the remote sensing of Earth’s surface. Besides brightness temperature and backscattering coefficient, microwave indices, defined as a combination of data collected at different frequencies and polarizations, revealed a good sensitivity to hydrological cycle parameters such as surface soil moisture, vegetation water content, and snow depth and its water equivalent. The differences between microwave backscattering and emission at more frequencies and polarizations have been well established in relation to these parameters, enabling operational retrieval algorithms based on microwave indices to be developed. This Special Issue aims at providing an overview of microwave signal capabilities in estimating the main land parameters of the hydrological cycle, e.g., soil moisture, vegetation water content, and snow water equivalent, on both local and global scales, with a particular focus on the applications of microwave indices

Monitoring crops water needs at high spatio-temporal resolution by synergy of optical/thermal and radar observations

L'optimisation de la gestion de l'eau en agriculture est essentielle dans les zones semi-arides afin de préserver les ressources en eau qui sont déjà faibles et erratiques dues à des actions humaines et au changement climatique. Cette thèse vise à utiliser la synergie des observations de télédétection multispectrales (données radar, optiques et thermiques) pour un suivi à haute résolution spatio-temporelle des besoins en eau des cultures. Dans ce contexte, différentes approches utilisant divers capteurs (Landsat-7/8, Sentinel-1 et MODIS) ont été developpées pour apporter une information sur l'humidité du sol (SM) et le stress hydrique des cultures à une échelle spatio-temporelle pertinente pour la gestion de l'irrigation. Ce travail va parfaitement dans le sens des objectifs du projet REC "Root zone soil moisture Estimates at the daily and agricultural parcel scales for Crop irrigation management and water use impact: a multi-sensor remote sensing approach" (http://rec.isardsat.com/) qui visent à estimer l'humidité du sol dans la zone racinaire (RZSM) afin d'optimiser la gestion de l'eau d'irrigation. Des approches innovantes et prometteuses sont mises en place pour estimer l'évapotranspiration (ET), RZSM, la température de surface du sol (LST) et le stress hydrique de la végétation à travers des indices de SM dérivés des observations multispectrales à haute résolution spatio-temporelle. Les méthodologies proposées reposent sur des méthodes basées sur l'imagerie, la modélisation du transfert radiatif et la modélisation du bilan hydrique et d'énergie et sont appliquées dans une région à climat semi-aride (centre du Maroc). Dans le cadre de ma thèse, trois axes ont été explorés. Dans le premier axe, un indice de RZSM dérivé de LST-Landsat est utilisé pour estimer l'ET sur des parcelles de blé et des sols nus. L'estimation par modélisation de ET a été explorée en utilisant l'équation de Penman-monteith modifiée obtenue en introduisant une relation empirique simple entre la résistance de surface (rc) et l'indice de RZSM. Ce dernier est estimé à partir de la température de surface (LST) dérivée de Landsat, combinée avec les températures extrêmes (en conditions humides et sèches) simulée par un modèle de bilan d'énergie de surface piloté par le forçage météorologique et la fraction de couverture végétale dérivée de Landsat. La méthode utilisée est calibrée et validée sur deux parcelles de blé situées dans la même zone près de Marrakech au Maroc. Dans l'axe suivant, une méthode permettant de récupérer la SM de la surface (0-5 cm) à une résolution spatiale et temporelle élevée est développée à partir d'une synergie entre données radar (Sentinel-1) et thermique (Landsat) et en utilisant un modèle de bilan d'énergie du sol. L'approche développée a été validée sur des parcelles agricoles en sol nu et elle donne une estimation précise de la SM avec une différence quadratique moyenne en comparant à la SM in situ, égale à 0,03 m3 m-3. Dans le dernier axe, une nouvelle méthode est développée pour désagréger la MODIS LST de 1 km à 100 m de résolution en intégrant le SM proche de la surface dérivée des données radar Sentinel-1 et l'indice de végétation optique dérivé des observations Landsat. Le nouvel algorithme, qui inclut la rétrodiffusion S-1 en tant qu'entrée dans la désagrégation, produit des résultats plus stables et robustes au cours de l'année sélectionnée. Dont, 3,35 °C était le RMSE le plus bas et 0,75 le coefficient de corrélation le plus élevé évalués en utilisant le nouvel algorithme.Optimizing water management in agriculture is essential over semi-arid areas in order to preserve water resources which are already low and erratic due to human actions and climate change. This thesis aims to use the synergy of multispectral remote sensing observations (radar, optical and thermal data) for high spatio-temporal resolution monitoring of crops water needs. In this context, different approaches using various sensors (Landsat-7/8, Sentinel-1 and MODIS) have been developed to provide information on the crop Soil Moisture (SM) and water stress at a spatio-temporal scale relevant to irrigation management. This work fits well the REC "Root zone soil moisture Estimates at the daily and agricultural parcel scales for Crop irrigation management and water use impact: a multi-sensor remote sensing approach" (http://rec.isardsat.com/) project objectives, which aim to estimate the Root Zone Soil Moisture (RZSM) for optimizing the management of irrigation water. Innovative and promising approaches are set up to estimate evapotranspiration (ET), RZSM, land surface temperature (LST) and vegetation water stress through SM indices derived from multispectral observations with high spatio-temporal resolution. The proposed methodologies rely on image-based methods, radiative transfer modelling and water and energy balance modelling and are applied in a semi-arid climate region (central Morocco). In the frame of my PhD thesis, three axes have been investigated. In the first axis, a Landsat LST-derived RZSM index is used to estimate the ET over wheat parcels and bare soil. The ET modelling estimation is explored using a modified Penman-Monteith equation obtained by introducing a simple empirical relationship between surface resistance (rc) and a RZSM index. The later is estimated from Landsat-derived land surface temperature (LST) combined with the LST endmembers (in wet and dry conditions) simulated by a surface energy balance model driven by meteorological forcing and Landsat-derived fractional vegetation cover. The investigated method is calibrated and validated over two wheat parcels located in the same area near Marrakech City in Morocco. In the next axis, a method to retrieve near surface (0-5 cm) SM at high spatial and temporal resolution is developed from a synergy between radar (Sentinel-1) and thermal (Landsat) data and by using a soil energy balance model. The developed approach is validated over bare soil agricultural fields and gives an accurate estimates of near surface SM with a root mean square difference compared to in situ SM equal to 0.03 m3 m-3. In the final axis a new method is developed to disaggregate the 1 km resolution MODIS LST at 100 m resolution by integrating the near surface SM derived from Sentinel-1 radar data and the optical-vegetation index derived from Landsat observations. The new algorithm including the S-1 backscatter as input to the disaggregation, produces more stable and robust results during the selected year. Where, 3.35 °C and 0.75 were the lowest RMSE and the highest correlation coefficient assessed using the new algorithm

Prediction of topsoil properties at field-scale by using C-band SAR data

Designing and validating digital soil mapping (DSM) techniques can facilitate precision agriculture implementation. This study generates and validates a technique for the spatial prediction of soil properties based on C-band radar data. To this end, (i) we focused on working at farm-field scale and conditions, a fact scarcely reported; (ii) we validated the usefulness of Random Forest regression (RF) to predict soil properties based on C-band radar data; (iii) we validated the prediction accuracy of C-band radar data according to the coverage condition (for example: crop or fallow); and (iv) we aimed to find spatial relationship between soil apparent electrical conductivity and C-band radar. The experiment was conducted on two agricultural fields in the southern Argentine Pampas. Fifty one Sentinel 1 Level-1 GRD (Grid) products of C-band frequency (5.36 GHz) were processed. VH and VV polarizations and the dual polarization SAR vegetation index (DPSVI) were estimated. Soil information was obtained through regular-grid sample scheme and apparent soil electrical conductivity (ECa) measurements. Soil properties predicted were: texture, effective soil depth, ECa at 0-0.3m depth and ECa at 0-0.9m depth. The effect of water, vegetation and soil on the depolarization from SAR backscattering was analyzed. Complementary, spatial predictions of all soil properties from ordinary cokriging and Conditioned Latin hypercube sampling (cLHS) were evaluated using six different soil sample sizes: 20, 40, 60, 80, 100 and the total of the grid sampling scheme. The results demonstrate that the prediction accuracy of C-band SAR data for most of the soil properties evaluated varies considerably and is closely dependent on the coverage type and weather dynamics. The polarizations with high prediction accuracy of all soil properties showed low values of σVVo and σVHo, while those with low prediction accuracy showed high values of σVVo and low values of σVHo. The spatial patterns among maps of all soil properties using all samples and all sample sizes were similar. In conditions when summer crops demand large amount of water and there is soil water deficit backscattering showed higher prediction accuracy for most soil properties. During the fallow season, the prediction accuracy decreased and the spatial prediction accuracy was closely dependent on the number of validation samples. The findings of this study corroborates that DSM techniques at field scale can be achieved by using C-band SAR data. Extrapolation y applicability of this study to other areas remain to be tested.EEA BalcarceFil: Domenech, Marisa. Universidad Nacional del Sur. Departamento de Agronomía; Argentina.Fil: Amiottia, Nilda. Universidad Nacional del Sur. Departamento de Agronomía; Argentina.Fil: Amiottia, Nilda. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina.Fil: Costa, José Luis. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Balcarce; Argentina.Fil: Castro-Franco, Mauricio. Centro de Investigaciones de la Caña de Azúcar de Colombia. Estación Experimental Estación Experimental vía Cali-Florida; Colombia

On the use of temporal series of L-and X-band SAR data for soil moisture retrieval. Capitanata plain case study

This paper investigates the use of time series of ALOS/PALSAR-1 and COSMO-SkyMed data for the soil moisture retrieval (mv) by means of the SMOSAR algorithm. The application context is the exploitation of mv maps at a moderate spatial and temporal resolution for improving flood/drought monitoring at regional scale. The SAR data were acquired over the Capitanata plain in Southern Italy, over which ground campaigns were carried out in 2007, 2010 and 2011. The analysis shows that the mv retrieval accuracy is 5%-7% m^3/m^3 at L- and X band, although the latter is restricted to a use over nearly bare soil only

Surface soil moisture estimate from Sentinel-1 and Sentinel-2 data in agricultural fields in areas of high vulnerability to climate variations: the Marche region (Italy) case study

Surface soil moisture is a key hydrologic state variable that greatly influences the global environment and human society. Its significant decrease in the Mediterranean region, registered since the 1950s, and expected to continue in the next century, threatens soil health and crops. Microwave remote sensing techniques are becoming a key tool for the implementation of climate-smart agriculture, as a means for surface soil moisture retrieval that exploits the correlation between liquid water and the dielectric properties of soil. In this study, a workflow in Google Earth Engine was developed to estimate surface soil moisture in the agricultural fields of the Marche region (Italy) through Synthetic Aperture Radar data. Firstly, agricultural areas were extracted with both Sentinel-2 optical and Sentinel-1 radar satellites, investigating the use of Dual-Polarimetric Entropy-Alpha decomposition's bands to improve the accuracy of radar data classification. The results show that Entropy and Alpha bands improve the kappa index obtained from the radar data only by 4% (K = 0.818), exceeding optical accuracy in urban and water areas. However, they still did not allow to reach the overall optical accuracy (K = 0.927). The best classification results are reached with the total dataset (K = 0.949). Subsequently, Water Cloud and Tu Wien models were implemented on the crop areas using calibration parameters derived from literature, to test if an acceptable accuracy is reached without in situ observation. While the first model’s accuracy was inadequate (RMSD = 12.3), the extraction of surface soil moisture using Tu Wien change detection method was found to have acceptable accuracy (RMSD = 9.4)