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.

This pressure dataset is derived from three human activities datasets - Urban land use (on land) - Recreational boating and sports (updated layer for 2018 version, please see separate http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/8c30e828-1340-4162-b7f9-254586ae32b6) - Bathing sites These data are described in more detail in separate fact sheets. Urban land use data was first converted to 1 km grid cells and expanded with 1 km. Thus, coastal urban areas extended also to the sea. These areas were given value 1 and other sea areas, value 0. Bathing sites (points) were converted to 1km grid and given value 1, rest of the sea areas were given value 0. Normalized recreational boating data was converted to 1 km grid cells. These three layers were summed to produce the layer (values from 0 to 3), after that the layer was normalized. Hunting and recreational fishing data were excluded from human disturbance layer, as they are mostly reported per country and would have resulted in overestimation of the actual pressure.

The layer depicts the pressure of hazardous substances in the Baltic Sea, based on the data from the HOLAS 3 integrated hazardous substances assessment. The methodology utilizes the integrated status values available for each HELCOM assessment unit on level 3. The results are based on multiple hazardous substances groups integration, done through the CHASE tool. The integrated assessment assess the hazardous substances status in biota, water and sediment, and final result in based on the worst status. As the SPIA is carried out using a 1x1km grid and the Integrated hazardous substances is assessed on vector-based HELCOM assessment units, the vector data is rasterized. First, the vector data is rasterized to 100x100m resolution, and thereafter it is aggregated to 10x10km grid using a mean value. A 10 km grid is used in order to make the gradients between assessment units slightly smoother and finally values are converted to 1x1 km resolution. Please see "lineage" section below for further details.

Introduction of radionuclides is based on HELCOM MORS Discharge data from 2011 to 2014. The isotopes taken into account were: Cesium-137, Strontium-90, and Cobalt-60. The decay-corrected annual average of the sum of the radionuclide discharges (in Bq) were calculated for the pressure. 10 km buffer with linear decreasing function was used to represent the impact distance from the nuclear power plant outlets.

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.

Introduction of radionuclides is based on HELCOM MORS discharge data (2016-2020) . Annual averages of CO60, CS137 and SR90 from the period 2016-2020 per nuclear power plant. Gradual buffer around outlet to 10km distance (Type B decline). 10 km buffer with linear decline composed of 5 rings from discharges of radioactive substances (Type B decline)12.

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.

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.

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: Land claim - Area of polygon or 50 m buffer for points, 30m buffer for lines. Area of polygon - buffered line or point data, equals lost area. Watercourse modification - 50 m buffer. Area of polygon, buffered line or point data, equals lost area. Coastal defence and flood protection - 50 m buffer for lines, area of polygon. Area of polygon, buffered line or point data, equals lost area. Extraction of sand and gravel - Area of polygon. Area of polygon equals lost area. Dredging (capital) - Area of polygon or a 25/50 m buffer for <5000 m3 / >5000m3 sites. Area of polygon, buffered line or point data, equals lost area. Oil platforms - 25 m buffer. Buffered point data, equals lost area. Pipelines - 15 m buffer around cables with operational status. Area of polygon, buffered line or point data, equals lost area. Wind farms - 30 m buffer around each turbine with operational status. Buffered point data, equals lost area. Cables - 1.5 m buffer around cables with operational status. Buffered line data, equals lost area. Harbours - Polygon with 200 m buffer. Area of polygon, buffered line or point data, equals lost area. Marinas and leisure harbour - Point with 200 m buffer. Buffered point data, equals lost area. Bridges - 2 m buffer. Buffered line data, equals lost area. Finfish mariculture - 150 m buffer. Buffered point data, equals lost area. Shellfish mariculture - Area of polygon, 150 m buffer for points. Buffered point data, equals lost area. Activities are combined and potentially overlapping areas are removed. Dataset is clipped with coastline. Combined layer is intersected with 1 km grid to calculate % of area lost within a cell.