Input of heat pressure dataset contains delta heat values from warm water discharge of - Discharge of warm water from nuclear power plants - Fossil fuel energy production. The Discharge of warm water from nuclear power plants was requested from HELCOM Contracting Parties. Average temperature change of the cooling water (°C) and amount of cooling water (m3) was used to calculate the heat load in TWh. No data on heat load was available for the Leningrad nuclear power plant; therefore the average heat load (TWh) of discharge of warm water from nuclear power plants was given. No heat load information was available for fossil fuel energy production sites. Heat load of 1 (TWh) was given to all production sites, based on average heat load of individual production site in Helsinki from recent years. 1 km buffer was used for both datasets with steep decline around the outlet. Overlapping areas were summed. Heat load from both layers were summed and the 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.

The extraction of Sprat data set is based on: 1. http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/1fb1bd2d-8dff-493a-9ed3-a278aec8f371 for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle). Landing values were redistributed within each ICES rectangle by the c-square fishing effort data provided by ICES (all gears, 2011-2013). Tonnes / km² was calculated and the results were converted to 1 km x 1 km grid cells. The layer was log-transformed and normalised to produce the final pressure layer on extraction of Sprat. Please see "lineage" section below for further details on attributes, data source, data processing, etc.

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 occurrence of suitable nursery habitats is crucial for maintaining fish populations (Sundblad et al. 2013). Species distribution modelling studies have shown the importance of suitable environmental conditions for pikeperch recruitment. Due to lack of coherent data on pikeperch spawning and nursery areas across the Baltic Sea countries, environmental variables were used in delineating potential recruitment areas for pikeperch. The pikeperch recruitment area presented on the map is mainly delineated by selecting areas where depth < 5 m, logged exposure < 5, salinity < 7 PSU, Secchi depth < 2 m and distance to deep (10m) water < 4km. The threshold values have been obtained from literature (Veneranta et al. 2011, Bergström et al. 2013, Sundblad et al. 2013, Kallasvuo et al. 2016). Temperature, although important for pikeperch, was left out due to high variation in timing of suitable spawning temperatures across the Baltic Sea. In Finnish coastal waters, a national pikeperch model (Kallasvuo et al. 2016) has been used, with very suitable areas for pikeperch generalized to 1 km grid. In Sweden, the areas delineated by environmental variables have been complemented with information from national interview survey (Gunnartz et al. 2011) as well as expert opinion.

Boreal Baltic islets and small islands (according to Habitats Directive Annex I) are groups of skerries, islets or single small islands, mainly in the outer archipelago or offshore areas. They are important nesting sites for birds and resting sites for seals. The surrounding sublittoral vegetation is also included. The distribution map is based on data submission by HELCOM contracting parties. Only Sweden and Finland reported occurrences of boreal Baltic islets and small islands.

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.

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.

Submarine structures made by leaking gases (according to Habitats Directive Annex I) are also known as “bubbling reefs”. These formations support a zonation of diverse benthic communities consisting of algae and/or invertebrate specialists of hard marine substrates different to that of the surrounding habitat. The distribution map is based on data submission by HELCOM contracting parties. Only Sweden and Denmark reported occurrences of submarine structures made by leaking gases.

Reefs (according to Habitats Directive Annex I) are hard compact substrata (either biogenic or geogenic) on solid and soft bottoms, which arise from the seafloor in the sublittoral and littoral zone. Distribution of mapped Natura 2000 habitat “Reefs” based on data submission by HELCOM contracting parties. Most of the submitted data is based on modelling and limited ground-truthing. Data coverage, accuracy and the methods in obtaining the data vary between countries.