'''Short description:''' This product provides daily (nighttime), gap-free (Level-4, L4) maps of foundation Sea Surface Temperature (SST) - that is, the SST free from diurnal warming - over the Mediterranean Sea, at high (HR, 1/16°) and ultra-high (UHR, 1/100°) spatial resolutions, covering the period from 2008 to present. Each map represents nighttime SST values (centered at 00:00 UTC) and is produced by the Italian National Research Council – Institute of Marine Sciences (CNR-ISMAR). L4 maps are generated by selecting only the highest-quality SST observations from upstream Level-2 (L2) data acquired within a short local nighttime window, in order to minimize cloud contamination and avoid the effects of the diurnal cycle. The main L2 sources currently ingested include SLSTR from Sentinel-3A and -3B, VIIRS from NOAA-21, NOAA-20, and Suomi-NPP, AVHRR from Metop-B and -C, and SEVIRI. A two-step algorithm allows to interpolate SST data at high and ultra-high spatial resolution, applying statistical techniques (Buongiorno Nardelli et al., 2013; Buongiorno Nardelli et al., 2015). Additionally, from 2024 onwards, an improved first-guess field has been used in the generation of the MED UHR L4 data, enhancing the product's spatial resolution of SST features and the accuracy of SST gradients via machine learning techniques (Fanelli et al., 2024). '''DOI (product) :''' https://doi.org/10.48670/moi-00172

'''DEFINITION''' The global annual chlorophyll anomaly is computed by subtracting a reference climatology (1997-2014) from the annual chlorophyll mean, on a pixel-by-pixel basis and in log10 space. Both the annual mean and the climatology are computed employing ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al., 2018a) global products (i.e. using the standard OC-CCI chlorophyll algorithms, OCI) as distributed by CMEMS. '''CONTEXT''' Phytoplankton – and chlorophyll concentration as a proxy for phytoplankton – respond rapidly to changes in their physical environment. Some of those changes are seasonal and are determined by light and nutrient availability (Racault et al., 2012). By comparing annual mean values to a climatology, we effectively remove the seasonal signal, while retaining information on potential events during the year. Chlorophyll anomalies can be correlated to climate indexes in particular regions, such as the ENSO index in the equatorial Pacific (Behrenfeld et al. 2006; Racault et al., 2012) and the IOD index in the Indian Ocean (Brewin et al., 2012). It is important to study chlorophyll anomalies in consonance with sea surface temperature and sea level anomalies, as increases in chlorophyll are generally consistent with decreases in SST and sea level anomalies, suggesting an increase in mixing and vertical nutrient transport (von Schuckmann et al., 2016). '''CMEMS KEY FINDINGS''' The average global chlorophyll anomaly 2019 is -0.02 log10(mg m-3), with a maximum value of 1.7 log10(mg m-3) and a minimum value of -3.2 log10(mg m-3). That is to say that, in average, the annual 2019 mean value is slightly lower (96%) than the 1997-2014 climatological value. The positive signals reported in 2016 and 2017 (Sathyendranath et al., 2018b) in the southern Pacific Ocean could still be observed in the 2019 map, while the significant negative anomalies in the tropical waters of the northern Pacific Ocean were also detected to a lesser extent. Areas showing a change of anomaly sign from 2019 include the southern coast of Japan (no anomaly to positive) and the tropical Atlantic (anomalies close to zero for 2019). A marked increase in chlorophyll concentration was observed during 2019 in the Great Australian Bight, while negative anomalies became stronger in the Guatemala Basin and the region south of the Gulf of Guinea and, with values of chlorophyll reaching as low as 30% of the climatological value (anomaly < -0.5 log10(mg m-3)). The persistent positive anomalies in the higher latitudes of the North Atlantic (> 40°) match the cooling observed in the 2018 and previous years SST anomaly maps.

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The indicator of the Kuroshio extension phase variations is based on the standardized high frequency altimeter Eddy Kinetic Energy (EKE) averaged in the area 142-149°E and 32-37°N and computed from the DUACS (https://duacs.cls.fr) delayed-time (reprocessed version DT-2021, CMEMS SEALEVEL_GLO_PHY_L4_MY_008_047) and near real-time (CMEMS SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS_008_046) altimeter sea level gridded products. The change in the reprocessed version (previously DT-2018) and the extension of the mean value of the EKE (now 27 years, previously 20 years) induce some slight changes not impacting the general variability of the Kuroshio extension (correlation coefficient of 0.988 for the total period, 0.994 for the delayed time period only). '''CONTEXT''' The long-term mean and trends alone do not give a complete view of the likely changes in position of unstable western boundary current extensions (Kelly et al., 2010). The Kuroshio Extension is an eastward-flowing current in the subtropical western North Pacific after the Kuroshio separates from the coast of Japan at 35°N, 140°E. Being the extension of a wind-driven western boundary current, the Kuroshio Extension is characterized by a strong variability and is rich in large-amplitude meanders and energetic eddies (Niiler et al., 2003; Qiu, 2003, 2002). The Kuroshio Extension region has the largest sea surface height variability on sub-annual and decadal time scales in the extratropical North Pacific Ocean (Jayne et al., 2009; Qiu and Chen, 2010, 2005). Prediction and monitoring of the path of the Kuroshio are of huge importance for local economies as the position of the Kuroshio extension strongly determines the regions where phytoplankton and hence fish are located. '''CMEMS KEY FINDINGS''' The different states of the Kuroshio extension phase have been presented and validated by Bessières et al. (2013) and further reported by Drévillon et al. (2018) in the Copernicus Ocean State Report #2. Two rather different states of the Kuroshio extension are observed: an ‘elongated state’ (also called ‘strong state’) corresponding to a narrow strong steady jet, and a ‘contracted state’ (also called ‘weak state’) in which the jet is weaker and more unsteady, spreading on a wider latitudinal band. When the Kuroshio Extension jet is in a contracted (elongated) state, the upstream Kuroshio Extension path tends to become more (less) variable and regional eddy kinetic energy level tends to be higher (lower). In between these two opposite phases, the Kuroshio extension jet has many intermediate states of transition and presents either progressively weakening or strengthening trends. In 2018, the indicator reveals an elongated state followed by a weakening neutral phase since then. '''DOI (product):''' https://doi.org/10.48670/moi-00222

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR), in Rome operationally produces Level-4 product includes monthly averaged datasets of the diffuse attenuation coefficient of light at 490 nm (Kd490) for multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI at 300m of resolution) (at 1 km resolution) and Sentinel3A-OLCI observations (at 300m resolution). Kd490 is the diffuse attenuation coefficient of light at 490 nm, and is a measure of the turbidity of the water column, i.e., how visible light in the blue-green region of the spectrum penetrates the water column. It is directly related to the presence of absorbing and scattering matter in the water column and is estimated through the ratio between Rrs at 490 and 555 nm. For the multi-sensor dataset, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by the GOS. The QAA allows the inversion of the radiative transfer equations to compute the Inherent Optical Properties. Level-4 product includes monthly averages along with the standard deviation and the number of observations in the period of integration. '''Processing information:''' Multi-sensor products are constituted by MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS and Sentinel3A-OLCI. For consistency with NASA L2 dataset, BRDF correction was applied to Sentinel3A-OLCI prior to band shifting and multi sensor merging. Hence, the single sensor OLCI data set is also distributed after BRDF correction. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution (300m for Sentinel3A-OLCI) using cylindrical equirectangular projection. Afterwards, single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged with a technique aimed at smoothing the differences among different sensors. This technique is developed by The Global Ocean Satellite monitoring and marine ecosystem study group (GOS) of the Italian National Research Council (CNR, Rome). Then geophysical fields (i.e. chlorophyll, kd490, bbp, aph and adg) are estimated via state-of-the-art algorithms for better product quality. Time averages are computed on the delayed-time data. '''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. '''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 :''' *dataset-oc-med-opt-multi-l4-kd490_1km_monthly-rt-v02 *dataset-oc-med-opt-olci-l4-kd490_300m_monthly-rt '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00117

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Global Ocean- in-situ reprocessed Carbon observations. This product contains observations and gridded files from two up-to-date carbon and biogeochemistry community data products: Surface Ocean Carbon ATlas SOCATv2021 and GLobal Ocean Data Analysis Project GLODAPv2.2021. The SOCATv2021-OBS dataset contains >25 million observations of fugacity of CO2 of the surface global ocean from 1957 to early 2021. The quality control procedures are described in Bakker et al. (2016). These observations form the basis of the gridded products included in SOCATv2020-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.2021-OBS contains >1 million observations from individual seawater samples of temperature, salinity, oxygen, nutrients, dissolved inorganic carbon, total alkalinity and pH from 1972 to 2019. 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 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. '''DOI (product) :''' https://doi.org/10.48670/moi-00035

