# Seascapes of fear and competition shape regional seabird movement ecology

Nicolas Courbin, Lorien Pichegru, Mduduzi Seakamela, Azwianewi Makhado, Michael Meÿer, Pieter Kotze, Steven Mc Cue, Clara Péron, David Grémillet

► To cite this version:

Nicolas Courbin, Lorien Pichegru, Mduduzi Seakamela, Azwianewi Makhado, Michael Meÿer, et al. "Seascapes of fear and competition shape regional seabird movement ecology." Communications Biology, 2022, 5 (1), pp.208. 10.1038/s42003-022-03151-z. hal-03650515

HAL Id: hal-03650515

https://hal.science/hal-03650515

Submitted on 3 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Seascapes of fear and competition shape regional seabird movement ecology

Nicolas Courbin ID 1,2,8✉, Lorien Pichegru³, Mduduzi Seakamela⁴, Azwianewi Makhado⁴,⁵, Michael Meÿer⁴, Pieter G. H. Kotze⁴, Steven A. Mc Cue⁴, Clara Péron⁶ & David Grémillet ID 5,7,8✉

Fear effects of predators on prey distributions are seldom considered in marine environments, especially over large spatial scales and in conservation contexts. To fill these major gaps, we tested the Seascape of Fear Hypothesis in the Benguela marine ecosystem off South Africa. Using electronic tracking data, we showed that Cape gannets and their predator, the Cape fur seal, co-occurred in daytime and competed with fisheries within coastal areas. At night, gannets are particularly vulnerable to seals, and 28% of the birds returned to the safety of their breeding colony. The remaining 72% slept at the sea surface, but shifted to offshore areas with lower seal attendance, reducing predation risk by 25%. Overall, our integrative study demonstrates how fear and competition shape the seascape of threatened Cape gannets within a marine environment perturbed by climate change and overfishing. Such knowledge has strong implications for the design of marine protected areas.

<footnote n="1">¹Centre d'Ecologie Fonctionnelle et Evolutive (CEFE), UMR 5175 CNRS - Université de Montpellier - Université Paul Valéry Montpellier 3 - EPHE - IRD, Montpellier, France. ²Laboratoire d'Écologie Alpine (LECA), UMR 5553 Université Grenoble Alpes - Université Savoie Mont-Blanc - CNRS, Le Bourget-du-Lac, France. ³Institute for Coastal and Marine Research, Nelson Mandela University, Gqeberha, South Africa. ⁴Department of Forestry, Fisheries and the Environment, Branch Ocean and Coasts, Cape Town, South Africa. ⁵Percy FitzPatrick Institute of African Ornithology, University of Cape Town, Rondebosch, South Africa. ⁶Laboratoire de Biologie des Organismes et Ecosystèmes Aquatiques (UMR BOREA) MNHN, CNRS, IRD, SU, UCN, UA, Paris, France. ⁷Centre d'Etudes Biologiques de Chizé (CEBC), UMR 7372 CNRS - La Rochelle Université, Villiers-en-Bois, France. ⁸These authors contributed equally: Nicolas Courbin, David Grémillet. ✉email: ncourbin@gmail.com; david.gremillet@cebc.cnrs.fr</footnote>
Seascape ecology is on the rise. Decades after the advent of landscape ecology, its aquatic counterpart has taken a giant footstep¹. Indeed, this fusion of geography and ecology benefited from geospatial revolution, with the rapid development of global data acquisition and mapping tools. Those allow assessing the effect of structural environmental conditions upon the biogeography of marine organisms. In this context, a major constraint is the fear effect, i.e., the spatial variation in prey perception of predation risk, termed the Landscape of Fear (LoF)², and later extended to the Seascape of Fear (SoF)³. Fear effects modify prey behavior, induce diel migratory patterns⁴ and entail energetic costs, with cascading effects onto population dynamics and entire ecosystems⁵,⁶. The fearscape concept has gained substantial interest across the last decade, with renewed focus on temporal dynamics⁴,⁵. Surprisingly, such fear effects are not addressed by Pittman¹ in their seminal work on seascape ecology.

