The main objective of the SWIM Level 3 (L3) product is to provide global gridded wave statistics from the along-track directional wave spectra measures, as observed with SWIM instrument onboard CFOSAT. CFOSAT (Chinese French Ocean SATellite) is a french-chinese mission launched in 2018, whose aim is to provide wind (SCAT instrument) and wave (SWIM instrument) measurements over the sea surface. The processing uses the integrated parameters (wave length, significant wave height, direction) of the wave systems associated with the individual partitions retrieved from SWIM wave spectra, computing monthly statistics such as distribution, min, max, mean and median over a global 2 degree resolution grid. Three sub-products are generated separately for each SWIM beam (from 4° to 10°) from respectively the SWIM Level 2 spectra produced by the CWWIC center (CNES), the SWIM Level 2S spectra produced by the IWWOC center (Ifremer) and forecast spectra produced from WaveWatch3 wave model at Ifremer as an independent source for comparison. The SWIM level 3 dataset is generated in delayed mode. It is currently mainly meant as a validation dataset aiming at inter-comparing the different SWIM Level 2 between each other, with an independent model (here WaveWatch3) or other external similar wave statistics products. The SWIM L3 dataset is generated and distributed by Ifremer / CERSAT in the frame of the Ifremer Wind and Wave Operation Center (IWWOC) co-funded by Ifremer and CNES and dedicated to the processing of the delayed mode data of CFOSAT mission.

The Level-2 Ka-band Radar Interferometer (KaRIn) low rate (LR, ocean) sea surface height (SSH) data product from the Surface Water and Ocean Topography (SWOT) mission, also referenced by the short name L2_LR_SSH, provides ocean topography measurements from the low rate ocean data stream of the KaRIn instrument, spanning 60 km on either side of the nadir altimeter with a nadir gap. The L2_LR_SSH product is available continuously and globally, although different versions of the product may be produced at different latencies and/or through different reprocessing with refined input data. Note that L2_LR_SSH does not include SSH data from the SWOT nadir altimeter. The SWOT L2_LR_SSH product is organized in four files, the L2_LR_SSH ['Expert'] is described in this metadata sheet. The 3 other file types (['Basic'], ['WindWave'], ['Unsmoothed']) are described by 3 different metadata sheets that can be accessed via the links below. The ['Expert'] file is intended for expert users who are interested in the details of how the KaRIn measurements were derived and who may use detailed information for their own custom processing. The ['Basic'] file is intended for users who are interested in SSH measurements and who will use the KaRIn measurements as provided. The ['WindWave'] file is intended for users interested in wind and wave information. The ['Unsmoothed'] file, also intended for expert users, is provided on a finer 'native' grid of 250-m (with minimal smoothing applied), and has a significantly larger data volume than the other files. The ['Expert'] L2_LR_SSH includes copies of all variables in the Basic and WindWave files (sea surface height (SSH), sea surface height anomaly (SSHA), data quality flags, geophysical reference fields, height correction information, significant wave height (SWH), normalized radar cross section (NRCS or backscatter cross section or sigma0), wind speed derived from sigma0 and SWH, wind and wave model information, quality flags on a 2 km geographically fixed grid), plus more detailed information on the KaRIn instrument and environmental corrections, radiometer data, and geophysical models on a 2 km geographically fixed grid. May 2025: v3.0 (version D) Production and distribution of the L2_LR_SSH version D products: - PID0 for forward-processed version D products: from May 6, 2025 onward, - PGD0 for reprocessed version D products: from March 30 to July 10, 2023 (phase CalVal) and from July 26, 2023 to May 19, 2025 (phase Science) is ongoing. August 2024: v2.0 (version D) L2_LR_SSH version C products declared as validated by the SWOT project. March 2024: v2.0 (version C) Production and distribution of the pre-validated L2_LR_SSH version C products: - PIC0 for forward-processed version C products: November 23, 2023 to present, - PGC0 for reprocessed version C products: from March 30 to July 10, 2023 (phase CalVal) and from July 26, 2023 to January 25, 2024 (phase Science) November 2023: v1.0 The beta pre-validated L2_LR_SSH version 1.0 product (summer 2023 reprocessing release) is available only for the 1-day CalVal orbit phase, from March 29 to July 10, 2023, and the 21-day Science orbit phase from September 7 to November 21, 2023.

Fronts are ubiquitous discrete features of the global ocean often associated with enhanced vertical velocities, in turn boosting primary production and so forth. Fronts thus form dynamical and ephemeral ecosystems where numerous species meet across all trophic levels. Fronts are also targeted by fisheries. Capturing ocean fronts and studying their long-term variability in relation with climate change is thus key for marine resource management and spatial planning. The Mediterranean Sea and the Southwest Indian Ocean are natural laboratories to study front-marine life interactions due to their energetic flow at sub-to-mesoscales, high biodiversity (including endemic and endangered species) and numerous conservation initiatives. Based on remotely-sensed Sea Surface Temperature and Height, we compute thermal fronts (2003-2020) and attracting Lagrangian Coherent Structures (1994-2020), in both regions over several decades. We advocate for the combined use of both thermal fronts and attracting Lagrangian Coherent Structures to study front-marine life interactions. The resulting front database differs from other alternatives by its high spatio-temporal resolution, long time coverage, and relevant thresholds defined for ecological provinces.

