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

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 )

GOSUD aims at assembling in-situ observations of the world ocean surface collected by a variety of ships and at distributing quality controlled datasets.  At present time the variables considered by GOSUD are temperature and salinity. The GOSUD data are mostly collected using thermosalinographs (TSG) installed on research vessels, on commercial ships and in some cases on sailing exploration ships GOSUD manages both near-real time data and delayed mode (reprocessed) data.

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

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

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

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.

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

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.