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.

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.

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.

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.

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.

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.

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.

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.

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.

Dataset represents the radioactive discharges from nuclear facilities in the Baltic Sea area. Data includes isotopes CS137, CO60 and SR90 Aquatic discharges in 2011-2014 with decay corrections.