This visualization product displays the cigarette related items abundance of marine macro-litter (> 2.5cm) per beach per year from Marine Strategy Framework Directive (MSFD) monitoring surveys without UNEP-MARLIN data. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of beach litter have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols and reference lists used on a European scale. Preliminary processings were necessary to harmonize all the data: - Exclusion of OSPAR 1000 protocol: in order to follow the approach of OSPAR that it is not including these data anymore in the monitoring; - Selection of MSFD surveys only (exclusion of other monitoring, cleaning and research operations); - Exclusion of beaches without coordinates; - Selection of cigarette related items only. The list of selected items is attached to this metadata. This list was created using EU Marine Beach Litter Baselines, the European Threshold Value for Macro Litter on Coastlines and the Joint list of litter categories for marine macro-litter monitoring from JRC (these three documents are attached to this metadata); - Exclusion of surveys referring to the UNEP-MARLIN list: the UNEP-MARLIN protocol differs from the other types of monitoring in that cigarette butts are surveyed in a 10m square. To avoid comparing abundances from very different protocols, the choice has been made to distinguish in two maps the cigarette related items results associated with the UNEP-MARLIN list from the others; - Normalization of survey lengths to 100m & 1 survey / year: in some case, the survey length was not exactly 100m, so in order to be able to compare the abundance of litter from different beaches a normalization is applied using this formula: Number of cigarette related items of the survey (normalized by 100 m) = Number of cigarette related items of the survey x (100 / survey length) Then, this normalized number of cigarette related items is summed to obtain the total normalized number of cigarette related items for each survey. Finally, the median abundance of cigarette related items for each beach and year is calculated from these normalized abundances of cigarette related items per survey. Sometimes the survey length was null or equal to 0. Assuming that the MSFD protocol has been applied, the length has been set at 100m in these cases. Percentiles 50, 75, 95 & 99 have been calculated taking into account cigarette related items from MSFD monitoring data (excluding UNEP-MARLIN protocol) for all years. More information is available in the attached documents. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area.

This visualization product displays the density of floating micro-litter per net normalized per km² per year from research and monitoring protocols. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Before 2021, there was no coordinated effort at the regional or European scale for micro-litter. Given this situation, EMODnet Chemistry proposed to adopt the data gathering and data management approach as generally applied for marine data, i.e., populating metadata and data in the CDI Data Discovery and Access service using dedicated SeaDataNet data transport formats. EMODnet Chemistry is currently the official EU collector of micro-litter data from Marine Strategy Framework Directive (MSFD) National Monitoring activities (descriptor 10). A series of specific standard vocabularies or standard terms related to micro-litter have been added to SeaDataNet NVS (NERC Vocabulary Server) Common Vocabularies to describe the micro-litter. European micro-litter data are collected by the National Oceanographic Data Centres (NODCs). Micro-litter map products are generated from NODCs data after a test of the aggregated collection including data and data format checks and data harmonization. A filter is applied to represent only micro-litter sampled according to research and monitoring protocols as MSFD monitoring. Densities were calculated for each net using the following calculation: Density (number of particles per km²) = Micro-litter count / (Sampling effort (km) * Net opening (cm) * 0.00001) When information about the sampling effort (km) was lacking and point coordinates were known (start and end of the sampling), the sampling effort was calculated using the PostGIS ST_DistanceSpheroid function with a WGS84 measurement spheroid. When the number of micro-litters or the net opening was not filled, it was not possible to calculate the density. Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the National Oceanographic Data Centre (NODC) for this area.

