The pressure layer represents biological pressure caused by introduction of non-indigenous species. The data is obtained from core indicator Trend in the arrival of new non-indigenous species (BSEP 129b: http://www.helcom.fi/Lists/Publications/BSEP129B.pdf). For the Baltic Sea Impact Index, 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 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.

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.

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.

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

Potential cumulative impacts on benthic habitats is based on the same method than http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/9477be37-94a9-4201-824a-f079bc27d097, but is focused on physical pressures and benthic habitats. The dataset was created based on separate analysis for potential cumulative impacts on only the benthic habitats, as these are particularly affected by physical pressures. In this case the evaluation was based on pressure layers representing http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/ea0ef0fa-0517-40a9-866a-ce22b8948c88 and http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/05e325f3-bc30-44a0-8f0b-995464011c82, combined with information on the distribution of eight broad benthic habitat types and five habitat-forming species (http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/363cb353-46da-43f4-9906-7324738fe2c3, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/f9cc7b2c-4080-4b19-8c38-cac87955cb91, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/264ed572-403c-43bd-9707-345de8b9503c, http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/822ddece-d96a-4036-9ad8-c4b599776eca and http://metadata.helcom.fi/geonetwork/srv/eng/catalog.search#/metadata/ca327bb1-d3cb-46c2-8316-f5f62f889090). The potential cumulative impacts has been estimated based on currently best available data, but spatial and temporal gaps may occur in underlying datasets. Please scroll down to "Lineage" and visit http://stateofthebalticsea.helcom.fi/cumulative-impacts/ for more info.

The seals' distribution maps show the distribution and abundance of grey, harbour and ringed seals across the Baltic Sea. The ecosystem component maps on mammals' distribution were drafted by EG MAMA harbour porpoise and seal distribution teams. The maps were prepared as expert-derived distribution categories to be used in the HELCOM Third Holistic Assessment of the Ecosystem health of the Baltic Sea.

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.

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.

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.