Ocean Salinity Research Articles

How Can Seasonality Influence the Performance of Recent Microwave Satellite Soil Moisture Products?

In addition to technical issues related to the instruments used, differences between soil moisture (SM) measured using ground-based methods and microwave remote sensing (RS) can be related to the main features of the study areas, which are intricately connected to hydraulic–hydrological conditions and soil properties. When long-term analysis is performed, these discrepancies are mitigated by the contribution of SM seasonality and are only evident when high-frequency variations (i.e., signal anomalies) are investigated. This study sought to examine the responsiveness of SM to seasonal variations in terrestrial ecoregions located in areas covered by the in situ Romanian Soil Moisture Network (RSMN). To achieve this aim, several remote sensing-derived retrievals were considered: (i) NASA’s Soil Moisture Active and Passive (SMAP) L4 V5 model assimilated product data; (ii) the European Space Agency’s Soil Moisture and Ocean Salinity INRA–CESBIO (SMOS-IC) V2.0 data; (iii) time-series data extracted from the H115 and H116 SM products, which are derived from the analysis of Advanced Scatterometer (ASCAT) data acquired via MetOp satellites; (iv) Copernicus Global Land Service SSM 1 km data; and (v) the “combined” European Space Agency’s Climate Change Initiative for Soil Moisture (ESA CCI SM) product v06.1. An initial assessment of the performance of these products was conducted by checking the anomaly of long-term fluctuations, quantified using the Absolute Variation of Local Change of Environment (ALICE) index, within a time frame spanning 2015 to 2020. These correlations were then compared with those based on raw data and anomalies computed using a moving window of 35 days. Prominent correlations were observed with the SMAP L4 dataset and across all ecoregions, and the Balkan mixed forests (646) exhibited strong concordance regardless of the satellite source (with a correlation coefficient RALICE > 0.5). In contrast, neither the Central European mixed forests (No. 654) nor the Pontic steppe (No. 735) were adequately characterized by any satellite dataset (RALICE < 0.5). Subsequently, the phenological seasonality and dynamic behavior of SM were computed to investigate the effects of the wetting and drying processes. Notably, the Central European mixed forests (654) underwent an extended dry phase (with an extremely low p-value of 2.20 × 10−16) during both the growth and dormancy phases. This finding explains why the RSMN showcases divergent behavior and underscores why no satellite dataset can effectively capture the complexities of the ecoregions covered by this in situ SM network.

Real-Time Optical Fiber Salinity Interrogator Based on Time-Domain Demodulation and TPMF Incorporated Sagnac Interferometer.

A novel demodulation scheme for a point-type fiber sensor is designed for salinity concentration monitoring based on a Sagnac interferometer (SI) composed of a tapered polarization-maintaining fiber (TPMF) and optical time stretching technology. The SI, constructed using a PMF with a taper region of 5.92 μm and an overall length of 30 cm, demonstrated a notable enhancement in the evanescent field, which intensifies the interaction between the light field and external salinity. This enhancement allows for a direct assessment of salinity concentration changes by analyzing the variations in the SI reflection spectra and the experimental results indicate that the sensitivity of the sensor is 0.151 nm/‱. In contrast to traditional fiber optic sensors that depend on spectral demodulation with slower response rates, this work introduces a new approach where the spectral shift is translated to the time domain, utilizing a dispersion compensation fiber (DCF) with the demodulation rate reaching up to 50 MHz. The experimental outcomes reveal that the sensor exhibits a sensitivity of -0.15 ns/‱ in the time domain. The designed sensor is anticipated to play a pivotal role in remote, real-time monitoring of ocean salinity.

Generation of state‐dependent ensemble perturbations based on time‐varying seawater density for GloSea5 initialization

