'''DEFINITION''' The CMEMS IBI_OMI_seastate_extreme_var_swh_mean_and_anomaly OMI indicator is based on the computation of the annual 99th percentile of Significant Wave Height (SWH) from model data. Two different CMEMS products are used to compute the indicator: The Iberia-Biscay-Ireland Multi Year Product (IBI_MULTIYEAR_WAV_005_006) and the Analysis product (IBI_ANALYSIS_FORECAST_WAV_005_005). 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 in the whole period (1993-2021). • Anomaly of the 99th percentile in 2022: The 99th percentile of the year 2022 is computed from the Analysis product. The anomaly is obtained by subtracting the mean percentile to the percentile in 2022. This indicator is aimed at monitoring the extremes of annual significant wave height and evaluate the spatio-temporal variability. The use of percentiles instead of annual maxima, makes this extremes study less affected by individual data. This approach was first successfully applied to 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 Álvarez-Fanjul et al., 2019). Further details and in-depth scientific evaluation can be found in the CMEMS Ocean State report (Álvarez- Fanjul et al., 2019). '''CONTEXT''' The sea state and its related spatio-temporal variability affect dramatically maritime activities and the physical connectivity between offshore waters and coastal ecosystems, impacting therefore on the biodiversity of marine protected areas (González-Marco et al., 2008; Savina et al., 2003; Hewitt, 2003). Over the last decades, significant attention has been devoted to extreme wave height events since their destructive effects in both the shoreline environment and human infrastructures have prompted a wide range of adaptation strategies to deal with natural hazards in coastal areas (Hansom et al., 2019). Complementarily, there is also an emerging question about the role of anthropogenic global climate change on present and future extreme wave conditions. The Iberia-Biscay-Ireland region, which covers the North-East Atlantic Ocean from Canary Islands to Ireland, is characterized by two different sea state wave climate regions: whereas the northern half, impacted by the North Atlantic subpolar front, is of one of the world’s greatest wave generating regions (Mørk et al., 2010; Folley, 2017), the southern half, located at subtropical latitudes, is by contrast influenced by persistent trade winds and thus by constant and moderate wave regimes. The North Atlantic Oscillation (NAO), which refers to changes in the atmospheric sea level pressure difference between the Azores and Iceland, is a significant driver of wave climate variability in the Northern Hemisphere. The influence of North Atlantic Oscillation on waves along the Atlantic coast of Europe is particularly strong in and has a major impact on northern latitudes wintertime (Martínez-Asensio et al. 2016; Bacon and Carter, 1991; Bouws et al., 1996; Bauer, 2001; Wolf et al., 2002; Gleeson et al., 2017). Swings in the North Atlantic Oscillation index produce changes in the storms track and subsequently in the wind speed and direction over the Atlantic that alter the wave regime. When North Atlantic Oscillation index is in its positive phase, storms usually track northeast of Europe and enhanced westerly winds induce higher than average waves in the northernmost Atlantic Ocean. Conversely, in the negative North Atlantic Oscillation phase, the track of the storms is more zonal and south than usual, with trade winds (mid latitude westerlies) being slower and producing higher than average waves in southern latitudes (Marshall et al., 2001; Wolf et al., 2002; Wolf and Woolf, 2006). Additionally a variety of previous studies have uniquevocally determined the relationship between the sea state variability in the IBI region and other atmospheric climate modes such as the East Atlantic pattern, the Arctic Oscillation, the East Atlantic Western Russian pattern and the Scandinavian pattern (Izaguirre et al., 2011, Martínez-Asensio et al., 2016). In this context, long‐term statistical analysis of reanalyzed model data is mandatory not only to disentangle other driving agents of wave climate but also to attempt inferring any potential trend in the number and/or intensity of extreme wave events in coastal areas with subsequent socio-economic and environmental consequences. '''CMEMS KEY FINDINGS''' The climatic mean of 99th percentile (1993-2021) reveals a north-south gradient of Significant Wave Height with the highest values in northern latitudes (above 8m) and lowest values (2-3 m) detected southeastward of Canary Islands, in the seas between Canary Islands and the African Continental Shelf. This north-south pattern is the result of the two climatic conditions prevailing in the region and previously described. The 99th percentile anomalies in 2023 show that during this period, the central latitudes of the domain (between 37 ºN and 50 ºN) were affected by extreme wave events that exceeded up to twice the standard deviation of the anomalies. These events impacted not only the open waters of the Northeastern Atlantic but also European coastal areas such as the west coast of Portugal, the Spanish Atlantic coast, and the French coast, including the English Channel. Additionally, the impact of significant wave extremes exceeding twice the standard deviation of anomalies was detected in the Mediterranean region of the Balearic Sea and the Algerian Basin. This pattern is commonly associated with the impact of intense Tramontana winds originating from storms that cross the Iberian Peninsula from the Gulf of Biscay. '''Figure caption''' Iberia-Biscay-Ireland Significant Wave Height extreme variability: Map of the 99th mean percentile computed from the Multi Year Product (left panel) and anomaly of the 99th percentile in 2022 computed from the Analysis product (right panel). Transparent grey areas (if any) represent regions where anomaly exceeds the climatic standard deviation (light grey) and twice the climatic standard deviation (dark grey). '''DOI (product):''' https://doi.org/10.48670/moi-00249

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

