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

Distribution of Fucus sp. based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Fucus 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. All data (Fucus points and the raster presenting predicted presence of Fucus) were generalized to 5km x 5km grid cells.

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.

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.

This map presents the Special Protection Areas (SPAs) with reported breeding areas for birds. The spatial data on Special Protection Areas were gathered from the HELCOM contracting parties by Lund University, Sweden. In the data, the countries also indicated whether the sites were designated mainly due to wintering or breeding birds in the area. For Denmark, the information was obtained from standard forms for Natura 2000 sites. For Denmark, the data was updated after review process 20 February 2017. For Germany, the areas that were reported as “NA”(=information not available) were included in both breeding and wintering area maps. Many of the SPAs are both wintering and breeding areas. For the Baltic Sea Impact Index, the data was converted to 1 km x 1km 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.

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.

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