AbstractIn this study, we developed a flow‐dependent oceanic initialization system for initializing the oceanic temperature and salinity in the Global Seasonal forecast system version 5 (GloSea5). Our algorithm overcomes the limitation of stationary perturbations for Ensemble Optimal Interpolation (EnOI) by spreading observed information along isopycnal lines to create three‐dimensional snapshot density states. The proposed algorithm, which we call state‐dependent ensemble‐based EnOI (SD‐EnOI), takes into account changes in the background error covariance over time without relying on ensemble model simulations. To evaluate the quality of the oceanic initial conditions (ICs) produced by SD‐EnOI, we compared them with those generated by the Global Ocean Data Assimilation and Prediction System version 1 (GODAPS1) operated by the Korea Meteorological Agency (KMA) throughout January 2017 to December 2017. Our findings show that the thermal construction of the SD‐EnOI ICs is more realistic than that of GODAPS1, particularly in the tropical Pacific region. The strong warm bias in sea surface temperature (SST) and the shallow mixed‐layer depth bias observed in the GODAPS1 ICs are not shown in SD‐EnOI. Due to the more realistic oceanic thermal structure present in the SD‐EnOI ICs, their use in retrospective forecast experiments resulted in a systematic reduction in climatological SST drift in the central‐eastern Pacific for forecasts up to four lead months compared to using GODAPS1 ICs. This demonstrates the significant impact of the initialization process on the quality of dynamical seasonal forecasts.

A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data

High-resolution three-dimensional (3D) variations in ocean temperature and salinity fields are of great significance for ocean environment monitoring. Currently, AI-based 3D temperature and salinity field predictions rely on expensive 3D data, and as the prediction period increases, the stacking of high-resolution 3D data greatly increases the difficulty of model training. This paper transforms the prediction of 3D temperature and salinity into the prediction of sea surface elements and the inversion of subsurface temperature and salinity using sea surface elements, by leveraging the relationship between sea surface factors and subsurface temperature and salinity. This method comprehensively utilizes multi-source ocean data to avoid the issue of data volume caused by stacking high-resolution historical data. Specifically, the model first utilizes 1/4° low-resolution satellite remote sensing data to construct prediction models for sea surface temperature (SST) and sea level anomaly (SLA), and then uses 1/12° high-resolution temperature and salinity data as labels to build an inversion model of subsurface temperature and salinity based on SST and SLA. The prediction model and inversion model are integrated to obtain the final high-resolution 3D temperature and salinity prediction model. Experimental results show that the 20-day prediction results in the two sea areas of the coastal waters of China and the Northwest Pacific show good performance, accurately predicting ocean temperature and salinity in the vast majority of layers, and demonstrate higher resource utilization efficiency.

P-band radiometry for enhanced vegetation optical depth (VOD) and soil moisture retrieval in dense crop canopies

P-band (0.3–1 GHz) radiometry shows promise in enhancing the next generation of radiometer missions for global soil moisture and vegetation optical depth (VOD) monitoring, particularly in moderately to densely vegetated areas, which remain challenging for current L-band (1.4 GHz) missions such as Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP). However, the superiority of P-band over L-band in this regard has not been experimentally verified, and the quantitative differences remain unexplored. To address this gap, the P-band Radiometer Inferred Soil Moisture (PRISM) tower-based experiment was conducted, pioneering a comprehensive evaluation of P- (0.75 GHz) and L-band (1.4 GHz) radiometry in densely vegetated areas with vegetation water content (VWC) up to 20 kg/m2. The results revealed that P-band VOD was consistently lower than L-band VOD, suggesting reduced vegetation attenuation and thus a higher VWC threshold for reliable soil moisture sensing at P-band. The VWC thresholds for P- and L-band were estimated at 7.7 kg/m2 and 4.0 kg/m2, respectively, suggesting that P-band radiation could go through the vegetation approximately twice as effectively as L-band radiation. Additionally, the Microwave Polarization Difference Index (MPDI) was identified as a strong predictor of retrieval performance because it established a good link between retrieval errors and VWC. Regarding VOD retrieval, robust logarithmic relationships (R2 > 0.9) between VOD and VWC were established for both bands within the 0–20 kg/m2 VWC range, challenging the previously assumed linear relationship. Moreover, P-band demonstrated greater retrieval robustness as its dual-polarized TB observations contain more independent information than those of L-band. This research establishes the foundation for a successful passive microwave mission at P-band, aiming to improve the quality of global VOD and soil moisture datasets.

