Bayesian unsupervised machine learning approach to segment arctic sea ice using SMOS data

Microwave radiometry at L-band is sensitive to sea ice thickness (SIT) up to ~ 60 cm. Current methods to infer SIT depend on ice-physical properties and data provided by the ESA’s Soil Moisture and Ocean Salinity (SMOS) mission. However, retrieval accuracy is limited due to seasonally and regionally variable surface conditions during the formation and melting of sea ice. In this work, Arctic sea ice is segmented using a Bayesian unsupervised learning algorithm aiming to recognize spatial patterns by harnessing multi-incidence angle brightness temperature observations. The approach considers both statistical characteristics and spatial correlations of the observations. The temporal stability and separability of classes are analyzed to distinguish ambiguous from well-determined regions. Model uncertainty is quantified from class membership probabilities using information entropy. The presented approach opens up a new scope to improve current SIT retrieval algorithms, and can be particularly beneficial to investigate merged satellite products.This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No.713673. It was also funded through the award “Unidad de Excelencia María de Maeztu” MDM-2016-0600, by the Spanish Ministry of Science and Innovation through the project “L-band” ESP2017-89463-C3-2-R, and the project “Sensing with Pioneering Opportunistic Techniques (SPOT)” RTI2018-099008-B-C21/AEI/10.13039/501100011033.Peer ReviewedPostprint (published version

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

Variability and uncertainty of satellite sea surface salinity in the subpolar North Atlantic (2010-2019)

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Yu, L. Variability and uncertainty of satellite sea surface salinity in the subpolar North Atlantic (2010-2019). Remote Sensing, 12(13), (2020): 2092, doi:10.3390/rs12132092.Satellite remote sensing of sea surface salinity (SSS) in the recent decade (2010–2019) has proven the capability of L-band (1.4 GHz) measurements to resolve SSS spatiotemporal variability in the tropical and subtropical oceans. However, the fidelity of SSS retrievals in cold waters at mid-high latitudes has yet to be established. Here, four SSS products derived from two satellite missions were evaluated in the subpolar North Atlantic Ocean in reference to two in situ gridded products. Harmonic analysis of annual and semiannual cycles in in situ products revealed that seasonal variations of SSS are dominated by an annual cycle, with a maximum in March and a minimum in September. The annual amplitudes are larger (>0.3 practical salinity scale (pss)) in the western basin where surface waters are colder and fresher, and weaker (~0.06 pss) in the eastern basin where surface waters are warmer and saltier. Satellite SSS products have difficulty producing the right annual cycle, particularly in the Labrador/Irminger seas where the SSS seasonality is dictated by the influx of Arctic low-salinity waters along the boundary currents. The study also found that there are basin-scale, time-varying drifts in the decade-long SMOS data records, which need to be corrected before the datasets can be used for studying climate variability of SSSThis research was funded by NASA Ocean Salinity Science Team (OSST) activities through Grant 80NSSC18K1335

2000 days of SMOS at the Barcelona Expert Centre: a tribute to the work of Jordi Font

Soil Moisture and Ocean Salinity (SMOS) is the first satellite mission capable of measuring sea surface salinity and soil moisture from space. Its novel instrument (the L-band radiometer MIRAS) has required the development of new algorithms to process SMOS data, a challenging task due to many processing issues and the difficulties inherent in a new technology. In the wake of SMOS, a new community of users has grown, requesting new products and applications, and extending the interest in this novel brand of satellite services. This paper reviews the role played by the Barcelona Expert Centre under the direction of Jordi Font, SMOS co-principal investigator. The main scientific activities and achievements and the future directions are discussed, highlighting the importance of the oceanographic applications of the mission.Peer ReviewedPostprint (published version

First SMOS Sea Surface Salinity dedicated products over the Baltic Sea

26 pages, 24 figures, 4 tables.-- Data availability: Access to the data is provided by the Barcelona Expert Center, through its FTP service. The DOI of the L3 product is https://doi.org/10.20350/digitalCSIC/13859 (González-Gambau et al., 2021a). The DOI of the L4 product is https://doi.org/10.20350/digitalCSIC/13860 (González-Gambau et al., 2021b). Seasonal averaged L4 SSS products are also available in the HELCOM catalogue (https://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/9d979033-1136-4dd1-a09b-7ee9e512ad14, BEC team, 2021b), and they can be visualized in the HELCOM Map and Data service (https://maps.helcom.fi/website/mapservice/?datasetID=9d979033-1136-4dd1-a09b-7ee9e512ad14, last access: 9 November 2021).-- This work is a contribution to the CSIC Thematic Interdisciplinary Platform TeledetectThis paper presents the first Soil Moisture and Ocean Salinity (SMOS) Sea Surface Salinity (SSS) dedicated products over the Baltic Sea. The SSS retrieval from L-band brightness temperature (TB) measurements over this basin is really challenging due to important technical issues, such as the land–sea and ice–sea contamination, the high contamination by radio-frequency interference (RFI) sources, the low sensitivity of L-band TB at SSS changes in cold waters, and the poor characterization of dielectric constant models for the low SSS range in the basin. For these reasons, exploratory research in the algorithms used from the level 0 up to level 4 has been required to develop these dedicated products. This work has been performed in the framework of the European Space Agency regional initiative Baltic+ Salinity Dynamics. Two Baltic+ SSS products have been generated for the period 2011–2019 and are freely distributed: the Level 3 (L3) product (daily generated 9 d maps in a 0.25∘ grid; https://doi.org/10.20350/digitalCSIC/13859, González-Gambau et al., 2021a) and the Level 4 (L4) product (daily maps in a 0.05∘ grid; https://doi.org/10.20350/digitalCSIC/13860, González-Gambau et al., 2021b)​​​​​​​, which are computed by applying multifractal fusion to L3 SSS with SST maps. The accuracy of L3 SSS products is typically around 0.7–0.8 psu. The L4 product has an improved spatiotemporal resolution with respect to the L3 and the accuracy is typically around 0.4 psu. Regions with the highest errors and limited coverage are located in Arkona and Bornholm basins and Gulfs of Finland and Riga. The impact assessment of Baltic+ SSS products has shown that they can help in the understanding of salinity dynamics in the basin. They complement the temporally and spatially very sparse in situ measurements, covering data gaps in the region, and they can also be useful for the validation of numerical models, particularly in areas where in situ data are very sparseThis work has been carried out as part of the Baltic+ Salinity Dynamics project (4000126102/18/I-BG), funded by the European Space Agency. It has been also supported in part by the Spanish R&D project INTERACT (PID2020-114623RB-C31), which is funded by MCIN/AEI/10.13039/501100011033. We also received funding from the Spanish government through the “Severo Ochoa Centre of Excellence” accreditation (CEX2019-000928-S)Peer reviewe

Improved BEC SMOS Arctic Sea Surface Salinity product v3.1

17 pages, 13 figures, 1 table.-- Data availability: The product (Martínez et al., 2019) is freely distributed on the BEC (Barcelona Expert Center) web page (http://bec.icm.csic.es/, last access: 25 January 2022) with the DOI number https://doi.org/10.20350/digitalCSIC/12620 (Martínez et al., 2019) and on the Digital CSIC server: https://digital.csic.es/handle/10261/219679 (last access: 25 January 2022). Data can be downloaded from the FTP service: http://bec.icm.csic.es/bec-ftp-service/ (last access: 25 January 2022). The maps are distributed in the standard grid EASE-Grid 2.0, which has a spatial resolution of 25 km. In addition to the product validated in this work (L3 with temporal resolution of 9 d), L3 products having a temporal resolution of 3 and 18 d and the L2 product are available. These Arctic SSS products cover the period from 2011 to 2019.-- This work represents a contribution to the CSIC Thematic Interdisciplinary Platform PTI Teledetect and PolarCSIC. Argo data were collected and made freely available by the International Argo program and the national programs that contribute to it (https://argo.ucsd.edu, https://www.ocean-ops.org, last access: 25 January 2022). The Argo program is part of the Global Ocean Observing SystemMeasuring salinity from space is challenging since the sensitivity of the brightness temperature (TB) to sea surface salinity (SSS) is low (about 0.5 K psu−1), while the SSS range in the open ocean is narrow (about 5 psu, if river discharge areas are not considered). This translates into a high accuracy requirement of the radiometer (about 2–3 K). Moreover, the sensitivity of the TB to SSS at cold waters is even lower (0.3 K psu−1), making the retrieval of the SSS in the cold waters even more challenging. Due to this limitation, the ESA launched a specific initiative in 2019, the Arctic+Salinity project (AO/1-9158/18/I-BG), to produce an enhanced Arctic SSS product with better quality and resolution than the available products. This paper presents the methodologies used to produce the new enhanced Arctic SMOS SSS product (Martínez et al., 2019) . The product consists of 9 d averaged maps in an EASE 2.0 grid of 25 km. The product is freely distributed from the Barcelona Expert Center (BEC, http://bec.icm.csic.es/, last access: 25 January 2022) with the DOI number https://doi.org/10.20350/digitalCSIC/12620 (Martínez et al., 2019). The major change in this new product is its improvement of the effective spatial resolution that permits better monitoring of the mesoscale structures (larger than 50 km), which benefits the river discharge monitoringThis work has been carried out as part of the ESA Arctic+Salinity project (AO/1-9158/18/I-BG), which permitted the production of the database, and the Ministry of Economy and Competitiveness, Spain, through the National R&D Plan under L-BAND project ESP2017-89463-C3-1-R. [...] With the funding support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), of the Spanish Research Agency (AEI)Peer reviewe

Forecasting Salinity in the Laguna Madre Using Deep Learning

Salinity is an important metric in the Laguna Madre for establishing the long term health of the local ecological population. By utilizing Deep Learning (DL) techniques, the predicted and forecasted salinity in the Laguna Madre is generated from data provided by the Moderate Resolution Imaging Spectroradiometer (MODIS)-Aqua satellite. Currently, only one other DL model has been used to forecast Sea Surface Salinity (SSS), being a Recurrent Neural Network (RNN). However, the RNN model requires the prediction of a full area of salinity to function. As such, several model architectures were tested, with the best one, being a Multi-input MPNN, utilized to evaluate the feasibility of forecasting utilizing simpler DL models. The results show that a one-day forecast is plausible, while three and five-day forecasts would require a data-rich environment, unlike that of the Laguna Madre

Spatial predictions on physically constrained domains: Applications to Arctic sea salinity data

In this paper, we predict sea surface salinity (SSS) in the Arctic Ocean based on satellite measurements. SSS is a crucial indicator for ongoing changes in the Arctic Ocean and can offer important insights about climate change. We particularly focus on areas of water mistakenly flagged as ice by satellite algorithms. To remove bias in the retrieval of salinity near sea ice, the algorithms use conservative ice masks, which result in considerable loss of data. We aim to produce realistic SSS values for such regions to obtain more complete understanding about the SSS surface over the Arctic Ocean and benefit future applications that may require SSS measurements near edges of sea ice or coasts. We propose a class of scalable nonstationary processes that can handle large data from satellite products and complex geometries of the Arctic Ocean. Barrier Overlap-Removal Acyclic directed graph GP (BORA-GP) constructs sparse directed acyclic graphs (DAGs) with neighbors conforming to barriers and boundaries, enabling characterization of dependence in constrained domains. The BORA-GP models produce more sensible SSS values in regions without satellite measurements and show improved performance in various constrained domains in simulation studies compared to state-of-the-art alternatives. The R package is available on https://github.com/jinbora0720/boraGP

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