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 )

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

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 data set contains the gridded hydrographic and transport data for the biennial Go-Ship A25 Greenland–Portugal OVIDE section from 2002 to 2012. The properties and transports are mapped on a 7km x 1m grid. Using a common grid facilitates the comparison between the different occupations of the line and the averaging. This data set was used in Daniault et al. (2016, Progress in Oceanography) to which the reader is referred for a description of the gridding method.

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.

This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2023 and GLobal Ocean Data Analysis Project GLODAPv2.2023. The SOCATv2023-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2023. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2023-GRIDDED: monthly, yearly and decadal averages of fCO2 over a 1x1 degree grid over the global ocean, and a 0.25x0.25 degree, monthly average for the coastal ocean. GLODAPv2.2023-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2021. These data were subjected to an extensive quality control and bias correction described in Olsen et al. (2020). GLODAPv2-GRIDDED contains global climatologies for temperature, salinity, oxygen, nitrate, phosphate, silicate, dissolved inorganic carbon, total alkalinity and pH over a 1x1 degree horizontal grid and 33 standard depths using the observations from the previous major iteration of GLODAP, GLODAPv2. SOCAT and GLODAP are based on community, largely volunteer efforts, and the data providers will appreciate that those who use the data cite the corresponding articles (see References below) in order to support future sustainability of the data products.

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.