Assimilating multivariate remote sensing data into a fully coupled subsurface-land surface hydrological model

Hydrological models play a crucial role in tracking and predicting terrestrial water storage, yet they face challenges due to uncertainties and inaccuracies caused by various factors such as meteorological processes and data limitations. To refine these models, data assimilation has emerged as a valuable tool, utilizing a new source of data to update model states while considering associated uncertainties, thereby enhancing our comprehension and predictive capabilities in hydrological processes. In this context, satellite data are receiving increasing attention because they can cover large areas and are useful in detecting spatial and temporal variability of water.In most existing studies related to satellite data assimilation in hydrological models, one source satellite data is used in the analysis. This study focuses on improving subsurface water storage model accuracy by assimilating data from different satellite sources. In particular, we used data from the Soil Moisture and Ocean Salinity (SMOS) satellite and terrestrial water storage data from the Gravity Recovery and Climate Experiment (GRACE). The data are assimilated into a fully coupled subsurface-surface hydrological model, developed with ParFlow-CLM (PARallel FLOW- Community Land Model). The investigation is conducted in Iran. Employing an Ensemble Kalman Filter, three assimilation scenarios are explored: (i) GRACE, (ii) SMOS, and (iii) the combined assimilation of both GRACE and SMOS data (joint). Findings are validated against the Soil Moisture Active Passive (SMAP) and in-situ groundwater data by using a novel probabilistic reliability framework, demonstrating the advantages of joint data assimilation.The study highlights the influence of assimilated remote sensing data type on the effectiveness of data assimilation. Assimilating GRACE data enhances groundwater level estimations. However, SMOS data positively impacts topsoil moisture estimations, but adversely affects groundwater level estimates. Importantly, the assimilation of both GRACE and SMOS data through multivariate (joint) data assimilation significantly improves accuracy for both soil moisture estimation and groundwater level estimation.

Impact of sampling depth on CO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} flux estimates

The exchange of trace gases between the atmosphere and the ocean plays a key role in the Earth’s climate. Fluxes at the air-sea interface are affected mainly by wind blowing over the ocean and seawater temperature and salinity changes. This study aimed to quantify the use of CO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} partial pressure (pCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}) measurements at different depths (1, 5, and 10 m) in ocean surface layers to determine CO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} fluxes (FCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}) and to investigate the impacts of wind-sheltered and wind-exposed regions on the carbon budget. Vertical profiles of temperature, salinity, and pCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} were considered during a daily cycle. pCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} profiles exhibited relatively high values during sunny hours, associated with relatively high sea temperatures. However, the largest FCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} corresponded with higher wind speeds. Estimated fluxes between measurements at 1 and 10 m depths decreased by 71% in the sheltered region and 44% in the exposed region. According to the SOCAT dataset, at a depth of 5 m, the Atlantic basin emits approximately 0.29 Tg month-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document} of CO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} to the atmosphere; nevertheless, our estimates suggest that FCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document} at the surface is 12.02 Tg month-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document}, which is 97.6% greater than that at 5 m depth. Therefore, future studies should consider sampling depth to adequately estimate the FCO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}.

Evaluating soil moisture retrieval in Arctic and sub-Arctic environments using passive microwave satellite data

