# Indicator-based assessment of marine biological diversity – lessons from 10 case studies across the European Seas

Laura Uusitalo, Hugues Blanchet, Jesper H. Andersen, Olivier Beauchard, Torsten Berg, Silvia Bianchelli, Annalucia Cantafaro, Jacob Carstensen, Laura Carugati, Sabine Cochrane, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Laura Uusitalo, Hugues Blanchet, Jesper H. Andersen, Olivier Beauchard, Torsten Berg, et al. "Indicator-based assessment of marine biological diversity – lessons from 10 case studies across the European Seas. Frontiers in Marine Science, 2016, 3, pp.159. 10.3389/fmars.2016.00159. hal-04165566

HAL Id: hal-04165566

https://hal.science/hal-04165566

Submitted on 19 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

CrossMark

# Indicator-Based Assessment of Marine Biological Diversity—Lessons from 10 Case Studies across the European Seas

Laura Uusitalo¹*, Hugues Blanchet²,³, Jesper H. Andersen⁴, Olivier Beauchard⁵,
Torsten Berg⁶, Silvia Bianchelli⁷, Annalucia Cantafaro⁷, Jacob Carstensen⁸,
Laura Carugati⁷, Sabine Cochrane²,⁹, Roberto Danovaro⁷,¹⁰, Anna-Stiina Heiskanen¹,
Ville Karvinen¹, Snejana Moncheva¹¹, Ciaran Murray⁴,⁸, João M. Neto¹², Henrik Nygård¹,
Maria Pantazi¹³, Nadia Papadopoulou¹⁴, Nomiki Simboura¹⁵, Greta Srébalienė¹⁶,
Maria C. Uyarra¹⁷ and Angel Borja¹⁷

OPEN ACCESS

Edited by:

Michael Elliott,
University of Hull, UK

Reviewed by:

Jan Marcin Weslawski,
Institute of Oceanology of the Polish
Academy of Sciences, Poland
Bernardo Antonio Perez Da Gama,
Federal Fluminense University, Brazil

*Correspondence:

Laura Uusitalo
laura.uusitalo@ymparisto.fi;
laura.uusitalo@iki.fi

Specialty section:

This article was submitted to
Marine Ecosystem Ecology,
a section of the journal
Frontiers in Marine Science

Received: 08 June 2016

Accepted: 19 August 2016

Published: 06 September 2016

Citation:

Uusitalo L, Blanchet H, Andersen JH,
Beauchard O, Berg T, Bianchelli S,
Cantafaro A, Carstensen J,
Carugati L, Cochrane S, Danovaro R,
Heiskanen A-S, Karvinen V,
Moncheva S, Murray C, Neto JM,
Nygård H, Pantazi M,
Papadopoulou N, Simboura N,
Srèbalienè G, Uyarra MC and Borja A
(2016) Indicator-Based Assessment of
Marine Biological Diversity—Lessons
from 10 Case Studies across the
European Seas.
Front. Mar. Sci. 3:159.
doi: 10.3389/fmars.2016.00159

¹ Marine Research Centre, Finnish Environment Institute (SYKE), Helsinki, Finland, ² SALT Lofoten AS, Svolvaer, Norway, ³ UMR 5805 EPOC, University of Bordeaux, Talence, Bordeaux, France, ⁴ NIVA Denmark Water Research, Copenhagen, ⁵ Estuarine and Delta Systems Department, Royal Netherlands Institute for Sea Research, Yerseke, Netherlands, ⁶ MariLim Aquatic Research GmbH, Schönkirchen, Germany, ⁷ Department of Life and Environmental Sciences, Marche Polytechnic University, Ancona, Italy, ⁸ Department of Bioscience, Aarhus University, Roskilde, Denmark, ⁹ Akvaplan-niva, FRAM - High North Centre for Climate and the Environment, Tromsø, Norway, ¹⁰ Stazione Zoologica Anton Dohm Napoli, Napoli, Italy, ¹¹ Marine Biology and Ecology, Institute Of Oceanology - Bulgarian Academy of Sciences, Varna, Bulgaria, ¹² Department of Life Sciences, IMAR - Institute of Marine Research, University of Coimbra, Coimbra, Portugal, ¹³ Hellenic Centre for Marine Research, Institute of Marine Biological Resources and Inland Waters, Athens, Greece, ¹⁴ Hellenic Centre for Marine Research, Institute of Marine Biological Resources and Inland Waters, Heraklion, Greece, ¹⁵ Institute of Oceanography, Hellenic Centre for Marine Research, Attica, Greece, ¹⁶ Marine Science and Technology Centre, Klaipėda University, Klaipėda, Lithuania, ¹⁷ Marine Research Division, AZTI Tecnalia, Gipuzkoa, Spain

The Marine Strategy Framework Directive requires the environmental status of European marine waters to be assessed using biodiversity as 1 out of 11 descriptors, but the complexity of marine biodiversity and its large span across latitudinal and salinity gradients have been a challenge to the scientific community aiming to produce approaches for integrating information from a broad range of indicators. The Nested Environmental status Assessment Tool (NEAT), developed for the integrated assessment of the status of marine waters, was applied to 10 marine ecosystems to test its applicability and compare biodiversity assessments across the four European regional seas. We evaluate the assessment results as well as the assessment designs of the 10 cases, and how the assessment design, particularly the choices made regarding the area and indicator selection, affected the results. The results show that only 2 out of the 10 case study areas show more than 50% probability of being in good status in respect of biodiversity. No strong pattern among the ecosystem components across the case study areas could be detected, but marine mammals, birds, and benthic vegetation indicators tended to indicate poor status while zooplankton indicators indicated good status when included into the assessment. The analysis shows that the assessment design, including the selection of indicators, their target values, geographical resolution
and habitats to be assessed, has potentially a high impact on the result, and the assessment structure needs to be understood in order to make an informed assessment. Moreover, recommendations are provided for the best practice of using NEAT for marine status assessments.

**Keywords: biodiversity, assessment tool, MSFD, environmental status, spatial aggregation, integration, indicator sensitivity**

## INTRODUCTION

Biological diversity is widely recognized as one of the cornerstones of healthy ecosystems (e.g., *Worm et al., 2006*). Diversity may safeguard ecosystems against undesired regime shifts (*Folke et al., 2004*) and guarantee the continued delivery of ecosystem goods and services (*Duarte, 2000; Beaumont et al., 2007*). The need to maintain biodiversity is also recognized by international legislation (e.g., *Convention of Biological Diversity; UNEP, 1992*); to European Union (EU) level, the Marine Strategy Framework Directive (MSFD; *European Union, 2008*) requires its member states to assess the status of marine biodiversity and take action to guarantee that it remains at, or is restored to, Good Environmental Status (GES). A definition of what can be interpreted as good status can be consulted in *Borja et al. (2013)*.

In order to conduct an assessment of status, and to determine the effectiveness of any implemented remedial measures, we need a clear definition of biodiversity and a unified approach for its assessment. In the marine assessments like MSFD, biodiversity is defined on the level of species, communities, habitats, and ecosystems, as well as in the genetic level (*Cochrane et al., 2010*). Indicators that show the ecosystem response to human pressures form the basis of the tool kit with which we can describe environmental status (*Borja et al., 2016*). Based on qualitative environmental objectives, targets are set for each indicator which allow policy makers to implement management measures should these not be reached (*Borja et al., 2012*).

One of the challenges faced during the first round of MSFD initial assessments is the diverging data availability for biodiversity across highly variable systems, but yet an overarching need to conduct compatible assessments across European regional seas (*Hummel et al., 2015*). European marine ecosystems comprise a complexity and variability both in space and time, ranging from fully saline systems such as in Mediterranean and Atlantic waters to the brackish Baltic Sea, and exposed open water systems such as in the northern Norwegian and Barents seas to fully enclosed systems such as the Black Sea. The levels of available knowledge and data within these systems vary, as well as the biological parameters and indicators used for assessments (*Hummel et al., 2015*).

The conclusions of the European Commission, in their evaluation of the EU member states' reports on the initial assessment carried out in 2010–2012 was that there is an apparent lack of coherence and comparability in the indicators used and in the final evaluation of the overall status, between the countries and within all regional seas (*Palialexis et al., 2014*). Therefore, there is an urgent need for coherent frameworks and methodologies to allow consistent approach in biodiversity

status assessment across the European Regional Seas. This would also be needed in order to allow coherence in the biodiversity assessments for the EU Birds and Habitats directives and the EU Biodiversity Strategy 2020.

While we could argue that we cannot compare studies if we do not have directly comparable datasets, in practice this is rarely possible, and certainly not at large spatial scales, or involving multiple research institutes and member states. Since there is no single way of describing biodiversity that fits all purposes, and since regional seas have intrinsic differences, we need a pragmatic selection of indicators which are appropriate to the specific questions asked, as well as a flexible and transparent indicator-based tool for assessment of biodiversity status. There is a large number of operational indicators, which have been used to describe the status in different types of aquatic systems (*Birk et al., 2012; Borja et al., 2016*). As biological diversity is multifaceted, including different taxonomic and functional groups, it cannot be expressed with a single indicator. Consequently, sets of different indicators are needed to cover the broad aspects of biological diversity and it is their combination into a single assessment that becomes a challenge (*Borja et al., 2014; Probst and Lynam, 2016*). In order to obtain a single overall assessment value, or conclusion, the results of the multiple indicators used in the assessment need to be aggregated, depending on the purpose of the assessment; e.g., if the aim is to inform different stakeholders and to set overall targets for the improvement of the marine environment, or depending on the assessment scale (*Borja et al., 2014*). Clear and transparent aggregation and integration rules are needed to interpret indicator information onto an environmental status assessment (see *Borja et al., 2014* for a review on integration methods).

A variety of assessment tools enabling the integration of indicators already exists (see e.g., *HELCOM, 2009a; Andersen et al., 2014; Borja et al., 2016*). However, only few of them have treated biological diversity in a comprehensive way, have been tested broadly (i.e., outside the region in which they have been developed), or consider the complexity at an adequate level of detail for the spatial scale for which they are applied. To overcome these issues, in the context of the EU funded project DEVOTES (DEVelopment Of innovative Tools for understanding marine biodiversity and assessing GES), the Nested Environmental status Assessment Tool (NEAT; *Berg et al., 2016; Borja et al., 2016*) has been developed to assess biodiversity status of marine waters under the MSFD. NEAT uses a combination of high-level integration of habitats and spatial units, and averaging approach (*Borja et al., 2014*), allowing for specification on structural and spatial levels, applicable to any geographical scale.
In this contribution NEAT is applied to the assessment of marine biological diversity in 10 different case studies distributed across the European regional seas (Figure 1). The assessment results are discussed, but the main focus of the paper is on: (i) analyzing the outcome of these assessments in light of the practical choices that have to be made to apply this tool, and (ii) proposing best practices for marine biological diversity assessment using this tool.

# MATERIALS AND METHODS

## Case Study Areas

The case study areas were selected to represent a wide range of marine systems (Figure 1), with different climatic and hydrographic characteristics as well as exposure to different

human activities and management challenges (Table 1). These areas represent a wide range of marine biogeographical areas from subtropical waters to temperate and Arctic, covering the four European regional seas (i.e., Mediterranean, Atlantic, Black, and Baltic Seas). The surface areas of these case studies varied from <3000 km² in Saronikos Gulf (Greece) to >820,000 km² in the Barents Sea (Norway; Table 1). Detailed descriptions of the case study areas, with relevant references, can be found in Supplementary Material (S1-S10).

## NEAT

NEAT is a structured, hierarchical tool for making marine status assessments (Berg et al., 2016; Borja et al., 2016), and freely available at www.devotes-project.eu/neat. In NEAT, the study area can be divided into hierarchical spatial assessment units

