Gridded product for visualization of Water body chlorophyll-a in the Western part of Black Sea generated by DIVA 4.6.11 using all EMODNET Chemistry data from 2000 to 2014. Depth range (IODE standard depths): 0, -5, -10, –20m. 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. Analysis: logaritmic transformation. 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: mg/m^3

Gridded product for visualization of Water body ammonium in the Black Sea - West Slope generated by DIVA 4.6.9 using all EMODNET Chemistry data from 1980 to 2013. Depth range (IODE standard depths): 0 – 75m. 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

Gridded product for visualization of Water body nitrate plus nitrite in the Black Sea generated by DIVA 4.6.9 using all EMODNET Chemistry data from 1976 to 2013. Depth range (IODE standard depths): 0 – 150m. 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

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 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

This visualization product displays plastic bags density per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during bottom trawl surveys. In cases where the wingspread and/or number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were calculated on each trawl and year using the following computation: Density of plastic bags (number of items per km²) = ∑Number of plastic bags related items / Swept area (km²) Percentiles 50, 75, 95 & 99 were calculated taking into account data for all years. The list of selected items for this product is attached to this metadata. Information on data processing and calculation is detailed in the attached methodology document. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area.

This visualization product displays the density of seafloor litter per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during bottom trawl surveys. In cases where the wingspread and/or the number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were calculated on each trawl and year using the following computation: Density (number of items per km²) = ∑Number of items / Swept area (km²) Percentiles 50, 75, 95 & 99 were calculated taking into account data for all years. More information on data processing and calculation are detailed in the attached document. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area.

This visualization product displays the spatial distribution of fishing related items density per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during bottom trawl surveys. In cases where the wingspread and/or number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were calculated on each trawl and year using the following computation: Density of fishing related items (number of items per km²) = ∑Number of fishing related items / Swept area (km²) Then a grid with 30km x 30km cells was used to calculate the weighted mean of densities in each cell from the formula : Weighted mean (number of items per km²) = ∑ (Distance (km) * Density (number of items per km²)) / ∑ Distance (km) Percentiles 50, 75, 95 & 99 were calculated taking into account data for all years. More information on data processing and calculation are detailed in the attached methodology document. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area. This work is based on the work presented in the following scientific article: O. Gerigny, M. Brun, M.C. Fabri, C. Tomasino, M. Le Moigne, A. Jadaud, F. Galgani, Seafloor litter from the continental shelf and canyons in French Mediterranean Water: Distribution, typologies and trends, Marine Pollution Bulletin, Volume 146, 2019, Pages 653-666, ISSN 0025-326X, https://doi.org/10.1016/j.marpolbul.2019.07.030.

This visualization product displays the spatial distribution of plastic bags density per trawl. EMODnet Chemistry included the collection of marine litter in its 3rd phase. Since the beginning of 2018, data of seafloor litter collected by international fish-trawl surveys have been gathered and processed in the EMODnet Chemistry Marine Litter Database (MLDB). The harmonization of all the data has been the most challenging task considering the heterogeneity of the data sources, sampling protocols (OSPAR and MEDITS protocols) and reference lists used on a European scale. Moreover, within the same protocol, different gear types are deployed during bottom trawl surveys. In cases where the wingspread and/or number of items were/was unknown, it was not possible to use the data because these fields are needed to calculate the density. Data collected before 2011 are concerned by this filter. When the distance reported in the data was null, it was calculated from: - the ground speed and the haul duration using the following formula: Distance (km) = Haul duration (h) * Ground speed (km/h); - the trawl coordinates if the ground speed and the haul duration were not filled in. The swept area was calculated from the wingspread (which depends on the fishing gear type) and the distance trawled: Swept area (km²) = Distance (km) * Wingspread (km) Densities were calculated on each trawl and year using the following computation: Density of plastic bags (number of items per km²) = ∑Number of plastic bags related items / Swept area (km²) Then a grid with 30km x 30km cells was used to calculate the weighted mean of densities in each cell from the formula : Weighted mean (number of items per km²) = ∑ (Distance (km) * Density (number of items per km²)) / ∑ Distance (km) Percentiles 50, 75, 95 & 99 were calculated taking into account data for all years. More information on data processing and calculation are detailed in the attached methodological document. Warning: the absence of data on the map does not necessarily mean that they do not exist, but that no information has been entered in the Marine Litter Database for this area. This work is based on the work presented in the following scientific article: O. Gerigny, M. Brun, M.C. Fabri, C. Tomasino, M. Le Moigne, A. Jadaud, F. Galgani, Seafloor litter from the continental shelf and canyons in French Mediterranean Water: Distribution, typologies and trends, Marine Pollution Bulletin, Volume 146, 2019, Pages 653-666, ISSN 0025-326X, https://doi.org/10.1016/j.marpolbul.2019.07.030.