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 (48 parameters), and covers the Black Sea with 21504 CDI records divided per matrices: 7 biota profiles, 19677 water profiles, 1820 sediment profiles. Vertical profiles temporal range is from 1974-08-24 to 2017-10-06. Data were aggregated and quality controlled by ‘National Institute for Marine Research and Development "Grigore Antipa"’ from Romania. 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 48 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

Gridded product for visualization of Water body dissolved oxygen concentration in the Black Sea generated by DIVA 4.6.11 using all EMODNET Chemistry data from 1970 to 2014. Depth range (IODE standard depths): 0 – 250m. DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L), and both parameters vertically filtered. Background field: the data mean value is subtracted from the data. Detrending of data: no. Advection constraint applied: no. Every year of the time dimension corresponds to a 10-year centred average for each season: - winter: December - February, - spring: March - May, - summer: June - August, - autumn: September – November Units: umol/l

EMODnet Chemistry aims to provide access to marine chemistry datasets and derived data products concerning eutrophication, acidity and contaminants. The importance of the selected substances and other parameters relates to the Marine Strategy Framework Directive (MSFD). This aggregated dataset contains all unrestricted EMODnet Chemistry data on potential hazardous substances, despite the fact that some data might not be related to pollution (e.g. collected by deep corer). Temperature, salinity and additional parameters are included when available. It covers the Black Sea. Data were harmonised and validated by the ‘National Institute for Marine Research and Development "Grigore Antipa" in Romania. The dataset contains water, sediment and biota profiles. The temporal coverage is 1974–2021 for water measurements, 1990–2022 for sediment measurements and 2008–2016 for biota measurements. Regional datasets concerning contaminants are automatically harvested and the resulting collections are harmonised and validated using ODV Software and following a common methodology for all sea regions ( https://doi.org/10.6092/8b52e8d7-dc92-4305-9337-7634a5cae3f4). Parameter names are based on P01 vocabulary, which relates to BODC Parameter Usage Vocabulary and is available at: https://vocab.nerc.ac.uk/search_nvs/P01/. The harmonised dataset can be downloaded as as an ODV spreadsheet, which is composed of a metadata header followed by tab separated values. This spreadsheet can be imported into ODV Software for visualisation (more information can be found at: https://www.seadatanet.org/Software/ODV). In addition, the same dataset is offered also as a txt 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 into subcomponents (substance, CAS number, matrix...).This transposed format is more adapted to worksheet applications (e.g. LibreOffice Calc).

EMODnet Chemistry aims to provide access to marine chemistry datasets and derived data products concerning eutrophication, acidity and contaminants. The importance of the selected substances and other parameters relates to the Marine Strategy Framework Directive (MSFD). This aggregated dataset contains all unrestricted EMODnet Chemistry data on eutrophication and acidity, and covers the Black Sea, Marmara Sea and Sea of Azov. Data were aggregated and quality controlled by 'National Institute for Marine Research and Development 'Grigore Antipa' in Romania. ITS-90 water temperature and water body salinity variables have also been included ('as are') to complete the eutrophication and acidity data. If you use these variables for calculations, please refer to SeaDataNet for the quality flags: https://www.seadatanet.org/Products/Aggregated-datasets. Regional datasets concerning eutrophication and acidity are automatically harvested, and the resulting collections are aggregated and quality controlled using ODV Software and following a common methodology for all sea regions ( https://doi.org/10.13120/8xm0-5m67). Parameter names are based on P35 vocabulary, which relates to EMODnet Chemistry aggregated parameter names and is available at: https://vocab.nerc.ac.uk/search_nvs/P35/. When not present in original data, water body nitrate plus nitrite was calculated by summing all nitrate and nitrite parameters. The same procedure was applied for water body dissolved inorganic nitrogen (DIN), which was calculated by summing all nitrate, nitrite, and ammonium parameters. Concentrations per unit mass were converted to a unit volume using a constant density of 1.25 kg/L. The aggregated dataset can also be downloaded as an ODV collection and spreadsheet, which is composed of a 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).

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.2, 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.2, 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.2, 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_silicate in the Black Sea for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: from 1970-1975 until 2013-2018 Depth range (IODE standard depths): -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. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Log transformation was used in the analyses. 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/

Moving 6-year analysis of Water_body_phosphate in the Black Sea for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: from 1970-1975 until 2013-2018 Depth range (IODE standard depths): -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. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Log transformation was used in the analyses. Background field: the data mean value is subtracted from the data. Detrending of data: no, Advection constraint applied: no. Units: umol/l.

Moving 6-year analysis of Water_body_dissolved_oxygen_concentration in the Black Sea for each season: spring (March-May), summer (June-August), autumn (September-November), winter (December-February). Every year of the time dimension corresponds to the 6-years centred average of each season. 6-years periods span: from 1970-1975 until 2015-2020 Depth range (IODE standard depths): -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. Analysed field masked using relative error threshold 0.3 and 0.5 DIVA settings: signal-to-noise ratio and correlation length were estimated using data mean distance as a minimum (for L) and vertically filtered. Log transformation was used in the analyses. 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: https://emodnet.ec.europa.eu/geoviewer/