ABSTRACT Soil Moisture (SM) is a key parameter in northern Arctic and sub-Arctic (A-SA) environments that are highly vulnerable to climate change. We evaluated six SM satellite passive microwave datasets using thirteen ground-based SM stations across Northwestern America. The best agreement was obtained with SMAP (Soil Moisture Active Passive) products with the lowest RMSD (Root Mean Square Difference) (0.07 m 3 m − 3 ) and the highest R (0.55). ESA CCI (European Space Agency Climate Change Initiative) also performed well in terms of correlation with a similar R (0.55) but showed a strong variation among sites. Weak results were obtained over sites with high water body fractions. This study also details and evaluates a dedicated retrieval of SM from SMOS (Soil Moisture and Ocean Salinity) brightness temperatures based on the τ − ω model. Two soil dielectric models (Mironov and Bircher) and a dedicated soil roughness and single scattering albedo parameterization were tested. Water body correction in the retrieval shows limited improvement. The metrics of our retrievals (RMSD = 0.08 m 3 m − 3 and R = 0.41) are better than SMOS but outperformed by SMAP. Passive microwave satellite remote sensing is suitable for SM retrieval in the A-SA region, but a dedicated approach should be considered.

Climatic Effects of Ocean Salinity on M Dwarf Exoplanets

Ocean salinity is known to dramatically affect the climates of Earth-like planets orbiting Sun-like stars, with high salinity leading to less ice and higher surface temperature. However, how ocean composition impacts climate under different conditions, such as around different types of stars or at different positions within the habitable zone, has not been investigated. We used ROCKE-3D, an ocean-atmosphere general circulation model, to simulate how planetary climate responds to ocean salinities for planets with G-star versus M dwarf hosts at several stellar fluxes. We find that increasing ocean salinity from 20 to 100 g kg−1 in our model results in nonlinear ice reduction and warming on G-star planets, sometimes causing abrupt transitions to different climate states. Conversely, sea ice on M dwarf planets responds more gradually and linearly to increasing salinity. Moreover, reductions in sea ice on M dwarf planets are not accompanied by significant surface warming as on G-star planets. High salinity can modestly bolster the resilience of M dwarf planets against snowball glaciation and allow these planets to retain surface liquid water further from their host star, but the effects are muted compared to G-star planets that experience snowball bifurcation and climate hysteresis due to the ice-albedo feedback.

Controls on Upper Ocean Salinity Variability in the Eastern Subpolar North Atlantic During 1992–2017

AbstractThe upper ocean salinity in the eastern subpolar North Atlantic undergoes decadal fluctuations. A large fresh anomaly event occurred during 2012–2016. Using the ECCOv4r4 state estimate, we diagnose and compare mechanisms of this low salinity event with those of the 1990s fresh anomaly event. To avoid issues related to the choice of reference salinity values in the freshwater budget, we perform a salt mass content budget analysis of the eastern subpolar North Atlantic. It shows that the recent low salt content anomaly occurs due to the circulation of anomalous salinity by mean currents entering the eastern subpolar basin from its western boundary via the North Atlantic Current. This is in contrast to the early 1990s, when the dominant mechanism governing the low salt content anomaly was the transport of the mean salinity field by anomalous currents.

Seasonal-scale intercomparison of SMAP and fused SMOS-SMAP soil moisture products

