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%

Evaluation of a global soil moisture product from finer spatial resolution sar data and ground measurements at Irish sites

In the framework of the European Space Agency Climate Change Initiative, a global, almost daily, soil moisture (SM) product is being developed from passive and active satellite microwave sensors, at a coarse spatial resolution. This study contributes to its validation by using finer spatial resolution ASAR Wide Swath and in situ soil moisture data taken over three sites in Ireland, from 2007 to 2009. This is the first time a comparison has been carried out between three sets of independent observations from different sensors at very different spatial resolutions for such a long time series. Furthermore, the SM spatial distribution has been investigated at the ASAR scale within each Essential Climate Variable (ECV) pixel, without adopting any particular model or using a densely distributed network of in situ stations. This approach facilitated an understanding of the extent to which geophysical factors, such as soil texture, terrain composition and altitude, affect the retrieved ECV SM product values in temperate grasslands. Temporal and spatial variability analysis provided high levels of correlation (p < 0.025) and low errors between the three datasets, leading to confidence in the new ECV SM global product, despite limitations in its ability to track the driest and wettest conditions

Dynamic Programming and Bayesian Inference

Dynamic programming and Bayesian inference have been both intensively and extensively developed during recent years. Because of these developments, interest in dynamic programming and Bayesian inference and their applications has greatly increased at all mathematical levels. The purpose of this book is to provide some applications of Bayesian optimization and dynamic programming

Surface Soil Moisture Retrievals from Remote Sensing:Current Status, Products & Future Trends

Advances in Earth Observation (EO) technology, particularly over the last two decades, have shown that soil moisture content (SMC) can be measured to some degree or other by all regions of the electromagnetic spectrum, and a variety of techniques have been proposed to facilitate this purpose. In this review we provide a synthesis of the efforts made during the last 20 years or so towards the estimation of surface SMC exploiting EO imagery, with a particular emphasis on retrievals from microwave sensors. Rather than replicating previous overview works, we provide a comprehensive and critical exploration of all the major approaches employed for retrieving SMC in a range of different global ecosystems. In this framework, we consider the newest techniques developed within optical and thermal infrared remote sensing, active and passive microwave domains, as well as assimilation or synergistic approaches. Future trends and prospects of EO for the accurate determination of SMC from space are subject to key challenges, some of which are identified and discussed within. It is evident from this review that there is potential for more accurate estimation of SMC exploiting EO technology, particularly so, by exploring the use of synergistic approaches between a variety of EO instruments. Given the importance of SMC in Earth’s land surface interactions and to a large range of applications, one can appreciate that its accurate estimation is critical in addressing key scientific and practical challenges in today’s world such as food security, sustainable planning and management of water resources. The launch of new, more sophisticated satellites strengthens the development of innovative research approaches and scientific inventions that will result in a range of pioneering and ground-breaking advancements in the retrievals of soil moisture from space

Power Sensitivity Analysis of Multi-Frequency, Multi-Polarized, Multi-Temporal SAR Data for Soil-Vegetation System Variables Characterization

The knowledge of spatial and temporal variability of soil water content and others soil-vegetation variables (leaf area index, fractional cover) assumes high importance in crop management. Where and when the cloudiness limits the use of optical and thermal remote sensing techniques, synthetic aperture radar (SAR) imagery has proven to have several advantages (cloud penetration, day/night acquisitions and high spatial resolution). However, measured backscattering is controlled by several factors including SAR configuration (acquisition geometry, frequency and polarization), and target dielectric and geometric properties. Thus, uncertainties arise about the more suitable configuration to be used. With the launch of the ALOS Palsar, Cosmo-Skymed and Sentinel 1 sensors, a dataset of multi-frequency (X, C, L) and multi-polarization (co- and cross-polarizations) images are now available from a virtual constellation; thus, significant issues concerning the retrieval of soil-vegetation variables using SAR are: (i) identifying the more suitable SAR configuration; (ii) understanding the affordability of a multi-frequency approach. In 2006, a vast dataset of both remotely sensed images (SAR and optical/thermal) and in situ data was collected in the framework of the AgriSAR 2006 project funded by ESA and DLR. Flights and sampling have taken place weekly from April to August. In situ data included soil water content, soil roughness, fractional coverage and Leaf Area Index (LAI). SAR airborne data consisted of multi-frequency and multi-polarized SAR images (X, C and L frequencies and HH, HV, VH and VV polarizations). By exploiting this very wide dataset, this paper, explores the capabilities of SAR in describing four of the main soil-vegetation variables (SVV). As a first attempt, backscattering and SVV temporal behaviors are compared (dynamic analysis) and single-channel regressions between backscattering and SVV are analyzed. Remarkably, no significant correlations were found between backscattering and soil roughness (over both bare and vegetated plots), whereas it has been noticed that the contributions of water content of soil underlying the vegetation often did not influence the backscattering (depending on canopy structure and SAR configuration). Most significant regressions were found between backscattering and SVV characterizing the vegetation biomass (fractional cover and LAI). Secondly, the effect of SVV changes on the spatial correlation among SAR channels (accounting for different polarization and/or frequencies) was explored. An inter-channel spatial/temporal correlation analysis is proposed by temporally correlating two-channel spatial correlation and SVV. This novel approach allowed a widening in the number of significant correlations and their strengths by also encompassing the use of SAR data acquired at two different frequencie

