Input of impulsive anthropogenic sound includes impulsive events from 2011-2016 • Seismic surveys (HELCOM-OSPAR Registry; national data call submissions as lines in the folder of data) • Explosions (HELCOM-OSPAR Registry) • Pile driving (HELCOM-OSPAR Registry) • Airguns (HELCOM-OSPAR Registry) For the different event types, numeric intensity value was used to represent the pressure as categorized in HELCOM-OSPAR Impulsive noise registry. All nationally reported seismic surveys were given intensity values “Very low” (0.25) - Very low (0.25) - Low (0.5) - Medium (0.75) - High (1) The impact distance has not been taken into account due to the different nature of separate datasets used for the pressure layer. We acknowledge that e.g. pile driving and airguns may impact up to 20 km from the source event. The spread of the sound wave depends on the sound frequency, water salinity, temperature and density.

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.

This map presents the Special Protection Areas (SPAs) with reported breeding areas for birds. The spatial data on Special Protection Areas were gathered from the HELCOM contracting parties by Lund University, Sweden. In the data, the countries also indicated whether the sites were designated mainly due to wintering or breeding birds in the area. For Denmark, the information was obtained from standard forms for Natura 2000 sites. For Denmark, the data was updated after review process 20 February 2017. For Germany, the areas that were reported as “NA”(=information not available) were included in both breeding and wintering area maps. Many of the SPAs are both wintering and breeding areas. For the Baltic Sea Impact Index, the data was converted to 1 km x 1km grid cells.

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.

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.

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

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 “infralittoral sand” includes classes “Sand” and “Muddy sand” 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.

Potential effect of continuous noise to mobile species assesses the cumulative potential effect of continuous noise on mobile species and their presence in the HELCOM area. The evaluation is based on the pressure layer on input of continuous noise, combined with information on the distribution of 15 mobile species and their habitats. The thematic analyses on the potential effect of continuous noise to mobile species 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.