Combined product of Water body silicate based on DIVA 4D 10-year analysis on five regions : Northeast Atlantic Ocean, North Sea, Baltic Sea, Black Sea, Mediterranean Sea. The boundaries and overlapping zones between these five regions were filtered to avoid any unrealistic spatial discontinuities. This combined water body silicate product is masked using the relative error threshold 0.5. Units: umol/l

The analysis was performed per season using DIVA software tool (Data-Interpolating Variational Analysis). The analyses products are stored as NetCDF CF files and made available as WMS layers for easy browsing and adding. Every step of the time dimension corresponds to a 6-year moving average from 1983 to 2016. The depth dimension spans from surface to 1000 m, with 21 vertical levels. The boundaries and overlapping zones between these regions were filtered to avoid any unrealistic spatial discontinuities. This combined water body dissolved oxygen concentration product is masked using the relative error threshold 0.5. Units: µmol/l Created by 'University of Liège, GeoHydrodynamics and Environment Research (ULiège-GHER)'. The data used as input for DIVA have been extracted from the EMODnet Chemistry Download Service: https://emodnet-chemistry.maris.nl/search Intermediate regional data products: Mediterranean Sea - DIVA 4D 6-year analysis of Water body silicate 1965/2016 v2018, Arctic Ocean - DIVA 4D 6-year analysis of Water body silicate 1980/2017 v2018, North Sea - DIVA 4D 6-year analysis of Water body silicate 1980/2017 v2018, Black Sea - DIVA 4D 6-year analysis of Water body silicate 1990/2016 v2018, North East Atlantic Ocean - DIVA 4D 6-year analysis of Water body silicate 1971/2017 v2018

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

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.

This gridded product visualizes 1960 - 2014 water body silicate concentration (umol/l) in the North Sea domain, for each season (winter: December – February; spring: March – May; summer: June – August; autumn: September – November). It is produced as a Diva 4D analysis, version 4.6.9: a reference field of all seasonal data between 1960-2014 was used; results were logit transformed 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. Every step of the time dimension corresponds to a 10-year moving average for each season. The depth dimension allows visualizing the gridded field at various depths.

Moving 10-years analysis of Silicate at Northeast Atlantic Ocean for each season. - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December Every year of the time dimension corresponds to the 10-year centred average of each season. Decades span - from 1963-1972 until 2005-2014 (winter) - from 1979-1988 until 2005-2014 (spring) - from 1990-1999 until 2005-2014 (summer) - from 1964-1973 until 2005-2014 (autumn) Observational data span from 1962 to 2014. Depth range (IODE standard depths): -3000.0, -2500.0, -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. Signal to noise ratio and correlation length were optimized and filtered vertically and a seasonally-averaged profile was used. 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 consraint applied: no. Units: umol/l

Moving 10-years analysis of silicates at 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 10-year centered average of each season. Decades span from 1960-1969 until 2004-2013. Observational data span from 1960 to 2013. 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. Analyzed 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 3. Correlation length was optimized and filtered vertically and a seasonally-averaged profile was used. 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. 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. Units: umol/l

EMODnet Chemistry aims to provide access to marine chemistry data sets and derived data products concerning eutrophication, acidity and contaminants. The chemicals chosen reflect importance to the Marine Strategy Framework Directive (MSFD). ITS-90 water temperature and Water body salinity variables have been also included (as-is) to complete the Eutrophication and Acidity data. If you use these variables for calculations, please refer to SeaDataNet for having the quality flags: https://www.seadatanet.org/Products/Aggregated-datasets . This aggregated dataset contains all unrestricted EMODnet Chemistry data on Eutrophication and Acidity (14 parameters with quality flag indicators), and covers the Mediterranean Sea with 67271 CDI records (60876 Vertical profiles and 6395 Time series). Vertical profiles temporal range is from 1911-08-17 to 2017-10-25. Time series temporal range is from 1974-06-04 to 2015-12-22. Data were aggregated and quality controlled by 'Hellenic Centre for Marine Research, Hellenic National Oceanographic Data Centre (HCMR/HNODC)' from Greece. Regional datasets concerning eutrophication and acidity are automatically harvested and resulting collections are aggregated and quality controlled using ODV Software and following a common methodology for all Sea Regions ( https://doi.org/10.6092/9f75ad8a-ca32-4a72-bf69-167119b2cc12). When not present in original data, Water body nitrate plus nitrite was calculated by summing up the Nitrates and Nitrites. Same procedure was applied for Water body dissolved inorganic nitrogen (DIN) which was calculated by summing up the Nitrates, Nitrites and Ammonium. Parameter names are based on P35, EMODnet Chemistry aggregated parameter names vocabulary, which is available at: https://www.bodc.ac.uk/resources/vocabularies/vocabulary_search/P35/. Detailed documentation is available at: https://doi.org/10.6092/ec8207ef-ed81-4ee5-bf48-e26ff16bf02e The aggregated dataset can be downloaded as ODV spreadsheet, which is composed of metadata header followed by tab separated values. This worksheet can be imported to ODV Software for visualisation (More information can be found at: https://www.seadatanet.org/Software/ODV ). The original datasets can be searched and downloaded from EMODnet Chemistry Download Service: https://emodnet-chemistry.maris.nl/search

EMODnet Chemistry aims to provide access to marine chemistry data sets and derived data products concerning eutrophication, acidity and contaminants. The chemicals chosen reflect importance to the Marine Strategy Framework Directive (MSFD). ITS-90 water temperature and Water body salinity variables have been also included (as-is) to complete the Eutrophication and Acidity data. If you use these variables for calculations, please refer to SeaDataNet for having the quality flags: https://www.seadatanet.org/Products/Aggregated-datasets . This aggregated dataset contains all unrestricted EMODnet Chemistry data on Eutrophication and Acidity (15 parameters with quality flag indicators), and covers the Mediterranean Sea with 187722 CDI records (180932 Vertical profiles and 6790 Time series). Vertical profiles temporal range is from 1911-08-17 to 2020-10-16. Time series temporal range is from 1974-06-04 to 2018-02-20. Data were aggregated and quality controlled by 'Hellenic Centre for Marine Research, Hellenic National Oceanographic Data Centre (HCMR/HNODC)' from Greece. Regional datasets concerning eutrophication and acidity are automatically harvested and resulting collections are aggregated and quality controlled using ODV Software and following a common methodology for all Sea Regions ( https://doi.org/10.6092/9f75ad8a-ca32-4a72-bf69-167119b2cc12). When not present in original data, Water body nitrate plus nitrite was calculated by summing up the Nitrates and Nitrites. Same procedure was applied for Water body dissolved inorganic nitrogen (DIN) which was calculated by summing up the Nitrates, Nitrites and Ammonium. Parameter names are based on P35, EMODnet Chemistry aggregated parameter names vocabulary, which is available at: https://www.bodc.ac.uk/resources/vocabularies/vocabulary_search/P35/. Detailed documentation is available at: https://dx.doi.org/10.6092/4e85717a-a2c9-454d-ba0d-30b89f742713 Explore and extract data at: https://emodnet-chemistry.webodv.awi.de/eutrophication%3EMediterranean The aggregated dataset can also be downloaded as ODV collection and spreadsheet, which is composed of metadata header followed by tab separated values. This spreadsheet can be imported to ODV Software for visualisation (More information can be found at: https://www.seadatanet.org/Software/ODV ). The original datasets can be searched and downloaded from EMODnet Chemistry Chemistry CDI Data and Discovery Access Service: https://emodnet-chemistry.maris.nl/search