<figure><img src="image_4.png" /><figcaption>FIGURE 1 | The case study areas. For the area codes, see Table 1. More detailed case study maps can be found in Supplementary Material (S1-S10).</figcaption></figure>
<table><caption>TABLE 1 | Characteristics of the case study areas.</caption><thead><tr><th>Case study name</th><th>Code</th><th>Regional Sea, subregion</th><th>Salinity range</th><th>Seas Surface temperature range (°C)</th><th>Mean/max depth (m)</th><th>Special features</th><th>Main pressures</th><th>Biodiversity components in poorest status</th><th>Biodiversity components in best status</th><th>Observed gaps in biodiversity components</th><th>Observed gaps in habitats</th><th>Key references</th></tr></thead><tbody><tr><td>Norwegian Barents Sea–Lofoten</td><td>BSL</td><td>North-East Atlantic, Barents Sea</td><td>34.4–35</td><td>-1.5–7.5, permanent/ seasonal ice cover/ice-free areas</td><td>230/500</td><td>Very large ecosystem; polar region</td><td>Climate change, fishing</td><td>Benthic vegetation, birds</td><td>Phytoplankton, Benthic fauna, zooplankton, fish</td><td>Lack of indicators for microbes, reptiles, large benthic organisms, cephalopods, alien species</td><td>Rocky and mixed substratum and biogenic reef (ex: Lophelia reefs)</td><td>(Wassmann et al., 2006)</td></tr><tr><td>Gulf of Finland</td><td>FIN</td><td>Baltic Sea</td><td>0–7</td><td>0–20, seasonal ice cover</td><td>37/100</td><td>Estuary; Heavy shipping traffic incl. oil tankers; Large drainage basin; Important resting and breeding area for migratory birds</td><td>Nutrient inputs, hazardous substances</td><td>Birds, water column habitat, benthic fauna, mammals</td><td>Zooplankton</td><td>Hard bottom fauna, microbes, alien species; habitat-level indicators</td><td>Hard bottoms</td><td>Ojaveer et al., 2010</td></tr><tr><td>Lithuanian marine waters</td><td>LIT</td><td>Baltic Sea</td><td>5–18</td><td>9.6–10.4</td><td>50/125</td><td>High wave exposure</td><td>Tourism, recreation, ports, shipping, urban, industry, agriculture, fisheries</td><td>Phytoplankton, birds, benthic vegetation</td><td>Benthic habitat, benthic fauna, benthic vegetation, zooplankton</td><td>Mammals</td><td></td><td></td></tr><tr><td>Kattegat</td><td>KAT</td><td>Baltic Sea, North Sea</td><td>18–32 (surface)</td><td>2–19 (surface) 5–14 (bottom)</td><td>23.9/80</td><td>Transition zone between the Baltic Sea and the North Sea.</td><td>Nutrient inputs, hazardous substances, fishing, physical modification</td><td>Birds, phytoplankton</td><td>Zooplankton, fish</td><td>Non-commercial fish, alien species, microbes</td><td>Hard bottoms</td><td>Naturstyrelsens, 2012; Andersen et al., 2013</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>Case study name</th><th>Code</th><th>Regional Sea, subregion</th><th>Salinity range</th><th>Seas Surface temperature range (°C)</th><th>Mean/max depth (m)</th><th>Special features</th><th>Main pressures</th><th>Biodiversity components in poorest status</th><th>Biodiversity components in best status</th><th>Observed gaps in biodiversity components</th><th>Observed gaps in habitats</th><th>Key references</th></tr></thead><tbody><tr><td>Dutch North Sea</td><td>DNS</td><td>North Sea</td><td>27.6-34.8</td><td>4-18</td><td>6/54</td><td>Shallow and highly productive area, soft bottom sediments, contrasted benthic habitats (stressed in the southern part, rarely disturbed in the northern one)</td><td>Beam-trawl fishing, organic loading, harbor infrastructures</td><td>Fish</td><td>Birds, mammals, benthic fauna</td><td>Benthic vegetation</td><td>Herman et al., 2014</td><td></td></tr><tr><td>Basque Coast</td><td>BAS</td><td>North-East Atlantic</td><td>32.5-35.5</td><td>11-23</td><td>Range 0-6000</td><td>High wave exposure, long fetch</td><td>Fishing/ shipping</td><td>Mammals</td><td>Phytoplankton, zooplankton, benthic habitat, water column habitat</td><td>Microbes</td><td>Deep and abyssal habitats</td><td>Borja and Collins, 2004</td></tr><tr><td>Portuguese continental subdivision</td><td>POR</td><td>North-East Atlantic</td><td>35-36</td><td>13-23</td><td>Range 0-6000</td><td>High wave exposure and seasonal upwelling at north-west coast, smoother conditions at south and south-west coast</td><td>Coastal trawling (also dredging, sediment disposal, harbors)</td><td>Fish, birds, pelagic fauna</td><td>Zooplankton</td><td>Reptiles, mammals, microbes. Species distribution and coverage area. The current status of all stocks has been assessed, although using secondary indicators</td><td>Deep and abyssal habitats and open-sea areas (>20 nm)</td><td>MAMAOT, 2012 (https://dl.dropbox usercontent.com/u/ 103729442/Estrateg iaMarinha_subdv_Co ntinente.pdf) Dupont et al., 2014</td></tr><tr><td>Black Sea coast (Large Varna Bay)</td><td>BSC</td><td>Black Sea</td><td>11-18</td><td>0-28</td><td>30</td><td>Hydrologically dynamic area</td><td>Tourism industry, urbanization, port activity, aquaculture</td><td>Fish, phytoplankton, zooplankton</td><td>Water column habitat, benthic habitat</td><td>Microbes, alien species, food-web indicators; insufficient fish monitoring</td><td>No monitoring in the open sea areas</td><td>National Initial Assessment GES Report, 2013</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>Case study name</th><th>Code</th><th>Regional Sea, subregion</th><th>Salinity range</th><th>Seas Surface temperature range (°C)</th><th>Mean/max depth (m)</th><th>Special features</th><th>Main pressures</th><th>Biodiversity components in poorest status</th><th>Biodiversity components in best status</th><th>Observed gaps in biodiversity components</th><th>Observed gaps in habitats</th><th>Key references</th></tr></thead><tbody><tr><td>Saronikos Gulf</td><td>SAR</td><td>Mediterranean (Eastern Mediterranean/ Aegean Sea)</td><td>38–39</td><td>12–28</td><td>100/450</td><td>It is the natural marine gateway of the city of Athens and Piraeus harbor. It receives the effluents of the central sewage outfall of Athens through a deep underwater outlet.</td><td>Activities: Waste Water Treatment Plant (WWTP), Shipping, tourism, fishing.<br>Pressures: input of organics/nutrients, contaminants, habitat loss, resource exploitation, alien species, litter</td><td>Mammals, fish</td><td>Reptiles</td><td>Birds, zooplankton</td><td></td><td>Simboura et al., 2014, 2015; SoHeIFME, 2005; SoHeIFI, 2007</td></tr><tr><td>Adriatic Sea</td><td>ADR</td><td>Mediterranean</td><td>37–39</td><td>18–25</td><td>35/1200</td><td>Cyclonic circulation; input from the Po river</td><td>Organic loading; phytoplankton overexploitation; bottom trawling</td><td>Mammals, phytoplankton</td><td>Fish, birds, zooplankton</td><td>Reef and mud habitats</td><td></td><td>Artegiani et al., 1997a; Artegiani et al., 1997b; (UNEP(DEPI)/MED WG.408/Inf.14, 2015).</td></tr></tbody></table>
(SAU) and habitat types (HBT); e.g., SAU “archipelago zone” could include “inner archipelago” and “outer archipelago” as lower-level SAUs, and they, in turn, could include, e.g., water bodies as yet lower-level SAUs. Similarly, the HBT “seafloor” could include HBTs “soft bottom” and “hard bottom,” which again could be further sub-divided (Figure 2). NEAT classifies the status of each SAU based on indicators that have been defined for that SAU; if one SAU has indicators describing different HBTs, the status of each HBT within a SAU is assessed first, and each HBT is then given equal weight in assessing the status of the SAU. The overall assessment is an average of the SAUs, weighted by their surface areas (km²). Other weighting schemes can be applied, if desired.

Each indicator must be explicitly linked to a SAU and a HBT— the same indicator, e.g., “the maximum depth of seaweed,” can be included multiple times for multiple SAUs and HBTs if it has been assessed for multiple areas. These instances of indicators are called “indicator values” in this paper, while the indicators describing a certain ecological concept, e.g., the growth depth of a macrophyte species, or the reproduction rate of a bird species, are called “unique indicators.”

In order to aggregate indicators by weighted average, it is necessary to transform all indicators to a common scale. In NEAT, indicators are transformed into values that range from 0 to 1 using a continuous piecewise linear function. On this scale, the value of 0.6 corresponds to the boundary between

good (>0.6) and not good (<0.6) status. Transformation to this scale is defined by specifying the values of the indicator in the original measurement scale, which corresponds to the transformed values of 0, 0.2, 0.4, 0.6, 0.8, and 1.0. Though the transformation function is piecewise linear, the definition of 5 segments allows a reasonable approximation to non-linear functions. These five segments are also used here for illustrative purposes, and they are called bad/poor/moderate/good/high classes, although it is recognized that the boundary between GES and non-GES lies between the “moderate” and “good” classes.

## Indicator Selection and Specification

The indicators used for this assessment represent the best available data and expertise for the six biological descriptors of the MSFD [i.e., D1 (biodiversity), D2 (non-indigenous species), D3 (commercially important species), D4 (food webs), D5 (eutrophication), and D6 (sea floor integrity)] in each case study area. These indicators include the national and regional indicators used for the MSFD assessment, and indicators derived from scientific literature and expertise. They have been selected to be representative of various biodiversity components, habitats, and geographical areas relevant for each case study area; however it is possible that no indicators exist to be used for some relevant components. The list of indicators included in each case study is available in Supplementary Material S11.

<figure><img src="image_5.png" /></figure>
Each indicator is associated to an ecosystem component class that describes the ecosystem component that the indicator describes. In this study, 12 ecosystem components were defined in order to accommodate all indicators used in all of the case studies. These components were phytoplankton, zooplankton, fish, reptiles, marine mammals, birds, benthic fauna, benthic vegetation, pelagic fauna (composite indicators consisting of data from multiple pelagic fauna groups), all taxa (composite indicators consisting of data from multiple taxa), benthic habitat, and water column habitat. The latter two components gathered indicators related to physico-chemical conditions of the habitat, necessary to maintain life (e.g., oxygen or nutrients), whilst the "all taxa," benthic fauna, and pelagic fauna groups included composite indicators encompassing many species groups; the other nine ecosystem components were taxonomic groups.

## Biodiversity Status

The status of the biological diversity was assessed for each case study area using NEAT. The analysis provides an overall assessment for each case study area and a separate assessment for each of the ecosystem components included in the assessment. The final value has an associated uncertainty value, which is the probability of being in a determinate class status (GES/non-GES). This uncertainty was determined by the standard error linked to the indicator values (Carstensen and Lindegarth, 2016).

## Evaluation of Assessment Design and Its Effects on the Status Assessment

The application of NEAT to a broad range of marine regions provides an opportunity to test and compare the NEAT assessment approaches and evaluate the consequences of design choice for the general environmental status assessment. How the available data are combined within the tool might have consequences on the results of the status assessment of biodiversity (Borja et al., 2014; Probst and Lynam, 2016). Therefore, one of our aims is to evaluate the consequences of the way the assessment was designed on the general assessment result.

NEAT gives a framework to organize the assessment, but it does not prescribe the number of assessment components, i.e., indicators, SAUs, HBTs, or ecosystem components to be used in an assessment. The user has the option to organize the different components of NEAT depending on the case, e.g., the morphological characteristics of the area, availability and resolution of data, and how the selected local indicators are defined.

In order to describe the assessment design, the following key components were summarized for each case study: (i) the total number of SAUs and how many hierarchical SAU levels there are, (ii) the total number of HBTs and their hierarchical levels, (iii) the number of ecosystem components covered by the indicators, (iv) the number of unique indicators (i.e., not repetition of the same indicator on a different spatial unit), as well as (v) the quantity of data, defined as the number of different indicator values (e.g., if the same indicator is defined separately for five different SAUs, they would comprise five indicator values).

NEAT assigns weights to the indicators based on the SAU and HBT that they represent (see Section Evaluation of the Assessment Results). The SAUs are weighted according to their surface area and the HBTs are weighed equally within a SAU. Therefore, the indicator values contribute to the assessment with different weights, the highest weight being assigned to an indicator representing a large SAU with a small number of indicators, and within it a HBT with a small number of indicators. The relative weights of the indicator values were used to identify the indicators that contribute 90% of the weight of the final assessment. In addition, the relative weight of each ecosystem component in each case study assessment was calculated. These summary statistics highlighted differences in aggregating information among case studies.

To test the sensitivity of the case study assessments to the selection and number of indicator values, a sensitivity analysis was performed by running the assessment using randomly selected indicator values. The number of indicator values included into the assessment varied from 1 to the maximum number of indicators in the case study minus one. This process was repeated 100 times for each number of indicator values. For example, take a case study with 120 indicator values. First, one random indicator value is selected and the assessment is done using only that indicator. This procedure is repeated 100 times. Then, two indicator values are picked at random, and the assessment is run using them; this again is repeated 100 times. This procedure is repeated for all numbers of indicator values up to 119. This results in a large number of values whose divergence can be analyzed to see if any patterns can be identified.

# RESULTS

## Assessment Design

The number of SAUs as well as how many hierarchical levels were used in these varied widely between the case studies. The number of SAUs included in the Gulf of Finland and Portugal continental sub-division cases were much higher (>60) than in all other case studies which included, on average, 9 different SAUs. Excluding these two case studies, larger areas were usually assessed using more SAUs. The number of hierarchical SAU levels varied between 1 and 5, but in 7 out of 10 cases, there were 3 or 4 levels (Table 2, Figure 2). The total number of HBTs included in the assessment varied between 3 and 9, and 9 out of 10 case studies had 2 or 3 hierarchical HBT levels (Table 2).

Not all SAUs necessarily included all habitat types, and indicators or data may not exist for all defined HBT types for each SAU. The number of SAU-HBT combinations that were assessed by at least one indicator value, varied between 6 and 132 (Table 2).

The number of ecosystem components included in the analyses varied between 5 and 9, with an average of 7.3 (Table 2). It has to be noted that all ecosystem components identified in this study were not applicable to all areas; an example being reptiles that do not occur in most of the study sites.

The number of unique indicators applied in each case study area varied between 11 and 116 (Table 2, Supplementary Material S11). The number of indicator values varied greatly with 466
<table><caption>TABLE 2 | Synthesis of the structure used by the different case studies for the nested assessment.</caption><thead><tr><th>Case study name</th><th>Area (km²)</th><th>Number of SAU levels (total number of SAU)</th><th>Number of HBT levels (total number of HBT)</th><th>Number of SAU*HBT combinations with data</th><th>Number of ecosystem components included</th><th>Number of unique indicators (see Supplement 11)</th><th>Number of indicator values</th></tr></thead><tbody><tr><td>Norwegian Barents Sea-Lofoten</td><td>821 478</td><td>4 (13)</td><td>3 (9)</td><td>21</td><td>7</td><td>40</td><td>74</td></tr><tr><td>Gulf of Finland</td><td>22 482</td><td>5 (60)</td><td>2 (3)</td><td>103</td><td>8</td><td>25</td><td>147</td></tr><tr><td>Lithuanian marine waters</td><td>6 426</td><td>2 (4)</td><td>2 (7)</td><td>6</td><td>9</td><td>27</td><td>50</td></tr><tr><td>Kattegat</td><td>17 440</td><td>3 (11)</td><td>2 (7)</td><td>21</td><td>8</td><td>31</td><td>69</td></tr><tr><td>Dutch North Sea</td><td>57 000</td><td>1 (1)</td><td>3 (6)</td><td>6</td><td>6</td><td>15</td><td>31</td></tr><tr><td>Basque Coast</td><td>10 794</td><td>3 (8)</td><td>3 (6)</td><td>22</td><td>9</td><td>48</td><td>109</td></tr><tr><td>Portuguese continental subdivision</td><td>268 645</td><td>4 (61)</td><td>4 (7)</td><td>132</td><td>7</td><td>14</td><td>466</td></tr><tr><td>Black Sea coast (Large Varna Bay)</td><td>1 434</td><td>3 (7)</td><td>2 (4)</td><td>15</td><td>7</td><td>35</td><td>112</td></tr><tr><td>Saronikos Gulf</td><td>2 907</td><td>3 (4)</td><td>3 (6)</td><td>10</td><td>7</td><td>17</td><td>29</td></tr><tr><td>Adriatic Sea</td><td>138 600</td><td>3 (10)</td><td>2 (5)</td><td>17</td><td>5</td><td>116</td><td>177</td></tr><tr><td>Mean</td><td>134 721</td><td>3.1 (17.9)</td><td>2.6 (6.0)</td><td>35.3</td><td>7.3</td><td>39.2</td><td>126.4</td></tr><tr><td>Stdev</td><td>255 807</td><td>1.1 (22.7)</td><td>0.7 (1.7)</td><td>44.3</td><td>1.3</td><td>30.7</td><td>128.8</td></tr></tbody></table>

SAU, spatial assessment unit; HBT, habitats. The case studies are ordered according to their latitude.

values at the higher end in Portugal continental sub-division and between 20 and 200 values in all other case studies (Table 2).

## Biological Diversity Status

The summary of the test NEAT assessments of the 10 case study areas is presented in Figure 3. The assessment resulted in GES for the Basque EEZ and the Barents Sea-Lofoten, with 100 and 66% confidence, respectively, the remaining eight case studies presented non-GES (i.e., bad, poor, or moderate; Figure 3). Lithuanian coast has the potential for being in GES, but with a low confidence of 20% (Figure 3). For the other case studies, this probability of achieving GES was <1% (Figure 3).

The different ecosystem components showed different status in the case study areas (Figure 4). No strong pattern among the ecosystem components could be detected, but some commonalities were found: Indicators based on marine mammals generally indicated degraded situation in 6 cases out of 7 (Figure 4). When included, birds and benthic vegetation indicators as well as water column indicators of physico-chemical status also indicated degraded situation in 5 cases out of 7. Indicators encompassing several ecosystem components (“AT,” on Figure 4) always indicated degraded situations. On the other hand, indicators of benthic habitats’ physico-chemical status and of zooplankton community status indicated GES when they were included in the assessment (Table 1, Figure 4).

## Relative Contribution of Indicator Values and Biodiversity Components

The indicator values contributed differently to the final assessment result (Figure 5); indicator values defined for larger SAUs tend to have more weight, particularly if there are only few indicators defined for these SAUs. In 7 out of the 10 case studies,

<10 indicator values already contributed to more than 50% of the final assessment result. For 9 case studies, <50 indicator values contributed to >90% of the final assessment. This 90% of the final assessment was reached with <20 indicator values in five case studies (Figure 5). The five indicator values that made the highest contribution to the final assessments of each case study are listed in Table 3. These indicator values were dominated by mammal, bird, fish, and benthic fauna indicators.

The 12 different ecosystem components’ contribution to the final assessment result did not correspond to the number of indicator values defined for each component (Table 4). For example, most case studies had a large proportion of benthic fauna indicator values (average: 22.4% of indicators values), which ultimately did not reflect proportionally in the final assessment (average contribution: 11.7%). In contrast, the proportion of fish and marine mammals indicator values were lower, but these components contributed to a higher proportion of the final assessment. In five case studies (i.e., Barents-Lofoten, Gulf of Finland, Dutch North Sea, Saronikos Gulf, and Adriatic Sea), “Benthic fauna” was the component with the highest proportion of indicator values (Table 4); the other five case studies each had a different component with the highest number of indicator values. However, in none of the case studies, benthic fauna was the component with highest contribution to the final assessment (Table 4); in five (i.e., Gulf of Finland, Dutch North Sea, Basque coast, Portuguese continental sub-division, and Black Sea coast) and two case studies (i.e., Barents Sea-Lofoten and Adriatic Sea) respectively, fish and mammals were the components carrying the highest weight to the final assessment (Table 4). However, other ecosystem components, that overall did not contribute to many case study assessments, were very relevant for specific case studies (e.g., the composite
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Probabilities for the five environmental status classes for each of the 10 case study assessments. Good environmental status is assumed attained if the cumulative probability of "Good" and "High" is higher than the cumulative probability of "Moderate," "Poor," and "Bad." If opposite, the Good environmental status is not attained. For case study codes see **Table 1**.</figcaption></figure>

group "all taxa" in the Saronikos Gulf and benthic habitat in the Lithuanian coast).

## Sensitivity Analysis

The sensitivity analysis shows that there are major differences in how much the result varies if only a subset of the indicator values is included in the assessment (Figure 6). For example, if only a small number (close to 0) indicators were included, the assessment results in all studies could be anywhere between high and bad status, except in Barents Sea and Portuguese continental subdivision, where they could range from poor to high status. As more indicator values are added, the range of outcomes narrows down. However, how steeply that happens when indicator values are added varies between the case study areas (Figure 6).

## DISCUSSION

The current NEAT-based assessment demonstrates a large-scale marine biodiversity assessment, providing a feasible solution to the apparent problem pointed out by the European Commission, in their evaluation of the EU member states' reports on the MSFD initial assessments carried out in 2010–2012 (Palialexis et al., 2014). This problem was the apparent lack of coherence and comparability in indicators used and in the final evaluation of the overall status between the countries and within all regional seas (Palialexis et al., 2014). Despite the available guidance and Commission Decision (European Union, 2010) on GES descriptors, criteria and indicators, the overall picture in assessments was patchy and non-coherent (European Commission, 2014). The use of NEAT, and its validation in different regional seas and case study areas, is a crucial contribution from the DEVOTES project to provide a harmonized approach and methodology for a coherent and comparable environmental status assessment across the European regional seas. It also shows that although the regional seas have different characteristics and human pressures impacting those (Claudet and Fraschetti, 2010; Micheli et al., 2013a; Andersen et al., 2015), a coherent assessment framework can be employed to evaluate differences in the environmental status and the ecological components that are impacted by different pressures.

The study and the comparison of the case studies brought into light several issues that need attention in order to improve the coherent and comparable "biodiversity status" assessments of the European regional seas. These issues are related to the data and indicator availability, how the assessments are structured, how the integrative assessment should be structured, and how this structure should be taken into account when defining the spatial resolution and indicator selection of the assessments. The current study revealed that while these assessments could be carried out, there are two major problems in achieving the objectives of the MSFD assessments: (i) there are still multiple gaps in the availability and coverage of indicators in the various areas, and (ii) comparability of the status assessments across different regions would benefit from a more unified assessment framework, even if indicators suitable for each area remained different. NEAT provides a general framework that could be accompanied with guidelines for the selection of SAUs, HBTs, and indicators.

Each of the case studies was initially designed with the best available selection of spatial units, habitats, and indicators, adhering to the NEAT methodology but without specific guidelines for the indicator selection, target level setting, etc. This situation resembles the situation where the new users would start using NEAT on their area. For the purposes of this study, the assessments were evaluated and harmonized to some degree, e.g., if the same indicator appeared in multiple case studies, it was ensured that it was associated to the same biodiversity component (e.g., chlorophyll a levels would be assigned to phytoplankton). Despite this harmonizing, there were major differences in how the case studies were constructed in terms of spatial resolution, habitats, and indicator definition. The current assessment is based on best available data and evaluation of the experts participating within this exercise, and the biodiversity status results of this study should be considered as indicative, not definitive.

The indicators selected for the assessments are designed or adapted for each area separately, including the geographical and habitat specification and the target level, i.e., which values are considered good and which less than good in any given area and habitat. This means that the "good" status is scaled according to the area: In areas with a naturally low biodiversity, lower
<figure><img src="image_7.png" /><figcaption>FIGURE 4 | Probabilities for the five environmental status classes separated by ecosystem component for each of the 10 case studies (for codes see Table 1). Ecosystem components considered are: Phytoplankton (PPL), zooplankton (ZPL), benthic vegetation (BV) benthic fauna (BF), fish, birds, mammals (Mamm.), reptiles (Rep.), pelagic fauna (PF), all taxa (AT), water column habitat (WCH), benthic habitat (BH).</figcaption></figure>

biodiversity is also considered “good” than in areas with naturally high diversity. This makes the assessment relevant for each area, and the result must be interpreted to be in relation to undisturbed condition of that area rather than in absolute terms of diversity.

According to a categorization of rules or methods for combining or aggregating indicators or criteria within a given descriptor (Prins et al., 2013; Borja et al., 2014), NEAT is classified as a high-level integration method which reduces the risks associated to the “one out, all out” principle of the Water Framework Directive approach (Borja and Rodríguez,

2010) while giving an overall and specific (to descriptors and components) assessment.

According to the relevant guidance document for the MSFD (Prins et al., 2013), the spatial scales are not the same for all indicators within the biodiversity descriptor, where depending on the species or habitat a different spatial scale may be used. It is also recommended to address uncertainties and assess confidence of the classification result (as a secondary assessment). In our study, the NEAT software treats equally all assessment elements assigning equal weights, but gives more weight in cases of larger spatial coverage, with higher data representativeness, in that way
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Cumulative contribution to the final assessment (cumulative weight) in relation to the number of indicator values. For each case study, values have been ranked from the most important indicator value (highest weight) to the least important value (lowest weight). The X-axis has been cut; Portuguese case study has a total of 466 indicator values. For case study codes, see **Table 1**.</figcaption></figure>

incorporating the spatial scales issue and the confidence level into the assessment. This could be the reason for which some ecosystem components (e.g., seabirds, mammals, and fishes) have more weight in the final assessment, since they are normally assessed at large scale spatial areas, which have more weight when aggregating (e.g., Saronikos gulf). However, NEAT also includes the possibility to weight indicators differently.

## Implications of the Assessment Design

Most of the case study areas lacked indicators regarding one or several biodiversity components and habitats (Table 1, Figure 4), even those that were deemed important in the area. The lack of indicators stemmed either from lack of monitoring data regarding the area or biological diversity component (e.g., birds, reptiles, pelagic fauna), or from obstacles in the indicator development, including the lack of expert time to develop indicators, or insufficient knowledge about the target levels due to lack of long-term or reference condition data (Hummel et al., 2015). In some cases, more basic ecological research is needed in order to understand the ecological processes well enough to develop indicators. In fact, most of the assessments undertaken until now by member states is more qualitative than quantitative (Hummel et al., 2015), representing a challenge for the assessment.

The habitats and biodiversity components for which no indicators are available potentially affect the final assessment result. It is entirely possible that adding even one indicator that would represent a poorly-represented, large area or habitat, would change the overall assessment for better or for worse. Therefore, in order to make a reliable assessment of the status of the biological diversity, the

critical gaps in each assessment case need to be evaluated for their potential to affect the overall result. If such high-leverage gaps exist, the assessment result must be taken with caution.

Different indicator values and spatial assessment units had varying weights in the final assessment result in all of the cases (Table 3, Figure 2). The differences in the indicator value weights stem from the fact that the default NEAT assessment first assesses the result for each SAU, giving equal weight to each HBT with similar hierarchy, and combines these SAUs hierarchically so that each SAU is given weight according to its area. Therefore, if a SAU has a large surface area and only a small number of indicators per one or several of its habitat types, these indicator values end up contributing strongly to the final assessment.

This emphasizes the importance of the balanced nature of the indicator set, and particularly the reliable assessment of indicators that are used to assess the status of large areas, and particularly their habitats with only few indicators (Feary et al., 2014). Therefore, particular attention should be paid to both the observed value, the boundary values between the classes, and the uncertainty estimation of these most influential indicator values.

The fact that the SAUs are weighted according to their surface area in the default mode of NEAT also emphasizes the need for careful consideration of the definition of the SAUs. Ideally, the SAUs should be defined in the manner that an indicator value defined for a SAU can be expected to reasonably represent all of the SAU. On the other hand, if the assessment area is split into several sub-SAUs and only a fraction of them actually has indicator data, their value will be generalized to represent the whole super-area in the hierarchical assessment anyway.
<table><caption>TABLE 3 | List of the top-five indicator values contributing the most to the overall assessment for each case study.</caption><thead><tr><th>Spatial Assessment Unit (SAU)</th><th>Habitat</th><th>Ecosystem component</th><th>Indicator name</th><th>% contribution to assessment</th></tr></thead><tbody><tr><td><strong>BARENTS SEA LOFOTEN</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>56</strong></td></tr><tr><td>Norwegian BARENTS SEA-LOFOTEN</td><td>sea-ice</td><td>Mammals</td><td>Harp seal, <em>Pagophilus groenlandicus</em> stock size</td><td>16.6</td></tr><tr><td>Ice-free zone</td><td>pelagic</td><td>Birds</td><td>Kittiwake, <em>Rissa tridactyla</em> breeding success over last 5 year</td><td>16.4</td></tr><tr><td>Seasonal-ice zone</td><td>sea-ice</td><td>Mammals</td><td>Proportion of non-threatened pagophiletic mammals</td><td>15</td></tr><tr><td>Offshore (ice-free zone)</td><td>shelf (muddy sediments)</td><td>Benthic Fauna</td><td>4 indicators related to macrobenthic fauna: ES100, abundance level, evenness and AMBI</td><td>3.7 (each)</td></tr><tr><td><strong>GULF OF FINLAND</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>60</strong></td></tr><tr><td>Open sea</td><td>Benthic</td><td>Benthic Fauna</td><td>Average regional diversity</td><td>18.4</td></tr><tr><td>Gulf of Finland</td><td>Pelagic</td><td>Fish</td><td>Abundance of salmon spawners and smolt</td><td>12.5</td></tr><tr><td>Gulf of Finland</td><td>Pelagic</td><td>Fish</td><td>Herring, Spawning stock biomass</td><td>12.5</td></tr><tr><td>Gulf of Finland</td><td>Marine (general)</td><td>Mammals</td><td>3 indicators: Gray seal population growth rate, Gray seal pregnancy rate, Ringed seal population growth rate</td><td>8.3 (each)</td></tr><tr><td><strong>LITHUANIAN COAST</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>47</strong></td></tr><tr><td>Territorial sea</td><td>Benthic</td><td>Benthic habitat</td><td>Extent of the seabed significantly affected by human activities</td><td>31</td></tr><tr><td>Territorial sea</td><td>Marine (general)</td><td>Others</td><td>Biopollution level index (invasive species)</td><td>4</td></tr><tr><td>Territorial sea</td><td>Marine (general)</td><td>Birds</td><td>Abundance of wintering populations of seabirds: (1) Red-throated Diver + Black-throated Loon, (2) Great Crested Grebe (<em>Podiceps cristatus</em>), (3) Common Merganser (<em>Megus merganser</em>), (4) Velvet Scotter (<em>Melanitta fusca</em>), (5) Long-tailed Duck (<em>Clangula hyemalis</em>), and (6) Common Goldeneye (<em>Bucephala clangula</em>)</td><td>4 (each)</td></tr><tr><td><strong>KATTEGAT DK</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>30</strong></td></tr><tr><td>KATTEGAT, central parts</td><td>Marine (general)</td><td>Birds</td><td>Fulmar winter abundance (encounter rate)</td><td>6.6</td></tr><tr><td>KATTEGAT, central parts</td><td>Marine (general)</td><td>Birds</td><td>Kittiwake winter abundance (encounter rate)</td><td>6.6</td></tr><tr><td>KATTEGAT, central parts</td><td>Marine (general)</td><td>Birds</td><td>Guillemot winter abundance (encounter rate)</td><td>6.6</td></tr><tr><td>KATTEGAT, central parts</td><td>Benthic</td><td>Benthic Fauna</td><td>BQI</td><td>5</td></tr><tr><td>KATTEGAT, central parts</td><td>Benthic</td><td>Benthic Fauna</td><td>DKI</td><td>5</td></tr><tr><td><strong>DUTCH NORTH SEA</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>21</strong></td></tr><tr><td>Dutch EEZ</td><td>Benthic</td><td>Benthic Fauna</td><td>Benthic invertebrates total number of species</td><td>16.7</td></tr><tr><td>Dutch EEZ</td><td>Muddy deep bottom</td><td>Benthic Fauna</td><td>3 indicators related to typological group sensitive to seafloor physical impact, based on (1) density, (2) biomass, (3) number of species</td><td>1.1 (each)</td></tr><tr><td>Dutch EEZ</td><td>Muddy deep bottom</td><td>Benthic Fauna</td><td>4 indicators related to typological group highly sensitive to seafloor physical impact, based on (1) density, (2) biomass, (3) number of species</td><td>1.1 (each)</td></tr><tr><td><strong>BASQUE EEZ</strong></td><td></td><td></td><td><strong>Cumulative contrib. of top-five indicators</strong></td><td><strong>33</strong></td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Benthic</td><td>Benthic habitat</td><td>Seabed affected by human activities</td><td>9.9</td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Pelagic</td><td>Birds</td><td>Biological value Seabirds</td><td>6.6</td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Pelagic</td><td>Mammals</td><td>Biological value Mammals</td><td>6.6</td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Pelagic</td><td>Phytoplankton</td><td>Eutrophication indicator: Chlorophyll a, 90th percentile</td><td>6.6</td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Sedimentary</td><td>Benthic Fauna</td><td>2 indicators: (1) M-AMBI, (2) AMBI</td><td>3.3 (each)</td></tr><tr><td>Offshore waters (>200 m depth)</td><td>Sedimentary</td><td>Fish</td><td>Biological value Demersal Fish</td><td>3.3</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 3 | Continued</caption><thead><tr><th>Spatial Assessment Unit (SAU)</th><th>Habitat</th><th>Ecosystem component</th><th>Indicator name</th><th>% contribution to assessment</th></tr></thead><tbody><tr><td colspan="5"><strong>PORTUGUESE CONTINENTAL SUBDIVISION</strong></td></tr><tr><td>Continental_A2_600</td><td>Marine (general)</td><td>Birds</td><td><strong>Cumulative contrib. of top-five indicators</strong><br>Biological Value Marine Birds</td><td>17<br>3.9</td></tr><tr><td>Continental_A2_600</td><td>Pelagic</td><td>Zooplankton</td><td>Biological Value Zooplankton</td><td>3.9</td></tr><tr><td>Continental_B4_600</td><td>Sedimentary</td><td>Benthic Fauna</td><td>Biological Value Benthic communities</td><td>3.7</td></tr><tr><td>Continental_B4_600</td><td>Marine (general)</td><td>Birds</td><td>Biological Value Marine Birds</td><td>3.7</td></tr><tr><td>Continental_B5_600</td><td>Marine (general)</td><td>Birds</td><td>Biological Value Marine Birds</td><td>2.2</td></tr><tr><td colspan="5"><strong>COASTAL BLACK SEA</strong></td></tr><tr><td>Black Sea coastal</td><td>Pelagic</td><td>Fish</td><td><strong>Cumulative contrib. of top-five indicators</strong><br>4 indicators: (1) Mean length of Sprattus sprattus, (2) Catch/biomass ratio of S. sprattus, Biomass of S. sprattus, (4) Sexually mature specimen of S. sprattus</td><td>57<br>11.4 (each)</td></tr><tr><td>Black Sea coastal</td><td>Benthic</td><td>Fish</td><td>4 indicators: (1) Mean length of Scophtalmus maximus, (2) Catch/biomass ratio of S. maximus, (3) Biomass of S. maximus, (4) Sexually mature specimen of S. maximus</td><td>11.4 (each)</td></tr><tr><td colspan="5"><strong>SARONIKOS GULF</strong></td></tr><tr><td>Saronikos Gulf</td><td>Marine (general)</td><td>All Taxa</td><td><strong>Cumulative contrib. of top-five indicators</strong><br>CIMPAL index (alien species)</td><td>62<br>44</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Reptiles</td><td>% loss of spawning areas of sea turtle Caretta caretta</td><td>4.4</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Fish</td><td>4 indicators: Fisging mortality for (1) Engraulis encrasicolus, (2) Sardina pilchardus, (3) Merluccius merluccius, (4) Mullus barbatus</td><td>4.4 (each)</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Mammals</td><td>% Threatened mammals</td><td>4.4</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Fish</td><td>% Threatened sharks</td><td>4.4</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Fish</td><td>% of stocks that meet GES based on fishing mortality</td><td>4.4</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Fish</td><td>% of stocks that meet GES based on reproductive capacity</td><td>4.4</td></tr><tr><td>Saronikos Gulf</td><td>Pelagic</td><td>Fish</td><td>% of stocks that meet GES based on reproductive capacity and biomass indices</td><td>4.4</td></tr><tr><td colspan="5"><strong>ADRIATIC SEA</strong></td></tr><tr><td>Adriatic sea</td><td>Pelagic</td><td>Mammals</td><td><strong>Cumulative contrib. of top-five indicators</strong><br>Tursiops truncatus, distributional range</td><td>58<br>11.7</td></tr><tr><td>Adriatic sea</td><td>Pelagic</td><td>Mammals</td><td>Tursiops truncatus, distributional range</td><td>11.7</td></tr><tr><td>Adriatic sea</td><td>Pelagic</td><td>Mammals</td><td>Grampus griseus, distributional range</td><td>11.7</td></tr><tr><td>Adriatic sea</td><td>Pelagic</td><td>Mammals</td><td>Ziphius cavirostris, distributional range</td><td>11.7</td></tr><tr><td>Adriatic sea</td><td>Pelagic</td><td>Mammals</td><td>Tursiops truncatus, abundance</td><td>11.7</td></tr></tbody></table>

In case of equal contribution of several indicator values, all the indicator values are given. The contribution to the overall assessment (in %) of each indicator value is given. Numerical values are rounded. ES100, expected number of species in 100 individuals; AMBI, AZTI's Marine Biotic Index; BQI, Benthic Quality Index; DKI, Danish Index; CIMPAL, Cumulative IMPacts of invasive ALien species; M-AMBI, multivariate AMBI; GES, good environmental status.

In NEAT, it is possible to weight the SAUs according to their perceived ecological relevance instead of their surface area; for example, biodiversity hotspots, important reproduction areas, marine protected areas, etc., could be given a higher weight than their area alone would imply. In this study, this option was not used in any of the case studies.

Uncertainty of the results is assessed based on Monte Carlo simulations, using the observed value as mean and the standard error value as the standard deviations, assuming a Gaussian distribution (Carstensen and Lindegarth, 2016). Based on these simulations, NEAT determines how often the sampled value falls into each of the five classes, and this distribution is reported. Therefore, the standard error values assigned to the indicators

play a major role in the uncertainty associated with the final assessment result. This emphasizes the importance of careful evaluation of the standard deviation, particularly with indicators that have a high weight in the assessment.

## Evaluation of the Assessment Results

There are other tools to assess the status of marine systems, e.g., the Ocean Health Index (OHI; Halpern et al., 2012). This index has different concept and a much broader spatial scale, and a comparison between NEAT and OHI results (BD values presented in Table S6 in Selig et al., 2013) shows that the results are quite different (Table 5).
<table><caption>TABLE 4 | Number of indicator values per ecosystem component, and the relative weight of ecosystem components in each case study.</caption><thead><tr><th rowspan="2"></th><th colspan="12">Number of indicator values</th><th colspan="6">% of total weight</th></tr><tr><th>Phytoplankton</th><th>Zooplankton</th><th>Benthic Veg</th><th>Benthic Fauna</th><th>Fish</th><th>Birds</th><th>Mammals</th><th>Reptiles</th><th>Pelagic Fauna</th><th>All Taxa</th><th>Water column habitat</th><th>Benthic habitat</th><th>Phytoplankton (%)</th><th>Zooplankton (%)</th><th>Benthic Veg (%)</th><th>Benthic Fauna (%)</th><th>Mammals (%)</th><th>Reptiles (%)</th><th>Pelagic Fauna (%)</th><th>All Taxa (%)</th><th>Water column habitat (%)</th><th>Benthic habitat (%)</th></tr></thead><tbody><tr><td>Norwegian Barents</td><td>1</td><td>4</td><td>3</td><td>37</td><td>9</td><td>14</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>0</td><td>20</td><td>14</td><td>28</td><td>34</td><td></td><td></td><td></td></tr><tr><td>Sea-Lofoten</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gulf of Finland</td><td>53</td><td>2</td><td>34</td><td>46</td><td>4</td><td>2</td><td>3</td><td>0</td><td>0</td><td>0</td><td>3</td><td>0</td><td>10</td><td>4</td><td>2</td><td>20</td><td>30</td><td>3</td><td>25</td><td></td><td></td><td>6</td></tr><tr><td>Lithuanian marine waters</td><td>3</td><td>6</td><td>1</td><td>1</td><td>6</td><td>19</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>8</td><td>3</td><td>4</td><td>9</td><td>0</td><td>0</td><td>1</td><td>29</td><td></td><td>5</td><td>20</td><td>32</td></tr><tr><td>Kattegat</td><td>5</td><td>1</td><td>25</td><td>9</td><td>11</td><td>3</td><td>1</td><td>0</td><td>0</td><td>0</td><td>14</td><td>0</td><td>2</td><td>0</td><td>26</td><td>14</td><td>16</td><td>20</td><td>1</td><td></td><td></td><td>21</td><td></td></tr><tr><td>Dutch North Sea</td><td>2</td><td>2</td><td></td><td>19</td><td>6</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>7</td><td>7</td><td></td><td></td><td>33</td><td>45</td><td>4</td><td>4</td><td></td><td></td><td></td></tr><tr><td>Basque Coast</td><td>8</td><td>5</td><td>5</td><td>19</td><td>42</td><td>6</td><td>6</td><td>0</td><td>0</td><td>0</td><td>12</td><td>6</td><td>8</td><td>1</td><td>0</td><td>10</td><td>54</td><td>8</td><td>8</td><td>8</td><td></td><td>25</td><td>1</td><td>11</td></tr><tr><td>Portuguese continental subdivision</td><td>44</td><td>42</td><td></td><td>27</td><td>108</td><td>44</td><td>0</td><td>0</td><td>177</td><td>24</td><td>0</td><td>0</td><td>8</td><td>9</td><td></td><td></td><td>8</td><td>29</td><td>21</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Black Sea coast (Large Varna Bay)</td><td>36</td><td>24</td><td>4</td><td>12</td><td>12</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>20</td><td>4</td><td>2</td><td>1</td><td>2</td><td>2</td><td>92</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Saronikos Gulf</td><td>0</td><td>0</td><td>2</td><td>12</td><td>8</td><td>0</td><td>1</td><td>1</td><td>0</td><td>1</td><td>4</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Adriatic Sea</td><td>3</td><td>0</td><td>0</td><td>100</td><td>0</td><td>0</td><td>30</td><td>0</td><td>0</td><td>0</td><td>37</td><td>7</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The OHI tends to give a more reduced range of status values (74–97) than those provided by NEAT (0.37–0.69) for these areas. The OHI does not provide a GES/non-GES status, but in general provides higher values than those by NEAT. The OHI study (Selig et al., 2013) has been applied globally, and includes a large variety of worldwide cases with great differences in setting and problems. In that context, e.g., the Mediterranean and the Baltic Sea seem to be in a (seemingly more homogenous) better state than e.g., waters around Africa or Indonesia and Philippines.

An interesting observation is that there is a negative rather than a positive correlation between these results, and those areas ranked low in NEAT (such as the Gulf of Finland and Kattegat) get high scores in OHI, while the best-scoring area in NEAT (Basque EEZ) gets lowest score in OHI (Table 5). This discrepancy is partly due to the fact that the OHI scores are given by country, thus covering larger areas than the case studies assessed here with NEAT. Therefore, the local status of a case study area may be masked by the results from the rest of the country in OHI. The NEAT results are reported here for the entirety of each of the case study areas, but where the case study area includes smaller SAUs, the results can be viewed for each of them separately as well, yielding even a more detailed geographical resolution.

Another factor possibly contributing to this discrepancy is the use of different indicators; the OHI assessment used publicly available data with little local/regional detail, which can vary the final assessment when applying to regional scales (Halpern et al., 2014), while the current NEAT assessment used indicators specifically designed for marine status assessment. The species scores of OHI focused on the extinction risk of marine species (Selig et al., 2013), while the indicators in the NEAT assessments included a wider spectrum of indicators of species status. The OHI habitat scores were based on condition estimates of mangroves, coral reefs, seagrass beds, salt marshes, sea ice, and subtidal soft-bottom (Selig et al., 2013) while the NEAT assessments were tailored for each area.

The NEAT assessment results were in most cases in line with previous regional/local assessments, understanding, or known pressure gradients (Table 1, Figure 4). For example, The Baltic Sea biodiversity has been assessed by HELCOM (2009a, 2010) to be in poor to bad status in all of the three Baltic case study areas included in this analysis (Gulf of Finland, Lithuanian marine waters, Kattegat), being similar to the NEAT results but not to the OHI assessments. The difference between the NEAT and OHI results in these cases is probably largely due to eutrophication, which is documented to be major pressure threatening the ecosystem functioning of the Baltic Sea (HELCOM, 2009b, 2010). While it is reflected in the status of phytoplankton and water column habitats, and also affects the higher trophic levels of the food web (Österblom et al., 2007) and the seafloor (Karlson et al., 2002), it is not likely to be strongly reflected in the extinction threat of marine species (used in OHI), although it does affect the habitat scores, particularly seagrasses (Table S1 in Selig et al., 2013). Another factor affecting the discrepancy in the case of Finland is that the Gulf of Finland area has poorer biodiversity status than the Finnish marine waters on average (HELCOM, 2010).
<figure><img src="image_9.png" /></figure>

**FIGURE 6** | Variation of the overall assessment results in the case studies, if only a subset of the indicator values is taken into account. The x axis indicates how many indicators are included into the assessment, and the gray area spans the assessment result values that emerged. The black line shows the mean assessment result across the 100 runs conducted for each number of indicator values. Note that the x axis goes only up to 200 indicator values; the Portuguese continental subdivision study included 466 indicator values in total.

In the North Sea, fishing is considered the main pressure, and the results show fish to be the ecosystem component in poorest status; the other assessed ecosystem components (birds, mammals, benthic fauna, and phytoplankton) were assessed to be in GES, with the exception of zooplankton that showed sub-GES (moderate) status (Figure 4). The Black Sea Coast case results obtained in this study also corresponded very well to known pressure gradients, such as nutrient enrichment affecting the status of the plankton community (Figure 4). Phytoplankton and benthic vegetation assessments correspond to category “poor” in the Varna Bay itself (Dencheva and Doncheva, 2014; Moncheva

et al., 2015) as the most affected by anthropogenic pressure among the BSC sub-SAUs (Shtereva et al., 2012). The lowest benthic fauna score is also found there, which is fully in compliance with recently published results (National Report on the State and Protection of the Environment in Bulgaria, 2014). Similarly, the Basque area, which was previously assessed as being in good status, using a different methodology (Borja et al., 2011) also results in good status after applying NEAT; only mammals were assessed to be in sub-GES status (Figure 4).

In Saronikos Gulf the assessment results correspond to the ecological status categorization according to the WFD which is
<table><caption>TABLE 5 | Comparison of the biodiversity assessments obtained using the Ocean Health Index (OHI; data from Selig et al., 2013) and the Nested Environmental status Assessment Tool (NEAT) (this study) in the countries for which NEAT has case studies.</caption><thead><tr><td>Country (Case study)</td><td>Status (NEAT)</td><td>Biodiversity score (OHI)</td></tr></thead><tbody><tr><td>Norway (Barents Sea)</td><td>0.646</td><td>90</td></tr><tr><td>Finland (Finland)</td><td>0.401</td><td>97</td></tr><tr><td>Lithuania (Lithuania)</td><td>0.583</td><td>89</td></tr><tr><td>Denmark (Kattegat)</td><td>0.394</td><td>93</td></tr><tr><td>Netherlands (North Sea)</td><td>0.523</td><td>85</td></tr><tr><td>Spain (Basque Country)</td><td>0.689</td><td>74</td></tr><tr><td>Portugal</td><td>0.566</td><td>83</td></tr><tr><td>Bulgaria (Bulgaria)</td><td>0.495</td><td>87</td></tr><tr><td>Greece (Saronikos Gulf)</td><td>0.520</td><td>91</td></tr><tr><td>Italy/Croatia (Adriatic)</td><td>0.370</td><td>86/89</td></tr></tbody></table>

poor in the sewage outfall area and moderate in the inner central gulf (Simboura et al., 2014, 2015, 2016). Aliens, fish including threatened sharks, and mammals contributed to the moderate status seen for the outer Saronikos and overall Saronikos. In general, the respective assessment results, although not definitive, are in line with pertinent studies (Frantzis, 2009; Katsanevakis et al., 2013; Papaconstantinou, 2014; Vasilakopoulos et al., 2015; Zenetos et al., 2015; Simboura et al., 2016) regarding the Greek marine waters. The Saronikos Gulf result obtained in this analysis was lower than the OHI assessment of the Greek waters, which was to be expected, as the Gulf is intensely exploited.

Results from the Norwegian part of the Barents sea indicated a general good status, which is in accordance with indicators of fish status on exploited large marine ecosystems (Kleisner et al., 2014; Coll et al., 2015), the report on the Barents Sea management plan (Sunnana et al., 2010) and the work from Certain et al. (2011). Nevertheless, several indicators indicated potentially degraded situations both in the coastal area and in the area of seasonal ice presence: (1) Along northern Norway coast, the current extent of kelp forest, an important component of fjords ecosystem and coastal landscape, cannot be considered as good in northern Norway. Kelp forests along the Norwegian and Russian coast were indeed dramatically grazed during the early 1970s and replaced by barren grounds dominated by sea urchins (Norderhaug and Christie, 2009). Though a progressive northward recovery of kelp forests extent is observed, its recovery status is still partial in northern Norway (Sivertsen, 2006; Rinde et al., 2014). (2) In northernmost part of the Barents sea, sea-ice extent is undergoing a particularly dramatic decrease (Parkinson et al., 1999) with a significant decrease rate of -3.5% per decade of winter ice extent (Sorteberg and Kvingedal, 2006) as a response to climate warming (Boitsov et al., 2014). This dramatic loss of habitat has consequences on the associated communities (Kovacs et al., 2011) as well as in the functioning of the Barents sea ecosystem as a whole (Wassmann et al., 2006). The growing evidence of impacts of climate change on this area rises the issue of exogenic unmanaged pressures on this system and the issue of shifting baselines for the definition of target values. In addition,

there are still no indicators of the impact of trawling activities included in this assessment (see however Jørgensen et al., 2016).

For the Portuguese coast, the initial assessment officially provided in the scope of the MSFD (MAMAOT, 2012), presented a general environmental quality status higher than the NEAT results calculated in this study. This may be partly due to the fact that the present assessment did not include some special areas with a higher degree of protection (such as Berlengas' Marine Reserve and Professor Luiz Saldanha's Marine Park or Goringe Seafloor). These areas, which have restricted access by the public, are important for marine high trophic level species (e.g., marine birds, mammals), some of which were not included in the present assessment. Due to inconsistencies in the data (now being improved by projects such as MARPRO—Conservation of Marine Protected Species in mainland Portugal, http://marprolife.org), marine mammals, reptiles and benthic vegetation were not included in the current NEAT assessment, which may also contribute to the lower environmental quality results achieved by NEAT. The higher result reported by the OHI may be related to the methodology used for the scores' calculation, and may reflect more specifically the trend than the present environmental status.

An exception to the good correspondence between the current and previous assessments is the Adriatic Sea, where the assessment provided by NEAT appears too low considering the current trends, also reported in the scientific literature, and available information from expert opinions (Coll et al., 2010; Bastari et al., 2016). Despite the historical impacts on this shallow water basin, the Adriatic Sea is still characterized by a wide diversity of habitats, including rocky and soft bottoms, large estuaries and lagoons, seagrass meadows and in, its southern part, also deep-water environments. The habitat richness is reflected by a high biodiversity (Coll et al., 2012; Micheli et al., 2013b), with approximately 49% of the species described for the Mediterranean Sea (Boudouresque et al., 2009; UNEP, 2015) and a variety of endemic species (e.g., 18% of the endemic fish species of the Mediterranean; UNEP/MAP-RAC/SPA, 2015). Human activities and multiple stressors, and in particular bottom trawling, hydraulic dredging and habitat loss, are certainly still impacting the Adriatic Sea (Micheli et al., 2013a; Pusceddu et al., 2014). However, the overall environmental condition is not worsening with respect to the past decade. Eutrophication and dystrophic crises, related to the high nutrient discharge from the Po River combined with an alteration in water circulation, have caused hypoxia, anoxia and massive mucilage events, with consequent mortality of the benthic organisms, but the frequency of these events decreased significantly (or even disappeared) in the last decade (Degobbis et al., 2000; Danovaro et al., 2009). Thus, we hypothesize that the assessment of the environmental status obtained by using NEAT can be affected by the number and typology of data included in the specific exercise. An improvement of the number and type of the biological indicators (e.g., species or ecosystem functioning) could be crucial to obtain a more realistic classification of the marine environmental health of the Adriatic Sea.

Birds and mammals were found to be in poor status in many of the case study areas. This reflects the fact that seabirds are
indeed considered as more threatened than any other comparable
groups of bird species in general and display a faster trend of
decline than other bird species during the last decades (Croxall
et al., 2012). In addition, using IUCN Red list categories, it has
been evidenced that, among seabirds, pelagic species of seabirds
are disproportionately more threatened than coastal resident
or coastal non-breeding visitor species (Croxall et al., 2012).
Pelagic seabirds are particularly sensitive to disturbance as most
species lay only a single egg, adults do not reproduce every
year and usually reproduce several years after reaching sexual
maturity (Furness and Camphuysen, 1997). Most seabird species
display very large home range and thus integrate the state of the
environment and impacts of pressures over larger scale.

The conservation status of marine mammals is of particular
concern with an estimated proportion of threatened species
ranging worldwide between 23 and 61% of species (Schipper
et al., 2008). The North Atlantic region, which includes several of
the cases studied here, is one of the areas where the proportion of
threatened marine mammals is the highest, as shown by the low
quality values in Barents Sea, Kattegat, and Basque case studies
(Figure 4). The main reported threats explaining the bad status
of marine mammals are a long history of harvesting, accidental
mortalities through bycatch and collisions with vessel as well
as a very large panel of pollutions (from sound pollution to
contaminants and marine debris) and climate change (Schipper
et al., 2008). The sensitivity of these species to changes in their
environment might be related to their very slow population
dynamics, low densities in correlation with their large body-
size (Cardillo et al., 2008). Those life traits are also related with
relatively large home range. As a consequence, indicators of
marine mammals are usually measured over large scale, and
they are difficult to monitor with precision, leading to higher
uncertainty on many indicators (Taylor et al., 2007).

In two of the areas (Lithuania and Basque coasts), the
indicator contributing the most to the final assessment was “the
extent of the seabed significantly affected by human activities,”
which is a direct indicator of pressure. This is interesting since
some authors (e.g., Borja et al., 2013) have supported the use of
pressures instead of assessing the environmental status, if there
are not enough indicators. This should be done under the premise
that if an area has no obvious pressures then any changes in the
area must be due to natural changes which are outside the control
of management and vice versa.

Sensitivity Analysis

The sensitivity analysis results show differences among the case
studies in terms of how many indicator values are needed
before the assessment results will show approximately the same
results regardless of which indicator values are selected into the
assessment (Figure 6). This implies that there is no universally
sufficient number of indicator values needed to make a reliable
assessment, but that the number varies among case studies. No
clear patterns could be found among the 10 cases evaluated
in this study that would indicate a number of indicator values
of biodiversity components that can be considered sufficient
regardless of the case study and its structure.

The variation in the assessment result depends on the set
of indicator values that is available for the assessment. If the
indicator values are close to each other, i.e., all indicating similar
status, the variation in the results is naturally smaller. In contrast,
if the different indicator values indicate very different status, e.g.,
some areas or biodiversity components are in good status while
others are in bad, this naturally incurs a larger variation when
a subset of these variables are selected, as e.g., in the Gulf of
Finland.

These observations lead to the conclusions that if there is
variation among the status of the geographical or biodiversity
components in the study area, all of them should be covered by
indicators if possible. Particularly the inclusion of high-leverage
indicator values, i.e., those that have high weight and whose value
differs from the overall mean, can change the assessment result.
Therefore, the careful evaluation of the value and class limits of
these indicators should be a priority.

CONCLUSIONS

The structured assessment forces us to critically evaluate
the available indicator set in terms of ecological and spatial
representativeness of each indicator. This framework highlights
the gaps in the assessment as well as those parts that are well-
represented by current monitoring and available indicators. This,
in turn, helps in determining the best way to improve the quality
of the assessment: (i) via developing additional indicators to fill
in the gaps within the ecosystem approach (i.e., if not all the
important trophic levels of key species/ groups are covered in
the existing indicator set), (ii) working to determine the optimal
SAU for different categories of indicators that are targeted to
assess various trophic levels and functions in the food web, as
well as the HBT classification for each area, and (iii) working
toward improving specificity, robustness, and pressure relevance
of the indicators and enabling estimation of their standard
errors.

The development of NEAT and this extensive testing with
10 case studies in very different European marine areas offers
insight both to the status of the marine waters and to the state-
of-the art of the available indicator assemblages as well as the
development needs of the marine biological diversity assessment.
The application of the tool will make the improvement and
harmonization needs of the assessments visible and pave the way
toward a harmonized assessment across large geographical scales.

In conclusion, we propose the following recommendations
for the best practice in performing the environmental status
assessment using NEAT:

- Careful attention needs to be paid particularly to the current status and class boundaries of the indicators that cover large geographical areas (such as mobile birds and mammals), as they tend to carry a lot of weight in the final assessment.

- In order to make the assessment comparable between the different sub-regions and areas in the regional seas and provide a harmonized assessment among the regional seas, the design of the assessment needs to be harmonized. Attention must be paid to the selection of ecosystem components, and definition
of size and hierarchy of the spatial assessment units as well as
the definition of habitats.

- Consider the possibility of using different weighting for the individual indicator values, if that is ecologically more justified than using the weight based on the spatial area and habitat weighting.

- Contextualize the outputs on the basis of existing data. Different ecosystem components may present quite different data coverage, frequency, and data quality for the evaluation, and that may be reflected in the results. Consider carefully the standard deviation assigned to the indicators, but also consider how well the available indicators represent the ecosystem component and/or area as a whole.

- Consider not only the overall assessment, but the partial assessments (e.g., biological components or MSFD descriptors), as partial assessments can contribute to increased understanding of results and defining management measures for specific issues or areas.

AUTHOR CONTRIBUTIONS

LU, HB, and AB conceived the paper. The following partners provided the case studies: HB, SC (Barents Sea–Lofoten), LU, VK, HN (Gulf on Finland), GS (Lithuanian marine waters), CM, JA, JC (Kattegat), OB (Dutch North Sea), AB, MU (Basque Coast), JN (Portuguese continental subdivision), SM (Black Sea Coast), NP, MP, NS (Saronikos Gulf), RD, LC, AC, SB (Adriatic Sea). LU wrote the first draft. TB and CM contributed to the

calculations and sensitivity analyses. LU, HB, AB, SC, AH, and
MU contributed largely to the introduction and discussion. All
authors contributed to the last draft and to the discussions.

FUNDING

This manuscript is a result of DEVOTES (DEVelopment Of innovative Tools for understanding marine biodiversity and assessing good Environmental Status) project, funded by the European Union under the 7th Framework Programme, “The Ocean of Tomorrow” Theme (grant agreement no. 308392), www.devotes-project.eu. MU is partially funded through the Spanish programme for Talent and Employability in R+D+I “Torres Quevedo.” Moreover, the monitoring of Saronikos Gulf was financed by the Athens Water Supply and Sewerage Company (EYDAP SA).

ACKNOWLEDGMENTS

The authors would like to thank DEVOTES Advisory Board
members Paul Snelgrove and Simon Greenstreet for their
constructive comments at the early phase of the work.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: http://journal.frontiersin.org/article/10.3389/fmars.
2016.00159

REFERENCES

Andersen, J. H., Dahl, K., Göke, C., Hartvig, M., Murray, C., Rindorf, A., et al. (2014). Integrated assessment of marine biodiversity status using a prototype indicator-based assessment tool. *Front. Mar. Sci.* 1:55. doi: 10.3389/fmars.2014.00055

Andersen, J. H., Halpern, B. S., Korpinen, S., Murray, C., and Reker, J. (2015). Baltic Sea biodiversity status vs. cumulative human pressures. *Estuar. Coast. Shelf Sci.* 161, 88–92. doi: 10.1016/j.ecss.2015.05.002

Andersen, J.H., and Stock, A., Heinänen, S. M., Mannerla, and Vinther, M. (2013). *Human uses, Pressures and Impacts in the Eastern North Sea. Aarhus University, DCE-Danish Centre for Environment and Energy*, Vol. 18. Technical Report from DCE-Danish Centre for Environment and Energy, 134.

Artegiani, A., Bregant, D., Paschini, E., Pinardi, N., Raicich, F., and Russo, A. (1997a). The Adriatic Sea general circulation. Part I: air-sea interactions and water mass structure. *J. Phys. Oceanogr.* 27, 1492–1514.

Artegiani, A., Bregant, D., Paschini, E., Pinardi, N., Raicich, F., and Russo, A. (1997b). The Adriatic Sea general circulation. Part II: baroclinic circulation structure. *J. Phys. Oceanogr.* 1515–1532. [Epub ahead of print].

Bastari, A., Micheli, F., Ferretti, F., Pusceddu, A., and Cerrano, C. (2016). Large marine protected areas (LMPAs) in the Mediterranean Sea: the opportunity of the Adriatic Sea. *Mar. Policy* 68, 165–177. doi: 10.1016/j.marpol.2016.03.010

Beaumont, N. J., Austen, M. C., Atkins, J. P., Burdon, D., Degraer, S., Dentinho, T. P., et al. (2007). Identification, definition and quantification of goods and services provided by marine biodiversity: Implications for the ecosystem approach. *Mar. Pollut. Bull.* 54, 253–265. doi: 10.1016/j.marpolbul.2006.12.003

Berg, T., Murray, C., Carstensen, J., and Andersen, J.H. (2016). *Manual, Guidelines and Software for Biodiversity Assessment*. DEVOTES Deliverable D6.3. Available online at: http://www.devotes-project.eu/download/25439/

Birk, S., Bonne, W., Borja, A., Brucet, S., Courrat, A., Poikane, S., et al. (2012). Three hundred ways to assess Europe's surface waters: an almost complete

overview of biological methods to implement the water framework directive.
*Ecol. Indic.* 18, 31–41. doi: 10.1016/j.ecolind.2011.10.009

Boitsov, V. D., Karsakov, A. L., and Trofimov, A. G. (2014). Atlantic water temperature and climate in the Barents Sea, 2000–2009. *ICES J. Mar. Sci.* 69, 833–840. doi: 10.1093/icesjms/fss075

Borja, A., and Collins, M. (2004). *Oceanography and Marine Environment of the Basque Country*, Vol. 70. Elsevier Oceanography Series, 616.

Borja, Á., and Rodriguez, J. G. (2010). Problems associated with the ‘one-out, all-out’ principle, when using multiple ecosystem components in assessing the ecological status of marine waters. *Mar. Pollut. Bull.* 60, 1143–1146. doi: 10.1016/j.marpolbul.2010.06.026

Borja, Á., Dauer, D. M., and Grémare, A. (2012). The importance of setting targets and reference conditions in assessing marine ecosystem quality. *Ecol. Indic.* 12, 1–7. doi: 10.1016/j.ecolind.2011.06.018

Borja, A., Elliott, M., Andersen, J. H., Berg, T., Carstensen, J., Halpern, S., et al. (2016). Overview of integrative assessment of marine systems: the ecosystem approach in practice. *Front. Mar. Sci.* 3:20. doi: 10.3389/fmars.2016.00020

Borja, A., Elliott, M., Andersen, J. H., Cardoso, A. C., Carstensen, J., Ferreira, J. G., et al. (2013). Good Environmental Status of marine ecosystems: What is it and how do we know when we have attained it? *Mar. Pollut. Bull.* 76, 16–27. doi: 10.1016/j.marpolbul.2013.08.042

Borja, A., Galparsoro, I., Irigoien, X., Iriondo, A., Menchaca, I., Muxika, I., et al. (2011). Implementation of the European marine strategy framework directive: a methodological approach for the assessment of environmental status, from the Basque Country (Bay of Biscay). *Mar. Pollut. Bull.* 62, 889–904. doi: 10.1016/j.marpolbul.2011.03.031

Borja, A., Prins, T., Simboura, N., Andersen, J. H., Berg, T., Marques, J. C., et al. (2014). Tales from a thousand and one ways to integrate marine ecosystem components when assessing the environmental status. *Front. Mar. Sci.* 1:72. doi: 10.3389/fmars.2014.00072
Boudouresque, C., Bernard, G., Pergent, G., Pergent, G., Shili, A., and Verlaque, M. (2009). Regression of Mediterranean seagrasses caused by natural processes and anthropogenic disturbances and stress: a critical review. *Bot. Mar.* 52, 395-418. doi: 10.1515/BOT.2009.057

Cardillo, M., Mace, G. M., Gittleman, J. L., Jones, K. E., Bielby, J., and Purvis, A. (2008). The predictability of extinction: biological and external correlates of decline in mammals. *Proc. Biol. Sci.* 275, 1441-1448. doi: 10.1098/rspb.2008.0179

Carstensen, J., and Lindegarth, M. (2016). Confidence in ecological indicators: A framework for quantifying uncertainty components from monitoring data. *Ecol. Indic.* 67, 306-317. doi: 10.1016/j.ecolind.2016.03.002 1470-160X

Certain, G., Skarpaas, O., Bjerke, J. W., Framstad, E., Lindholm, M., Nilsen, J. E., et al. (2011). The nature index: a general framework for synthesizing knowledge on the state of biodiversity. *PLoS ONE* 6:e18930. doi: 10.1371/journal.pone.0018930

Claudet, J., and Fraschetti, S. (2010). Human-driven impacts on marine habitats: a regional meta-analysis in the Mediterranean Sea. *Biol. Conserv.* 143, 2195-2206. doi: 10.1016/j.biocon.2010.06.004

Cochrane, S. K. J., Connor, D. W., Nilsson, P., Mitchell, I., Reker, J., Franco, J., et al. (2010). *Marine Strategy Framework Directive – Task Group 1 Report Biological Diversity. EUR 24337 EN-Joint Research Centre.* Office for Official Publications of the European Communities, 110.

Coll, M., Piroddi, C., Albouy, C., Lasram, F. B. R., Cheung, W. W. L., Christensen, V., et al. (2012). The Mediterranean Sea under siege: spatial overlap between marine biodiversity, cumulative threats and marine reserves. *Glob. Ecol. Biogeogr.* 21, 465-480. doi: 10.1111/gb.1466-8238.2011.00697.x

Coll, M., Piroddi, C., Steenbeek, J., Kaschner, K., Ben Rais Lasram, F., Aguzzi, J., et al. (2010). The biodiversity of the Mediterranean Sea: estimates, patterns, and threats. *PLoS ONE* 5:e11842. doi: 10.1371/journal.pone.0011842

Coll, M., Shannon, L. J., Kleisner, K. M., Juan-jordá, M. J., Bundy, A., Akoglu, A. G., et al. (2015). Ecological indicators to capture the effects of fishing on biodiversity and conservation status of marine ecosystems. *Ecol. Indic.* 60, 947-962. doi: 10.1016/j.ecolind.2015.08.048

Croxall, J. P., Butchart, S. H. M., Lascelles, B., Stattersfield, A. J., Sullivan, B., Symes, A., et al. (2012). Seabird conservation status, threats and priority actions: a global assessment. *Bird Conserv. Int.* 22, 1-34. doi: 10.1017/S0959270912000020

Danovaro, R., Fonda Umani, S., and Pusceddu, A. (2009). Climate change and the potential spreading of marine mucilage and microbial pathogens in the Mediterranean Sea. *PLoS ONE* 4:e7006. doi: 10.1371/journal.pone. 0007006

Degobbis, D., Precali, R., Ivancic, I., Smodlaka, N., Fuks, D., and Kveder, S. (2000). Long-term Changes in the northern Adriatic ecosystem related to anthropogenic eutrophication. *Int. J. Environ. Pollut.* 13, 495-533. doi: 10.1504/IJEP.2000.002332

Dencheva, K., and Doncheva, V. (2014). "Ecological Index (EI) - tool for estimation of ecological status in coastal and transitional waters in compliance with European Water Framework Directive," in *Proceedings of Twelfth International Conference on Marine Sciences and technologies (Varna)*, 219-226.

Duarte, C. M. (2000). Marine biodiversity and ecosystem services: an elusive link. *J. Exp. Mar. Biol. Ecol.* 250, 117-131. doi: 10.1016/S0022-0981(00)00194-5

Dupont, C., Belin, A., Moreira, G., and Vermonden, B. (2014). *Article 12 - Technical Assessment of the MSFD 2012 Obligations*. National Report, Portugal Milieu Ltd; Consortium. Brussels, 35.

European Commission (2014). *Annex Accompanying the Document Commission Report to the Council and the European Parliament The First Phase of Implementation of the Marine Strategy Framework Directive (2008/56/EC)-The European Commission's Assessment And Guidance*. Brussels.

European Union (2008). *Directive 2008/56/EC of the European Parliament and of the Council of 17 June 2008 Establishing a Framework for Community Action in the Field of Marine Environmental Policy (Marine Strategy Framework Directive)*. Available online at: http://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32008L0056&from=en

European Union (2010). *COMMISSION DECISION of 1 September 2010 on Criteria and Methodological Standards on Good Environmental Status of Marine Waters*. Official Journal of the European Union. Available online at: http://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2010: 232:0014:0024:EN:PDF

Feary, D. A., Fowler, A. M., and Ward, T. J. (2014). Developing a rapid method for undertaking the World Ocean Assessment in data-poor regions - a case study using the South China Sea Large Marine Ecosystem. *Ocean Coast. Manag.* 95, 129-137. doi: 10.1016/j.ocecoaman.2014.04.006

Folke, C., Carpenter, S., Walker, B., Scheffer, M., Elmqvist, T., Gunderson, L., et al. (2004). Regime shifts, resilience, and biodiversity in ecosystem management. *Annu. Rev. Ecol. Evol. Syst.* 35, 557-581. doi: 10.1146/annurev.ecolsys.35. 021103.105711

Frantzis, A. (2009). *Cetaceans in Greece: Present Status of Knowledge*. Initiative for the Conservation of Cetaceans in Greece, Athens, 94.

Furness, R. W., and Camphuysen, C. J. (1997). Seabirds as monitors of the marine environment. *ICES J. Mar. Sci.* 54, 726-737.

Halpern, B. S., Longo, C., Hardy, D., McLeod, K. L., Samhoury, J. F., Katona, S. K., et al. (2012). An index to assess the health and benefits of the global ocean. *Nature* 488, 615-620. doi: 10.1038/nature11397

Halpern, B. S., Longo, C., Scarborough, C., Hardy, D., Best, B. D., Doney, S. C., et al. (2014). Assessing the health of the U.S. West Coast with a regional-scale application of the Ocean health index. *PLoS ONE* 9:e98995. doi: 10.1371/journal.pone.0098995

HELCOM (2009a). "Biodiversity in the Baltic Sea - An integrated thematic assessment on biodiversity and nature conservation in the Baltic Sea," in *Baltic Sea Environment Proceeding No. 116B*. Available online at: http://helcom.fi/ Lists/Publications/BSEP116B.pdf

HELCOM (2009b). "Eutrophication in the Baltic Sea - an integrated thematic assessment of the effects of nutrient enrichment and eutrophication in the Baltic Sea region Balt," in *Baltic Sea Environment Proceeding No. 115B*. Available online at: http://www.helcom.fi/Lists/Publications/BSEP115b.pdf

HELCOM (2010). "Ecosystem health of the Baltic Sea 2003-2007: HELCOM Initial Holistic Assessment," in *Baltic Sea Environment Proceeding No. 122*. Available online at: http://www.helcom.fi/lists/publications/bsep122.pdf.

Herman, P., Peter Beauchard, O., and van Duren, L. (2014). *De Staat van de Noordzee*. Available online at: http://www.noordzeedagen.nl/web/file?uuid=c49a717f-f648-4ef5-a8e6-8154ccac1b9d&owner=e80c7a5-1430d-4d1d-b34d-c4fd6e0f278d

Hummel, H., Frost, M., Juanes, J. A., Kochmann, J., Castellanos Perez Bolde, C. F., Aneiros, F., et al. (2015). A comparison of the degree of implementation of marine biodiversity indicators by European countries in relation to the Marine Strategy Framework Directive (MSFD). *J.Mar. Biol. Assoc. U.K.* 95, 1519-1531. doi: 10.1017/S0025315415000235

Jørgensen, L. L., Planque, B., Thangstad, T. H., and Certain, G. (2016). Vulnerability of megabenthic species to trawling in the Barents Sea. *ICES J. Mar. Sci.* 73(Suppl. 1), i84-i97. doi: 10.1093/icesjms/fsv107

Karlson, K., Rosenberg, R., and Bonsdorff, E. (2002). Temporal and spatial large-scale effects of eutrophication and oxygen deficiency on benthic fauna in Scandinavian and Baltic waters - a review. *Oceanogr. Mar. Biol.* 40, 427-489.

Katsanevakis, S., Zenetos, A., Poursanidis, D., Nunes, A. L., Deriu, I., Bogucarskis, K., et al. (2013). ELNAIS meets EASIN: distribution of marine alien species in Greece using EASIN mapping services and ELNAIS spatial data. *Mediterr. Mar. Sci.* 14, 95-98. doi: 10.12681/mms.329

Kleisner, K. M., Coll, M., Lynam, C. P., Bundy, A., Shannon, L., Shin, Y. J., et al. (2014). Evaluating changes in marine communities that provide ecosystem services through comparative assessments of community indicators. *Ecosyst. Serv.* 16, 413-429. doi: 10.1016/j.ecoser.2015.02.002

Kovacs, K. M., Lydersen, C., Overland, J. E., and Moore, S. E. (2011). Impacts of changing sea-ice conditions on Arctic marine mammals. *Mar. Biodivers.* 41, 181-194. doi: 10.1007/s12526-010-0061-0

MAMAOT (2012). *Estratégia Marinha para a Subdivisão do Continente*. Diretiva Quadro Estratégia Marinha. Ministério da Agricultura, do Mar, do Ambiente e do Ordenamento do Território. Outubro de 2012.

Micheli, F., Halpern, B. S., Walbridge, S., Ciriaco, S., Ferretti, F., Fraschetti, S., et al. (2013a). Cumulative human impacts on Mediterranean and Black Sea marine ecosystems: assessing current pressures and opportunities. *PLoS ONE* 8:e79889. doi: 10.1371/journal.pone.0079889

Micheli, F., Levin, N., Giakoumi, S., Katsanevakis, S., Abdulla, A., Coll, M., et al. (2013b). Setting priorities for regional conservation planning in the Mediterranean Sea. *PLoS ONE* 8:e59038. doi: 10.1371/journal.pone.0059038

Moncheva, S., Stefanova, K., Doncheva, V., Hristova, O., Dzhurova, B., and Racheva, E. (2015). "Plankton indicators to inform eutrophication
management,” in *Proceedings of the Twelfth International Conference on the Mediterranean Coastal Environment MEDCOAST’15*, ed E. Ozhan (Varna; Dalyan: Mediterranean Coastal Foundation), 351–362.

National Initial Assessment and GES Report (2013). Available online at: http://bsbd.org/UserFiles/File/Initial%20Assessment.pdf

National Report on the State and Protection of the Environment in Bulgaria (2014). State of Marine Waters. Sofia: Ministry of Environment and Water; Executive Environmental Agency.

Naturstyrelsen (2012). Danmarks Havstrategi. Basisanalyse, 100 (In Danish).

Norderhaug, K. M., and Christie, H. C. (2009). Sea urchin grazing and kelp re-vegetation in the NE Atlantic. *Mar. Biol. Res.* 5, 515–528. doi: 10.1080/17451000902932985

Ojaveer, H., Jaanus, A., MacKenzie, B. R., Martin, G., Olenin, S., Radziejewska, T., et al. (2010). Status of biodiversity in the Baltic Sea. *PLoS ONE* 5:e12467. doi: 10.1371/journal.pone.0012467

Österblom, H., Hansson, S., Larsson, U., Hjerne, O., Wulff, F., and Elmgren, R., et al. (2007). Human-induced trophic cascades and ecological regime shifts in the Baltic Sea. *Ecosystems* 10, 877–889. doi: 10.1007/s10021-007-9069-0

Palialexis, A., Tornero, V., Barbone, E., Gonzalez, D., Hanke, G., Cardoso, A. C., et al. (2014). *In-depth Assessment of the EU Member States' Submissions for the Marine Framework Strategy Framework Directive under Articles 8, 9 and 10 JRC Scientific and Technical Reports*. Publications Office of the European Union. Available online at: http://publications.jrc.ec.europa.eu/repository/bitstream/111111111/30749/1/lbna26473enn.pdf

Papaconstantinou, C. (2014). *Fauna Graeciae. An Updated Checklist of the Fishes in the Hellenic Seas, Monographs on Marine Sciences*, Vol. 7, Athens, 340.

Parkinson, C. L., Cavalieri, D. J., Gloersen, P., Zwally, H. J., and Comiso, J. C. (1999). Arctic sea ice extent, areas, and trends, 1978–1996. *J. Geophys. Res.* 104, 20837–20856. doi: 10.1029/1999JC900082

Prins, T. C., Borja, A., Simboura, N., Tsangaris, C., Van der Meulen, M. D., Boon, A. R., et al. (2013). *Coherent Geographic Scales and Aggregation Rules for Environmental Status Assessment within the Marine Strategy Framework Directive*. Towards a draft guidance. Deltares/AZTI/HCMR, Report 1207879-000-ZKS-0014 to the European Commission (Delft), 56.

Probst, W. N., and Lynam, C. P. (2016). Integrated assessment results depend on aggregation method and framework structure - A case study within the European marine strategy framework directive. *Ecol. Indic.* 61, 871–881. doi: 10.1016/j.ecolind.2015.10.040

Pusceddu, A., Bianchelli, S., Martin, J., Puig, P., Palanques, A., Masqué, P., et al. (2014). Chronic and intensive bottom trawling impairs deep-sea biodiversity and ecosystem functioning. *Proc. Natl. Acad. Sci. U.S.A.* 111, 8861–8866. doi: 10.1073/pnas.1405454111

Rinde, E., Christie, H., Fagerli, C. W., Bekkby, T., Gundersen, H., Norderhaug, K. M., et al. (2014). The influence of physical factors on kelp and sea urchin distribution in previously and still grazed areas in the NE Atlantic. *PLoS ONE* 9:e100222. doi: 10.1371/journal.pone.0100222

Schipper, J., Chanson, J. S., Chiozza, F., Cox, N. A., Hoffmann, M., Katariya, V., et al. (2008). The status of the world's land and marine mammals: diversity, threat, and knowledge. *Science* 322, 225–230. doi: 10.1126/science.1165115

Selig, E. R., Longo, C., Halpern, B. S., Best, B. D., Hardy, D., Elfes, C. T., et al. (2013). Assessing global marine biodiversity status within a coupled socio-ecological perspective. *PLoS ONE* 8:e60284. doi: 10.1371/journal.pone.0060284

Shtereva, G., Hristova, O., Dzhurova, B., and Krastev, A. (2012). *Hydrochemistry of Varna Bay (2000–2008)*, JEPE 13, No 3A, 1695–1702. Available online at: http://www.jepe-journal.info/journal-content/vol13-no-3a

Simboura, N., Pavlidou, A., Bald, J., Tsapakis, M., Pagou, P., Ch. Zeri, A., et al. (2016). Response of ecological indices to nutrient and chemical contaminant

stress factors in Eastern Mediterranean coastal waters. *Ecol. Indic.* 70, 89–105. doi: 10.1016/j.ecolind.2016.05.018

Simboura, N., Tsapakis, M., Pavlidou, A., Assimakopoulou, G., Pagou, K., Kontoyiannis, H., et al. (2015). Assessment of the environmental status in Hellenic coastal waters (Eastern Mediterranean): from the water framework directive to the marine strategy water framework directive. *Mediterr. Mar. Sci.* 16, 46–64. doi: 10.12681/mms.960

Simboura, N., Zenetos, A., and Pancucci-Papadopoulou, M. A. (2014). Benthic community indicators over a long period of monitoring (2000–2012) of the Saronikos Gulf, Greece, Eastern Mediterranean. *Environ. Monit. Assess.* 186, 3809–3821. doi: 10.1007/s10661-014-3659-z

Sivertsen, K. (2006). Overgrazing of kelp beds along the coast of Norway. *J. Appl. Phycol.* 18, 599–610. doi: 10.1007/s10811-006-9064-4

SoHelFI (2007). *State of Hellenic Fisheries*. Athens: HCMR Publication.

SoHelME (2005). *State of the Hellenic Marine Environment*. HCMR Publications.

Sorteberg, A., and Kvingedal, B. (2006). Atmospheric forcing on the Barents Sea winter ice extent. *J. Clim.* 19, 4772–4784. doi: 10.1175/JCLI 3885.1

Sunnana, K., Fossheim, M., and og Olseng, C.D. (red.) (2010). *Forvaltningsplan Barentshavet – Rapport fra Overvakingsgruppen 2010*. Fisken og havet, særn, Norway: Institute of Marine Research, 1b-2010.

Taylor, B. L., Martinez, M., Gerrodette, T., Barlow, J., and Hrovat, Y. N. (2007). Lessons from monitoring trends in abundance of marine mammals. *Mar. Mamm. Sci.* 23, 157–175. doi: 10.1111/j.1748-7692.2006.00092.x

UNEP (The United Nations Convention on Biological Diversity) (1992). Available online at: https://www.cbd.int/doc/legal/cbd-en.pdf

UNEP (2015). "United Nation Environment Program Mediterranean Action plan 2015," *Twelfth Meeting of Focal Points for Specially Protected Areas* (Athens).

UNEP/MAP-RAC/SPA (2015). *Adriatic Sea: Description of the Ecology and Identification of the Areas That May Deserve to be Protected*, Edited by C. Cerrano. Tunis: RAC/SPA, 92.

Vasilakopoulos, P., Maravelias, C. D., Tserpes, G. (2015). The alarming decline of mediterranean fish stocks *Curr. Biol.* 24, 1643–1648. doi: 10.1016/j.cub.2014.05.070

Wassmann, P., Reigstad, M., Haug, T., Rudels, B., Carroll, M. L., Hop, H., et al. (2006). Food webs and carbon flux in the Barents Sea. *Prog. Oceanogr.* 71, 232–287. doi: 10.1016/j.pocean.2006.10.003

Worm, B., Barbier, E. B., Beaumont, N., Duffy, J. E., Folke, C., Halpern, B. S., et al. (2006). Impacts of biodiversity loss on ocean ecosystem services. *Science* 314, 787–790. doi: 10.1126/science.1132294

ZenetoΣ, A., Akel, E. H. KH., Apostolidis, C., Bilecenoglu, M., Bitar, G., Buchet, V., et al. (2015). New Mediterranean biodiversity records (April 2015). *Mediterr. Mar. Sci.* 16, 266–284. doi: 10.12681/mms.1292

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2016 Uusitalo, Blanchet, Andersen, Beauchard, Berg, Bianchelli, Cantafaro, Carstensen, Carugati, Cochrane, Danovaro, Heiskanen, Karvinen, Moncheva, Murray, Neto, Nygård, Pantazi, Papadopoulou, Simboura, Srèbalienè, Uyarra and Borja. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.