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 hard substrate” includes classes “Rock and other hard substrate” and “Coarse substrate” 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.

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 “infralittoral hard substrate” includes classes “Rock and other hard substrate” and “Coarse substrate” of the original data, in the infralittoral 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.

Introduction of radionuclides is based on HELCOM MORS Discharge data from 2011 to 2014. The isotopes taken into account were: Cesium-137, Strontium-90, and Cobalt-60. The decay-corrected annual average of the sum of the radionuclide discharges (in Bq) were calculated for the pressure. 10 km buffer with linear decreasing function was used to represent the impact distance from the nuclear power plant outlets.

The occurrence of suitable nursery habitats is crucial for maintaining fish populations (Sundblad et al. 2013). For perch, species distribution modelling studies (Snickars et al. 2010, Bergström et al. 2013, Sundblad et al. 2013) have shown the importance of suitable environmental conditions for reproduction. Due to lack of coherent data on perch spawning and nursery areas across the Baltic Sea countries, environmental variables were used in delineating potential recruitment areas for perch. The distribution area or perch recruitment is delineated by selecting areas where depth < 4 m (For Danish waters < 3 m), logged exposure < 5 (exposure model described in Isæus 2004), and salinity < 10 PSU. The threshold values have been obtained from literature (Snickars et al. 2010, Bergström et al. 2013, Skovrind et al. 2013, Sundblad et al. 2013). Relatively “loose” thresholds have been used, to rather overestimate than underestimate the recruitment area (precautionary approach). Along the Finnish coastline a national model has been used (Kallasvuo et al. 2016), with suitable environments for perch recruitment generalized to 1 km x 1 km grid.

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.

Potential cumulative impacts on benthic habitats is based on the same method than http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/9477be37-94a9-4201-824a-f079bc27d097, but is focused on physical pressures and benthic habitats. The dataset was created based on separate analysis for potential cumulative impacts on only the benthic habitats, as these are particularly affected by physical pressures. In this case the evaluation was based on pressure layers representing http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/ea0ef0fa-0517-40a9-866a-ce22b8948c88 and http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/05e325f3-bc30-44a0-8f0b-995464011c82, combined with information on the distribution of eight broad benthic habitat types and five habitat-forming species (http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/363cb353-46da-43f4-9906-7324738fe2c3, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/f9cc7b2c-4080-4b19-8c38-cac87955cb91, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/264ed572-403c-43bd-9707-345de8b9503c, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/822ddece-d96a-4036-9ad8-c4b599776eca and http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/ca327bb1-d3cb-46c2-8316-f5f62f889090). The potential cumulative impacts has been estimated based on currently best available data, but spatial and temporal gaps may occur in underlying datasets. Please scroll down to "Lineage" and visit http://stateofthebalticsea.helcom.fi/cumulative-impacts/ for more info.

The Baltic Sea Pressure Index is a calculation of quantity and spatial distribution of potential cumulative impacts on the Baltic Sea. The BSII is based on georeferenced datasets of human activities (36 datasets), pressures (18 datasets) and ecosystem components (36 datasets), and on sensitivity estimates of ecosystem components (so-called sensitivity scores) that combine the pressure and ecosystem component layers, created in http://www.helcom.fi/helcom-at-work/projects/holas-ii project. The assessment can be applied with a focus on pressures only by using the Baltic Sea Pressure Index (BSPI) which shows the anthropogenic pressures/human activities in the defined assessment units without including ecosystem components. The BSPI however includes a weighting component in order to grade the effect of the pressures on the ecosystem in a generalized perspective. Cumulative impacts are calculated for each assessment unit (1 km2 grid cells) by summing all impacts occurring in the unit. All datasets and methodologies used in the index calculations are approved by all HELCOM Contracting Parties in review and acceptance processes. This dataset covers the time period 2011-2016. Please scroll down to "Lineage" and visit http://stateofthebalticsea.helcom.fi/cumulative-impacts/ for more info.

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

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.