Running 6-year analysis of Water body dissolved oxygen concentration in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1960-2021. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis in n dimensions) tool. GEBCO_08 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the horizontal correlation length is fixed to 80 km and decreasing towards the coastline, the vertical correlation length is varying with depth. Signal to noise ratio is fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.

Moving 6-year analysis of Water body dissolved oxygen concentration in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1980-1985 and 2013-2018. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. 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 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the Correlation length in the horizontal direction were fixed to 80 km and in the vertical direction to 25 m. Signal to noise ratio were fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.

Running 6-year analysis of Water body phosphate in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1975-2021. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis in n dimensions) tool. GEBCO_08 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the horizontal correlation length is fixed to 80 km and decreasing towards the coastline, the vertical correlation length is varying with depth. Signal to noise ratio is fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is 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. This aggregated dataset contains all unrestricted EMODnet Chemistry data on Eutrophication and Acidity (14 parameters with quality flag indicators), and covers the Baltic Sea with 175834 CDI stations (175778 Vertical profiles and 56 Time series). Vertical profiles temporal range is from 1902-08-05 to 2020-10-10. Time series temporal range is from 2010-01-12 to 2016-02-10. Data were aggregated and quality controlled by "Swedish Meteorological and Hydrological Institute (SMHI)" from Sweden. 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%3EBaltic 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

Moving 6-year analysis and visualization of Water body dissolved oxygen concentration in the North Sea, four seasons: winter (dec-feb), spring (mar-may), summer (jun-aug), autumn (sep-nov). 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.2. 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.

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.

Running 6-year analysis of Water body dissolved inorganic nitrogen in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1975-2021. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. Data Sources: observational data from SeaDataNet/EMODnet Chemistry Data Network. Description of DIVA analysis: Geostatistical data analysis by DIVAnd (Data-Interpolating Variational Analysis in n dimensions) tool. GEBCO_08 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the horizontal correlation length is fixed to 80 km and decreasing towards the coastline, the vertical correlation length is varying with depth. Signal to noise ratio is fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.

Moving 6-year analysis of Water body phosphate in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1980-1985 and 2013-2018. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. 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 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the Correlation length in the horizontal direction were fixed to 80 km and in the vertical direction to 25 m. Signal to noise ratio were fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.

Moving 6-year analysis of Water body dissolved inorganic nitrogen (DIN) in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1980-1985 and 2013-2018. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. 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 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the Correlation length in the horizontal direction were fixed to 80 km and in the vertical direction to 25 m. Signal to noise ratio were fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.

Moving 6-year analysis of Water body silicate in the Baltic Sea. Four seasons (March-May, June-August, September-November, December-February). Every year of the time dimension corresponds to a 6-year centred average. Periods span between 1980-1985 and 2013-2018. Analyses for depths (m) (HELCOM standard depths): 0, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 225, 250, 275, 300. 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 Grid (30 arc-seconds) topography is used for the contouring preparation. Files contain analysed fields, error fields and combined field with the deepest value for each grid point selected. Also pre-masked fields using relative error threshold 0.3 and 0.5 are included. In the analyses the Correlation length in the horizontal direction were fixed to 80 km and in the vertical direction to 25 m. Signal to noise ratio were fixed to 1.0. Background fields were created using data for the given time period and season. Log transformation was used in the analyses. No detrending, advection constraints or weighting are applied. Unit is umol/l.