Seasonal Climatology of Dissolved Inorganic Nitrogen (DIN) for Loire River for the period 1990-2019 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Observational data span from 1990 to 2019. Depth levels (m): -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -25.0, -20.0, -15.0, -10.0, -8.0, -6.0, -4.0, -2.0, -0.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.4, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. Signal to noise ratio was fixed to 1 for vertical profiles and to 0.1 for time series to account for the redundancy in the time series observations. 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.

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

Seasonal climatology of Water body dissolved inorganic nitrogen for Loire river for the period 1971-2021 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1971 to 2021. 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.4, 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.Anam.loglin) was applied to the data prior to the analysis to avoid unrealistic negative values. Background field: the vertically-filtered data mean profile is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Moving 6-year analysis of Silicate 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 1972-1977 until 2015-2020 (winter), - from 1984-1989 until 2015-2020 (spring), - from 1972-1977 until 2015-2020 (summer), - from 1971-1976 until 2015-2020 (autumn). Observational data span from 1971 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

This product displays the stations where fluoranthene has been measured and the values present in EMODnet Chemistry infrastructure are always below the limit of detection or quantification (LOD/LOQ), i.e quality values found in EMODnet validated dataset can be equal to 6 or Q. 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 (above, below or above/below). 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

Seasonal Climatology of Phosphate for Loire River for the period 1964-2019 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Observational data span from 1964 to 2019. Depth levels (m): -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -25.0, -20.0, -15.0, -10.0, -8.0, -6.0, -4.0, -2.0, -0.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.4, using GEBCO 30sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees. 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. The weight of time series was decreased by a factor of 10 relative to the weight of the profiles to account for the redundancy in the time series observations. Detrending of data: no, Advection constraint applied: no. Units: umol/l.

Water body dissolved oxygen concentration - Monthly Climatology for the European Seas for the period 1960-2020 on the domain from longitude -45.0 to 70.0 degrees East and latitude 24.0 to 83.0 degrees North. 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.9, using GEBCO 30sec topography for the spatial connectivity of water masses. Horizontal correlation length and vertical correlation length vary spatially depending on the topography and domain. Depth range: 0.0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0, 55.0, 60.0, 65.0, 70.0, 75.0, 80.0, 85.0, 90.0, 95.0, 100.0, 125.0, 150.0, 175.0, 200.0, 225.0, 250.0, 275.0, 300.0, 325.0, 350.0, 375.0, 400.0, 425.0, 450.0, 475.0, 500.0, 550.0, 600.0, 650.0, 700.0, 750.0, 800.0, 850.0, 900.0, 950.0, 1000.0, 1050.0, 1100.0, 1150.0, 1200.0, 1250.0, 1300.0, 1350.0, 1400.0, 1450.0, 1500.0, 1550.0, 1600.0, 1650.0, 1700.0, 1750.0, 1800.0, 1850.0, 1900.0, 1950.0, 2000.0, 2100.0, 2200.0, 2300.0, 2400.0, 2500.0, 2600.0, 2700.0, 2800.0, 2900.0, 3000.0, 3100.0, 3200.0, 3300.0, 3400.0, 3500.0, 3600.0, 3700.0, 3800.0, 3900.0, 4000.0, 4100.0, 4200.0, 4300.0, 4400.0, 4500.0, 4600.0, 4700.0, 4800.0, 4900.0, 5000.0, 5100.0, 5200.0, 5300.0, 5400.0, 5500.0 m. Units: umol/l. The horizontal resolution of the produced DIVAnd analysis is 0.25 degrees.

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 1965-1970 until 2014-2019. Observational data span from 1960 to 2019. Depth range (IODE standard depths): 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, -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. Profiles were interpolated at standard depths using weighted parabolic interpolation algorithm (Reiniger and Ross, 1968). GEBCO 1min topography is used for the contouring preparation. Analysed filed masked using relative error threshold 0.3 and 0.5. DIVA settings: A constant value for signal-to-noise ratio was used equal to 1. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. 'log(data)-exp(analysis' 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. The entire set of related maps can be found in the viewing service: http://ec.oceanbrowser.net/emodnet/ . 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.

Moving 6-year analysis and visualization of Water body silicate in the North Sea. Four seasons (December-February, March-May, June-August, September-November). Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis) tool, version 2.7.9. results were subjected to the minfield option in DIVAnd to avoid negative/underestimated values in the interpolated results; error threshold masks L1 (0.3) and L2 (0.5) are included as well as the unmasked field. The depth dimension allows visualizing the gridded field at various depths.

Seasonal Climatology of Silicate for Loire River for the period 1965-2019 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Observational data span from 1965 to 2019. Depth levels (m): -125.0, -100.0, -75.0, -50.0,-40.0, -30.0, -25.0, -20.0, -15.0, -10.0, -8.0, -6.0, -4.0, -2.0, -0.0 Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Description of DIVAnd 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. The horizontal resolution of the produced DIVAnd maps grids is 0.01 degrees. 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. The weight of time series was decreased by a factor of 10 relative to the weight of the profiles to account for the redundancy in the time series observations. Detrending of data: no, Advection constraint applied: no. Units: umol/l.