'''Short description:''' This product provides daily (nighttime), merged multi-sensor (Level-3S, L3S) 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 and is produced by the Italian National Research Council – Institute of Marine Sciences (CNR-ISMAR). L3S 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. These L3S data serve as input to an optimal interpolation procedure used to generate gap-free Level-4 (L4) SST fields, as implemented in product 010_004 (Buongiorno Nardelli et al., 2013). '''DOI (product) :''' https://doi.org/10.48670/moi-00171

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

'''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 distributes Remote Sensing Reflectances (Rrs) and diffuse attenuation coefficient of light at 490 nm (kd490) data. These datasets derived from Rrs multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI) spectra at the state-of-the-art algorithms for multi-sensor merging. 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. Reprocessed (multi-year) products are consistent and homogeneous in terms of format, algorithms and processing software. Rrs is defined as the ratio of upwelling radiance and downwelling irradiance at any wavelength (412, 443, 490, 555, and 670 nm). Kd490 is defined as 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 within the water column. It is directly related to the presence of scattering particles in the water column and is estimated through the ratio between Rrs at 490 and 555 nm. Kd490 is achieved via Mediterranean regional algorithm developed by GOS on the basis of MedBiOp in situ dataset (Volpe et al., 2019). The current day data temporal consistency is evaluated as Quality Index (QI): QI=(CurrentDataPixel-ClimatologyDataPixel)/STDDataPixel where QI is the difference between current data and the relevant climatological field as a signed multiple of climatological standard deviations (STDDataPixel). Inherent Optical Properties (aph443, adg443 and bbp443 at 443nm) are derived via QAAv6 model. '''Processing information:''' Multi-sensor product is 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. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution 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 and kd490, bbp, aph and adg) are estimated via state-of-the-art algorithms for better product quality. The entire data set is consistent and processed in one-shot mode (with an unique software version and identical configurations). '''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-l3-rrs412_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs490_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs510_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs555_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-rrs670_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-kd490_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-bbp443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-adg443_1km_daily-rep-v02 * dataset-oc-med-opt-multi-l3-aph443_1km_daily-rep-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00116

