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

This dataset provide a times series of daily mean fields of Sea Surface Temperature (SST) foundation at ultra-high resolution (UHR) on a 0.02 x 0.02 degree grid (approximately 2 x 2 km) for the Mediterranean Sea, 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. 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. The processing is the same as for the Atlantic Near Real Time (NRT) L4 dataset available on Copernicus Marine Service [SST_ATL_SST_L4_NRT_OBSERVATIONS_010_025 dataset] and users can refer to the user manual and quality documents available there for more details. This dataset was developed in the frame of European Space Agency (ESA)'s Medspiration project.

This daily High-Resolution (HR) Level 3 gridded wind product is derived from Copernicus Sentinel-1 SAR (Synthetic Aperture Radar) observations, over the North Western Atlantic ("ATL" area). It is based on the European Space Agency (ESA) Level-2 OCN products at the highest available resolution. Although L2-OCN products already contain wind vectors, those are calculated using the CMOD5.n Geophysical Model Function (GMF) applied to the co-polarized (co-pol) VV channel (emitting in Vertical polarization and receiving in Vertical polarization). This VV GMF was mapped from scatterometer sensors (Hersbach et al., 2007) which are only able to use co-pol measurements. However, these co-pol GMF are known to lose sensitivity for wind above 20 m/s. Therefore, wind based on such GMF alone, are known to under-estimate wind speed (Polverari et al., 2022). For the L3 products winds based on SAR, we take advantage of the available cross-polarized (cross-pol) VH channel (emitting in Vertical polarization and receiving in Horizontal polarization) for which GMF were specifically derived based on C-Band SAR (Mouche et al., 2017, Mouche et al., 2019). Winds estimated from the combination of both the co-pol and cross-pol channels are referred to as dual-polarization (or dual-pol) winds. As shown in Mouche et al. (2019), taking advantage of the dual polarization strongly improves the wind estimation for high wind conditions thanks to the much greater VH channel sensitivity compared to VV. These new wind estimations are then gridded with a 0.012 degree resolution (between 0.5 and 1.2 km in zonal direction depending on the latitude and 1.3 km in meridional direction) using a cylindrical equidistant projection, independently for ascending and descending satellite passes and for each satellite (so 4 wind fields are available per day for two satellites). This dataset is generated over all Sentinel-1 mission time series starting from March 2018 and updated in delayed mode with a 4-months delay. It is also produced for 4 other different European areas. This dataset is produced and disseminated in the frame of Copernicus Marine Service.

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

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

The Level 4 merged microwave wind product is a complete set of hourly global 10-m wind maps on a 0.25x0.25 degree latitude-longitude grid, spanning 1 Jan 2010 through the end of 2020. The product combines background neutral equivalent wind fields from ERA5, daily surface current fields from CMEMS, and stress equivalent winds obtained from several microwave passive and active sensors to produce hourly surface current relative stress equivalent wind analyses. The satellite winds include those from recently launched L-band passive sensors capable of measuring extreme winds in tropical cyclones, with little or no degradation from precipitation. All satellite winds used in the analyses have been recalibrated using a large set of collocated satellite-SFMR wind data in storm-centric coordinates. To maximize the use of the satellite microwave data, winds within a 24-hour window centered on the analysis time have been incorporated into each analysis. To accomodate the large time window, satellite wind speeds are transformed into deviations from ERA5 background wind speeds interpolated to the measurement times, and then an optical flow-based morphing technique is applied to these wind speed increments to propagate them from measurement to analysis time. These morphed wind speed increments are then added to the background wind speed at the analysis time to yield a set of total wind speeds fields for each sensor at the analysis time. These individual sensor wind speed fields are then combined with the background 10-m wind direction to yield vorticity and divergence fields for the individual sensor winds. From these, merged vorticity and divergence fields are computed as a weighted average of the individual vorticity and divergence fields. The final vector wind field is then obtained directly from these merged vorticity and divergence fields. Note that one consequence of producing the analyses in terms of vorticity and divergence is that there are no discontinuities in the wind speed fields at the (morphed) swath edges. There are two important points to be noted: the background ERA5 wind speed fields have been rescaled to be globally consistent with the recalibrated AMSR2 wind speeds. This rescaling involves a large increase in the ERA5 background winds beyond about 17 m/s. For example, an ERA5 10 m wind speed of 30 m/s is transformed into a wind speed of 41 m/s, and a wind speed of 34 m/s is transformed into a wind speed of about 48 m/s. Besides the current version of the product is calibrated for use within tropical cyclones and is not appropriate for use elsewhere. This dataset was produced in the frame of ESA MAXSS project. The primary objective of the ESA Marine Atmosphere eXtreme Satellite Synergy (MAXSS) project is to provide guidance and innovative methodologies to maximize the synergetic use of available Earth Observation data (satellite, in situ) to improve understanding about the multi-scale dynamical characteristics of extreme air-sea interaction.

