Quantifying ecosystem-level consequences of ocean acidification and warming by scaling up individual level responses of a predator snail and its trophic interactions

Understanding long‐term, ecosystem‐level impacts of climate change is challenging because experimental research frequently focuses on short‐term, individual‐level impacts in isolation. We address this shortcoming first through an interdisciplinary ensemble of novel experimental techniques to investigate the impacts of 14‐month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterize a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual‐level responses, while acidification had a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual‐level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large‐scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local‐environmental conditions and resource availability. Such changes in macro‐scale distributions cannot be predicted by investigating individual‐level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long‐term, multiscale responses to multiple stressors, in an ecosystem context.

Scaling up experimental ocean acidification and warming research: from individuals to the ecosystem

Understanding long-term, ecosystem-level impacts of climate change is challenging because experimental research frequently focuses on short-term, individual-level impacts in isolation. We address this shortcoming first through an inter-disciplinary ensemble of novel experimental techniques to investigate the impacts of 14-month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterise a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual level responses, while acidification has a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large-scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local environmental conditions and resource availability. Such changes in macro-scale distributions cannot be predicted by investigating individual level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long-term, multi-scale responses to multiple stressors, in an ecosystem context.

Scaling up experimental ocean acidification and warming research: from individuals to the ecosystem

Understanding long-term, ecosystem-level impacts of climate change is challenging because experimental research frequently focuses on short-term, individual-level impacts in isolation. We address this shortcoming first through an inter-disciplinary ensemble of novel experimental techniques to investigate the impacts of 14-month exposure to ocean acidification and warming (OAW) on the physiology, activity, predatory behaviour and susceptibility to predation of an important marine gastropod (Nucella lapillus). We simultaneously estimated the potential impacts of these global drivers on N. lapillus population dynamics and dispersal parameters. We then used these data to parameterise a dynamic bioclimatic envelope model, to investigate the consequences of OAW on the distribution of the species in the wider NE Atlantic region by 2100. The model accounts also for changes in the distribution of resources, suitable habitat and environment simulated by finely resolved biogeochemical models, under three IPCC global emissions scenarios. The experiments showed that temperature had the greatest impact on individual level responses, while acidification has a similarly important role in the mediation of predatory behaviour and susceptibility to predators. Changes in Nucella predatory behaviour appeared to serve as a strategy to mitigate individual level impacts of acidification, but the development of this response may be limited in the presence of predators. The model projected significant large-scale changes in the distribution of Nucella by the year 2100 that were exacerbated by rising greenhouse gas emissions. These changes were spatially heterogeneous, as the degree of impact of OAW on the combination of responses considered by the model varied depending on local environmental conditions and resource availability. Such changes in macro-scale distributions cannot be predicted by investigating individual level impacts in isolation, or by considering climate stressors separately. Scaling up the results of experimental climate change research requires approaches that account for long-term, multi-scale responses to multiple stressors, in an ecosystem context.

Acquisition of continuous plankton recorder data

The Continuous Plankton Recorder has been deployed on a seasonal basis in the north Pacific since 2000, accumulating a database of abundance measurements for over 290 planktonic taxa in over 3,500 processed samples. There is an additional archive of over 10,000 samples available for further analyses. Exxon Valdez Oil Spill Trustee Council financial support has contributed to about half of this tally, through four projects funded since 2002. Time series of zooplankton variables for sub-regions of the survey area are presented together with abstracts of eight papers published using data from these projects. The time series covers a period when the dominant climate signal in the north Pacific, the Pacific Decadal Oscillation (PDO), switched with unusual frequency between warm/positive states (pre-1999 and 2003-2006) and cool/negative states (1999-2002 and 2007). The CPR data suggest that cool negative years show higher biomass on the shelf and lower biomass in the open ocean, while the reverse is true in warm (PDO positive) years with lower shelf biomass (except 2005) and higher oceanic biomass. In addition, there was a delay in plankton increase on the Alaskan shelf in the colder spring of 2007, compared to the warmer springs of the preceding years. In warm years, smaller species of copepods which lack lipid reserves are also more common. Availability of the zooplankton prey to higher trophic levels (including those that society values highly) is therefore dependent on the timing of increase and peak abundance, ease of capture and nutritional value. Previously published studies using these data highlight the wide-ranging applicability of CPR data and include collaborative studies on; phenology in the key copepod species Neocalanus plumchrus, descriptions of distributions of decapod larvae and euphausiid species, the effects of hydrographic features such as mesoscale eddies and the North Pacific Current on plankton populations and a molecularbased investigation of macro-scale population structure in N. cristatus. The future funding situation is uncertain but the value of the data and studies so far accumulated is considerable and sets a strong foundation for further studies on plankton dynamics and interactions with higher trophic levels in the northern Gulf of Alaska.

