'''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 Global Ocean- In-situ observation yearly delivery in delayed mode. The In Situ delayed mode product designed for reanalysis purposes integrates the best available version of in situ data for temperature and salinity measurements. These data are collected from main global networks (Argo, GOSUD, OceanSITES, World Ocean Database) completed by European data provided by EUROGOOS regional systems and national system by the regional INS TAC components. It is updated on a yearly basis. This version is a merged product between the previous verion of CORA and EN4 distributed by the Met Office for the period 1950-1990. '''DOI (product) :''' https://doi.org/10.17882/46219

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

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

'''Short description:''' DTU Space produces polar covering Near Real Time gridded ice displacement fields obtained by MCC processing of Sentinel-1 SAR, Envisat ASAR WSM swath data or RADARSAT ScanSAR Wide mode data . The nominal temporal span between processed swaths is 24hours, the nominal product grid resolution is a 10km. '''DOI (product) :''' https://doi.org/10.48670/moi-00135

'''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (MY) product as distributed by CMEMS. 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. '''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 trend average for the 1997-2023 period is slightly negative (-0.73±0.65% per year) emphasising the results obtained from previous release (1997-2022). This result is in contrast with the analysis of Sathyendranath et al. (2018) that reveals an increasing trend in chlorophyll concentration in all the European Seas. Starting from 2010-2011, except for 2018-2019, the decrease of chlorophyll concentrations is quite evident in the deseasonalized timeseries (in green), and in the maxima of the observations (grey line), starting from 2015. This attenuation of chlorophyll values of the last decade, results in an overall negative trend for 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 '''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. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00195

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The time series are derived from the regional chlorophyll reprocessed (REP) product as distributed by CMEMS. This dataset, derived from multi-sensor (SeaStar-SeaWiFS, AQUA-MODIS, NOAA20-VIIRS, NPP-VIIRS, Envisat-MERIS and Sentinel3A-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. (2021). 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. '''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). '''CMEMS KEY FINDINGS''' In the Mediterranean Sea, the trend average for the 1997-2020 period is slightly negative (-0.580.62% per year). Due to the change in processing techniques and chlorophyll retrieval, this trend estimate cannot be compared directly to those previously reported. The observations time series (in grey) shows minima values have been quite constant until 2015 and then there is a little decrease up to 2020, when an absolute minimum occurs with values lower than 0.04 mg m-3. Throughout the time series, maxima are variable year by year (with absolute maximum in 2015, >0.14 mg m-3), showing an evident reduction since 2016. In the last years of the series, the decrease of chlorophyll concentrations is also observed in the deseasonalized timeseries (in green) with a marked step in 2020. This attenuation of chlorophyll values in the last years results in an overall negative trend for the Mediterranean Sea. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''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:''' For the North Atlantic and Arctic oceans, the ESA Ocean Colour CCI Remote Sensing Reflectance (merged, bias-corrected Rrs) data are used to compute surface Chlorophyll (mg m-3, 1 km resolution) using the regional OC5CCI chlorophyll algorithm. The Rrs are generated by merging the data from SeaWiFS, MODIS-Aqua, MERIS, VIIRS and OLCI-3A sensors and realigning the spectra to that of the MERIS sensor. The algorithm used is OC5CCI - a variation of OC5 (Gohin et al., 2002) developed by IFREMER in collaboration with PML. As part of this development, an OC5CCI look up table was generated specifically for application over OC-CCI merged daily remote sensing reflectances. The resulting OC5CCI algorithm was tested and selected through an extensive calibration exercise that analysed the quantitative performance against in situ data for several algorithms in these specific regions. Processing information: PML's Remote Sensing Group has the capability to automatically receive, archive, process and map global data from multiple polar-orbiting sensors in both near-real time and delayed time. OLCI products are downloaded at level-2 from CODA, the Copernicus Hub and/or via EUMETCAST. These products are remapped at nominal 300m and 1 Km spatial resolution using cylindrical equirectangular projection. 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. By comparing reflectances at different wavelengths and calibrating the result against in situ measurements, an estimate of chlorophyll content can be derived. '''Processing information:''' ESA OC-CCI Rrs raw data are provided by Plymouth Marine Laboratory, currently at 4km resolution globally. These are processed to produce chlorophyll concentration using the same in-house software as in the operational processing. The entire CCI data set is consistent and processing is done in one go. Both OC CCI and the REP product are versioned. Standard masking criteria for detecting clouds or other contamination factors have been applied during the generation of the Rrs, i.e., land, cloud, sun glint, atmospheric correction failure, high total radiance, large solar zenith angle (70deg), large spacecraft zenith angle (56deg), coccolithophores, negative water leaving radiance, and normalized water leaving radiance at 560 nm 0.15 Wm-2 sr-1 (McClain et al., 1995). For the regional products, a variant of the OC-CCI chain is run to produce high resolution data at the 1km resolution necessary. A detailed description of the ESA OC-CCI processing system can be found in OC-CCI (2014e). '''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. By comparing reflectances at different wavelengths and calibrating the result against in-situ measurements, an estimate of chlorophyll content can be derived. '''Quality / Accuracy / Calibration information:''' Detailed description of cal/val is given in the relevant QUID, associated validation reports and quality documentation. '''DOI (product) :''' https://doi.org/10.48670/moi-00070

''' 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). [https://help.marine.copernicus.eu/en/articles/4444611-how-to-cite-or-reference-copernicus-marine-products-and-services How to cite] '''DOI (product) :''' https://doi.org/10.48670/moi-00170