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

'''Short description:''' For the Mediterranean Sea - the CNR diurnal sub-skin Sea Surface Temperature (SST) product provides daily gap-free (L4) maps of hourly mean sub-skin SST at 1/16° (0.0625°) horizontal resolution over the CMEMS Mediterranean Sea (MED) domain, by combining infrared satellite and model data (Marullo et al., 2014). The implementation of this product takes advantage of the consolidated operational SST processing chains that provide daily mean SST fields over the same basin (Buongiorno Nardelli et al., 2013). The sub-skin temperature is the temperature at the base of the thermal skin layer and it is equivalent to the foundation SST at night, but during daytime it can be significantly different under favorable (clear sky and low wind) diurnal warming conditions. The sub-skin SST L4 product is created by combining geostationary satellite observations aquired from SEVIRI and model data (used as first-guess) aquired from the CMEMS MED Monitoring Forecasting Center (MFC). This approach takes advantage of geostationary satellite observations as the input signal source to produce hourly gap-free SST fields using model analyses as first-guess. The resulting SST anomaly field (satellite-model) is free, or nearly free, of any diurnal cycle, thus allowing to interpolate SST anomalies using satellite data acquired at different times of the day (Marullo et al., 2014). '''DOI (product) :''' https://doi.org/10.48670/moi-00170

'''Short description:''' For the European Ocean, the L4 multi-sensor daily satellite product is a 2km horizontal resolution subskin sea surface temperature analysis. This SST analysis is run by Meteo France CMS and is built using the European Ocean L3S products originating from bias-corrected European Ocean L3C mono-sensor products at 0.02 degrees resolution. This analysis uses the analysis of the previous day at the same time as first guess field. '''DOI (product) :''' https://doi.org/10.48670/moi-00161

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

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

'''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 [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. '''DOI (product) :''' https://doi.org/10.48670/moi-00140

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

'''Short description:''' The High-Resolution Ocean Colour (HR-OC) Consortium (Brockmann Consult, Royal Belgian Institute of Natural Sciences, Flemish Institute for Technological Research) distributes Remote Sensing Reflectances (RRS, expressed in sr-1), Turbidity (TUR, expressed in FNU), Solid Particulate Matter Concentration (SPM, expressed in mg/l), spectral particulate backscattering (BBP, expressed in m-1) and chlorophyll-a concentration (CHL, expressed in µg/l) for the Sentinel 2/MSI sensor at 100m resolution for a 20km coastal zone. The products are delivered on a geographic lat-lon grid (EPSG:4326). To limit file size the products are provided in tiles of 600x800 km². RRS and BBP are delivered at nominal central bands of 443, 492, 560, 665, 704, 740, 783, 865 nm. The primary variable from which it is virtually possible to derive all the geophysical and transparency products is the spectral RRS. This, together with the spectral BBP, constitute the category of the 'optics' products. The spectral BBP product is generated from the RRS products using a quasi-analytical algorithm (Lee et al. 2002). The 'transparency' products include TUR and SPM). They are retrieved through the application of automated switching algorithms to the RRS spectra adapted to varying water conditions (Novoa et al. 2017). The GEOPHYSICAL product consists of the Chlorophyll-a concentration (CHL) retrieved via a multi-algorithm approach with optimized quality flagging (O'Reilly et al. 2019, Gons et al. 2005, Lavigne et al. 2021). The NRT products are generally provided withing 24 hours after end of the day.The RRS product is accompanied by a relative uncertainty estimate (unitless) derived by direct comparison of the products to corresponding fiducial reference measurements provided through the AERONET-OC network. The current day data temporal consistency is evaluated as Quality Index (QI) for TUR, SPM and CHL: 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:''' The HR-OC processing system is deployed on Creodias where Sentinel 2/MSI L1C data are available. The production control element is being hosted within the infrastructure of Brockmann Consult. The processing chain consists of: * Resampling to 60m and mosaic generation of the set of Sentinel-2 MSI L1C granules of a single overpass that cover a single UTM zone. * Application of a coastal mask with 20km water + 20km land. The result is a L1C mosaic tile with data just in the coastal area optimized for compression. * Level 2 processing with pixel identification (IdePix), atmospheric correction (C2RCC and ACOLITE or iCOR), in-water processing and merging (HR-OC L2W processor). The result is a 60m product with the same extent as the L1C mosaic, with variables for optics, transparency, and geophysics, and with data filled in the water part of the coastal area. * Daily L3 aggregation merges all Level 2 mosaics of a day intersecting with a target tile. All valid water pixels are included in the 20km coastal stripes; all other values are set to NaN. There may be more than a single overpass a day, in particular in the northern regions. The main contribution usually is the mosaic of the zone, but also adjacent mosaics may overlap. This step comprises resampling to the 100m target grid. * Monthly L4 aggregation combines all Level 3 products of a month and a single tile. The output is a set of 3 NetCDF datasets for optics, transparency, and geophysics respectively, for the tile and month. * Gap filling combines all daily products of a period and generates (partially) gap-filled daily products again. The output of gap filling are 3 datasets for optics (BBP443 only), transparency, and geophysics per day. '''Description of observation methods/instruments:''' Ocean colour technique exploits the emerging electromagnetic radiation from the sea surface in different wavelengths. The spectral variability of this signal defines the so-called ocean colour which is affected by the presence of phytoplankton. '''Quality / Accuracy / Calibration information:''' A detailed description of the calibration and validation activities performed over this product can be found on the CMEMS web portal and in CMEMS-BGP_HR-QUID-009-201to212. '''Suitability, Expected type of users / uses:''' This product is meant for use for educational purposes and for the managing of the marine safety, marine resources, marine and coastal environment and for climate and seasonal studies. '''Dataset names: ''' *cmems_obs_oc_ibi_bgc_geophy_nrt_l3-hr_P1D-v01 *cmems_obs_oc_ibi_bgc_transp_nrt_l3-hr_P1D-v01 *cmems_obs_oc_ibi_bgc_optics_nrt_l3-hr_P1D-v01 '''Files format:''' *netCDF-4, CF-1.7 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00109