We tested the SoF hypothesis at a regional scale in the ecologically perturbed Benguela upwelling ecosystem off South Africa. There, Cape gannets (*Morus capensis*) are at risk from predation by Cape fur seals (*Arctocephalus pusillus pusillus*), specifically when resting at the water surface during the night when they are most vulnerable to underwater predator attacks. Such predation events have been documented on land and at sea⁷,⁸ (Supplementary Data 1), but their impacts on Cape gannet at-sea behavior has never been tested, thereby blurring our understanding of non-consumptive effects at a regional scale. Moreover, gannets compete with seals and regional fisheries for the consumption of diminishing stocks of small pelagic fish (sardines *Sardinops sagax* and anchovies *Engraulis encrasicolus*), but are provisioned with fishery wastes dumped at sea by trawlers targeting hake (*Merluccius spp.*)⁹. In this context, we hypothesized that spatiotemporal variations in fear and competition contribute to shaping the spatial ecology of gannets in the Benguela upwelling ecosystem. Specifically, we predicted that (1) birds would avoid seals at all times, especially during the riskiest nighttime period. This was predicted to induce diel gannet spatial shifts at sea, towards safer areas at night. However, during the less

risky daytime period, gannets have to balance predation risk by seals with their foraging needs. We predicted that (2) birds would compete with seals and fisheries while feeding at sea during daytime. Alternatively, we hypothesized that thermoregulatory constraints also shaped gannet at sea distributions. We took advantage of a large multi-year data set, which included electronic tracking of at-sea movements for seals and gannets in the Southern Benguela upwelling ecosystem, as well as the distribution of fishing activities and sea surface temperatures in this region. On the basis of this information, we assessed the at-sea spatial ecology of the different marine predators involved in the Benguela seascape.

## Results

Testing prediction 1: Seals at sea were active night and day, with similar habitat selection patterns and a decrease in occurrence with increasing distance from the coast (Fig. 1, Supplementary Note 1, Supplementary Fig. 1, Supplementary Table 1). Gannets displayed two tactics when at sea: 1) Some (28%) came back to the colony at night after daytime foraging. 2) Most (72%) remained at sea in-between daytime foraging bouts, and spent 76% of nighttime resting at the sea surface (Fig. 1, Table 1, Supplementary Note 2). In agreement with our first prediction, all gannets, regardless of whether they remained at sea at night, or came back to the colony, selected foraging areas that decreased their risk of encountering seals during daytime (Fig. 2, Supplementary Note 3, Supplementary Table 2). Crucially, birds which rested at the sea surface overnight displayed an offshore spatial shift (13.1 km) between locations dedicated to daytime foraging and nighttime resting (Fig. 3, Table 1, Supplementary Note 4, Supplementary Table 3). This shift allowed resting gannets to decrease their risk of encountering seals by 25% at night, on average, whereas nighttime resting waters were on average only 0.3 °C warmer than daytime coastal foraging areas (Table 1, Supplementary Note 4). Indeed, the likelihood of nighttime resting over daytime foraging continuously increased with the decrease in seal occurrence, except in areas with high encounter

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Seascape of fear induced by Cape fur seals for Cape gannets. **a** Daytime seascape of fear (yellow to red index of seal occurrence) and daytime foraging locations for gannets returning to the colony at night (gray dots) or resting at sea all night (green dots). **b** Nighttime seascape of fear (yellow to red index of seal occurrence) and gannet nighttime at-sea resting locations (blue dots). Seals (n = 25 individuals) were Argos-tracked from three colonies (black triangles) and gannets (n = 197 individuals) were GPS-tracked from Malgas Island (black dot).</figcaption></figure>
<table><caption>Table 1 Summary statistics of habitat use for 197 Cape gannets.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="2">Daytime foraging, spent night at colony<br>N = 6890 locations, 53 birds</th><th colspan="2">Daytime foraging, spent night at sea<br>N = 47953 locations, 144 birds</th><th colspan="2">Nighttime resting, spent night at sea<br>N = 85142 locations, 142 birds</th></tr><tr><th>Mean (95% CI)</th><th>Min;max</th><th>Mean (95% CI)</th><th>Min;max</th><th>Mean (95% CI)</th><th>Min;max</th></tr></thead><tbody><tr><td>Distance to shore (km)</td><td>41.9 (37.5;46.5)</td><td>0.0;95.6</td><td>42.1 (37.3;47.2)</td><td>0.2;157.4</td><td>57.0 (48.5;63.2)</td><td>0.5;142.0</td></tr><tr><td>Distance to the colony (km)</td><td>63.7 (50.3;78.0)</td><td>0.7;154.8</td><td>88.6 (71.9;105.7)</td><td>0.3;278.1</td><td>96.6 (81.8;110.5)</td><td>18.3;270.6</td></tr><tr><td>Bathymetry (m)</td><td>-238.0 (-284.0; -190.4)</td><td>-1480.9; -6.0</td><td>-242.5 (-265.5; -220.8)</td><td>-1465.6; -2.6</td><td>-358.7 (-397.9; -320.3)</td><td>-2253.9; -8.4</td></tr><tr><td>Seal occurrence index</td><td>6.6 (5.9;7.4)</td><td>1;10</td><td>6.4 (6.0;6.9)</td><td>1;10</td><td>4.8 (4.3;5.4)</td><td>1;10</td></tr><tr><td>Log(purse-seiner catch (tonnes)+1)</td><td>2.9 (1.4;4.4)</td><td>0.0;9.9</td><td>3.2 (2.3;4.2)</td><td>0.0;10.2</td><td>1.8 (1.0;2.6)</td><td>0.0;9.4</td></tr><tr><td>Log(trawler catch (tonnes)+1)</td><td>2.9 (1.9;3.8)</td><td>0.0;6.6</td><td>2.6 (2.4;2.9)</td><td>0.0;6.6</td><td>3.6 (3.2;4.0)</td><td>0.0;6.8</td></tr><tr><td>Sea surface temperature (°C)</td><td>16.6 (16.4;16.8)</td><td>13;19</td><td>16.6 (16.4;16.8)</td><td>13;19</td><td>16.9 (16.6;17.1)</td><td>13;19</td></tr></tbody></table>