This dataset contains the high-frequency total horizontal current at 15m depth on a global grid at 1/4° resolution. It is composed by the addition of two components, the first is the Geostrophic current derived by Altimetry, from the DT-2018 CMEMS database, and the second is the unsteady-Ekman ageostrophic component forced by the wind. All the details about the algorithm and the physical content of this ageostrophy component are given in the ATBD. The data are available through HTTP and FTP; access to the data is free and open. This dataset was generated by Datlas Ocean and is distributed by Ifremer / CERSAT in the frame of the World Ocean Circulation (WOC) project funded by the European Space Agency (ESA).

Level 2 skin Sea Surface Temperature derived from IASI on Metop-B, global and provided in full-resolution swath (12 km at nadir to 40 km), in GHRSST compliant netCDF format. SST is retrieved 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.

“SMOS Near Real Time Level 2 swath wind speed” (SMOS L2WS NRT) products, which are SMOS retrieved surface wind speed gridded maps with a spatial sampling of 1/4°x 1/4° and consisting of orbital segments (containing parts of ascending and descending half orbits) following the granularity of the SMOS Level 1B near real time (NRT) data products. SMOS L2WS NRT products are available within 4 to 6 hours from sensing and are generated in NetCDF-4 format as described in the "SMOS Wind Data Service Product Description Document" (http://www.smosstorm.org/Document-tools/SMOS-Wind-Data-Service-Documentation).

This daily High-Resolution (HR) Level 3 gridded wind product is derived from Copernicus Sentinel-1 SAR (Synthetic Aperture Radar) observations, over the Mediterranean Sea ("MED" area). It is based on the European Space Agency (ESA) Level-2 OCN products at the highest available resolution. Although L2-OCN products already contain wind vectors, those are calculated using the CMOD5.n Geophysical Model Function (GMF) applied to the co-polarized (co-pol) VV channel (emitting in Vertical polarization and receiving in Vertical polarization). This VV GMF was mapped from scatterometer sensors (Hersbach et al., 2007) which are only able to use co-pol measurements. However, these co-pol GMF are known to lose sensitivity for wind above 20 m/s. Therefore, wind based on such GMF alone, are known to under-estimate wind speed (Polverari et al., 2022). For the L3 products winds based on SAR, we take advantage of the available cross-polarized (cross-pol) VH channel (emitting in Vertical polarization and receiving in Horizontal polarization) for which GMF were specifically derived based on C-Band SAR (Mouche et al., 2017, Mouche et al., 2019). Winds estimated from the combination of both the co-pol and cross-pol channels are referred to as dual-polarization (or dual-pol) winds. As shown in Mouche et al. (2019), taking advantage of the dual polarization strongly improves the wind estimation for high wind conditions thanks to the much greater VH channel sensitivity compared to VV. These new wind estimations are then gridded with a 0.012 degree resolution (between 0.5 and 1.2 km in zonal direction depending on the latitude and 1.3 km in meridional direction) using a cylindrical equidistant projection, independently for ascending and descending satellite passes and for each satellite (so 4 wind fields are available per day for two satellites). This dataset is generated over all Sentinel-1 mission time series starting from March 2018 and updated in delayed mode with a 4-months delay. It is also produced for 4 other different European areas. This dataset is produced and disseminated in the frame of Copernicus Marine Service.