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.

Input of hazardous substances pressure layer is interpolated from CHASE Assessment tool concentration component. The contamination ratio values were calculated with CHASE Assessment tool for hazardous substances monitored in water, sediment and biota. Classified mean contamination ratio was used in the interpolation. Classification is based on the http://stateofthebalticsea.helcom.fi/about-helcom-and-the-assessment/downloads-and-data/. The points were interpolated to cover the entire Baltic Sea with Spline with barriers interpolation method. Please see "lineage" section below for further details on attributes, data source, data processing, etc.

This dataset is built from the following Human activities dataset: • Hunting of seals The number of hunted seals (see separate metadata on hunting of seals) were averaged over 2011-2014 separately for grey seals, ringed seals and harbour seals (e.g. number of hunted grey seals / year). In Sweden the numbers of hunted grey seals in 2011 (74) were reported for the whole Swedish territorial waters), but here the numbers were set only to Swedish Gulf of Bothnia, as corresponding numbers were reported there in 2013 (75) and 2014 (65). The area of the reporting unit was used to calculate the number of hunted seals / km2 and the data was converted to 1km x 1km grid. For the Baltic Sea Impact Index, the values were normalized. Normalized value 0.5 was set to the level of quota for hunting of seal species in the Baltic Sea. The following quotas for hunting were used: Grey seal: 2000, Ringed seal: 350, Harbour seal 230.

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.

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.

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.

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

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.

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.

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.