We considered gannet habitat use during the daytime foraging and nighttime resting periods and in relation to their nighttime resting behavior (rest at colony or rest at sea) during October-November between 2008 and 2015, estimated with mixed linear models accounting for repeated measures within individual and year.

<figure><img src="image_4.png" /><figcaption>Fig. 2 Daytime foraging habitat selection of Cape gannets. Habitat selection coefficients of foraging behavior during daytime, depending on nighttime resting tactics for 197 Cape gannets from Malgas Island during October-November between 2008 and 2015. Error bars represent the 95% confidence interval. Confidence interval including 0 means no selection.</figcaption></figure>

risk, but gannets remained more likely to forage than rest in such areas (Fig. 4a, Supplementary Note 4, Supplementary Tables 4, 5). In contrast, gannets behaved equally with an increase in sea surface temperature, until high-temperature values where resting became more likely than foraging (Fig. 4b). Therefore, the large decrease in predation risk outperformed a small thermoregulatory benefit, on average. Nighttime periods spent at the sea surface may also allow adult gannets to digest a first food load, before catching another the next morning, to be brought back to their chick¹⁰.

Testing prediction 2: During the day, gannets ventured into areas with higher risk of seal encounter compared to nighttime (Fig. 4a, Table 1), but still selected foraging areas that reduced seal co-occurrence (Fig. 2, Supplementary Note 3). Gannets competed with fisheries during the day, and were more likely to forage in areas with an increase in average catch of small pelagic fish by purse-seiners, relative to nighttime (Fig. 4c, Supplementary Note 4, Supplementary Table 5). Conversely, gannets were less likely to forage during the day, than rest at night, in areas with high average catch of hake by trawlers (Fig. 4d, Supplementary Table 5). On average, daytime gannet foraging areas had higher average catch of small pelagic fish by purse-seiners and lower average catch of hake by trawlers, compared to gannet nighttime resting areas (Table 1). However, at a daytime scale, foraging habitat selection of gannets did not rely on fisheries (Fig. 2, Supplementary Note 3).

Discussion

Using a multi-year, large electronic tracking data set of gannet and seal at-sea movements and fisheries' log-books, we validated our hypothesis: both fear and competition contribute to shaping regional gannet spatial ecology in the Benguela upwelling ecosystem. Gannets foraged in risky and competitive areas during the day. They also used movement tactics aiming to reduce spatial overlap with seals, especially during the riskiest night period. This involved coming back to the colony or performing a marked
<figure><img src="image_5.png" /><figcaption>**Fig. 3** *Diel spatial shift of 144 Cape gannets which stayed at sea at night during October-November between 2008 and 2015.* **a** Distribution of gannet nighttime resting locations (dark gray) and daytime foraging locations (light gray) in relation to distance to shore. **b** At-sea spatial distribution of daytime foraging (green dot) and nighttime resting locations (blue dot) of gannets. Triangles show the three monitored seals colonies. Map generated with the marmap⁴¹ package for R.</figcaption></figure>

