Hostname: page-component-7b9c58cd5d-nzzs5 Total loading time: 0 Render date: 2025-03-15T19:31:03.375Z Has data issue: false hasContentIssue false

Essence of the patterns of cover and richness of intertidal hard bottom communities: a pan-European study

Published online by Cambridge University Press:  04 October 2016

Jonne Kotta*
Affiliation:
Estonian Marine Institute, University of Tartu, Tallinn, Estonia
Helen Orav-Kotta
Affiliation:
Estonian Marine Institute, University of Tartu, Tallinn, Estonia
Holger Jänes
Affiliation:
Estonian Marine Institute, University of Tartu, Tallinn, Estonia
Herman Hummel
Affiliation:
Monitor Taskforce, Royal Netherlands Institute for Sea Research (NIOZ), Yerseke, the Netherlands
Christos Arvanitidis
Affiliation:
Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research, Heraklion, Crete, Greece
Pim Van Avesaath
Affiliation:
Monitor Taskforce, Royal Netherlands Institute for Sea Research (NIOZ), Yerseke, the Netherlands
Guy Bachelet
Affiliation:
Arcachon Marine Station, CNRS, Université de Bordeaux, EPOC, Arcachon, France
Lisandro Benedetti-Cecchi
Affiliation:
Department of Biology, University of Pisa, Pisa, Italy
Natalia Bojanić
Affiliation:
Institute of Oceanography and Fisheries, Split, Croatia
Serena Como
Affiliation:
CNR, Institute for Coastal Marine Environment, Torregrande, Oristano, Italy
Stefania Coppa
Affiliation:
CNR, Institute for Coastal Marine Environment, Torregrande, Oristano, Italy
Jennifer Coughlan
Affiliation:
School of Biology and Environmental Science and Earth Institute, University College Dublin, Ireland
Tasman Crowe
Affiliation:
School of Biology and Environmental Science and Earth Institute, University College Dublin, Ireland
Martina Dal Bello
Affiliation:
Department of Biology, University of Pisa, Pisa, Italy
Steven Degraer
Affiliation:
Royal Belgian Institute of Natural Sciences, Operational Directorate Natural Environment, Aquatic and Terrestrial Ecology, Marine Ecology and Management, Brussels and Oostende, Belgium
Jose Antonio Juanes De La Pena
Affiliation:
Environmental Hydraulics Institute, Universidad de Cantabria, Santander, Spain
Valentina Kirienko Fernandes De Matos
Affiliation:
MARE – Marine and Environmental Sciences Centre and IMAR – Institute of Marine Research, University of the Azores, 9901-862 Horta, Portugal
Free Espinosa
Affiliation:
Laboratorio de Biología Marina, Universidad de Sevilla, Sevilla, Spain
Sarah Faulwetter
Affiliation:
Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research, Heraklion, Crete, Greece
Matt Frost
Affiliation:
Marine Biological Association, Plymouth, UK
Xabier Guinda
Affiliation:
Environmental Hydraulics Institute, Universidad de Cantabria, Santander, Spain
Emilia Jankowska
Affiliation:
Institute of Oceanology, Polish Academy of Sciences, Sopot, Poland
Jérôme Jourde
Affiliation:
Observatoire de la biodiversité (OBIONE), UMR 7266 LIttoral ENvironnement et Sociétés, CNRS/University of La Rochelle, France
Francis Kerckhof
Affiliation:
Royal Belgian Institute of Natural Sciences, Operational Directorate Natural Environment, Aquatic and Terrestrial Ecology, Marine Ecology and Management, Brussels and Oostende, Belgium
Nicolas Lavesque
Affiliation:
Arcachon Marine Station, CNRS, Université de Bordeaux, EPOC, Arcachon, France
Jean-Charles Leclerc
Affiliation:
Sorbonne Universités, UPMC Univ. Paris 6, Station Biologique, Place Georges Teissier, 29680 Roscoff, France CNRS, UMR 7144, Station Biologique, Place Georges Teissier, 29680 Roscoff, France
Paolo Magni
Affiliation:
CNR, Institute for Coastal Marine Environment, Torregrande, Oristano, Italy
Christina Pavloudi
Affiliation:
Institute of Marine Biology, Biotechnology and Aquaculture, Hellenic Centre for Marine Research, Heraklion, Crete, Greece
Maria Luiza Pedrotti
Affiliation:
Sorbonne Universités, UPMC Univ. Paris 6, UMR 7093, LOV, Villefranche-sur-mer, France
Ohad Peleg
Affiliation:
National Institute of Oceanography, Israel Oceanographic and Limnological Research, Haifa, Israel
Angel Pérez-Ruzafa
Affiliation:
Department of Ecology and Hydrology, Regional Campus of International Excellence ‘Campus Mare Nostrum’, University of Murcia, Spain
Araceli Puente
Affiliation:
Environmental Hydraulics Institute, Universidad de Cantabria, Santander, Spain
Pedro Ribeiro
Affiliation:
MARE – Marine and Environmental Sciences Centre and IMAR – Institute of Marine Research, University of the Azores, 9901-862 Horta, Portugal
Gil Rilov
Affiliation:
National Institute of Oceanography, Israel Oceanographic and Limnological Research, Haifa, Israel
Maria Rousou
Affiliation:
Marine and Environmental Research Lab Ltd, Nicosia, Cyprus
Tomas Ruginis
Affiliation:
Open Access Centre for Marine Research, Klaipeda University, Lithuania
Teresa Silva
Affiliation:
Laboratório de Ciências do Mar, Marine and Environmental Sciences Centre (MARE), Universidade de Évora, Sines, Portugal
Nathalie Simon
Affiliation:
Sorbonne Universités, UPMC Univ. Paris 6, Station Biologique, Place Georges Teissier, 29680 Roscoff, France CNRS, UMR 7144, Station Biologique, Place Georges Teissier, 29680 Roscoff, France
Isabel Sousa-Pinto
Affiliation:
Centre for Marine and Environmental Research, CIIMAR, Porto, Portugal
Jesús Troncoso
Affiliation:
ECIMAT, Station of Marine Sciences of Toralla, Department of Ecology and Animal Biology, University of Vigo, Spain
Jan Warzocha
Affiliation:
Marine Fisheries Institute, Gdynia, Poland
Jan Marcin Weslawski
Affiliation:
Institute of Oceanology, Polish Academy of Sciences, Sopot, Poland
*
Correspondence should be addressed to: J. Kotta, Estonian Marine Institute, University of Tartu, Tallinn, Estonia Email: jonne@sea.ee
Rights & Permissions [Opens in a new window]

Abstract

Coastal ecosystems are highly complex and driven by multiple environmental factors. To date we lack scientific evidence for the relative contribution of natural and anthropogenic drivers for the majority of marine habitats in order to adequately assess the role of different stressors across the European seas. Such relationship can be investigated by analysing the correlation between environmental variables and biotic patterns in multivariate space and taking into account non-linearities. Within the framework of the EMBOS (European Marine Biodiversity Observatory System) programme, hard bottom intertidal communities were sampled in a standardized way across European seas. Links between key natural and anthropogenic drivers and hard bottom communities were analysed using Boosted Regression Trees modelling. The study identified strong interregional variability and showed that patterns of hard bottom macroalgal and invertebrate communities were primarily a function of tidal regime, nutrient loading and water temperature (anomalies). The strength and shape of functional form relationships varied widely however among types of organisms (understorey algae composing mostly filamentous species, canopy-forming algae or sessile invertebrates) and aggregated community variables (cover or richness). Tidal regime significantly modulated the effect of nutrient load on the cover and richness of understorey algae and sessile invertebrates. In contrast, hydroclimate was more important for canopy algae and temperature anomalies and hydroclimate separately or interactively contributed to the observed patterns. The analyses also suggested that climate-induced shifts in weather patterns may result in the loss of algal richness and thereby in the loss of functional diversity in European hard bottom intertidal areas.

Type
Research Article
Copyright
Copyright © Marine Biological Association of the United Kingdom 2016 

INTRODUCTION

