Distribution of Charophytes (Chara spp., Nitella spp., Nitellopsis spp., Tolypella spp.) mainly based on data submission by HELCOM contracting parties. Submitted point data was originally gathered in national mapping and monitoring campaigns, or for scientific research. Also scientific publications were used to complement the data (in Curonian, Vistula and Szczechin lagoons, see reference list). Polygon data from Poland was digitized based on Polish Marine Atlas. From Estonian waters, a predictive model was used (200m resolution), that was converted to presence/absence using minimized difference threshold (MDT) criteria. All data (points, polygon and the raster presenting predicted presence of Charophytes) were generalized to 5km x 5km grid cells.

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.

Input of heat pressure dataset contains delta heat values from warm water discharge of - Discharge of warm water from nuclear power plants - Fossil fuel energy production. The Discharge of warm water from nuclear power plants was requested from HELCOM Contracting Parties. Average temperature change of the cooling water (°C) and amount of cooling water (m3) was used to calculate the heat load in TWh. No data on heat load was available for the Leningrad nuclear power plant; therefore the average heat load (TWh) of discharge of warm water from nuclear power plants was given. No heat load information was available for fossil fuel energy production sites. Heat load of 1 (TWh) was given to all production sites, based on average heat load of individual production site in Helsinki from recent years. 1 km buffer was used for both datasets with steep decline around the outlet. Overlapping areas were summed. Heat load from both layers were summed and the layer was normalized.

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 mixed substrate” includes classes “mixed sediment” 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 mud” includes classes “Fine mud”, “Mud to sandy mud” and “Sandy mud” 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.

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.

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.

This dataset contains modelled small vessel fuel consumption. This describes the geographical distribution of the fuel used by small boats. The total fuel consumption was modelled in SHEBA project to study emissions from pleasure boats. The model is based on locations and berths in marinas and leisure harbours, AIS information, statistics on fuel sale and extensive survey. For 2018 version the layer is weighted with depth, log-transformed and normalised (please see below). This dataset was also used on HOLAS 3.

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.

The occurrence of suitable nursery habitats is crucial for maintaining fish populations (Sundblad et al. 2013). For perch, species distribution modelling studies (Snickars et al. 2010, Bergström et al. 2013, Sundblad et al. 2013) have shown the importance of suitable environmental conditions for reproduction. Due to lack of coherent data on perch spawning and nursery areas across the Baltic Sea countries, environmental variables were used in delineating potential recruitment areas for perch. The distribution area or perch recruitment is delineated by selecting areas where depth < 4 m (For Danish waters < 3 m), logged exposure < 5 (exposure model described in Isæus 2004), and salinity < 10 PSU. The threshold values have been obtained from literature (Snickars et al. 2010, Bergström et al. 2013, Skovrind et al. 2013, Sundblad et al. 2013). Relatively “loose” thresholds have been used, to rather overestimate than underestimate the recruitment area (precautionary approach). Along the Finnish coastline a national model has been used (Kallasvuo et al. 2016), with suitable environments for perch recruitment generalized to 1 km x 1 km grid.