The 4D Marine Heatwaves (MHW) atlas contains 4D (x, y, z, t) **daily temperature and marine heatwaves categories** for global region [82.875°S-89.875°N, 0.125°E-359.875°E], from 0 to 300m depth and a spatial resolution of 1/8°. It covers the period 1993-2022. The MHW atlas has been computed from the temperature 4D fields of the ARMOR3D global product delivered in the Copernicus Marine Service (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012 - https://doi.org/10.48670/moi-00052 ). The MHW categories are derived from the Hobday’s method [Hobday et al.,2018]. Each MHW event is classified among four categories (moderate to extreme), identified in terms of multiples of the local difference between the 90th percentile and climatological values, and defined as moderate (1-2×, Category I), strong (2-3×, Category II), severe (3-4×, Category III), and extreme (>4×, Category IV). When the category is zero, this means that there is no MHW. The period 1993-2021 is used as a baseline for defining the climatology to be as close as possible to the 30-year period suggested by Hobday. This choice is motivated by the need of altimetry data to constrain the vertical temperature reconstruction, which is required for most ocean reanalyses as well. Additionally, ancillary data are provided together with the data. It consists of 4D daily **temperature climatology** and **90 percentiles of the temperature**. These fields have been used to compute the MHW categories. They are delivered over the same domain as the MHW atlas. ARMOR3D **temperature uncertainties** are also supplied as they can help users to select only the most reliable events in the database. This dataset was generated by CLS (Collecte Localisation satellite) and is distributed by Ifremer /CERSAT in the frame of the CAREHeat project (CAREHeat Website) funded by the European Space Agency (ESA).

In recent years, large datasets of in situ marine carbonate system parameters (partial pressure of CO2 (pCO2), total alkalinity, dissolved inorganic carbon and pH) have been collated. These carbonate system datasets have highly variable data density in both space and time, especially in the case of pCO2, which is routinely measured at high frequency using underway measuring systems. This variation in data density can create biases when the data are used, for example for algorithm assessment, favouring datasets or regions with high data density. A common way to overcome data density issues is to bin the data into cells of equal latitude and longitude extent. This leads to bins with spatial areas that are latitude and projection dependent (eg become smaller and more elongated as the poles are approached). Additionally, as bin boundaries are defined without reference to the spatial distribution of the data or to geographical features, data clusters may be divided sub-optimally (eg a bin covering a region with a strong gradient). To overcome these problems and to provide a tool for matching in situ data with satellite, model and climatological data, which often have very different spatiotemporal scales both from the in situ data and from each other, a methodology has been created to group in situ data into ‘regions of interest’, spatiotemporal cylinders consisting of circles on the Earth’s surface extending over a period of time. These regions of interest are optimally adjusted to contain as many in situ measurements as possible. All in situ measurements of the same parameter contained in a region of interest are collated, including estimated uncertainties and regional summary statistics. The same grouping is done for each of the other datasets, producing a dataset of matchups. About 35 million in situ datapoints were then matched with data from five satellite sources and five model and re-analysis datasets to produce a global matchup dataset of carbonate system data, consisting of 287,000 regions of interest spanning 54 years from 1957 to 2020. Each region of interest is 100 km in diameter and 10 days in duration. An example application, the reparameterisation of a global total alkalinity algorithm, is shown. This matchup dataset can be updated as and when in situ and other datasets are updated, and similar datasets at finer spatiotemporal scale can be constructed, for example to enable regional studies. This dataset was funded by ESA Satellite Oceanographic Datasets for Acidification (OceanSODA) project which aims at developing the use of satellite Earth Observation for studying and monitoring marine carbonate chemistry.

This dataset provide a times series of daily mean fields of Sea Surface Temperature (SST) foundation at high resolution (UHR) on a 0.10 x 0.10 degree grid (approximately 10 x 10 km) for the Global Ocean, every 24 hours. An Optimal interpolation (OI) technique is used to combine coincident swath measures of SST from different types satellite sensors and to fill gaps where no observations are available or obstructed by clouds. This multi-sensor compositing and interpolation process categorizes this dataset as a Level 4 product. Each daily mean field is available the next day before 12:00. Whereas along swath observation data essentially represent the skin or sub-skin SST, the L4 SST product is defined to represent the SST foundation (SSTfnd). SSTfnd is defined within GHRSST-PP as the temperature at the base of the diurnal thermocline. It is so named because it represents the foundation temperature on which the diurnal thermocline develops during the day. SSTfnd changes only gradually along with the upper layer of the ocean, and by definition it is independent of skin SST fluctuations due to wind- and radiation-dependent diurnal stratification or skin layer response. It is therefore updated at intervals of 24 hrs. SSTfnd corresponds to the temperature of the upper mixed layer which is the part of the ocean represented by the top-most layer of grid cells in most numerical ocean models. It is never observed directly by satellites, but it comes closest to being detected by infrared and microwave radiometers during the night, when the previous day's diurnal stratification can be assumed to have decayed. The processing combines the observations of multiple polar orbiting and geostationary satellites, embedding infrared of microwave radiometers. All these sources are intercalibrated with eachother before merging. This dataset is generated daily within a 24 delay and is therefore suitable for assimilation into operational models. It is produced in the frame of Copernicus Marine Service and the data available through various tools and protocols with a simple user registration from this service (product identifier: SST_GLO_PHY_L4_NRT_010_043) at: https://data.marine.copernicus.eu/product/SST_GLO_PHY_L4_NRT_010_043

Level 3 hourly sub-skin Sea Surface Temperature derived from Meteosat at 0° longitude, covering 60S-60N and 60W-60E and re-projected on a 0.05° regular grid, in GHRSST compliant netCDF format. The satellite input data has successively come from Meteosat level 1 data processed at EUMETSAT. SST is retrieved from SEVIRI using a multi-spectral algorithm and a cloud mask. Atmospheric profiles of water vapor and temperature from a numerical weather prediction model, Sea Surface Temperature from an analysis, together with a radiative transfer model, are used to correct the multispectral algorithm for regional and seasonal biases due to changing atmospheric conditions.The quality of the products is monitored regularly by daily comparison of the satellite estimates against buoy measurements. The product format is compliant with the GHRSST Data Specification (GDS) version 2. Users are advised to use data only with quality levels 3, 4 and 5.

Level 3, four times a day, sub-skin Sea Surface Temperature derived from AVHRR on Metop satellites and VIIRS or AVHRR on NOAA and NPP satellites, over North Atlantic and European Seas and re-projected on a polar stereographic at 2 km resolution, in GHRSST compliant netCDF format. This catalogue entry presents NOAA-20 North Atlantic Regional Sea Surface Temperature. SST is retrieved from infrared channels using a multispectral algorithm and a cloud mask. Atmospheric profiles of water vapor and temperature from a numerical weather prediction model, Sea Surface Temperature from an analysis, together with a radiative transfer model, are used to correct the multispectral algorithm for regional and seasonal biases due to changing atmospheric conditions. The quality of the products is monitored regularly by daily comparison of the satellite estimates against buoy measurements. The product format is compliant with the GHRSST Data Specification (GDS) version 2.Users are advised to use data only with quality levels 3,4 and 5.

Level 3 twice-daily sub-skin Sea Surface Temperature derived from AVHRR on Metop-A, global and re-projected on a 0.05° regular grid, in GHRSST compliant netCDF format. The satellite input data has successively come from Metop-A and Metop-B level 1 data processed at EUMETSAT. SST is retrieved from AVHRR infrared channels (3.7, 10.8 and 12.0 µm) using a multi-spectral algorithm and a cloud mask. Atmospheric profiles of water vapor and temperature from a numerical weather prediction model, Sea Surface Temperature from an analysis, together with a radiative transfer model, are used to correct the multispectral algorithm for regional and seasonal biases due to changing atmospheric conditions.The quality of the products is monitored regularly by daily comparison of the satellite estimates against buoy measurements. The product format is compliant with the GHRSST Data Specification (GDS) version 2. Users are advised to use data only with quality levels 3, 4 and 5.

The ESA Sea State Climate Change Initiative (CCI) project has produced global daily merged multi-sensor time-series of along-track satellite altimeter significant wave height data (referred to as Level 3 (L3) data) with a particular focus for use in climate studies. This dataset contains the Version 3 Remote Sensing Significant Wave Height product, which provides along-track data at approximately 6 km spatial resolution. It has been generated from upstream Sea State CCI L2P products, edited and merged into daily products, retaining only valid and good quality measurements from all altimeters over one day, with simplified content (only a few key parameters). This is close to what is delivered in Near-Real Time by the CMEMS (Copernicus - Marine Environment Monitoring Service) project. It covers the date range from 2002-2021. The altimeter data used in the Sea State CCI dataset v3 come from multiple satellite missions (Envisat, CryoSat-2, Jason-1, Jason-2, Jason-3, SARAL, Sentinel-3A), therefore spanning over a shorter time range than version 1.1. Unlike version 1.1, this version 3 involved a complete and consistent retracking of all the included altimeters. Many altimeters are bi-frequency (Ku-C or Ku-S) and only measurements in Ku band were used, for consistency reasons, being available on each altimeter but SARAL (Ka band).

This dataset provides the meridional and zonal components of both the stress-equivalent wind (U10S) and wind stress (Tau) vectors. The ERA* product is a correction of the ECMWF Fifth Reanalysis (ERA5) output by means of geo-located scatterometer-ERA5 differences over a 15-day temporal window. The product also contains ERA5 U10S and Tau. The data are available through HTTP and FTP; access to the data is free and open. In order to be informed about changes and to help us keep track of data usage, we encourage users to register at: https://forms.ifremer.fr/lops-siam/access-to-esa-world-ocean-circulation-project-data/ This dataset was generated by ICM (Institute of Marine Sciences) / CSIC (Consejo Superior de Investigaciones Científicas) and is distributed by Ifremer / CERSAT in the frame of the World Ocean Circulation (WOC) project funded by the European Space Agency (ESA).

The ESA Sea State Climate Change Initiative (CCI) project has produced global multi-sensor time-series of along-track satellite altimeter significant wave height data (referred to as Level 2P (L2P) data) with a particular focus for use in climate studies. This dataset contains the Version 3 Remote Sensing Significant Wave Height product, which provides along-track data at approximately 6 km spatial resolution, separated per satellite and pass, including all measurements with flags, corrections and extra parameters from other sources. These are expert products with rich content and no data loss. The altimeter data used in the Sea State CCI dataset v3 come from multiple satellite missions spanning from 2002 to 2022021 (Envisat, CryoSat-2, Jason-1, Jason-2, Jason-3, SARAL, Sentinel-3A), therefore spanning over a shorter time range than version 1.1. Unlike version 1.1, this version 3 involved a complete and consistent retracking of all the included altimeters. Many altimeters are bi-frequency (Ku-C or Ku-S) and only measurements in Ku band were used, for consistency reasons, being available on each altimeter but SARAL (Ka band).

Data record (2004-2012) of level 3 hourly sub-skin Sea Surface Temperature derived from Meteosat at 0° longitude, covering 60S-60N and 60W-60E and re-projected on a 0.05° regular grid, in GHRSST compliant netCDF format. The satellite input data come from the imager SEVIRI on MSG satellites (Meteosat-8 and Meteosat-9). SST is retrieved from SEVIRI infrared channels (10.8 and 12.0 µm) using a multispectral algorithm and the cloud mask from CM SAF. NWP outputs (temperature and humidity profiles), OSTIA Sea Surface Temperature re-analysis and analysis, together with a radiative transfer model (RTTOV), are used to correct the multispectral algorithm for regional and seasonal biases due to changing atmospheric conditions. The product format is compliant with the GHRSST Data Specification (GDS) version 2. Users are advised to use data only with quality levels 3,4 and 5.