'''DEFINITION''' The sea level ocean monitoring indicator has been presented in the Copernicus Ocean State Report #8. The ocean monitoring indicator on mean sea level is derived from the DUACS delayed-time (DT-2024 version, “my” (multi-year) dataset used when available) sea level anomaly maps from satellite altimetry based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and by the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The time series of area averaged anomalies correspond to the area average of the maps in the Northeast Atlantic Ocean and adjacent seas Sea weighted by the cosine of the latitude (to consider the changing area in each grid with latitude) and by the proportion of ocean in each grid (to consider the coastal areas). The time series are corrected from regional mean GIA correction (weighted GIA mean of a 27 ensemble model following Spada et Melini, 2019). The time series are adjusted for seasonal annual and semi-annual signals and low-pass filtered at 6 months. Then, the trends/accelerations are estimated on the time series using ordinary least square fit. Uncertainty is provided in a 90% confidence interval. It is calculated as the weighted mean uncertainties in the region from Prandi et al., 2021. This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation depending on the period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not considered. ""CONTEXT "" Change in mean sea level is an essential indicator of our evolving climate, as it reflects both the thermal expansion of the ocean in response to its warming and the increase in ocean mass due to the melting of ice sheets and glaciers (WCRP Global Sea Level Budget Group, 2018). At regional scale, sea level does not change homogenously. It is influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). The adverse effects of floods, storms and tropical cyclones, and the resulting losses and damage, have increased as a result of rising sea levels, increasing people and infrastructure vulnerability and food security risks, particularly in low-lying areas and island states (IPCC, 2022a). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022b). In this region, sea level variations are influenced by the North Atlantic Oscillation (NAO) (e.g. Delworth and Zeng, 2016) and the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). Hermans et al., 2020 also reported the dominant influence of wind on interannual sea level variability in a large part of this area. This region encompasses the Mediterranean, IBI, North-West shelf, Black Sea and Baltic regions with different sea level dynamics detailed in the regional indicators. ""KEY FINDINGS"" Over the [1999/02/20 to 2025/10/18] period, the area-averaged sea level in the Northeast Atlantic Ocean and adjacent seas area rises at a rate of 4.0 ± 0.8 mm/yr with an acceleration of 0.23 ± 0.06 mm/yr². This trend estimation is based on the altimeter measurements corrected from regional GIA correction (Spada et Melini, 2019) to consider the ongoing movement of land. The TOPEX-A is no longer included in the computation of regional mean sea level parameters (trend and acceleration) with version 2024 products due to potential drifts, and ongoing work aims to develop a new empirical correction. Calculation begins in February 1999 (the start of the TOPEX-B period). '''DOI (product):''' https://doi.org/10.48670/mds-00335