'''Short description:''' For the NWS/IBI Ocean- Sea Surface Temperature L3 Observations . This product provides daily foundation sea surface temperature from multiple satellite sources. The data are intercalibrated. This product consists in a fusion of sea surface temperature observations from multiple satellite sensors, daily, over a 0.05° resolution grid. It includes observations by polar orbiting from the ESA CCI / C3S archive . The L3S SST data are produced selecting only the highest quality input data from input L2P/L3P images within a strict temporal window (local nightime), to avoid diurnal cycle and cloud contamination. The observations of each sensor are intercalibrated prior to merging using a bias correction based on a multi-sensor median reference correcting the large-scale cross-sensor biases. '''DOI (product) :''' https://doi.org/10.48670/moi-00311

'''DEFINITION''' The OMI_EXTREME_SL_IBI_slev_mean_and_anomaly_obs indicator is based on the computation of the 99th and the 1st percentiles from in situ data (observations). It is computed for the variable sea level measured by tide gauges along the coast. The use of percentiles instead of annual maximum and minimum values, makes this extremes study less affected by individual data measurement errors. The annual percentiles referred to annual mean sea level are temporally averaged and their spatial evolution is displayed in the dataset omi_extreme_sl_ibi_slev_mean_and_anomaly_obs, jointly with the anomaly in the target year. This study of extreme variability was first applied to sea level variable (Pérez Gómez et al 2016) and then extended to other essential variables, sea surface temperature and significant wave height (Pérez Gómez et al 2018). '''CONTEXT''' Sea level (SLEV) is one of the Essential Ocean Variables most affected by climate change. Global mean sea level rise has accelerated since the 1990’s (Abram et al., 2019, Legeais et al., 2020), due to the increase of ocean temperature and mass volume caused by land ice melting (WCRP, 2018). Basin scale oceanographic and meteorological features lead to regional variations of this trend that combined with changes in the frequency and intensity of storms could also rise extreme sea levels up to one meter by the end of the century (Vousdoukas et al., 2020, Tebaldi et al., 2021). This will significantly increase coastal vulnerability to storms, with important consequences on the extent of flooding events, coastal erosion and damage to infrastructures caused by waves (Boumis et al., 2023). The increase in extreme sea levels over recent decades is, therefore, primarily due to the rise in mean sea level. Note, however, that the methodology used to compute this OMI removes the annual 50th percentile, thereby discarding the mean sea level trend to isolate changes in storminess. The Iberian Biscay Ireland region shows positive sea level trend modulated by decadal-to-multidecadal variations driven by ocean dynamics and superposed to the long-term trend (Chafik et al., 2019). '''COPERNICUS MARINE SERVICE KEY FINDINGS''' The completeness index criteria is fulfilled by 57 stations in 2021, two more than those available in 2021 (55), recently added to the multi-year product INSITU_GLO_PHY_SSH_DISCRETE_MY_013_053. The mean 99th percentiles reflect the great tide spatial variability around the UK and the north of France. Minimum values are observed in the Irish eastern coast (e.g.: 0.66 m above mean sea level in Arklow Harbour) and the Canary Islands (e.g.: 0.93 and 0.96 m above mean sea level in Gomera and Hierro, respectively). Maximum values are observed in the Bristol and English Channels (e.g.: 6.26, 5.58 and 5.17 m above mean sea level in Newport, St. Malo and St. Helier, respectively). The annual 99th percentiles standard deviation reflects the south-north increase of storminess, ranging between 1-2 cm in the Canary Islands to 12 cm in Newport (Bristol Channel). Although less pronounced and general than in 2021, negative or close to zero anomalies of 2022 99th percentile still prevail throughout the region this year reaching up to -14 cm in St.Helier (Jersey Island, Channel Islands), or -12 cm in St. Malo. Positive anomalies of 2022 99th percentile are found in the northern part of the region (Irish eastern coast and west Scotland coast) and at a couple of stations in Southern England, with values reaching 9 cm in Bangor (Northern Ireland) and 6 cm in Portsmouth (South England). '''DOI (product):''' https://doi.org/10.48670/moi-00253

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''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

