EMODnet Chemistry aims to provide access to marine chemistry data sets and derived data products concerning eutrophication, ocean acidification and contaminants. The chemicals chosen reflect importance to the Marine Strategy Framework Directive (MSFD). This regional aggregated dataset contains all unrestricted EMODnet Chemistry data on contaminants (20 parameters), and covers the Arctic Ocean with 548 CDI records divided per matrices: 235 water profiles, 313 sediment profiles. Vertical profiles temporal range is from 1974-03-25 to 2015-06-15. Data were aggregated and quality controlled by ‘Institute of Marine Research - Norwegian Marine Data Centre (NMD)’ from Norway. Regional datasets concerning contaminants are automatically harvested. Parameter names in these datasets are based on P01, BODC Parameter Usage Vocabulary, which is available at: http://seadatanet.maris2.nl/bandit/browse_step.php . Each measurement value has a quality flag indicator. The resulting data collections for each Sea Basin are harmonised, and the collections are quality controlled by EMODnet Chemistry Regional Leaders using ODV Software and following a common methodology for all Sea Regions. Harmonisation means that: (1) unit conversion is carried out to express contaminant concentrations with a limited set of measurement units (according to EU directives 2013/39/UE; Comm. Dec. EU 2017/848) and (2) merging of variables described by different “local names” ,but corresponding exactly to the same concepts in BODC P01 vocabulary. The harmonised dataset can be downloaded as ODV spreadsheet (TXT file), 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 same dataset is offered also as XLSX file in a long/vertical format, in which each P01 measurement is a record line. Additionally, there are a series of columns that split P01 terms in subcomponents (measure, substance, CAS number, matrix...).This transposed format is more adapted to worksheet applications users (e.g. LibreOffice Calc). The 20 parameter names in this metadata record are based on P02, SeaDataNet Parameter Discovery Vocabulary, which is available at: http://seadatanet.maris2.nl/v_bodc_vocab_v2/vocab_relations.asp?lib=P02 . Detailed documentation will be published soon. The original datasets can be searched and downloaded from EMODnet Chemistry Download Service: https://emodnet-chemistry.maris.nl/search

Moving 10-years analysis of pH 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 1959-1968 until 2005-2014. Observational data span from 1930 to 2014. 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 1. 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: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o 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 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o 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: pH units

Moving 10-years analysis of dissolved_oxygen_concentration 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 1971-1980 until 2006-2015. Observational data span from 1911 to 2015. Depth range (IODE standard depths):-3000.0, -2500.0, -2000.0, -1750.0, -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 1. 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: o Brunetti Fabio (OGS) o Cardin Vanessa, Bensi Manuel doi:10.6092/36728450-4296-4e6a-967d-d5b6da55f306 o Cardin Vanessa, Bensi Manuel, Ursella Laura, Siena Giuseppe doi:10.6092/f8e6d18e-f877-4aa5-a983-a03b06ccb987 o Cataletto Bruno (OGS) o Cinzia Comici Cinzia (OGS) o Civitarese Giuseppe (OGS) o DeVittor Cinzia (OGS) o Giani Michele (OGS) o Kovacevic Vedrana (OGS) o Mosetti Renzo (OGS) o 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 o Celio Massimo (ARPA FVG) o Malaguti Antonella (ENEA) o Fonda Umani Serena (UNITS) o Bignami Francesco (ISAC/CNR) o Boldrini Alfredo (ISMAR/CNR) o Marini Mauro (ISMAR/CNR) o Miserocchi Stefano (ISMAR/CNR) o Zaccone Renata (IAMC/CNR) o 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

Water body dissolved inorganic nitrogen (DIN) - 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.

Water body phosphate - 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.

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.

Water body silicate - 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 dissolved oxygen concentration 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, 2500.0, 3000.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 phosphate in the Arctic Ocean, for each season in the period 1970-2017. Every year of the time dimension corresponds to the 6-year centered average for each season. Winter: December-February, Spring: March-May, Summer: June-August, Autumn: September-November. Depth range (IODE standard depths): 0, 5, 10, 20, 30, 40, 50, 75, 100, 125, 150, 200, 250, 300, 400, ..., 1500, 1750, 2000, 2500m. Units: umol/l. Description of DIVA analysis: The computation was done with 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 0.1 degrees. Signal-to-noise ratio was fixed to 3.0, horizontal correlation length varying from 45 km near the coast to 150 km, and vertical correlation length varying between 25 and 1000 m. Logarithmic transformation is applied to the data prior to the analysis. Background field: analysis with signal-to-noise ratio = 10, horizontal correlation length 60-200 km, and vertical correlation length 25-1000 m.

Moving 6-year analysis of silicate in the Arctic Ocean, for each season in the period 1970-2017. Every year of the time dimension corresponds to the 6-year centered average for each season. Winter: December-February, Spring: March-May, Summer: June-August, Autumn: September-November. Depth range (IODE standard depths): 0, 5, 10, 20, 30, 40, 50, 75, 100, 125, 150, 200, 250, 300, 400, ..., 1500, 1750, 2000, 2500m. Units: umol/l. Description of DIVA analysis: The computation was done with 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 0.1 degrees. Signal-to-noise ratio was fixed to 3.0, horizontal correlation length varying from 45 km near the coast to 150 km, and vertical correlation length varying between 25 and 1000 m. Logarithmic transformation is applied to the data prior to the analysis. Background field: analysis with signal-to-noise ratio = 10, horizontal correlation length 60-200 km, and vertical correlation length 25-1000 m.