This dataset comprises the global frequency, classification and distribution of marine heat waves (MHWs) from 1996-2020, in Chauhan et al. 2023 (https://doi.org/10.3389/fmars.2023.1177571). The classification was done based on their attributes and using different baselines. Daily SST values were extracted from the NOAA-OISST v2 high-resolution (0.25°) dataset from 1982-2020. MHWs were detected using the method presented by Hobday et al. 2016 (https://doi.org/10.1016/j.pocean.2015.12.014), and by using the 95th percentile of the accumulated temperature distribution to flag the extreme events. A shifting baseline of 8-year rolling period was selected between the years 1982-1996, since this period shows relatively stable maximum values of temperature across different ocean regions. The shifting baseline aims to account for the decadal changes of westerly winds, temperatures and ocean gyres circulations. The classification was done using the KMeans clustering algorithm to identify the relevant features of MHWs and classify them into separate groups based on feature similarities. This algorithm takes MHW features, namely duration, maximum intensity, rate onset and rate decline, as input vectors and applies clustering in the 4-dimensional feature space where each data point represents an MHW event. Note that all the MHWs features are standardized because unequal variances can put more weight on variables with smaller variances. This record comprehends the geospatial datasets of: Average number of MHW days per year (i.e., the sum of all MHW days divided by the total number of years, 1996-2020). Average cumulative intensity per year (i.e., the sum of cumulative intensity divided by the total number of years, 1996-2020). Total number of MHW events across the different periods averaged on the total number of years (1989-2020). The period 1982-1988 was only used as an initial baseline without calculating MHWs. Spatial distribution of three MHW categories: moderate MHWs, abrupt and Intense MHWs and extreme MHWs; displaying the total number of MHW days normalized by the number of years considered (i.e., 1989-2020). Distribution of Extreme MHWs across the different periods (A) 1989-1996, (B) 1997-2004, (C) 2005-2012, (D) 2013-2020. The relative frequency (γ) is a ratio of extreme MHWs in a specific period and all extreme MHWs in the same cluster for all periods.

The anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) populations in the Bay of Biscay are jointly surveyed each year in May since 2000 and in September since 2003 by means of acoustic surveys. The integrated survey PELGAS (Doray et al., 2018) is run in May by France and covers the French shelf of the Bay of Biscay. Its objectives are to monitor the Bay of Biscay pelagic ecosystem in springtime and assess the biomass of its small pelagic fish species, including sardine and anchovy. Amongst many information on the ecosystem, the survey PELGAS provides knowledge on the adults of anchovy and sardine during their spawning in spring. The survey JUVENA (Boyra et al., 2013) is run in September by Spain. It has a larger spatial coverage than PELGAS, including part of the Spanish coast and open ocean outside the shelf because it targets juvenile anchovy. It also provides knowledge on sardine as well as other pelagic species. Both surveys are coordinated by the ICES Working Group on Acoustic and Egg Surveys for Small Pelagic Fish (WGACEGG), together with other pelagic surveys in ICES areas 7, 8 and 9. Survey protocoles are detailed in Doray et al. (2021). Briefly, fish backscatter data are recorded along the survey transect lines and pelagic trawl hauls are undertaken opportunistically to identify the echotraces to species and collect fish samples for biometric data. The trawl haul catches have provided the anchovy and sardine length data, from which the maps presented here are derived. At each trawl haul, the catch is sorted by species and weighted. A subsample by species is measured to estimate the species’ length distribution. The four maps presented here correspond to the average maps of anchovy and sardine length distributions in May and September, derived from the PELGAS and JUVENA trawl haul data series. The maps were obtained by kriging, following the procedure explained in Petitgas et al. (2011) for mapping functions instead of variables. For each species, the experimental length distribution at each haul was fitted by a linear combination of Legendre polynomials, the coefficients of which were co-kriged. The number of polynomials varied from 15 to 22 depending on the survey and species, with a higher number for sardine and in autumn. The length histogram at each grid node was then deduced from the mapped coefficients. When the length distribution at a given haul was estimated with less than 40 individual fish, the haul was not taken into account for mapping. This threshold defined presence and absence of the species in the haul data sets. The trawl hauls from 2000 to 2019 were pooled for the PELGAS series (1965 stations) and from 2003 to 2020 for JUVENA (852 stations). The mapping was performed on the same grid for both PELGAS and JUVENA and both species, and with similar moving kriging neighbourhoods. The grid has a mesh size of 0.25 x 0.25 decimal degree square. In addition to mapping the length distribution, presence/absence was mapped for each species by ordinary kriging on the same grid and with the same neighbourhoods as previously. The computations were performed in R (version 4.0.5) with the RGeostats package (version 13.0.1) freely available at http://rgeostats.free.fr. The map data files comprise for each species the following information: the geographical coordinates of the grid points, the probability of presence and the probability of each length class of width 0.5 cm ranging from 2.5 to 27 cm.

Moving 6-year analysis and visualization of Water body dissolved oxygen concentration in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

'''Short description:''' For the Global Ocean- Sea Surface Temperature L3 Observations . This product provides daily foundation sea surface temperature from multiple satellite sources. The data are intercalibrated. This product consists in a fusion of sea surface temperature observations from multiple satellite sensors, daily, over a 0.05° resolution grid. It includes observations by polar orbiting from the ESA CCI / C3S archive . The L3S SST data are produced selecting only the highest quality input data from input L2P/L3P images within a strict temporal window (local nightime), to avoid diurnal cycle and cloud contamination. The observations of each sensor are intercalibrated prior to merging using a bias correction based on a multi-sensor median reference correcting the large-scale cross-sensor biases. '''DOI (product) :''' https://doi.org/10.48670/mds-00329

'''This product has been archived''' This dataset provide a times series of gap free map of Sea Surface Temperature (SST) foundation at high resolution on a 0.10 x 0.10 degree grid (approximately 10 x 10 km) for the Global Ocean, every 24 hours. Whereas along swath observation data essentially represent the skin or sub-skin SST, the Level 4 SST product is defined to represent the SST foundation (SSTfnd). SSTfnd is defined within GHRSST as the temperature at the base of the diurnal thermocline. It is so named because it represents the foundation temperature on which the diurnal thermocline develops during the day. SSTfnd changes only gradually along with the upper layer of the ocean, and by definition it is independent of skin SST fluctuations due to wind- and radiation-dependent diurnal stratification or skin layer response. It is therefore updated at intervals of 24 hrs. SSTfnd corresponds to the temperature of the upper mixed layer which is the part of the ocean represented by the top-most layer of grid cells in most numerical ocean models. It is never observed directly by satellites, but it comes closest to being detected by infrared and microwave radiometers during the night, when the previous day's diurnal stratification can be assumed to have decayed. The processing combines the observations of multiple polar orbiting and geostationary satellites, embedding infrared of microwave radiometers. All these sources are intercalibrated with each other before merging. A ranking procedure is used to select the best sensor observation for each grid point. An optimal interpolation is used to fill in where observations are missing. '''DOI (product) :''' https://doi.org/10.48670/mds-00321

Moving 6-year analysis and visualization of Water body silicate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.