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.

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.

Potential cumulative impacts of eutrophication and hazardous substances assesses the cumulative potential effect of eutrophication and hazardous substances over all ecosystem components. The evaluation is based on the pressure layer on eutrophication and hazardous substances, combined with information on all ecosystem components (57 layers) included in SPIA for HOLAS 3. The thematic analyses is calculated for each assessment unit (1 km2 grid cells) and the data set covers the time period 2016-2021. Spatial Pressure and Impact Assessment (SPIA) is the framework for assessing spatial and cumulative pressures and impacts in HOLAS 3, and this analyses present a thematic assessment including only a certain subset of layers. The framework also includes results for the Baltic Sea Impact Index (full cumulative impact assessment), Baltic Sea Pressure Index (full cumulative pressure assessment), and other thematic assessments where a subset of pressure and ecosystem layers are used. For more info please - visit the HOLAS 3 website (http://stateofthebalticsea.helcom.fi/) - download the report thematic assessment of spatial distribution of pressures and impacts 2016-2021 (https://helcom.fi/post_type_publ/holas3_spa) - or check out the HELCOM SPIA online tool to make calculations for any desired combination of pressures and ecosystem layers (https://maps.helcom.fi/website/bsii/). Please scroll down to "Lineage" for a more detailed description of the methodology.

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.

The extraction of herring data set is based on: 1. http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/a3b67a55-7c1e-469e-b692-58c4e7b79279 for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle). Landing values were redistributed within each ICES rectangle by the c-square fishing effort data provided by ICES (all gears, 2011-2013). Tonnes / km² was calculated and the results were converted to 1 km x 1 km grid cells. The layer was log-transformed and normalised to produce the final pressure layer on extraction of herring. Please see "lineage" section below for further details on attributes, data source, data processing, etc.

The Baltic Sea Pressure index (BSPI) assesses the potential cumulative pressures in the Baltic Sea. The BSPI is based on georeferenced datasets of human activities (28 datasets), pressures (17 datasets) and uses the average sensitivity of each pressure layer to all ecosystems to weigth the pressure. Cumulative pressures are calculated for each assessment unit (1 km2 grid cells) and the data set covers the time period 2016-2021. Spatial Pressure and Impact Assessment (SPIA) is the framework for assessing spatial and cumulative pressures and impacts in HOLAS 3, where the BSPI presents the full cumulative pressure assessment where all pressures are included. The framework also includes results for the Baltic Sea Impact Index (full cumulative impact assessment), and other thematic assessments where a subset of pressure and ecosystem layers are used. For more info please - visit the HOLAS 3 website (http://stateofthebalticsea.helcom.fi/) - download the report thematic assessment of spatial distribution of pressures and impacts 2016-2021 (https://helcom.fi/post_type_publ/holas3_spa) - or check out the HELCOM SPIA online tool to make calculations for any desired combination of pressures and ecosystem layers (https://maps.helcom.fi/website/bsii/). Please scroll down to "Lineage" for a more detailed description of the methodology.

Essential fish habitat (EFH) map on Potential spawning areas for sprat was prepared in PanBalticScope project (co-founded by the European Maritime and Fisheries Fund of the European Union) http://www.panbalticscope.eu/ Sprat (Sprattus sprattus) occurs in the entire Baltic Sea, and mainly in open sea areas. It is assessed as a single stock in the Baltic Sea within fisheries management. Sprat eggs are pelagic, and sprat spawning is well known from the deep basins in the central Baltic, where it typically occurs from February to August. Further north, spawning starts later in the year, and is less certain. Recent fisheries surveys indicate that sprat spawning does no longer occur in the Gulf of Finland. Sprat spawning areas were delineated using environmental variables due to lack of coherent field data across the Baltic Sea countries. “Potential sprat spawning areas” were delineated as areas with salinity > 6 and water depth > 30 m, but for the Arcona basin depth > 20 m was used (Grauman, 1980, Bauman et al. 2006, Voss et al. 2012). “High probability spawning areas” were delineated for areas deeper than 70 m. Stock: Sprat in subdivisions 22-32 (ICES) EFH type: Potential spawning areas Approach: Environmental envelope, corrected for areas 20-40 m south of Bornholm. Variables and thresholds: Potential spawning area: Depth > 30 m, Salinity > 6 (annual average) High probability spawning area: Depth >70 m, Salinity > 6 (annual average) Quality: The map is based on literature and environmental variables, not actual data on sprat spawning. The map might overestimate the spawning area west and north of Gotland. The data layers on environmental variables are based on modelling. Attribute information: Raster value representing no spawning (0), potential spawning area (0.5) and high probability spawning area (1). References: - Baumann, H, H Hinrichsen, C Mollmann, F Koster, A Malzahn, and A Temming (2006) Recruitment variability in Baltic Sea sprat (Sprattus sprattus) in tightly coupled to temperature and transport patterns affecting the larval and early juvenile stages. Canadian Journal of Fisheries and Aquatic Science 63:2191-2201 - Grauman GB (1980) Long term changes in the abundance data of eggs and larvae of sprat in the Baltic Sea. Fisheries research in the Baltic Sea, Riga. 15:138-150 (in Russian) - HELCOM (2018) Outcome of the regional expert workshop on essential fish habitats, organized by Pan Baltic Scope project and HELCOM (HELCOM Pan Baltic Scope EFH WS 1-2018) - Voss R, MA Peck, HH Hinrichsen, C Clemmesen, H Baumann, D Stepputis, M Bernreuther, JO Schmidt, A Temming, and FW Köster (2012) Recruitment processes in Baltic sprat - A re-evaluation of GLOBEC Germany hypotheses. Progress in Oceanography 107:61-79

The map of herring relative abundance is mainly based on Baltic International acoustic surveys (BIAS), years 2011-2016 (ICES WGBIFS reports), reported as millions of herring / ICES rectangle. Also herring landings data were 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-2016). MAX-value = millions of herring according to BIAS in the ICES rectangle with highest landings. ICES rectangles outside the BIAS survey area with no reported herring landings were given the value 0. The relative abundance values in each ICES rectangle were divided by the area of the rectangle to obtain values per 1km2. If the values in small coastal ICES rectangles (outside BIAS area) became unrealistically large due to high herring landings, the value of the neighboring rectangle was given. The final layer was converted to 1 km x 1km grid cells. Values were first log transformed and normalized.

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.

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.