offshore spatial shift to safer waters. Our study therefore adds substantial support to the idea that seascape ecology must account for the fear factor³,⁶,¹¹.

Interestingly, studies of zooplankton vertical migration first pointed to the LoF, but this concept has been mainly studied within terrestrial landscapes⁴,⁵. Indeed, so far only 4.3% (10/230, ISI Web of Science April 2021) of publications which addressed the LoF²,⁵ did so in a marine context. Most such studies focused on coastal reef habitats¹²,¹³, to the neglect of pelagic ecosystems³,¹⁴ and large spatial scales¹¹,¹⁵. We fill these major gaps, by providing, to the best of our knowledge, the first evidence for a SoF in an upwelling ecosystem, at a regional scale. Our analyses confirm that seabirds avoid seals¹⁶ whenever they can, but in the Benguela upwelling ecosystem other top predators add layers of complexity to the pelagic fearscape. Notably, great white sharks (*Carcharodon carcharias*) are known to target both seals and seabirds¹⁷, even though they are far less abundant than seals. Orcas (*Orcinus orca*) may also take seals and great white sharks. We are therefore calling for SoF studies considering trophic networks within pelagic ecosystems, to rigorously understand fear effects on the spatiotemporal distribution of aquatic animals⁶, and their population consequences⁵.

The Benguela ecosystem is being challenged by the combined effects of climate change and overfishing, threatening a series of endemic seabird species: Cape gannet, Cape and Bank Cormorants (*Phalacrocorax capensis* and *P. neglectus*), and African penguin (*Spheniscus demersus*). In this context, marine protected areas (MPAs) are an effective way to safeguard the integrity of ecosystems, and of threatened seabirds¹⁸. However, MPAs have so far been designed and implemented on the basis of knowledge of seabird critical foraging habitats, to the neglect of the SoF imposed onto seabirds by their potential predators. Considering our fearscape study in the Benguela upwelling ecosystem, we posit that seascape ecology and forthcoming marine management should consider the SoF when designing MPAs.

**Methods**

**Cape gannet movement tracking.** The study took place in the Western Cape, South Africa, where we studied chick-rearing Cape gannets from Malgas Island

(33.05° S, 17.93° E) during October–November from 2008 to 2015 (Fig. 1). We caught birds using a pole fitted with a loop and fitted 197 adult Cape gannets (22 in 2008, 16 in 2009, 38 in 2010, 11 in 2011, 29 in 2012, 29 in 2013, 23 in 2014, and 29 in 2015) with GPS-loggers (2008: GPS mass 65 g, i.e., 2.4 % adult body mass, Technosmart, Rom. 2009–2010: GPS mass 45 g, i.e., 1.7% adult body mass, Technosmart, Rom. From 2011: GPS mass 30 g, 1.1% of bird body mass, Catnip Technologies, Hong-Kong). Loggers were attached to the lower back with waterproof Tesa® tape and recorded position at a regular 30-s to 2-min intervals, reinterpolated over 1-min intervals. Devices were recovered after one foraging trip lasting a few hours to one week. Bird handling and tracking using these procedures do not have a measurable impact on foraging behavior¹⁹,²⁰. We caught adult birds at-random from the colony, and previous studies showed that this resulted in a well-balanced sex-ratio preventing confounding sex effects²¹. All experiments were performed under permit from South African National Parks with respect to animal ethics (N° RYAP/AGR/001-2002/V1).

**Cape gannet movement tactics and behavioral phases.** We identified two movement trip tactics for Cape gannets: After their daytime foraging activities, some birds returned to the colony at night (*rest at colony*) tactic) while others spent all the night at sea (*rest at sea*) tactic). Within the GPS tracks of gannets from these two categories, we discriminated resting, foraging, and commuting phases, with a segmentation-clustering method based on smoothed speed (i.e., speed smoothed over two steps before and after the focal location) and turning angle measured at constant step length. This corresponded to the angle between the focal location, the first location entering a circle of radius equal to the median step length, and the last location inside the circle²². We fitted behavioral identification with the *segclust2d* package²³ for the R software²⁴. See complete details on behavioral classification for Cape gannets tracks in Appendix 1 in Courbin et al.²⁵.

