Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7.

Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7

Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7.

Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7.

Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7

Units: umol/l. Method: spatial interpolation produced with DIVA (Data-Interpolating Variational Analysis). URL: http://modb.oce.ulg.ac.be/DIVA. Comment: Every year of the time dimension corresponds to a 10-year centred average for each season : - winter season (December-February), - spring (March-May), - summer (June-August), - autumn (September-November). Diva settings: Snr=1.0, CL=0.7.

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.

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.

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.