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

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

'''DEFINITION''' The CMEMS MEDSEA_OMI_tempsal_extreme_var_temp_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Sea Surface Temperature (SST) from model data. Two different CMEMS products are used to compute the indicator: The Iberia-Biscay-Ireland Multi Year Product (MEDSEA_MULTIYEAR_PHY_006_004) and the Analysis product (MEDSEA_ANALYSISFORECAST_PHY_006_013). Two parameters have been considered for this OMI: * Map of the 99th mean percentile: It is obtained from the Multi Year Product, the annual 99th percentile is computed for each year of the product. The percentiles are temporally averaged over the whole period (1987-2019). * Anomaly of the 99th percentile in 2020: The 99th percentile of the year 2020 is computed from the Near Real Time product. The anomaly is obtained by subtracting the mean percentile from the 2020 percentile. This indicator is aimed at monitoring the extremes of sea surface temperature every year and at checking their variations in space. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This study of extreme variability was first applied to the sea level variable (Pérez Gómez et al 2016) and then extended to other essential variables, such as sea surface temperature and significant wave height (Pérez Gómez et al 2018 and Alvarez Fanjul et al., 2019). More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Alvarez Fanjul et al., 2019). '''CONTEXT''' The Sea Surface Temperature is one of the Essential Ocean Variables, hence the monitoring of this variable is of key importance, since its variations can affect the ocean circulation, marine ecosystems, and ocean-atmosphere exchange processes. As the oceans continuously interact with the atmosphere, trends of sea surface temperature can also have an effect on the global climate. In recent decades (from mid ‘80s) the Mediterranean Sea showed a trend of increasing temperatures (Ducrocq et al., 2016), which has been observed also by means of the CMEMS SST_MED_SST_L4_REP_OBSERVATIONS_010_021 satellite product and reported in the following CMEMS OMI: MEDSEA_OMI_TEMPSAL_sst_area_averaged_anomalies and MEDSEA_OMI_TEMPSAL_sst_trend. The Mediterranean Sea is a semi-enclosed sea characterized by an annual average surface temperature which varies horizontally from ~14°C in the Northwestern part of the basin to ~23°C in the Southeastern areas. Large-scale temperature variations in the upper layers are mainly related to the heat exchange with the atmosphere and surrounding oceanic regions. The Mediterranean Sea annual 99th percentile presents a significant interannual and multidecadal variability with a significant increase starting from the 80’s as shown in Marbà et al. (2015) which is also in good agreement with the multidecadal change of the mean SST reported in Mariotti et al. (2012). Moreover the spatial variability of the SST 99th percentile shows large differences at regional scale (Darmariaki et al., 2019; Pastor et al. 2018). '''CMEMS KEY FINDINGS''' The Mediterranean mean Sea Surface Temperature 99th percentile evaluated in the period 1987-2019 (upper panel) presents highest values (~ 28-30 °C) in the eastern Mediterranean-Levantine basin and along the Tunisian coasts especially in the area of the Gulf of Gabes, while the lowest (~ 23–25 °C) are found in the Gulf of Lyon (a deep water formation area), in the Alboran Sea (affected by incoming Atlantic waters) and the eastern part of the Aegean Sea (an upwelling region). These results are in agreement with previous findings in Darmariaki et al. (2019) and Pastor et al. (2018) and are consistent with the ones presented in CMEMS OSR3 (Alvarez Fanjul et al., 2019) for the period 1993-2016. The 2020 Sea Surface Temperature 99th percentile anomaly map (bottom panel) shows a general positive pattern up to +3°C in the North-West Mediterranean area while colder anomalies are visible in the Gulf of Lion and North Aegean Sea . This Ocean Monitoring Indicator confirms the continuous warming of the SST and in particular it shows that the year 2020 is characterized by an overall increase of the extreme Sea Surface Temperature values in almost the whole domain with respect to the reference period. This finding can be probably affected by the different dataset used to evaluate this anomaly map: the 2020 Sea Surface Temperature 99th percentile derived from the Near Real Time Analysis product compared to the mean (1987-2019) Sea Surface Temperature 99th percentile evaluated from the Reanalysis product which, among the others, is characterized by different atmospheric forcing). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00266

'''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 surface chlorophyll of the European region by merging the daily chlorophyll regional products over the Atlantic Ocean, the Baltic Sea, the Black Sea, and the Mediterranean Sea. Single chlorophyll daily images are the Case I – Case II products, which are produced accounting for bio-optical differences in these two water types. The mosaic is built using the following datasets: • dataset-oc-atl-chl-multi_cci-l3-chl_1km_daily-rt-v01 for the North Atlantic Ocean • dataset-oc-bal-chl-modis_a-l3-nn_1km_daily-rt-v01 for the Baltic Sea • dataset-oc-bs-chl-multi-l3-chl_1km_daily-rt-v02 for the Black Sea • dataset-oc-med-chl-multi-l3-chl_1km_daily-rt-v02 for the Mediterranean Sea. '''Processing information:''' All details about the processing can be found in relevant product description: *OCEANCOLOUR_ATL_CHL_L3_NRT_OBSERVATIONS_009_036 *OCEANCOLOUR_BAL_CHL_L3_NRT_OBSERVATIONS_009_049 *OCEANCOLOUR_BS_CHL_L3_NRT_OBSERVATIONS_009_044 *OCEANCOLOUR_MED_CHL_L3_NRT_OBSERVATIONS_009_040 '''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-eur-chl-multi-l3-chl_1km_daily-rt-v02 '''DOI (product) :''' https://doi.org/10.48670/moi-00095

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