Two L-band passive microwave satellite sensors, onboard the Soil Moisture and Ocean Salinity (SMOS) launched in 2009 and Soil Moisture Active Passive (SMAP) launched in 2015, are specifically designed for surface soil moisture (SM) monitoring. The first global continuous fused L-band satellite SM product based on SMOS and SMAP observations (SMOS-SMAP-INRAE-BORDEAUX, the so-called Fused-IB) was recently released to the public. Currently, the performance of Fused-IB has only been evaluated collectively over the entire data records in the study period, without specific evaluation for individual seasons. To fill this gap, this study intercompared the Fused-IB and the enhanced SMAP-L3 version 6 (SMAP-E) SM products against in situ SM data from the International Soil Moisture Network (ISMN) from 2016 to 2020 regarding the whole period and different seasons. We aim to investigate the performance of these two products at different time scales and to explore the potential eco-hydrological factors (i.e., precipitation and vegetation) driving their seasonal variations. Results show that both SM products are in good agreement with the in situ measurements, demonstrating high median correlation (R) and low ubRMSD (median R = 0.70 and ubRMSD = 0.058 m3/m3 for Fused-IB vs. R = 0.68 and ubRMSD = 0.059 m3/m3 for SMAP-E) during 2016–2020. For most land use and land cover (LULC) types, Fused-IB outperformed SMAP-E with higher accuracy and lower errors, particularly in forests, partly due to the advantage of the robust SMAP-IB (SMAP-INRAE-BORDEAUX) algorithm used to generate Fused-IB in forests, which avoids the pronounced saturation effects of vegetation optical depth caused by relying on optical information. Besides, both products had superior performances across most LULC types in summer (JJA) and autumn (SON), yet increased uncertainties were observed in forests, grasslands, and croplands during spring (MAM) and winter (DJF). These uncertainties could be mainly attributed to the effects of vegetation growth in forests, grasslands and croplands, and the interception of water from rainfall events in grasslands. The results of this study can serve as a reference for algorithm developers to enhance the accuracy of SM and thus promote hydro-meteorological applications that benefit from L-band radiometer soil moisture products.

Ocean exploration-oriented temperature and salinity sensor based on bend-insensitive microfiber Mach-Zehnder interferometer.

Ocean exploration-oriented temperature and salinity (TS) sensor based on bend-insensitive microfiber Mach-Zehnder interferometer (MMZI) is proposed and demonstrated in a marine environment. To solve the demodulation problem induced by the narrow waveband of the spectrometer used in the sea trial, a MMZI sensor is calibrated and demodulated by machine learning method. Results show that even if the wavelength range used in demodulation is as small as tens of nanometers, a relatively accurate demodulation can still be achieved. Then a sea trial is performed based on the design of the test system on board. Compared with commercial conductivity-temperature-depth systems, mean relative errors of 2.28% and 2.29% for TS measurement are obtained. In addition, the influence of spectral range on the demodulation result and response of the spectrum to fiber bending is discussed. Repeatability and stability of the sensor are proved in repeated sea trials. The sensor demonstrated here provides a new optical method for measuring TS in seawater with the advantages of multi-parameter measurement, low cost, easy deployment, and high accuracy, which will hopefully become a useful complement to the existing electrical method used in ocean exploration.

A Robust Sparse Sensor Placement Strategy Based on Indicators of Noise for Ocean Monitoring

A well-performing data-driven sparse sensor deployment strategy is critical for marine monitoring systems, as it enables the optimal reconstruction of marine physical quantities with fewer sensors. However, ocean data typically contain substantial amounts of noise, including outliers (incomplete data) and inherent measurement noise, which heightens the complexity of sensor deployment. Therefore, this study optimizes the sparse sensor placement model by establishing noise indicators, including small noise weight and large noise weight, which are measured by entropy to minimize the prediction bias. Building on this, a robust sparse sensor placement algorithm is proposed, which utilizes the block coordinate update (BCU) iteration method to solve the function. During the iterative updating process, the proposed algorithm simultaneously updates the selection matrix, reconstruction matrix, and noise matrix. This allows for effective identification and mitigation of noise in the data through evaluation. Consequently, the deployed sensors achieve superior reconstruction performance compared to other deployment methods that do not incorporate noise evaluation. Experiments are also conducted on datasets of sea surface temperature (SST) and global ocean salinity, which demonstrate that our strategy significantly outperforms several other considered methods in terms of reconstruction accuracy while enabling autonomous sensor deployment under noisy conditions.

Influence of multi-stressor combinations of pCO2, temperature, and salinity on the toxicity of Heterosigma akashiwo (Raphidophyceae), a fish-killing flagellate.

