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.

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

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.

Concentration of phosphorus pressure layer is interpolated from annual seasonal average of total phosphorus measurements 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).

This layer is based on data from the BIAS project representing ambient underwater noise, modelled into a 0.5 km x 0.5 km grid, and representing sound pressure levels at 1/3 octave bands of 125 Hz exceeded at least 5% of the time. Measured and modelled acoustic data is provided as Sound Pressure Level (SPL). The time period for the data is annual values for year 2014. The selected depth interval is 0 m – bottom to represent the ambient underwater noise in the whole water column. The data were normalized setting level 0 at 92 db re 1µPa and level 1 at 127 db re 1µPa.

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 mud” includes classes “Fine mud”, “Sandy mud” and “Mud to sandy mud” 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 “circalittoral sand includes classes “Sand” and “Muddy sand” of the original data, in the circalittoral 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.

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.

Distribution of Furcellaria lumbricalis based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Furcellaria were submitted, originally gathered in national mapping and monitoring campaigns, or for scientific research purposes. From Estonian waters, a predictive model was used (200m resolution), that was converted to presence/absence using minimized difference threshold (MDT) criteria. For Poland, only confirmed occurrence of Furcellaria were included (Slupsk bansk, Rowy reef and reef at Orlowo cliff). All data (Furcellaria points and the raster presenting predicted presence of Furcellaria) were generalized to 5km x 5km grid cells.

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.