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.

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.

Carbon fluxes and pelagic ecosystem dynamics near two western Antarctic Peninsula Adélie penguin colonies: an inverse model approach

An inverse food-web model for the western Antarctic Peninsula (WAP) pelagic food web was constrained with data from Palmer Long Term Ecological Research (PAL-LTER) project annual austral summer sampling cruises. Model solutions were generated for 2 regions with Adelie penguin Pygoscelis adeliae colonies presenting different population trends (a northern and a southern colony) for a 12 yr period (1995-2006). Counter to the standard paradigm, comparisons of carbon flow through bacteria, microzooplankton, and krill showed that the diatom-krill-top predator food chain is not the dominant pathway for organic carbon exchanges. The food web is more complex, including significant contributions by microzooplankton and the microbial loop. Using both inverse model results and network indices, it appears that in the northern WAP the food web is dominated by the microbial food web, with a temporal trend toward its increasing importance. The dominant pathway for the southern WAP food web varies from year to year, with no detectable temporal trend toward dominance of microzooplankton versus krill. In addition, sensitivity analyses indicated that the northern colony of Adelie penguins, whose population size has been declining over the past 35 yr, appears to have sufficient krill during summer to sustain its basic metabolic needs and rear chicks, suggesting the importance of other processes in regulating the Adelie population decline.

DNA from historical and trophy samples provides insights into white shark population origins and genetic diversity

Characterizing genetic variation by retrospective genotyping of trophy or historical artifacts from endangered species is an important conservation tool. Loss of genetic diversity in top predators such as the white shark Carcharodon carcharias remains an issue, exacerbated in this species by declining, sometimes isolated philopatric populations. We successfully sequenced mitochondrial DNA (mtDNA) D-loop from osteodentine of contemporary South African white shark teeth (from 3 jaws), and from 34 to 129 yr old dried cartilage and skin samples from 1 Pacific Ocean and 5 Mediterranean sharks. Osteodentine-derived sequences from South African fish matched those derived from an individual’s finclips, but were generally of poorer quality than those from skin and cartilage of historical samples. Three haplotypes were identified from historical Mediterranean samples (n = 5); 2 individuals had unique sequences and 3 shared the contemporary Mediterranean haplotype. Placement of previously undescribed mtDNA haplotypes from historical material within both the Mediterranean and Pacific clades fits with the accepted intra-specific phylogeny derived from contemporary material, verifying our approaches. The utility of our methodology is in its provision of additional genetic resources from osteodentine (for species lacking tooth pulp) and cartilage of rare and endangered species held in often uncurated, contemporary and historical dry collections. Such material can usefully supplement estimates of connectivity, population history, and stock viability. We confirm the depauperate haplotype diversity of historical Mediterranean sharks, consistent with founding by a small number of Pacific colonizers. The consequent lack of diversity suggests serious challenges for the maintenance of this top predator and the Mediterranean ecosystem.

DNA from historical and trophy samples provides insights into white shark population origins and genetic diversity

Characterizing genetic variation by retrospective genotyping of trophy or historical artifacts from endangered species is an important conservation tool. Loss of genetic diversity in top predators such as the white shark Carcharodon carcharias remains an issue, exacerbated in this species by declining, sometimes isolated philopatric populations. We successfully sequenced mitochondrial DNA (mtDNA) D-loop from osteodentine of contemporary South African white shark teeth (from 3 jaws), and from 34 to 129 yr old dried cartilage and skin samples from 1 Pacific Ocean and 5 Mediterranean sharks. Osteodentine-derived sequences from South African fish matched those derived from an individual’s finclips, but were generally of poorer quality than those from skin and cartilage of historical samples. Three haplotypes were identified from historical Mediterranean samples (n = 5); 2 individuals had unique sequences and 3 shared the contemporary Mediterranean haplotype. Placement of previously undescribed mtDNA haplotypes from historical material within both the Mediterranean and Pacific clades fits with the accepted intra-specific phylogeny derived from contemporary material, verifying our approaches. The utility of our methodology is in its provision of additional genetic resources from osteodentine (for species lacking tooth pulp) and cartilage of rare and endangered species held in often uncurated, contemporary and historical dry collections. Such material can usefully supplement estimates of connectivity, population history, and stock viability. We confirm the depauperate haplotype diversity of historical Mediterranean sharks, consistent with founding by a small number of Pacific colonizers. The consequent lack of diversity suggests serious challenges for the maintenance of this top predator and the Mediterranean ecosystem.

