Distribution of Furcellaria lumbricalis based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Furcellaria 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. For Poland, only confirmed occurrence of Furcellaria were included (Slupsk bansk, Rowy reef and reef at Orlowo cliff). All data (Furcellaria points and the raster presenting predicted presence of Furcellaria) were generalized to 5km x 5km grid cells.

The dataset contains total landings of sprat for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle) under EU Joint Research Centre’s data collection framework for fisheries data. Russian data extracted from ICES annual reports.

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.

Amount of hunted birds (number of birds/area) per year per area (county) is given separately for each target species: common scooter (Melanitta nigra), velvet scoter (Melanitta fusca), eider (Somateri molissima) and long tailed duck (Clangula hymalis). The data was made available by HELCOM Contracting Parties in response to data request. The data was received from Denmark, Estonia, Finland and Sweden. The activity was declared as not relevant in Germany, Latvia, Lithuania and Poland. For each species, a total number of hunted birds during the time period and a calculated average (hunted birds/year), is given. Data includes a total number (sum) of all hunted birds during the time period per county (total number of hunted birds/ county) and an average for hunted birds annually (hunted individuals/year). Velvet scoter is protected species in Sweden and Finland, and not listed as a game in Estonia. Common scoter is also protected species in Finland, thus hunting data is not available. Attribute specification and units: Country: Country AreaCode: County’s national code Area: County, unit area TOTAL: Total number of hunted birds in 2011-2015 Average: An average of hunted birds in a year (hunted birds/year) 2011_Sco – 2015_Sco: Number of hunted common scoters in 2011-2015 SUM_Sco: Total number of hunted common scoters in 2011-2015 Mean_Sco: An average number of hunted common scoters in a year (hunted individuals/year) 2011_VSco – 2015_VSco: Number of hunted velvet scoters in 2011 - 2015 SUM_Vsco: Total number of hunted velvet scoters in 2011-2015 Mean_Vsco: An average number of hunted velvet scoters in a year (hunted individuals/year) 2011_Eider – 2015_Eider: Number of hunted eiders in 2011 - 2015 SUM_Eider: Total number of hunted eiders in 2011-2015 Mean_Eider: An average number of hunted eiders in a year (hunted individuals/year) 2011_LTDuc – 2015_LTDuc: Number of hunted long tailed ducks in 2011 – 2015 SUM_LTDuck: Total number of hunted long tailed ducks in 2011-2015 Mean_LTDuc: An average number of hunted long tailed ducks in a year (hunted individuals/year) Notes: Notes regarding the data

Pressure layer combines all human activities that cause changes to hydrological conditions. The human activities were presented as point data which were given spatial extents (given below). The pressure value was given as the proportion of the grid cell under the pressure. The following human activities were combined into the changes to hydrological conditions layer; - Hydropower dams (a 1km2 grid cell in the river estuary was selected) - Water course modification (1 km) - Wind turbines (operational, 0.3 km, linear decline) - Oil platforms (0.5 km, linear decline) The human activity datasets were first processed separately covering the whole Baltic Sea and then summed together and overlapping areas were dissolved to remove double counting. Attenuation gradients are assigned to each layer as described above. Area effected decreases when distance from avtivity increases. Layer was normalized.

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.

The dataset contains total landings of cod for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle) under EU Joint Research Centre’s data collection framework for fisheries data. Russian data extracted from ICES annual reports.

Physical loss pressure layer combines all human activities that cause physical loss of seabed. The pressure is given as area lost in each cell (km2). For the polygon datasets the area was assumed to be the lost area. For line and point datasets spatial extents were calculated with buffers (below in brackets). If no buffer extent is indicated, the data was reported as polygon. The human activities used for the physical loss pressure: - Bridges (2 m) - Cables (operational; 1,5 m) - Coastal defence and flood protection (area of polygon, 50 m for lines) - Dredging (capital dredging, Area of polygon or a 25/50 m buffer for <5000 m3 / >5000m3 points) - Extraction of sand and gravel - Finfish mariculture (150 m) - Harbours (polygon with 200 m buffer) - Land claim (area of polygon, 30m buffer for lines) - Marinas and leisure harbours (200 m) - Oil platforms (25 m) - Oil terminals and refineries (200 m) - Pipelines (operational; 15 m) - Shellfish mariculture (area of polygon, 150 m points) - Watercourse modification (50 m) - Wind turbines (operational; 30m point location of turbine) The datasets were first processed separately covering the whole Baltic Sea and then merged into one uniform data layer and minimizing the effect of overlapping areas. Polygon areas were clipped with coastline to remove buffered areas that reached to land.

Concentration of nitrogen pressure layer is interpolated from annual seasonal average of total nitrogen concentrations 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).

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