Power Sensitivity Analysis of Multi-Frequency, Multi-Polarized, Multi-Temporal SAR Data for Soil-Vegetation System Variables Characterization

The knowledge of spatial and temporal variability of soil water content and others soil-vegetation variables (leaf area index, fractional cover) assumes high importance in crop management. Where and when the cloudiness limits the use of optical and thermal remote sensing techniques, synthetic aperture radar (SAR) imagery has proven to have several advantages (cloud penetration, day/night acquisitions and high spatial resolution). However, measured backscattering is controlled by several factors including SAR configuration (acquisition geometry, frequency and polarization), and target dielectric and geometric properties. Thus, uncertainties arise about the more suitable configuration to be used. With the launch of the ALOS Palsar, Cosmo-Skymed and Sentinel 1 sensors, a dataset of multi-frequency (X, C, L) and multi-polarization (co- and cross-polarizations) images are now available from a virtual constellation; thus, significant issues concerning the retrieval of soil-vegetation variables using SAR are: (i) identifying the more suitable SAR configuration; (ii) understanding the affordability of a multi-frequency approach. In 2006, a vast dataset of both remotely sensed images (SAR and optical/thermal) and in situ data was collected in the framework of the AgriSAR 2006 project funded by ESA and DLR. Flights and sampling have taken place weekly from April to August. In situ data included soil water content, soil roughness, fractional coverage and Leaf Area Index (LAI). SAR airborne data consisted of multi-frequency and multi-polarized SAR images (X, C and L frequencies and HH, HV, VH and VV polarizations). By exploiting this very wide dataset, this paper, explores the capabilities of SAR in describing four of the main soil-vegetation variables (SVV). As a first attempt, backscattering and SVV temporal behaviors are compared (dynamic analysis) and single-channel regressions between backscattering and SVV are analyzed. Remarkably, no significant correlations were found between backscattering and soil roughness (over both bare and vegetated plots), whereas it has been noticed that the contributions of water content of soil underlying the vegetation often did not influence the backscattering (depending on canopy structure and SAR configuration). Most significant regressions were found between backscattering and SVV characterizing the vegetation biomass (fractional cover and LAI). Secondly, the effect of SVV changes on the spatial correlation among SAR channels (accounting for different polarization and/or frequencies) was explored. An inter-channel spatial/temporal correlation analysis is proposed by temporally correlating two-channel spatial correlation and SVV. This novel approach allowed a widening in the number of significant correlations and their strengths by also encompassing the use of SAR data acquired at two different frequencie

Neural Networks Applications for the Remote Sensing of Hydrological Parameters

The main artificial neural networks (ANN)‐based retrieval algorithms developed at the Institute of Applied Physics (IFAC) are reviewed here. These algorithms aim at retrieving the main hydrological parameters, namely the soil moisture content (SMC), the plant water content (PWC) of agricultural vegetation, the woody volume of forests (WV) and the snow depth (SD) or snow water equivalent (SWE), from data collected by active (SAR/scatterometers) and passive (radiometers) microwave sensors operating from space. Taking advantage of the fast computation, ANN are able to generate output maps of the target parameter at both local and global scales, with a resolution varying from hundreds of meters to tens of kilometres, depending on the considered sensor. A peculiar strategy adopted for the training, which has been obtained by combining satellite measurements with data simulated by electromagnetic models (based on the radiative transfer theory, RTT), made these algorithms robust and site independent. The obtained results demonstrated that ANN are a powerful tool for estimating the hydrological parameters at different spatial scales, provided that they have been trained with consistent datasets, made up by both experimental and theoretical data

Advanced methods for earth observation data synergy for geophysical parameter retrieval

