These monthly gridded climatology were produced using MBT, XBT, Profiling floats, Gliders, and ship-based CTD data from different database and carried out in the Med. between 1969 and 2013. The Mixed Layer Depth (MLD) is calculated with a delta T= 0.1 C criterion relative to 10m reference level on individual profiles. The Depth of the Bottom of the Seasonal Thermocline (DBST) is calculated on individual profiles as the maximum value from a vector composed of two elements: 1) the depth of the temperature minimum in the upper 200m; 2) the MLD. This double criterion for the calculation of DBST is necessary in areas where the mixed layer exceed 200m depth. DBST is the integration depth used in the calculation of the upper-ocean Heat Storage Rate. For more details about the data and the methods used, see: Houpert et al. 2015, Seasonal cycle of the mixed layer, the seasonal thermocline and the upper-ocean heat storage rate in the Mediterranean Sea derived from observations, Progress in Oceanography, http://doi.org/10.1016/j.pocean.2014.11.004

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.

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

Ensemble simulations of the ecosystem model Apecosm (https://apecosm.org) forced by the IPSL-CM6-LR climate model with the climate change scenario SSP5-8.5. The output files contain yearly mean biomass density for 3 communities (epipelagic, mesopelagic migratory and mesopelagic redidents) and 100 size classes (ranging from 0.12cm to 1.96m) The model grid file is also provided. Units are in J/m2 and can be converted in kg/m2 by dividing by 4e6. These outputs are associated with the "Assessing the time of emergence of marine ecosystems from global to local scales using IPSL-CM6A-LR/APECOSM climate-to-fish ensemble simulations" paper from the Earth's Future "Past and Future of Marine Ecosystems" Special Collection.

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.

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.

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).

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.

15 years of L-Band remote sensing Sea Surface Salinity (SSS) measurements have proven the capability of satellite SSS to resolve large scale to mesoscale SSS features in tropical to subtropical ocean. In mid to high latitude, L-Band measurements still suffer from large scale and time varying biases. Here, a simple method is proposed to mitigate the large scale and time varying biases. First, in order to estimate these biases, an Optimal Interpolation (OI) using a large correlation scale is used to map SMOS and SMAP L3 products and is compared to equivalent mapping of in situ observations. Then, a second mapping is performed on corrected SSS at scale of SMOS/SMAP resolution (~45 km). This procedure allows to correct and merge both products, and to increase signal to noise ratio of the absolute SSS estimates. Using thermodynamic equation of state (TEOS-10), the resulting L4 SSS product is combined with microwave satellite SST products to produce sea surface density and spiciness, useful to fully characterize the surface ocean water masses. The new L4 SSS products is validated against independent in situ measurements from low to high latitudes. The L4 products exhibits a significant improvement in mid-and high latitude in comparison to the existing SMOS and SMAP L3 products.

This delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for ocean surface currents and current vertical profiles. It concerns three delayed time datasets dedicated to near-surface currents measurements coming from three platforms (Lagrangian surface drifters, High Frequency radars and Argo floats) and velocity profiles within the water column coming from the Acoustic Doppler Current Profiler (ADCP, vessel mounted only). The latest version of Copernicus surface and sub-surface water velocity product is also distributed from Copernicus Marine catalogue.