Phenological trends and trophic mismatch across multiple levels of a North Sea pelagic food web

Differential phenological responses to climate among species are predicted to disrupt trophic interactions, but datasets to evaluate this are scarce. We compared phenological trends for species from 4 levels of a North Sea food web over 24 yr when sea surface temperature (SST) increased significantly. We found little consistency in phenological trends between adjacent trophic levels, no significant relationships with SST, and no significant pairwise correlations between predator and prey phenologies, suggesting that trophic mismatching is occurring. Finer resolution data on timing of peak energy demand (mid-chick-rearing) for 5 seabird species at a major North Sea colony were compared to modelled daily changes in length of 0-group (young of the year) lesser sandeels Ammodytes marinus. The date at which sandeels reached a given threshold length became significantly later during the study. Although the phenology of all the species except shags also became later, these changes were insufficient to keep pace with sandeel length, and thus mean length (and energy value) of 0-group sandeels at mid-chick-rearing showed net declines. The magnitude of declines in energy value varied among the seabirds, being more marked in species showing no phenological response (shag, 4.80 kJ) and in later breeding species feeding on larger sandeels (kittiwake, 2.46 kJ) where, due to the relationship between sandeel length and energy value being non-linear, small reductions in length result in relatively large reductions in energy. However, despite the decline in energy value of 0-group sandeels during chick-rearing, there was no evidence of any adverse effect on breeding success for any of the seabird species. Trophic mismatch appears to be prevalent within the North Sea pelagic food web, suggesting that ecosystem functioning may be disrupted.

Scale-dependent response diversity of seabirds to prey in the North Sea

Functional response diversity is defined as the diversity of responses to environmental change among species that contribute to the same ecosystem function. Because different ecological processes dominate on different spatial and temporal scales, response diversity is likely to be scale dependent. Using three extensive data sets on seabirds, pelagic fish, and zooplankton, we investigate the strength and diversity in the response of seabirds to prey in the North Sea over three scales of ecological organization. Two-stage analyses were used to partition the variance in the abundance of predators and prey among the different scales of investigation: variation from year to year, variation among habitats, and variation on the local patch scale. On the year-to-year scale, we found a strong and synchronous response of seabirds to the abundance of prey, resulting in low response diversity. Conversely, as different seabird species were found in habitats dominated by different prey species, we found a high diversity in the response of seabirds to prey on the habitat scale. Finally, on the local patch scale, seabirds were organized in multispecies patches. These patches were weakly associated with patches of prey, resulting in a weak response strength and a low response diversity. We suggest that ecological similarities among seabird species resulted in low response diversity on the year-to-year scale. On the habitat scale, we suggest that high response diversity was due to interspecific competition and niche segregation among seabird species. On the local patch scale, we suggest that facilitation with respect to the detection and accessibility of prey patches resulted in overlapping distribution of seabirds but weak associations with prey. The observed scale dependencies in response strength and diversity have implications for how the seabird community will respond to different environmental disturbances.

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.

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.

Spatial demography of Calanus finmarchicus in the Irminger Sea

Continuous Plankton Recorder data suggest that the Irminger Sea supports a major proportion of the surface-living population of the copepod Calanus finmarchicus in the northern North Atlantic, but there have been few studies of its population dynamics in the region. In this paper, we document the seasonal changes in the demographic structure of C finmarchicus in the Irminger Sea from a field programme during 2001/2002, and the associations between its developmental stages and various apparent bio-physical zones. Overwintering stages were found widely at depth (>500 m) across the Irminger Sea, and surviving females were widely distributed in the surface waters the following spring. However, recruitment of the subsequent generation was concentrated around the fringes of the Irminger Sea basin, along the edges of the Irminger and East Greenland Currents, and not in the central basin. In late summer animals were found descending back to overwintering depths in the Central Irminger Sea. The key factors dictating this pattern of recruitment appear to be (a) the general circulation regime, (b) predation on eggs in the spring, possibly by the surviving GO stock, and (c) mortality of first feeding naupliar stages in the central basin where food concentrations appear to be low throughout the year. We compared the demographic patterns in 2001/2002 with observations from the only previous major survey in 1963 and with data from the Continuous Plankton Recorder (CPR) surveys. In both previous data sets, the basic structure of GO ascent from the central basin and G1 recruitment around the fringes was a robust feature, suggesting that it is a recurrent phenomenon. The Irminger Sea is a complex mixing zone between polar and Atlantic water masses, and it has also been identified as a site of sporadic deep convection. The physical oceanographic characteristics of the region are therefore potentially sensitive to climate fluctuations. Despite this, the abundance of C finmarchicus in the region, as measured by the CPR surveys, appears not to have responded to climate factors linked to the North Atlantic Oscillation Index, in contrast with the stocks in eastern Atlantic areas. We speculate that this may because biological factors (production and mortality), rather than transport processes are the key factors affecting the population dynamics in the Irminger Sea. (c) 2007 Elsevier Ltd. All rights reserved.

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.

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.

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.

