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.

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 “circalittoral mud” includes classes “Fine mud”, “Sandy mud” and “Mud to sandy mud” of the original data, in the circalittoral zone. The original polygon maps have been converted to 1 km x 1 km 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.

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.

This dataset is the first dedicated SMOS Sea Surface Salinity (SSS) product for the Baltic basin to enhance the science capabilities in the Baltic region and help to fill the gaps and grand challenges identified by the scientific community. These new product has been created under the funded ESA project ITT Baltic+ Salinity dynamics (4000126102/18/I-BG). This basin is one of the most challenging regions for the satellite SSS retrieval. The available EO-based SSS products are quite limited in terms of spatio-temporal coverage and quality. This is mainly due to technical limitations that strongly affect the brightness temperatures (TB), such as the high contamination by interferences and the contamination close to land and ice edges. Moreover, the sensitivity of TB to SSS changes is very low and dielectric models present limitations in this low salinity regime. Baltic+ L4 SSS product comprises 9 years (2011-2019) of daily maps at 0.05 degrees. A detailed explanation of the product algorithms and validation can be found at http://bec.icm.csic.es/doc/BEC_PD_SSS_Baltic_L3_L4.pdf and in the publication: Gonzalez-Gambau et al., “First SMOS Sea Surface Salinity dedicated products over the Baltic Sea“, Earth System Science Data, 2021 We present here the seasonal averaged Baltic+ L4 SSS products for the period 2011-2019. The daily Baltic+ L4 SSS products can be downloaded from the BEC FTP service (sftp://becftp.icm.csic.es) in the directory OCEAN/SSS/SMOS/Baltic/v1.0/L4/daily/

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 “circalittoral hard substrate” includes classes “Rock and other hard substrate” and “Coarse substrate” of the original data, in the circalittoral zone. The original polygon maps have been converted to 1 km x 1 km 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.

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.

Essential fish habitat (EFH) map on Potential spawning areas for sprat was prepared in PanBalticScope project (co-founded by the European Maritime and Fisheries Fund of the European Union) http://www.panbalticscope.eu/ Sprat (Sprattus sprattus) occurs in the entire Baltic Sea, and mainly in open sea areas. It is assessed as a single stock in the Baltic Sea within fisheries management. Sprat eggs are pelagic, and sprat spawning is well known from the deep basins in the central Baltic, where it typically occurs from February to August. Further north, spawning starts later in the year, and is less certain. Recent fisheries surveys indicate that sprat spawning does no longer occur in the Gulf of Finland. Sprat spawning areas were delineated using environmental variables due to lack of coherent field data across the Baltic Sea countries. “Potential sprat spawning areas” were delineated as areas with salinity > 6 and water depth > 30 m, but for the Arcona basin depth > 20 m was used (Grauman, 1980, Bauman et al. 2006, Voss et al. 2012). “High probability spawning areas” were delineated for areas deeper than 70 m. Stock: Sprat in subdivisions 22-32 (ICES) EFH type: Potential spawning areas Approach: Environmental envelope, corrected for areas 20-40 m south of Bornholm. Variables and thresholds: Potential spawning area: Depth > 30 m, Salinity > 6 (annual average) High probability spawning area: Depth >70 m, Salinity > 6 (annual average) Quality: The map is based on literature and environmental variables, not actual data on sprat spawning. The map might overestimate the spawning area west and north of Gotland. The data layers on environmental variables are based on modelling. Attribute information: Raster value representing no spawning (0), potential spawning area (0.5) and high probability spawning area (1). References: - Baumann, H, H Hinrichsen, C Mollmann, F Koster, A Malzahn, and A Temming (2006) Recruitment variability in Baltic Sea sprat (Sprattus sprattus) in tightly coupled to temperature and transport patterns affecting the larval and early juvenile stages. Canadian Journal of Fisheries and Aquatic Science 63:2191-2201 - Grauman GB (1980) Long term changes in the abundance data of eggs and larvae of sprat in the Baltic Sea. Fisheries research in the Baltic Sea, Riga. 15:138-150 (in Russian) - HELCOM (2018) Outcome of the regional expert workshop on essential fish habitats, organized by Pan Baltic Scope project and HELCOM (HELCOM Pan Baltic Scope EFH WS 1-2018) - Voss R, MA Peck, HH Hinrichsen, C Clemmesen, H Baumann, D Stepputis, M Bernreuther, JO Schmidt, A Temming, and FW Köster (2012) Recruitment processes in Baltic sprat - A re-evaluation of GLOBEC Germany hypotheses. Progress in Oceanography 107:61-79

The occurrence of suitable nursery habitats is crucial for maintaining fish populations (Sundblad et al. 2013). Species distribution modelling studies have shown the importance of suitable environmental conditions for pikeperch recruitment. Due to lack of coherent data on pikeperch spawning and nursery areas across the Baltic Sea countries, environmental variables were used in delineating potential recruitment areas for pikeperch. The pikeperch recruitment area presented on the map is mainly delineated by selecting areas where depth < 5 m, logged exposure < 5, salinity < 7 PSU, Secchi depth < 2 m and distance to deep (10m) water < 4km. The threshold values have been obtained from literature (Veneranta et al. 2011, Bergström et al. 2013, Sundblad et al. 2013, Kallasvuo et al. 2016). Temperature, although important for pikeperch, was left out due to high variation in timing of suitable spawning temperatures across the Baltic Sea. In Finnish coastal waters, a national pikeperch model (Kallasvuo et al. 2016) has been used, with very suitable areas for pikeperch generalized to 1 km grid. In Sweden, the areas delineated by environmental variables have been complemented with information from national interview survey (Gunnartz et al. 2011) as well as expert opinion.

The Baltic Sea Pressure index (BSPI) assesses the potential cumulative pressures in the Baltic Sea. The BSPI is based on georeferenced datasets of human activities (28 datasets), pressures (17 datasets) and uses the average sensitivity of each pressure layer to all ecosystems to weigth the pressure. Cumulative pressures are 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, where the BSPI presents the full cumulative pressure assessment where all pressures are included. The framework also includes results for the Baltic Sea Impact Index (full cumulative impact 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.

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.