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.

Data shows the extent of land claim (permanent or temporary establishments of the sea) and the type of the construction. The data was made available by HELCOM Contracting Parties in response to data request. The data was received from Denmark, Finland, Sweden and Poland. The activity was declared as not relevant in Germany, Estonia, Latvia and Lithuanian. From Russia no data was reported. Attribute specification and units: Country: Country Type: Type of construction (land claim) Type_spec: More specified information about the type of land claim Year: Year of construction Estimated: Estimated year of construction from the identification information (environmental permit) given by the country in question Length: Length of the land reclamation (m) Area: Area (km2) of the land claim X_lon: Original Longitude coordinate point (for the data that has been transformed from points into lines) Y_Lat: Original latitude coordinate point (for the data that has been transformed from points into lines)

Distribution of blue mussel based on data submission by HELCOM contracting parties. Mainly pointwise occurrences of Mytilus spp. were submitted, originally gathered in national mapping and monitoring campaigns, or for scientific research. Point data from Poland was digitized based on Polish Marine Atlas. From Lithuania, a polygon delineating reefs was used to present Mytilus occurrence. For Germany, point data was complemented with a model describing Mytilus biomass in the German marine area (Darr et al. 2014), where predicted biomasses > 1g dw/ m2 were included as presence. 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 Mytilus) were generalized to 5km x 5km grid cells.

The dataset contains sand and gravel extraction activity during 2011–2015. The dataset is based on data submission by HELCOM Contacting Parties in response to data request. The dataset is quality assured and contains data from all the Baltic Sea countries. Common extraction sites with information on volume of extraction was included in the dataset. The other extraction sites (building sites, exclusive sites) had only summed data for all the sites, as volumes of extraction in specific sites was classified as confidential information. For these sites the volumes was calculated from reported sum of "other extraction areas". The amount was equally divided between corresponding extraction areas. For some common extraction sites a shared amount was reported.

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

Data set represents dredging activities around the Baltic Sea. The dataset contains information about the dredging activity and for some the type (capital, maintenance) and the year of activity as reported by HELCOM Contracting Parties in response to data request. The dredging data is missing from Denmark.

Baltic International Trawl Survey (BITS) data (2011-2016) from ICES DATRAS database was used as a base to create a map of cod relative abundance (quarter 1 data, CPUE values per ICES subdivision). Cod = 30cm was included. For ICES rectangles surveyed by BITS, values shown are the mean CPUE per ICES subdivision based on BITS data, average for 2011-2016. For ICES rectangles not surveyed by BITS, 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-2016). MAX-value = CPUE according to BITS in the ICES rectangle with highest landings. ICES rectangles outside the BITS survey area with no reported cod landings were given the value 0. Values were first log transformed and then 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.

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

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.