'''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:''' Mean Dynamic Topography that combines the global CNES-CLS-2022 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

'''DEFINITION''' The sea level ocean monitoring indicator has been presented in the Copernicus Ocean State Report #8. The ocean monitoring indicator on regional mean sea level is derived from the DUACS delayed-time (DT-2024 version, “my” (multi-year) dataset used when available, “myint” (multi-year interim) used after) sea level anomaly maps from satellite altimetry based on a stable number of altimeters (two) in the satellite constellation. These products are distributed by the Copernicus Climate Change Service and the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The time series of area averaged anomalies correspond to the area average of the maps in the Irish-Biscay-Iberian (IBI) Sea weighted by the cosine of the latitude (to consider the changing area in each grid with latitude) and by the proportion of ocean in each grid (to consider the coastal areas). The time series are corrected from regional mean GIA correction (weighted GIA mean of a 27 ensemble model following Spada et Melini, 2019). The time series are adjusted for seasonal annual and semi-annual signals and low-pass filtered at 6 months. Then, the trends/accelerations are estimated on the time series using ordinary least square fit.The trend uncertainty is provided in a 90% confidence interval. It is calculated as the weighted mean uncertainties in the region from Prandi et al., 2021. This estimate only considers errors related to the altimeter observation system (i.e., orbit determination errors, geophysical correction errors and inter-mission bias correction errors). The presence of the interannual signal can strongly influence the trend estimation considering to the altimeter period considered (Wang et al., 2021; Cazenave et al., 2014). The uncertainty linked to this effect is not considered. ""CONTEXT "" Change in mean sea level is an essential indicator of our evolving climate, as it reflects both the thermal expansion of the ocean in response to its warming and the increase in ocean mass due to the melting of ice sheets and glaciers (WCRP Global Sea Level Budget Group, 2018). At regional scale, sea level does not change homogenously. It is influenced by various other processes, with different spatial and temporal scales, such as local ocean dynamic, atmospheric forcing, Earth gravity and vertical land motion changes (IPCC WGI, 2021). The adverse effects of floods, storms and tropical cyclones, and the resulting losses and damage, have increased as a result of rising sea levels, increasing people and infrastructure vulnerability and food security risks, particularly in low-lying areas and island states (IPCC, 2022a). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022b). In IBI region, the RMSL trend is modulated by decadal variations. As observed over the global ocean, the main actors of the long-term RMSL trend are associated with anthropogenic global/regional warming. Decadal variability is mainly linked to the strengthening or weakening of the Atlantic Meridional Overturning Circulation (AMOC) (e.g. Chafik et al., 2019). The latest is driven by the North Atlantic Oscillation (NAO) for decadal (20-30y) timescales (e.g. Delworth and Zeng, 2016). Along the European coast, the NAO also influences the along-slope winds dynamic which in return significantly contributes to the local sea level variability observed (Chafik et al., 2019). ""KEY FINDINGS "" Over the [1999/02/21 to 2023/12/31] period, the area-averaged sea level in the IBI area rises at a rate of 4.7 ± 0.8 mm/yr with an acceleration of 0.30 ± 0.06 mm/yr². This trend estimation is based on the altimeter measurements corrected from global GIA correction (Spada et Melini, 2019) to consider the ongoing movement of land. T The TOPEX-A is no longer included in the computation of regional mean sea level parameters (trend and acceleration) with version 2024 products due to potential drifts, and ongoing work aims to develop a new empirical correction. Calculation begins in February 1999 (the start of the TOPEX-B period). '''DOI (product):''' https://doi.org/10.48670/moi-00252