The first part of the thesis focuses on the analysis of relevant factors to estimate the response time between satellite-based and in-situ soil moisture (SM) using a Dynamic Time Warping (DTW). DTW was applied to the SMOS L4 SM, and was compared to in-situ root-zone SM in the REMEDHUS network in Western Spain. The method was customized to control the evolution of time lag during wetting and drying conditions. Climate factors in combination with crop growing seasons were studied to reveal SM-related processes. The heterogeneity of land use was analyzed using high-resolution images of NDVI from Sentinel-2 to provide information about the level of spatial representativity of SMOS data to each in-situ station. The comparison of long-term precipitation records and potential evapotranspiration allowed estimation of SM seasons describing different SM conditions depending on climate and soil properties. The second part of the thesis focuses on data-driven methods for sea ice segmentation and parameter retrieval. A Bayesian framework is employed to segment sets of multi-source satellite data. The Bayesian unsupervised learning algorithm allows to investigate the ‘hidden link’ between multiple data. The statistical properties are accounted for by a Gaussian Mixture Model, and the spatial interactions are reflected using Hidden Markov Random Fields. The algorithm segments spatial data into a number of classes, which are represented as a latent field in physical space and as clusters in feature space. In a first application, a two-step probabilistic approach based on Expectation-Maximization and the Bayesian segmentation algorithm was used to segment SAR images to discriminate surface water from sea ice types. Information on surface roughness is contained in the radar backscattering images which can be - in principle - used to detect melt ponds and to estimate high-resolution sea ice concentration (SIC). In a second study, the algorithm was applied to multi-incidence angle TB data from the SMOS L1C product to harness the its sensitivity to thin ice. The spatial patterns clearly discriminate well-determined areas of open water, old sea ice and a transition zone, which is sensitive to thin sea ice thickness (SIT) and SIC. In a third application, SMOS and the AMSR2 data are used to examine the joint effect of CIMR-like observations. The information contained in the low-frequency channels allows to reveal ranges of thin sea ice, and thicker ice can be determined from the relationship between the high-frequency channels and changing conditions as the sea ice ages. The proposed approach is suitable for merging large data sets and provides metrics for class analysis, and to make informed choices about integrating data from future missions into sea ice products. A regression neural network approach was investigated with the goal to infer SIT using TB data from the Flexible Microwave Payload 2 (FMPL-2) of the FSSCat mission. Two models - covering thin ice up to 0.6m and the full-range of SIT - were trained on Arctic data using ground truth data derived from the SMOS and Cryosat-2. This work demonstrates that moderate-cost CubeSat missions can provide valuable data for applications in Earth observation.La primera parte de la tesis se centra en el análisis de los factores relevantes para estimar el tiempo de respuesta entre la humedad del suelo (SM) basada en el satélite y la in-situ, utilizando una deformación temporal dinámica (DTW). El DTW se aplicó al SMOS L4 SM, y se comparó con la SM in-situ en la red REMEDHUS en el oeste de España. El método se adaptó para controlar la evolución del desfase temporal durante diferentes condiciones de humedad y secado. Se estudiaron los factores climáticos en combinación con los períodos de crecimiento de los cultivos para revelar los procesos relacionados con la SM. La heterogeneidad del uso del suelo se analizó utilizando imágenes de alta resolución de NDVI de Sentinel-2 para proporcionar información sobre el nivel de representatividad espacial de los datos de SMOS a cada estación in situ. La comparación de los patrones de precipitación a largo plazo y la evapotranspiración potencial permitió estimar las estaciones de SM que describen diferentes condiciones de SM en función del clima y las propiedades del suelo. La segunda parte de esta tesis se centra en métodos dirigidos por datos para la segmentación del hielo marino y la obtención de parámetros. Se emplea un método de inferencia bayesiano para segmentar conjuntos de datos satelitales de múltiples fuentes. El algoritmo de aprendizaje bayesiano no supervisado permite investigar el “vínculo oculto” entre múltiples datos. Las propiedades estadísticas se contabilizan mediante un modelo de mezcla gaussiana, y las interacciones espaciales se reflejan mediante campos aleatorios ocultos de Markov. El algoritmo segmenta los datos espaciales en una serie de clases, que se representan como un campo latente en el espacio físico y como clústeres en el espacio de las variables. En una primera aplicación, se utilizó un enfoque probabilístico de dos pasos basado en la maximización de expectativas y el algoritmo de segmentación bayesiano para segmentar imágenes SAR con el objetivo de discriminar el agua superficial de los tipos de hielo marino. La información sobre la rugosidad de la superficie está contenida en las imágenes de backscattering del radar, que puede utilizarse -en principio- para detectar estanques de deshielo y estimar la concentración de hielo marino (SIC) de alta resolución. En un segundo estudio, el algoritmo se aplicó a los datos TB de múltiples ángulos de incidencia del producto SMOS L1C para aprovechar su sensibilidad al hielo fino. Los patrones espaciales discriminan claramente áreas bien determinadas de aguas abiertas, hielo marino viejo y una zona de transición, que es sensible al espesor del hielo marino fino (SIT) y al SIC. En una tercera aplicación, se utilizan los datos de SMOS y de AMSR2 para examinar el efecto conjunto de las observaciones tipo CIMR. La información contenida en los canales de baja frecuencia permite revelar rangos de hielo marino delgado, y el hielo más grueso puede determinarse a partir de la relación entre los canales de alta frecuencia y las condiciones cambiantes a medida que el hielo marino envejece. El enfoque propuesto es adecuado para fusionar grandes conjuntos de datos y proporciona métricas para el análisis de clases, y para tomar decisiones informadas sobre la integración de datos de futuras misiones en los productos de hielo marino. Se investigó un enfoque de red neuronal de regresión con el objetivo de inferir el SIT utilizando datos de TB de la carga útil de microondas flexible 2 (FMPL-2) de la misión FSSCat. Se entrenaron dos modelos - que cubren el hielo fino hasta 0.6 m y el rango completo del SIT - con datos del Ártico utilizando datos de “ground truth” derivados del SMOS y del Cryosat-2. Este trabajo demuestra que las misiones CubeSat de coste moderado pueden proporcionar datos valiosos para aplicaciones de observación de la Tierra.Postprint (published version