As the demand on terrestrial resources (energy, food and space) has intensified (e.g. Hubbert, Reference Hubbert1949; Hoogwijk et al., Reference Hoogwijk, Faaij, van den Broek, Berndes, Gielen and Turkenburg2003), there has been an increasing interest in utilizing marine resources. This is because oceans and seas cover over 70% of the Earth's surface and host an extraordinarily rich biodiversity. Nearshore areas in particular are among the most important biomes with all major taxonomic groups being represented and responsible for 90% of the world's marine primary production (Kaiser et al., Reference Kaiser, Attrill, Jennings, Thomas, Barnes, Brierley, Hiddink, Kaartokallio, Polunin and Raffaelli2011). The consequent increase and diversification in use of coastal natural resources has resulted in a necessity for stock-taking and possible protective measures for a variety of marine species (Halpern et al., Reference Halpern, Walbridge, Selkoe, Kappel, Micheli, D'Agrosa, Bruno, Casey, Ebert, Fox, Fujita, Heinemann, Lenihan, Madin, Perry, Selig, Spalding, Steneck and Watson2008). Obtaining information on the distribution of marine species has become an important goal; nevertheless, large-scale patterns of biodiversity and factors causing change in biodiversity are still not well understood, mostly due to the inconsistency of data sources and differences in methodology (e.g. Goulletquer et al., Reference Goulletquer, Gros, Boeuf and Weber2014).

Marine biodiversity varies over large scales whereas most studies focus on local and partly on regional scales. As direct mapping of large-scale patterns is extremely costly, knowledge-based management across large spatial scales requires alternative, low-cost sampling models as well as research strategies beyond the scope and vision of current national and institutional monitoring activities. There is also a requirement for a standard way of measuring local community patterns over large-scale environmental gradients. The recently built network of research locations in Europe (EMBOS) has addressed this issue by establishing a cost-effective monitoring programme for the diversity of hard bottom intertidal communities at a European scale. The principal aim of EMBOS is to assess broad-scale patterns in marine biodiversity in relation to natural and anthropogenic gradients. The assumption is that non-random patterns of variability in richness and cover of species in communities probably reflect the scale of underlining processes. Measuring variability at several spatial scales therefore enables the identification of the range of processes likely to be behind the observed differences among communities.

The observed patterns of each eco-region are not just defined by large-scale biogeographic processes, however (e.g. Roy et al., Reference Roy, Jablonski, Valentine and Rosenberg1998), but are also (and often interactively) a product of regional (Ramos et al., Reference Ramos, Puente and Juanes2016b) and local-scale variability in biogeochemical conditions (Dal Bello et al., Reference Dal Bello, Leclerc, Benedetti-Cecchi, Arvanitidis, van Avesaath, Bachelet, Bojanić, Como, Coppa, Coughlan, Crowe, Degraer, Espinosa, Faulwetter, Frost, Guinda, Ikauniece, Jankowska, Jourde, Kerckhof, Kotta, Lavesque, de Lucia, Magni, Fernandes de Matos, Orav-Kotta, Pavloudi, Pedrotti, Peleg, de la Pena, Puente, Ribeiro, Rigaut-Jalabert, Rilov, Rousou, Rubal, Ruginis, Pérez-Ruzafa, Silva, Simon, Sousa-Pinto, Troncoso, Warzocha, Weslawski and Hummel2016; Ramos et al., Reference Ramos, Puente, Guinda and Juanes2016a). The EMBOS programme has allowed a thorough examination of scale (see Dal Bello et al., Reference Dal Bello, Leclerc, Benedetti-Cecchi, Arvanitidis, van Avesaath, Bachelet, Bojanić, Como, Coppa, Coughlan, Crowe, Degraer, Espinosa, Faulwetter, Frost, Guinda, Ikauniece, Jankowska, Jourde, Kerckhof, Kotta, Lavesque, de Lucia, Magni, Fernandes de Matos, Orav-Kotta, Pavloudi, Pedrotti, Peleg, de la Pena, Puente, Ribeiro, Rigaut-Jalabert, Rilov, Rousou, Rubal, Ruginis, Pérez-Ruzafa, Silva, Simon, Sousa-Pinto, Troncoso, Warzocha, Weslawski and Hummel2016) but is also designed to measure the key environmental gradients at multiple spatial scales. The current study design matches with theoretical expectations that (1) local biological assemblages are defined by several interacting environmental gradients (Whittaker, Reference Whittaker1974); (2) local benthic diversity represents the interplay of historical evolutionary processes (i.e. the number of species in the regional species pool) combined with separate and interactive effects of local environmental gradients (e.g. Ricklefs & Schluter, Reference Ricklefs, Schluter, Ricklefs and Schluter1993; Holyoak et al., Reference Holyoak, Leibold, Moquet, Holt, Hoopes, Holyoak, Leibold, Moquet and Holt2005; Ramos et al., Reference Ramos, Puente, Guinda and Juanes2016a); and (3) these relationships are scale dependent i.e. local communities assemble from the macro- to the microscale and within these constraints, environmental stress plays a central role in shaping community diversity (Menge & Sutherland, Reference Menge and Sutherland1987).

In general, three main local gradients/stresses have been recognized in intertidal hard bottom habitats: wetness, exposure to wave action and salinity (Kaiser et al., Reference Kaiser, Attrill, Jennings, Thomas, Barnes, Brierley, Hiddink, Kaartokallio, Polunin and Raffaelli2011). When moving away from the waterline toward the highest part of the shore the environment becomes progressively drier. However, this gradient is amplified by waves and tides. Specifically, gently sloping areas tend to dissipate the wave energy arriving at the shore while steeper sloping shores experience much greater physical impact. In addition the wetness gradient is expected to be more pronounced on macrotidal shores (North Atlantic) compared with microtidal shores (Mediterranean and Baltic Seas). Nevertheless, features of rocky shore zonation patterns are similar worldwide independent of biogeographic region suggesting that similar abiotic forces (e.g. a combination of wave and tide effect) can generate similar biotic patterns (Stephenson & Stephenson, Reference Stephenson and Stephenson1972; Bird et al., Reference Bird, Franklin, Smith and Toonen2013 but see also Underwood, Reference Underwood1978). The salinity gradient occurs in estuaries and/or semi-enclosed seas and is often characterized by a significant loss of marine taxa and/or lowered diversity along a gradient of salinity reduction (e.g. Bonsdorff & Pearson, Reference Bonsdorff and Pearson1999).

In addition to natural environmental gradients, anthropogenic impacts are becoming increasingly important. To date, human activities have globally been linked to a diverse range of ecological changes in marine ecosystems including profound changes in coastal biodiversity (Goulletquer et al., Reference Goulletquer, Gros, Boeuf and Weber2014). Among the most commonly reported effects, climate change is expected to induce shifts in air and seawater temperature and affect benthic communities globally but often these effects are very context-specific and depend on the local environmental setting (Lima et al., Reference Lima, Ribeiro, Queiroz, Hawkins and Santos2007a, Reference Lima, Ribeiro, Queiroz, Xavier, Tarroso, Hawkins and Santosb, Kotta et al., Reference Kotta, Möller, Orav-Kotta and Pärnoja2014; Thorner et al., Reference Thorner, Kumar and Smith2014). In addition eutrophication is considered to be a significant concern in some regional seas (e.g. the Baltic, North Sea) (Conley et al., Reference Conley, Carstensen, Aigars, Axe, Bonsdorff, Eremina, Haahti, Humborg, Jonsson, Kotta, Lännegren, Larsson, Maximov, Rodriguez Medina, Lysiak-Pastuszak, Remeikaité-Nikiené, Walve, Wilhelms and Zillén2011) but locally also in the Atlantic Ocean and Mediterranean Sea (Howarth et al., Reference Howarth, Billen, Swaney, Townsend, Jaworski, Lajtha, Downing, Elmgren, Caraco, Jordan, Berendse, Freney, Kudeyarov, Murdoch and Zhu1996; Karydis & Kitsiou, Reference Karydis and Kitsiou2012). In order to assess the relative ecological importance of human impacts, it is necessary to understand how marine communities respond to various natural and anthropogenic ecological disturbances either separately or as interacting drivers of change. Such knowledge is of the utmost importance for understanding biodiversity patterns and assessing the consequences of various human activities on ecosystem integrity which in turn informs the management of natural resources.

Intertidal community structure is shaped by natural as well as by anthropogenic variables (e.g. Zacharias & Roff, Reference Zacharias and Roff2001; Crowe et al., Reference Crowe, Cusson, Bulleri, Davoult, Arenas, Aspden, Benedetti-Cecchi, Bevilacqua, Davidson, Defew, Fraschetti, Golléty, Griffin, Herkül, Kotta, Migné, Molis, Nicol, Noël, Sousa Pinto, Valdivia, Vaselli and Jenkins2013), however, only a handful of studies have been conducted that examine and compare the importance of these variables across distant aquatic ecosystems (Rees et al., Reference Rees, Pendle, Waldock, Limpenny and Boyd1999; Danovaro et al., Reference Danovaro, Canals, Gambi, Heussner, Lampadariou and Vanreusel2009; Narayanaswamy et al., Reference Narayanaswamy, Coll, Danovaro, Davidson, Ojaveer and Renaud2013; Schiele et al., Reference Schiele, Darr, Zettler, Berg, Blomqvist, Daunys, Jermakovs, Korpinen, Kotta, Nygård, von Weber, Voss and Warzocha2016). In the current study, we identified the role of multiple environmental variables on the patterns of cover and richness of intertidal hard-bottom communities throughout Europe and established the functional form relationships between these drivers and biotic patterns. The following hypotheses were put forward: (1) interregional differences account for a large part of the variability in hard bottom intertidal communities; (2) higher shore levels are characterized by reduced richness of macroalgae and invertebrates compared with lower levels; (3) elevated wave exposure and increased tidal range will result in systematically higher benthic richness compared with sheltered habitats; (4) elevated regional-scale loads of nutrients are expected to increase both cover and richness of understorey algae compared with cover and richness of canopy-forming algae; (5) temperature anomalies (higher than regular temperatures) are expected to reduce the growth of canopy forming algae and facilitate the growth of alternative understorey algae.

MATERIALS AND METHODS

Study area

Europe covers an extensive range of coastal areas spanning from subtropical to subpolar climatic regions. Such a wide climatic gradient leads to a great variety of environmental conditions impacting the diversity and stability of intertidal benthic macrophyte and invertebrate communities (e.g. Bulleri et al., Reference Bulleri, Benedetti-Cecchi, Cusson, Maggi, Arenas, Aspden, Bertocci, Crowe, Davoult, Eriksson, Fraschetti, Golléty, Griffin, Jenkins, Kotta, Kraufvelin, Molis, Sousa Pinto, Terlizzi, Valdivia and Paterson2012). A remarkable interregional variability exists: for example the semi-enclosed Mediterranean Sea with its oligotrophic conditions, high biodiversity and rates of endemism is very different to the Baltic Sea, which is characterized by eutrophic conditions and low biodiversity. Furthermore, the coasts of the North-East Atlantic Ocean, Arctic Ocean, Mediterranean, North and Baltic Seas differ greatly in terms of temperature, light availability, wave exposure and tidal regime.

Sampling and supporting environmental variables

The large-scale network of locations established through the implementation of EMBOS (http://embos.info/) was utilized for the study. In each study location the representative field sites were chosen based on commonalities and integrity of intertidal habitat as well as the lack of mobile sediment impacting bedrock. The network of locations covered all major European basins (Figure 1).

Fig. 1. Map of the sampling stations in the study area. The annual average SST isotherms (°C) are also indicated.

Sampling was carried out early in biological spring for each region in 2014 (i.e. in late March/early April in the Mediterranean, April in the temperate NE Atlantic, end of April to mid-May in the Baltic, and about one month after the loss of ice-cover in the Arctic). At each location two stations were sampled. Each station was 5–15 m wide and separated by a distance of 50–100 m. The stations were placed in representative areas of the shore where commonly occurring species were observed i.e. avoiding rockpools and other atypical features.

At each station (N = 39), the assessment was carried out at mid intertidal and lower intertidal levels. The mid intertidal was defined operationally as approximately 25% of the vertical extent of the shore centred on Mean Tidal Level. The lower intertidal was defined as 25% of the vertical extent of the shore upwards from Mean Spring Low Water. In most of the sites, however, high intertidal and upper subtidal were also sampled. In subtidal systems, stations were placed in areas dominated by macro-algae, rather than in urchin barren habitats.

During sampling the biological characterization involved the assessment of per cent coverage of different taxa of canopy-forming and understorey algae as well as sessile invertebrates within each of five quadrats placed haphazardly within a site at each shore height. The understorey algae mostly consisted of fleshy and filamentous species. In macrotidal shores (e.g. NE Atlantic) the quadrat size was 0.5 × 0.5 m (0.25 m2) and in microtidal systems (e.g. Mediterranean) this was 0.2 × 0.2 m (0.04 m2). The size of quadrats differed between the two regions to account for differences in tidal amplitude and size of organisms. The studied micro-tidal habitats are often very narrow and larger quadrats would not have allowed the differentiation between the low-shore and mid-shore habitats in this micro-tidal system. Smaller quadrats in the NE Atlantic would not have been appropriate to sample usually large-sized fucoid and kelp algae. The gathered taxonomic information was optimized for up-scaling, so that a minimum amount of information was lost when datasets were coupled and integrated. This improved the efficiency and performance of large-scale approaches and analyses.

Geographic coordinates (as latitude and longitude) were recorded during sampling. Tidal range was obtained from tide tables. Regional-scale proxies for weather, eutrophication and climate change variables were obtained from different online databases. A mean climatology for wind, solar energy and cloud cover was obtained from the Climatic Research Unit at the University of East Anglia (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1073). The climatology collection was compiled from existing data sources and algorithms, and was designed to satisfy the needs of modellers and investigators of the global carbon, water and energy cycle. The site data were interpolated as a function of latitude, longitude and elevation using thin-plate splines (Stackhouse & Gupta, Reference Stackhouse, Gupta, Hall, Collatz, Meeson, Los, De Colstoun and Landis2012).

Sea ice concentration (expressed as a mean percentage of ocean area covered by sea ice) was obtained from the Goddard Space Flight Center, Sea Ice Concentrations from Nimbus-7 Scanning Multichannel Microwave Radiometer and the Defense Meteorological Satellites Program Special Sensor Microwave/Imager passive microwave data (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=981). These original data were re-gridded by the National Snow and Ice Data Center from their original 25 km spatial resolution and EASE-Grid into equal angle Earth grids with quarter degree spatial resolutions in latitude/longitude (Armstrong & Knowles, Reference Armstrong, Knowles, Hall, Collatz, Meeson, Los, De Colstoun and Landis2010).

Nutrient pollution, sea surface temperatures and temperature anomalies data were obtained from Halpern et al. (Reference Halpern, Walbridge, Selkoe, Kappel, Micheli, D'Agrosa, Bruno, Casey, Ebert, Fox, Fujita, Heinemann, Lenihan, Madin, Perry, Selig, Spalding, Steneck and Watson2008). In their study nutrient pollution data came from Food and Agriculture Organization (FAO) national statistics (http://faostat.fao.org) on average annual use of fertilizers (nutrients) for the years 1993–2002 and were distributed across landscapes in agricultural lands with dasymetric techniques. Sea surface temperature (SST) data from 1985 to 2005 have been constructed from Advanced Very High Resolution Radiometer Pathfinder Version 5.0 by NOAA's National Oceanographic Data Center and the University of Miami's Rosenstiel School of Marine and Atmospheric Science (http://pathfinder.nodc.noaa.gov). Then it was assessed how many times the SST anomaly exceeded the standard deviation of SSTs for each 4 × 4 km grid cell and week of the year. By incorporating the standard deviation, this threshold-based approach accounts for natural variability at a given location, which can vary widely from place to place. Finally, a SST change metric (later referred to as temperature anomalies) was developed by subtracting the number of non-zero positive anomalies in the early period (1985–1990) from the number in the recent period (2000–2005). As the anthropogenic driver data (nutrient load and temperature anomalies) had extreme left-skewed distributions the data were log[X + 1]-transformed and then rescaled between 0–1, with the highest log-transformed value set = 1.

Boosted regression tree modelling

Environmental gradients are expected to affect biota interactively and when exploring potential influences on the diversity of intertidal hard bottom communities, the analytical framework should have capabilities to quantify the functional form of the relationships between separate environmental gradients and the biota, as well as their interactive effects. Modelling is the most appropriate tool, and several refined statistical approaches have already been applied in the field (Herkül et al., Reference Herkül, Kotta, Kutser and Vahtmäe2013; Sandman et al., Reference Sandman, Wikström, Blomqvist, Kautsky and Isaeus2013). The main focus of spatial modelling has often been directed towards predicting species distribution ranges (Peterson, Reference Peterson2003; Soberon & Peterson, Reference Soberon and Peterson2005) without investigating underlying mechanisms. Modelling can however improve our understanding of the relationship between the environment and biota and inform the theoretical understanding of biotic patterns and their driving forces (Kotta et al., Reference Kotta, Möller, Orav-Kotta and Pärnoja2014).

A novel predictive modelling technique, Boosted Regression Trees (BRT) does not use any predefined data model but instead uses an algorithm to analyse the relationship between the biota and environmental variables. BRT copes with different non-linear relationships which are common in ecological data but difficult to analyse using more traditional methods. The BRT also avoids overfitting the data, thereby providing robust estimates. Most importantly from an ecological perspective, BRT automatically detects and models interactive effects between predictors. Due to its strong predictive performance, BRT is increasingly used in ecology (Elith et al., Reference Elith, Leathwick and Hastie2008; Kotta et al., Reference Kotta, Kutser, Teeveer, Vahtmäe and Pärnoja2013).

Multicollinearity can be an issue with modelling when explaining relationships between environment and the biota. Thus, prior to modelling, a correlation analysis was conducted for environmental variables. Among the studied environmental variables only water temperature, solar radiation and cloudiness showed any correlation (r between 0.82 and 0.90, P < 0.01). Thus, in order to avoid multicollinearity issues, the final models included only water temperature.

The contribution of different environmental variables to the cover and richness of hard bottom intertidal macroalgae and invertebrates was explored using the BRT modelling. The BRT technique iteratively develops a large ensemble of small regression trees constructed from random subsets of the data. Each successive tree predicts the residuals from the previous tree to gradually boost the predictive performance of the overall model (Elith et al., Reference Elith, Leathwick and Hastie2008). In fitting a BRT the learning rate and the tree complexity must be specified. The learning rate determines the contribution of each successive tree to the final model, as it proceeds through the iterations. The selected value of a tree complexity determines whether only main effects (tree complexity = 1) or interactions are also included (tree complexity >1). Ultimately, the learning rate and tree complexity combined determine the total number of trees in the final model. Following the suggestions by Elith et al. (Reference Elith, Leathwick and Hastie2008) the model learning rate was kept at 0.1 and tree complexity at 5 for all models. It was also checked that the final models had more than 1000 trees. Nevertheless, a selection of model parameters had only marginal impact on model performance with optimal models improving predictions less than 1%. Model performance was evaluated using the cross validation statistics calculated during model fitting (Hastie et al., Reference Hastie, Tibshirani and Friedman2009). Thus, when running models, a random selection of 80% of the data was used for training the model and the rest of the data i.e. 20% was assigned for testing model accuracy. The BRT modelling was done in the statistical software R using the gbm package (RDC Team, 2013).

RESULTS

Outputs of the BRT models described a significant proportion of variability in the cover and richness of macroalgae and invertebrates inhabiting hard bottom intertidal habitats. For the majority of models the explained variance was higher than 80%. The model explained 71% of overall variability when only considering cover of canopy algae.

The results of the BRT models including region (i.e. basin) as an explanatory variable showed that interregional variability was very important for hard bottom intertidal macrophytes and invertebrates. There was a large difference however in the role of regional variability among types of organisms (understorey algae, canopy-forming algae or sessile invertebrates) and aggregated community variables (cover or richness). In general, regional variability was a better explanation for algal patterns rather than invertebrate patterns and total cover rather than richness. Regional variability explained 60% of the variability in canopy algae cover and 53% of the understorey algae cover but only 37% of the sessile invertebrate cover. In the richness models, the same values were about one-third lower than the levels in the cover models. Understorey algae had the highest cover in the English Channel, the Irish Sea and the Strait of Gibraltar and highest richness in the Bay of Biscay. Canopy algae had the highest cover in the Irish Sea, the Mediterranean Sea, the English Channel and the Baltic Sea and highest richness in the Bay of Biscay, the English Channel, the Mediterranean Sea and the Mid North Atlantic coasts. Both the cover and richness of sessile invertebrates were highly variable with no clear spatial patterns (Figure 2).

Fig. 2. Average values of the cover and richness of understorey algae, canopy algae and sessile invertebrates in different regions. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal.

As our main goal was to quantify roles of different environmental gradients on intertidal algae and invertebrates and regional variability largely incorporates changes in various natural and anthropogenic forcing, the BRT models below did not include the variable ‘region’ among the predictors. Such BRT models identified tidal level, tidal range, nutrient loading and water temperature (anomalies) to be the most significant predictors of both cover and richness of macroalgae and invertebrates (Table 1).

Table 1. Relative contribution of environmental predictors in the BRT models of total cover and richness of macroalgae and invertebrates.

The cover and richness of understorey algae was higher in subtidal and low intertidal areas compared with medium and high intertidal areas. Moreover, higher tidal range was associated with elevated values of cover and richness of understorey algae. Nutrient load was also among the most important variables predicting the understorey algal patterns. The cover of understorey algae increased with nutrient load whereas the richness of understorey algae was lowest at intermediate loads (Figure 3).

Fig. 3. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of understorey algae in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

As for the understorey algae, cover and richness of canopy algae was higher in subtidal and low intertidal areas compared with medium and high intertidal areas. Elevated tidal range reduced the cover of canopy algae but increased the associated taxonomic richness. Water temperature strongly contributed to the patterns of canopy algae with cover and richness being higher at reduced water temperatures. In contrast to understorey algae, the effect of nutrient load for patterns of canopy algae was moderate. Specifically, elevated nutrient loads were associated with higher cover of canopy algae whereas highest richness values were found at intermediate nutrient loads (Figure 4).

Fig. 4. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of canopy algae in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

The richness of canopy algae and sessile invertebrates showed similar patterns along the tidal levels and temperature gradient. However, no such similarities emerged for the cover of sessile invertebrates. Instead, the cover of sessile invertebrates was higher in high and medium intertidal areas compared with low intertidal and subtidal areas. Nutrient loads and tidal range highly contributed to the cover of sessile invertebrates with higher cover values expected at intermediate nutrient loads and tidal range (Figure 5).

Fig. 5. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of sessile invertebrates in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the section of methods for further information on environmental variables.

The BRT analyses also identified a few important interactive effects. Specifically, the effect of nutrient loads on understorey algae and sessile invertebrates varied widely among tidal levels. Although the cover of understorey algae increased with nutrient load across tidal levels, the effect was stronger in subtidal and low intertidal areas compared with medium and high intertidal areas. Similarly, the cover and richness of sessile invertebrates increased with nutrient loads in contrast to understorey algae where the stronger effects were found in high and medium intertidal areas compared with low intertidal and subtidal areas. The richness of understorey algae decreased with nutrient loads with the effect being proportionally stronger in subtidal and low intertidal areas compared with medium and high intertidal areas (Figure 6).

Fig. 6. Partial dependence plots of the BRT model for the cover and richness of understorey algae and sessile invertebrates at different tidal levels in the study area. The values of nutrient loads are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

Water temperature interacted with other environmental variables in the models of canopy algae. The cover of canopy algae decreased at elevated temperatures; however, this effect was much stronger in subtidal and low intertidal areas compared with medium and high intertidal areas. The richness of canopy algae was reduced by elevated ice scour. However, the effect disappeared in areas characterized by higher temperature anomalies (Figure 7).

Fig. 7. A partial dependence plot of the BRT model for the cover of canopy algae at different tidal levels and a three-dimensional partial dependence plot in the BRT model for the richness of canopy algae along the gradients of water temperature anomalies and ice cover in the study area. The values of temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

DISCUSSION

Based on the results of the analyses, initially proposed hypotheses were supported. The results confirmed that interregional variability accounted for about half of the variability in intertidal hard bottom macroalgal and invertebrate communities at the pan-European scale. Interregional variability displayed a systematically stronger effect on algal communities than on invertebrate communities and had a significantly stronger effect on cover patterns compared with richness of species. It was found that higher shore levels were characterized by reduced algal and invertebrate richness although elevated wave exposure and increased tidal range mitigated the effect. Elevated nutrient loads increased both cover and richness of understorey species but had a moderate effect on canopy algae and increases in seawater temperature and associated weather anomalies did not seem to have any effect on the cover of canopy algae whereas a very strong negative effect on canopy richness was observed. On the other hand, the patterns of understorey algae were essentially uncoupled with those of elevated temperature anomalies.

There exist broad differences in how biogeographic history affects local patterns of species cover and richness and these differences may explain why regional variability was more important for species cover than richness. Specifically, the number of species in the regional pool sets the maximum limit of species in a local community. At small spatial scales, however, local communities become quickly saturated by species and the number of species in the species pool is much higher than local communities can actually accommodate (Kotta & Witman, Reference Kotta, Witman and Wahl2009). Earlier studies have found clear evidence for saturation at the studied quadrat scale with the saturation effect being amplified towards the coastline characterized by high environmental stress (Russell et al., Reference Russell, Wood, Allison and Menge2006). Therefore, local species richness of the hard bottom intertidal habitats is defined either by limiting environmental factors or short-term biotic mechanisms that leads to the exclusion of some of the species from a community rather than regional speciation processes. As observed in this study, intertidal hard bottom habitats are often highly patchy, especially at the smallest spatial scales with different microhabitats being suitable for different sets of species. Such microhabitat heterogeneity also suggests the existence of a strong small-scale component of variability and a weak regional pool effects in intertidal species richness.

Total cover values, in contrast to species richness, are determined by environmental stressors such as desiccation, wave exposure, ice scour and nutrient loads that often play a central role in defining community structure (Menge & Sutherland, Reference Menge and Sutherland1987; Smale et al., Reference Smale, Burrows, Evans, King, Sayer, Yunnie and Moore2016). Such environmental stressors do not only operate at a regional scale, but can even introduce significant variability within the same region (Smale et al., Reference Smale, Burrows, Evans, King, Sayer, Yunnie and Moore2016) and thus synchronize cover values at multiple spatial scales, particularly where they are having a direct effect on recruitment or mortality (Dal Bello et al., Reference Dal Bello, Leclerc, Benedetti-Cecchi, Arvanitidis, van Avesaath, Bachelet, Bojanić, Como, Coppa, Coughlan, Crowe, Degraer, Espinosa, Faulwetter, Frost, Guinda, Ikauniece, Jankowska, Jourde, Kerckhof, Kotta, Lavesque, de Lucia, Magni, Fernandes de Matos, Orav-Kotta, Pavloudi, Pedrotti, Peleg, de la Pena, Puente, Ribeiro, Rigaut-Jalabert, Rilov, Rousou, Rubal, Ruginis, Pérez-Ruzafa, Silva, Simon, Sousa-Pinto, Troncoso, Warzocha, Weslawski and Hummel2016). For example, stressors such as wave exposure, ice scour and desiccation at upper intertidal zones have a notable impact on the abundance patterns of many sessile intertidal species (Dayton, Reference Dayton1971; Heaven & Scrosati, Reference Heaven and Scrosati2008). When the intensity or the severity of disturbances outweighs colonization, then the inhabited area remains smaller in size than potentially possible (Johnston et al., Reference Johnston, Keough and Qian2002). Regional variability in productivity and consumption may also have a fundamental influence on patterns of species cover. It has been shown that high productivity regions are more prone to overgrazing (functional types of grazers which cause barrens of diversity) than to nutrient loads, and low productivity regions are less affected by grazing (no barrens) and more affected by nutrients (Connell & Irving, Reference Connell and Irving2008).

We also expected that environmental gradients are major underlying components of the observed interregional variability. Apart from the potential effects of the regional species pool and species biogeography not resolved in this study, our BRT models showed that most variability in cover and richness of hard bottom macroalgae and invertebrates across Europe was explained by site-specific abiotic characteristics such as tidal regime and climate variables coupled with anthropogenic eutrophication.

More specifically, we expected higher shores to be characterized by reduced richness of macroalgae and invertebrates compared with lower shores. The BRT analyses confirmed this expectation as both cover and richness of understorey and canopy algae increased at lower intertidal and subtidal zones compared with mid- and high intertidal zones. Such findings confirm well-known principles of intertidal ecology whereby higher cover and richness of algal communities is found at lower intertidal and subtidal zones due to reduced risk of desiccation and thermal stress (Dayton, Reference Dayton1975) but in some cases probably due to reduced grazing or a combination of both (Underwood & Jernakoff, Reference Underwood and Jernakoff1981; Hawkins & Hartnoll, Reference Hawkins and Hartnoll1985). Species richness of sessile invertebrates also increased in line with the previously observed zonation patterns for algal communities which can be linked to physiological tolerances of species, foundation effects due to canopy-forming seaweeds and broader availability of suitable habitats as space, these being some of the most limiting factors for intertidal communities (Tomanek & Helmuth, Reference Tomanek and Helmuth2002). On the other hand, the cover of sessile invertebrates displayed the complete opposite i.e. cover increased at high and medium intertidal areas compared with low intertidal and subtidal areas. This most probably relates to the fact that competition and predation pressure decreases at higher intertidal zones allowing an increase in the cover of species that are able to effectively colonize high intertidal areas (Dayton, Reference Dayton1971).

Elevated wave exposure and increased tidal range were also expected to support higher benthic richness compared with sheltered and/or microtidal habitats. The BRT analyses showed that tidal range accounted for a large variability of the cover and richness patterns of benthic intertidal species. High-energy shores were characterized by high cover and richness of understorey algae; thus, an elevated tidal range probably reduced the desiccation stress of small non-corticated filamentous macroalgae, otherwise common in microtidal shores. The tide height also relates to the tidal stream strength and wave energy reaching the shore so lower cover of larger organisms (canopy algae and some invertebrates) is associated with higher tidal range possibly due to mechanical disturbance related to elevated disturbance by waves (Denny et al., Reference Denny, Miller, Stokes, Hunt and Helmuth2003; Helmuth & Denny, Reference Helmuth and Denny2003).

Elevated regional-scale loads of nutrients were expected to increase both cover and richness of understorey algae species, whereas relationship between nutrient loads and cover and richness of canopy-forming species is bell-shaped. Observed nutrient loads displayed a significant separate effect on the cover and richness of the intertidal communities. It was found that elevated nutrient loads increased both cover and richness of understorey species but had a moderate role on canopy algae. Such observations might relate to functional traits of the studied organisms. For example, understorey algae and sessile invertebrates are both characterized by a short lifespan, thus, elevated nutrient availability might be beneficial throughout their life cycle. Beneficial impacts of nutrients on filamentous understorey algae correlates well with previous findings suggesting shifts in floral dominance patterns from perennial to ephemeral species in case of prolonged events of nutrient influx (Gorostiaga & Diez, Reference Gorostiaga and Diez1996; Diez et al., Reference Diez, Secilla, Santolaria and Gorostiaga1999). However, the highest richness values of understorey algae were observed when nutrient loads were the lowest, minimum richness at intermediate loads and moderate richness at highest loads. This could potentially reflect shifts in taxonomic composition along the eutrophication gradient as communities under pristine conditions differ significantly from those under eutrophic conditions.

The increase in nutrient loads was positively correlated with the cover of canopy algae, but the slope of the function was not as steep when compared with understorey algae. Canopy algae have internal nutrient reserves and might therefore be less influenced by external nutrient availability (Pedersen & Borum, Reference Pedersen and Borum1996). Due to relative stability in internal nutrient reserves the canopy algae are favoured over smaller ephemeral species at low nutrient environments. On the other hand, the highest canopy richness was observed at intermediate nutrient levels that resulted in a bell-shaped relationship. In oligotrophic conditions an increment of nutrients promotes richness by directly improving algal growth conditions and indirectly increasing habitat heterogeneity (Field et al., Reference Field, Behrenfeld, Randerson and Falkowski1998; Kotta et al., Reference Kotta, Möller, Orav-Kotta and Pärnoja2014) whereas nutrient overload is detrimental for many canopy-forming algae and a steep decrease in cover and richness values can occur. The latter phenomenon might be a result of algal blooms or increased epiphyte overgrowth that reduces light availability for macroalgae and/or induces supplemental herbivory (Valiela et al., Reference Valiela, Mcclelland, Hauxwell, Behr, Hersh and Foreman1988; Orav-Kotta & Kotta, Reference Orav-Kotta and Kotta2004; Kraufvelin et al., Reference Kraufvelin, Ruuskanen, Nappu and Kiirikki2007). Ultimately, this leads to competitive exclusion where only the best competitors for light are favoured (Abrams, Reference Abrams1995).

Cover of sessile invertebrates was highest at intermediate nutrient load. Here, several mechanisms can generate the observed functional form relationship but the exact mechanism can only be validated experimentally. Most plausibly, filamentous understorey macroalgae outcompete invertebrates for space in areas of high nutrient loads.

Finally, it was expected that elevated temperatures and increased temperature anomalies would facilitate understorey algal species but hinder the growth of canopy-forming algae. The BRT models demonstrated that temperature described a large proportion of the variability of intertidal communities, especially those of canopy algae. It was found that climate-induced increases in seawater temperature and associated weather anomalies did not seem to have any effect on the cover of canopy algae whereas a very strong negative effect on canopy richness was observed. This suggests that the maintenance of macroalgae-related functions of understorey assemblages in coastal environments (Britton-Simmons, Reference Britton-Simmons2006; Wikström & Kautsky, Reference Wikström and Kautsky2007; Golléty et al., Reference Golléty, Migné and Davoult2008) may be threatened by climate change as species richness is often considered as a reliable measure to reflect the state of an ecosystem (Maestre et al., Reference Maestre, Quero, Gotelli, Escudero, Ochoa, Delgado-Baquerizo, García-Gómez, Bowker, Soliveres, Escolar, García-Palacios, Berdugo, Valencia, Gozalo, Gallardo, Aguilera, Arredondo, Blones, Boeken, Bran, Conceição, Cabrera, Chaieb, Derak, Eldridge, Espinosa, Florentino, Gaitán, Gatica, Ghiloufi, Gómez-González, Gutiérrez, Hernández, Huang, Huber-Sannwald, Jankju, Miriti, Monerris, Mau, Morici, Naseri, Ospina, Polo, Prina, Pucheta, Ramírez-Collantes, Romão, Tighe, Torres-Díaz, Val, Veiga, Wang and Zaady2012) and is regularly found as a basis of management advice. It is important to note that species richness incorporates taxonomic identity of species in an ecosystem that is also linked to functional traits of organisms. If the richness of canopy algae should decrease then it is highly questionable whether sustained rates of high canopy cover could provide similar functions required by the ecosystem as well as whether full recovery is likely. For example, in temperate Western Australia a canopy dominated temperate reef system switched to a turf dominated state due to an exceptional heatwave in 2011 and the recovery to the previous state is further hindered by high grazing rates of herbivores (Smale & Wernberg, Reference Smale and Wernberg2013; Wernberg et al., Reference Wernberg, Smale, Tuya, Thomsen, Langlois, de Bettignies, Bennett and Rousseaux2013). Increased frequency of disturbance events and rising temperatures are expected to cause species range shifts (Lima et al., Reference Lima, Ribeiro, Queiroz, Hawkins and Santos2007a) and pose a significant threat to the distribution of canopy-forming algae and jeopardize the stability of algal-dominated marine regions as well as the overall food-web stability (Wernberg et al., Reference Wernberg, Thomsen, Tuya, Kendrick, Staehr and Toohey2010; Smale et al., Reference Smale, Burrows, Moore, O'Connor and Hawkins2013). Unexpectedly, the patterns of understorey algae were essentially uncoupled with those of elevated temperature anomalies. One likely explanation for such a phenomenon could be that filamentous algal species dominate in understorey communities and rapidly replace one another when environmental conditions become unfavourable to some taxa. The major loss of canopy algae highlights the need to understand whether filamentous algae could compensate for such a loss and provide similar functions to an ecosystem as canopy algae. For example, it has been demonstrated that filamentous algae can be more productive than canopy-forming algae on a biomass basis whereas the annual net production per area is two to seven times lower than for canopy-forming algae (Copertino et al., Reference Copertino, Connell and Cheshire2005). Likewise, it has been shown that filamentous algae may host similar invertebrate species as perennial macroalgal species (Wikström & Kautsky, Reference Wikström and Kautsky2007). However, due to their seasonal occurrence, ephemeral species are not likely to replace the majority of functions provided by perennial species such as habitat formation and associated diversity (Råberg & Kautsky, Reference Råberg and Kautsky2007; Arnold et al., Reference Arnold, Teagle, Brown and Smale2016).

All components of the environment and biota are interlinked and depend on each other. Thus, any model that resolves only separate effects of the environment on biota is only a simplified representation of reality. In this study, we assessed the importance of different environmental interactions and found that, surprisingly, only a few interactive effects emerged as being important. For the ephemeral species, the effect of nutrient loads varied along a gradient of tidal level. This may be partly due to the fact that the sensitivity of habitats to eutrophication decreases with the intensity of water exchange (e.g. Tett et al., Reference Tett, Gilpin, Svendsen, Erlandsson, Larsson, Kratzer, Fouilland, Janzen, Lee, Grenz, Newton, Ferreira, Fernandes and Scory2003) and areas higher in the shore have lower exposure to water enriched with nutrients. Alternatively, elevated intensities of other stressors in the higher shore levels (including elevated rates of desiccation or ice scour, the effects of nutrient loads) are expected to be weaker (Jutterström et al., Reference Jutterström, Andersson, Omstedt and Malmaeus2014). For the canopy-forming species, tidal level moderated the response of temperature to species cover. The high inter-tidal level–temperature interaction had the weakest effect on the cover of canopy algae as these algae can be highly adapted to elevated temperatures/desiccation events (Guenther & Martone, Reference Guenther and Martone2014). Similarly, our study showed that elevated temperature anomalies may counterbalance the effects of ice scour on the richness of canopy algae, promoting an elevated richness of canopy algae at high latitudes (especially in mid and high-intertidal areas). This correlates with recent documentation of the current effects of seawater warming and ice retreat on littoral macrophytes at high latitudes (Weslawski et al., Reference Weslawski, Wiktor and Kotwicki2010). On the other hand, elevated temperature anomalies reduced algal cover at low latitudes (especially in low intertidal areas).

Correlative species distribution models are frequently applied as a sophisticated tool to improve our understanding of the relationship between environment and biota (Elith et al., Reference Elith, Graham, Anderson, Dudík, Ferrier, Guisan, Hijmans, Huettmann, Leathwick, Lehmann, Li, Lohmann, Loiselle, Manion, Moritz, Nakamura, Nakazawa, Mc Overton, Peterson Townsend, Phillips, Richardson, Scachetti-Pereira, Schapire, Soberón, Williams, Wisz and Zimmermann2006, Reference Elith, Kearney and Phillips2010). It is important that the analytical framework has the capacity to quantify the functional form of separate and interactive relationships in order to apply potential reasoning for causations. The machine learning methodology used in this study provides a robust framework to seek for the mechanisms that determine spatial distribution patterns of intertidal species under current environmental conditions. The established functions between environmental gradients and anthropogenic stressors can potentially be used to predict the cover and diversity of intertidal hard bottom communities at the pan-European scale and seek for potential biodiversity hotspots. Such knowledge provides better starting conditions to understand various aquatic ecosystems and contributes towards environmental policies such as marine ecosystem-based management that seeks to manage marine resources in accordance with ecosystem health while providing socio-economic services needed by people (Crowder & Norse, Reference Crowder and Norse2008; Espinosa-Romero et al., Reference Espinosa-Romero, Chan, McDaniels and Dalmer2011; Long et al., Reference Long, Charles and Stephenson2015). However, care should be taken when interpreting cause-effect relationships based on the results of this study as further experimental research is required to validate the obtained results. This is mainly because we were only able to include one temporal replicate of sampled sites which poses a notable limitation for interpretation of results, but is understandable given the scale of this study. Nevertheless, extreme weather conditions such as heavy winter storms in the NE Atlantic Ocean or explicitly warm summers in the Mediterranean Sea could result in atypically sparse algal canopies. Regional deviations from typical disturbance history may therefore affect the credibility of data obtained in one sampling season as such large-scale environmental stochasticity may override the effect of environmental variability operating at smaller spatial scales.

In this study we demonstrate how the distribution patterns of benthic hard bottom algal and invertebrate communities at pan-European scale are affected by various environmental gradients and human-induced stressors. For the majority of large-scale study initiatives the main difficulty is the amount of time required in order to overcome restrictions to access required data and/or to harmonize different databases. This amplifies the need for global science to become more unified and opened through data sharing as data compiled in this study is virtually based on the voluntary initiative of scientists and good will. Projects such as EMBOS, that go beyond the scope and vision of current national and institutional monitoring activities, promote a simple and standardized sampling framework required in order to address complications associated with large-scale studies that examine biological mechanisms of life on earth.

ACKNOWLEDGEMENT

We are thankful to Enrique Ostalé-Valriberas for his kind support during the sampling programme and Caroline Broudin, Fanny Noisette, Aline Migné, Céline Houbin, Laurent Lévêque, Dominique Davoult, François Vandenbosch, Pascal Riera, Elena Maggi, Ofrat Raveh and Elvira Ramos for their help with sampling.

FINANCIAL SUPPORT

This article is based upon work from COST Action ES1003 Development and implementation of a pan-European Marine Biodiversity Observatory System (EMBOS), supported by COST (European Cooperation in Science and Technology). Jonne Kotta, Holger Jänes and Helen Orav-Kotta were supported by Institutional research funding IUT02-20 of the Estonian Research Council and the BONUS project BAMBI, the joint Baltic Sea research and development programme (Art 185), funded jointly from the European Union's Seventh Programme for research, technological development and demonstration and from the Estonian Research Council. Pedro Ribeiro was supported by FCT through post-doctoral grant ref. SFRH/BPD/69232/2010 funded through QREN and COMPETE, and the strategic project UID/MAR/04292/2013 granted to MARE. Jérōme Jourde was financially supported by the Région Poitou-Charentes through CPER funding, La Rochelle University and CNRS. Martina Dal Bello and Lisandro Benedetti-Cecchi were supported by a grant from the Italian Ministry for Research and Education (PRIN project ‘Biocostruzioni costiere: struttura, funzione, e gestione’). Gil Rilov was partly supported by the Israeli Science Foundation (grant number 1217/10) and the Israel Ministry of Environmental Protection. Valentina Kirienko Fernandes de Matos was supported by the Portuguese Science Foundation (FCT) through a doctoral grant (ref. SFRH/BD/86390/2012). Emilia Jankowska and Jan Marcin Węsławski were supported by the Statutory Funds of Institute of Oceanology Polish Academy of Sciences.

CONFLICT OF INTEREST

None.

ETHICAL STANDARDS

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional guides.

ETHICS POLICY

Authors complied with best practice and a high standard of ethics. Authors declare that the work has not been previously published or submitted for publication elsewhere, appropriate ethics conduct and approval were adhered to where field campaigns were conducted such as there were no destruction of habitat or killing of species, none of the text has been copied or paraphrased from other publications.

CONTRIBUTORS

The authors of this article are partners within the COST Action EMBOS on a European Marine Biodiversity Observatory System. First author of the paper has performed statistical analyses and led the writing. Co-authors have made a significant contribution to the conception, design, execution, or interpretation of the reported study. All authors have commented on drafts of the manuscript and have approved the final article to be true and included in the disclosure.

References

REFERENCES

Abrams, P.A. (1995) Monotonic or unimodal diversity-productivity gradients: what does competition theory predict? Ecology 76, 20192027.Google Scholar
Armstrong, R.L. and Knowles, K. (2010) ISLSCP II global sea ice concentration. In Hall, F.G., Collatz, G., Meeson, B., Los, S., De Colstoun, E.B. and Landis, D. (eds) ISLSCP initiative II collection. Data set. Available online [http://daac.ornl.gov/] from Oak Ridge, Tennessee, USA: Oak Ridge National Laboratory Distributed Active Archive Center. doi: 10.3334/ORNLDAAC/981.Google Scholar
Arnold, M., Teagle, H., Brown, M.P. and Smale, D.A. (2016) The structure of biogenic habitat and epibiotic assemblages associated with the global invasive kelp Undaria pinnatifida in comparison to native macroalgae. Biological Invasions 18, 661676.Google Scholar
Bird, C.E., Franklin, E.C., Smith, C.M. and Toonen, R.J. (2013) Between tide and wave marks: a unifying model of physical zonation on littoral shores. PeerJ 1, e154.Google Scholar
Bonsdorff, E. and Pearson, T.H. (1999) Variation in the sublittoral macrozoobenthos of the Baltic Sea along environmental gradients: a functional group approach. Australian Journal of Ecology 24, 312326.CrossRefGoogle Scholar
Britton-Simmons, K.H. (2006) Functional group diversity, resource preemption and the genesis of invasion resistance in a community of marine algae. Oikos 113, 395401.CrossRefGoogle Scholar
Bulleri, F., Benedetti-Cecchi, L., Cusson, M., Maggi, E., Arenas, F., Aspden, R., Bertocci, I., Crowe, T.P., Davoult, D., Eriksson, B.K., Fraschetti, S., Golléty, C., Griffin, J., Jenkins, S.R., Kotta, J., Kraufvelin, P., Molis, M., Sousa Pinto, I., Terlizzi, A., Valdivia, N. and Paterson, D.M. (2012) Temporal stability of European rocky shore assemblages: variation across a latitudinal gradient and the role of habitat-formers. Oikos 121, 18011809.Google Scholar
Conley, D.J., Carstensen, J., Aigars, J., Axe, P., Bonsdorff, E., Eremina, T., Haahti, B.M., Humborg, C., Jonsson, P., Kotta, J., Lännegren, C., Larsson, U., Maximov, A., Rodriguez Medina, M., Lysiak-Pastuszak, E., Remeikaité-Nikiené, N., Walve, J., Wilhelms, S. and Zillén, L. (2011) Hypoxia is increasing in the coastal zone of the Baltic Sea. Environmental Science and Technology 45, 67776783.Google Scholar
Connell, S.D. and Irving, A.D. (2008) Integrating ecology with biogeography using landscape characteristics: a case study of subtidal habitat across continental Australia. Journal of Biogeography 35, 16081621.CrossRefGoogle Scholar
Copertino, M.S., Connell, S.D. and Cheshire, A. (2005) Prevalence and production of turf-forming algae on a temperate subtidal coast. Phycologia 43, 241248.Google Scholar
Crowder, L. and Norse, E. (2008) Essential ecological insights for marine ecosystem-based management and marine spatial planning. Marine Policy 32, 772778.Google Scholar
Crowe, T.P., Cusson, M., Bulleri, F., Davoult, D., Arenas, F., Aspden, R., Benedetti-Cecchi, L., Bevilacqua, S., Davidson, I., Defew, E., Fraschetti, S., Golléty, C., Griffin, J.N., Herkül, K., Kotta, J., Migné, A., Molis, M., Nicol, S.K., Noël, L.M.-L.J., Sousa Pinto, I., Valdivia, N., Vaselli, S. and Jenkins, S.R. (2013) Large-scale variation in combined impacts of canopy loss and disturbance on community structure and ecosystem functioning. PLoS ONE 8, e66238.Google Scholar
Dal Bello, M., Leclerc, J.-C., Benedetti-Cecchi, L., Arvanitidis, C., van Avesaath, P., Bachelet, G., Bojanić, N., Como, S., Coppa, S., Coughlan, J., Crowe, T., Degraer, S., Espinosa, F., Faulwetter, S., Frost, M., Guinda, X., Ikauniece, A., Jankowska, E., Jourde, J., Kerckhof, F., Kotta, J., Lavesque, N., de Lucia, A., Magni, P., Fernandes de Matos, V.K., Orav-Kotta, H., Pavloudi, C., Pedrotti, M.L., Peleg, O., de la Pena, J.A.J., Puente, A., Ribeiro, P., Rigaut-Jalabert, F., Rilov, G., Rousou, M., Rubal, M., Ruginis, T., Pérez-Ruzafa, A., Silva, T., Simon, N., Sousa-Pinto, I., Troncoso, J., Warzocha, J., Weslawski, J.M. and Hummel, H. (2016) Consistent patterns of spatial variability between Atlantic and Mediterranean rocky shores. Journal of the Marine Biological Association of the United Kingdom, In press.Google Scholar
Danovaro, R., Canals, M., Gambi, C., Heussner, S., Lampadariou, N. and Vanreusel, A. (2009) Exploring benthic biodiversity patterns and hotspots on European margin slopes. Oceanography 22, 1625.Google Scholar
Dayton, P.K. (1971) Competition, disturbance, and community organization: the provision and subsequent utilization of space in a rocky intertidal community. Ecological Monographs 41, 351389.CrossRefGoogle Scholar
Dayton, P.K. (1975) Experimental evaluation of ecological dominance in a rocky intertidal algal community. Ecological Monographs 45, 137159.Google Scholar
Denny, M.W., Miller, L.P., Stokes, M.D., Hunt, L.J.H. and Helmuth, B.S.T. (2003) Extreme water velocities: topographical amplification of wave-induced flow in the surf zone of rocky shores. Limnology and Oceanography 48, 18.Google Scholar
Diez, I., Secilla, A., Santolaria, A., and Gorostiaga, J.M. (1999) Phytobenthic intertidal community structure along an environmental pollution gradient. Marine Pollution Bulletin 38, 463472.Google Scholar
Elith, J., Graham, C.H., Anderson, R.P., Dudík, M., Ferrier, S., Guisan, A., Hijmans, R.J., Huettmann, F., Leathwick, J.R., Lehmann, A., Li, J., Lohmann, L.G., Loiselle, B.A., Manion, G., Moritz, C., Nakamura, M., Nakazawa, Y., Mc Overton, J.C.M., Peterson Townsend, A., Phillips, S.J., Richardson, K., Scachetti-Pereira, R., Schapire, R.E., Soberón, J., Williams, S., Wisz, M.S. and Zimmermann, N.E. (2006) Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129151.Google Scholar
Elith, J., Kearney, M. and Phillips, S. (2010) The art of modelling range-shifting species. Methods in Ecology and Evolution 1, 330342.Google Scholar
Elith, J., Leathwick, J.R. and Hastie, T. (2008) A working guide to boosted regression trees. Journal of Animal Ecology 77, 802813.CrossRefGoogle ScholarPubMed
Espinosa-Romero, M.J., Chan, K.M.A., McDaniels, T. and Dalmer, D.M. (2011) Structuring decision-making for ecosystem-based management. Marine Policy 35, 575583.CrossRefGoogle Scholar
Field, C.B., Behrenfeld, M.J., Randerson, J.T. and Falkowski, P. (1998) Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237240.CrossRefGoogle ScholarPubMed
Golléty, C., Migné, A. and Davoult, D. (2008) Benthic metabolism on a sheltered rocky shore: role of the canopy in the carbon budget. Journal of Phycology 44, 11461153.Google Scholar
Gorostiaga, J.M. and Diez, I. (1996) Changes in the sublittoral benthic marine macroalgae in the polluted area of Abra de Bilbao and proximal coast (Northern Spain). Marine Ecology Progress Series 130, 157167.Google Scholar
Goulletquer, P., Gros, P., Boeuf, G. and Weber, J. (2014) Biodiversity in the marine environment. Amsterdam: Springer.CrossRefGoogle Scholar
Guenther, R.J. and Martone, P.T. (2014) Physiological performance of intertidal coralline algae during a simulated tidal cycle. Journal of Phycology 50, 310321.CrossRefGoogle ScholarPubMed
Halpern, B.S., Walbridge, S., Selkoe, K.A., Kappel, C.V., Micheli, F., D'Agrosa, C., Bruno, J.F., Casey, K.S., Ebert, C., Fox, H.E., Fujita, R., Heinemann, D., Lenihan, H.S., Madin, E.M.P., Perry, M.T., Selig, E.R., Spalding, M., Steneck, R. and Watson, R. (2008) A global map of human impact on marine ecosystems. Science 319, 948952.Google Scholar
Hastie, T., Tibshirani, R. and Friedman, J.H. (2009) The elements of statistical learning: data mining, inference, and prediction. New York, NY: Springer-Verlag.Google Scholar
Hawkins, S.J. and Hartnoll, R.G. (1985) Factors determining the upper limits of intertidal canopy-forming algae. Marine Ecology Progress Series 20, 265271.Google Scholar
Heaven, C.S. and Scrosati, R.A. (2008) Benthic community composition across gradients of intertidal elevation, wave exposure, and ice scour in Atlantic Canada. Marine Ecology Progress Series 369, 1323.Google Scholar
Helmuth, B. and Denny, M.W. (2003) Predicting wave exposure in the rocky intertidal zone: do bigger waves always lead to larger forces? Limnology and Oceanography 48, 13381345.CrossRefGoogle Scholar
Herkül, K., Kotta, J., Kutser, T. and Vahtmäe, E. (2013) Relating remotely sensed optical variability to marine benthic biodiversity. PLoS ONE 8, e55624.Google Scholar
Holyoak, M., Leibold, M.A., Moquet, N., Holt, R.D. and Hoopes, M.F. (2005) Metacommunities: a framework for large-scale community ecology. In Holyoak, M., Leibold, M.A., Moquet, N. and Holt, R.D. (eds) Metacommunities: spatial dynamics and ecological communities. Chicago, IL: University of Chicago Press, pp. 134.Google Scholar
Hoogwijk, M., Faaij, A., van den Broek, R., Berndes, G., Gielen, D. and Turkenburg, W. (2003) Exploration of the ranges of the global potential of biomass for energy. Biomass and Bioenergy 25, 119133.CrossRefGoogle Scholar
Howarth, R.W., Billen, G., Swaney, D., Townsend, A., Jaworski, N., Lajtha, K., Downing, J.A., Elmgren, R., Caraco, N., Jordan, T., Berendse, F., Freney, J., Kudeyarov, V., Murdoch, P. and Zhu, Z-L. (1996) Regional nitrogen budgets and riverine N & P fluxes for the drainages to the North Atlantic Ocean: natural and human influences. Biogeochemistry 35, 75139.Google Scholar
Hubbert, M.K. (1949) Energy from fossil fuels. Science 109, 103109.Google Scholar
Johnston, E.L., Keough, M. and Qian, P. (2002) Maintenance of species dominance through pulse disturbances to a sessile marine invertebrate assemblage in Port Shelter, Hong Kong. Marine Ecology Progress Series 226, 103114.Google Scholar
Jutterström, S., Andersson, H.C., Omstedt, A. and Malmaeus, J.M. (2014) Multiple stressors threatening the future of the Baltic Sea–Kattegat marine ecosystem: implications for policy and management actions. Marine Pollution Bulletin 86, 468480.CrossRefGoogle ScholarPubMed
Kaiser, M.J., Attrill, M.J., Jennings, S., Thomas, D.N., Barnes, D.K.A., Brierley, A.S., Hiddink, J.G., Kaartokallio, H., Polunin, N.V.C. and Raffaelli, D.G. (2011) Marine ecology: processes, systems, and impacts. Oxford: Oxford University Press.Google Scholar
Karydis, M. and Kitsiou, D. (2012) Eutrophication and environmental policy in the Mediterranean Sea: a review. Environmental Monitoring and Assessment 184, 49314984.Google Scholar
Kotta, J., Kutser, T., Teeveer, K., Vahtmäe, E. and Pärnoja, M. (2013) Predicting species cover of marine macrophyte and invertebrate species combining hyperspectral remote sensing, machine learning and regression techniques. PLoS ONE 8, e63946.Google Scholar
Kotta, J., Möller, T., Orav-Kotta, H. and Pärnoja, M. (2014) Realized niche width of a brackish water submerged aquatic vegetation under current environmental conditions and projected influences of climate change. Marine Environmental Research 102, 88101.CrossRefGoogle ScholarPubMed
Kotta, J. and Witman, J. (2009) Regional-scale patterns. In Wahl, M. (ed.) Marine hard bottom communities. Berlin: Springer-Verlag, pp. 8999.Google Scholar
Kraufvelin, P., Ruuskanen, A.T., Nappu, N. and Kiirikki, M. (2007) Winter colonisation and succession of filamentous macroalgae on artificial substrates and possible relationships to Fucus vesiculosus settlement in early summer. Estuarine Coastal and Shelf Science 72, 665674.CrossRefGoogle Scholar
Lima, F.P., Ribeiro, P.A., Queiroz, N., Hawkins, S.J. and Santos, A.M. (2007a) Do distributional shifts of northern and southern species of algae match the warming pattern? Global Change Biology 13, 25922604.Google Scholar
Lima, F.P., Ribeiro, P.A., Queiroz, N., Xavier, R., Tarroso, P., Hawkins, S.J. and Santos, A.M. (2007b) Modelling past and present geographical distribution of the marine gastropod Patella rustica as a tool for exploring responses to environmental change. Global Change Biology 13, 20652077.Google Scholar
Long, R.D., Charles, A. and Stephenson, R.L. (2015) Key principles of marine ecosystem-based management. Marine Policy 57, 5360.Google Scholar
Maestre, F.T., Quero, J.L., Gotelli, N.J., Escudero, A., Ochoa, V., Delgado-Baquerizo, M., García-Gómez, M., Bowker, M.A., Soliveres, S., Escolar, C., García-Palacios, P., Berdugo, M., Valencia, E., Gozalo, B., Gallardo, A., Aguilera, L., Arredondo, T., Blones, J., Boeken, B., Bran, D., Conceição, A.A., Cabrera, O., Chaieb, M., Derak, M., Eldridge, D.J., Espinosa, C.I., Florentino, A., Gaitán, J., Gatica, M.G., Ghiloufi, W., Gómez-González, S., Gutiérrez, J.R., Hernández, R.M., Huang, X., Huber-Sannwald, E., Jankju, M., Miriti, M., Monerris, J., Mau, R.L., Morici, E., Naseri, K., Ospina, A., Polo, V., Prina, A., Pucheta, E., Ramírez-Collantes, D.A., Romão, R., Tighe, M., Torres-Díaz, C., Val, J., Veiga, J.P., Wang, D. and Zaady, E. (2012) Plant species richness and ecosystem multifunctionality in global drylands. Science 335, 214218.Google Scholar
Menge, B.A. and Sutherland, J.P. (1987) Community regulation: variation in disturbance, competition and predation in relation to environmental stress and recruitment. American Naturalist 130, 730757.Google Scholar
Narayanaswamy, B.E., Coll, M., Danovaro, R., Davidson, K., Ojaveer, H. and Renaud, P.E. (2013) Synthesis of knowledge on marine biodiversity in European Seas: from census to sustainable management. PLoS ONE 8, e58909.Google Scholar
Orav-Kotta, H. and Kotta, J. (2004) Food and habitat choice of the isopod Idotea baltica in the northeastern Baltic Sea. Hydrobiologia 514, 7985.Google Scholar
Pedersen, M.F. and Borum, J. (1996) Nutrient control of algal growth in estuarine waters. Nutrient limitation and the importance of nitrogen requirements and nitrogen storage among phytoplankton and species of macroalgae. Marine Ecology Progress Series 142, 261272.Google Scholar
Peterson, A.T. (2003) Predicting the geography of species’ invasions via ecological niche modeling. Quarterly Review of Biology 78, 419433.Google Scholar
Råberg, S. and Kautsky, L. (2007) A comparative biodiversity study of the associated fauna of perennial fucoids and filamentous algae. Estuarine, Coastal and Shelf Science 73, 249258.Google Scholar
Ramos, E., Puente, A., Guinda, X. and Juanes, J.A. (2016a) A hierarchical classification system along the NE Atlantic coast: focusing on the local scale (Cantabria, N Spain). European Journal of Phycology, in press. doi: 10.1080/09670262.2016.1221469.Google Scholar
Ramos, E., Puente, A. and Juanes, J. (2016b) An ecological classification of rocky shores at a regional scale: a predictive tool for management of conservation values. Marine Ecology 37, 311328.Google Scholar
RDC Team (2013) R: a language and environment for statistical computing. Vienna: The R Foundation for Statistical Computing. Available at: http://www.r-project.org/.Google Scholar
Rees, H.L., Pendle, M.A., Waldock, R., Limpenny, D.S. and Boyd, S.E. (1999) A comparison of benthic biodiversity in the North Sea, English Channel, and Celtic Seas. ICES Journal of Marine Science 56, 228246.CrossRefGoogle Scholar
Ricklefs, R.E. and Schluter, D. (1993) Species diversity: regional and historical influences. In Ricklefs, R.E. and Schluter, D. (eds) Species diversity in ecological communities: historical and geographical perspectives. Chicago, IL: University of Chicago Press, pp. 350364.Google Scholar
Roy, K., Jablonski, D., Valentine, J.W. and Rosenberg, G. (1998) Marine latitudinal diversity gradients: tests of causal hypotheses. Proceedings of the National Academy of Sciences USA 95, 36993702.Google Scholar
Russell, R., Wood, S.A., Allison, G. and Menge, B.A. (2006) Scale, environment and trophic status: the context dependency of community saturation in rocky intertidal communities. American Naturalist 167, E158E170.Google Scholar
Sandman, A.N., Wikström, S.A., Blomqvist, M., Kautsky, H. and Isaeus, M. (2013) Scale-dependent influence of environmental variables on species distribution: a case study on five coastal benthic species in the Baltic Sea. Ecography 36, 354363.Google Scholar
Schiele, K., Darr, A., Zettler, M.L., Berg, T., Blomqvist, M., Daunys, D., Jermakovs, V., Korpinen, S., Kotta, J., Nygård, H., von Weber, M., Voss, J. and Warzocha, J. (2016) Rating species sensitivity throughout gradient systems – a consistent approach for the Baltic Sea. Ecological Indicators 61, 447455. doi: 10.1016/j.ecolind.2015.09.046.Google Scholar
Smale, D.A., Burrows, M.T., Evans, A.J., King, N., Sayer, M.D.J., Yunnie, A.L.E. and Moore, P.J. (2016) Linking environmental variables with regional-scale variability in ecological structure and standing stock of carbon within kelp forests in the United Kingdom. Marine Ecology Progress Series 542, 7995.Google Scholar
Smale, D.A., Burrows, M.T., Moore, P., O'Connor, N. and Hawkins, S.J. (2013) Threats and knowledge gaps for ecosystem services provided by kelp forests: a northeast Atlantic perspective. Ecology and Evolution 3, 40164038.CrossRefGoogle ScholarPubMed
Smale, D.A. and Wernberg, T. (2013) Extreme climatic event drives range contraction of a habitat-forming species. Proceedings of the Royal Society B 280, 20122829.Google Scholar
Soberon, J. and Peterson, A.T. (2005) Interpretation of models of fundamental ecological niches and species’ distributional areas. Biodiversity Informatics 2, 110.Google Scholar
Stackhouse, P.W. and Gupta, S.K. (2012) ISLSCP II cloud and meteorology parameters. In Hall, F.G., Collatz, G., Meeson, B., Los, S., De Colstoun, E.B. and Landis, D. (eds) ISLSCP Initiative II Collection. Data set. Available online [http://daac.ornl.gov/] from Oak Ridge, Tennessee, USA: Oak Ridge National Laboratory Distributed Active Archive Center. doi: 10.3334/ORNLDAAC/1073.Google Scholar
Stephenson, T.A. and Stephenson, A. (1972) Life between tidemarks on rocky shores. San Francisco, CA: WH Freeman.Google Scholar
Tett, P., Gilpin, L., Svendsen, H., Erlandsson, C.P., Larsson, U., Kratzer, S., Fouilland, E., Janzen, C., Lee, J.-Y., Grenz, C., Newton, A., Ferreira, J.G., Fernandes, T. and Scory, S. (2003) Eutrophication and some European waters of restricted exchange. Continental Shelf Research 23, 16351671.Google Scholar
Thorner, J., Kumar, L. and Smith, S.D.A. (2014) Impacts of climate-change-driven sea level rise on intertidal rocky reef habitats will be variable and site specific. PLoS ONE 9, e86130.Google Scholar
Tomanek, L. and Helmuth, B. (2002) Physiological ecology of rocky intertidal organisms: a synergy of concepts. Integrative and Comparative Biology 42, 771775.Google Scholar
Underwood, A.J. (1978) A refutation of critical tidal levels as determinants of the structure of intertidal communities on British shores. Journal of Experimental Marine Biology and Ecology 33, 261276.Google Scholar
Underwood, A.J. and Jernakoff, P. (1981) Effects of interactions between algae and grazing gastropods on the structure of a low-shore intertidal algal communities. Oecologia 48, 221233.Google Scholar
Valiela, I., Mcclelland, J., Hauxwell, J., Behr, P.J., Hersh, D. and Foreman, K. (1988) Macroalgal blooms in shallow estuaries: controls and ecophysiological and ecosystem consequences. Limnology and Oceanography 42, 11051118.Google Scholar
Wernberg, T., Smale, D.A., Tuya, F., Thomsen, M.S., Langlois, T.J., de Bettignies, T., Bennett, S. and Rousseaux, C.S. (2013) An extreme climatic event alters marine ecosystem structure in a global biodiversity hotspot. Nature Climate Change 3, 7882.CrossRefGoogle Scholar
Wernberg, T., Thomsen, M.S., Tuya, F., Kendrick, G.A., Staehr, P.A. and Toohey, B.D. (2010) Decreasing resilience of kelp beds along a latitudinal temperature gradient: potential implications for a warmer future. Ecology Letters 13, 685694.Google Scholar
Weslawski, J.M., Wiktor, J. Jr and Kotwicki, L. (2010) Increase in biodiversity in the arctic rocky littoral, Sorkappland, Svalbard, after 20 years of climate warming. Marine Biodiversity 40, 123130.Google Scholar
Whittaker, R.O. (1974) Communities and ecosystems. New York, NY: Macmillan.Google Scholar
Wikström, S.A. and Kautsky, L. (2007) Structure and diversity of invertebrate communities in the presence and absence of canopy-forming Fucus vesiculosus in the Baltic Sea. Estuarine Coastal and Shelf Science 72, 168176.Google Scholar
Zacharias, M.A. and Roff, J.C. (2001) Explanations of patterns of intertidal diversity at regional scales. Journal of Biogeography 28, 471483.Google Scholar
Figure 0

Fig. 1. Map of the sampling stations in the study area. The annual average SST isotherms (°C) are also indicated.

Figure 1

Fig. 2. Average values of the cover and richness of understorey algae, canopy algae and sessile invertebrates in different regions. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal.

Figure 2

Table 1. Relative contribution of environmental predictors in the BRT models of total cover and richness of macroalgae and invertebrates.

Figure 3

Fig. 3. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of understorey algae in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

Figure 4

Fig. 4. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of canopy algae in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

Figure 5

Fig. 5. Standardized functional-form relationships showing effects of the four most important environmental variables on the cover and richness of sessile invertebrates in the study area, whilst all other variables are held at their means. The variables are ordered by their relative contribution in the BRT model (shown in brackets). Upward tickmarks on x-axis show the frequency of distribution of data along this axis. The code of tidal levels is as follows: H – high intertidal, M – mid intertidal, L – low intertidal, S – subtidal. The values of nutrient loads and temperature anomalies are log-transformed and rescaled between 0–1. See the section of methods for further information on environmental variables.

Figure 6

Fig. 6. Partial dependence plots of the BRT model for the cover and richness of understorey algae and sessile invertebrates at different tidal levels in the study area. The values of nutrient loads are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.

Figure 7

Fig. 7. A partial dependence plot of the BRT model for the cover of canopy algae at different tidal levels and a three-dimensional partial dependence plot in the BRT model for the richness of canopy algae along the gradients of water temperature anomalies and ice cover in the study area. The values of temperature anomalies are log-transformed and rescaled between 0–1. See the Methods section for further information on environmental variables.