Moving 6-year analysis of Water_body_silicate in the Black Sea for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: from 1970-1975 until 2013-2018 Depth range (IODE standard depths): -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Log transformation was used in the analyses. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l. The entire set of related maps can be found in the viewing service: http://ec.oceanbrowser.net/emodnet/

This product displays the stations present in EMODnet validated dataset where DDT 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

Water body dissolved oxygen concentration - Seasonal Climatology for Po River for the period 1960-2019 on the domain: Lon E12.0-14.0 deg.E, Lat N44.0-46.00 deg.N. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30sec topography for the spatial connectivity of water masses. A relative correlation length used based on topography gradients and the horizontal correlation length varies from 25 km to 12.5 km. The vertical correlation length (in meters) is: [4.,4.,4.,4.,4.,10.,10.,10.,10.,20.,20.,20.,20.,20.]. Depth range: [0., 2., 4., 6., 8., 10., 15., 20., 25., 30., 40., 50., 60., 70.]. Units: umol/l. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.009 degrees (around 975m and 708m accordingly). Additional analysis parameters: The error variance (epsilon2) was set equal to 1. Data weight was reduced near the coasts. Data weighting was applied with the following threshold values: lon=0.03, lat=0.03, depth=2. Threshold value for log transformation=200. Seasons definition: Winter=Jan-Mar, Spring=Apr-Jun, Summer=Jul-Sep, Autumn=Oct-Dec. DOI: 10.6092/tasv-9066

Water body dissolved inorganic nitrogen (DIN) - - Seasonal Climatology for Danube Delta (NW Black Sea) for the period 1980-2018 on the domain: Lon E28.5-30.5°E, Lat N43.7-45.6°N. (winter: December-February, - spring: March-May, - summer: June-August, - autumn: September-November). Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.5, using GEBCO 30sec topography for the spatial connectivity of water masses. Horizontal correlation length of 47.730km and vertical correlation length of 45m were applied. Depth range: 0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 40.0, 50.0. Units: umol/l. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees.

Water body chlorophyll-a - Seasonal Climatology for Danube Delta (NW Black Sea) for the period 2001-2018 on the domain: Lon E28.5-30.5°E, Lat N43.7-45.6°N. (winter: December-February, - spring: March-May, - summer: June-August, - autumn: September-November) Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.5, using GEBCO 30sec topography for the spatial connectivity of water masses. Horizontal correlation length of 47.730km and vertical correlation length of 45m were applied. Depth range: 0, 2.0, 4.0, 6.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0 m. Units: mg/m^3. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees.

This product displays the stations present in EMODnet validated dataset where fluoranthene levels have been measured in water. 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

Moving 6-year analysis of Water_body_dissolved_inorganic_nitrogen_(DIN) in the Black Sea.for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: - from 1990-1995 until 2013-2018 (spring), - from 1990-1995 until 2013-2018 (summer), - from 1990-1995 until 2013-2018 (autumn), - from 1990-1995 until 2011-2016 (winter). Depth range (IODE standard depths): -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l. The entire set of related maps can be found in the viewing service: http://ec.oceanbrowser.net/emodnet/ .

Water body silicate - Seasonal Climatology for Gulf of Riga for the period 1994-2019 on the domain: Lon 22.3-25.0°E, Lat 56.8-58.4°N. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: The computation was done with the DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7, using GEBCO 30sec topography for the spatial connectivity of water masses. Horizontal correlation length of 25 km and vertical correlation length of 10 m. Depth range: 0.0, 2.5, 5, 10, 15, 20, 25, 30, 40, 50 m. Units: µmol/l. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees.

Moving 6-year analysis of Dissolved Inorganic Nitrogen at Atlantic 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 centred average of each season. 6-year periods span - from 1988-1993 until 2015-2020 (winter), - from 1988-1993 until 2015-2020 (spring), - from 1986-1991 until 2015-2020 (summer), - from 1984-1989 until 2015-2020 (autumn). Observational data span from 1984 to 2020. Depth range (IODE standard depths): -2000.0, -1750, -1500.0, -1400.0, -1300.0, -1200.0, -1100.0, -1000.0, -900.0, -800.0, -700.0, -600.0, -500.0, -400.0, -300.0, -250.0, -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analyzed filed masked using relative error threshold 0.3 and 0.5 DIVA settings. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Signal to noise ratio was fixed to 1. Logarithmic transformation applied to the data prior to the analysis. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l

Moving 6-year analysis of Water_body_phosphate in the Black Sea for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: from 1970-1975 until 2013-2018 Depth range (IODE standard depths): -200.0, -150.0, -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -20.0, -10.0, -5.0, -0.0. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVA (Data-Interpolating Variational Analysis) tool. GEBCO 1min topography is used for the contouring preparation. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Log transformation was used in the analyses. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l. The entire set of related maps can be found in the viewing service: http://ec.oceanbrowser.net/emodnet/