Observations of Sea surface temperature and salinity are now obtained from voluntary sailing ships using medium or small size sensors. They complement the networks installed on research vessels or commercial ships. The delayed mode dataset proposed here is upgraded annually as a contribution to GOSUD (http://www.gosud.org )

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 dataset comprises two netcdf files. The first file contains the six global two-dimensional maps necessary to implement the tidal mixing parameterization presented in de Lavergne et al. (2020). Four power fields (E_wwi, E_sho, E_cri and E_hil) represent depth-integrated internal tide energy dissipation, with units of Watts per square meter. Each power field corresponds to a specific dissipative process and associated vertical structure of turbulence production. The two remaining fields, H_cri and H_bot, are decay heights (with units of meters) that enter the vertical structures of the E_cri and E_hil components, respectively. The second file contains three-dimensional fields of turbulence production (with units of Watts per kilogram) obtained by application of the parameterization to the WOCE global hydrographic climatology. The file includes the total turbulence production (epsilon_tid), its four components (epsilon_wwi, epsilon_sho, epsilon_cri, epsilon_hil), and the underlying hydrographic fields, as a function of longitude, latitude and depth. All maps have a horizontal resolution of 0.5º. Detailed documentation of the parameterization can be found in the following publication: de Lavergne, C., Vic, C., Madec, G., Roquet, F., Waterhouse, A.F., Whalen, C.B., Cuypers, Y., Bouruet-Aubertot, P., Ferron, B., Hibiya, T. A parameterization of local and remote tidal mixing. Journal of Advances in Modeling Earth Systems, 12, e2020MS002065 (2020). https://doi.org/10.1029/2020MS002065

A quantitative understanding of the integrated ocean heat content depends on our ability to determine how heat is distributed in the ocean and what are the associated coherent patterns. This dataset contains the results of the Maze et al., 2017 (Prog. Oce.) study demonstrating how this can be achieved using unsupervised classification of Argo temperature profiles. The dataset contains: - A netcdf file with classification~results (labels and probabilities) and coordinates (lat/lon/time) of 100,684 Argo temperature profiles in North Atlantic. - A netcdf file with a Profile Classification Model (PCM) that can be used to classify new temperature profiles from observations or numerical models. The classification method used is a Gaussian Mixture Model that decomposes the Probability Density Function of the dataset into a weighted sum of Gaussian modes. North Atlantic Argo temperature profiles between 0 and 1400m depth were interpolated onto a regular 5m grid, then compressed using Principal Component Analysis and finally classified using a Gaussian Mixture Model. To use the netcdf PCM file to classify new data, you can checkout our PCM Matlab and Python toolbox here: https://github.com/obidam/pcm

This dataset contains the dynamical outputs of a global ocean simulation coupling dynamics and biogeochemistry at ¼° over the year 2019. The simulation has been performed using the coupled circulation/ecosystem model NEMO/PISCES (https://www.nemo-ocean.eu/), which is here enhanced to perform an ensemble simulation with explicit simulation of modeling uncertainties in the physics and in the biogeochemistry. This dataset is one of the 40 members of the ensemble simulation. This study was part of the Horizon Europe project SEAMLESS (https://seamlessproject.org/Home.html), with the general objective of improving the analysis and forecast of ecosystem indicators.   See Popov et al. (https://os.copernicus.org/articles/20/155/2024/) for more details on the study.

This dataset provides a World Ocean Atlas of Argo inferred statistics. The primary data are exclusively Argo profiles. The statistics are done using the whole time range covered by the Argo data, starting in July 1997. The atlas is provided with a 0.25° resolution in the horizontal and 63 depths from 0 m to 2,000 m in the vertical. The statistics include means of Conservative Temperature (CT), Absolute Salinity, compensated density, compressiblity factor and vertical isopycnal displacement (VID); standard deviations of CT, VID and the squared Brunt Vaisala frequency; skewness and kurtosis of VID; and Eddy Available Potential Energy (EAPE). The compensated density is the product of the in-situ density times the compressibility factor. It generalizes the virtual density used in Roullet et al. (2014). The compressibility factor is defined so as to remove the dependency with pressure of the in-situ density. The compensated density is used in the computation of the VID and the EAPE.

This data set provides a monthly time series of the upper limb of the Meridional Overturning Circulation (MOC) intensity at the A25 Greenland-Portugal OVIDE line from 1993 to 2015. The MOC was derived by combining AVISO altimetry with ISAS temperature and salinity data. The reader is referred to Mercier et al. (2015, Progress in Oceanography) for a full description of the method.

This product integrates sea level observations aggregated and validated from the Regional EuroGOOS consortium (Arctic-ROOS, BOOS, NOOS, IBI-ROOS, MONGOOS) and Black Sea GOOS as well as from the Global telecommunication system (GTS) used by the Met Offices. The latest version of Copernicus delayed-mode Sea level product is also distributed from Copernicus Marine catalogue.

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.

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.