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.

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.

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.

Distribution of Charophytes (Chara spp., Nitella spp., Nitellopsis spp., Tolypella spp.) mainly based on data submission by HELCOM contracting parties. Submitted point data was originally gathered in national mapping and monitoring campaigns, or for scientific research. Also scientific publications were used to complement the data (in Curonian, Vistula and Szczechin lagoons, see reference list). Polygon data from Poland was digitized based on Polish Marine Atlas. 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 Charophytes) were generalized to 5km x 5km grid cells.

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.

Sandbanks (according to Habitats Directive Annex I) are areas elevated from their surroundings that consist mainly of sand, but where cobbles and boulders can occur. Distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on modelling, GIS analysis and only limited ground-truthing has been carried out. Data coverage, accuracy and the methods in obtaining the data vary between countries.

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.

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 mixed substrate” includes classes “mixed sediment” of the original data, in the infralittoral zone. The original polygon maps have been converted to 1 km x 1km 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.

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.

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.