'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Products ADG (volume absorption coefficient of radiative flux in sea water due to dissolved organic matter and non algal particles), APH (volume absorption coefficient of radiative flux in sea water due to phytoplankton) and ATOT (volume absorption coefficient of radiative flux in sea water) are described in the PML Inherent Optical Property model (Smyth, T.J., Moore, G.F., Hirata, T. Aiken, J. (2006), a semi-analytic model for the derivation of ocean color inherent optical properties. The RRS product is defined as the spectral ratio of upwelling radiance and downwelling irradiance which can also be expressed as the ratio of normalized water leaving Radiance (nLw) and the extra-terrestrial solar irradiance (F0). The KD490 product identifies the turbidity of the water column, i.e., how visible light in the blue-green region of the spectrum penetrates within the water column. It is directly related to the presence of scattering particles in the water column. Inorganic Suspended Particulate Matter (SPM) is defined as all inorganic matter that stays on a glass fibre filter with an approximate pore size of 0.7 micrometres. Heavy metals and various organic micropollutants adsorb to SPM, the transport of which can affect the ecosystem. High concentrations of SPM cause turbidity which in turn affects the underwater light conditions, thus influencing primary production by phytoplankton and other algae in coastal waters. Products derived from OLCI are Rrs (400, 412, 443, 490, 510, 560, 620, 665, 674, 681, 709) and KD490. From the CCI multiple-sensor product are derived Rrs (410, 443, 490, 510, 560, 665nm), the Inherent Optical Properties, IOPs (ADG, APH, ATOT) and SPM. These products are remapped at nominal 300m (OLCI) and 1 Km spatial resolution using cylindrical equirectangular projection. '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of in water absorption parameters can be derived. '''DOI (product) :''' https://doi.org/10.48670/moi-00076

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.

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the Global Ocean- In-situ observation delivered in delayed mode. This In Situ delayed mode product integrates the best available version of in situ oxygen, chlorophyll / fluorescence and nutrients data '''DOI (product) :''' https://doi.org/10.17882/86207

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the European Ocean- The L3 multi-sensor (supercollated) product is built from bias-corrected L3 mono-sensor (collated) products at the resolution 0.02 degrees. If the native collated resolution is N and N < 0.02 the change (degradation) of resolution is done by averaging the best quality data. If N > 0.02 the collated data are associated to the nearest neighbour without interpolation nor artificial increase of the resolution. A synthesis of the bias-corrected L3 mono-sensor (collated) files remapped at resolution R is done through a selection of data based on the following hierarchy: AVHRR_METOP_B, VIIRS_NPP, SLSTRA, SEVIRI, AVHRRL-19, MODIS_A, MODIS_T, AMSR2. This hierarchy can be changed in time depending on the health of each sensor. '''DOI (product) :''' https://doi.org/10.48670/moi-00163

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the North Atlantic and Arctic oceans, the ESA Ocean Colour CCI Remote Sensing Reflectance (merged, bias-corrected Rrs) data are used to compute surface Chlorophyll (mg m-3, 1 km resolution) using the regional OC5CCI chlorophyll algorithm. The Rrs are generated by merging the data from SeaWiFS, MODIS-Aqua, MERIS, VIIRS and OLCI-3A sensors and realigning the spectra to that of the MERIS sensor. The algorithm used is OC5CCI - a variation of OC5 (Gohin et al., 2002) developed by IFREMER in collaboration with PML. As part of this development, an OC5CCI look up table was generated specifically for application over OC-CCI merged daily remote sensing reflectances. The resulting OC5CCI algorithm was tested and selected through an extensive calibration exercise that analysed the quantitative performance against in situ data for several algorithms in these specific regions. Processing information: PML's Remote Sensing Group has the capability to automatically receive, archive, process and map global data from multiple polar-orbiting sensors in both near-real time and delayed time. OLCI products are downloaded at level-2 from CODA, the Copernicus Hub and/or via EUMETCAST. These products are remapped at nominal 300m and 1 Km spatial resolution using cylindrical equirectangular projection. Description of observation methods/instruments: Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in situ measurements, an estimate of chlorophyll content can be derived. '''Processing information:''' ESA OC-CCI Rrs raw data are provided by Plymouth Marine Laboratory, currently at 4km resolution globally. These are processed to produce chlorophyll concentration using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (70deg), large spacecraft zenith angle (56deg), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 560 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). For the regional products, a variant of the OC-CCI chain is run to produce high resolution data at the 1km resolution necessary. A detailed description of the ESA OC-CCI processing system can be found in OC-CCI (2014e). '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean colour which is affected by the presence of phytoplankton. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''DOI (product) :''' https://doi.org/10.48670/moi-00070

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' The low resolution ocean physics analysis and forecast for the North-West European Shelf is produced using a forecasting ocean assimilation model, with tides, at 7 km horizontal resolution. The ocean model is NEMO (Nucleus for European Modelling of the Ocean), using the 3DVar NEMOVAR system to assimilate observations. These are surface temperature, vertical profiles of temperature and salinity, and along track satellite sea level anomaly data. The model is forced by lateral boundary conditions from the UK Met Office North Atlantic Ocean forecast model and by the CMEMS Baltic forecast product [https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=BALTICSEA_ANALYSISFORECAST_PHY_003_006 BALTICSEA_ANALYSISFORECAST_PHY_003_006]. The atmospheric forcing is given by the operational UK Met Office Global Atmospheric model. The river discharge is from a daily climatology. Further details of the model, including the product validation are provided in the [http://catalogue.marine.copernicus.eu/documents/QUID/CMEMS-NWS-QUID-004-001.pdf CMEMS-NWS-QUID-004-001]. Products are provided as hourly instantaneous and daily 25-hour, de-tided, averages. The datasets available are temperature, salinity, horizontal currents, sea level, mixed layer depth, and bottom temperature. Temperature, salinity and currents, as multi-level variables, are interpolated from the model 51 hybrid s-sigma terrain-following system to 24 standard geopotential depths (z-levels). Grid-points near to the model boundaries are masked. The product is updated daily, providing a 6-day forecast and the previous 2-day assimilative hindcast. See [http://catalogue.marine.copernicus.eu/documents/PUM/CMEMS-NWS-PUM-004-001_002.pdf CMEMS-NWS-PUM-004-001_002] for further details. '''Associated products:''' This model is coupled with a biogeochemistry model (ERSEM) available as CMEMS product [https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=NWSHELF_ANALYSISFORECAST_BGC_004_002 NWSHELF_ANALYSISFORECAST_BGC_004_002] A reanalysis product is available from: [https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=NWSHELF_MULTIYEAR_PHY_004_009 NWSHELF_MULTIYEAR_PHY_004_009]. '''DOI (product) :''' https://doi.org/10.48670/moi-00057

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The ocean monitoring indicator on regional mean sea level is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Irish-Biscay-Iberian (IBI) region is derived from the average of the gridded sea level maps weighted by the cosine of the latitude. The annual and semi-annual periodic signals are removed (least square fit of sinusoidal function) and the time series is low-pass filtered (175 days cut-off). The curve is corrected for the regional mean effect of the Glacial Isostatic Adjustment (GIA) using the ICE5G-VM2 GIA model (Peltier, 2004). During 1993-1998, the Global men sea level (hereafter GMSL) has been known to be affected by a TOPEX-A instrumental drift (WCRP Global Sea Level Budget Group, 2018; Legeais et al., 2020). This drift led to overestimate the trend of the GMSL during the first 6 years of the altimetry record (about 0.04 mm/y at global scale over the whole altimeter period). A correction of the drift is proposed for the Global mean sea level (Legeais et al., 2020). Whereas this TOPEX-A instrumental drift should also affect the regional mean sea level (hereafter RMSL) trend estimation, currently this empirical correction is currently not applied to the altimeter sea level dataset and resulting estimated for RMSL. Indeed, the pertinence of the global correction applied at regional scale has not been demonstrated yet and there is no clear consensus achieved on the way to proceed at regional scale. Additionally, the estimation of such a correction at regional scale is not obvious, especially in areas where few accurate independent measurements (e.g. in situ)- necessary for this estimation - are available. The trend uncertainty is provided in a 90% confidence interval (Prandi et al., 2021). This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation considering to the altimeter period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not taken into account. '''CONTEXT''' The indicator on area averaged sea level is a crucial index of climate change, and individual components contribute to sea level rise, including expansion due to ocean warming and melting of glaciers and ice sheets (WCRP Global Sea Level Budget Group, 2018). According to the recent IPCC 6th assessment report, global mean sea level (GMSL) increased by 0.20 (0.15 to 0.25) m over the period 1901 to 2018 with a rate 25 of rise that has accelerated since the 1960s to 3.7 (3.2 to 4.2) mm yr-1 for the period 2006–2018. Human activity was very likely the main driver of observed GMSL rise since 1970 (IPCC WGII, 2021). The weight of the different contributions evolves with time and in the recent decades the mass change has increased, contributing to the on-going acceleration of the GMSL trend (IPCC, 2022a; Legeais et al., 2020; Horwath et al., 2022). At regional scale, sea level does not change homogenously, and RMSL rise can also be influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). Rising sea level can strongly affect population and infrastructures in coastal areas, increase their vulnerability and risks for food security, particularly in low lying areas and island states. Adverse impacts from floods, storms and tropical cyclones with related losses and damages have increased due to sea level rise, and increase their vulnerability and increase risks for food security, particularly in low lying areas and island states (IPCC, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). In IBI region, the RMSL trend is modulated by decadal variations. As observed over the global ocean, the main actors of the long-term RMSL trend are associated with anthropogenic global/regional warming. Decadal variability is mainly linked to the strengthening or weakening of the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). The latest is driven by the North Atlantic Oscillation (NAO) (e.g. Delworth and Zeng, 2016). Along the European coast, the NAO also influences the along-slope winds dynamic which in return significantly contributes to the local sea level variability observed (Chafik et al., 2019). '''CMEMS KEY FINDINGS''' Over the [1993/01/01, 2021/08/02] period, the basin-wide RMSL in the IBI area rises at a rate of 3.8  0.82 mm/year. '''DOI (product):''' https://doi.org/10.48670/moi-00252

"'Short description:''' The High-Resolution Ocean Colour (HR-OC) Consortium (Brockmann Consult, Royal Belgian Institute of Natural Sciences, Flemish Institute for Technological Research) distributes Level 4 (L4) Turbidity (TUR, expressed in FNU), Solid Particulate Matter Concentration (SPM, expressed in mg/l), particulate backscattering at 443nm (BBP443, expressed in m-1) and chlorophyll-a concentration (CHL, expressed in µg/l) for the Sentinel 2/MSI sensor at 100m resolution for a 20km coastal zone. The products are delivered on a geographic lat-lon grid (EPSG:4326). BBP443, constitute the category of the 'optics' products. The BBP443 product is generated from the L3 RRS products using a quasi-analytical algorithm (Lee et al. 2002). he 'tur_tsm_chl' products include TUR, SPM and CHL. They are retrieved through the application of automated switching algorithms to the RRS spectra adapted to varying water conditions (Novoa et al. 2017). The GEOPHYSICAL product consists of the Chlorophyll-a concentration (CHL) retrieved via a multi-algorithm approach with optimized quality flagging (O'Reilly et al. 2019, Gons et al. 2005, Lavigne et al. 2021). Monthly products (P1M) are temporal aggregates of the daily L3 products. Daily products contain gaps in cloudy areas and where there is no overpass at the respective day. Aggregation collects the non-cloudy (and non-frozen) contributions to each pixel. Contributions are averaged per variable. While this does not guarantee data availability in all pixels in case of persistent clouds, it provides a more complete product compared to the sparsely filled daily products. The Monthly L4 products (P1M) are generally provided withing 4 days after the last acquisition date of the month. Daily gap filled L4 products (P1D) are generated using the DINEOF (Data Interpolating Empirical Orthogonal Functions) approach which reconstructs missing data in geophysical datasets by using a truncated Empirical Orthogonal Functions (EOF) basis in an iterative approach. DINEOF reconstructs missing data in a geophysical dataset by extracting the main patterns of temporal and spatial variability from the data. While originally designed for low resolution data products, recent research has resulted in the optimization of DINEOF to handle high resolution data provided by Sentinel-2 MSI, including cloud shadow detection (Alvera-Azcárate et al., 2021). These types of L4 products are generated and delivered one month after the respective period. '''Processing information:''' The HR-OC processing system is deployed on Creodias where Sentinel 2/MSI L1C data are available. The production control element is being hosted within the infrastructure of Brockmann Consult. The processing chain consists of: * Resampling to 60m and mosaic generation of the set of Sentinel-2 MSI L1C granules of a single overpass that cover a single UTM zone. * Application of a glint correction taking into account the detector viewing angles * Application of a coastal mask with 20km water + 20km land. The result is a L1C mosaic tile with data just in the coastal area optimized for compression. * Level 2 processing with pixel identification (IdePix), atmospheric correction (C2RCC and ACOLITE or iCOR), in-water processing and merging (HR-OC L2W processor). The result is a 60m product with the same extent as the L1C mosaic, with variables for optics, transparency, and geophysics, and with data filled in the water part of the coastal area. * invalid pixel identification takes into account corrupted (L1) pixels, clouds, cloud shadow, glint, dry-fallen intertidal flats, coastal mixed-pixels, sea ice, melting ice, floating vegetation, non-water objects, and bottom reflection. * Daily L3 aggregation merges all Level 2 mosaics of a day intersecting with a target tile. All valid water pixels are included in the 20km coastal stripes; all other values are set to NaN. There may be more than a single overpass a day, in particular in the northern regions. This step comprises resampling to the 100m target grid. * Monthly L4 aggregation combines all Level 3 products of a month and a single tile. The output is a set of 3 NetCDF datasets for (1) optics and (2) turbidity, suspended matter and chlorophyll concentration, respectively for the month. * Gap filling combines all daily products of a period and generates (partially) gap-filled daily products again. The output of gap filling are 2 datasets for (1) optics (BBP443 only) and (2) turbidity, suspended mattr and chlorophyll concentration per day. '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so-called ocean colour which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of the calibration and validation activities performed over this product can be found on the CMEMS web portal and in CMEMS-BGP_HR-QUID-009-201_to_212. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''Dataset names: ''' *cmems_obs_oc_med_bgc_tur_spm_chl_nrt_l4-hr-mosaic_P1M-v01 *cmems_obs_oc_med_bgc_optics_nrt_l4-hr-mosaic_P1M-v01 *cmems_obs_oc_med_bgc_tur_spm_chl_nrt_l4-hr-mosaic_P1D-v01 *cmems_obs_oc_med_bgc_optics_nrt_l4-hr-mosaic_P1D-v01 '''Files format:''' *netCDF-4, CF-1.7 *INSPIRE compliant." '''DOI (product) :''' https://doi.org/10.48670/moi-00110

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' We have derived an annual eutrophication and eutrophication indicator map for the North Atlantic Ocean using satellite-derived chlorophyll concentration. Using the satellite-derived chlorophyll products distributed in the regional North Atlantic CMEMS REP Ocean Colour dataset (OC- CCI), we derived P90 and P10 daily climatologies. The time period selected for the climatology was 1998-2017. For a given pixel, P90 and P10 were defined as dynamic thresholds such as 90% of the 1998-2017 chlorophyll values for that pixel were below the P90 value, and 10% of the chlorophyll values were below the P10 value. To minimise the effect of gaps in the data in the computation of these P90 and P10 climatological values, we imposed a threshold of 25% valid data for the daily climatology. For the 20-year 1998-2017 climatology this means that, for a given pixel and day of the year, at least 5 years must contain valid data for the resulting climatological value to be considered significant. Pixels where the minimum data requirements were met were not considered in further calculations. We compared every valid daily observation over 2020 with the corresponding daily climatology on a pixel-by-pixel basis, to determine if values were above the P90 threshold, below the P10 threshold or within the [P10, P90] range. Values above the P90 threshold or below the P10 were flagged as anomalous. The number of anomalous and total valid observations were stored during this process. We then calculated the percentage of valid anomalous observations (above/below the P90/P10 thresholds) for each pixel, to create percentile anomaly maps in terms of % days per year. Finally, we derived an annual indicator map for eutrophication levels: if 25% of the valid observations for a given pixel and year were above the P90 threshold, the pixel was flagged as eutrophic. Similarly, if 25% of the observations for a given pixel were below the P10 threshold, the pixel was flagged as oligotrophic. '''CONTEXT''' Eutrophication is the process by which an excess of nutrients – mainly phosphorus and nitrogen – in a water body leads to increased growth of plant material in an aquatic body. Anthropogenic activities, such as farming, agriculture, aquaculture and industry, are the main source of nutrient input in problem areas (Jickells, 1998; Schindler, 2006; Galloway et al., 2008). Eutrophication is an issue particularly in coastal regions and areas with restricted water flow, such as lakes and rivers (Howarth and Marino, 2006; Smith, 2003). The impact of eutrophication on aquatic ecosystems is well known: nutrient availability boosts plant growth – particularly algal blooms – resulting in a decrease in water quality (Anderson et al., 2002; Howarth et al.; 2000). This can, in turn, cause death by hypoxia of aquatic organisms (Breitburg et al., 2018), ultimately driving changes in community composition (Van Meerssche et al., 2019). Eutrophication has also been linked to changes in the pH (Cai et al., 2011, Wallace et al. 2014) and depletion of inorganic carbon in the aquatic environment (Balmer and Downing, 2011). Oligotrophication is the opposite of eutrophication, where reduction in some limiting resource leads to a decrease in photosynthesis by aquatic plants, reducing the capacity of the ecosystem to sustain the higher organisms in it. Eutrophication is one of the more long-lasting water quality problems in Europe (OSPAR ICG-EUT, 2017), and is on the forefront of most European Directives on water-protection. Efforts to reduce anthropogenically-induced pollution resulted in the implementation of the Water Framework Directive (WFD) in 2000. '''CMEMS KEY FINDINGS''' Some coastal and shelf waters, especially between 30 and 400N showed active oligotrophication flags for 2020, with some scattered offshore locations within the same latitudinal belt also showing oligotrophication. Eutrophication index is positive only for a small number of coastal locations just north of 40oN, and south of 30oN. In general, the indicator map showed very few areas with active eutrophication flags for 2019 and for 2020. The Third Integrated Report on the Eutrophication Status of the OSPAR Maritime Area (OSPAR ICG-EUT, 2017) reported an improvement from 2008 to 2017 in eutrophication status across offshore and outer coastal waters of the Greater North Sea, with a decrease in the size of coastal problem areas in Denmark, France, Germany, Ireland, Norway and the United Kingdom. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00195

'''This product has been archived''' '''DEFINITION''' Estimates of Ocean Heat Content (OHC) are obtained from integrated differences of the measured temperature and a climatology along a vertical profile in the ocean (von Schuckmann et al., 2018). The regional OHC values are then averaged from 60°S-60°N aiming i) to obtain the mean OHC as expressed in Joules per meter square (J/m2) to monitor the large-scale variability and change. ii) to monitor the amount of energy in the form of heat stored in the ocean (i.e. the change of OHC in time), expressed in Watt per square meter (W/m2). Ocean heat content is one of the six Global Climate Indicators recommended by the World Meterological Organisation for Sustainable Development Goal 13 implementation (WMO, 2017). '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the ocean shapes our perspectives for the future (von Schuckmann et al., 2020). Variations in OHC can induce changes in ocean stratification, currents, sea ice and ice shelfs (IPCC, 2019; 2021); they set time scales and dominate Earth system adjustments to climate variability and change (Hansen et al., 2011); they are a key player in ocean-atmosphere interactions and sea level change (WCRP, 2018) and they can impact marine ecosystems and human livelihoods (IPCC, 2019). '''CMEMS KEY FINDINGS''' Since the year 2005, the upper (0-700m) near-global (60°S-60°N) ocean warms at a rate of 0.6 ± 0.1 W/m2. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00234