Scale-dependent foraging ecology of a marine top predator modelled using passive acoustic data

1.Understanding which environmental factors drive foraging preferences is critical for the development of effective management measures, but resource use patterns may emerge from processes that occur at different spatial and temporal scales. Direct observations of foraging are also especially challenging in marine predators, but passive acoustic techniques provide opportunities to study the behaviour of echolocating species over a range of scales. 2.We used an extensive passive acoustic data set to investigate the distribution and temporal dynamics of foraging in bottlenose dolphins using the Moray Firth (Scotland, UK). Echolocation buzzes were identified with a mixture model of detected echolocation inter-click intervals and used as a proxy of foraging activity. A robust modelling approach accounting for autocorrelation in the data was then used to evaluate which environmental factors were associated with the observed dynamics at two different spatial and temporal scales. 3.At a broad scale, foraging varied seasonally and was also affected by seabed slope and shelf-sea fronts. At a finer scale, we identified variation in seasonal use and local interactions with tidal processes. Foraging was best predicted at a daily scale, accounting for site specificity in the shape of the estimated relationships. 4.This study demonstrates how passive acoustic data can be used to understand foraging ecology in echolocating species and provides a robust analytical procedure for describing spatio-temporal patterns. Associations between foraging and environmental characteristics varied according to spatial and temporal scale, highlighting the need for a multi-scale approach. Our results indicate that dolphins respond to coarser scale temporal dynamics, but have a detailed understanding of finer-scale spatial distribution of resources.

Identifying predictable foraging habitats for a wide-ranging marine predator using ensemble ecological niche models

Aim: Ecological niche modelling can provide valuable insight into species' environmental preferences and aid the identification of key habitats for populations of conservation concern. Here, we integrate biologging, satellite remote-sensing and ensemble ecological niche models (EENMs) to identify predictable foraging habitats for a globally important population of the grey-headed albatross (GHA) Thalassarche chrysostoma. Location: Bird Island, South Georgia; Southern Atlantic Ocean. Methods: GPS and geolocation-immersion loggers were used to track at-sea movements and activity patterns of GHA over two breeding seasons (n = 55; brood-guard). Immersion frequency (landings per 10-min interval) was used to define foraging events. EENM combining Generalized Additive Models (GAM), MaxEnt, Random Forest (RF) and Boosted Regression Trees (BRT) identified the biophysical conditions characterizing the locations of foraging events, using time-matched oceanographic predictors (Sea Surface Temperature, SST; chlorophyll a, chl-a; thermal front frequency, TFreq; depth). Model performance was assessed through iterative cross-validation and extrapolative performance through cross-validation among years. Results: Predictable foraging habitats identified by EENM spanned neritic (<500 m), shelf break and oceanic waters, coinciding with a set of persistent biophysical conditions characterized by particular thermal ranges (3–8 °C, 12–13 °C), elevated primary productivity (chl-a > 0.5 mg m−3) and frequent manifestation of mesoscale thermal fronts. Our results confirm previous indications that GHA exploit enhanced foraging opportunities associated with frontal systems and objectively identify the APFZ as a region of high foraging habitat suitability. Moreover, at the spatial and temporal scales investigated here, the performance of multi-model ensembles was superior to that of single-algorithm models, and cross-validation among years indicated reasonable extrapolative performance. Main conclusions: EENM techniques are useful for integrating the predictions of several single-algorithm models, reducing potential bias and increasing confidence in predictions. Our analysis highlights the value of EENM for use with movement data in identifying at-sea habitats of wide-ranging marine predators, with clear implications for conservation and management.

