Concentration of nitrogen pressure layer is interpolated from annual seasonal average of total nitrogen concentrations 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).

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.

Pressure layer combines all human activities that cause physical disturbance or damage to seabed. For several human activity datasets, spatial extents were given (table below). Buffers with decreasing value rates were applied to represent the impact distance of physical disturbance. The following human activities were combined into the physical disturbance layer; - Cables (under construction, 1 km buffer) - Coastal defence and flood protection (under construction, 500 m buffer) - Deposit of dredged material (500 m buffer for points and areas) - Dredging (maintenance) (500 m buffer for points and areas) - Extraction of sand and gravel (500 m buffer) - Finfish mariculture (1 km buffer) - Fishing intensity 2011-2016 average (subsurface swept area ratio) - Furcellaria harvesting - Pipelines (0,3 km buffer) - Recreational boating and sports - Shellfish mariculture - Shipping density - Wind farms (under construction) (1 km buffer) - Wind farms (operational) (0,1 km buffer) The human activity data sets were first processed separately covering the whole Baltic Sea and then summed together. In this integration, some data layers were down-weighted to arrive at a balanced pressure layer, as described below. High pressure intensity and/or slow recovery (weighting factor 1): Coastal defence and flood protection, Deposit of dredged material, Dredging, Extraction of sand and gravel and Fishing intensity Moderate to high (Weighting factor 0,8): Pipelines and Shipping density Moderate (Weighting factor 0,6): Finfish mariculture, Shellfish mariculture and Wind farms (under construction) Low to moderate (Weighting factor 0,4): Cables Low (Weighting factor 0,2): Maerl and Furcellaria harvesting, Recreational boating and sports and Wind farms (operational) Harbours and marinas were left out from the physical disturbance pressure to avoid double counting due to their representation in the shipping density and recreational boating and sports data sets.

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.

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

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.

Summary This dataset shows model results for the average bottom temperature in the Baltic region in the plant growth season from April to September. Description This dataset shows model results for the average bottom temperature in the Baltic region in the plant growth season from April to September.

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.

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.

Raster grid of the Baltic Sea bathymetry computed with ArcGIS Spatial Analyst (KRIGING) from the original Digital Topography of the Baltic Sea (IOWTOPO) database produced by the Baltic Sea Research Institute of Warnemunde. Output resolution of the grid is 250 m, data is projected into ERTS89_LAEA CRS (Lambert Azimuthal Equal Area projection, ETRS89 datum), file format is Erdas Imagine (IMG), data format is continuous, float.