'''DEFINITION''' The Mediterranean water mass formation rates are evaluated in 4 areas as defined in the Ocean State Report issue 2 (OSR2, von Schuckmann et al., 2018) section 3.4 (Simoncelli and Pinardi, 2018): (1) the Gulf of Lions for the Western Mediterranean Deep Waters (WMDW); (2) the Southern Adriatic Sea Pit for the Eastern Mediterranean Deep Waters (EMDW); (3) the Cretan Sea for Cretan Intermediate Waters (CIW) and Cretan Deep Waters (CDW); (4) the Rhodes Gyre, the area of formation of the so-called Levantine Intermediate Waters (LIW) and Levantine Deep Waters (LDW). Annual water mass formation rates have been computed using daily mixed layer depth estimates (density criteria Δσ = 0.01 kg/m3, 10 m reference level) considering the annual maximum volume of water above mixed layer depth with potential density within or higher the specific thresholds specified in Table 1 then divided by seconds per year. Annual mean values are provided using the Mediterranean 1/24o eddy resolving reanalysis (Escudier et al. 2020, 2021). Time spans from 1987 to the year preceding the current one [-1Y], operationally extended yearly. '''CONTEXT''' The formation of intermediate and deep water masses is one of the most important processes occurring in the Mediterranean Sea, being a component of its general overturning circulation. This circulation varies at interannual and multidecadal time scales and it is composed of an upper zonal cell (Zonal Overturning Circulation) and two main meridional cells in the western and eastern Mediterranean (Pinardi and Masetti 2000). The objective is to monitor the main water mass formation events using the eddy resolving Mediterranean Sea Reanalysis (MEDSEA_MULTIYEAR_PHY_006_004, Escudier et al. 2020, 2021) and considering Pinardi et al. (2015) and Simoncelli and Pinardi (2018) as references for the methodology. The Mediterranean Sea Reanalysis can reproduce both Eastern Mediterranean Transient and Western Mediterranean Transition phenomena and catches the principal water mass formation events reported in the literature. This will permit constant monitoring of the open ocean deep convection process in the Mediterranean Sea and a better understanding of the multiple drivers of the general overturning circulation at interannual and multidecadal time scales. Deep and intermediate water formation events reveal themselves by a deep mixed layer depth distribution in four Mediterranean areas: Gulf of Lions, Southern Adriatic Sea Pit, Cretan Sea and Rhodes Gyre. '''KEY FINDINGS''' The Western Mediterranean Deep Water (WMDW) formation events in the Gulf of Lion appear to be larger after 1999 consistently with Schroeder et al. (2006, 2008) related to the Eastern Mediterranean Transient event. This modification of WMDW after 2005 has been called Western Mediterranean Transition. WMDW formation events are consistent with Somot et al. (2016) and the event in 2009 is also reported in Houpert et al. (2016). The Eastern Mediterranean Deep Water (EMDW) formation in the Southern Adriatic Pit region displays a period of water mass formation between 1988 and 1993, in agreement with Pinardi et al. (2015), in 1996, 1999 and 2000 as documented by Manca et al. (2002). Weak deep water formation in winter 2006 is confirmed by observations in Vilibić and Šantić (2008). An intense deep water formation event is detected in 2012-2013 (Gačić et al., 2014). Last years are characterized by large events starting from 2017 (Mihanovic et al., 2021). Cretan Intermediate Water formation rates present larger peaks between 1989 and 1993 with the ones in 1992 and 1993 composing the Eastern Mediterranean Transient phenomena. The Cretan Deep Water formed in 1992 and 1993 is characterized by the highest densities of the entire period in accordance with Velaoras et al. (2014). The Levantine Deep Water formation rate in the Rhode Gyre region presents the largest values between 1992 and 1993 in agreement with Kontoyiannis et al. (1999). '''DOI (product):''' https://doi.org/10.48670/mds-00318

'''Short description:''' The Reprocessed (REP) Mediterranean (MED) dataset provides a stable and consistent long-term Sea Surface Temperature (SST) time series over the Mediterranean Sea (and the adjacent North Atlantic box) developed for climate applications. This product consists of daily (nighttime), merged multi-sensor (L3S), satellite-based estimates of the foundation SST (namely, the temperature free, or nearly-free, of any diurnal cycle) at 0.05° resolution grid covering the period from 1st January 1981 to present (approximately one month before real time). The MED-REP-L3S product is built from a consistent reprocessing of the collated level-3 (merged single-sensor, L3C) climate data record (CDR) v.3.0, provided by the ESA Climate Change Initiative (CCI) and covering the period up to 2021, and its interim extension (ICDR) that allows the regular temporal extension for 2022 onwards. '''DOI (product) :''' https://doi.org/10.48670/moi-00314

'''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 Reprocessed (REP) Mediterranean (MED) dataset provides a stable and consistent long-term Sea Surface Temperature (SST) time series over the Mediterranean Sea (and the adjacent North Atlantic box) developed for climate applications. This product consists of daily (nighttime), optimally interpolated (L4), satellite-based estimates of the foundation SST (namely, the temperature free, or nearly-free, of any diurnal cycle) at 0.05° resolution grid covering the period from 1st January 1981 to present (approximately one month before real time). The MED-REP-L4 product is built from a consistent reprocessing of the collated level-3 (merged single-sensor, L3C) climate data record (CDR) v.3.0, provided by the ESA Climate Change Initiative (CCI) and covering the period up to 2021, and its interim extension (ICDR) that allows the regular temporal extension for 2022 onwards. '''DOI (product) :''' https://doi.org/10.48670/moi-00173

'''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

'''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

'''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:''' 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 at 300m of resolution) (at 1 km resolution) and the monthly averaged chlorophyll concentration for the multi-sensor (at 1 km resolution) and Sentinel-OLCI Level-3 (at 300m 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. The other 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. '''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. 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. '''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 pre+D2sence 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-l4-chl_1km_monthly-rt-v02 *dataset-oc-med-chl-multi-l4-interp_1km_daily-rt-v02 *dataset-oc-med-chl-olci-l4-chl_300m_monthly-rt-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00113

'''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 surface chlorophyll concentration (mg m-3) derived from multi-sensor (MODIS-AQUA, NOAA20-VIIRS, NPP-VIIRS, and Sentinel3A-OLCI at 300m of resolution) (at 1 km resolution) and Sentinel3A-OLCI (at high resolution, 300m) Rrs spectra. Chlorophyll datasets are obtained by means of the Mediterranean Ocean Colour 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 at pixel-by-pixel level 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 1 and Case 2 information is performed estimating the Mahalanobis distance between the observed and reference spectra and using it as weight for the final merged value. This product identifies the average chlorophyll content of the surface layer as defined by the first optical depth (roughly one fifth of the euphotic depth). For multi-sensor observations, 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. 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). '''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 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. '''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-rt-v02 *dataset-oc-med-chl-olci-l3-chl_300m_daily-rt-v02 '''Files format:''' *CF-1.4 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00111