Predator–prey reversal: A possible mechanism for ecosystem hysteresis in the North Sea?

Removal of large predatory fishes from marine ecosystems has resulted in persistent ecosystem shifts, with collapsed predator populations and super-abundant prey populations. One explanation for these shifts is reversals of predator–prey roles that generate internal feedbacks in the ecosystems. Pelagic forage fish are often predators and competitors to the young life stages of their larger fish predators. I show that cod recruitment in the North Sea has been negatively related to the spawning-stock biomass of herring for the last 44 years. Herring, together with the abundance of Calanus finmarchicus, the major food for cod larvae, were the main predictors of cod recruitment. These predictors were of equivalent importance, worked additively, and explained different parts of the dynamics in cod recruitment. I suggest that intensive harvesting of cod has released herring from predator control, and that a large population of herring suppresses cod recruitment through predation on eggs and larvae. This feedback mechanism can promote alternative stable states and therefore cause hysteresis to occur under changing conditions; however, harvesting of herring might at present prevent a shift in the ecosystem to a herring-dominated state.

Regional and annual variation in black-legged kittiwake breeding productivity is related to sea surface temperature

Spatiotemporal variation in seabird demographic parameters is often pronounced and may be an important source of information on the state of marine ecosystems. Black-legged kittiwakes Rissa tridactyla in Britain and Ireland show strong regional structure in breeding productivity, and both temporal and spatial variation are probably related to abundance of the principal prey of breeding kittiwakes, the lesser sandeel Ammodytes marinus. Annual regional estimates of sandeel abundance do not exist, prohibiting direct tests of this hypothesis. We examined relationships between kittiwake breeding productivity and 2 potential proxies of sandeel abundance, winter sea surface temperature (SST) and abundance of Calanus copepods, within and among 6 regions in Britain and Ireland from 1986 to 2004. Means and trends in winter SST differed among regions, with higher means and less pronounced increasing trends in western (Atlantic) regions than in eastern (North Sea) regions. A negative relationship between breeding productivity and winter SST in the previous year was found within 2 regions (East Scotland and Orkney), as well as in a cross-regional analysis. Results were inconclusive for Calanus abundance, with a positive relationship in East Scotland and negative in Orkney. These results demonstrate that although a single environmental driver (SST) is related to both within- and between-region variation in a key demographic parameter, regional heterogeneity in SST trends as well as the importance of other factors may lead to highly variable responses. Understanding this heterogeneity is critical for predicting long-term effects of climate change or other anthropogenic drivers on marine ecosystems.

A synthesis of large-scale patterns in the planktonic prey of larval and juvenile cod (Gadus morhua)

Data from 40 published studies of the diet composition of larval and juvenile cod (Gadus morhua) from around the northern North Atlantic were summarized to assess generic patterns in ontogenetic and regional variability in the key prey. The results showed that larvae at the northern edge of the latitudinal range of cod depend primarily on development stages of the copepod Calanus finmarchicus, whilst those at the southern edge depend on Para- and Pseudocalanus species. Juvenile cod preyed on a wider range of taxa than larvae, but euphausiids were the main target prey. Analysis of regional variations in the relative abundances of C. finmarchicus and Para/Pseudocalanus spp. in the plankton, as estimated by the continuous plankton recorder (CPR) surveys, showed a similar geographical pattern to the larval cod stomach contents. Comparison of CPR data from the 1960s and 70s with data from the 1990s showed that the boundary between C. finmarchicus and Para/Pseudocalanus spp. dominance has shifted northwards on both sides of the Atlantic, whilst the abundance of euphausiids in the southern cod stock regions has declined. The results are discussed in relation to regional differences in the response of cod stocks to climate variability.

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.

Empirically modelling the potential effects of changes in temperature and prey availability on the growth of cod larvae in UK shelf seas

It has been hypothesized that changes in zooplankton community structure over the past four decades led to reduced growth and survival of prerecruit Atlantic cod (Gadus morhua) and that this was a key factor underlying poor year classes, contributing to stock collapse, and inhibiting the recovery of stocks around the UK. To evaluate whether observed changes in plankton abundance, species composition and temperature could have led to periods of poorer growth of cod larvae, we explored the effect of prey availability and temperature on early larval growth using an empirical trophodynamic model. Prey availability was parameterized using species abundance data from the Continuous Plankton Recorder. Our model suggests that the observed changes in plankton community structure in the North Sea may have had less impact on cod larval growth, at least for the first 40 days following hatching, than previously suggested. At least in the short term, environmental and prey conditions should be able to sustain growth of cod larvae and environmental changes acting on this early life stage should not limit stock recovery.

Macro-ecology of plankton–seabird associations in the North Pacific Ocean