Global latitudinal variations in marine copepod diversity and environmental factors

Latitudinal gradients in diversity are among the most striking features in ecology. For terrestrial species, climate (i.e. temperature and precipitation) is believed to exert a strong influence on the geographical distributions of diversity through its effects on energy availability. Here, we provide the first global description of geographical variation in the diversity of marine copepods, a key trophic link between phytoplankton and fish, in relation to environmental variables. We found a polar-tropical difference in copepod diversity in the Northern Hemisphere where diversity peaked at subtropical latitudes. In the Southern Hemisphere, diversity showed a tropical plateau into the temperate regions. This asymmetry around the Equator may be explained by climatic conditions, in particular the influence of the Inter-Tropical Convergence Zone, prevailing mainly in the northern tropical region. Ocean temperature was the most important explanatory factor among all environmental variables tested, accounting for 54 per cent of the variation in diversity. Given the strong positive correlation between diversity and temperature, local copepod diversity, especially in extra-tropical regions, is likely to increase with climate change as their large-scale distributions respond to climate warming.

Placing biodiversity in ecosystem models without getting lost in translation

A key challenge to progressing our understanding of biodiversity’s role in the sustenance of ecosystem function is the extrapolation of the results of two decades of dedicated empirical research to regional, global and future landscapes. Ecosystem models provide a platform for this progression, potentially offering a holistic view of ecosystems where, guided by the mechanistic understanding of processes and their connection to the environment and biota, large-scale questions can be investigated. While the benefits of depicting biodiversity in such models are widely recognized, its application is limited by difficulties in the transfer of knowledge from small process oriented ecology into macro-scale modelling. Here, we build on previous work, breaking down key challenges of that knowledge transfer into a tangible framework, highlighting successful strategies that both modelling and ecology communities have developed to better interact with one another. We use a benthic and a pelagic case-study to illustrate how aspects of the links between biodiversity and ecosystem process have been depicted in marine ecosystem models (ERSEM and MIRO), from data, to conceptualisation and model development. We hope that this framework may help future interactions between biodiversity researchers and model developers by highlighting concrete solutions to common problems, and in this way contribute to the advance of the mechanistic understanding of the role of biodiversity in marine (and terrestrial) ecosystems.

The meso-scale response of subarctic North Pacific seabird community structure to lower trophic level abundance and diversity

