- Open Access
Activity seascapes highlight central place foraging strategies in marine predators that never stop swimming
© The Author(s). 2018
- Received: 16 February 2018
- Accepted: 28 May 2018
- Published: 21 June 2018
Central place foragers (CPF) rest within a central place, and theory predicts that distance of patches from this central place sets the outer limits of the foraging arena. Many marine ectothermic predators behave like CPF animals, but never stop swimming, suggesting that predators will incur ‘travelling’ costs while resting. Currently, it is unknown how these CPF predators behave or how modulation of behavior contributes to daily energy budgets. We combine acoustic telemetry, multi-sensor loggers, and hidden Markov models (HMMs) to generate ‘activity seascapes’, which combine space use with patterns of activity, for reef sharks (blacktip reef and grey reef sharks) at an unfished Pacific atoll.
Sharks of both species occupied a central place during the day within deeper, cooler water where they were less active, and became more active over a larger area at night in shallower water. However, video cameras on two grey reef sharks revealed foraging attempts/success occurring throughout the day, and that multiple sharks were refuging in common areas. A simple bioenergetics model for grey reef sharks predicted that diel changes in energy expenditure are primarily driven by changes in swim speed and not body temperature.
We provide a new method for simultaneously visualizing diel space use and behavior in marine predators, which does not require the simultaneous measure of both from each animal. We show that blacktip and grey reef sharks behave as CPFs, with diel changes in activity, horizontal and vertical space use. However, aspects of their foraging behavior may differ from other predictions of traditional CPF models. In particular, for species that never stop swimming, patch foraging times may be unrelated to patch travel distance.
- Hidden Markov models
- Coral reefs
Central place foraging (CPF) is a ubiquitous behavior seen across animal groups ranging from insects, to birds, and humans . Unlike random movements within a home range, CPF behavior consists of periodic and predictable movements to and from a central place, often with multiple individuals sharing the central place . CPF animals tend to rest at the central place, with their energy costs increasing as they travel greater distances from this location [1, 2]. As such, the costs associated with travel distance to the patch should define the limits of the animals foraging range from the central place [2, 3]. CPF behavior can lead to heterogeneity in habitat or prey distribution as the animal’s foraging rates will likely vary with distance from the central place [4, 5]. As patch distance to the central place increases, travel costs also increase, and the animal should spend more time foraging at the patch [1, 4].
A key assumption of CPF theory is that an animal rests at the central place, and foraging costs increase with travel distance to a feeding patch. Yet, there are many species of marine predators that exhibit CPF-like behavior, but never stop swimming, and never truly rest. For these animals, energy costs may be independent of travel distance to the patch, and simply a function of swim speeds. Such predators include large coral reef-associated fishes (sharks and teleosts), which swim within a relatively small core area during the day, and move over an expanded area at night, with periodicity of movements related to diel or tidal cycles [6–9]. Marine animals also move in a three-dimensional (3D) environment and CPF behavior can include a vertical component as well as a horizontal one, with individuals performing diel vertical migrations (DVM) relative to the central place [10–12]. While few studies have measured actual activity and swim speeds, some tropical reef sharks display patterns of activity that also vary with diel and tidal cycles [13–17]. However, why predators that do not use a shelter or ever stop swimming require a central place, is unclear.
CPF animals that never stop swimming are almost exclusively ectotherms, so metabolic rates are sensitive to changes in ambient temperature. At any particular time, routine metabolic rates should be a function of body temperature, movement speed, and other aspects of the movement process (e.g. turning costs more than straight line swimming ). While the animal may not stop swimming, they can establish the central place in cooler waters where metabolic rates are reduced. If the animal simultaneously maintains low activity in the central place while cooler, then energy costs may be essentially similar to ‘resting’ . In addition to changes in body temperature, routine metabolic rates can be modulated via changes in swim speed. Hence, the energetics of CPF in these animals must consider body temperature and movement rates.
Accelerometers have become popular for measuring both the activity and energy requirements of free-ranging marine animals [14, 19–21]. Accelerometer data can be combined with information about the geographic location of the animal to generate a spatial representation of the animal’s energy costs [19–21]. However, it is difficult to separate areas of high-energy expenditure (e.g. traversing through an energetically expensive habitat) from areas of high animal activity related to specific behaviors (e.g. foraging) within the landscape, especially for animals whose behaviors cannot be easily defined from sensor measurements (e.g. continuously swimming fish). For our general case of CPF foragers, we are interested in how the predator’s foraging activity varies spatially in relation to the central place. This challenge is further complicated in fishes because space use and activity are often measured at different temporal scales, due to limitations of tagging technology. Movements can be recorded via telemetry over multiple years, while fine-scale behaviors from accelerometers are recorded over time-frames of days to weeks. The diel behavior of CPF animals is likely to be predictable, suggesting that measurements made over shorter time periods are representative of long-term behavior . We took advantage of this predictability to develop a new spatial representation of activity we term an ‘activity seascape’, for marine predators which show CPF behavior. The activity seascape combines long-term patterns of space use derived from acoustic telemetry data with the diel probability of being more or less active, which is based on statistical time-series models inferred from short deployment accelerometers . The seascapes allow the locations/times where high activity may be strongest/weakest to be identified without requiring space use and activity of individuals to be measured simultaneously.
We use multi-sensor loggers and statistical time-series models to build activity seascapes for two species of reef shark (blacktip reef sharks (Carcharhinus melanopterus), and grey reef sharks (C. amblyrhynchos)), whose movements appear similar to a CPF as they use smaller areas during the day than at night [7–9, 12, 14, 23]. We also used these activity seascapes to better understand the dynamics of CPF in animals that never stop swimming. Specifically, we predict that sharks use a central place during the day while in deeper and cooler water where they decrease activity and reduce energy costs. We predict they move over a larger area at night into shallower water, where they become more active. We then build a simple bioenergetics model for grey reef sharks to predict the relative contribution of body temperature and swim speed to daily energy budget. We further predict that grey reef sharks will maintain lower body temperatures within their central place as a tactic to reduce ‘travel costs’ within their refuge.
Research was conducted at Palmyra Atoll (5°54’ N; 162°05’ W), located in the central Pacific Ocean. Palmyra has been a U.S. federal wildlife refuge since 2001, with only a research station housing a small number of scientists and staff. The atoll consists of two lagoons, surrounded by sandflats and backreefs, which transition to forereef habitats (see Additional file 1: Appendix S1). Due to its protected status, large shark populations are found at the atoll . Further details of the site can be found in [9, 14].
Patterns of activity
Details of sharks tagged with data-loggers or acoustic transmitters. BLT (blacktip reef shark), GR (grey reef shark), TL (Total Length)
TL range (cm)
TL range (cm)
2) Acoustic transmitters
Analysis of shark activity
We filtered the static contribution of gravity from raw acceleration data, and then calculated overall dynamic body acceleration (ODBA, ). ODBA was used as a measure of activity as it incorporates tail beat frequency as well as activity along other body axes. However, being time series data, ODBA values are highly correlated such that inferences based on statistical models that do not take this key feature into account will usually be invalid. Furthermore, the autocorrelation structure itself will provide interesting behavioral information as the probability of the animal being active is likely to be a function of how active it was previously. Understanding the biological importance of changes in ODBA in animals that swim continuously is difficult as there is no true ‘resting’ period. We were interested in periods of increased activity, particularly ‘bursts’ which could be indicative of foraging behavior.
Hidden Markov models (HMMs) are stochastic time series models where observed data (e.g. travel speed, depth, ODBA) are assumed to be driven by an underlying hidden process. We assume that the hidden process can be in either of N = 2 states, roughly corresponding to behavioral states, which we label ‘relatively high activity’ and ‘relatively low activity’, respectively. The observations can be considered noisy measurements of the behavioral state, which cannot be explicitly observed [22, 26]. Traditionally, HMMs have inferred behavioral states from the movement process itself (e.g. rate of movement, turning angles), but more recently have been used with behavioral data and ODBA specifically . We developed a 2-state HMM, where sharks were in either a low state of activity (state 1) or high state of activity (State 2), based on ODBA data from accelerometers. We could then compute probabilities of sharks changing or remaining within behavioral states based on time of day, tidal state, or swimming depth and water temperature. ODBA values were averaged over 1 s intervals before applying the HMMs, based on observations of behaviors from video footage (see below). We removed the first 4 h of data from each animal so that we did not infer behavior when the animals may still have been highly stressed. All HMMs were built in the statistical environment R using customized code. All HMM details can be found in Additional file 2: Appendix S2.
The HMMs provide a data-driven, objective approach to analyzing acceleration data, but we still cannot identify which specific behavior (e.g. feeding, predator avoidance) sharks perform while within the various states. We had the unique opportunity to correlate the results of the HMM with the simultaneously collected 22 h of video obtained from two grey reef sharks. Video data allowed us to see what the sharks were doing (and in which habitats) when inferred by the HMM to be in certain states, providing validation for our interpretation of the model. We first observed all 22 h of video noting times of foraging, increased activity (both of the individual shark and conspecifics) as well as the presence of other sharks. For those specific time periods, we could then compare observations with the behavioral state predicted by the HMM.
We also analyzed swim speed data in grey reef sharks as a direct measure of energy expenditure. To explore diel changes in swim speed, we constructed generalized additive models (GAMs) for each individual grey reef shark (N = 3) using hourly mean swim speed data. Speed will still likely suffer from serial correlation, which we accounted for by including an AR(1) (first-order auto-regressive) process with time as the position variable. The correlation at lag = 1 was included in the model to specify the correlation structure. The GAM was constructed with a Gaussian error distribution and time of day was modeled with a cyclic smooth spline. Model fit was assessed by examining residual diagnostic plots, and Akaike’s information criterion  (AIC) was used to assess model performance against a null model (intercept only), with improved model fit indicated by a minimum ΔAIC value > 3 . GAM analyses were conducted in R using the mgcv package.
Patterns of space use
We quantified reef shark space use patterns using acoustic telemetry methods. Between 2010 and 2014, sharks were caught on hook and line and surgically implanted with an acoustic transmitter (V16, 69 kHz, Vemco ltd, Nova Scotia). A small incision was made on the ventral surface, the transmitter was inserted into the body cavity and a single suture was used to close the incision. The sharks were measured, sexed, externally dart tagged and released. These individuals were not the same sharks as those tagged with data-loggers above, although there was temporal overlap between periods of data collection (i.e. individuals were being tracked during the same time when sharks were carrying loggers). We acoustically tagged a total of 20 blacktip reef sharks, and 43 grey reef sharks. Of those, 9 blacktip reef sharks and 13 grey reef sharks were tagged with V16 PT transmitters that also measured pressure (depth) and body temperature, and 6 blacktip reef sharks were tagged with V13 AP transmitters that measured 3D acceleration in addition to depth. Transmitters were detected by an array of up to 70 subsurface omni-directional acoustic receivers (VR2W) deployed throughout the atoll . Every time a shark swam within range of the receiver, the date/time of detection, along with swimming depth and body temperature (for those individuals with PT tags) were recorded. We downloaded receivers annually. Receiver detection range can vary by habitat and range tests of a subsample of receivers showed an approximate range of 250 m on backreefs and close to 500 m on forereef habitats. Note that blacktip reef shark acceleration, depth, and body temperatures were analyzed in .
We calculated shark spatial utilization distributions (UD) using a Brownian bridge movement model (BBMM), where consecutive acoustic detections between receivers by moving sharks are linked by conditional random walks [29, 30]. The BBMM also incorporates measurement error (we set this to 300 m as the average receiver detection range) to provide a more realistic representation of the possible space used by the animal . CPF behavior will consist of periodic excursions from the central place, and we were interested in the distances sharks moved throughout the diel period from this core location. We calculated the central place for each individual separately. For blacktip reef sharks, we defined the central place as the 50% UD predicted by the BBMM. However, grey reef sharks were detected on far fewer receivers, making it difficult to calculate 2-dimensional BBMMs. As such, for each individual we defined the core area as the receiver where 80% of detections were made, with a 1000 m buffer (assuming a detection range of 500 m on either side of forereef receivers). For both species of shark, we then quantified the distance of movement from this core area throughout the diel cycle. That is, for each individual we created a time series of distances to its own central place and then in the last step we averaged these distances into a time series of average displacement for each shark species. Analyses were performed in R using the adehabitat LT and adehabitat HR packages.
To explore diel patterns of vertical habitat use by grey reef sharks, we constructed generalized additive mixed models (GAMMs) using hourly mean depth and body temperature data from the acoustic transmitters. Specifically, we estimated the effect of time of day on swimming depth and temperature, with individual shark modeled as a random effect. All other model components were the same as described above for swim speed. Diel changes in depth and temperature for blacktips were reported in .
Patterns of space use and activity seascapes
The acceleration data and HMMs provided probabilities of sharks being in a relatively high active state throughout the diel cycle. We then combined these with three-dimensional UDs calculated over months to years of movement from the telemetry data, to generate an overall visual representation of the ‘activity seascape’ of individual animals within their environment. We used the telemetry data to produce space-time cubes (STC) where the bottom two dimensions represent the geographic space over which the movement occurs and the third dimension is time . In a 2D kernel density UD, the surface representing this density is divided into square grid cells or pixels, the value of which represents the probability of movement in that particular location. Analogously, in a 3D UD, the space-time volume is divided into so-called voxels, i.e. three-dimensional grid cubes, with each voxel assigned a probability of movement in that particular location in space and time. Our three-dimensional BBMM generalizations allow us to use time as part of the BBMM calculation, and show time visually on the third axis of the space-time density volume . We built 3D UDs for each shark by aggregating telemetry data by day. The space-time UDs were then combined with the diel activity probabilities from HMMs in order to visually emphasize areas in the space-time UDs when sharks were the most active. We did this by multiplying each voxel in the space-time UD volume with the probability at that particular moment in time of the shark being in an active state, as determined by the HMMs. We removed any individuals from analysis that had < 100 detections. The activity seascape algorithm is in the process of being published as R package. In the mean time, a preliminary version of the R code will be placed to https://github.com/udemsar upon publication.
To predict the relative contribution of swimming speed and body temperature on daily energy budgets, we built a simple bioenergetics model for grey reef sharks. We calculated shark mass using W = 0.0045 L3.21 where L is total length in cm and W is weight in kg (http://fishbase.org). Routine metabolic rates (M) were then estimated using the equation for ectothermic sharks, Log10M = 0.79*Log10W + 2.31 . Metabolic rates were corrected for changes in body temperature assuming a Q10 of 1.65 and 3.0, to include the range of Q10 values seen within tropical sharks [32, 33]. To predict the effect of changing swim speed, we used the model developed for requiem sharks in . Briefly, we assume that the average swimming speed for each animal represents its optimal travel speed that minimizes its cost of transport. The percentage increase or decrease in swim speeds from this average leads to an equivalent change in metabolic rate (e.g. 1% increase in speed = 1% increase in active metabolic rate ). We calculated hourly changes in routine metabolic rate for a 38.3 kg (average grey reef shark body mass at Palmyra) individual using the observed diel changes in body temperature, and then assuming the animal maintained constant body temperature throughout the diel cycle (average 28.0 °C). Diel changes in swim speed were inferred from the GAM results from the swim speed sensor data, described above. Model results were compared with a paired 2 sample t-test. Our goal was not to estimate daily energy costs, but rather to determine if observed diel changes in body temperature were biologically relevant (at least regarding energy costs) and the relative contribution of shark swim speed to the daily energy budget.
Patterns of activity
Grey reef sharks (154–159 cm TL) swam at average speeds of 0.51–0.63 m/s with maximum speeds of 3.6 m/s. Although the average speeds differed between the three individuals, they all showed the same diel pattern with increased speeds at 02:00–03:00 and lowest speeds from 12:00–18:00 (Additional file 6: Appendix S6). Average swim speeds increased by 5–10% during the night (GAM AIC values = − 13,561.7.27 (A), − 10,896.4 (C), − 27,936.8 (D) compared to null models AIC values = − 8423.6 (A), − 7730.9 (C), − 16,824.0 (D), Additional file 5: Appendix S5).
Patterns of space use and activity seascapes
Shark vertical habitat use appeared to vary relative to the diel cycle and the area of core use. Our previous analysis of blacktip reef shark vertical movements and body temperature (from acoustic telemetry data), revealed very moderate diel changes, with these sharks using water 2–3 m deeper during the day . Grey reef sharks displayed more dramatic diel shifts in swimming depths, using deeper depths during the day (average 53 m) than at night (average 23 m, GAMM, AIC value = 193,034.1, null model AIC value = 206,293.0, Additional file 7: Appendix S7). Body temperatures in grey reef sharks also varied with the diel period, with lower body temperatures during the day (average 27.8 °C) when the sharks were deeper, than at night (average 28.1 °C, GAMM, AIC value = 11,059.5, null model AIC value =15,518.0, Additional file 7: Appendix S7), although there was more variability between individuals compared to diving patterns.
Due to our limited sample size of data-logger equipped sharks, and to ensure that shark movement data (from telemetry) at least overlapped temporally with the period of data-logger deployment, we only generated activity seascapes for blacktip reef sharks tagged in the backreef and grey reef sharks tagged along the SW portion of the atoll (locations where data-logger equipped sharks were caught and released). As such, activity seascape data consists of 17 blacktip reef sharks tagged on the backreefs and 17 grey reef sharks tagged within and just outside the channel (Additional file 1: Appendix S1). Activity seascapes for both species showed a high concentration of space use/activity in the daytime central place, with activity becoming spatially more diffuse at night (Figs. 6 and 7). For both species, multiple individuals used the same daytime core area, suggestive of refuging behavior (Figs. 6 and 7). As changes in diel activity were more pronounced in blacktip reef sharks, activity seascapes showed greater variation from simple space use UDs. The importance of daytime activity was greatly reduced within the central place, while the early evening use of the main channel was highlighted as a location where activity of blacktip reef sharks may be high (Fig. 6).
As predicted, both species of shark used a smaller central place in deeper water during the day when they were less active, before expanding movements into shallow water at night with a concurrent increase in activity. Blacktip reef sharks showed more dramatic diel changes in activity (but much more subtle changes in swimming depth) and were most active during the early evening, while grey reef sharks tended to be most active in the middle of the night. The early evening peak in activity seen in blacktip reef sharks was previously suggested to be a function of increased foraging success (due to falling or low light levels) and predator/prey thermal physiology . Aspects of CPF behavior clearly vary between the species, as blacktip reef sharks use the central place day and night; they just expand the total area they use at night. Grey reef sharks were rarely detected at their central place at night so were not reusing the same area throughout the diel cycle. Furthermore, there is strong spatial separation between the species at Palmyra, with blacktip reef shark activity spaces located in the backreef, lagoons, and shallow forereef, while grey reefs are situated on the deeper sections of the forereefs . Activity cycles are based on small samples sizes, although there is evidence that diel activity rhythms in tropical reef sharks fit a more general pattern [23, 36]. For example, an acoustic camera placed in the main channel at Palmyra documented a peak in number of sharks (assumed to be blacktip reef sharks based on nighttime fishing) seen at the exact time and place we predict maximal blacktip reef shark activity (early evening in the channel ).
Changes in activity were concurrent with diel vertical migrations to and from the central place. Blacktip reef sharks perform small shifts in swimming depth, but tides played a greater role in swimming depth, with individuals moving close to the surface during daytime low tide periods, where they are directly warmed by the sun [14, 38, 39]. Grey reef sharks maintained cooler body temperatures during the day by occasionally diving below the thermocline and showed a greater diel shift in depth. However, their diel variation in body temperature was small and dives below the thermocline during the day were short and tended to occur while they were in a high activity state. Due to thermal inertia effects, grey reef sharks would have warmer muscle temperatures than mesophotic teleost prey residing below the thermocline, as long as dive durations are short. The increased muscle performance from being warmer could provide sharks with a hunting advantage during these periodic deeper excursions. Although behavior below the thermocline was consistent between individual grey reef sharks, there was also clearly some variability in the depths selected by individuals, with one shark staying considerably deeper than the others. The reasons for this variability are unknown, but could reduce intra-specific competition which is also thought to drive horizontal separation between grey reef shark groups at Palmyra .
Activity seascapes and video footage suggest that some foraging may occur in the central place of grey reef sharks at Palmyra. If foraging has some significance to CPF behavior in grey reef sharks, then we would expect the location of the central place to be in an area of high diurnal prey biomass. We recently conducted a mark-recapture study to determine spatially dependent densities for grey reef sharks at Palmyra . All fishing was performed during the day so technically these distribution patterns represent daytime central place locations. We also used diver surveys to estimate spatial patterns of biomass for reef fishes of lengths 10–60 cm, as a proxy of potential prey distribution (, NOAA Coral Reef Ecosystem Program, Additional file 8: Appendix S8). Qualitatively, there was broad overlap between shark density hotspots and prey biomass, with highest levels of both metrics at the eastern and western tips and lower values along the southern and northern portions of the atoll (Additional file 8: Appendix S8). While abiotic conditions may explain these patterns (e.g. areas of high current flow), we at least provide evidence that grey reef shark central place locations are in areas where prey is most abundant. Similarly, seabirds establish central place locations dependent on local productivity and prey availability, although of course seabird prey are not located directly at the central place .
Although our results match our general predictions, our bioenergetics model for grey reef sharks suggests that body temperatures are having an insignificant effect on routine metabolic rates, at least at the diel scale. Instead, diel changes in routine metabolic rates are driven by changes in swim speed. Similarly, activity made a much larger contribution to the short-term field metabolic rates of sea snakes than changes in temperature . These results suggest that CPF behavior of grey reef sharks may involve some different mechanisms than those seen in other animals, at least for the Palmyra population (i.e. individuals at other locations may display more dramatic diel changes in body temperatures). Firstly, we cannot assume that foraging does not take place within the central place, which can have implications for patterns of intra-specific competition in refuging animals. Secondly, energy costs will not be a function of travel distance to the central place, but instead related to movement speeds and likely tortuosity . Unlike seabirds and other colonial CPF, access to foraging patches by sharks should not be a function of distance to the patch as they still have to swim within the central place; sharks simply modulate speed . This would also call into question a key prediction of traditional CPF models; in reef sharks foraging time within a patch should be independent of patch to central place distance, in species that never stop swimming . Of course, this assumes the only currency of value to the sharks is energy, when time itself may also be important, especially if coral reefs are heterogeneous with regards habitat quality. In this case, the time taken to swim to distant patches may come at lost opportunities of foraging at more reliable patches close to the central place.
Why would sharks, which do not use a nest/shelter or stop swimming, need a central place? CPF would require some degree of memory although the animal may only have to use basic path integration and memorization of simple landmarks to return to a central place, which should optimize foraging rates [42, 43]. Therefore, CPF behavior may optimize foraging success via improved navigation to foraging areas, while minimizing memory requirements (which come with costs) when prey communities are relatively abundant [42, 43]. CPF behavior may also foster social associations between reef sharks, as ‘meeting’ would require a common daytime habitat. We found that multiple individuals used the same core daytime areas, and video footage showed large numbers of conspecifics sharing the same areas and at times even foraging simultaneously. It is becoming increasingly apparent that some species of shark, including grey and blacktip reef sharks, are capable of forming social associations with conspecifics, and ‘refuging’ during the day may be a mechanism to allow such associations to form [6, 15, 44, 45].
Ideally, one would record spatial position and high resolution activity from each individual simultaneously, and over long time periods, improving our ability to identify ‘activity hotspots’ . Doing so with fishes will require technological advancements to existing satellite or acoustic telemetry techniques or the development of completely new technologies [20, 21]. While cyclical patterns of activity are likely reasonably consistent in CPF animals that show high residency, this may change throughout the year and requires testing. For example, a semi-captive study with lemon sharks suggested that the extent of diel activity differences could be a function of sex and social position . New logging technologies may also enable feeding events and meal size to be directly measured in free-ranging sharks, rather than inferred from acceleration signals . Analytically, state space models and HMMs can now incorporate multiple variables simultaneously from individual animals (e.g. diving depth, acceleration, stomach temperature, turn angles) to infer and validate behavioral states and produce population-level activity budgets [47, 48]. Future studies could also then address individual variability in behavior rather than homogenizing it within population scale analyses.
Understanding the small-scale dynamics of predation requires approaching behavior from a spatial and temporal context in concert. When considering predator effects on prey, we must consider patterns of predator habitat selection, but also predator behavior within those habitats. A recent study suggested that even trophic cascades in marine systems may be confined to spatial ‘hotspots’ used by predators (e.g. ). Our activity seascape approach simultaneously displays the space used and diel activity of marine predators and could further identify locations of ecological importance, taking into account space use and behaviour. These advancements could significantly improve our understanding of how marine predators may influence lower trophic levels and further explain the evolution of CPF behavior in marine animals.
We develop a new method for simultaneously displaying diel space use and activity, in marine animals where space use and behavior are often measured over different time scales. We use this method to show that reef sharks behave similar to central place foragers, occupying deeper water during the day where they are less active, and being more active in shallower water at night. However, unlike other CPF animals, grey reef sharks at least, showed some evidence of foraging in the central place. Furthermore, diel changes in energy expenditure appear to be regulated by changes in swim speed and not body temperature. The implications are that traditional predictions from CPF theory regarding patch foraging times, may not apply for CPF animals that never stop swimming.
We would like to thank K. Pollock, A. Pollock, H. Lynch and all the TNC staff who helped accommodate our stay at Palmyra. We thank I. Williams and the Coral Reef Ecosystem Program for reef fish biomass data. This is publication number PARC-0129 from the Palmyra Atoll Research Consortium.
This work was supported by the Marisla Foundation, the Grants-in-Aids for scientific research from the Japan Society for the Promotion of Science (25850138), the Hawaii Undersea Research Laboratory/NOAA Coral Reef Conservation Program (NA05OAR4301108, NA09OAR43002), and the Pelagic Fisheries Research Program/Joint Institute for Marine and Atmospheric Research (NA09OAR4320075).
Availability of data and materials
Data used in this study is available upon reasonable request to the lead author.
YPP, YYW, and JEC conceived the study. Statistical analysis was performed by RL, VL, and DB, while spatial analysis was performed by UD. YPP, DB, CGL, AMF, JEC, KW, and YYW performed fieldwork. All authors contributed to writing the manuscript. All authors read and approved the final manuscript.
Shark tagging was permitted under UCSB IACUC #856, and US Fish and Wildlife permit # (Permit numbers #12533–14,011, #12533–13,011, #12533–12,011, #12533–11,007, #12533–10,011, #12533–09010, #12533–08011, and #12533–07006).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hamilton WJ, Watt KEF. Refuging. Ann Rev Ecol Syst. 1970;1:263–86.View ArticleGoogle Scholar
- Olsson O, Bolin A. A model for habitat selection and species distribution derived from central place foraging theory. Oecologia. 2014;175:537–48.View ArticlePubMedGoogle Scholar
- Boyd C, Punt AE, Weimerskirch H, Bertrand S. Movement models provide insights into variation in the foraging effort of central place foragers. Ecol Model. 2014;286:13–25.View ArticleGoogle Scholar
- Kacelnik A. Central place refuging in starlings (Sturnus vulgaris). I. Patch residence time. J Anim Ecol. 1984;53:283–99.View ArticleGoogle Scholar
- Chase JM. Central-place forager effects on food web dynamics and spatial pattern in northern California meadows. Ecology. 1998;79:1236–45.View ArticleGoogle Scholar
- Klimley AP, Nelson DR. Diel movement patterns of the scalloped hammerhead shark (Sphyrna lewini) in relation to El Bajo Espiritu Santo: a refuging-central position social system. Behav Ecol Sociol. 1984;15:45–54.View ArticleGoogle Scholar
- McKibben JN, Nelson DR. Patterns of movement and grouping of gray reef sharks, Carcharhinus amblyrhynchos , at Enewetak, Marshall Islands. Bull Mar Sci. 1986;38:89–110.Google Scholar
- Field IC, Meekan MG, Speed CW, White W, Bradshaw CJA. Quantifying movement patterns for shark conservation at remote coral atolls of the Indian Ocean. Coral Reefs. 2011;30:61–71.View ArticleGoogle Scholar
- Papastamatiou YP, Friedlander AM, Caselle JE, Lowe CG. Long-term movement patterns and trophic ecology of blacktip reef sharks (Carcharhinus melanopterus) at Palmyra atoll. J Exp Mar Biol Ecol. 2010;386:94–102.View ArticleGoogle Scholar
- Sims DW, Weamouth VJ, Southall EJ, Hill JM, Moore P, Rawlinson K, Hutchinson N, Budd GC, Righton D, Metcalfe JD, Nash JP, Morritt D. Hunt warm, rest cool: bioenergetic strategy underlying diel vertical migration of a benthic shark. J Anim Ecol. 2006;75:176–90.View ArticlePubMedGoogle Scholar
- Chapman DD, Pikitch EK, Babcock EA, Shivji MS. Deep-diving and diel changes in vertical habitat use by Caribbean reef sharks, Carcharhinus perezi. Mar Ecol Prog Ser. 2007;344:271–5.View ArticleGoogle Scholar
- Heupel MR, Simpfendorfer CA. Long-term movement patterns of a coral reef predator. Coral Reefs. 2015;34:679–91.View ArticleGoogle Scholar
- Whitney NM, Papastamatiou YP, Holland KN, Lowe CG. Use of an acceleration data logger to measure diel activity patterns in captive whitetip reef sharks, Triaenodon obesus. Aquat Liv Res. 2007;20:299–305.View ArticleGoogle Scholar
- Papastamatiou YP, Watanabe YY, Bradley D, Dee LE, Weng K, Lowe CG, Caselle JE. Drivers of daily routines in an ectothermic marine predator: hunt warm, rest warmer? PLoS One. 2015;10:e0127807.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson ADM, Brownscombe JW, Krause J, Krause S, Gutowsky LFG, Brooks EJ, Cooke SJ. Integrating network analysis, sensor tags, and observations to understand shark ecology and behavior. Behav Ecol. 2015;26:1577–86.View ArticleGoogle Scholar
- Barnett A, Payne NL, Semmens JM, Fitzpatrick R. Ecotourism increases the field metabolic rate of whitetip reef sharks. Biol Conserv. 2016;199:132–6.View ArticleGoogle Scholar
- Shipley ON, Brownscombe JW, Danylchuk AJ, Cooke SJ, O’Shea OR, Brooks EK. Fine-scale movement and activity patterns of Caribbean reef sharks (Carcharhinus perezi) in the Bahamas. Environ Biol Fish. 2018; In PressGoogle Scholar
- Wilson RP, Griffiths IW, Legg PA, Friswell MI, Bidder OR, Halsey LG, Lambertucci SA, Shepard ELC. Turn costs change the value of animal search paths. Ecol Lett. 2013;16:1145–50.View ArticlePubMedGoogle Scholar
- Wilson RP, Quintana F, Hobson VJ. Construction of energy landscapes can clarify the movement and distribution of foraging animals. Proc Roy Soc B. 2011; https://doi.org/10.1098/rspb.2011.1544.
- Brownscombe JW, Cooke SJ, Danylchuk AJ. Spatiotemporal drivers of energy expenditure in a coastal marine fish. Oecologia. 2017;183:689–99.View ArticlePubMedGoogle Scholar
- Udyawer V, Simpfendorfer CA, Heupel MR, Clark TD. Temporal and spatial activity-associated energy partitioning in free-swimming sea snakes. Funct Ecol. 2017;31:1739–49.View ArticleGoogle Scholar
- Leos-Barajas V, Photopoulou T, Langrock R, Patterson TA, Watanabe YY, Murgatroyd M, Papastamatiou YP. Analysis of animal accelerometer data using hidden Markov models. Meth Ecol Evol. 2017;8:161–73.View ArticleGoogle Scholar
- Vanderklift MA, Boschetti F, Roubertie C, Pillans RD, Haywood MDE, Babcock RC. Density of reef sharks estimated by applying an agent-based model to video surveys. Mar Ecol Prog Ser. 2014;508:201–9.View ArticleGoogle Scholar
- Bradley D, Conklin E, Papastamatiou YP, McCauley DJ, Pollock K, Pollock A, Kendall BE, Gaines SD, Caselle JE. Resetting predator baselines in coral reef ecosystems. Sci Rep. 2017;7:43131.View ArticlePubMedPubMed CentralGoogle Scholar
- Gleiss AC, Wilson RP, Shepard ELC. Making overall dynamic body acceleration work: on the theory of acceleration as a proxy for energy expenditure. Meth Ecol Evol. 2011;2:23–33.View ArticleGoogle Scholar
- Langrock R, King R, Matthiopoulos J, Thomas L, Morales JM. Flexible and practical modelling of animal telemetry data: hidden Markov models and extensions. Ecology. 2012;93:2336–42.View ArticlePubMedGoogle Scholar
- Akaike H. Information theory and an extension of the maximum likelihood principle. Proceedings of the 2nd International Symposium on Information Theory. Budapest: Akademiai Kiado; 1973. p. 267–81.Google Scholar
- Richards SA. Testing ecological theory using the information-theoretic approach: examples and cautionary results. Ecology. 2005;86(10):2805–14.Google Scholar
- Demšar U, Buchin K, van Loon EE, Shamoun-Baranes J. Stacked space-time densities: a geovisualisation approach to explore dynamics of space use over time. GeoInformatica. 2015;19:85–1185.View ArticleGoogle Scholar
- Horne JS, Garton EO, Krone SM, Lewis JS. Analyzing animal movements using Brownian bridges. Ecology. 2007;88:2354–63.View ArticlePubMedGoogle Scholar
- Payne NL, Snelling EP, Fitzpatrick R, Seymour J, Courtney R, Barnett A, Watanabe YY, Sims DW, Squire Jr L, Semmens JM. A new method for resolving uncertainty of energy requirements in large water breathers: the ‘mega-flume’ seagoing swim-tunnel respirometer. Meth Ecol Evol 2015, 6: 668–677.Google Scholar
- White CR, Phillips NF, Seymour RS. The scaling and temperature dependence of vertebrate metabolism. Biol Lett. 2006;2:125–7.View ArticlePubMedGoogle Scholar
- Dowd WW, Brill RW, Bushnell PG. Standard and routine metabolic rates of juvenile sandbar sharks (Carcharhinus plumbeus), including the effects of body mass and acute temperature change. Fish Bull. 2006;104:323–31.Google Scholar
- Iosilevskii G, Papastamatiou YP. Relations between morphology, buoyancy and energetics of requiem sharks. Roy Soc Op Sci. 2016;3:160406.View ArticleGoogle Scholar
- Papastamatiou YP, Bodey TW, Friedlander AM, Lowe CG, Bradley D, Weng K, Priestley V, Caselle JE. Spatial separation without territoriality in shark communities. Oikos. 2018, In Press.Google Scholar
- Mourier J, Maynard J, Parravicini V, Ballesta L, Clua E, Domeier ML. Extreme inverted trophic pyramid of reef sharks supported by spawning groupers. Curr Biol. 2016;26:2011–6.Google Scholar
- McCauley DJ, DeSalles PA, Young HS, Gardner JPA, Micheli F. Use of high-resolution acoustic cameras to study reef shark behavioral ecology. J Exp Mar Biol Ecol. 2016;482:128–33.View ArticleGoogle Scholar
- Speed CW, Meekan MG, Field IC, McMahon CR, Bradshaw CJA. Heat-seeking sharks: support for behavioural thermoregulation in reef sharks. Mar Ecol Prog Ser. 2012;463:231–44.View ArticleGoogle Scholar
- Rasher DB, Hoey AS, Hay ME. Cascading predator effects in a Fijian coral reef ecosystem. Sci Rep. 2017;7:15684.View ArticlePubMedPubMed CentralGoogle Scholar
- McCoy K, Heenan A, Asher J, Ayotte P, Gorospe K, Gray A, Lino K, Zamzow J, Williams I. Pacific Reef Assessment and Monitoring Program. 2015, NOAA Coral Reef Ecosystem Program.Google Scholar
- Sandvick H, et al. Modelled drift patterns of fish larvae link coastal morphology to seabird colony distribution. Nat Commun. 2016;7:11599.View ArticleGoogle Scholar
- Barraquand F, Inchausti P, Bretagnolle V. Cognitive abilities of a central place forager interact with spatial aggregation in their effect on intake rate. Anim Behav. 2009;78:505–14.View ArticleGoogle Scholar
- Fagan WF, Lewis MA, Auger-Methe M, Avgar T, Benhamou S, Breed G, LaDage L, Schlagel UE, Tang W, Papastamatiou YP, Forester J, Mueller T. Spatial memory and animal movement. Ecol Lett. 2013;16:1316–29.View ArticlePubMedGoogle Scholar
- Jacoby DMP, Papastamatiou YP, Freeman R. Inferring animal social networks and leadership: applications for passive monitoring arrays. J Roy Soc Inter. 2016;13:20160676.View ArticleGoogle Scholar
- Mourier J, Vercelloni J, Planes S. Evidence of social communities in a spatially structured network of a free-ranging shark species. Anim Behav. 2012;83:389–401.View ArticleGoogle Scholar
- Brewster LR, Dale JJ, Guttridge TL, Gruber SH, Hansell AC, Elliott M, Cowx IG, Whitney NM, Gleiss AC. Development and application of a machine learning algorithm for classification of elasmobranch behavior from accelerometry data. Mar Biol. 2018;165:62.View ArticlePubMedPubMed CentralGoogle Scholar
- DeRuiter SL, Langrock R, Skirbutas T, Goldbogen JA, Calambokidis J, Friedlander AS, Southall BL. A multivariate mixed hidden Markov model for blue whale behaviour and responses to sound exposure. Ann Appl Stat. 2017; In PressGoogle Scholar
- McClintock BT, Russell DJF, Matthiopoulos J, King R. Combining individual animal movement and ancillary biotelemetry data to investigate population-level activity budgets. Ecology. 2013;94:838–49.View ArticleGoogle Scholar