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.

Lagoons are expanses of shallow coastal waters, wholly or partially separated from the sea by sandbanks or shingle, or by rocks. Salinity may vary from brackish water to hypersalinity depending on rainfall, evaporation and addition of fresh seawater from storms, temporary flooding, or tidal exchange. The distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on modelling and/or GIS analysis. 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 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.

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.

Summary The following 6 categories of annual mean salinity were applied delineating the Kattegat and the Baltic Sea into regions with differences in salinity regime (fig. 15): I. Oligohaline I (< 5psu). II. Oligohaline II (5 - 7.5psu). III. Mesohaline I (7.5 - 11psu). IV. Mesohaline II (11 - 18psu). V. Polyhaline (18 - 30psu). VI. Euhaline (>30psu). Description This dataset was produced by NERI, Denmark, for the BSR INTERREG IIIB project BALANCE. Due to the stratification in the Baltic Sea it was decided to use bottom salinity for the development of the benthic marine landscapes and difference in surface to bottom salinity for the pelagic landscapes. The following 6 categories of annual mean salinity were applied delineating the Kattegat and the Baltic Sea into regions with differences in salinity regime (fig. 15): I. Oligohaline I (< 5psu). II. Oligohaline II (5 - 7.5psu). III. Mesohaline I (7.5 - 11psu). IV. Mesohaline II (11 - 18psu). V. Polyhaline (18 - 30psu). VI. Euhaline (>30psu).

The pressure layer represents biological pressure caused by introduction of non-indigenous species. The data is obtained from core indicator Trend in the arrival of new non-indigenous species (BSEP 129b: http://www.helcom.fi/Lists/Publications/BSEP129B.pdf). For the Baltic Sea Impact Index, the 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.

Baltic International Trawl Survey (BITS) data (2011-2016) from ICES DATRAS database was used as a base to create a map of cod relative abundance (quarter 1 data, CPUE values per ICES subdivision). Cod = 30cm was included. For ICES rectangles surveyed by BITS, values shown are the mean CPUE per ICES subdivision based on BITS data, average for 2011-2016. For ICES rectangles not surveyed by BITS, 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-2016). MAX-value = CPUE according to BITS in the ICES rectangle with highest landings. ICES rectangles outside the BITS survey area with no reported cod landings were given the value 0. Values were first log transformed and then normalized.

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.

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.