Input of impulsive anthropogenic sound includes impulsive events from 2011-2016 • Seismic surveys (HELCOM-OSPAR Registry; national data call submissions as lines in the folder of data) • Explosions (HELCOM-OSPAR Registry) • Pile driving (HELCOM-OSPAR Registry) • Airguns (HELCOM-OSPAR Registry) For the different event types, numeric intensity value was used to represent the pressure as categorized in HELCOM-OSPAR Impulsive noise registry. All nationally reported seismic surveys were given intensity values “Very low” (0.25) - Very low (0.25) - Low (0.5) - Medium (0.75) - High (1) The impact distance has not been taken into account due to the different nature of separate datasets used for the pressure layer. We acknowledge that e.g. pile driving and airguns may impact up to 20 km from the source event. The spread of the sound wave depends on the sound frequency, water salinity, temperature and density.

Large shallow inlets bays (according to Habitats Directive Annex I) are large, shallow indentations of the coast, sheltered from wave action and where, in contrast to estuaries, the influence of freshwater is generally limited. The distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on GIS analysis and modelling, but also field inventories and ground-truthing has been carried out in some areas. Data coverage, accuracy and the methods in obtaining the data vary between countries.

Broad-scale habitat maps for the Baltic Sea have been produced in the EUSeaMap project in 2016. For German and Estonian marine areas, national (more accurate) datasets were used. German data included both substrate and light information (division into infralittoral/circalittoral). Estonian data included only substrate and the division into light regimes was obtained from the EuSeaMap data. Here, the habitat class “circalittoral mixed substrate” includes classes “mixed sediment” of the original data, in the circalittoral zone. The original polygon maps have been converted to 1 km x 1 km grid. The scale of the substrate data used in broad-scale habitat maps varies from 1:250 000 to 1:1M (data from EMODnet Geology). Coarser resolution data has been used in areas, where 1: 250 000 substrate data has not been available. Due to different scales used, the habitat classes may show different sized patterns in different areas.

Concentration of phosphorus pressure layer is interpolated from annual seasonal average of total phosphorus measurements from surface waters (0-10 m) extracted from ICES’s oceanographic database, database of Swedish Meteorological and Hydrological Institute, EEA’s Eionet database and Data from Gulf of Finland year 2014. The points were interpolated to cover the entire Baltic Sea with Spline with barriers interpolation method. Values were log-transformed and normalised (more detailed description below).

The extraction of cod pressure layer is based on two datasets: 1. http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/7a1389b3-382a-487f-8888-ac45c94c5a97 for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle). 2. http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/debeafcd-948b-4455-88ae-7a3d1618f5a8 from ICES recreational fisheries reports for 2011-2016, reported per country (only coastal areas included). Landing values were redistributed within each ICES rectangle by the c-square fishing effort data provided by ICES (all gears, 2011-2013). Tonnes / km² were calculated for both data sets and the results were converted to 1 km x 1 km grid cells. The layers were summed together, log-transformed and normalised to produce the final pressure layer on extraction of cod. Please see "lineage" section below for further details on attributes, data source, data processing, etc.

Distribution of Charophytes (Chara spp., Nitella spp., Nitellopsis spp., Tolypella spp.) mainly based on data submission by HELCOM contracting parties. Submitted point data was originally gathered in national mapping and monitoring campaigns, or for scientific research. Also scientific publications were used to complement the data (in Curonian, Vistula and Szczechin lagoons, see reference list). Polygon data from Poland was digitized based on Polish Marine Atlas. From Estonian waters, a predictive model was used (200m resolution), that was converted to presence/absence using minimized difference threshold (MDT) criteria. All data (points, polygon and the raster presenting predicted presence of Charophytes) were generalized to 5km x 5km grid cells.

The pressure oil slicks and spills is combination of following datasets: • Illegal oil discharges • Polluting ship accidents Illegal oil discharge data is based on airborne surveillance with remote sensing equipment in the Baltic Sea Area. The area of the detected spills in 2011–2016 was used to represent the pressure. The value of spills under 1km2 were directly given to grid cell, spills over 1km2 were buffered based on estimate spill area. For polluting ship accidents the reported oil spill volumes (m3) in years 2011-2015 were used for the pressure. Some polluting ship accidents spills were missing spilled oil volume, thus a mean of reported volumes was given to accidents with missing oil volume. Datasets were handled separately. Both layers were normalized, summed and normalized again to produce the “oil slicks and spills” pressure layer. Please see below for further details.

Broad-scale habitat maps for the Baltic Sea have been produced in the EUSeaMap project in 2016. For German and Estonian marine areas, national (more accurate) datasets were used. German data included both substrate and light information (division into infralittoral/circalittoral). Estonian data included only substrate and the division into light regimes was obtained from the EuSeaMap data. Here, the habitat class “circalittoral mud” includes classes “Fine mud”, “Sandy mud” and “Mud to sandy mud” of the original data, in the circalittoral zone. The original polygon maps have been converted to 1 km x 1 km grid. The scale of the substrate data used in broad-scale habitat maps varies from 1:250 000 to 1:1M (data from EMODnet Geology). Coarser resolution data has been used in areas, where 1: 250 000 substrate data has not been available. Due to different scales used, the habitat classes may show different sized patterns in different areas.

Distribution of blue mussel based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Mytilus spp. were submitted, originally gathered in national mapping and monitoring campaigns, or for scientific research. Point data from Poland was digitized based on Polish Marine Atlas. From Lithuania, a polygon delineating reefs was used to present Mytilus occurrence. For Germany, point data was complemented with a model describing Mytilus biomass in the German marine area (Darr et al. 2014), where predicted biomasses > 1g dw/ m2 were included as presence. From Estonian waters, a predictive model was used (200m resolution), that was converted to presence/absence using minimized difference threshold (MDT) criteria. All data (points, polygon and the raster presenting predicted presence of Mytilus) were generalized to 5km x 5km grid cells.

Distribution of Fucus sp. based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Fucus were submitted, originally gathered in national mapping and monitoring campaigns, or for scientific research purposes. From Estonian waters, a predictive model was used (200m resolution), that was converted to presence/absence using minimized difference threshold (MDT) criteria. All data (Fucus points and the raster presenting predicted presence of Fucus) were generalized to 5km x 5km grid cells.