'''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, distributes Level-4 product including the daily interpolated chlorophyll field with no data voids starting from the multi-sensor (MODIS-Aqua, NOAA-20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI), the monthly averaged chlorophyll concentration for the multi-sensor and climatological fields, all at 1 km resolution. Chlorophyll field are obtained by means of the Mediterranean regional algorithms: an updated version of the MedOC4 (Case 1 waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 waters, Berthon and Zibordi, 2004). Discrimination between the two water types is performed by comparing the satellite spectrum with the average water type spectral signature from in situ measurements for both water types. Reference insitu dataset is MedBiOp (Volpe et al., 2019) where pure Case II spectra are selected using a k-mean cluster analysis (Melin et al., 2015). Merging of Case I and Case II information is performed estimating the Mahalanobis distance between the observed and reference spectra and using it as weight for the final merged value. The interpolated gap-free Level-4 Chl concentration is estimated by means of a modified version of the DINEOF algorithm by GOS (Volpe et al., 2018). DINEOF is an iterative procedure in which EOF are used to reconstruct the entire field domain. As first guess, it uses the SeaWiFS-derived daily climatological values at missing pixels and satellite observations at valid pixels. Monthly Level-4 dataset is the time averages of the L3 fields and includes the standard deviation and the number of observations in the monthly period of integration. SeaWiFS daily climatology provides reference for the calculation of Quality Indices (QI) for Chl observations. '''Processing information:''' Multi-sensor product is 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. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution 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) are estimated via state-of-the-art algorithms for better product quality. Level-4 includes both monthly time averages and the daily-interpolated fields. Time averages are computed on the delayed-time data. The interpolated product starts from the L3 products at 1 km resolution. At the first iteration, DINEOF procedure uses, as first guess for each of the missing pixels the relative daily climatological pixel. A procedure to smooth out spurious spatial gradients is applied to the daily merged image (observation and climatology). From the second iteration, the procedure uses, as input for the next one, the field obtained by the EOF calculation, using only a number of modes: that is, at the second round, only the first two modes, at the third only the first three, and so on. At each iteration, the same smoothing procedure is applied between EOF output and initial observations. The procedure stops when the variance explained by the current EOF mode exceeds that of noise. The entire data set is consistent and processed in one-shot mode (with an unique software version and identical configurations). For the climatology Rrs data are derived from the latest SeaWiFS NASA reprocessing (R2018.0) and converted to chlorophyll concentrations with the same algorithm as the one used for other L4 and/or L3 products. '''Description of observation methods/instruments:''' Ocean color technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so called ocean color which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of both the cal/val and a more in depth view of the method is given in QUID-009-038to045-071-073-078-079-095-096.pdf over Copernicus 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-chl-multi-l4-chl_1km_monthly-rep-v02 *dataset-oc-med-chl-multi-l4-interp_1km_daily-rep *dataset-oc-med-chl-seawifs-l4-chl_1km_daily-climatology-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant. '''DOI (product) :''' https://doi.org/10.48670/moi-00114

'''DEFINITION:''' The regional 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 the regional products as distributed by CMEMS, derived by application of the regional chlorophyll algorithms over remote sensing reflectances (Rrs) provided by the ESA Ocean Colour Climate Change Initiative (ESA OC-CCI, Sathyendranath et al., 2018a). '''CONTEXT:''' Phytoplankton – and chlorophyll concentration as their proxy – respond rapidly to changes in their physical environment. In the North Atlantic region these changes present a distinct seasonality and are mostly determined by light and nutrient availability (González Taboada et al., 2014). By comparing annual mean values to a climatology, we effectively remove the seasonal signal at each grid point, while retaining information on potential events during the year (Gregg and Rousseaux, 2014). In particular, North Atlantic anomalies can then be correlated with oscillations in the Northern Hemisphere Temperature (Raitsos et al., 2014). Chlorophyll anomalies also provide information on the status of the North Atlantic oligotrophic gyre, where evidence of rapid gyre expansion has been found for the 1997-2012 period (Polovina et al. 2008, Aiken et al., 2017, Sathyendranath et al., 2018b). '''CMEMS KEY FINDINGS:''' The average chlorophyll anomaly in the North Atlantic is -0.02 log10(mg m-3), with a maximum value of 1.0 log10(mg m-3) and a minimum value of -1.0 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. A moderate increase in chlorophyll concentration was observed in 2019 over the Bay of Biscay and regions close to Iceland and Greenland, such as the Irminger Basin and the Denmark Strait. In particular, the annual average values for those areas are around 160% of the 1997-2014 average (anomalies > 0.2 log10(mg m-3)). While the significant negative anomalies reported for 2016-2017 (Sathyendranath et al., 2018c) in the area west of the Ireland and Scotland coasts continued to manifest, the Irish and North Seas returned to their normative regime during 2019, with anomalies close to zero. A change in the anomaly sign (positive to negative) was also detected for the West European Basin, with annual values as low as 60% of the 1997-2014 average. This reduction in chlorophyll might be matched with negative anomalies in sea level during the period, indicating a dominance of upwelling factors over stratification.

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-2014) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 700 m depth: OHC=∫_(z_1)^(z_2)ρ_0 c_p (T_yr-T_clim )dz [1] with a reference density of = 1030 kgm-3 and a specific heat capacity of cp = 3980 J kg-1 °C-1 (e.g. von Schuckmann et al., 2009). Time series of annual mean values area averaged ocean heat content is provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and is evaluated for topography deeper than 300m. '''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 oceans shape our perspectives for the future. The quality evaluation of MEDSEA_OMI_OHC_area_averaged_anomalies is based on the “multi-product” approach as introduced in the second issue of the Ocean State Report (von Schuckmann et al., 2018), and following the MyOcean’s experience (Masina et al., 2017). Six global products and a regional (Mediterranean Sea) product have been used to build an ensemble mean, and its associated ensemble spread. The reference products are: • The Mediterranean Sea Reanalysis at 1/24 degree horizontal resolution (MEDSEA_MULTIYEAR_PHY_006_004, DOI: https://doi.org/10.25423/CMCC/MEDSEA_MULTIYEAR_PHY_006_004_E3R1, Escudier et al., 2020) • Four global reanalyses at 1/4 degree horizontal resolution (GLOBAL_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The ensemble mean ocean heat content anomaly time series over the Mediterranean Sea shows a continuous increase in the period 1993-2022 at rate of 1.38±0.08 W/m2 in the upper 700m. After 2005 the rate has clearly increased with respect the previous decade, in agreement with Iona et al. (2018). '''DOI (product):''' https://doi.org/10.48670/moi-00261

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''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 near-real 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, HY-2B). The system exploits the most recent datasets available based on the enhanced OGDR/NRT+IGDR/STC production. All the missions are homogenized with respect to a reference mission. Part of the processing is fitted to the European Sea area. (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 http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_NOISE_L4_NRT_OBSERVATIONS_008_032 [http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=SEALEVEL_GLO_PHY_NOISE_L4_STATIC_008_033] 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-00140

'''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 distributes reprocessed surface chlorophyll concentration (Chl) and phytoplankton functional types (PFT). Input Rrs multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, Sentinel3A-OLCI) spectra at the state-of-the-art algorithms for multi-sensor merging. Single sensor Rrs fields are band-shifted, over the SeaWiFS native bands (using the QAAv6 model, Lee et al., 2002) and merged. Reprocessed (multi-year) products are consistent and homogeneous in terms of format, algorithms and processing software. Chl is obtained by means of the Mediterranean regional algorithms: an updated version of the MedOC4 (Volpe et al., 2019) and AD4 (Berthon and Zibordi, 2004). Discrimination between the two water types is performed by comparing the satellite spectrum with the average spectrum from in situ measurements. Reference insitu dataset is MedBiOp (Volpe et al., 2019) where Case II spectra are selected with a k-mean cluster analysis (Melin et al., 2015). Merging of Case I and Case II information is performed estimating the Mahalanobis distance between observed and reference spectra and using it as weight for the final value. The PFT provides estimates of Chl concentration of 9 phytoplankton groups: Micro, Nano, Pico, Diato, Dino, Crypto, Hapto, Green and Prokar. Micro consists of Diato and Dino, Nano includes Crypto and Hapto and Pico is referred to Green and Prokar with the adjustment of Brewin et al. (2010) in the ultra-oligotrophic water for Pico and Nano. These classes are estimated via empirical regional functions, correlating Chl concentration with each in-situ PFT fraction computed by a regional diagnostic pigment analysis (Di Cicco et al. 2017). '''Processing information:''' Multi-sensor product is 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. Single sensor NASA Level-2 data are destriped and then all Level-2 data are remapped at 1 km spatial resolution 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 and kd490) are estimated via state-of-the-art algorithms for better product quality. The entire data set is consistent and processed in one-shot mode (with an unique software version and identical configurations). '''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-chl-multi-l3-chl_1km_daily-rep-v02 * dataset-oc-med-pft-multi-l3-pft_1km_daily-rep-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00112

'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''Short description:''' Experimental altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean with a 5Hz (~1.3km) sampling. All the missions are homogenized with respect to a reference mission (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 Atmosphic Correction, Ocean Tides, Long Wavelength Errors, Internal tide, …) that can be used to change the physical content for specific needs This product was generated as experimental products in a CNES R&D context. It was processed by the DUACS multimission altimeter data processing system. '''DOI (product) :''' https://doi.org/10.48670/moi-00137

'''DEFINITION''' Time mean meridional Eulerian streamfunctions are computed using the velocity field estimate provided by the Copernicus Marine Mediterranean Sea reanalysis over the period from 1987 to the year preceding the current one [-1Y], operationally extended yearly. The Eulerian meridional streamfunction is evaluated by integrating meridional velocity daily data first in a vertical direction, then in a meridional direction, and finally averaging over the reanalysis period. The Mediterranean overturning indices are derived for the eastern and western Mediterranean Sea by computing the annual streamfunction in the two areas separated by the Strait of Sicily around 36.5°N, and then considering the associated maxima. In each case a geographical constraint focused the computation on the main region of interest. For the western index, we focused on deep-water formation regions, thus excluding both the effect of shallow physical processes and the Gibraltar net inflow. For the eastern index, we investigate the Levantine and Cretan areas corresponding to the strongest meridional overturning cell locations, thus only a zonal constraint is defined. Time series of annual mean values is provided for the Mediterranean Sea using the Mediterranean 1/24o eddy resolving reanalysis (Escudier et al., 2020, 2021). More details can be found in the Copernicus Marine Ocean State Report issue 4 (OSR4, von Schuckmann et al., 2020) Section 2.4 (Lyubartsev et al., 2020) and in the QUID. '''CONTEXT''' The western and eastern Mediterranean clockwise meridional overturning circulation is connected to deep-water formation processes. The Mediterranean Sea 1/24o eddy resolving reanalysis (MEDSEA_MULTIYEAR_PHY_006_004, Escudier et al., 2020, 2021) is used to show the interannual variability of the Meridional Overturning Index. Details on the product are delivered in the PUM and QUID of this OMI. The Mediterranean Meridional Overturning Index is defined here as the maxima of the clockwise cells in the eastern and western Mediterranean Sea and is associated with deep and intermediate water mass formation processes that occur in specific areas of the basin: Gulf of Lion, Southern Adriatic Sea, Cretan Sea and Rhodes Gyre (Pinardi et al., 2015). As in the global ocean, the overturning circulation of the western and eastern Mediterranean are paramount to determine the stratification of the basins (Cessi, 2019). In turn, the stratification and deep water formation mediate the exchange of oxygen and other tracers between the surface and the deep ocean (e.g., Johnson et al., 2009; Yoon et al., 2018). In this sense, the overturning indices are potential gauges of the ecosystem health of the Mediterranean Sea, and in particular they could instruct early warning indices for the Mediterranean Sea to support the Sustainable Development Goal (SDG) 13 Target 13.3. '''CMEMS KEY FINDINGS''' The western and eastern Mediterranean overturning indices (WMOI and EMOI) are synthetic indices of changes in the thermohaline properties of the Mediterranean basin related to changes in the main drivers of the basin scale circulation. The western sub-basin clockwise overturning circulation is associated with the deep-water formation area of the Gulf of Lion, while the eastern clockwise meridional overturning circulation is composed of multiple cells associated with different intermediate and deep-water sources in the Levantine, Aegean, and Adriatic Seas. On average, the EMOI shows higher values than the WMOI indicating a more vigorous overturning circulation in eastern Mediterranean. The difference is mostly related to the occurrence of the eastern Mediterranean transient (EMT) climatic event, and linked to a peak of the EMOI in 1992 (Roether et al. 1996, 2014, Gertman et al. 2006). In 1999, the difference between the two indices started to decrease because EMT water masses reached the Sicily Strait flowing into the western Mediterranean Sea (Schroeder et al., 2016). The western peak in 2006 is discussed to be linked to anomalous deep-water formation during the Western Mediterranean Transition (Smith, 2008; Schroeder et al., 2016). Thus, the WMOI and EMOI indices are a useful tool for long-term climate monitoring of overturning changes in the Mediterranean Sea. '''DOI (product):''' https://doi.org/10.48670/mds-00317