'''DEFINITION''' The Strong Wave Incidence index is proposed to quantify the variability of strong wave conditions in the Iberia-Biscay-Ireland regional seas. The anomaly of exceeding a threshold of Significant Wave Height is used to characterize the wave behavior. A sensitivity test of the threshold has been performed evaluating the differences using several ones (percentiles 75, 80, 85, 90, and 95). From this indicator, it has been chosen the 90th percentile as the most representative, coinciding with the state-of-the-art. Two CMEMS products are used to compute the Strong Wave Incidence index: • IBI-WAV-MYP: IBI_REANALYSIS_WAV_005_006 • IBI-WAV-NRT: IBI_ANALYSIS_FORECAST_WAV_005_005 The Strong Wave Incidence index (SWI) is defined as the difference between the climatic frequency of exceedance (Fclim) and the observational frequency of exceedance (Fobs) of the threshold defined by the 90th percentile (ThP90) of Significant Wave Height (SWH) computed on a monthly basis from hourly data of IBI-WAV-MYP product: SWI = Fobs(SWH > ThP90) – Fclim(SWH > ThP90) Since the Strong Wave Incidence index is defined as a difference of a climatic mean and an observed value, it can be considered an anomaly. Such index represents the percentage that the stormy conditions have occurred above/below the climatic average. Thus, positive/negative values indicate the percentage of hourly data that exceed the threshold above/below the climatic average, respectively. '''CONTEXT''' Ocean waves have a high relevance over the coastal ecosystems and human activities. Extreme wave events can entail severe impacts over human infrastructures and coastal dynamics. However, the incidence of severe (90th percentile) wave events also have valuable relevance affecting the development of human activities and coastal environments. The Strong Wave Incidence index based on the CMEMS regional analysis and reanalysis product provides information on the frequency of severe wave events. The IBI-MFC covers the Europe’s Atlantic coast in a region bounded by the 26ºN and 56ºN parallels, and the 19ºW and 5ºE meridians. The western European coast is located at the end of the long fetch of the subpolar North Atlantic (Mørk et al., 2010), one of the world’s greatest wave generating regions (Folley, 2017). Several studies have analyzed changes of the ocean wave variability in the North Atlantic Ocean (Bacon and Carter, 1991; Kursnir et al., 1997; WASA Group, 1998; Bauer, 2001; Wang and Swail, 2004; Dupuis et al., 2006; Wolf and Woolf, 2006; Dodet et al., 2010; Young et al., 2011; Young and Ribal, 2019). The observed variability is composed of fluctuations ranging from the weather scale to the seasonal scale, together with long-term fluctuations on interannual to decadal scales associated with large-scale climate oscillations. Since the ocean surface state is mainly driven by wind stresses, part of this variability in Iberia-Biscay-Ireland region is connected to the North Atlantic Oscillation (NAO) index (Bacon and Carter, 1991; Hurrell, 1995; Bouws et al., 1996, Bauer, 2001; Woolf et al., 2002; Tsimplis et al., 2005; Gleeson et al., 2017). However, later studies have quantified the relationships between the wave climate and other atmospheric climate modes such as the East Atlantic pattern, the Arctic Oscillation pattern, the East Atlantic Western Russian pattern and the Scandinavian pattern (Izaguirre et al., 2011, Matínez-Asensio et al., 2016). The Strong Wave Incidence index provides information on incidence of stormy events in four monitoring regions in the IBI domain. The selected monitoring regions (Figure 1.A) are aimed to provide a summarized view of the diverse climatic conditions in the IBI regional domain: Wav1 region monitors the influence of stormy conditions in the West coast of Iberian Peninsula, Wav2 region is devoted to monitor the variability of stormy conditions in the Bay of Biscay, Wav3 region is focused in the northern half of IBI domain, this region is strongly affected by the storms transported by the subpolar front, and Wav4 is focused in the influence of marine storms in the North-East African Coast, the Gulf of Cadiz and Canary Islands. More details and a full scientific evaluation can be found in the CMEMS Ocean State report (Pascual et al., 2020). '''CMEMS KEY FINDINGS''' The analysis of the index in the last decades do not show significant trends of the strong wave conditions over the period 1992-2021 with 99% confidence. The maximum wave event reported in region WAV1 (B) occurred in February 2014, producing an increment of 25% of strong wave conditions in the region. Two maximum wave events are found in WAV2 (C) with an increment of 15% of high wave conditions in November 2009 and February 2014. As in regions WAV1 and WAV2, in the region WAV3 (D), a strong wave event took place in February 2014, this event is one of the maximum events reported in the region with an increment of strong wave conditions of 20%, two months before (December 2013) there was a storm of similar characteristics affecting this region, other events of similar magnitude are detected on October 2000 and November 2009. The region WAV4 (E) present its maximum wave event in January 1996, such event produced a 25% of increment of strong wave conditions in the region. Despite of each monitoring region is affected by independent wave events; the analysis shows several past higher-than-average wave events that were propagated though several monitoring regions: November-December 2010 (WAV3 and WAV2); February 2014 (WAV1, WAV2, and WAV3); and February-March 2018 (WAV1 and WAV4). The analysis of the NRT period (January 2022 onwards) depicts a significant event that occurred in November 2022, which affected the WAV2 and WAV3 regions, resulting in a 15% and 25% increase in maximum wave conditions, respectively. In the case of the WAV3 region, this event was the strongest event recorded in this region. In the WAV4 region, an event that occurred in February 2024 was the second most intense on record, showing an 18% increase in strong wave conditions in the region. In the WAV1 region, the NRT period includes two high-intensity events that occurred in February 2024 (21% increase in strong wave conditions) and April 2022 (18% increase in maximum wave conditions). '''Figure caption''' (A) Mean 90th percentile of Sea Wave Height computed from IBI_REANALYSIS_WAV_005_006 product at an hourly basis. Gray dotted lines denote the four monitoring areas where the Strong Wave Incidence index is computed. (B, C, D, and E) Strong Wave Incidence index averaged in monitoring regions WAV1 (A), WAV2 (B), WAV3 (C), and WAV4 (D). Panels show merged results of two CMEMS products: IBI_REANALYSIS_WAV_005_006 (blue), IBI_ANALYSIS_FORECAST_WAV_005_005 (orange). The trend and 99% confidence interval of IBI-MYP product is included (bottom right). '''DOI (product):''' https://doi.org/10.48670/moi-00251

'''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 MY 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 2021 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''' The coastal and shelf waters, especially between 30 and 400N that showed active oligotrophication flags for 2020 have reduced in 2021 and a reversal to eutrophic flags can be seen in places. Again, the eutrophication index is positive only for a small number of coastal locations just north of 40oN in 2021, however south of 40oN there has been a significant increase in eutrophic flags, particularly around the Azores. In general, the 2021 indicator map showed an increase in oligotrophic areas in the Northern Atlantic and an increase in eutrophic areas in the Southern Atlantic. 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