EMODnet Chemistry aims to provide access to marine chemistry data sets and derived data products concerning eutrophication, ocean acidification, contaminants and litter. The chosen parameters are relevant for the Marine Strategy Framework Directive (MSFD), in particular for descriptors 5, 8, 9 and 10. The dataset contains standardized, harmonized and validated data collections from beach litter (monitoring and other sources). Datasets concerning beach and seafloor litter data are loaded in a central database after a semi-automated validation phase. Once loaded, a data assessment is performed in order to check data consistency and potential errors are corrected thanks to a feedback loop with data originators. For beach litter, the harmonized datasets contain all unrestricted EMODnet Chemistry data on beach litter, including monitoring data, data from cleaning surveys and data from research. A relevant part of the monitoring data has been considered for assessment purposes by the European institutions and therefore is tagged as MSFD_monitoring. EMODnet beach litter data and databases are hosted and maintained by 'Istituto Nazionale di Oceanografia e di Geofisica Sperimentale, Division of Oceanography (OGS/NODC)' from Italy. Data are formatted following Guidelines and forms for gathering marine litter data, which can be found at: https://doi.org/10.6092/15c0d34c-a01a-4091-91ac-7c4f561ab508 The updated vocabularies of admitted values are available at: https://nodc.ogs.it/marinelitter/vocab The harmonized datasets can be downloaded as EMODnet Beach litter data format Version 7.0, which is a spreadsheet file composed of 4 sheets: beach metadata, survey metadata, animals and litter. Local_CDI field in the survey metadata sheet allows to retrieve the original data.

This visualization product displays the fishing & aquaculture related plastic items abundance of marine macro-litter (> 2.5cm) per beach per year from Marine Strategy Framework Directive (MSFD) monitoring surveys. 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 plastic bags 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); - 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 plastic bags related items of the survey (normalized by 100 m) = Number of plastic bags related items of the survey x (100 / survey length) Then, this normalized number of plastic bags related items is summed to obtain the total normalized number of plastic bags related items for each survey. Finally, the median abundance of plastic bags related items for each beach and year is calculated from these normalized abundances of plastic bags 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 plastic bags related items from MSFD data 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.

'''DEFINITION''' The temporal evolution of thermosteric sea level in an ocean layer is obtained from an integration of temperature driven ocean density variations, which are subtracted from a reference climatology to obtain the fluctuations from an average field. The products used include three global reanalyses: GLORYS, C-GLORS, ORAS5 (GLOBAL_MULTIYEAR_PHY_ENS_001_031) and two in situ based reprocessed products: CORA5.2 (INSITU_GLO_PHY_TS_OA_MY_013_052) , ARMOR-3D (MULTIOBS_GLO_PHY_TSUV_3D_MYNRT_015_012). Additionally, the time series based on the method of von Schuckmann and Le Traon (2011) has been added. The regional thermosteric sea level values are then averaged from 60°S-60°N aiming to monitor interannual to long term global sea level variations caused by temperature driven ocean volume changes through thermal expansion as expressed in meters (m). '''CONTEXT''' The global mean sea level is reflecting changes in the Earth’s climate system in response to natural and anthropogenic forcing factors such as ocean warming, land ice mass loss and changes in water storage in continental river basins. Thermosteric sea-level variations result from temperature related density changes in sea water associated with volume expansion and contraction (Storto et al., 2018). Global thermosteric sea level rise caused by ocean warming is known as one of the major drivers of contemporary global mean sea level rise (Cazenave et al., 2018; Oppenheimer et al., 2019). '''CMEMS KEY FINDINGS''' Since the year 2005 the upper (0-2000m) near-global (60°S-60°N) thermosteric sea level rises at a rate of 1.3±0.3 mm/year. Note: The key findings will be updated annually in November, in line with OMI evolutions. '''DOI (product):''' https://doi.org/10.48670/moi-00240

Moving 6-year analysis of Water body phosphate in the NorthEast Atlantic 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 centred average of each season. 6-year periods span from 1950/1955 until 2019/2024. Observation data span from 1950 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.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 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 250km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. 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: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

'''Short description:''' The NWSHELF_ANALYSISFORECAST_PHY_LR_004_001 is produced by a coupled hydrodynamic-biogeochemical model system with tides, implemented over the North East Atlantic and Shelf Seas at 7 km of horizontal resolution and 24 vertical levels. The product is updated daily, providing 7-day forecast for temperature, salinity, currents, sea level and mixed layer depth. Products are provided at quarter-hourly, hourly, daily de-tided (with Doodson filter), and monthly frequency. '''DOI (product) :''' https://doi.org/10.48670/mds-00367

Moving 6-year analysis of Water body dissolved inorganic nitrogen (DIN) 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 1990-1995 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 21 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.]. 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.]. 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 (1990-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.

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

Moving 6-year analysis of Water body dissolved inorganic nitrogen in the NorthEast Atlantic 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 centred average of each season. 6-year periods span from 1971/1976 until 2019/2024. Observation data span from 1971 to 2024. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.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 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. 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: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

This visualization product displays the density of floating micro-litter per net normalized in grams per km² 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 (weight of particles per km²) = Micro-litter weight / (Sampling effort (km) * Net opening (cm) * 0.00001) When the weight of microlitters 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.

'''DEFINITION''' Important note to users: These data are not to be used for navigation. The data is 100 m resolution and as high quality as possible. It has been produced with state-of-the-art technology and validated to the best of the producer’s ability and where sufficient high-quality data were available. These data could be useful for planning and modelling purposes. The user should independently assess the adequacy of any material, data and/or information of the product before relying upon it. Neither Mercator Ocean International/Copernicus Marine Service nor the data originators are liable for any negative consequences following direct or indirect use of the product information, services, data products and/or data. Product overview: This is a satellite derived bathymetry product covering the global coastal area (where data retrieval is possible), with 100 m resolution, based on Sentinel-2. This global coastal product has been developed based on 3 methodologies: Intertidal Satellite-Derived Bathymetry; Physics-based optical Satellite-Derived Bathymetry from RTE inversion; and Wave Kinematics Satellite-Derived Bathymetry from wave dispersion. There is one dataset for each of the methods (including a quality index based on uncertainty) and an additional one where the three datasets were merged (also includes a quality index). Using their expertise and special techniques the consortium tried to achieve an optimal balance between coverage and data quality. '''DOI (product):''' https://doi.org/10.48670/mds-00364