'''Short description:''' Altimeter satellite along-track sea surface heights anomalies (SLA) computed with respect to a twenty-year [1993, 2012] mean. 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) that can be used to change the physical content for specific needs This product is processed by the DUACS multimission altimeter data processing system. It serves in near-real time the main operational oceanography and climate forecasting centers in Europe and worldwide. It processes data from all altimeter missions: Jason-3, Sentinel-3A, HY-2A, Saral/AltiKa, Cryosat-2, Jason-2, Jason-1, T/P, ENVISAT, GFO, ERS1/2. It provides a consistent and homogeneous catalogue of products for varied applications, both for near real time applications and offline studies. To produce maps of SLA (Sea Level Anomalies) in near-real time, the system exploits the most recent datasets available based on the enhanced OGDR+IGDR production. The system acquires and then synchronizes altimeter data and auxiliary data; each mission is homogenized using the same models and corrections. The Input Data Quality Control checks that the system uses the best altimeter data. The multi-mission cross-calibration process removes any residual orbit error, or long wavelength error (LWE), as well as large scale biases and discrepancies between various data flows; all altimeter fields are interpolated at crossover locations and dates. After a repeat-track analysis, a mean profile, which is peculiar to each mission, or a Mean Sea Surface (MSS) (when the orbit is non repetitive) is subtracted to compute sea level anomaly. The MSS is available via the Aviso+ dissemination (http://www.aviso.altimetry.fr/en/data/products/auxiliary-products/mss.html [http://www.aviso.altimetry.fr/en/data/products/auxiliary-products/mss.html]). Data are then cross validated, filtered from residual noise and small scale signals, and finally sub-sampled (sla_filtered variable). The ADT (Absolute Dynamic Topography, adt_filtered variable) can computed as follows: adt_filtered=sla_filtered+MDT where MDT. The Mean Dynamic Topography distributed by Aviso+ (http://www.aviso.altimetry.fr/en/data/products/auxiliary-products/mdt.html [http://www.aviso.altimetry.fr/en/data/products/auxiliary-products/mdt.html]). '''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] 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.

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

'''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (MY) product as distributed by CMEMS (OCEANCOLOUR_MED_BGC_L3_NRT_009_141). This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3-OLCI) Rrs spectra produced by CNR using an in-house processing chain, is obtained by means of the Mediterranean Ocean Colour regional algorithms: an updated version of the MedOC4 (Case 1 (off-shore) waters, Volpe et al., 2019, with new coefficients) and AD4 (Case 2 (coastal) waters, Berthon and Zibordi, 2004). The processing chain and the techniques used for algorithms merging are detailed in Colella et al. (2023). Monthly regional mean values are calculated by performing the average of 2D monthly mean (weighted by pixel area) over the region of interest. The deseasonalized time series is obtained by applying the X-11 seasonal adjustment methodology on the original time series as described in Colella et al. (2016), and then the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are subsequently applied to obtain the magnitude of trend. This OMI has been introduced since the 2nd issue of Ocean State Report in 2017. '''CONTEXT''' Phytoplankton and chlorophyll concentration as a proxy for phytoplankton respond rapidly to changes in environmental conditions, such as light, temperature, nutrients and mixing (Colella et al. 2016). The character of the response depends on the nature of the change drivers, and ranges from seasonal cycles to decadal oscillations (Basterretxea et al. 2018). Therefore, it is of critical importance to monitor chlorophyll concentration at multiple temporal and spatial scales, in order to be able to separate potential long-term climate signals from natural variability in the short term. In particular, phytoplankton in the Mediterranean Sea is known to respond to climate variability associated with the North Atlantic Oscillation (NAO) and El Niño Southern Oscillation (ENSO) (Basterretxea et al. 2018, Colella et al. 2016). '''KEY FINDINGS''' In the Mediterranean Sea, the average chlorophyll trend for the 1997–2024 period is slightly negative, at -0.77 ± 0.59% per year, reinforcing the findings of the previous releases. This result contrasts with the analysis by Sathyendranath et al. (2018), which reported increasing chlorophyll concentrations across all European seas. From around 2010–2011 onward, excluding the 2018–2019 period, a noticeable decline in chlorophyll levels is evident in the deseasonalized time series (green line) and in the observed maxima (grey line), particularly from 2015. This sustained decline over the past decade contributes to the overall negative trend observed in the Mediterranean Sea. '''DOI (product):''' https://doi.org/10.48670/moi-00259

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