Encounter success of free-ranging marine predator movements across a dynamic prey landscape

Movements of wide-ranging top predators can now be studied effectively using satellite and archival telemetry. However, the motivations underlying movements remain difficult to determine because trajectories are seldom related to key biological gradients, such as changing prey distributions. Here, we use a dynamic prey landscape of zooplankton biomass in the north-east Atlantic Ocean to examine active habitat selection in the plankton-feeding basking shark Cetorhinus maximus. The relative success of shark searches across this landscape was examined by comparing prey biomass encountered by sharks with encounters by random-walk simulations of ‘model’ sharks. Movements of transmitter-tagged sharks monitored for 964 days (16754km estimated minimum distance) were concentrated on the European continental shelf in areas characterized by high seasonal productivity and complex prey distributions. We show movements by adult and sub-adult sharks yielded consistently higher prey encounter rates than 90% of random-walk simulations. Behavioural patterns were consistent with basking sharks using search tactics structured across multiple scales to exploit the richest prey areas available in preferred habitats. Simple behavioural rules based on learned responses to previously encountered prey distributions may explain the high performances. This study highlights how dynamic prey landscapes enable active habitat selection in large predators to be investigated from a trophic perspective, an approach that may inform conservation by identifying critical habitat of vulnerable species.

From plankton to top predators: bottom-up control of a marine food web across four trophic levels

1. Abundant mid-trophic pelagic fish often play a central role in marine ecosystems, both as links between zooplankton and top predators and as important fishery targets. In the North Sea, the lesser sandeel occupies this position, being the main prey of many bird, mammal and fish predators and the target of a major industrial fishery. However, since 2003, sandeel landings have decreased by > 50%, and many sandeel-dependent seabirds experienced breeding failures in 2004. 2. Despite the major economic implications, current understanding of the regulation of key constituents of this ecosystem is poor. Sandeel abundance may be regulated 'bottom-up' by food abundance, often thought to be under climatic control, or 'top-down' by natural or fishery predation. We tested predictions from these two hypotheses by combining unique long-term data sets (1973–2003) on seabird breeding productivity from the Isle of May, SE Scotland, and plankton and fish larvae from the Continuous Plankton Recorder survey. We also tested whether seabird breeding productivity was more tightly linked to sandeel biomass or quality (size) of individual fish. 3. The biomass of larval sandeels increased two- to threefold over the study period and was positively associated with proxies of the abundance of their plankton prey. Breeding productivity of four seabirds bringing multiple prey items to their offspring was positively related to sandeel larval biomass with a 1-year lag, indicating dependence on 1-year-old fish, but in one species bringing individual fish it was strongly associated with the size of adult sandeels. 4. These links are consistent with bottom-up ecosystem regulation and, with evidence from previous studies, indicate how climate-driven changes in plankton communities can affect top predators and potentially human fisheries through the dynamics of key mid-trophic fish. However, the failing recruitment to adult sandeel stocks and the exceptionally low seabird breeding productivity in 2004 were not associated with low sandeel larval biomass in 2003, so other mechanisms (e.g. predation, lack of suitable food after metamorphosis) must have been important in this case. Understanding ecosystem regulation is extremely important for predicting the fate of keystone species, such as sandeels, and their predators.

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.

Scaling laws of ambush predator 'waiting' behaviour are tuned to a common ecology

The decisions animals make about how long to wait between activities can determine the success of diverse behaviours such as foraging, group formation or risk avoidance. Remarkably, for diverse animal species, including humans, spontaneous patterns of waiting times show random ‘burstiness’ that appears scale-invariant across a broad set of scales. However, a general theory linking this phenomenon across the animal kingdom currently lacks an ecological basis. Here, we demonstrate from tracking the activities of 15 sympatric predator species (cephalopods, sharks, skates and teleosts) under natural and controlled conditions that bursty waiting times are an intrinsic spontaneous behaviour well approximated by heavy-tailed (power-law) models over data ranges up to four orders of magnitude. Scaling exponents quantifying ratios of frequent short to rare very long waits are species-specific, being determined by traits such as foraging mode (active versus ambush predation), body size and prey preference. A stochastic–deterministic decision model reproduced the empirical waiting time scaling and species-specific exponents, indicating that apparently complex scaling can emerge from simple decisions. Results indicate temporal power-law scaling is a behavioural ‘rule of thumb’ that is tuned to species’ ecological traits, implying a common pattern may have naturally evolved that optimizes move–wait decisions in less predictable natural environments.

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.