Understanding the mechanisms that structure communities and influence biodiversity are fundamental goals of ecology. To test the hypothesis that the abundance and diversity of upper-trophic level predators (seabirds) is related to the underlying abundance and diversity of their prey (zooplankton) and ecosystem-wide energy availability (primary production), we initiated a monitoring program in 2002 that jointly and repeatedly surveys seabird and zooplankton populations across a 7,500 km British Columbia-Bering Sea-Japan transect. Seabird distributions were recorded by a single observer (MH) using a strip-width technique, mesozooplankton samples were collected with a Continuous Plankton Recorder, and primary production levels were derived using the appropriate satellite parameters and the Vertically Generalized Production Model (Behrenfeld and Falkowski 1997). Each trophic level showed clear spatio-temporal patterns over the course of the study. The strongest relationship between seabird abundance and diversity and the lower trophic levels was observed in March/April ('spring') and significant relationships were also found through June/July ('summer'). No discernable relationships were observed during the September/October ('fall') months. Overall, mesozooplankton abundance and biomass explained the dominant portion of seabird abundance and diversity indices (richness, Simpson's Index, and evenness), while primary production was only related to seabird richness. These findings underscore the notion that perturbations of ocean productivity and lower trophic level ecosystem constituents influenced by climate change, such as shifts in timing (phenology) and synchronicity (match-mismatch), could impart far-reaching consequences throughout the marine food web.

Basin scale distribution of zooplankton in the Ligurian Sea (north-western Mediterranean) in late autumn

Mesozooplankton biomass and abundance were evaluated in epipelagic waters at 59 stations covering the Italian sector of the Ligurian Sea (north-western Mediterranean) in December 1990. This region is characterised by a cyclonic circulation which encloses a central divergence zone and is associated with a main thermohaline front offshore the western Ligurian coast. At the end of autumn, mesozooplankton biomass (range: 0.80–4.24 mg DW m−3) and the abundance (range: 83.8–932 ind. m−3) were lower in the divergence zone. On the contrary, in the Ligurian frontal zone at the periphery of the divergence and on the eastern continental shelf the greatest values of biomass and abundance were recorded. Copepods and appendicularians dominated the mesozooplankton community, the main taxa being the copepods Clausocalanus spp. (46% of total zooplankton) and Oithona spp. (15%) and the appendicularian Fritillaria spp. (12%). Three hydrological sub-regions, i.e. the divergence, the eastern continental shelf and the periphery of the divergence, were characterised by different zooplankton communities and characteristic species. Environmental differences between the three zones were mainly related to changes in bottom topography, sea surface temperatures and quantity of particulate organic matter. Vertical mesozooplankton abundance and taxa distribution from the surface to 1,900 m depth were also examined in one station. The results showed that the bulk of the community was concentrated in the upper 200 m, small copepods being dominant particularly in the upper 50 m. The copepod community was more diversified in sub-superficial waters, with a maximum observed in the 200–400 m layer. The distributions of main zooplankton taxa described in epipelagic waters in the eastern Ligurian Sea in autumn were compared with their distribution at surface in the north-western Mediterranean obtained by sampling performed with the Continuous Plankton Recorder in 1997–1999. The analysis of the zooplankton community in CPR samples confirms the dominance of small copepods (Paracalanus spp., Clausocalanus spp., Oithona spp.) and appendicularians in the north-western Mediterranean in late autumn-winter and shows that their distribution is mainly related to the main mesoscale hydrographic features characterising this basin.

Quantifying pursuit diving seabirds' associations with fine-scale physical features in a high tidal energy environment

The rapid increase in the number of tidal stream turbine arrays will create novel and unprecedented levels of anthropogenic activity within habitats characterized by horizontal current speeds exceeding 2 ms−1. However, the potential impacts on pursuit‐diving seabirds exploiting these tidal stream environments remain largely unknown. Identifying similarities between the fine‐scale physical features (100s of metres) suitable for array installations, and those associated with foraging pursuit‐diving seabirds, could identify which species are most vulnerable to either collisions with moving components, or displacement from these installations. A combination of vessel‐based observational surveys, Finite Volume Community Ocean Model outputs and hydroacoustic seabed surveys provided concurrent measures of foraging distributions and physical characteristics at a fine temporal (15 min) and spatial (500 m) resolution across a tidal stream environment suitable for array installations, during both breeding and non‐breeding seasons. These data sets were then used to test for associations between foraging pursuit‐diving seabirds (Atlantic puffins Fratercula arctica, black guillemots Cepphus grylle, common guillemots Uria aalge, European shags Phalacrocorax aristotelis) and physical features. These species were associated with areas of fast horizontal currents, slow horizontal currents, high turbulence, downward vertical currents and also hard–rough seabeds. The identity and strength of associations differed among species, and also within species between seasons, indicative of interspecific and intraspecific variations in habitat use. However, Atlantic puffins were associated particularly strongly with areas of fast horizontal currents during breeding seasons, and European shags with areas of rough–hard seabeds and downward vertical currents during non‐breeding seasons. Synthesis and applications. Atlantic puffins’ strong association with fast horizontal current speeds indicates that they are particularly likely to interact with installations during breeding seasons. Any post‐installation monitoring and mitigation measures should therefore focus on this species and season. The multi‐species associations with high turbulence and downward vertical currents, which often coincide with fast horizontal current speeds, also highlight useful pre‐installation mitigation measures via the omission of devices from these areas, reducing the overall likelihood of interactions. Environmental impact assessments (EIA) generally involve once‐a‐month surveys across 2‐year periods. However, the approaches used in this study show that more focussed surveys can greatly benefit management strategies aiming to reduce the likelihood of negative impacts by facilitating the development of targeted mitigation measures. It is therefore recommended that these approaches contribute towards EIA within development sites.

Broad-scale patterns of sex ratios in Patella spp.: a comparison of range edge and central range populations in the British Isles and Portugal

Sex change, or sequential hermaphroditism, occurs in the plant and animal kingdoms and often determines a predominance of the first sex. Our aim was to explore changes in sex ratios within the range of the species studied: Patella vulgata and Patella depressa. The broad-scale survey of sex with size of limpets covered a range of latitudes from Zambujeira do Mar (southern Portugal) to the English Channel. Indirect evidence was found for the occurrence of protandry in P. vulgata populations from the south of England, with females predominating in larger size-classes; cumulative frequency distributions of males and females were different; sex ratios were biased towards males and smallest sizes of males were smaller than the smallest sizes of females. In contrast in Portugal females were found in most size-classes of P. vulgata. In P. depressa populations from the south coast of England and Portugal females were interspersed across most size-classes; size distributions of males and females and size at first maturity of males and females did not differ. P. depressa did, however, show some indications of the possibility of slight protandry occurring in Portugal. The test of sex ratio variation with latitude indicated that P. vulgata sex ratios might be involved in determining the species range limit, particularly at the equatorward limit since the likelihood of being male decreased from the south coast of England to southern Portugal. Thus at the southern range limit, sperm could be in short supply due to scarcity of males contributing to an Allee effect.

Broad-scale patterns of sex ratios in Patella spp.: a comparison of range edge and central range populations in the British Isles and Portugal

Sex change, or sequential hermaphroditism, occurs in the plant and animal kingdoms and often determines a predominance of the first sex. Our aim was to explore changes in sex ratios within the range of the species studied: Patella vulgata and Patella depressa. The broad-scale survey of sex with size of limpets covered a range of latitudes from Zambujeira do Mar (southern Portugal) to the English Channel. Indirect evidence was found for the occurrence of protandry in P. vulgata populations from the south of England, with females predominating in larger size-classes; cumulative frequency distributions of males and females were different; sex ratios were biased towards males and smallest sizes of males were smaller than the smallest sizes of females. In contrast in Portugal females were found in most size-classes of P. vulgata. In P. depressa populations from the south coast of England and Portugal females were interspersed across most size-classes; size distributions of males and females and size at first maturity of males and females did not differ. P. depressa did, however, show some indications of the possibility of slight protandry occurring in Portugal. The test of sex ratio variation with latitude indicated that P. vulgata sex ratios might be involved in determining the species range limit, particularly at the equatorward limit since the likelihood of being male decreased from the south coast of England to southern Portugal. Thus at the southern range limit, sperm could be in short supply due to scarcity of males contributing to an Allee effect.