'''DEFINITION''' The sea level ocean monitoring indicator is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and are also available in the CMEMS catalogue (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the North-West Shelf region is derived from the average of the gridded sea level maps weighted by the cosine of the latitude. The annual and semi-annual periodic signals are adjusted and the time series is low-pass filtered. The curve is corrected for the effect of the Glacial Isostatic Adjustment using the ICE5G-VM2 GIA model (Peltier, 2004). '''CONTEXT''' Mean sea level evolution has a direct impact on coastal areas and is a crucial index of climate change since it reflects both the amount of heat added in the ocean and the mass loss due to land ice melt (e.g. IPCC, 2013; Dieng et al., 2017). Long-term and inter-annual variations of the sea level are observed at global and regional scales. They are strongly related to the internal variability observed at basin scale and these variations can strongly affect population living in coastal areas. '''CMEMS KEY FINDINGS''' Using the latest reprocessed altimeter sea level products, it is possible to estimate the sea level rise in the North-West Shelf region since 01/1993 (see the proposed figure of the indicator for the updated trend value). The associated uncertainty is provided in a 90% confidence interval and only errors related to the altimeter observation system have been considered in the sea level trend uncertainty (Prandi et al., 2021). The uncertainty due to the sea level internal variability of the observed ocean is not included and should be considered additionally. The North-West Shelf region is a large, deep and open ocean area, which explains the relatively reduced inter annual variability observed in the sea level record compared to semi-enclosed basin with shallow bathymetry (Baltic and Black seas). '''DOI (product):''' https://doi.org/10.48670/moi-00271

'''DEFINITION''' The subsurface salinity trends have been derived from regional reanalysis and forecast modelling results of the CMEMS BAL MFC group for the Baltic Sea (product references BALTICSEA_REANALYSIS_PHY_003_011). The salinity trend has been obtained through a linear fit for each time series of horizontally averaged (13 °E - 31 °E and 53 °N - 66 °N; excluding the Skagerrak strait) annual salinity and at each depth level (Mulet et al., 2018). '''CONTEXT''' The Baltic Sea is a brackish semi-enclosed sea in North-Eastern Europe. The surface salinity varies horizontally from ~10 near the Danish Straits down to ~2 at the northernmost and easternmost sub-basins of the Baltic Sea. The halocline, a vertical layer with rapid changes of salinity with depth that separates the well-mixed surface layer from the weakly stratified layer below, is located at the depth range of 60-80 meters (Matthäus, 1984). The bottom layer salinity below the halocline depth varies from 15 in the south down to 3 in the northern Baltic Sea (Väli et al., 2013). The long-term salinity is determined by net precipitation and river discharge as well as saline water inflows from the North Sea (BACCII Author Team, 2015). Long-term salinity decrease may reduce the occurrence and biomass of the Fucus vesiculosus - Idotea balthica association/symbiotic aggregations (Kotta et al., 2019). Changes in salinity and oxygen content affect the survival of the Baltic cod eggs (Raudsepp et al, 2019; von Dewitz et al., 2018). '''CMEMS KEY FINDINGS''' The subsurface salinity over the 1993-2020 period shows no trend in the surface layer of 40-m. Then, the trend increases monotonically until the depth of 80-m, where it reaches the value of about 0.04 per year. Below that depth the salinity trend is uniform within the uncertainty limits. '''DOI (product):''' https://doi.org/10.48670/moi-00207

'''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_bal_bgc_geophy_nrt_l3-hr_P1D-v01 *cmems_obs_oc_bal_bgc_transp_nrt_l3-hr_P1D-v01 *cmems_obs_oc_bal_bgc_optics_nrt_l3-hr_P1D-v01 '''Files format:''' *netCDF-4, CF-1.7 *INSPIRE compliant. '''DOI (product) :''' https://doi.org/10.48670/moi-00079

'''DEFINITION''' The sea level ocean monitoring indicator is derived from the DUACS delayed-time (DT-2021 version) altimeter gridded maps of sea level anomalies based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and are also available in the CMEMS catalogue (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Black Sea is derived from the average of the gridded sea level maps weighted by the cosine of the latitude. The annual and semi-annual periodic signals are adjusted and the time series is low-pass filtered. The curve is corrected for the effect of the Glacial Isostatic Adjustment using the ICE5G-VM2 GIA model (Peltier, 2004). '''CONTEXT''' Mean sea level evolution has a direct impact on coastal areas and is a crucial index of climate change since it reflects both the amount of heat added in the ocean and the mass loss due to land ice melt (e.g. IPCC, 2013; Dieng et al., 2017). Long-term and inter-annual variations of the sea level are observed at global and regional scales. They are strongly related to the internal variability observed at basin scale and these variations can strongly affect population living in coastal areas. '''CMEMS KEY FINDINGS''' Using the latest reprocessed altimeter sea level products, it is possible to estimate the sea level rise in the Black Sea since 01/1993 (see the proposed figure of the indicator for the updated trend value). The associated uncertainty is provided in a 90% confidence interval and only errors related to the altimeter observation system have been considered in the sea level trend uncertainty (Prandi et al., 2021). The uncertainty due to the sea level internal variability of the observed ocean is not included and should be considered additionally. The Black Sea is a relatively small semi-enclosed basin with shallow bathymetry, which explains the high level of inter annual variability observed in the sea level record compared to large, deeper and open ocean areas (North West Shelf and Iberian-Biscay-Ireland regions). '''DOI (product):''' https://doi.org/10.48670/moi-00215

'''Short description:''' Mean Dynamic Topography that combines the global CNES-CLS18 MDT, the Black Sea CMEMS2020 MDT and the Med Sea CMEMS2020 MDT. It is an estimate of the mean over the 1993-2012 period of the sea surface height above geoid. This is consistent with the reference time period also used in the DUACS products '''DOI (product) :''' https://doi.org/10.48670/moi-00150

'''Short description:''' Baltic Sea ice concentration, extent, and classification daily time series since the winter 1981-1982 '''DOI (product) :''' https://doi.org/10.48670/moi-00131

'''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:''' 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_bal_bgc_geophy_nrt_l4-hr_P1M-v01 *cmems_obs_oc_bal_bgc_transp_nrt_l4-hr_P1M-v01 *cmems_obs_oc_bal_bgc_optics_nrt_l4-hr_P1M-v01 *cmems_obs_oc_bal_bgc_geophy_nrt_l4-hr_P1D-v01 *cmems_obs_oc_bal_bgc_transp_nrt_l4-hr_P1D-v01 *cmems_obs_oc_bal_bgc_optics_nrt_l4-hr_P1D-v01 '''Files format:''' *netCDF-4, CF-1.7 *INSPIRE compliant '''DOI (product) :''' https://doi.org/10.48670/moi-00080

'''DEFINITION''' This product includes the Baltic Sea satellite chlorophyll trend map from 1997 to 2020 based on regional chlorophyll reprocessed (REP) product as distributed by CMEMS OC-TAC which, in turn, result from the application of the regional chlorophyll algorithms over remote sensing reflectances (Rrs) provided by the Plymouth Marine Laboratory (PML) using the ESA Ocean Colour Climate Change Initiative processor (ESA OC-CCI, Sathyendranath et al., 2018a). The chlorophyll product is derived from a Multi Layer Perceptron neural-net (MLP) developed on field measurements collected within the BiOMaP program of JRC/EC (Zibordi et al., 2011). The algorithm is an ensemble of different MLPs that use Rrs at different wavelengths as input. The processing chain and the techniques used to develop the algorithm are detailed in Brando et al. (2021a; 2021b). 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 in the Baltic Sea depends on the nature of the change drivers, and ranges from seasonal cycles to decadal oscillations (Kahru and Elmgren 2014) and anthropogenic climate change. Eutrophication is one of the most important issue for the Baltic Sea (HELCOM, 2018), therefore 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''' The average Baltic Sea trend for the 1997-2020 period is 0.5% per year. A positive trend characterizes the central area of the basin from Northern Baltic Proper to Southern part, throughout Eastern and Western Gotland Basin. This result is in accordance to those of Sathyendranath et al. (2018b), that reveal an increasing trend in chlorophyll concentration in most of the European Seas. Weak negative trend is observable in the eastern sector of Gulf of Finland, Bothnian Bay and over the Gulf of Riga. Generally, along the coast of the basin the trend is no significant. Finally, in the 1997-2020 time window, the Bothnian Bay does not show a specific positive or negative trend, with percentage close to zero. '''DOI (product):''' https://doi.org/10.48670/moi-00198

'''Short description:''' For the Baltic Sea- The DMI Sea Surface Temperature L3S aims at providing daily multi-sensor supercollated data at 0.03deg. x 0.03deg. horizontal resolution, using satellite data from infra-red radiometers. Uses SST satellite products from these sensors: NOAA AVHRRs 7, 9, 11, 14, 16, 17, 18 , Envisat ATSR1, ATSR2 and AATSR. '''DOI (product) :''' https://doi.org/10.48670/moi-00154