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

The Southern Ocean plays a fundamental role in regulating the global climate. This ocean also contains a rich and highly productive ecosystem, potentially vulnerable to climate change. Very large national and international efforts are directed towards the modeling of physical oceanographic processes to predict the response of the Southern Ocean to global climate change and the role played by the large-scale ocean climate processes. However, these modeling efforts are greatly limited by the lack of in situ measurements, especially at high latitudes and during winter months. The standard data that are needed to study ocean circulation are vertical profiles of temperature and salinity, from which we can deduce the density of seawater. These are collected with CTD (Conductivity-Temperature-Depth) sensors that are usually deployed on research vessels or, more recently, on autonomous Argo profilers. The use of conventional research vessels to collect these data is very expensive, and does not guarantee access to areas where sea ice is found at the surface of the ocean during the winter months. A recent alternative is the use of autonomous Argo floats. However, this technology is not easy to use in glaciated areas. In this context, the collection of hydrographic profiles from CTDs mounted on marine mammals is very advantageous. The choice of species, gender or age can be done to selectively obtain data in particularly under-sampled areas such as under the sea ice or on continental shelves. Among marine mammals, elephant seals are particularly interesting. Indeed, they have the particularity to continuously dive to great depths (590 ± 200 m, with maxima around 2000 m) for long durations (average length of a dive 25 ± 15 min, maximum 80 min). A Conductivity-Temperature-Depth Satellite Relay Data Logger (CTD-SRDLs) has been developed in the early 2000s to sample temperature and salinity vertical profiles during marine mammal dives (Boehme et al. 2009, Fedak 2013). The CTD-SRDL is attached to the seal on land, then it records hydrographic profiles during its foraging trips, sending the data by satellite ARGOS whenever the seal goes back to the surface.While the principle intent of seal instrumentation was to improve understanding of seal foraging strategies (Biuw et al., 2007), it has also provided as a by-product a viable and cost-effective method of sampling hydrographic properties in many regions of the Southern Ocean (Charrassin et al., 2008; Roquet et al., 2013).

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

The glider operations in the MOOSE network started to be deployed regularly in 2010 in the North Western Mediterranean Sea, thanks to the setup of national glider facilities at DT-INSU/Ifremer (http://www.dt.insu.cnrs.fr/gliders/gliders.php) and with the support of the European project FP7-PERSEUS. Two endurance lines are operated: MooseT00 (Nice-Calvi; Ligurian Sea) and MooseT02 (Marseille-Menorca; Gulf of Lion). The all dataset here corresponds to raw data in the EGO format.

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 dataset is an aggregation of all availale in situ data from Coriolis and Copernicus in situ data centres, observed in the French DCSMM area. It contains 5167 NetCDF CF files from 1903 to 2017. Each file contains the observations of a specific platform (e.g. vessel, mooring site, sea level station). Observed parameters are temperature, salinity, pressure, oxygen, nitrate, chlorophyll (and other bio-geo-chemicals), current, wave, sea level, river flow.  

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.