This visualization product displays the cigarette related items abundance of marine macro-litter (> 2.5cm) per beach per year from Marine Strategy Framework Directive (MSFD) monitoring surveys without UNEP-MARLIN data. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of beach litter have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols and reference lists used on a European scale. Preliminary processings were necessary to harmonize all the data: - Exclusion of OSPAR 1000 protocol: in order to follow the approach of OSPAR that it is not including these data anymore in the monitoring; - Selection of MSFD surveys only (exclusion of other monitoring, cleaning and research operations); - Exclusion of beaches without coordinates; - Selection of cigarette related items only. The list of selected items is attached to this metadata. This list was created using EU Marine Beach Litter Baselines, the European Threshold Value for Macro Litter on Coastlines and the Joint list of litter categories for marine macro-litter monitoring from JRC (these three documents are attached to this metadata); - Selection of surveys referring to the UNEP-MARLIN list: the UNEP-MARLIN protocol differs from the other types of monitoring in that cigarette butts are surveyed in a 10m square. To avoid comparing abundances from very different protocols, the choice has been made to distinguish in two maps the cigarette related items results associated with the UNEP-MARLIN list from the others; - Normalization of survey lengths to 100m & 1 survey / year: in some case, the survey length was not exactly 100m, so in order to be able to compare the abundance of litter from different beaches a normalization is applied using this formula: Number of cigarette related items of the survey (normalized by 100 m) = Number of cigarette related items of the survey x (100 / survey length) Then, this normalized number of cigarette related items is summed to obtain the total normalized number of cigarette related items for each survey. Finally, the median abundance of cigarette related items for each beach and year is calculated from these normalized abundances of cigarette related items per survey. Sometimes the survey length was null or equal to 0. Assuming that the MSFD protocol has been applied, the length has been set at 100m in these cases. Percentiles 50, 75, 95 & 99 have been calculated taking into account cigarette related items from MSFD monitoring data (excluding UNEP-MARLIN protocol) for all years. More information is available in the attached documents. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area.

EMODnet Chemistry aims to provide access to marine chemistry datasets and derived data products concerning eutrophication, acidity and contaminants. The importance of the selected substances and other parameters relates to the Marine Strategy Framework Directive (MSFD). This aggregated dataset contains all unrestricted EMODnet Chemistry data on potential hazardous substances, despite the fact that some data might not be related to pollution (e.g. collected by deep corer). Temperature, salinity and additional parameters are included when available. It covers the Northeast Atlantic Ocean (40W). Data were harmonised and validated by '‘IFREMER / IDM / SISMER - Scientific Information Systems for the SEA’ in France. The dataset contains water (profiles), sediment (profiles and timeseries) and biota (timeseries). The temporal coverage is 1974–2018 for water measurements, 1966–2022 for sediment measurements and 1979–2023 for biota measurements. Regional datasets concerning contaminants are automatically harvested and the resulting collections are harmonised and validated using ODV Software and following a common methodology for all sea regions ( https://doi.org/10.6092/8b52e8d7-dc92-4305-9337-7634a5cae3f4 ). Parameter names are based on P01 vocabulary, which relates to BODC Parameter Usage Vocabulary and is available at: https://vocab.nerc.ac.uk/search_nvs/P01/ . The harmonised dataset can be downloaded as as an ODV spreadsheet, which is composed of a metadata header followed by tab separated values. This spreadsheet can be imported into ODV Software for visualisation (more information can be found at: https://www.seadatanet.org/Software/ODV ). In addition, the same dataset is offered also as a txt file in a long/vertical format, in which each P01 measurement is a record line. Additionally, there are a series of columns that split P01 terms into subcomponents (substance, CAS number, matrix...).This transposed format is more adapted to worksheet applications (e.g. LibreOffice Calc).

Moving 6-year analysis of Water body silicate in the Mediterranean Sea for each season: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Every year of the time dimension corresponds to the 6-year centered average of the season. 6-years periods span from 1970-1975 until 2018-2023. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 25 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.,1200.,1300.,1400.,1500.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2023) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. Originators of Italian data sets-List of contributors: - Brunetti Fabio (OGS) - Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 - Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 - Cataletto Bruno (OGS) - Cinzia Comici Cinzia (OGS) - Civitarese Giuseppe (OGS) - DeVittor Cinzia (OGS) - Giani Michele (OGS) - Kovacevic Vedrana (OGS) - Mosetti Renzo (OGS) - Solidoro C.,Beran A.,Cataletto B.,Celussi M.,Cibic T.,Comici C.,Del Negro P.,De Vittor C.,Minocci M.,Monti M.,Fabbro C.,Falconi C.,Franzo A.,Libralato S.,Lipizer M.,Negussanti J.S.,Russel H.,Valli G., doi:10.6092/e5518899-b914-43b0-8139-023718aa63f5 - Celio Massimo (ARPA FVG) - Malaguti Antonella (ENEA) - Fonda Umani Serena (UNITS) - Bignami Francesco (ISAC/CNR) - Boldrini Alfredo (ISMAR/CNR) - Marini Mauro (ISMAR/CNR) - Miserocchi Stefano (ISMAR/CNR) - Zaccone Renata (IAMC/CNR) - Lavezza, R., Dubroca, L. F. C., Ludicone, D., Kress, N., Herut, B., Civitarese, G., Cruzado, A., Lefèvre, D.,Souvermezoglou, E., Yilmaz, A., Tugrul, S., and Ribera d'Alcala, M.: Compilation of quality controlled nutrient profiles from the Mediterranean Sea, doi:10.1594/PANGAEA.771907, 2011.

