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

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 )

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

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.

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

Argo is a global array of 3,000 free-drifting profiling floats that measures the temperature and salinity of the upper 2000 m of the ocean. This allows, for the first time, continuous monitoring of the temperature, salinity, and velocity of the upper ocean, with all data being relayed and made publicly available within hours after collection. The array provides 100,000 temperature/salinity profiles and velocity measurements per year distributed over the global oceans at an average of 3-degree spacing. Some floats provide additional bio-geo parameters such as oxygen or chlorophyll. All data collected by Argo floats are publically available in near real-time via the Global Data Assembly Centers (GDACs) in Brest (France) and Monterey (California) after an automated quality control (QC), and in scientifically quality controlled form, delayed mode data, via the GDACs within six months of collection. The BGC-Argo Sprof snapshot is a subset of the global Argo data snapshot.  It is created to ease BGC-Argo data usage.  The content is the same if you are to download the global Argo data snapshot, and then select all the BGC-Argo Sprof files.  Please use the same DOI and citation as the global Argo data snapshot.

Ensemble simulations of the ecosystem model Apecosm (https://apecosm.org) forced by the IPSL-CM6-LR climate model with the climate change scenario SSP1-2.6. 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.

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