This dataset is built from following Human activities datasets: • http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/81c917ea-492d-48e2-9f00-e1bb7fe3e4fc • http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/4fcd51dd-b8be-4e83-8cad-37c566782e8f The game hunting of seabirds data (see separate metadata): The total number of hunted seabirds were averaged over 2011-2015 (number of hunted seabirds / year). The area of the reporting unit was used to calculate the number of hunted seabirds / km2 and the data was converted to 1km x 1km grid. The predator control of seabirds data (see separate metadata): The total number of hunted cormorants were averaged over 2011-2015 (number of hunted cormorants / year). The area of the reporting unit was used to calculate the number of hunted cormorants / km2 and the data was converted to 1km x 1km grid. The two datasets were first separately log transformed and then summed, to get the total value for each grid cell. Zero values were given to all grid cells with no reported seabird hunting activity. The layer was normalized.

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

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.

Estuaries (according to Habitats Directive Annex I) are coastal inlets that are strongly influenced by freshwater. The distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on modelling, GIS analysis and/or aerial photos. Data coverage, accuracy and the methods in obtaining the data vary between countries.

This map presents the Special Protection Areas (SPAs) with reported breeding areas for birds. The spatial data on Special Protection Areas were gathered from the HELCOM contracting parties by Lund University, Sweden. In the data, the countries also indicated whether the sites were designated mainly due to wintering or breeding birds in the area. For Denmark, the information was obtained from standard forms for Natura 2000 sites. For Denmark, the data was updated after review process 20 February 2017. For Germany, the areas that were reported as “NA”(=information not available) were included in both breeding and wintering area maps. Many of the SPAs are both wintering and breeding areas. For the Baltic Sea Impact Index, the data was converted to 1 km x 1km grid cells.

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 map of sprat relative abundance is mainly based on Baltic International acoustic surveys (BIAS), years 2011-2016, (ICES WGBIFS reports), reported as millions of sprat per ICES rectangle. The BIAS surveys cover almost the whole area where sprat is commonly encountered. Outside BIAS area, sprat landings data was used to complement the data. For ICES rectangles surveyed by BIAS, values shown are the mean values per ICES rectangle based on BIAS data, average for 2011-2016. For ICES rectangles not surveyed by BIAS, values are calculated as: MAX-value x Weighting factor. The weighting factor is specific to each ICES rectangle, calculated as the ratio between the commercial landings in that rectangle and the commercial landings in the ICES rectangle with highest landings (based on averages for 2011-2015). MAX-value = millions of sprat according to BIAS in the ICES rectangle with highest landings. ICES rectangles outside the BIAS survey area with no reported sprat landings were given the value 0. The abundance values / ICES rectangle were divided by the area of the rectangle to obtain values per 1km2, and then converted to 1 km x 1km grid cells. Values were first log transformed and then normalised.

'Availability of deep water habitat, based on occurrence of H2S' layer describes the suitability of the bottom areas for the Baltic Sea biota, with regard to oxygen conditions of the near bottom waters. The data used to produce the layer was received from Leibniz-Institut für Ostseeforschung Warnemünde (IOW): - areas (polygons) with hydrogen sulfide (H2S) based on point measurements and modelling. Five time periods / year, for years 2011-2016 (altogether 30 layers). The polygons were converted to raster layers in a way, that for each time period (6 years, 5 time periods each year), areas with H2S got a value 0, other areas got the value 1. All layers were summed, (representing 6 years, 5 time periods each year, maximum value 30) and data was normalised. For more detailed information on the data used, please see Feistel et al. 2016.

Pressure layer combines all human activities that cause physical disturbance or damage to seabed. For several human activity datasets, spatial extents were given (table below). Buffers with decreasing value rates were applied to represent the impact distance of physical disturbance. The following human activities were combined into the physical disturbance layer; - Cables (under construction, 1 km buffer) - Coastal defence and flood protection (under construction, 500 m buffer) - Deposit of dredged material (500 m buffer for points and areas) - Dredging (maintenance) (500 m buffer for points and areas) - Extraction of sand and gravel (500 m buffer) - Finfish mariculture (1 km buffer) - Fishing intensity 2011-2016 average (subsurface swept area ratio) - Furcellaria harvesting - Pipelines (0,3 km buffer) - Recreational boating and sports - Shellfish mariculture - Shipping density - Wind farms (under construction) (1 km buffer) - Wind farms (operational) (0,1 km buffer) The human activity data sets were first processed separately covering the whole Baltic Sea and then summed together. In this integration, some data layers were down-weighted to arrive at a balanced pressure layer, as described below. High pressure intensity and/or slow recovery (weighting factor 1): Coastal defence and flood protection, Deposit of dredged material, Dredging, Extraction of sand and gravel and Fishing intensity Moderate to high (Weighting factor 0,8): Pipelines and Shipping density Moderate (Weighting factor 0,6): Finfish mariculture, Shellfish mariculture and Wind farms (under construction) Low to moderate (Weighting factor 0,4): Cables Low (Weighting factor 0,2): Maerl and Furcellaria harvesting, Recreational boating and sports and Wind farms (operational) Harbours and marinas were left out from the physical disturbance pressure to avoid double counting due to their representation in the shipping density and recreational boating and sports data sets.

This dataset shows sea bottom risk areas for mines sunk in the World War II. The big areas in Danish and German areas as well as in the Gdansk Bay are British flight mine areas. This dataset was created by the HELCOM Expert Group on Environmental Risks of Hazardous Submerged Objects (SUBMERGED). SUBMERGED works to compile and assess information about all kinds of hazardous objects and assess the associated risks. The dataset was provided by Gunnar Möller (Mine Warfare Data Center (C MWDC), 4th Naval Warfare Flottilla, Berga, Sweden) for the HELCOM Maritime Assessment published in 2018.