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

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

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 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 DBCP – Data Buoy Cooperation Panel - is an international program coordinating the use of autonomous data buoys to observe atmospheric and oceanographic conditions, over ocean areas where few other measurements are taken. DBCP coordinates the global array of 1 600 active drifting buoys (August 2020) and historical observation from 14 000 drifting buoys. Data and metadata collected by drifting buoys are publically available in near real-time via the Global Data Assembly Centers (GDACs) in Coriolis-Ifremer (France) and MEDS (Canada) after an automated quality control (QC). In long term, scientifically quality controlled delayed mode data will be distributed on the GDACs. Disclaimer: the DB-GDAC is under construction. It is currently (January 2020) aggregating data from the Coriolis DAC (E-Surfmar, Canada). Additional DACs are considered. An interim provision from GTS real-time data to GDAC may be provided from Coriolis DAC.  

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

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

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 dataset provides a global Look-Up Table (LUT) of physiological ratios for the real-time adjustment of chlorophyll-a fluorescence measured by biogeochemical Argo (BGC-Argo) profiling floats. The physiological ratios aim to account for the global variability in the relationship between fluorescence and chlorophyll-a concentration, as influenced by phytoplankton physiology. The LUT was developed using two different gap-filled observational Argo-based products (SOCA machine learning-based methodology ; Sauzède et al., 2016; Sauzède et al., 2024). The first product provides gap-filled chlorophyll-a data derived from fluorescence corrected for dark signal and non-photochemical quenching (NPQ) following Schmechtig et al. (2023), while the second product provides chlorophyll-a concentrations derived from light attenuation. The latter is based on the downward irradiance at 490 nm (ED490) derived from the SOCA-light method (Renosh et al., 2023). From this, the diffuse attenuation coefficient (KD490) is computed, which is subsequently used to estimate the chlorophyll-a concentration through the bio-optical relationships described by Morel et al. (2007). These two products, based on fluorescence and radiometry, enable the derivation of spatially varying correction factors, or physiological ratios. These ratios provide a validated grounded framework for adjusting real-time fluorescence observations from OneArgo floats into chlorophyll-a concentrations. The LUT is distributed in NetCDF format and is provided on a regular 1°×1° latitude–longitude grid covering the global ocean. Each grid cell contains the temporal mean, averaged over the water column (from the surface to 1.5 times the euphotic depth), of the physiological ratio. The file also includes metadata describing variable definitions, units, and other relevant information. Variables included: - physiological_ratio — fluorescence-to-radiometry-based chlorophyll correction factor (dimensionless) - physiological_ratio_sd — temporal standard deviation (over the twelve months) of the fluorescence-to-radiometry-based chlorophyll correction factor (dimensionless) - lat, lon — spatial coordinates (degrees north/east) - Global attributes — dataset description, reference citation, and contact information