Seasonal climatology of Water body dissolved oxygen concentration for Loire river for the period 1950-2021 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1950 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): 150000 (winter), 100000 (spring), 70000 (summer), 115000 (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.

Seasonal climatology of Water body phosphate for Loire river for the period 1950-2021 and for the following seasons: - winter: January-March, - spring: April-June, - summer: July-September, - autumn: October-December. Observation data span from 1950 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): 230000 (winter), 264000 (spring), 140000 (summer), 135000 (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 Water body dissolved inorganic nitrogen in the NorthEast Atlantic 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 1971/1976 until 2016/2021. Observation data span from 1971 to 2021. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 400km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. 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: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

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 1970-1975 until 2017-2022. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Units: umol/l. 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. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 25 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.,1200.,1300.,1400.,1500.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2022) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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 Water body phosphate 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 1970-1975 until 2017-2022. Data Sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Units: umol/l. 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. The horizontal resolution of the produced DIVAnd maps grids is dx=dy=0.125 degrees (around 13.5km and 10.9km accordingly). The vertical resolution is 25 depth levels: [0.,5.,10.,20.,30.,50.,75.,100.,125.,150.,200.,250.,300.,400.,500.,600.,700.,800.,900.,1000.,1100.,1200.,1300.,1400.,1500.]. The horizontal correlation length is 200km. The vertical correlation length (in meters) was set twices the vertical resolution: [10.,10.,20.,20.,40.,50.,50.,50.,50.,100.,100.,100.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.,200.]. Duplicates check was performed using the following criteria for space and time: dlon=0.001deg., dlat=0.001deg., ddepth=1m, dtime=1hour, dvalue=0.1. The error variance (epsilon2) was set equal to 1 for profiles and 10 for time series to reduce the influence of close data near the coasts. An anamorphosis transformation was applied to the data (function DIVAnd.Anam.loglin) to avoid unrealistic negative values: threshold value=200. A background analysis field was used for all years (1970-2022) with correlation length equal to 600km and error variance (epsilon2) equal to 20. Quality control of the observations was applied using the interpolated field (QCMETHOD=3). Residuals (differences between the observations and the analysis (interpolated linearly to the location of the observations) were calculated. Observations with residuals outside the minimum and maximum values of the 99% quantile were discarded from the analysis. 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 Water body phosphate in the NorthEast Atlantic 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 1950/1955 until 2016/2021. Observation data span from 1950 to 2021. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 250km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. 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: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

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.

Moving 6-year analysis of Water body silicate in the NorthEast Atlantic 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 1950/1955 until 2016/2021. Observation data span from 1950 to 2021. Depth levels (IODE standard depths): [0.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 75.0, 100.0, 125.0, 150.0, 200.0, 250.0, 300.0, 400.0, 500.0, 600.0, 700.0, 800.0, 900.0, 1000.0, 1100.0, 1200.0, 1300.0, 1400.0, 1500.0, 1750.0, 2000.0]. Data sources: observational data from SeaDataNet/EMODNet Chemistry Data Network. Descrption of DIVAnd analysis: the computation was done with DIVAnd (Data-Interpolating Variational Analysis in n dimensions), version 2.7.4, using GEBCO 30 sec topography for the spatial connectivity of water masses. The horizontal resolution of the produced DIVAnd maps is 0.1 degrees. Horizontal correlation length varies from 300km in open sea regions to 50km at the coast. Vertical correlation length is defined as twice the vertical resolution. 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: a vertically-filtered profile of the seasonal data mean value (including all years) is substracted from the data. Detrending of data: no, advection constraint applied: no. Units: umol/l.

Moving 6-year analysis and visualization of Water body dissolved oxygen concentration 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.

Moving 6-year analysis and visualization of Water body phosphate 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.