Climate change and global warming have led to more frequent harmful algal blooms in the last decade. Among these blooms, Heterosigma akashiwo, a golden-brown phytoflagellate, is one of the 40 species with a high potential to form harmful blooms, leading to significant fish mortality. Climate change leads to rising atmospheric and ocean temperatures. These changes, along with altered rainfall patterns and meltwater input, can cause fluctuations in ocean salinity. Elevated atmospheric carbon dioxide (CO2) levels increase water acidity as oceans absorb CO2. This study investigated the effects of temperature, salinity, and CO2 levels on lipid production, hemolytic activity, and toxicity of H. akashiwo using the design of experiment approach, which can be used to investigate the effect of two or more factors on the same response simultaneously in a precise manner with fewer experiments and materials but in a larger region of the factor space. The lipid content was measured using a high-throughput Nile Red method, and the highest level of lipid content was detected at 25°C, a salinity of 30, and a CO2 concentration of 400 ppm. Hemolytic activity was assessed using rabbit blood erythrocytes in a 96-well plate, and the optimal conditions for achieving the highest hemolytic activity were determined at 15°C, a salinity of 10, and a CO2 concentration of 400 ppm. As the chemical structure of the toxin is not known, we used the toxicity against the cell line RTgill-W1 as the cell toxicity proxy. The maximum toxicity was identified at 15°C, a salinity of 10, and a CO2 level of 700 ppm.

A machine learning-based framework for spatio-temporal extension and filling of SMOS surface soil moisture observations over Canada

ABSTRACT Advancements in global satellite missions have revolutionized the assessment of Surface Soil Moisture (SSM) at global to local scales. However, spatio-temporal data discontinuities in specific regions remain a challenge. This study proposes a Machine Learning (ML)-based framework to extend the Soil Moisture Ocean Salinity (SMOS) SSM product both in the spatial and temporal domains, over Canada. In the first phase of the proposed framework, ML models based on Random Forest (RF) and Convolutional Neural Networks (CNN) are trained and validated with SMOS SSM as target and SSM-relevant climatic variables and geophysical variables, obtained from fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data (ERA5), for the 2011–2020 period, as predictors. Developed models, when tested on unseen data for the years 2021–2022, suggest slightly better performance of the RF model compared to CNN, with root mean square error (RMSE) of 0.033 and 0.056 respectively; prediction biases mostly noted for regions with large inter-annual variability. The spatial filling of SSM for grid cells that were excluded during the training process, with similar land types as those in the SMOS training data, yields reasonable performance, with RF (RMSE = 0.013) performing better than CNN (RMSE = 0.064). In the second phase, the RF model is selected to extend the SMOS dataset for the 2008–2010 period. The temporal correlation between extended SMOS and ASCAT (Advanced Scatterometer) SSM demonstrates a reasonable association, with correlation coefficient exceeding 0.6. Additionally, spatial correlation analysis reveals similar patterns between the two datasets, with smaller values for the summer season owing to the importance of local processes on SSM during this period. However, extending SMOS SSM spatially for surface types that were not included in the training process such as peatlands, remains a challenge warranting additional studies. The developed framework is robust and can address spatio-temporal discontinuities in other SSM products.

Estimating the uncertainties of satellite derived soil moisture at global scale

