'''This product has been archived'''                For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' The ocean monitoring indicator of regional mean sea level 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 the Copernicus Marine Service (SEALEVEL_GLO_PHY_CLIMATE_L4_MY_008_057). The mean sea level evolution estimated in the Mediterranean 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 removed (least square fit of sinusoidal function) and the time series is low-pass filtered (175 days cut-off). The curve is corrected for the regional mean effect of the Glacial Isostatic Adjustment (GIA) using the ICE5G-VM2 GIA model (Peltier, 2004). During 1993-1998, the Global men sea level (hereafter GMSL) has been known to be affected by a TOPEX-A instrumental drift (WCRP Global Sea Level Budget Group, 2018; Legeais et al., 2020). This drift led to overestimate the trend of the GMSL during the first 6 years of the altimetry record (about 0.04 mm/y at global scale over the whole altimeter period). A correction of the drift is proposed for the Global mean sea level (Legeais et al., 2020). Whereas this TOPEX-A instrumental drift should also affect the regional mean sea level (hereafter RMSL) trend estimation, this empirical correction is currently not applied to the altimeter sea level dataset and resulting estimated for RMSL. Indeed, the pertinence of the global correction applied at regional scale has not been demonstrated yet and there is no clear consensus achieved on the way to proceed at regional scale. Additionally, the estimate of such a correction at regional scale is not obvious, especially in areas where few accurate independent measurements (e.g. in situ)- necessary for this estimation - are available. The trend uncertainty is provided in a 90% confidence interval (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 taken into account. '''CONTEXT''' The indicator on area averaged sea level is a crucial index of climate change, and individual components contribute to sea level rise, including expansion due to ocean warming and melting of glaciers and ice sheets (WCRP Global Sea Level Budget Group, 2018). According to the recent IPCC 6th assessment report, global mean sea level (GMSL) increased by 0.20 (0.15 to 0.25) m over the period 1901 to 2018 with a rate 25 of rise that has accelerated since the 1960s to 3.7 (3.2 to 4.2) mm yr-1 for the period 2006–2018. Human activity was very likely the main driver of observed GMSL rise since 1970 (IPCC WGII, 2021). The weight of the different contributions evolves with time and in the recent decades the mass change has increased, contributing to the on-going acceleration of the GMSL trend (IPCC, 2022a; Legeais et al., 2020; Horwath et al., 2022). At regional scale, sea level does not change homogenously, and RMSL rise can also be 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). Rising sea level can strongly affect population and infrastructures in coastal areas, increase their vulnerability and risks for food security, particularly in low lying areas and island states. Adverse impacts from floods, storms and tropical cyclones with related losses and damages have increased due to sea level rise, and increase their vulnerability and increase risks for food security, particularly in low lying areas and island states (IPCC, 2022b). Adaptation and mitigation measures such as the restoration of mangroves and coastal wetlands, reduce the risks from sea level rise (IPCC, 2022c). Beside a clear long-term trend, the regional mean sea level variation in the Mediterranean Sea shows an important interannual variability, with a high trend observed before 1999 and lower values afterward. This variability is associated with a variation of the different forcing. Steric effect has been the most important forcing before 1999 (Fenoglio-Marc, 2002; Vigo et al., 2005). Important change of the deep-water formation site also occurred in 1995. The latest is preconditioned by an important change of the sea surface circulation observed in the Ionian Sea in 1997-1998 (e.g. Gačić et al., 2011), under the influence of the North Atlantic Oscillation (NAO) and negative Atlantic Multidecadal Oscillation (AMO) phases (Incarbona et al., 2016). They may also impact the sea level trend in the basin (Vigo et al., 2005). In 2010-2011, high regional mean sea level has been related to enhanced water mass exchange at Gibraltar, under the influence of wind forcing during the negative phase of NAO (Landerer and Volkov, 2013). '''CMEMS KEY FINDINGS''' Over the [1993/01/01, 2021/08/02] period, the basin-wide RMSL in the Mediterranean Sea rises at a rate of 2.7  0.83 mm/year. '''DOI (product):''' https://doi.org/10.48670/moi-00264

Output of the 2019 EUSeaMap broad-scale predictive model, produced by EMODnet Seabed Habitats and aggregated into the Benthic Broad Habitat Types of the Marine Strategy Framework Directive (as defined in the Commission Decision 17 May 2017). The extent of the mapped area includes the Mediterranean Sea, Black Sea, Baltic Sea, and areas of the North Eastern Atlantic extending from the Canary Islands in the south to the Barents Sea in the north. The map was produced using a "top-down" modelling approach using classified habitat descriptors to determine a final output habitat. Habitat descriptors differ per region but include: - Biological zone - Energy class - Oxygen regime - Salinity regime - Seabed substrate - Riverine input Habitat descriptors (excepting Substrate) are calculated using underlying physical data and thresholds derived from statistical analyses or expert judgement on known conditions.

'''DEFINITION''' Based on daily, global climate sea surface temperature (SST) analyses generated by the Copernicus Climate Change Service (C3S) (product SST-GLO-SST-L4-REP-OBSERVATIONS-010-024). Analysis of the data was based on the approach described in Mulet et al. (2018) and is described and discussed in Good et al. (2020). The processing steps applied were: 1. The daily analyses were averaged to create monthly means. 2. A climatology was calculated by averaging the monthly means over the period 1991 - 2020. 3. Monthly anomalies were calculated by differencing the monthly means and the climatology. 4. The time series for each grid cell was passed through the X11 seasonal adjustment procedure, which decomposes a time series into a residual seasonal component, a trend component and errors (e.g., Pezzulli et al., 2005). The trend component is a filtered version of the monthly time series. 5. The slope of the trend component was calculated using a robust method (Sen 1968). The method also calculates the 95% confidence range in the slope. '''CONTEXT''' Sea surface temperature (SST) is one of the Essential Climate Variables (ECVs) defined by the Global Climate Observing System (GCOS) as being needed for monitoring and characterising the state of the global climate system (GCOS 2010). It provides insight into the flow of heat into and out of the ocean, into modes of variability in the ocean and atmosphere, can be used to identify features in the ocean such as fronts and upwelling, and knowledge of SST is also required for applications such as ocean and weather prediction (Roquet et al., 2016). '''CMEMS KEY FINDINGS''' Warming trends occurred over most of the globe between 1982 and 2024, with the strongest warming in the Northern Pacific and Atlantic Oceans. However, there were cooling trends in parts of the Southern Ocean and the South-East Pacific Ocean. '''DOI (product):''' https://doi.org/10.48670/moi-00243

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''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). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00260

This product displays the stations present in EMODnet validated dataset where tributyltin levels have been measured in sediment. EMODnet Chemistry has included the gathering of contaminants data since the beginning of the project in 2009. For the maps for EMODnet Chemistry Phase III, it was requested to plot data per matrix (water,sediment, biota), per biological entity and per chemical substance. The series of relevant map products have been developed according to the criteria D8C1 of the MSFD Directive, specifically focusing on the requirements under the new Commission Decision 2017/848 (17th May 2017). The Commission Decision points to relevant threshold values that are specified in the WFD, as well as relating how these contaminants should be expressed (units and matrix etc.) through the related Directives i.e. Priority substances for Water. EU EQS Directive does not fix any threshold values in sediments. On the contrary Regional Sea Conventions provide some of them, and these values have been taken into account for the development of the visualization products. To produce the maps the following process has been followed: 1. Data collection through SeaDataNet standards (CDI+ODV) 2. Harvesting, harmonization, validation and P01 code decomposition of data 3. SQL query on data sets from point 2 4. Production of map with each point representing at least one record that match the criteria The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols. Preliminary processing were necessary to harmonize all the data : • For water: contaminants in the dissolved phase; • For sediment: data on total sediment (regardless of size class) or size class < 2000 μm • For biota: contaminant data will focus on molluscs, on fish (only in the muscle), and on crustaceans • Exclusion of data values equal to 0

This product displays the stations present in EMODnet validated dataset where cadmium levels have been measured in sediment. EMODnet Chemistry has included the gathering of contaminants data since the beginning of the project in 2009. For the maps for EMODnet Chemistry Phase III, it was requested to plot data per matrix (water,sediment, biota), per biological entity and per chemical substance. The series of relevant map products have been developed according to the criteria D8C1 of the MSFD Directive, specifically focusing on the requirements under the new Commission Decision 2017/848 (17th May 2017). The Commission Decision points to relevant threshold values that are specified in the WFD, as well as relating how these contaminants should be expressed (units and matrix etc.) through the related Directives i.e. Priority substances for Water. EU EQS Directive does not fix any threshold values in sediments. On the contrary Regional Sea Conventions provide some of them, and these values have been taken into account for the development of the visualization products. To produce the maps the following process has been followed: 1. Data collection through SeaDataNet standards (CDI+ODV) 2. Harvesting, harmonization, validation and P01 code decomposition of data 3. SQL query on data sets from point 2 4. Production of map with each point representing at least one record that match the criteria The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols. Preliminary processing were necessary to harmonize all the data : • For water: contaminants in the dissolved phase; • For sediment: data on total sediment (regardless of size class) or size class < 2000 μm • For biota: contaminant data will focus on molluscs, on fish (only in the muscle), and on crustaceans • Exclusion of data values equal to 0

This product displays the stations where fluoranthene has been measured in biota and the values present in EMODnet Chemistry infrastructure are not compliant with EQSD, i.e for the substance, in that station, quality values found in EMODnet validated dataset can be equal to 6 or Q and values are above the 30% of the value established by EQSD for fluoranthene in biota. It is necessary to take into account that LOD/LOQ can change with time. These products aggregate data by station, producing only one final value for each station. EMODnet Chemistry has included the gathering of contaminants data since the beginning of the project in 2009. For the maps for EMODnet Chemistry Phase III, it was requested to plot data per matrix (water,sediment, biota), per biological entity and per chemical substance. The series of relevant map products have been developed according to the criteria D8C1 of the MSFD Directive, specifically focusing on the requirements under the new Commission Decision 2017/848 (17th May 2017). The Commission Decision points to relevant threshold values that are specified in the WFD, as well as relating how these contaminants should be expressed (units and matrix etc.) through the related Directives i.e. Priority substances for Water. EU EQS Directive does not fix any threshold values in sediments. On the contrary Regional Sea Conventions provide some of them, and these values have been taken into account for the development of the visualization products. To produce the maps the following process has been followed: 1. Data collection through SeaDataNet standards (CDI+ODV) 2. Harvesting, harmonization, validation and P01 code decomposition of data 3. SQL query on data sets from point 2 4. Production of map with each point representing at least one record that match the criteria The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols. Preliminary processing were necessary to harmonize all the data : • For water: contaminants in the dissolved phase; • For sediment: data on total sediment (regardless of size class) or size class < 2000 μm • For biota: contaminant data will focus on molluscs, on fish (only in the muscle), and on crustaceans • Exclusion of data values equal to 0

This product displays the stations present in EMODnet validated dataset where benzo[A]pyrene levels have been measured in sediment. EMODnet Chemistry has included the gathering of contaminants data since the beginning of the project in 2009. For the maps for EMODnet Chemistry Phase III, it was requested to plot data per matrix (water,sediment, biota), per biological entity and per chemical substance. The series of relevant map products have been developed according to the criteria D8C1 of the MSFD Directive, specifically focusing on the requirements under the new Commission Decision 2017/848 (17th May 2017). The Commission Decision points to relevant threshold values that are specified in the WFD, as well as relating how these contaminants should be expressed (units and matrix etc.) through the related Directives i.e. Priority substances for Water. EU EQS Directive does not fix any threshold values in sediments. On the contrary Regional Sea Conventions provide some of them, and these values have been taken into account for the development of the visualization products. To produce the maps the following process has been followed: 1. Data collection through SeaDataNet standards (CDI+ODV) 2. Harvesting, harmonization, validation and P01 code decomposition of data 3. SQL query on data sets from point 2 4. Production of map with each point representing at least one record that match the criteria The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols. Preliminary processing were necessary to harmonize all the data : • For water: contaminants in the dissolved phase; • For sediment: data on total sediment (regardless of size class) or size class < 2000 μm • For biota: contaminant data will focus on molluscs, on fish (only in the muscle), and on crustaceans • Exclusion of data values equal to 0

This raster dataset represents the physical loss of the seabed in the European seas. Physical habitat loss is an extreme pressure on the marine ecosystem which was defined to include all impacts on the seabed which take more than 12 years to recover, a time span influencing even long living marine mammals and seabirds. Habitat is lost if its substrate, morphology or topography is permanently altered. The dataset is an aggregation of several different human activities that cause physical loss in the seabed: dredging and dumping of dredged material, oil and gas rigs, offshore installations, ports, sand and gravel extraction and windfarms. The resulting dataset is a raster (10km grid cell) derived from EMODnet, MED-IAMER and 4C Offshore datasets and with temporal reference from 2012 to 2017. Higher values mean more loss causing activity in the area. Due to its resolution, loss causing activity in a cell does not usually mean that all the seabed habitat in the area is lost. This dataset has been prepared for the calculation of the combined effect index, produced for the ETC/ICM Report 4/2019 "Multiple pressures and their combined effects in Europe's seas" available on: https://www.eionet.europa.eu/etcs/etc-icm/etc-icm-report-4-2019-multiple-pressures-and-their-combined-effects-in-europes-seas-1.

'''This product has been archived''' For operationnal and online products, please visit https://marine.copernicus.eu '''DEFINITION''' Oligotrophic subtropical gyres are regions of the ocean with low levels of nutrients required for phytoplankton growth and low levels of surface chlorophyll-a whose concentration can be quantified through satellite observations. The gyre boundary has been defined using a threshold value of 0.15 mg m-3 chlorophyll for the Atlantic gyres (Aiken et al. 2016), and 0.07 mg m-3 for the Pacific gyres (Polovina et al. 2008). The area inside the gyres for each month is computed using monthly chlorophyll data from which the monthly climatology is subtracted to compute anomalies. A gap filling algorithm has been utilized to account for missing data. Trends in the area anomaly are then calculated for the entire study period (September 1997 to December 2020). '''CONTEXT''' Oligotrophic gyres of the oceans have been referred to as ocean deserts (Polovina et al. 2008). They are vast, covering approximately 50% of the Earth’s surface (Aiken et al. 2016). Despite low productivity, these regions contribute significantly to global productivity due to their immense size (McClain et al. 2004). Even modest changes in their size can have large impacts on a variety of global biogeochemical cycles and on trends in chlorophyll (Signorini et al. 2015). Based on satellite data, Polovina et al. (2008) showed that the areas of subtropical gyres were expanding. The Ocean State Report (Sathyendranath et al. 2018) showed that the trends had reversed in the Pacific for the time segment from January 2007 to December 2016. '''CMEMS KEY FINDINGS''' The trend in the North Atlantic gyre area for the 1997 Sept – 2020 December period was positive, with a 0.39% year-1 increase in area relative to 2000-01-01 values. This trend has decreased compared with the 1997-2019 trend of 0.45%, and is statistically significant (p<0.05). During the 1997 Sept – 2020 December period, the trend in chlorophyll concentration was positive (0.24% year-1) inside the North Atlantic gyre relative to 2000-01-01 values. This time series extension has resulted in a reversal in the rate of change, compared with the -0.18% trend for the 1997-209 period and is statistically significant (p<0.05). Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00226