The present study assesses the performance of state-of-the-art atmospheric correction (AC) algorithms applied to Sentinel-2-MultiSpectral Instrument (S2-MSI) and Sentinel-3-Ocean and Land Color Instrument (S3-OLCI) data recorded over moderately to highly turbid estuarine waters, considering the Gironde Estuary (SW France) as a test site. Three spectral bands of water-leaving reflectance ( $Rhow$ ) are considered: green (560 nm), red (655 or 665 nm) and near infrared (NIR) (865 nm), required to retrieve the suspended particulate matter (SPM) concentrations in clear to highly turbid waters (SPM ranging from 1 to 2000 mg/L). A previous study satisfactorily validated Acolite short wave infrared (SWIR) AC algorithm for Landsat-8-Operational Land Imager (L8-OLI) in turbid estuarine waters. The latest version of Acolite Dark Spectrum Fitting (DSF) is tested here and shows very good agreement with Acolite SWIR for OLI data. L8-OLI satellite data corrected for atmospheric effects using Acolite DSF are then used as a reference to assess the validity of atmospheric corrections applied to other satellite data recorded over the same test site with a minimum time difference. Acolite DSF and iCOR (image correction for atmospheric effects) are identified as the best performing AC algorithms among the tested AC algorithms (Acolite DSF, iCOR, Polymer and C2RCC (case 2 regional coast color)) for S2-MSI. Then, the validity of six different AC algorithms (OLCI Baseline Atmospheric Correction (BAC), iCOR, Polymer, Baseline residual (BLR), C2RCC-V1 and C2RCC-V2) applied to OLCI satellite data is assessed based on comparisons with OLI and/or MSI Acolite DSF products recorded on a same day with a minimum time lag. Results show that all the AC algorithms tend to underestimate $Rhow$ in green, red and NIR bands except iCOR in green and red bands. The iCOR provides minimum differences in green (slope = 1.0 ± 0.15, BIAS = 1.9 ± 4.5% and mean absolute percentage error (MAPE) = 12 ± 5%) and red (slope = 1.0 ± 0.17, BIAS = −9.8 ± 9% and MAPE = 28 ± 20%) bands with Acolite DSF products from OLI and MSI data. For the NIR band, BAC provides minimum differences (slope = 0.7 ± 0.13, BIAS = −33 ± 17% and MAPE = 55 ± 20%) with Acolite DSF products from OLI and MSI data. These results based on comparisons between almost simultaneous satellite products are supported by match-ups between satellite-derived and field-measured SPM concentrations provided by automated turbidity stations. Further validation of satellite products based on rigorous match-ups with in-situ $Rhow$ measurements is still required in highly turbid waters.

Because of the large diversity of case 2 waters ranging from extremely absorbing to extremely scattering waters and the complexity of light transfer due to external terrestrial inputs, retrieving main biogeochemical parameters such as chlorophyll-a or suspended particulate matter concentration in these waters is still challenging. By providing optical and biogeochemical parameters for 180 sampling stations with turbidity and chlorophyll-a concentration ranging from 1 to 700 FNU and from 0.9 to 180 mg m⁻³ respectively, the HYPERMAQ dataset will contribute to a better description of marine optics in optically complex water bodies and can help the scientific community to develop algorithms. The HYPERMAQ dataset provides biogeochemical parameters (i.e. turbidity, pigment and chlorophyll-a concentration, suspended particulate matter), apparent optical properties (i.e. water reflectance from above water measurements) and inherent optical properties (i.e. absorption and attenuation coefficients) from six different study areas. These study areas include large estuaries (i.e. the Rio de la Plata in Argentina, the Yangtze estuary in China, and the Gironde estuary in France), inland (i.e. the Spuikom in Belgium and Chascomùs lake in Argentina), and coastal waters (Belgium). The dataset is available from Lavigne et al. (2022) at https://doi.org/10.1594/PANGAEA.944313.

Hydro-sedimentary numerical models have been widely employed to derive suspended particulate matter (SPM) concentrations in coastal and estuarine waters. These hydro-sedimentary models are computationally and technically expensive in nature. Here we have used a computationally less-expensive, well-established methodology of self-organizing maps (SOMs) along with a hidden Markov model (HMM) to derive profiles of suspended particulate inorganic matter (SPIM). The concept of the proposed work is to benefit from all available data sets through the use of fusion methods and machine learning approaches that are able to process a growing amount of available data. This approach is applied to two different data sets entitled “Hidden” and “Observable”. The hidden data are composed of 15 months (27 September 2007 to 30 December 2008) of hourly SPIM profiles extracted from the Regional Ocean Modeling System (ROMS). The observable data include forcing parameter variables such as significant wave heights ( $Hs$ and $Hs50$ (50 days)) from the Wavewatch 3-HOMERE database and barotropic currents ( $Ubar$ and $Vbar$ ) from the Iberian–Biscay–Irish (IBI) reanalysis data. These observable data integrate hourly surface samples from 1 February 2002 to 31 December 2012. The time-series profiles of the SPIM have been derived from four different stations in the English Channel by considering 15 months of output hidden data from the ROMS as a statistical representation of the ocean for ≈11 years. The derived SPIM profiles clearly show seasonal and tidal fluctuations in accordance with the parent numerical model output. The surface SPIM concentrations of the derived model have been validated with satellite remote sensing data. The time series of the modeled SPIM and satellite-derived SPIM show similar seasonal fluctuations. The ranges of concentrations for the four stations are also in good agreement with the corresponding satellite data. The high accuracy of the estimated 25 h average surface SPIM concentrations (normalized root-mean-square error— $NRMSE$ of less than 16%) is the first step in demonstrating the robustness of the method.

Satellite remote sensing observations allow the ocean surface to be sampled synoptically over large spatio-temporal scales. The images provided from visible and thermal infrared satellite observations are widely used in physical, biological, and ecological oceanography. The present work proposes a method to understand the multi-scaling properties of satellite products such as the Chlorophyll-a (Chl-a), and the Sea Surface Temperature (SST), rarely studied. The specific objectives of this study are to show how the small scale heterogeneities of satellite images can be characterised using tools borrowed from the fields of turbulence. For that purpose, we show how the structure function, which is classically used in the frame of scaling time series analysis, can be used also in 2D. The main advantage of this method is that it can be applied to process images which have missing data. Based on both simulated and real images, we demonstrate that coarse-graining (CG) of a gradient modulus transform of the original image does not provide correct scaling exponents. We show, using a fractional Brownian simulation in 2D, that the structure function (SF) can be used with randomly sampled couple of points, and verify that 1 million of couple of points provides enough statistics.