Identifying predictable foraging habitats for a wide-ranging marine predator using ensemble ecological niche models

Aim: Ecological niche modelling can provide valuable insight into species' environmental preferences and aid the identification of key habitats for populations of conservation concern. Here, we integrate biologging, satellite remote-sensing and ensemble ecological niche models (EENMs) to identify predictable foraging habitats for a globally important population of the grey-headed albatross (GHA) Thalassarche chrysostoma. Location: Bird Island, South Georgia; Southern Atlantic Ocean. Methods: GPS and geolocation-immersion loggers were used to track at-sea movements and activity patterns of GHA over two breeding seasons (n = 55; brood-guard). Immersion frequency (landings per 10-min interval) was used to define foraging events. EENM combining Generalized Additive Models (GAM), MaxEnt, Random Forest (RF) and Boosted Regression Trees (BRT) identified the biophysical conditions characterizing the locations of foraging events, using time-matched oceanographic predictors (Sea Surface Temperature, SST; chlorophyll a, chl-a; thermal front frequency, TFreq; depth). Model performance was assessed through iterative cross-validation and extrapolative performance through cross-validation among years. Results: Predictable foraging habitats identified by EENM spanned neritic (<500 m), shelf break and oceanic waters, coinciding with a set of persistent biophysical conditions characterized by particular thermal ranges (3–8 °C, 12–13 °C), elevated primary productivity (chl-a > 0.5 mg m−3) and frequent manifestation of mesoscale thermal fronts. Our results confirm previous indications that GHA exploit enhanced foraging opportunities associated with frontal systems and objectively identify the APFZ as a region of high foraging habitat suitability. Moreover, at the spatial and temporal scales investigated here, the performance of multi-model ensembles was superior to that of single-algorithm models, and cross-validation among years indicated reasonable extrapolative performance. Main conclusions: EENM techniques are useful for integrating the predictions of several single-algorithm models, reducing potential bias and increasing confidence in predictions. Our analysis highlights the value of EENM for use with movement data in identifying at-sea habitats of wide-ranging marine predators, with clear implications for conservation and management.

Feeding rates and prey selectivity of planktonic decapod larvae in the Western English Channel

Meroplankton are seasonally important contributors to the zooplankton, particularly at inshore sites, yet their feeding ecology is poorly known relative to holoplankton. While several studies have measured feeding in decapod larvae, few studies have examined the feeding rates of decapod larvae on natural prey assemblages throughout the reproductive season. We conducted 8 feeding experiments with Necora puber, Liocarcinus spp. and Upogebia spp. zoea larvae collected from the L4 monitoring site off Plymouth (50°15.00′N, 4°13.02′W) during spring–summer 2009 and 2010. This period spanned moderate-to-high food availability (0.5–1.6 µg chl-a L−1), but a great range in food composition with small cells <20 µm dominating in 2010. Daily rations averaged 17, 60 and 22 % of body C for the 3 respective decapod species. Clearance rates differed according to prey type, and all 3 decapod genera showed evidence of selection of dinoflagellates. Importantly, small cells including nano- and pico-plankton were ingested, this being demonstrated independently by flow cytometric analysis of the feeding experiments and molecular analysis. PCR-based analysis of the haptophyte portion of the diet revealed ingestion of Isochrysis galbana by decapod larvae in the bottle incubations and Isochrysis galbana and Phaeocystis globosa by decapod larvae collected directly from the field. This study has shown that pico- and nano-sized plankton form an important supplement to the diverse and variable diet of decapod larvae.