Environmental context explains Lévy and Brownian movement patterns of marine predators

Journal name:
Nature
Volume:
465,
Pages:
1066–1069
Date published:
(24 June 2010)
DOI:
doi:10.1038/nature09116
Received
Accepted
Published online

An optimal search theory, the so-called Lévy-flight foraging hypothesis1, predicts that predators should adopt search strategies known as Lévy flights where prey is sparse and distributed unpredictably, but that Brownian movement is sufficiently efficient for locating abundant prey2, 3, 4. Empirical studies have generated controversy because the accuracy of statistical methods that have been used to identify Lévy behaviour has recently been questioned5, 6. Consequently, whether foragers exhibit Lévy flights in the wild remains unclear. Crucially, moreover, it has not been tested whether observed movement patterns across natural landscapes having different expected resource distributions conform to the theory’s central predictions. Here we use maximum-likelihood methods to test for Lévy patterns in relation to environmental gradients in the largest animal movement data set assembled for this purpose. Strong support was found for Lévy search patterns across 14 species of open-ocean predatory fish (sharks, tuna, billfish and ocean sunfish), with some individuals switching between Lévy and Brownian movement as they traversed different habitat types. We tested the spatial occurrence of these two principal patterns and found Lévy behaviour to be associated with less productive waters (sparser prey) and Brownian movements to be associated with productive shelf or convergence-front habitats (abundant prey). These results are consistent with the Lévy-flight foraging hypothesis1, 7, supporting the contention8, 9 that organism search strategies naturally evolved in such a way that they exploit optimal Lévy patterns.

Figures at a glance

  1. Figure 1: Examples of good fits to power-law and truncated power-law distributions.
    Examples of good fits to power-law and truncated power-law distributions.

    a, Synthetic power-law and truncated power-law (Pareto) distributions with upper truncations set to 50, 250, 5,000. bf, Empirical power-law and truncated power-law fits to dive data from individual blue sharks (Prionace glauca; b, d) and an ocean sunfish (Mola mola, e), together with the diving time series for the individual in b (over ~8d; c) and the individual in e (over ~4d; f). The red line indicates a synthetic power law in a, a power law in b and truncated power-law MLE model fits to empirical data in d and e.

  2. Figure 2: Behavioural switching between Lévy and Brownian motion in relation to habitat type.
    Behavioural switching between Lévy and Brownian motion in relation to habitat type.

    a–e, Split moving-window analysis showing significant discontinuities in the dive time series of blue shark 10. Red lines indicate points where the time series was divided into sections (SEC1–SEC5). f–j, MLE analysis with μ values for sections best fitting a truncated power-law distribution: black circles, observed step lengths; red lines, best-fit truncated power law; blue lines, best-fit exponential distribution. k–o, Depth profiles of sea temperature recorded using electronic tags. p, q, Geo-referenced track sections of blue shark 10 overlaid on chlorophyll a concentrations (p) and bathymetry (q). Section numbers correspond to those in a–e and different data-point colours correspond to different sections: SEC1, black (higher latitude); SEC2, white (higher latitude); SEC3, grey; SEC4, black (lower latitude); SEC5, white (lower latitude).

  3. Figure 3: Spatial occurrence of Lévy and Brownian behaviour types.
    Spatial occurrence of Lévy and Brownian behaviour types.

    Frequencies of behaviour types in productive (frontal/shelf) and less productive (off-shelf) habitats in the northeast Atlantic (a), and in productive (frontal) and less productive (stratified) habitats in the central eastern Pacific (b). Tests of two predictions of the LFF hypothesis (Lévy behaviour where prey is sparse; Brownian movement where prey is abundant and not sparsely distributed) were performed on frequency data (not per cent frequency data). See main text for details of the statistical tests.

Author information

Affiliations

  1. Marine Biological Association of the United Kingdom, The Laboratory, Citadel Hill, Plymouth PL1 2PB, UK

    • Nicolas E. Humphries,
    • Nuno Queiroz,
    • Jennifer R. M. Dyer,
    • Nicolas G. Pade,
    • Victoria J. Wearmouth,
    • Emily J. Southall &
    • David W. Sims
  2. Marine Biology and Ecology Research Centre, Marine Institute, School of Marine Sciences and Engineering, University of Plymouth, Drake Circus, Plymouth PL4 8AA, UK

    • Nicolas E. Humphries &
    • David W. Sims
  3. CIBIO – Universidade do Porto, Campus Agrário de Vairão, Rua Padre Armando Quintas, 4485-668 Vairão, Portugal

    • Nuno Queiroz
  4. Institute of Biological and Environmental Sciences, School of Biological Sciences, University of Aberdeen, Tillydrone Avenue, Aberdeen AB24 2TZ, UK

    • Nuno Queiroz,
    • Nicolas G. Pade,
    • Catherine S. Jones &
    • Leslie R. Noble
  5. Joint Institute for Marine and Atmospheric Research, University of Hawaii at Manoa, Kewalo Research Facility/NOAA Fisheries, 1125-B Ala Mona Boulevard, Honolulu, Hawaii 96814, USA

    • Michael K. Musyl
  6. Inter-American Tropical Tuna Commission, 8604 La Jolla Shores Drive, La Jolla, California 92037-1508, USA

    • Kurt M. Schaefer &
    • Daniel W. Fuller
  7. ETH Zurich, Raemistrasse 101, CH-8092 Zurich, Switzerland

    • Juerg M. Brunnschweiler
  8. Coastal and Marine Resources Centre, ERI, University College Cork, Glucksman Marine Facility, Naval Base, Haulbowline, Cobh, Cork, Ireland

    • Thomas K. Doyle
  9. School of Biological Sciences, Queen’s University Belfast, Medical Biology Centre, 97 Lisburn Road, Belfast BT9 7BL, UK

    • Jonathan D. R. Houghton
  10. Department of Pure and Applied Ecology, Institute of Environmental Sustainability, Swansea University, Singleton Park, Swansea SA2 8PP, UK

    • Graeme C. Hays

Contributions

D.W.S. designed the study. N.E.H. and D.W.S. completed data analysis with contributions from N.Q. and J.R.M.D. N.E.H. designed and developed the software for MLE and split moving-window analyses. D.W.S. and N.E.H. wrote the paper and all authors contributed to subsequent drafts. Field data were collected by D.W.S., E.J.S., N.Q., N.G.P., M.K.M., K.M.S., D.W.F., J.M.B., T.K.D., J.D.R.H., G.C.H. and V.J.W.

Competing financial interests

The authors declare no competing financial interests.

Corresponding author

Correspondence to:

Author Details

Supplementary information

PDF files

  1. Supplementary Information (5.3M)

    This file contains Supplementary Methods, Supplementary Results, Supplementary Tables S1-S4, Supplementary Figures S1-S9 with legends and References.

Additional data