This study attempts to derive the uncertainty of the soil moisture estimation from passive microwave satellite mission at global scale. To do so, the approach is based on the sensitivity of the Soil Moisture and Ocean Salinity (SMOS) soil moisture retrieval quality to the land surface characteristics within its footprint (presence of forest, topography, open water bodies, sand, clay, bulk density and soil organic carbon content). First, we performed a global assessment of SMOS using in situ measurements from the International Soil Moisture Network (ISMN) as reference, with more than 1900 ISMN stations and 10 years of SMOS data. This assessment shows that the ubRMSD scores vary greatly between locations (with a mean of 0.074 m3m−3 and an interquartile range of 0.030 m3m−3). Second, the scores are analyzed for different surface conditions within the satellite footprint. The best agreement between the ground measurement and SMOS time series are obtained for low forest cover, low topographic complexity, and marginal presence of open water bodies within the SMOS footprint. Soil parameters also have an impact, with better scores for sandier soils with a high bulk-density and low soil organic carbon content. Finally, we propose to extrapolate the obtained relationships, using a multiple linear regression, in order to derive a global map of SMOS uncertainties based on surface conditions. This map of predicted uncertainties show a diverse range of ubRMSD values across the globe (with a mean of 0.076 m3m−3 and an interquartile range of 0.031 m3m−3) depending on the surface characteristics. At the ISMN site location, the predicted ubRMSD shows similar results than the comparison between SMOS and the in situ measurements. The map of predicted SMOS ubRMSD represents an upper bound estimate of the SMOS uncertainty, as it includes the uncertainties of the in situ sensor measurements and the scale mismatch. Further investigations will focus on the different components of this uncertainty budget to obtain a better assessment of the absolute uncertainties of SMOS soil moisture retrievals across the globe.

SMOS-derived Antarctic thin sea ice thickness: data description and validation in the Weddell Sea

Abstract. Accurate satellite measurements of the thickness of Antarctic sea ice are urgently needed but pose a particular challenge. The Antarctic data presented here were produced using a method to derive the sea ice thickness from 1.4 GHz brightness temperatures previously developed for the Arctic, with only modified auxiliary data. The ability to observe the thickness of thin sea ice using this method is limited to cold conditions, meaning it is only reasonable during the freezing period, typically March to October. The Soil Moisture and Ocean Salinity (SMOS) level-3 sea ice thickness product contains estimates of the sea ice thickness and its uncertainty up to a thickness of about 1 m. The sea ice thickness is provided as a daily average on a polar stereographic projection grid with a sample resolution of 12.5 km, while the SMOS brightness temperature data used have a footprint size of about 35–40 km in diameter. Data from SMOS have been available since 2010, and the mission's operation has been extended to continue until at least the end of 2025. Here we compare two versions of the SMOS Antarctic sea ice thickness product which are based on different level-1 input data (v3.2 based on SMOS L1C v620 and v3.3 based on SMOS L1C 724). A validation is performed to generate a first baseline reference for future improvements of the retrieval algorithm and synergies with other sensors. Sea ice thickness measurements to validate the SMOS product are particularly rare in Antarctica, especially during the winter season and for the valid range of thicknesses. From the available validation measurements, we selected datasets from the Weddell Sea that have varying degrees of representativeness: Helicopter-based EM Bird (HEM), Surface and Under-Ice Trawl (SUIT), and stationary Upward-Looking Sonars (ULS). While the helicopter can measure hundreds of kilometres, SUIT's use is limited to distances of a few kilometres and thus only captures a small fraction of an SMOS footprint. Compared to SMOS, the ULS are point measurements and multi-year time series are necessary to enable a statistically representative comparison. Only four of the ULS moorings have a temporal overlap with SMOS in the year 2010. Based on selected averaged HEM flights and monthly ULS climatologies, we find a small mean difference (bias) of less than 10 cm and a root mean square deviation of about 20 cm with a correlation coefficient R > 0.9 for the valid sea ice thickness range between 0 and about 1 m. The SMOS sea ice thickness showed an underestimate of about 40 cm with respect to the less representative SUIT validation data in the marginal ice zone. Compared with sea ice thickness outside the valid range, we find that SMOS strongly underestimates the real values, which underlines the need for combination with other sensors such as altimeters. In summary, the overall validity of the SMOS sea ice thickness for thin sea ice up to a thickness of about 1 m has been demonstrated through validation with multiple datasets. To ensure the quality of the SMOS product, an independent regional sea ice extent index was used for control. We found that the new version, v3.3, is slightly improved in terms of completeness, indicating fewer missing data. However, it is worth noting that the general characteristics of both datasets are very similar, also with the same limitations. Archived data are available in the PANGAEA repository at https://doi.org/10.1594/PANGAEA.934732 (Tian-Kunze and Kaleschke, 2021) and operationally at https://doi.org/10.57780/sm1-5ebe10b (European Space Agency, 2023).