Seasonal climatology of Water body dissolved inorganic nitrogen for Loire river for the period 1971-2024 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1971 to 2024. Depth levels (m): [0.0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0, 130.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 15 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.01 degrees. Horizontal correlation length is defined seasonally (in meters): 72000 (winter), 43000 (spring), 215000 (summer), 150000 (autumn). Vertical correlation length was optimized and vertically filtered and a seasonally-averaged profile was used (DIVAnd.fitvertlen). Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

'''DEFINITION''' Ocean heat content (OHC) is defined here as the deviation from a reference period (1993-2014) and is closely proportional to the average temperature change from z1 = 0 m to z2 = 700 m depth: OHC=∫_(z_1)^(z_2)ρ_0 c_p (T_yr-T_clim )dz [1] with a reference density of = 1030 kgm-3 and a specific heat capacity of cp = 3980 J kg-1 °C-1 (e.g. von Schuckmann et al., 2009). Time series of annual mean values area averaged ocean heat content is provided for the Mediterranean Sea (30°N, 46°N; 6°W, 36°E) and is evaluated for topography deeper than 300m. '''CONTEXT''' Knowing how much and where heat energy is stored and released in the ocean is essential for understanding the contemporary Earth system state, variability and change, as the oceans shape our perspectives for the future. The quality evaluation of MEDSEA_OMI_OHC_area_averaged_anomalies is based on the “multi-product” approach as introduced in the second issue of the Ocean State Report (von Schuckmann et al., 2018), and following the MyOcean’s experience (Masina et al., 2017). Six global products and a regional (Mediterranean Sea) product have been used to build an ensemble mean, and its associated ensemble spread. The reference products are: • The Mediterranean Sea Reanalysis at 1/24 degree horizontal resolution (MEDSEA_MULTIYEAR_PHY_006_004, DOI: https://doi.org/10.25423/CMCC/MEDSEA_MULTIYEAR_PHY_006_004_E3R1, Escudier et al., 2020) • Four global reanalyses at 1/4 degree horizontal resolution (GLOBAL_MULTIYEAR_PHY_ENS_001_031): GLORYS, C-GLORS, ORAS5, FOAM • Two observation based products: CORA (INSITU_GLO_PHY_TS_OA_MY_013_052) and ARMOR3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Details on the products are delivered in the PUM and QUID of this OMI. '''CMEMS KEY FINDINGS''' The ensemble mean ocean heat content anomaly time series over the Mediterranean Sea shows a continuous increase in the period 1993-2022 at rate of 1.38±0.08 W/m2 in the upper 700m. After 2005 the rate has clearly increased with respect the previous decade, in agreement with Iona et al. (2018). '''DOI (product):''' https://doi.org/10.48670/moi-00261

Seasonal climatology of Water body dissolved oxygen concentration for Loire river for the period 1950-2024 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1950 to 2024. Depth levels (m): [0.0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0, 130.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.12, using GEBCO 15 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.01 degrees. Horizontal correlation length is defined seasonally (in meters): 150000 (winter), 100000 (spring), 70000 (summer), 115000 (autumn). Vertical correlation length was optimized and vertically filtered and a seasonally-averaged profile was used (DIVAnd.fitvertlen). Signal-to-noise ratio was fixed to 1 for vertical profiles and 0.1 for time series to account for the redundancy in the time series observations. A logarithmic transformation (DIVAnd.loglin) was applied to the data prior to the analysis. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

'''DEFINITION''' Significant wave height (SWH), expressed in metres, is the average height of the highest third of waves. This OMI provides global maps of the seasonal mean and trend of significant wave height (SWH), as well as time series in three oceanic regions of the same variables and their trends from 2002 to 2020, calculated from the reprocessed global L4 SWH product (WAVE_GLO_PHY_SWH_L4_MY_014_007). The extreme SWH is defined as the 95th percentile of the daily maximum SWH for the selected period and region. The 95th percentile is the value below which 95% of the data points fall, indicating higher than normal wave heights. The mean and 95th percentile of SWH (in m) are calculated for two seasons of the year to take into account the seasonal variability of waves (January, February and March, and July, August and September). Trends have been obtained using linear regression and are expressed in cm/yr. For the time series, the uncertainty around the trend was obtained from the linear regression, while the uncertainty around the mean and 95th percentile was bootstrapped. For the maps, if the p-value obtained from the linear regression is less than 0.05, the trend is considered significant. '''CONTEXT''' Grasping the nature of global ocean surface waves, their variability, and their long-term interannual shifts is essential for climate research and diverse oceanic and coastal applications. The sixth IPCC Assessment Report underscores the significant role waves play in extreme sea level events (Mentaschi et al., 2017), flooding (Storlazzi et al., 2018), and coastal erosion (Barnard et al., 2017). Additionally, waves impact ocean circulation and mediate interactions between air and sea (Donelan et al., 1997) as well as sea-ice interactions (Thomas et al., 2019). Studying these long-term and interannual changes demands precise time series data spanning several decades. Until now, such records have been available only from global model reanalyses or localised in situ observations. While buoy data are valuable, they offer limited local insights and are especially scarce in the southern hemisphere. In contrast, altimeters deliver global, high-quality measurements of significant wave heights (SWH) (Gommenginger et al., 2002). The growing satellite record of SWH now facilitates more extensive global and long-term analyses. By using SWH data from a multi-mission altimetric product from 2002 to 2020, we can calculate global mean SWH and extreme SWH and evaluate their trends, regionally and globally. '''KEY FINDINGS''' From 2002 to 2020, positive trends in both Significant Wave Height (SWH) and extreme SWH are mostly found in the southern hemisphere (a, b). The 95th percentile of wave heights (q95), increases faster than the average values, indicating that extreme waves are growing more rapidly than average wave height (a, b). Extreme SWH’s global maps highlight heavily storms affected regions, including the western North Pacific, the North Atlantic and the eastern tropical Pacific (a). In the North Atlantic, SWH has increased in summertime (July August September) but decreased in winter. Specifically, the 95th percentile SWH trend is decreasing by 2.1 ± 3.3 cm/year, while the mean SWH shows a decrease of 2.2 ± 1.76 cm/year. In the south of Australia, during boreal winter, the 95th percentile SWH is increasing at 2.6 ± 1.5 cm/year (c), with the mean SWH increasing by 0.5 ± 0.66 cm/year (d). Finally, in the Antarctic Circumpolar Current, also in boreal winter, the 95th percentile SWH trend is 3.2 ± 2.14 cm/year (c) and the mean SWH trend is 1.7 ± 0.84 cm/year (d). These patterns highlight the complex and region-specific nature of wave height trends. Further discussion is available in A. Laloue et al. (2024). '''DOI (product):''' https://doi.org/10.48670/mds-00352

This visualization product displays the density of floating micro-litter per net normalized per m³ per year from specific protocols different from research and monitoring protocols. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Before 2021, there was no coordinated effort at the regional or European scale for micro-litter. Given this situation, EMODnet Chemistry proposed to adopt the data gathering and data management approach as generally applied for marine data, i.e., populating metadata and data in the CDI Data Discovery and Access service using dedicated SeaDataNet data transport formats. EMODnet Chemistry is currently the official EU collector of micro-litter data from Marine Strategy Framework Directive (MSFD) National Monitoring activities (descriptor 10). A series of specific standard vocabularies or standard terms related to micro-litter have been added to SeaDataNet NVS (NERC Vocabulary Server) Common Vocabularies to describe the micro-litter. European micro-litter data are collected by the National Oceanographic Data Centres (NODCs). Micro-litter map products are generated from NODCs data after a test of the aggregated collection including data and data format checks and data harmonization. A filter is applied to represent only micro-litter sampled according to a very specific protocol such as the Volvo Ocean Race (VOR) or Oceaneye. Densities were calculated for each net using the following calculation: Density (number of particles per m³) = Micro-litter count / Sampling effort (m³) When the number of micro-litters was not filled, it was not possible to calculate the density. Percentiles 50, 75, 95 & 99 have been calculated taking into account data for all years. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the National Oceanographic Data Centre (NODC) for this area.