Large shallow inlets bays (according to Habitats Directive Annex I) are large, shallow indentations of the coast, sheltered from wave action and where, in contrast to estuaries, the influence of freshwater is generally limited. The distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on GIS analysis and modelling, but also field inventories and ground-truthing has been carried out in some areas. Data coverage, accuracy and the methods in obtaining the data vary between countries.

Estuaries (according to Habitats Directive Annex I) are coastal inlets that are strongly influenced by freshwater. The distribution map is based on data submission by HELCOM contracting parties. Most of the submitted data is based on modelling, GIS analysis and/or aerial photos. Data coverage, accuracy and the methods in obtaining the data vary between countries.

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

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 mixed substrate” includes classes “mixed sediment” 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.

Springtime Chl-a concentration is here used as a proxy for productive surface waters. In the Baltic Sea Impact Index (BSII), areas with high springtime phytoplankton production will be given higher importance, as they are considered important areas for the Baltic Sea food web. In the current map, mean of springtime maximum weekly values (weeks 12-22, years 2003-2011) Chl-a concentration of the surface waters has been used, derived from satellite data (MERIS). Years 2003-2011 have been used, as there is no MERIS data available for years 2012-2016. The data for eastern Baltic Sea is provided by the Finnish Environment Institute (~300m resolution). Outside this high resolution data, MERIS-data downloaded from JRC-database has been used (~4 km resolution, to calculate average of maximum monthly values for April or May for 2003-2011). Both datasets were converted to 1 km x 1 km grid cells.

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.

Introduction of radionuclides is based on HELCOM MORS Discharge data from 2011 to 2014. The isotopes taken into account were: Cesium-137, Strontium-90, and Cobalt-60. The decay-corrected annual average of the sum of the radionuclide discharges (in Bq) were calculated for the pressure. 10 km buffer with linear decreasing function was used to represent the impact distance from the nuclear power plant outlets.

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.

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.

This map presents the Special Protection Areas (SPAs) with reported wintering areas for birds. The spatial data on SPAs 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.