Sea Surface Salinity Retrievals from Aquarius Using Neural Networks

Even though the Sea Surface Salinity (SSS) retrieved from Aquarius are generally very close to in-situ measurements, the level of similarity varies with the region and with the circumstances of the observations (wind speed, sea surface temperature, etc.). SSS is currently retrieved from the brightness temperatures measured by Aquarius and applying the current theoretical model for the propagation and emission of the natural thermal radiation. In this contribution we consider an alternative retrieval approach based on a Neural Network (NN) with the goal of improving the subsets of Aquarius SSS data that are in poorer agreement within-situ measurements. The subset considered here are the SSS retrieved at latitudes higher than 30 . The output of the NN approach are compared against in-situ measurements using four statistical metrics (correlation coefficient, bias, RMSD and 5% trimmed range). The output of the NN and the nominal Aquarius SSS are compared against SSS values from in-situ measurements and from ocean models. From these comparisons it appears that the output of the NN matches the in-situ measurements better than the nominal Aquarius SSS

An artificial neural network approach for soil moisture retrieval using passive microwave data

Soil moisture is a key variable that defines land surface-atmosphere (boundary layer) interactions, by contributing directly to the surface energy and water balance. Soil moisture values derived from remote sensing platforms only accounts for the near surface soil layers, generally the top 5cm. Passive microwave data at L-band (1.4 GHz, 21cm wavelength) measurements are shown to be a very effective observation for surface soil moisture retrieval. The first space-borne L-band mission dedicated to observing soil moisture, the European Space Agency's (ESA) Soil Moisture and Ocean Salinity (SMOS) mission, was launched on 2nd November 2009.Artificial Neural Network (ANN) methods have been used to empirically ascertain the complex statistical relationship between soil moisture and brightness temperature in the presence of vegetation cover. The current problem faced by this method is its inability to predict soil moisture values that are 'out-of-range' of the training data.In this research, an optimization model is developed for the Backpropagation Neural Network model. This optimization model utilizes the combination of the mean and standard deviation of the soil moisture values, together with the prediction process at different pre-determined, equal size regions to cope with the spatial and temporal variation of soil moisture values. This optimized model coupled with an ANN of optimum architecture, in terms of inputs and the number of neurons in the hidden layers, is developed to predict scale-to-scale and downscaling of soil moisture values. The dependency on the accuracy of the mean and standard deviation values of soil moisture data is also studied in this research by simulating the soil moisture values using a multiple regression model. This model obtains very encouraging results for these research problems.The data used to develop and evaluate the model in this research has been obtained from the National Airborne Field Experiments in 2005

Input parameters selection for soil moisture retrieval using an artificial neural network