'''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 (Faizal and Rafiuddin, 2011). 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_REANALYSIS_PHY_001_031): GLORYS, C-GLORS, ORAS5, FOAM Two observation based products: CORA (INSITU_GLO_TS_REP_OBSERVATIONS_013_001_b) and ARMOR3D (GLOBAL_REP_PHY_001_021). 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-2018 at rate of 1.5±0.2 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''' This product includes the Mediterranean Sea satellite chlorophyll trend map from 1997 to 2020 based on regional chlorophyll reprocessed (REP) product as distributed by CMEMS OC-TAC. This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3A-OLCI) (at 1 km resolution) 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. (2021). The trend map is obtained by applying Colella et al. (2016) methodology, where the Mann-Kendall test (Mann, 1945; Kendall, 1975) and Sens’s method (Sen, 1968) are applied on deseasonalized monthly time series, as obtained from the X-11 technique (see e. g. Pezzulli et al. 2005), to estimate, trend magnitude and its significance. The trend is expressed in % per year that represents the relative changes (i.e., percentage) corresponding to the dimensional trend [mg m-3 y-1] with respect to the reference climatology (1997-2014). Only significant trends (p < 0.05) are included. '''CONTEXT''' Phytoplankton are key actors in the carbon cycle and, as such, recognised as an Essential Climate Variable (ECV). 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). The Mediterranean Sea is an oligotrophic basin, where chlorophyll concentration decreases following a specific gradient from West to East (Colella et al. 2016). The highest concentrations are observed in coastal areas and at the river mouths, where the anthropogenic pressure and nutrient loads impact on the eutrophication regimes (Colella et al. 2016). The the use of long-term time series of consistent, well-calibrated, climate-quality data record is crucial for detecting eutrophication. Furthermore, chlorophyll analysis also demands the use of robust statistical temporal decomposition techniques, in order to separate the long-term signal from the seasonal component of the time series. '''CMEMS KEY FINDINGS''' Chlorophyll trend in the Mediterranean Sea, for the period 1997-2020, is negative over most of the basin. Positive trend areas are visible only in the southern part of the western Mediterranean basin, in the Gulf of Lion, Rhode Gyre and partially along the Croatian coast of the Adriatic Sea. On average the trend in the Mediterranean Sea is about -0.5% per year. Nevertheless, as shown by Salgado-Hernanz et al. (2019) in their analysis (related to 1998-2014 satellite observations), there is not a clear difference between western and eastern basins of the Mediterranean Sea. In the Ligurian Sea, the trend switch to negative values, differing from the positive regime observed in the trend maps of both Colella et al. (2016) and Salgado-Hernanz et al. (2019), referred, respectively, to 1998-2009 and 1998-2014 time period, respectively. The waters offshore the Po River mouth show weak negative trend values, partially differing from the markable negative regime observed in the 1998-2009 period (Colella et al., 2016), and definitely moving from the positive trend observed by Salgado-Hernanz et al. (2019). '''DOI (product):''' https://doi.org/10.48670/moi-00260