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.

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.

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.

The dataset contains total landings of cod for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle) under EU Joint Research Centre’s data collection framework for fisheries data. Russian data extracted from ICES annual reports.

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

'Availability of deep water habitat, based on occurrence of H2S' layer describes the suitability of the bottom areas for the Baltic Sea biota, with regard to oxygen conditions of the near bottom waters. The data used to produce the layer was received from Leibniz-Institut für Ostseeforschung Warnemünde (IOW): - areas (polygons) with hydrogen sulfide (H2S) based on point measurements and modelling. Five time periods / year, for years 2011-2016 (altogether 30 layers). The polygons were converted to raster layers in a way, that for each time period (6 years, 5 time periods each year), areas with H2S got a value 0, other areas got the value 1. All layers were summed, (representing 6 years, 5 time periods each year, maximum value 30) and data was normalised. For more detailed information on the data used, please see Feistel et al. 2016.

The dataset contains total landings of sprat for years 2011-2016 reported per ICES statistical rectangles (tonnes / ICES rectangle) under EU Joint Research Centre’s data collection framework for fisheries data. Russian data extracted from ICES annual reports.

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.