This dataset contains (1) outputs of Species Distribution Models (SDMs) for marine fish and (2) estimation of species richness using those outputs. SDMs use correlative algorithms to link presences of species to the environment recorded in place and time of their observation, calculate their environmental niche, estimate the geographical location suitable for them (habitat suitability) and in fine their geographical distribution. Here, we downloaded presences of marine fish from two open source databases, GBIF and OBIS and 13 environmental predictors known to be relevant in fish ecology (list below). We used the CEPHALOPOD pipeline, a framework allowing the user the compute a lot of species at the same time, with comparable methods and a verification of quality of inputs and outputs at every steps (Schickele et al., 2025). 3,642 fish made it to the final step and have the habitat suitability estimated for 12 month + annual mean, 10 bootstrap to quantify uncertainty and x algorithms. Those estimation are available in the “L2_marine_fish_*.nc” files, organized by water column position (bathydemersal, bathypelagic, benthopelagic, demersal, pelagic-neritic, pelagic-oceanic, reef-associated). We then used their annual mean to estimate their actual geographic distribution by applying (1) a bathymetric filtration and (2) a cutting procedure which removes isolated patches of high suitability (i.e., potential distribution) with no recorded presences (i.e., considered to be outside of the species dispersion range). Those geographic distribution were then stacked to estimate global species richness of every fish and each water column position.

The upper ocean pycnocline (UOP) monthly climatology is based on the ISAS20 ARGO dataset containing Argo and Deep-Argo temperature and salinity profiles on the period 2002-2020. Regardless of the season, the UOP is defined as the shallowest significant stratification peak captured by the method described in Sérazin et al. (2022), whose detection threshold is proportional to the standard deviation of the stratification profile. The three main characteristics of the UOP are provided -- intensity, depth and thickness -- along with hydrographic variables at the upper and lower edges of the pycnocline, the Turner angle and density ratio at the depth of the UOP. A stratification index (SI) that evaluates the amount of buoyancy required to destratify the upper ocean down to a certain depth, is also included. When evaluated at the bottom of the UOP, this gives the upper ocean stratification index (UOSI) as discussed in Sérazin et al. (2022). Three mixed layer depth variables are also included in this dataset, including the one using the classic density threshold of 0.03 kg.m-3, along with the minimum of these MLD variables. Several statistics of the UOP characteristics and the associated quantities are available in 2°×2° bins for each month of the year, whose results were smoothed using a diffusive gaussian filter with a 500 km scale. UOP characteristics are also available for each profile, with all the profiles sorted in one file per month.

This dataset is composed by the climatological seasonal field of the Ocean Salinity Stratification as defined from the Brunt-Vaisala frequency limited to the upper 300 m depth. The details are given in Maes, C., and T. J. O’Kane (2014), Seasonal variations of the upper ocean salinity stratification in the Tropics, J. Geophys. Res. Oceans, 119, 1706–1722, doi:10.1002/2013JC009366.

The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML). All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoys. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. This product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.

This dataset contains OAC-P results from application to Argo data in the World Ocean : - the 2000-2015 climatology of OAC-P results mapped onto a 0.5x0.5 grid with mapping error estimates; - the 2000-2015 probability density function of the permanent pycnocline potential density referenced to the sea surface vs Brunt-Väisälä frequency squared.OAC-P is an "Objective Algorithm for the Characterization of the permanent Pycnocline" developed to characterize subtropical gyre stratification features with both observed and modeled potential density profiles. OAC-P estimates the following properties: - for the permanent pycnocline: depth, upper and lower thicknesses, Brunt-Väisälä frequency squared, potential density, temperature and salinity; - for the surface mode water overlying the permanent pycnocline: depth, Brunt-Väisälä frequency squared, potential density, temperature and salinity. Argo data were download from Coriolis Argo GDAC on February, 8th 2016. Only Argo data with QC=1, 2, 5 or 8 were used.

The Coriolis Ocean Dataset for Reanalysis for the Ireland-Biscay-Iberia region (hereafter CORA-IBI) product is a regional dataset of in situ temperature and salinity measurements. The latest version of the product covers the period 1950-2014. The CORA-IBI observations comes from many different sources collected by Coriolis data centre in collaboration with the In Situ Thematic Centre of the Copernicus Marine Service (CMEMS INSTAC).  The observations integrated in the CORA-IBI product have been acquired both by autonomous platforms (Argo profilers, fixed moorings, gliders, drifters, sea mammals, fishery observing system from the RECOPESCA program), research or opportunity vessels ( CTDs, XBTs, ferrybox).  This CORA-IBI product has been controlled using an objective analysis (statistical tests) method and a visual quality control (QC). This QC procedure has been developed with the main objective to improve the quality of the dataset to the level required by the climate application and the physical ocean re-analysis activities. It provides T and S individual profiles on their original level with QC flags. The reference level of measurements is immersion (in meters) or pressure (in decibars). It is a subset on the IBI (Iberia-Bay-of-Biscay Ireland) of the CMEMS product referenced hereafter. The main new features of this regional product compared with previous global CORA products are the incorporation of coastal profiles from fishery observing system (RECOPESCA programme) in the Bay of Biscay and the English Channel as well as the use of an historical dataset collected by the Service hydrographique de la Marine (SHOM).

Satellite altimetry missions provide a quasi-global synoptic view of sea level over more than 25 years. The satellite altimetry constellation is used to build sea level maps and regional sea level indicators such as trends and accelerations. Estimating realistic uncertainties on these quantities is crucial to address some current climate science questions such as climate change detection and attribution or regional sea level budget closure for example. Previous studies have estimated the uncertainty for the global mean sea level (GMSL), but no uncertainty information is available at regional scales. In this study we estimate a regional satellite altimetry error budget and use it to derive maps of confidence intervals for local sea rise rates and accelerations. We analyze 27 years of satellite altimetry maps and derive the satellite altimetry error variance-covariance matrix at each grid point, prior to the estimation of confidence intervals on local trends and accelerations at the 90% confidence level using extended least squares estimators. Over 1993–2019, we find that the average local sea level trend uncertainty is 0.83 mm.yr-1 with local values ranging from 0.78 to 1.22 mm.yr-1. For accelerations, uncertainties range from 0.057 to 0.12 mm.yr-2, with a mean value of 0.063 mm.yr-2.   Change history: - 2020/07/08: initial dataset submission over 1993-2018 - 2020/10/21: 1993-2019 update and addition of error levels

This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2023 and GLobal Ocean Data Analysis Project GLODAPv2.2023. The SOCATv2023-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2023. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2023-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2023-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2021. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous major iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products.

Mesoscale eddy detection from 2000 to 2021 are computed using the AMEDA algorithm applied on AVISO L4 absolute dynamic topography at 1/8th degree. Eddy numbers correspond to tracks referenced in the DYNED atlas (https://doi.org/10.14768/2019130201.2). Detection is based on AVISO delyed-time product from 2000 to 2019 and on day+6 near-real-time altimetry from 2020 to 2021. Colocalisation is then made with available in situ profiles from Coriolis Ocean Dataset for Reanalysis (CORA) delayed-time data (113486 profiles) and Copernicus near-real-time profiles (43567).

This In Situ delayed mode product integrates the best available version of in situ oxygen, chlorophyll / fluorescence and nutrients data. The latest version of Copernicus delayed-mode BGC (bio-geo-chemical) product is also distributed from Copernicus Marine catalogue.