Factors other than soil moisture which influence the intensity of microwave emission from the soil include surface temperature, surface roughness, vegetation cover and soil texture which make this a non-linear and ill-posed problem. Artificial Neural Networks (ANNs) have been demonstrated to be good solutions to this type of problem. Since an ANN is a data driven model, proper input selection is a crucial step in its implementation as the presence of redundant or unnecessary inputs can severely impair the ability of the network to learn the target patterns. In this paper, the input parameters are chosen in combination with the brightness temperatures and are based on the use of incremental contributions of the variables towards soil moisture retrieval. Field experiment data obtained during the National Airborne Field Experiment 2005 (NAFE'05) are used. The retrieval accuracy with the input parameters selected is compared with the use of only brightness temperature as input and the use of brightness temperature in conjunction with a range of available parameters. Note that this research does not aim at selecting the best features for all ANN soil moisture retrieval problems using passive microwave. The paper shows that, depending on the problem and the nature of the data, some of the data available are redundant as the input of ANN for soil moisture retrieval. Importantly the results show that with the appropriate choice of inputs, the soil moisture retrieval accuracy of ANN can be significantly improved

Soil moisture deficit estimation using satellite multi-angle brightness temperature

Accurate soil moisture information is critically important for hydrological modelling. Although remote sensing soil moisture measurement has become an important data source, it cannot be used directly in hydrological modelling. A novel study based on nonlinear techniques (a local linear regression (LLR) and two feedforward artificial neural networks (ANNs)) is carried out to estimate soil moisture deficit (SMD), using the Soil Moisture and Ocean Salinity (SMOS) multi-angle brightness temperatures (Tbs) with both horizontal (H) and vertical (V) polarisations. The gamma test is used for the first time to determine the optimum number of Tbs required to construct a reliable smooth model for SMD estimation, and the relationship between model input and output is achieved through error variance estimation. The simulated SMD time series in the study area is from the Xinanjiang hydrological model. The results have shown that LLR model is better at capturing the interrelations between SMD and Tbs than ANNs, with outstanding statistical performances obtained during both training (NSE = 0.88, r = 0.94, RMSE = 0.008 m) and testing phases (NSE = 0.85, r = 0.93, RMSE = 0.009 m). Nevertheless, both ANN training algorithms (radial BFGS and conjugate gradient) have performed well in estimating the SMD data and showed excellent performances compared with those derived directly from the SMOS soil moisture products. This study has also demonstrated the informative capability of the gamma test in the input data selection for model development. These results provide interesting perspectives for data-assimilation in flood-forecasting

Use of Soil Moisture Variability in Artificial Neural Network Retrieval of Soil Moisture

Passive microwave remote sensing is one of the most promising techniques for soil moisture retrieval. However, the inversion of soil moisture from brightness temperature observations is not straightforward, as it is influenced by numerous factors such as surface roughness, vegetation cover, and soil texture. Moreover, the relationship between brightness temperature, soil moisture and the factors mentioned above is highly non-linear and ill-posed. Consequently, Artificial Neural Networks (ANNs) have been used to retrieve soil moisture from microwave data, but with limited success when dealing with data different to that from the training period. In this study, an ANN is tested for its ability to predict soil moisture at 1 km resolution on different dates following training at the same site for a specific date. A novel approach that utilizes information on the variability of soil moisture, in terms of its mean and standard deviation for a (sub) region of spatial dimension up to 40 km, is used to improve the current retrieval accuracy of the ANN method.A comparison between the ANN with and without the use of the variability information showed that this enhancement enables the ANN to achieve an average Root Mean Square Error (RMSE) of around 5.1% v/v when using the variability information, as compared to around 7.5% v/v without it. The accuracy of the soil moisture retrieval was further improved by the division of the target site into smaller regions down to 4 km in size, with the spatial variability of soil moisture calculated from within the smaller region used in the ANN. With the combination of an ANN architecture of a single hidden layer of 20 neurons and the dual-polarized brightness temperatures as input, the proposed use of variability and sub-region methodology achieves an average retrieval accuracy of 3.7% v/v. Although this accuracy is not the lowest as comparing to the research in this field, the main contribution is the ability of ANN in solving the problem of predicting “out-of-range” soil moisture values. However, the applicability of this method is highly dependent on the accuracy of the mean and standard deviation values within the sub-region, potentially limiting its routine application

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

Estimating Surface Soil Moisture from SMAP Observations Using a Neural Network Technique

A Neural Network (NN) algorithm was developed to estimate global surface soil moisture for April 2015 to March 2017 with a 2-3 day repeat frequency using passive microwave observations from the Soil Moisture Active Passive (SMAP) satellite, surface soil temperatures from the NASA Goddard Earth Observing System Model version 5 (GEOS-5) land modeling system, and Moderate Resolution Imaging Spectroradiometer-based vegetation water content. The NN was trained on GEOS-5 soil moisture target data, making the NN estimates consistent with the GEOS-5 climatology, such that they may ultimately be assimilated into this model without further bias correction. Evaluated against in situ soil moisture measurements, the average unbiased root mean square error (ubRMSE), correlation and anomaly correlation of the NN retrievals were 0.037 m(exp. 3)m(exp. -3), 0.70 and 0.66, respectively, against SMAP core validation site measurements and 0.026 m(exp. 3)m(exp. -3), 0.58 and 0.48, respectively, against International Soil Moisture Network (ISMN) measurements. At the core validation sites, the NN retrievals have a significantly higher skill than the GEOS-5 model estimates and a slightly lower correlation skill than the SMAP Level-2 Passive (L2P) product. The feasibility of the NN method was reflected by a lower ubRMSE compared to the L2P retrievals as well as a higher skill when ancillary parameters in physically-based retrievals were uncertain. Against ISMN measurements, the skill of the two retrieval products was more comparable. A triple collocation analysis against Advanced Microwave Scanning Radiometer 2 (AMSR2) and Advanced Scatterometer (ASCAT) soil moisture retrievals showed that the NN and L2P retrieval errors have a similar spatial distribution, but the NN retrieval errors are generally lower in densely vegetated regions and transition zones

Soil moisture estimation synergy using GNSS-R and L-Band microwave radiometry data from FSSCat/FMPL-2

The Federated Satellite System mission (FSSCat) was the winner of the 2017 Copernicus Masters Competition and the first Copernicus third-party mission based on CubeSats. One of FSSCat’s objectives is to provide coarse Soil Moisture (SM) estimations by means of passive microwave measurements collected by Flexible Microwave Payload-2 (FMPL-2). This payload is a novel CubeSat based instrument combining an L1/E1 Global Navigation Satellite Systems-Reflectometer (GNSS-R) and an L-band Microwave Radiometer (MWR) using software-defined radio. This work presents the first results over land of the first two months of operations after the commissioning phase, from 1 October to 4 December 2020. Four neural network algorithms are implemented and analyzed in terms of different sets of input features to yield maps of SM content over the Northern Hemisphere (latitudes above 45° N). The first algorithm uses the surface skin temperature from the European Centre of Medium-Range Weather Forecast (ECMWF) in conjunction with the 16 day averaged Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) to estimate SM and to use it as a comparison dataset for evaluating the additional models. A second approach is implemented to retrieve SM, which complements the first model using FMPL-2 L-band MWR antenna temperature measurements, showing a better performance than in the first case. The error standard deviation of this model referred to the Soil Moisture and Ocean Salinity (SMOS) SM product gridded at 36 km is 0.074 m3/m3. The third algorithm proposes a new approach to retrieve SM using FMPL-2 GNSS-R data. The mean and standard deviation of the GNSS-R reflectivity are obtained by averaging consecutive observations based on a sliding window and are further included as additional input features to the network. The model output shows an accurate SM estimation compared to a 9 km SMOS SM product, with an error of 0.087 m3/m3. Finally, a fourth model combines MWR and GNSS-R data and outperforms the previous approaches, with an error of just 0.063 m3/m3. These results demonstrate the capabilities of FMPL-2 to provide SM estimates over land with a good agreement with respect to SMOS SM.This work was supported by the 2017 ESA S3 challenge and Copernicus Masters overall winner award (“FSSCat” project). This work was (partially) sponsored by project SPOT: Sensing with Pioneering Opportunistic Techniques grant RTI2018-099008-B-C21 / AEI / 10.13039/501100011033, and by the Unidad de Excelencia Maria de Maeztu MDM-2016-0600. This work was also (partially) sponsored by the Spanish Ministry of Science and Innovation through the project ESP2017-89463-C3, by the Centro de Excelencia Severo Ochoa (CEX2019-000928-S), and by the CSIC Plataforma Temática Interdisciplinar de Teledetección (PTI-Teledetect). Joan Francesc Munoz-Martin received support from the grant for the recruitment of early-stage research staff FI-DGR 2018 of the AGAUR - Generalitat de Catalunya (FEDER), Spain; Christoph Herbert received the support of a fellowship from “la Caixa” Foundation (ID 100010434) with the fellowship code LCF/BQ/DI18/11660050 and funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie Grant Agreement No. 713673; David Llavería received support from an FPU fellowship from the Spanish Ministry of Education FPU18/06107.Peer ReviewedPostprint (published version

The impact of SMOS soil moisture data assimilation within the Operational Global Flood Awareness System (GloFAS)

In this study the impacts of Soil Moisture and Ocean Salinity (SMOS) soil moisture data assimilation upon the streamflow prediction of the operational Global Flood Awareness System (GloFAS) were investigated. Two GloFAS experiments were performed, one which used hydro-meteorological forcings produced with the assimilation of the SMOS data, the other using forcings which excluded the assimilation of the SMOS data. Both sets of experiment results were verified against streamflow observations in the United States and Australia. Skill scores were computed for each experiment against the observation datasets, the differences in the skill scores were used to identify where GloFAS skill may be affected by the assimilation of SMOS soil moisture data. In addition, a global assessment was made of the impact upon the 5th and 95th GloFAS flow percentiles to see how SMOS data assimilation affected low and high flows respectively. Results against in-situ observations found that GloFAS skill score was only affected by a small amount. At a global scale, the results showed a large impact on high flows in areas such as the Hudson Bay, central United States, the Sahel and Australia. There was no clear spatial trend to these differences as opposing signs occurred within close proximity to each other. Investigating the differences between the simulations at individual gauging stations showed that they often only occurred during a single flood event; for the remainder of the simulation period the experiments were almost identical. This suggests that SMOS data assimilation may affect the generation of surface runoff during high flow events, but may have less impact on baseflow generation during the remainder of the hydrograph. To further understand this, future work could assess the impact of SMOS data assimilation upon specific hydrological components such as surface and subsurface runoff