Convformer: A Model for Reconstructing Ocean Subsurface Temperature and Salinity Fields Based on Multi-Source Remote Sensing Observations

Observational data on ocean subsurface temperature and salinity are patently insufficient because in situ observations are complex and costly, while satellite remote-sensed measurements are abundant but mainly focus on sea surface data. To make up for the ocean interior data shortage and entirely use the abundant satellite data, we developed a data-driven deep learning model named Convformer to reconstruct ocean subsurface temperature and salinity fields from satellite-observed sea surface data. Convformer is designed by deeply optimizing Vision Transformer and ConvLSTM, consisting of alternating residual connections between multiple temporal and spatial attention blocks. The input variables consist of sea surface temperature (SST), sea surface salinity (SSS), sea surface height (SSH), and sea surface wind (SSW). Our results demonstrate that Convformer exhibits superior performance in estimating the temperature-salinity structure of the tropical Pacific Ocean. The all-depth average root mean square error (RMSE) of the reconstructed subsurface temperature (ST)/subsurface salinity (SS) is 0.353 °C/0.0695 PSU, with correlation coefficients (R²) of 0.98663/0.99971. In the critical thermocline, although the root mean square errors of ST and SS reach 0.85 °C and 0.121 PSU, respectively, they remain smaller compared to other models. Furthermore, we assessed Convformer’s performance from various perspectives. Notably, we also delved into the potential of Convformer to extract physical and dynamic information from a model mechanism perspective. Our study offers a practical approach to reconstructing the subsurface temperature and salinity fields from satellite-observed sea surface data.

DSE-NN: Discretized Spatial Encoding Neural Network for Ocean Temperature and Salinity Interpolation in the North Atlantic

The precise interpolation of oceanic temperature and salinity is crucial for comprehending the dynamics of marine systems and the implications of global climate change. Prior neural network-based interpolation methods face constraints related to their capacity to delineate the intricate spatio-temporal patterns that are intrinsic to ocean data. This research presents an innovative approach, known as the Discretized Spatial Encoding Neural Network (DSE-NN), comprising an encoder–decoder model designed on the basis of deep supervision, network visualization, and hyperparameter optimization. Through the discretization of input latitude and longitude data into specialized vectors, the DSE-NN adeptly captures temporal trends and augments the precision of reconstruction, concurrently addressing the complexity and fragmentation characteristic of oceanic data sets. Employing the North Atlantic as a case study, this investigation shows that the DSE-NN presents enhanced interpolation accuracy in comparison with a traditional neural network. The outcomes demonstrate its quicker convergence and lower loss function values, as well as the ability of the model to reflect the spatial and temporal distribution characteristics and physical laws of temperature and salinity. This research emphasizes the potential of the DSE-NN in providing a robust tool for three-dimensional ocean temperature and salinity reconstruction.

Brachiopods of the northern North Atlantic and Arctic, with a focus on Norwegian fauna

Ten species of Brachiopoda have been found living around Norway, including the Arctic Svalbard and Jan Mayen, while 26 are recorded here for the region extending from the Arctic Ocean in the North to the northernmost North Atlantic and Celtic Sea in the South. This paper provides an identification key and short descriptions for all species, including the new species Xenobrochus islandicus n. sp. from southwest of Iceland and the new East Atlantic species Dallina lusitanica n. sp. The distributions of all species have been reviewed and were generally found to follow the oceanic temperature, salinity and depth gradients of the region. Also organic enrichment had some influence on distribution. Some species have shown a possible ongoing northward shift of their biogeographic boundaries over the last two centuries, including the Arctic to subarctic species Hemithiris psittacea.