**Cape fur seal movement tracking and the seascape of fear.** We assessed the at-sea spatial distribution of Cape fur seals, a predator of Cape gannet fledglings⁷ and adults (Supplementary Data 1). We used Argos data collected from 25 lactating female seals before (2003 and 2004) and again concomitantly with gannet tracking (2012 and 2014). Seals were tracked during the same period of the year as gannets (i.e., September to November). Adult females nursing pups were selected at random and captured using a modified hoop net. Once restrained, anesthesia was induced using isoflurane gas delivered via a portable vaporizer (Stinger, Advanced Anesthesia Specialists, Gladesville, New South Wales, Australia). A satellite tag was glued to the guard hairs on the upper back. Individuals were allowed to recover from the anesthesia and resumed normal behavior within 45 min of capture. Throughout the process, the animals' breathing was closely monitored and their flippers were repeatedly flushed with seawater to prevent hyperthermia. Seals were equipped with Argos satellite transmitters at three colonies (Fig. 1): Kleinsee (29°35'09"S, 16°59'56"E) located ~400 km to the North of the gannet colony (n = 8 seals in 2003 and 2004); Vondeling Island (33°09'11"S, 17°58'57"E), ~12 km away from the gannet colony (n = 12 seals in 2012 and 2014); and Geyser Rock (34°41'19"S, 19°24'49"E) located ~230 km to the South of the gannet colony
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Probability of daytime foraging over nighttime resting for gannets in relation to environmental features (n = 142 birds). We predicted the effect of a seal occurrence, b sea surface temperature, c purse-seiner catches, and d trawler catches, with a mixed effect logistic regression and keeping all other covariates at their mean. Grey areas represent the 95% confidence interval. The x-axis is back-transformed (unscaled). Distribution and mean (dotted line) are shown for daytime foraging (green) and nighttime resting (blue) locations.</figcaption></figure>

(n = 5 seals in 2003). Seals at Vondeling Island were equipped with Argos-linked Spot-6 position transmitting tags (Wildlife Computers) following deployment procedures outlined in Kirkman et al.²⁶. Seals at Kleinsee and Geyser Rock were equipped with ST18 and ST20 satellite-linked platform terminal transmitters (Telonics, Mesa, USA), as detailed in Skern-Mauritzen et al.²⁷. Devices collected a well-balanced number of Argos locations during the day (n = 6080 locations) and at night (n = 6501 locations). See full details on seal tracking in Supplementary Table 6. All fieldwork was permitted by the Animal Ethics Committee of the Department of Environmental Affairs and Tourism's Marine and Coastal Management branch, which at the time was the management authority of South Africa's marine and coastal environment (Ref: DEAT2006-06-23).

We modeled both daytime and nighttime at-sea occurrences of seals for each colony with resource selection functions (RSF)²⁸,²⁹, a proxy of the fear effect for Cape gannets. RSF compared environmental features of seal's at-sea Argos positions (i.e., further 500 m than the colony) with five times more random

locations that captured the breadth of environmental conditions available to seals. We sampled random locations for each individual within the yearly area used by seals from each colony, delineated by the 95% kernel utilization distribution of the Argos locations of all seals of the colony. RSF were fitted with a generalized linear mixed model with a binomial distribution for errors. As environmental variables, we considered bathymetry (m), the slope of the bathymetry (°) and the distance to the colony (km) within the RSF. These variables were not highly correlated (|r| ≤ 0.61) and had low collinearity with a variance inflation factor VIF < 2³⁰.

All continuous predictors were centered and scaled. Following statistical recommendations in Muff et al.²⁹, we added random intercept for seal ID with a large fixed variance, a random slope for each predictor and weighted random locations by 1000. We assessed the robustness of the RSF using a leave-one-out cross validation with iteratively one of the individuals representing the testing set and the other seals representing the training set³¹,³². RSF with a high predictive power had a high average Spearman's rank correlation ($\bar{r}_s$) between the rank of the
RSF scores (relative probabilities of seal occurrence) split into ten bins and the area-adjusted frequency of the Argos locations³¹. We ran RSF with the *glmmTMB* package³³ for the R software²⁴.

Finally, we spatially predicted the binned RSF scores (10 bins) for each seal colony within our study area³⁴. We assessed the seascape of fear by overlapping the three colony-specific RSF maps and retaining the maximum value among the three maps for each pixel.

**Co-occurrence of Cape gannets and fisheries.** To assess the influence of fisheries on gannet movements we used vessel log-book records between 2008 and 2015. We mapped the yearly catch distribution of trawlers targeting hake (*Merluccius capensis*) with a 20 × 20 nautical mile resolution grid and of purse-seiners targeting anchovy (*Engraulis capensis*) and sardine (*Sardina pilchardus*) with a 10 × 10 nautical mile resolution grid. These data were made available by the branch: Fisheries Management of the Department of Forestry, Fisheries and the Environment of the Republic of South Africa.

**Daytime gannet foraging habitat selection.** We assessed whether Cape gannets experienced different risks of encountering seals and different levels of competition with fisheries, according to their movement tactics. For this analysis, we only kept daytime foraging locations, and estimated foraging habitat selection of gannets with a RSF for each tactic²⁸,²⁹. Thereby, we fitted a generalized linear mixed model with a binomial distribution for errors with random intercept for gannet ID (with a large fixed variance) and year, random slope for each predictor and weighted available locations by 1000. We determined availability using a design adapted for central place foragers in habitat selection studies, considering that individuals used areas close to the colony more frequently than elsewhere³⁵,³⁶. For each observed foraging trip, we simulated 10 tracks that started at the same location as the observed trip (i.e., at the colony) using a first-order vector autoregressive model³⁷. Simulated tracks considered no habitat preference, while respecting constraints on trip structure (duration and travel speed)³⁷. For each simulated track, we then determined available foraging locations using the segmentation-clustering method described above (see “Cape gannet movement tactics and behavioral phases” section). RSF included daytime seal encounter risk, as well as the presence of purse-seiners and trawlers, assessed through the spatial distribution of their yearly catches (log(tonnes + 1)). All continuous predictors were centered and scaled, and had low collinearity (|r| ≤ 0.62, VIF < 2)³⁰. We assessed RSF robustness as previously described, but using a k-fold cross validation with iteratively 80% of the individuals representing the training set and the 20% remaining birds representing the testing set³¹,³². We ran RSF with the *glmmTMB* package³³ for the R software²⁴.

**Cape gannet diel at-sea habitat use.** We tested whether Cape gannets adjusted their at-sea behavior according to the presence of seals and fishing boats between relatively safe daytime and risky nighttime. For this purpose, we only considered gannets that did not return to the colony at night (n = 142 individuals, 72% of our sample). We estimated the probability of daytime foraging and nighttime resting, depending on both fishing activities and seal encounter risk. Thereby, we calculated daytime seal encounter risk for daytime foraging locations and nighttime seal encounter risk for nighttime resting locations. We used a generalized linear mixed model with a binomial distribution for errors and individual ID and year as random intercepts. Models also included the average monthly sea surface temperature (SST, °C) to test whether by moving away from risky areas to rest, gannets may also benefit from better thermoregulation conditions, i.e., higher SST. SST data were extracted from Aqua MODIS satellite imagery, with a 4-km-resolution grid (NASA Goddard Space Flight Center, Ocean Ecology Laboratory, Ocean Biology Processing Group; (2018): Aqua MODIS Sensor Ocean Color Data, NASA OB.DAAC.). All continuous predictors were centered and scaled. We tested candidate models with linear or nonlinear effects (natural spline with df = 4) for predictors, and selected the best model using the Akaike's Information Criterion corrected for finite sample size. All candidate models did not include highly correlated variables (|r| < 0.4) and had low collinearity with VIF < 2³⁰. We fitted models with the *lme4* package³⁸ and model selection with the *MuMln* package³⁹ for the R software²⁴.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

Data are openly available in figshare at https://doi.org/10.6084/m9.fishshare.172990940⁴⁰, except Argos data for Cape fur seals that are under embargo until the end of 2023 as they are being used in the PhD thesis of Mduzudi Seakamela at the Department of Forestry, Fisheries and the Environment, Branch Ocean and Coasts, Cape Town, South Africa.

Received: 15 September 2021; Accepted: 8 February 2022;

Published online: 04 March 2022

## References

1. Pittman, S. J. (ed.). *Seascape Ecology* (John Wiley & Sons, 2017).

2. Laundré, J. W., Hernández, L. & Altendorf, K. B. Wolves, elk, and bison: reestablishing the 'landscape of fear' in Yellowstone National Park, USA. *Can. J. Zool.* **79**, 1401–1409 (2001).

3. Wiring, A. J., Heithaus, M. R., Frid, A. & Dill, L. M. Seascapes of fear: evaluating sublethal predator effects experienced and generated by marine mammals. *Mar. Mammal. Sci.* **24**, 1–15 (2008).

4. Courbin, N. et al. Zebra diel migrations reduce encounter risk with lions at night. *J. Anim. Ecol.* **88**, 92–101 (2018).

5. Gaynor, K. M., Brown, J. S., Middleton, A. D., Power, M. E. & Brashares, J. S. Landscapes of fear: spatial patterns of risk perception and response. *Trends Ecol. Evolution* **34**, 355–368 (2019).

6. Ainley, D. G. & Ballard, G. Non-consumptive factors affecting foraging patterns in Antarctic penguins: a review and synthesis. *Polar Biol.* **35**, 1–13 (2012).

7. Makhado, A. B., Crawford, R. J. & Underhill, L. G. Impact of predation by Cape fur seals *Arctocephalus pusillus pusillus* on Cape gannets *Morus capensis* at Malgas Island, Western Cape, South Africa. *Afr. J. Mar. Sci.* **28**, 681–687 (2006).

8. Crawford, R. J. M. et al. Trends in numbers of Cape gannets (*Morus capensis*) 1956/1957 – 2005/2006, with a consideration of the influence of food and other factors. *ICES J. Mar. Sci.* **64**, 169–177 (2007).

9. Cohen, L. A. et al. Changes in prey availability impact the foraging behaviour and fitness of Cape gannets over a decade. *Mar. Ecol. Prog. Ser.* **505**, 281–293 (2014).

10. Ropert-Coudert, Y. et al. A fine-scale time budget of Cape gannets provides insights into the foraging strategies of coastal seabirds. *Anim. Behav.* **67**, 985–992 (2004).

11. Matthews, C. J. D., Breed, G. A., LeBlanc, B. & Ferguson, S. H. Killer whale presence drives bowhead whale selection for sea ice in Arctic seascapes of fear. *Proc. Natl Acad. Sci. USA* **117**, 6590–6598 (2020).

12. Heithaus, M. R. et al. State-dependent risk-taking by green sea turtles mediates top-down effects of tiger shark intimidation in a marine ecosystem. *J. Anim. Ecol.* **76**, 837–844 (2007).

13. Catano, L. B. et al. Reefsccapes of fear: predation risk and reef heterogeneity interact to shape herbivore foraging behaviour. *J. Anim. Ecol.* **85**, 146–156 (2016).

14. Bishop, A. M., Brown, C. L., Sattler, R. & Horning, M. An integrative method for characterizing marine habitat features associated with predation: a case study on juvenile steller sea lions (*Eumetopias jubatus*). *Front. Mar. Sci.* **7**, 576716 (2020).

15. Hammerschlag, N. et al. Evaluating the landscape of fear between apex predatory sharks and mobile sea turtles across a large dynamic seascape. *Ecology* **96**, 2117–2126 (2015).

16. Masello, J. F., Kato, A., Sommerfeld, J., Mattern, T. & Quillfeldt, P. How animals distribute themselves in space: Variable energy landscapes. *Front. Zool.* **14**, 1–14 (2017).

17. Wcisel, M., O'Riain, M. J., de Vos, A. & Chivell, W. The role of refugia in reducing predation risk for Cape fur seals by white sharks. *Behav. Ecol. Sociobiol.* **69**, 127–138 (2015).

18. Pichegru, L., Grémillet, D., Crawford, R. J. M. & Ryan, P. G. Marine no-take zone rapidly benefits endangered penguin. *Biol. Lett.* **6**, 498–501 (2010).

19. Grémillet, D. et al. Offshore diplomacy, or how seabirds mitigate intra-specific competition. *Mar. Ecol. Prog. Ser.* **268**, 265–279 (2004).

20. Pichegru, L. et al. Foraging behaviour and energetics of Cape gannets *Morus capensis* feeding on live prey and fishery discards in the Benguela upwelling system. *Mar. Ecol. Prog. Ser.* **350**, 127–136 (2007).

21. Lewis, S. et al. Sex-specific foraging behaviour in a monomorphic seabird. *Proc. R. Soc. Lond.: B Biol. Sci.* **269**, 1687–1693 (2002).

22. Patin, R., Etienne, M.-P., Lebarbier, E., Chamaillé-Jammes, S. & Benhamou, S. Identifying stationary phases in multivariate time series for highlighting behavioural modes and home range settlements. *J. Anim. Ecol.* **89**, 44–56 (2020).

23. Patin, R., Etienne, M.-P., Lebarbier, E. & Benhamou, S. Segclust2d: Bivariate segmentation/clustering methods and tools. R package version 0.1.0. https://cran.r-project.org/web/packages/segclust2d (2018).

24. R Development Core Team. *R: A language and environment for statistical computing*. Version 3.6.2. http://www.R-project.org/ (R Foundation for Statistical Computing, 2019).

25. Courbin, N. et al. The dance of the Cape gannet may contain social information on foraging behaviour. *Anim. Behav.* **166**, 95–108 (2020).

26. Kirkman, S. P. et al. Dive behaviour and foraging effort of female Cape fur seals *Arctocephalus pusillus pusillus*. *R. Soc. Open Sci.* **6**, 191369 (2019).

27. Skern-Mauritzen, M. et al. Do inter-colony differences in Cape fur seal foraging behaviour reflect large-scale changes in the northern Benguela ecosystem? *Afr. J. Mar. Sci.* **31**, 399–408 (2009).
28. Manly, B., McDonald, L., Thomas, D. L., McDonald, T. L. & Erickson, W. P. *Resource Selection by Animals: Statistical Design and Analysis for Field Studies* (Springer Science and Business Media, 2002).

29. Muff, S., Signer, J. & Fieberg, J. Accounting for individual-specific variation in habitat-selection studies: efficient estimation of mixed-effects models using Bayesian or frequentist computation. *J. Anim. Ecol.* **89**, 80–92 (2020).

30. Dormann, C. F. et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. *Ecography* **36**, 27–46 (2013).

31. Boyce, M. S., Vernier, P. R., Nielsen, S. E. & Schmiegelow, F. K. A. Evaluating resource selection functions. *Ecol. Model.* **157**, 281–300 (2002).

32. Roberts, D. R. et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. *Ecography* **40**, 913–929 (2017).

33. Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zeroinflated generalized linear mixed modeling. *R. J.* **9**, 378–400 (2017).

34. Morris, L. R., Proffitt, K. M. & Blackburn, J. K. Mapping resource selection functions in wildlife studies: Concerns and recommendations. *Appl. Geogr.* **76**, 173–183 (2016).

35. Monsarrat, S. et al. How predictability of feeding patches affects home range and foraging habitat selection in avian social scavengers? *PLoS ONE* **8**, e53077 (2013).

36. Courbin, N. et al. Short-term prey field lability constrains individual specialisation in resource selection and foraging site fidelity in a marine predator. *Ecol. Lett.* **21**, 1043–1054 (2018).

37. Raymond, B. et al. Important marine habitat shift off east Antarctica revealed by two decades of multi-species predator tracking. *Ecography* **38**, 121–129 (2015).

38. Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48 (2015).

39. Barton, K. MuMln: Multi-Model Inference. R package version 1.43.17 (2020).

40. Courbin, N. et al. dataset/Seascapes_of_fear_and_competitionShape regional_seabird_movement_ecology. Figshare. https://doi.org/10.6084/ m9.fishshare.17299094 (2021).

41. Pante, E. & Simon-Bouhet, B. marmap: A package for importing, plotting and analyzing bathymetric and topographic data in R. *PLoS ONE* **8**, e73051 (2013).

## Acknowledgements

This work was funded by the Centre National de la Recherche Scientifique (CNRS), the University of Cape Town and Department of Forestry, Fisheries and the Environment of the Republic of South Africa (DFFE). We thank David G. Ainley for inspiration, as well as Pierre Nel and South African National Parks for accommodation and logistical support. We are also grateful to all persons who helped with the fieldwork, and all staff at the FitzPatrick Institute of African Ornithology and DFFE. Vessel log-book records were made available by the branch: Fisheries Management of the DFFE. We acknowledge

Anne Loison for her support. This is a contribution to the Excellence Chair ECOMM funded by the Region Nouvelle Aquitaine.

## Author contributions

N.C. and D.G. designed the research from an original idea of D.G. N.C., L.P., M.S., A.M., M.M., P.K., S.C., C.P., and D.G. gathered data. N.C. performed the analyses. N.C. and D.G. wrote the manuscript, with inputs from L.P., M.S., A.M., M.M., P.K., S.C., and C.P.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-022-03151-z.

**Correspondence** and requests for materials should be addressed to Nicolas Courbin or David Grémillet.

**Peer review information** *Communications Biology* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Quan-Xing Liu and Caitlin Karniski.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022