'''Short description:''' Altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 1Hz (~7km) sampling. It serves in delayed-time applications. This product is processed by the DUACS multimission altimeter data processing system. It processes data from all altimeter missions available (e.g. Sentinel-6A, Jason-3, Sentinel-3A, Sentinel-3B, Saral/AltiKa, Cryosat-2, Jason-1, Jason-2, Topex/Poseidon, ERS-1, ERS-2, Envisat, Geosat Follow-On, HY-2A, HY-2B, etc.). The system exploits the most recent datasets available based on the enhanced GDR/NTC production. All the missions are homogenized with respect to a reference mission. Part of the processing is fitted to the Global ocean. (see QUID document or http://duacs.cls.fr [http://duacs.cls.fr] pages for processing details). The product gives additional variables (e.g. Mean Dynamic Topography, Dynamic Atmospheric Correction, Ocean Tides, Long Wavelength Errors) that can be used to change the physical content for specific needs (see PUM document for details) '''Associated products''' A time invariant product https://resources.marine.copernicus.eu/product-detail/SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033/INFORMATION describing the noise level of along-track measurements is available. It is associated to the sla_filtered variable. It is a gridded product. One file is provided for the global ocean and those values must be applied for Arctic and Europe products. For Mediterranean and Black seas, one value is given in the QUID document. '''DOI (product)''': https://doi.org/10.48670/moi-00146

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' For the European Ocean - Sea Surface Temperature Mono-Sensor L3 Observations. One SST file per 24h per area and per sensor (bias corrected) closest to the original resolution: SLSTR-A, AMSR2, SEVIRI, AVHRR_METOP_B, AVHRR18_G, AVHRR_19L, MODIS_A, MODIS_T, VIIRS_NPP. One SST file per file window per area and per sensor (bias corrected) closest to the original resolution , while still manageable in terms volume over the processed area. '''Description of observation methods/instruments:''' The METOP_B derived SSTs are not bias corrected because METOP_B is used as the reference sensor for the correction method. '''DOI (product) :''' https://doi.org/10.48670/moi-00162

'''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''' Regional trends for the period 2005-2019 from the Copernicus Marine Service multi-ensemble approach show warming at rates ranging from the global mean average up to more than 8 W/m2 in some specific regions (e.g. northern hemisphere western boundary current regimes). There are specific regions where a negative trend is observed above noise at rates up to about -5 W/m2 such as in the subpolar North Atlantic, or the western tropical Pacific. These areas are characterized by strong year-to-year variability (Dubois et al., 2018; Capotondi et al., 2020). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00236

'''Short description: ''' For the Global Ocean - In-situ observation yearly delivery in delayed mode of Ocean surface currents. '''Detailed description: ''' The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for Ocean surface currents. The data are collected from the Surface Drifter Data Assembly Centre (SD-DAC at NOAA AOML) completed by European data provided by EUROGOOS regional systems and national systems by the regional INS TAC components. All surface drifters data have been processed to check for drogue loss. Drogued and undrogued drifting buoy surface ocean currents are provided with a drogue presence flag as well as a wind slippage correction for undrogued buoy. '''Processing information: ''' From the near real time INS TAC product validated on a daily and weekly basis for forecasting purposes, and from the SD-DAC quality controlled dataset a scientifically validated product is created . It s a """"reference product"""" updated on a yearly basis. This product has been processed using a method that checks for drogue loss. Altimeter and wind data have been used to extract the direct wind slippage from the total drifting buoy velocities. The obtained wind slippage values have then been analyzed to identify probable undrogued data among the drifting buoy velocities dataset. A simple procedure has then been applied to produce an updated dataset including a drogue presence flag as well as a wind slippage correction. '''Suitability, Expected type of users / uses: ''' The product is designed to be assimilated into or for validation purposes of operational models operated by ocean forecasting centers for reanalysis purposes or for research community. These users need data aggregated and quality controlled in a reliable and documented manner.

'''Short description:''' The Mean Dynamic Topography MDT-CMEMS_2024_EUR is an estimate of the mean over the 1993-2012 period of the sea surface height above geoid for the European Seas. This is consistent with the reference time period also used in the SSALTO DUACS products '''DOI (product) :''' https://doi.org/10.48670/mds-00337