In conjunction with the North Pacific Continuous Plankton Recorder program, we conducted surveys of seabirds from June 2002 to June 2007. Here, we tested the hypotheses of (i) east–west variations in coupled plankton and seabird abundance, and (ii) that surface-feeding and diving seabirds vary in their relationships to primary productivity and mesozooplankton species abundance and diversity. To test these hypotheses, we developed statistical models for 20 species of seabirds and 12 zooplankton taxonomic groups. Seabird density was highly variable between seasons, but was consistently higher in the western than eastern North Pacific. Seabird diversity was greater in the east. Zooplankton abundance did not differ between regions. We found associations at the “bulk” level between seabird density and net primary productivity, but only one association between seabirds and total zooplankton abundance or diversity. However, we found many relationships between seabird species and the abundance of different zooplankton summarized at the genus or family level. Some of these taxonomic relationships reflect direct predator–prey interactions, while others may reflect zooplankton that serve as ecological indicators of other prey, such as micronekton, upon which the birds may feed. Surface or near-surface feeding, mostly piscivorous seabirds, did not differ systematically from diving, mainly planktivorous seabirds in their zooplankton associations. Seabirds apparently respond to zooplankton taxonomic groupings more so than bulk zooplankton characteristics, such as abundance or diversity. Macro-ecological studies of remote marine ecosystems using zooplankton and seabirds as ecological indicators provide a framework for understanding and assessing spatial and temporal variations in these difficult-to-study pelagic environments.

Prey landscapes help identify potential foraging habitats for leatherback turtles in the NE Atlantic

Identifying key marine megavertebrate habitats has become ever more important as concern increases regarding global fisheries bycatch and accelerated climate change. This will be aided by a greater understanding of the patterns and processes determining the spatiotemporal distribution of species of conservation concern. We identify probable foraging grounds for leatherback turtles in the NE Atlantic using monthly landscapes of gelatinous organism distribution constructed from Continuous Plankton Recorder Survey data. Using sightings data (n = 2013 records, 1954 to 2003) from 9 countries (UK, Ireland, France, Belgium, The Netherlands, Denmark, Germany, Norway and Sweden), we show sea surface temperatures of approximately 10 to 12 degree C most likely indicate the lower thermal threshold for accessible habitats during seasonal foraging migrations to high latitudes. Integrating maps of gelatinous plankton as a possible indicator of prey distribution with thermal tolerance parameters demonstrates the dynamic (spatial and temporal) nature of NE Atlantic foraging habitats. We highlight the importance of body size- related thermal constraints in structuring leatherback foraging populations and demonstrate a latitudinal gradient in body size (Bergmann's rule) where smaller animals are excluded from higher latitude foraging areas. We highlight the marine area of the European continental shelf edge as being both thermally accessible and prey rich, and therefore potentially supporting appreciable densities of foraging leatherbacks, with some suitable areas not yet extensively surveyed.

Reply to comment: is climate change the most likely driver of range expansion of a critically endangered top predator in northeast Atlantic waters?

The comment by Votier et al. (2008) on our recently published article (Wynn et al. 2007) contains two main criticisms: (i) that our analytical approach is inappropriate and (ii) that we have failed to acknowledge other factors that may have contributed to the change in Balearic Shearwater numbers recorded throughout northwest European waters. We strongly disagree with both these criticisms.

Investigating the inter-relationships between water attenuated irradiance, primary production and DMS(P)

Both solar irradiance and primary production have been proposed as independent controls on seawater dimethyl sulphide (DMS) and dimethylsulphoniopropionate (DMSP) concentrations. However, irradiance also drives photosynthesis, and thus influences a complex set of inter-related processes that modulate marine DMS. We investigate the potential inter-relationships between the rate of primary production (carbon assimilation), water-attenuated irradiance and DMS/DMSP dynamics by applying correlation analysis to a high resolution, concurrently sampled in situ data set from a range of latitudes covering multiple biogeochemical provinces from 3 of the 4 Longhurst biogeochemical domains. The combination of primary production (PP) and underwater irradiance (Iz) within a multivariate regression model is able to explain 55% of the variance in DMS concentrations from all depths within the euphotic zone and 66% of the variance in surface DMS concentrations. Contrary to some previous studies we find a variable representing biological processes is necessary to better account for the variance in DMS. We find that the inclusion of Iz accounts for variance in DMS that is independent from the variance explained by PP. This suggests an important role for solar irradiance (beyond the influence of irradiance upon primary production) in mediating the relationship between the productivity of the ecosystem, DMS/DMSP production and ambient seawater DMS concentrations.

Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link

Exploring climate and anthropogenic impacts on marine ecosystems requires an understanding of how trophic components interact. However, integrative end-to-end ecosystem studies (experimental and/or modelling) are rare. Experimental investigations often concentrate on a particular group or individual species within a trophic level, while tropho-dynamic field studies typically employ either a bottom-up approach concentrating on the phytoplankton community or a top-down approach concentrating on the fish community. Likewise the emphasis within modelling studies is usually placed upon phytoplankton-dominated biogeochemistry or on aspects of fisheries regulation. In consequence the roles of zooplankton communities (protists and metazoans) linking phytoplankton and fish communities are typically under-represented if not (especially in fisheries models) ignored. Where represented in ecosystem models, zooplankton are usually incorporated in an extremely simplistic fashion, using empirical descriptions merging various interacting physiological functions governing zooplankton growth and development, and thence ignoring physiological feedback mechanisms. Here we demonstrate, within a modelled plankton food-web system, how trophic dynamics are sensitive to small changes in parameter values describing zooplankton vital rates and thus the importance of using appropriate zooplankton descriptors. Through a comprehensive review, we reveal the mismatch between empirical understanding and modelling activities identifying important issues that warrant further experimental and modelling investigation. These include: food selectivity, kinetics of prey consumption and interactions with assimilation and growth, form of voided material, mortality rates at different age-stages relative to prior nutrient history. In particular there is a need for dynamic data series in which predator and prey of known nutrient history are studied interacting under varied pH and temperature regimes.

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.