Introduction
Perception of predation risk induces physiological stress (Brown & Kotler, Reference Brown and Kotler2004; Clinchy et al., Reference Clinchy, Zanette, Boonstra, Winfield and Smith2004), stress-related behavioural responses such as avoidance/escape and vigilance (Childress & Lung, Reference Childress and Lung2003; Armitage, Reference Armitage2017), and changes in patterns of aggregation (Winnie & Creel, Reference Winnie and Creel2007) or movement (Sih & McCarthy, Reference Sih and Mccarthy2002). These responses are expected to conflict with other important fitness-related activities such as foraging, reinforcing energetic constraints on reproduction or survival (Sinclair & Arcese, Reference Sinclair and Arcese1995; Brown & Kotler, Reference Brown and Kotler2004). Because animals often perceive human presence as a form of predation risk, they react by displaying typical antipredator behavioural responses (Frid & Dill, Reference Frid and Dill2002). For example, tourism-related activities are known to influence vigilance patterns (Beauchamp, Reference Beauchamp2015). In marine animals that spend part of their biological cycle inshore as they feed, or reproduce near the shoreline, the potential for these disturbances increases with increased tourism. Following disturbance caused by marine tourism, animals shift from fundamental behavioural activities to vigilance or escape behaviour. For example, in bottlenose dolphins Tursiops spp. off the coast of Scotland, boat presence reduced foraging activity by 49% (Pirotta et al., Reference Pirotta, Merchant, Thompson, Barton and Lusseau2015). In Shark Bay, Australia, experimental boat approaches to dolphin Lagenorhynchus obscurus pods elicited more erratic travel direction compared to travel behaviour before and after approaches (Bejder et al., Reference Bejder, Samuels, Whitehead and Gales2006). Humpback whales Megaptera novaeangliae off the coast of Ecuador increased their swimming speed in the presence of boats (Scheidat et al., Reference Scheidat, Castro, González and Williams2004). Harp seal Phoca groenlandica pups spent less time resting, and the vigilance time of their mothers increased, in response to tourist presence, at the expense of suckling time for the pups (Kovacs & Innes, Reference Kovacs and Innes1990). A meta-analysis comparing cetacean responses to vessel approach across 15 studies (Senigaglia et al., Reference Senigaglia, Christiansen, Bejder, Gendron, Lundquist and Noren2016) found consistent evidence of disruptions of activity budgets and path directionalities. The viability of some cetacean populations has been affected, with behavioural disruptions leading to a decrease in female reproductive success and a consequent decrease in population size (Lusseau et al., Reference Lusseau, Slooten and Currey2006; Currey et al., Reference Currey, Dawson and Slooten2009).
Here, we focus on the whale shark Rhincodon typus, categorized as Endangered on the IUCN Red List (Pierce & Norman, Reference Pierce and Norman2016), in Bahia de Los Angeles, Gulf of Baja California, in northern Mexico. We experimentally mimicked the effect of ecotourism, allowing us to standardize the sources of disturbances (i.e. the approaches of vessels and swimmers) and investigate four aspects of whale shark–tourist interactions: (1) vigilance behaviour, (2) behaviours that reflect stress, (3) foraging probability, and (4) habituation to disturbance.
We investigated vigilance behaviour for two reasons. Firstly, monitoring vigilance is a more specific proxy of risk perception than, for example, a change in direction of a foraging individual, which may only reflect a change of foraging route following food distribution; the observation of the same individuals both with and without human stimulus allowed us to assess to what extent human approach caused changes in vigilance investment. Secondly, vigilance may be a reliable proxy of potential foraging costs (Blanchard & Fritz, Reference Blanchard and Fritz2007). We hypothesized that the probability of displaying behaviours expected to reflect stress (vigilance, changes in direction, diving or acceleration) should increase after an experimentally simulated tourist approach. To investigate whether tourist trips affect sharks foraging probability, we hypothesized that if tourist trips affect the perception of risk, the probability that whale sharks would be foraging should decrease following tourist approaches. As whale sharks are individually recognizable, we investigated whether they could habituate to disturbance: for individuals observed several times over the course of one season, we evaluated whether disturbance-related behaviours decreased in subsequent approaches compared to the first approach.
Study area
The study was carried out in Bahia de los Angeles, Gulf of Baja California, México (28°59′–29°03′N, 113°30′–113°26′W), a regular location for whale shark aggregation during June–November. This bay drains into the Gulf of California and is an important foraging habitat, with patches of high copepod densities (Nelson & Eckert, Reference Nelson and Eckert2007; Hernández-Nava & Alvarez-Borrego, Reference Hernández-Nava and Alvarez-Borrego2013). At the study site whale sharks can be observed approaching 2–10 m from the shore (authors, pers. obs.). Our study was conducted in September–November 2014, when whale shark sightings peak. This large planktivorous fish requires high density prey patches to make feeding energetically viable (Ronher et al., Reference Ronher, Armstrong, Pierce, Prebble, Cagua and Cochran2015). In the tropics, they face the challenge of finding sufficient food in an oligotrophic environment in which prey are sparsely distributed and their quantity and location are dynamic in space and time (Ronher et al., Reference Ronher, Armstrong, Pierce, Prebble, Cagua and Cochran2015). The occurrence of seasonal aggregations of whale sharks in the bay, following prey patches, has stimulated the development of ecotourism focused on these animals, with a potential cost for the sharks when disturbed during foraging. The Official Mexican Standard (Norma Official Mexicana: NMX-AA-142-SCFI-2008) provides guidelines for ecotourism activities, focusing on whale sharks while they are in Mexican waters. Touching or interrupting the activity of whale sharks is prohibited. However, in Bahia de los Angeles the presence of park rangers allowed to penalize these behaviours is restricted to a few days per month, and only outside the whale shark sighting season.
Methods
Experimental design
Whale sharks were located during the day in the area that tourism operators commonly use for trips. We used a 7 m panga (a skiff-like vessel with an outboard motor), driven by experienced tour operators, to approach whale sharks that were spotted at the water surface, typically while foraging. Vessels approached from behind and to one side of the whale shark, travelling at no more than 8 miles/h (c. 13 km/h), as established in the tourist operators regulations (Rodríguez-Dowdell et al., Reference Rodríguez-Dowdell, Enríquez-Andrade and Cárdenas-Torres2007). We only recorded encounters when whale sharks were at the surface and were first spotted from 10–20 m. We discarded encounters in which sharks appeared close to the vessel. A total of three different swimmers performed the trials, but only one swimmer participated per trial. The swimmer used snorkeling equipment for every interaction and was experienced in monitoring whale sharks. In each trial the swimmer always dived after the vessel approached, and swam from behind or approached on one side of the animal and at a low speed (slower than the vessel speed).
Whale sharks can be individually identified from the unique patterns of white dots and stripes on their dorsal surface. Photographs of whale sharks have been compiled since 2007 by the tourist management group PEJESAPO. We used this database to identify the individuals monitored. We used scan-sampling observations (Altmann, Reference Altmann1974) to record the behaviours of whale sharks in response to experimental approaches.
Whale shark behaviour was recorded by one observer, from the vessel, during four periods: immediately when the whale shark was sighted (i.e. at 10–20 m from the vessel) and before the vessel approached, after the vessel approached (until 4–5 m from the shark), after the swimmer approached (until 0.5–5 m from the shark), and (4) after human stimuli (i.e. after the swimmer returned to the vessel and the vessel was > 10 m from the whale shark). The observer recorded the following behaviours: foraging, vigilance, change of direction, diving and acceleration. Foraging individuals referred to animals gulping large amounts of water while swimming slowly or floating just below the water surface, following Ketchum et al. (Reference Ketchum, Galván-Magaña and Klimley2013) and Nelson & Eckert (Reference Nelson and Eckert2007). Vigilance was recorded when a whale shark was floating immobile at a shallow depth, with its mouth closed, sometimes orienting its head towards the vessel or swimmer (a behaviour referred to as alert by Cubero-Pardo et al., Reference Cubero-Pardo, Herron and González-Perez2011). Diving was recorded when individuals moved deeper in the water column. Acceleration was recorded when animals increased their horizontal swimming speed. We recorded date and time of sighting, location (using a global positioning system), sex (obtained by the swimmer diving under the whale shark, looking for presence or absence of claspers), and size of whale sharks observed. Size was determined by comparing the whale shark's total body length with the length of the vessel and swimmer (Nelson & Eckert, Reference Nelson and Eckert2007; Ketchum et al., Reference Ketchum, Galván-Magaña and Klimley2013). In September 2014 hurricane Odile hit the study area, significantly reducing tourism during the next 2 months, enabling us to perform the experimental trials without the confounding effect of approaches by tourists.
Statistical analyses
As the behaviour of each shark was recorded during four periods (see above) and as more than one individual could be observed on the same day, we included individual identity nested in observation date as random terms in all analyses, to avoid pseudoreplication. Mixed models were run using the lme4 package (Bates et al., Reference Bates, Mächler, Bolker and Walker2015) in R 3.2.3. (R Development Core Team, 2013). All the response variables were binary. We used general linear mixed models (GLMMs) with a binomial distribution and logit link function, and with the same explanatory variables: the period of observation (before human stimuli, after vessel approach, after swimmer approach, after overall human stimuli), sex, whale shark length, and whether the individual had or had not previously been observed in the area. Because we were interested in the effect of tourists on whale shark behaviour and because of the limited sample size, we only considered two-way interactions involving the period of observation. As we did not find any habituation effect we included all the observations recorded on the same individuals in the first three analyses (see below). In the first two analyses we focused on the individuals for which we had data for the four periods, whereas for the third analysis we considered the first and last period only, therefore increasing sample size. This is explained below. We discarded data for eight individuals of unknown sex. A total of 77 individuals were observed at least once.
Model selection began with a full model and terms were successively removed, beginning with interaction effects. We tested the change in deviance after removal of a term using a χ2 test. Whenever an interaction was tested, the main effects comprising the interaction were retained in the model. As probability values obtained from likelihood-ratio tests on GLMMs are not robust (Bolker et al., Reference Bolker, Brooks, Clark, Geange, Poulsen, Stevens and Whie2009), we retested all significant effects using a bootstrap method (with 200 replications): these probability values are indicated as Pboot. Model predictions and standard errors were obtained using the AICcmodavg package in R.
Are whale sharks more vigilant when approached?
The binary response variable was the state of vigilance of the individual: vigilant/not vigilant. However, because sharks were never vigilant before human stimuli and only one was vigilant after the vessel departed, models did not necessarily converge. Hence we first performed a χ2 test to investigate the relationship between being vigilant or not and the four observation periods, focusing only on the 44 individuals with data for all four periods (176 observations). We ran the GLMM on the two periods with sufficient variability in whale shark vigilance behaviour (after vessel and after swimmer; 130 observations).
Are whale sharks more stressed when approached?
The binary response variable was the display of behaviours expected to indicate stress (vigilance, change of direction, dive or acceleration): stress display/no stress display. Because of the absence of variability in shark stress before and after human stimuli (almost no sharks were stressed during these periods), we used the approach described above, and sample sizes are the same.
Do people affect whale shark foraging?
The binary response variable was the foraging behaviour of the individual: foraging/not foraging. We compared the observation period before and after overall human stimuli (i.e. two observations per individual per session; N = 68 individuals, 200 observations).
Do whale sharks exhibit within-season behavioural habituation?
We re-ran the above three analyses focusing on individuals observed at least twice in the season (twice: N = 22 individuals; three times: N = 10 individuals). We discarded individuals observed four and five times, as this applied to only two individuals. Given the small sample size, we only considered two-way interactions involving observation order. Because of convergence issues, we standardized all explanatory variables (except observation order, as it had three levels) using the arm package in R (Gelman, Reference Gelman2008).
Results
Are whale sharks more vigilant when approached?
No individuals (of 44) were vigilant before human stimuli and only one after, whereas vigilance was more common during approach by the vessel and by a swimmer (χ23 = 40.97, P < 0.0001; Fig. 1a). Whale sharks were 56% more likely to be vigilant after a swimmer's approach than after that of a vessel (estimate β 0.70 ± SE 0.38, Pboot = 0.044). They were also 106% more likely to be vigilant when they were new individuals at the location than when they had been seen previously (β 1.09 ± SE 0.45, Pboot = 0.010; Fig. 1b). Neither sex nor length had a significant effect (all P > 0.11) and none of the interactions with observation period (after vessel or after swimmer) significantly affected shark vigilance (all P > 0.84).
Probability of being vigilant (GLMM prediction ± SE) in relation to (a) observation period (after vessel or swimmer approach, together with data before and after human stimuli, for comparison), and (b) an individual's status in the population (previously observed or newly observed at the location).
Are whale sharks more stressed when approached?
As with vigilance behaviour, no individuals (of 44) were stressed before human stimuli and only four after, whereas stress was more common after encounters with a vessel or a swimmer (χ23 = 68.72, P < 0.0001; Fig. 2). Observation period was the only significant predictor (all other P > 0.15): sharks were 67% more likely to show stress after approach by a swimmer than after approach by a vessel (β 1.37 ± SE 0.46, Pboot = 0.010; Fig. 2). None of the interactions with observation period significantly affected shark stress (all P > 0.14).
Probability of exhibiting stress-related behaviour (see text for details; GLMM prediction ± SE) in relation to the observation period (after vessel or swimmer approach, together with data before and after human stimuli, for comparison).
Do people affect whale shark foraging behaviour?
Whale sharks were 24% more likely to forage before human stimuli than after (before: β 0.78 ± SE 0.08; after: β 0.63 ± 0.08; Pboot = 0.005 and there was a non-significant tendency for longer individuals to be observed foraging (β 0.64 ± SE 0.25; Pboot = 0.08). None of the interactions with observation period significantly affected the probability of observing whale sharks foraging (all P > 0.10).
Do whale sharks exhibit within-season behavioural habituation?
There was no evidence of behavioural habituation, regardless of the behaviour considered: none of the interactions between observation order and other variables were significant (all P > 0.37).
Discussion
Risk effects arise when prey alter their behaviour in response to predators, and this response carries costs (Creel & Christianson, Reference Creel and Christianson2008). In Bahia de los Angeles, mimicking ecotourism, in the form of vessels and swimmers approaching whale sharks, induced stress-related behavioural responses in both males and females. Both vessels and swimmers increased the probability of displaying vigilance, and stress-related behaviours more generally. Swimmers appeared to be more stressful to whale sharks than vessels. However, this could be a cumulative effect, as in our experimental design the swimmer stimulus always followed the vessel stimulus. It would be of value to separate each stimulus, to understand better the effects of each type of approach, but it is unlikely that in actual tourism these two stimuli would be presented in isolation or in a different order and therefore the sequence of presentation used here is the most relevant for an understanding of whale shark–tourist interactions in practice. Although vigilance may allow an animal to collect more information about a potential threat, changing direction, diving or accelerating suggests a motivation to move away from human stimulus. Moving away from boats was registered in dolphins Tursiops truncatus (Lusseau, Reference Lusseau2004): if dolphin–boat interactions were longer than a mean of 68 minutes, they modified their activity to escape from the area. A meta-analysis showed that variation in path sinuosity captured cetacean response to vessel presence, suggesting that this response may reflect an attempt to avoid vessels without leaving the area (Senigaglia et al., Reference Senigaglia, Christiansen, Bejder, Gendron, Lundquist and Noren2016).
We found that after a simulated tourist trip the probability that a whale shark would forage was decreased, with the potential for long-term negative effects on individuals as a result of reduced food intake. However, given that sharks were first spotted when at the surface, typically during foraging, we cannot exclude the possibility that some sharks naturally stopped foraging irrespective of human presence. For the dusky dolphin Lagenorhynchus obscurus, boat presence negatively affected long-term foraging efficiency, as they stopped feeding, to travel (Dans et al., Reference Dans, Crespo, Pedraza, Degrati and Garaffo2008). If copepod abundance is subject to dynamic and seasonal changes in Bahia de los Angeles (Lavaniegos, et al., Reference Lavaniegos, Heckel and Ladron de Guevara Porras2012), a likely cost of disturbance to whale sharks, particularly during foraging bouts, is lowered energy intake. This potentially translates into a reduction in reproductive output or in long-term survival (Lima, Reference Lima1998), with local consequences on population trends.
The analysis of intra-seasonal habituation suggested no such effect in whale sharks. However, 2–3 observations may be insufficient to detect such a pattern. Moreover, what we considered the first encounter could have been the first for some individuals but not for others, as whale sharks could have had encounters in other parts of the Gulf. Our results on the effect of previous sightings in the area on vigilance do not necessarily mean there is habituation (although this may be the case after < 3 experiences with tourists; i.e. at the end of the season), but rather could suggest that more sensitive individuals do not return to this area of high tourist pressure. Further studies are needed to discriminate between these possibilities. In Ningaloo Reef, Australia, there was no evidence that high exposure to tourism led to inter-annual decreases of whale sharks in the bay, which was interpreted as individuals having become accustomed to encounters with tourists (Sanzogni et al., Reference Sanzogni, Meekan and Meeuwig2015). Individuals exhibiting site fidelity could be vulnerable to human disturbance, and selection may favour those individuals with personality traits allowing them to habituate to encounters with tourists.
We found a non-significant trend for larger individuals to feed more, regardless of the period of observation (before or after the approach procedure). A meta-analysis showed that small cetacean species were less likely to rest in the presence of vessels than larger species (Senigaglia et al., Reference Senigaglia, Christiansen, Bejder, Gendron, Lundquist and Noren2016). In other vertebrates, the best predictor of the direction and magnitude of tolerance to human approaches is body size: for example, larger birds are more tolerant than smaller ones (Samia et al., Reference Samia, Blumstein, Stankowich and Cooper2016).
In conclusion, it remains to be determined whether short-term behavioural responses to perception of predation risk negatively affect fitness. Declines in whale shark sightings have been reported in several locations. In western Australia abundance declined by 40% from 1995 to 2004 (Mau & Wilson, Reference Mau, Wilson, Irvine and Keesing2005; Bradshaw et al., Reference Bradshaw, Fitzpatrick, Steinberg, Brook and Meekan2008), and anecdotal reports from Belize suggest that as tourism has increased, whale sharks sightings have decreased (Graham, Reference Graham2007). The effects of exploitation via non-lethal effects of ecotourism on the largest fish need to be measured both in terms of the economic benefits for local people and effects on the behaviour, physiology and, ultimately, fitness of the animals. Our results suggest that preventing swimmers from approaching, if whale sharks stop feeding when a vessel approaches, could avoid further increasing stress. This suggestion is also in the Mexican guidelines regulating ecotourism activities focusing on whale sharks (Norma Official Mexicana: NMX-AA-142-SCFI-2008). As whale sharks can also be observed close to the shore, tourist operators could instead develop walking tours to observe whale sharks with fixed telescopes, thereby reducing disturbance. Future studies could also investigate the way distance to vessel and swimmer shapes risk perception by whale sharks, which would facilitate guidelines for operators regarding minimum distances.
Acknowledgements
This work was funded by the Mexican Commission on Natural Protected Areas. We thank the Pejesapo monitoring group and PRONATURA Noroeste for logistical support. MOB is grateful to ECOLAB of University Paul Sabatier for its support during her sabbatical stay. The laboratory EDB is part of the Laboratoire d'Excellence TULIP (ANR -10-LABX-41). Lynna Kiere reviewed the English. José Arce-Smith, Ricardo Arce-Navarro and Joel Prieto-Villavicencio drove the vessel during trials. Ian Velazco-Espinoza and Ricardo Arce-Navarro participated in fieldwork. We are grateful to the Facultad de Ciencias Biológicas, Universidad Autónoma del Estado de Morelos, and to Elisabet Wehncke and Adalberto Aguilar for their critiques.
Author contributions
Conception and design of the experiment: MOB, AMQ; recording of the data: AMQ, AVH; data analysis: PB, TM; discussion of results: all authors; writing: MOB; revisions: PB and TM.
Conflicts of interest
None.
Ethical standards
We complied with the Oryx Code of Conduct for authors.
