Fine‐scale foraging habitat selection by two diving central place foragers in the Northeast Atlantic

Abstract Habitat selection and spatial usage are important components of animal behavior influencing fitness and population dynamic. Understanding the animal–habitat relationship is crucial in ecology, particularly in developing strategies for wildlife management and conservation. As this relationship is governed by environmental features and intra‐ and interspecific interactions, habitat selection of a population may vary locally between its core and edges. This is particularly true for central place foragers such as gray and harbor seals, where, in the Northeast Atlantic, the availability of habitat and prey around colonies vary at local scale. Here, we study how foraging habitat selection may vary locally under the influence of physical habitat features. Using GPS/GSM tags deployed at different gray and harbor seals’ colonies, we investigated spatial patterns and foraging habitat selection by comparing trip characteristics and home‐range similarities and fitting GAMMs to seal foraging locations and environmental data. To highlight the importance of modeling habitat selection at local scale, we fitted individual models to colonies as well as a global model. The global model suffered from issues of homogenization, while colony models showed that foraging habitat selection differed markedly between regions for both species. Despite being capable of undertaking far‐ranging trips, both gray and harbor seals selected their foraging habitat depending on local availability, mainly based on distance from the last haul‐out and bathymetry. Distance from shore and tidal current also influenced habitat preferences. Results suggest that local conditions have a strong influence on population spatial ecology, highlighting the relevance of processes occurring at fine geographical scale consistent with management within regional units.


| INTRODUC TI ON
Understanding species' distribution and relationships with habitat is central in ecology and to the development of strategies for wildlife management and conservation (Morris, 2003;Rhodes et al., 2005). This is particularly true in marine ecosystems as oceans face increasing threats from overexploitation and habitat destruction (Halpern et al., 2008;Jones et al., 2018). In this ecosystem, distribution of a species is shaped by interactions between internal (i.e., species' physiological tolerance, dispersal and reproductive strategies) and external factors (i.e., environmental features and regional richness).
Depending on these different pressures, individuals disproportionately use the available conditions and resources, defining habitat selection (Mayor et al., 2009). This habitat selection differs from use or association as it implies choice and is commonly measured as use relative to availability or as use versus non-use (Mayor et al., 2009).
Habitat selection is an important component of animal behavior and a key determinant of individual survival, reproductive success, and ultimately population dynamic. Habitat selection can be defined at different scales (Johnson, 1980). Firstly, a coarse-scale selection pertains the selection of species geographical range, determined by the dispersal and ability of species to relocate (Morris, 1992). At smaller (i.e., local) scale, habitat selection determines the use of various habitat characteristics of an individual or a social group within a home range. In the case of central place foragers, such as seabirds and pinnipeds, habitat selection at local scale is mainly performed around the colony. Then, in a metapopulation composed of different colonies in distinct geographical areas, habitat selection will vary locally due to the variation in physical habitat features and community structure. Furthermore, intra-and interspecific interactions might also drive spatial usage and consequently habitat selection through resource exploitation (i.e., prey depletion, Vance, 1984). The majority of studies on habitat selection of central place foragers focus on only one or two colonies within the metapopulation but few have compared the habitat selection across multiple colonies through a global analysis to establish a metapopulation-level framework (e.g., Wakefield et al., 2017). However, understanding the causes of variation in habitat selection for local populations and determining how trends can be organized in space and time represents a major challenge in ecology (Fortin et al., 2008). Such an approach allows the development of coherent explanatory frameworks that carry across discrete populations and provide insights into phenotypic plasticity and eco-evolutionary dynamics.
Selection of resources can be considered as the expression of different behavior forms (i.e., dispersal and migration) of an animal in a particular environment (Schoener, 1969). At a local scale, foraging behavior is perceived as the major behavior which ultimately influences reproductive success and survival rates (Breed et al., 2009;Morris, 1992) and was already taken into account in different studies focusing on habitat selection (Donazar et al., 1993;Duchamp et al., 2004;Monsarrat et al., 2013). Several studies on the ecology of marine central place foragers, such as pinnipeds and seabirds, used telemetry devices to incorporate foraging behavior in their analyses (Guinet et al., 1997;Hamer et al., 2001;Jonsen et al., 2007;Shiel et al., 1999).
In the Northeast Atlantic, gray seals (Halichoerus grypus) and harbor seals (Phoca vitulina) are two sympatric species occurring along the European continental, Irish and British coasts, with differing population trends (SCOS, 2017). Gray seal's and harbor seal's core populations are located in the UK with an estimated 141,000 and 43,500 seals, respectively (SCOS, 2017). In France and Ireland, colonies of both species are located at their southern and western limit of their European range, respectively. Seals move regularly between colonies and can remain at sea for a long period. However, they display a high degree of site fidelity and mainly forage around their haulout sites (Cronin et al., 2013;Huon et al., 2015;McConnell et al., 1999;Sjöberg & Ball, 2000). Both species are generalist and benthic feeders. Their diets vary regionally and seasonally (Breed et al., 2006;Hammond et al., 1994;McConnell et al., 1999;Spitz et al., 2010). Colonies of both species are located across a range of habitats and diet composition is likely related to prey availability and abundance surrounding the haulout region (Gosch et al., 2019).
This context represents an excellent case for comparing intra-and interspecific foraging habitat selection and spatial usage at local scale and to establish an interesting framework for seal species' metapopulations. Many gray and harbor seals have been tracked from different colonies in the Northeast Atlantic, with analysis focusing on habitat selection from only one or two sites (Aarts et al., 2008;Bailey et al., 2014;Huon et al., 2015); or at the population scale  pooling multiple datasets. In this study, we combined for the first time several datasets for the assessment of foraging habitat selection across the wider population, but at local scale. Analyzing the data in this way enables new insights into gray and harbor seal's ecology across the Northeast Atlantic. We aimed at 1) studying the foraging habitat selection of gray and harbor seals at a colony (i.e., local) scale and 2) investigating the relationship between population trends and physical habitat features on the seals' spatial patterns and foraging habitat selection.

| Study areas
Gray seals were tracked in 5 regions (Figure 1 to consider the proportion of seals at sea during the surveys. These conversion factors were 0.72 and 0.24 (SCOS, 2020), respectively, for harbor seals and gray seals. When seals haul out continuously along the shore, a radius of 50 km was chosen around the capture site to count the raw number of seals. These calculations gave an estimate of local gray seal numbers of 488 seals for ICS, 996 in the Irish Sea , 2,008 in the FoT (Morris et al., 2021), 542 in the Iroise Sea, and 688 in the EEC (Vincent et al., 2017). Local harbor seals numbers were estimated at 476 in the Kenmare Bay , 4,411 in the Inner Hebrides, 69 in the FoT (Morris et al., 2021), 107 in the BdM, 152 in the BdV, and 629 in the BdS (Vincent et al., 2017) seals, respectively. Within all regions where both species are hauling out (Inner Hebrides, FoT, and EEC), we obtained tracking data from both species at FoT and EEC, providing the opportunity to study potential influence of interspecific interactions on foraging habitat selection and spatial usage.   (Table 1).

| Breeding and molting period
During the breeding and molting periods, seals tend to strongly reduce their foraging activities, increasing the amount of time spent hauled out, or staying close to their haulout sites (Boness, 1984;Caudron et al., 2009;Lidgard et al., 2003). We aimed to focus on foraging activities outside these periods when seals need to replenish their body reserves and are less constrained to haulout sites. Data obtained during the breeding (September to December for gray F I G U R E 1 Map of the different study areas for gray seals (in orange): Irish Continental Shelf (ICS), Irish Sea, Firth of Tay (FoT), Iroise Sea, and Eastern English Channel (EEC); and harbor seals (in blue): Kenmare Bay, Inner Hebrides, Firth of Tay (FoT), Baie du Mont Saint michel (BdM), Baie des Veys (BdV), and Baie de Somme (BdS); including the tagging sites (black stars) seals, June to July for harbor seals) and early molting (January to February for gray seals and August to September for harbor seals) periods were excluded from the analyses.

| Return-trip selection
In the Northeast Atlantic, harbor seals undertake short movements from their haulout sites (10-20 Km) and show long-term site fidelity (Ries et al., 1997;Tollit et al., 1998;Vincent et al., 2010). In contrast, telemetry on gray seals showed frequent movements between colonies (McConnell et al., 1999). They can alternate return trips to their haulout site in specific areas (within areas where most of their foraging activities occur) but they also frequently travel over hundreds of kilometers to distinct haulout site (SCOS, 2017). To study habitat selection and spatial usage of gray seals, we only selected individuals' foraging return trips around the tagging colony (i.e., defined as haulout sites grouping), excluding travel trips and return trips from other areas (McConnell et al., 1999).
Bathymetry was obtained from the European Marine Observation and Data Network (EMODnet), with a grid size resolution of 0.125 × 0.125 min. 2 Sediment data were obtained from the MESH_ EUNIS model (Mapping European Seabed Habitat project), which predicts habitat types with a spatial resolution of 300 m. Sediment types were based on a simplified FOLK classification system (Folk, 1954) and limited to the most dominant types: rock, mud, sand, gravel, coarse, and mixed sediments. We used different tidal current datasets for Irish, Scottish, and French areas, scaled to similar resolutions for inter-site comparison. Datasets for the French study areas were obtained from Previmer (Lecornu & De Roeck, 2009) for the tracking period. These were created from the MARS 2D model with a resolution of 250m and were available at an hourly scale. The Irish Marine Institute provided tidal current data for the Irish Continental Shelf and the Irish Sea 3 . Data were obtained from a numerical model with a spatial resolution varying between 1.2 Km and 1.5 Km and corresponded to surface tidal current at 3 hr interval. This dataset did not cover the Kenmare bay (for which no tidal current data were available). We averaged model datasets for the French and Irish areas, respectively, in order to represent the tidal strength in space irrespective of instant tidal phases (ebb, slack, or rising tide). Tidal current data for Scottish sites were obtained from the Web vision renewable website and were calculated from the ABP mer model (Atlas of UK marine Renewable Energy Resources 2008 4 ). These data corresponded to the peak current speed of a mean spring tide (m.s -1 ), with a spatial resolution decreasing from 200 m to 5 Km from inshore to offshore areas. TA B L E 1 Details of gray and harbor seals fitted with GPS/GSM tags in the different study areas: Iroise Sea, Eastern English Channel ( Jones et al. (2015) The distance between each GPS location and the last haulout and the distance to the shore were also included as explanatory variables to describe accessibility to the environment (Aarts et al., 2008

| Identification of seals foraging dives
We identified foraging behavior using a vertical approach based on two diving criteria: the dive shape and vertical descent speed (Planque et al., 2020). First, we excluded dives with a maximum depth <3 m and a dive duration <30 s, considering that these very shallow and short dives are unlikely to be foraging dives. The two dive criteria indicating benthic foraging were determined at the individual level to consider potential inter-individual variability in diving due to physiological conditions, swimming capacities, and individual strategies (Austin et al., 2006;Beck et al., 2003), as suggested by Planque et al. (2020). Dive-shape was determined using the Time Allocation at Depth index (Fedak et al., 2001). TAD values vary from 0 to 1, where 0 correspond to dives close to the surface (i.e., when the animal spent most of its dive time at a shallower depth than the maximum depth) and 1 represents "U-shape" dives (i.e., when the animal spent most of its dive time at the maximum depth). Harbor and gray seals are generally considered benthic feeders and therefore mostly perform U-shaped dives when they forage. Following Planque et al. (2020), we selected the most U-shaped dives for each individual by selecting 25% of the highest TAD values. Because Ushaped dives with a very low vertical descent speed are more likely associated with resting/sleeping behavior, we excluded a further 10% of the most U-shaped with the lowest vertical descent speeds for each individual. We therefore selected ~22.5% of all dives that are more likely associated with foraging, and we used the interpolated location of these dives in habitat selection models.

| Use-availability design
Following the use-availability design (Keating & Cherry, 2004), we assessed the foraging habitat selection by comparing the environmental characteristics of foraging dive locations to those of randomly generated points (i.e., pseudo-absences), representing the habitat availability (Aarts et al., 2008;Johnson et al., 2006;Keating & Cherry, 2004;Lele & Keim, 2006). Two pseudo-absences per foraging dive point were created locally for each individual seal within the different study areas using the package sp in R. These random points were created in each study area within buffers three times the size of local population Minimum Convex Polygon (Burgman & Fox, 2003) of all seals' dive locations, limited by the continental shelf (seals do not travel further offshore). For the EEC and the FoT, where both species were present, one buffer was created for each species.

| Modeling analyses
We fitted generalized additive mixed models (GAMM) to the data, with the gam function mgcv R package. We used a binomial family argument with a logit-link function to estimate the parameters of an inverse-logit selection model based on seal foraging dives and random points (Johnson et al., 2006). Foraging dives and pseudoabsences were the response variable, taking the values 1 and 0, respectively. To consider intra-individual autocorrelation, we included individual as a random effect. Environmental variables were treated as fixed effects. Bathymetry, tidal current, distance from shore, and distance from the last haulout were included as continuous variables; sediments were treated as categorical variable. When one sediment type was over-represented, the model was forced to consider this sediment type as reference level (otherwise reference sediment type was included alphabetically). The multicollinearity between covariates was assessed using the VIF value (Kutner et al., 2004). The best model was selected using the AIC (Akaike, 1973). Firstly, we fitted one model per site for each species to focus at the local scale.
Secondly, to highlight the importance of modeling habitat selection around colonies, particularly when local habitat characteristics differ, we also fitted a global model for each species using pooled data from all colonies and included "site" as a random factor. We did not include Kenmare for harbor seals' global model, as tidal current dataset was not available for this area. However, as we had problem of convergence when running global model, we only included one diving point on three. For models fitted for each colony, we calculated the importance of each covariate using the prediction function of the GAMM, providing an index of the relative importance of each covariate in the chosen model. Maps of habitat selection predicted by the model were created with ArcGIS for all sites.

| Influence of intra-and interspecific interactions on spatial patterns and home-range segregation
Trip characteristics and measures of similarity between home ranges were used to evaluate the influence of intra-and interspecific interactions on spatial patterns, to get complementary information on foraging habitat selection. For each species, trips with duration lower than 3 hr were removed as they were considered to be in the vicinity of haulout sites (Cronin et al., 2013). We used trip duration and maximum distance from the haulout sites (values were log-transformed to correct for non-normal distribution). Interpolated tracks were used for these trip characteristics. To reduce sampling bias (between areas where seals spend more time diving or out of the water), we interpolated all GPS locations every 20 min using straight-line inter-  Bhattacharyya, 1943) to quantify spatial overlap in home range. This method quantifies the spatial overlap between two population spatial distribution (Fieberg et al., 2005) and provides a value ranging from 0 (i.e., no overlap or complete segregation) to 1 (i.e., complete overlap). We applied the BA on the 95% Kernel density of foraging dive locations between individuals of the same colony to study the influence of colony size (i.e., indirectly the density dependence), and between species when both species were tracked around the same colony (i.e., FoT and EEC) to study the influence of interspecific interactions. We used the Kerneloverlaphr function of the package adehabitatHR in R (Calenge, 2006;Fieberg, 2014).

| Foraging habitat selection
Only results obtained for models fitted at local scale are presented in this section. Global models highlighted an influence of colony in the foraging habitat selection (p < .05). Although overall the results obtained by the two types of models were mostly similar, some eco-  (Figure 2). In the Inner Hebrides and the FoT, harbor seals selected their foraging habitat in inshore and in distance to the shore.

| Influence of intra-and interspecific interactions on spatial usage and home-range segregation
For each site and species, the hypotheses of normality and homoscedasticity of trip duration and maximum extent were rejected (p < .05) leading to the use of the nonparametric Kruskal-Wallis test and Dunn test as a post hoc test. For both metrics and both species, the model containing only latitude was selected (Table 4) and its influence was not significant in all cases. Seal numbers at the colony did not influence the spatial pattern of each species.

| Gray seals' trip characteristics
Means of trip durations were significantly different between sites (p < .05, Figure 3, Appendix S4  To consider the differences in habitat availability between the different study areas for each species and make inter-site comparisons, we fitted several models at the local scale. To date, studies encompassing multiple sites with contrasting local habitats have used global models (e.g., Wakefield et al., 2017). We also fitted a global model for each species as it can give more ecological inferences, as here, Negative influence the ones of models at local scale, some ecological inconsistencies were highlighted. For example, the positive influence of bathymetry for harbor seal colonies in the English Channel, which is in contradiction with their ecology in this area as they remain very coastal.

| D ISCUSS I ON
Furthermore, the uncertainty obtained around the predictions (i.e., confidence intervals) for the global model, in particular for gray seals, is too important for one part of environmental characteristics.
This indicates that the uncertainty in the predictions is too large and might have an impact on the relevance of the results. These differences in predictions suggest problem of homogenization. This was also found by Paton and Matthiopoulos (2016), highlighting the limitations of global models when studying populations that respond to factors at the local scale. That is why, we chose to create a model for each haulout group in order to consider local habitat availability and difference of intraspecific interactions, to avoid this problem of homogenization. We did not include the distribution of seals' prey resources in the model, as these data were not available for most of the study areas. However, as gray and harbor seals are generalist benthic feeders, we used environmental characteristics that best matched the ecology of their potential prey.
Across models, distance from the last haulout was consistently the most important factor influencing the foraging habitat selection of gray seals in all study areas (from 45% to 76% of the explained deviance in the model respectively for the FoT and the Iroise Sea).
These results are consistent with previous studies throughout the range of gray seals in the Northeast Atlantic. Gray seals in the North Sea spent 43% of their time within 10Km of the haulout sites (McConnell et al., 1999) and preferentially selected habitat closer to haulout sites with a gradual decrease of habitat selection beyond tens of kilometers (Aarts et al., 2008). In the Baltic Sea (Sjöberg & Ball, 2000), noted short distance trips (from 10 to 15 Km), gray seals spending 75% of their time within a radius smaller than 50 Km around their haulout sites. Foraging habitat selection by gray seals was negatively influenced by bathymetry, but to a lesser extent than distance from haulout sites (varying from 10% to 40% in the explained deviance). Depending on study sites, the depth selection decreased until a depth of 50 m (Irish and Iroise Sea) and 100 m (FoT). These results are also consistent with previous studies noting usual dive depth between 10 to 80 m (Aarts et al., 2008;Tollit et al., 1998). Gray seals are generally considered as benthic feeders (Beck et al., 2003;Lydersen et al., 1994) and the influence of bathymetry on foraging habitat selection will presumably vary locally depending on the seabed topography and sediment type. In addition to the topography and depth accessibility, seals will search for seabed types favorable to their preys. Furthermore, distance from the shore and tidal current must also influence the behavior of seals, as they use them to orientate, to move, and to forage (Zamon, 2001(Zamon, , 2003. However, these two variables did not contribute much to the explained deviance of the models, and their influence varied among sites. Distribution of seals' prey resources was not included in the model, as the data were not available for some of the study areas. Results in gray seal spatial usage were in accordance with their habitat selection at the colony scale. There was a global concordance between the foraging habitats selected by gray seals in each site and their spatial pattern, as the prey found in their diet would also be found in the same habitat characteristics (Alheit & Hagen, 1997;Gosch et al., 2014;Hammond et al., 1994;Ridoux et al., 2007). Gray seals in the Firth of Tay and in the Eastern English Channel made longer trips but of shorter duration compared to gray seals in the other sites. In the North Sea, gray seals tend to travel long distances directly to offshore areas on specific sandbanks where sand eel availability is high (Hammond et al., 1994;McConnell et al., 1999;Wilson & Hammond, 2019). This offshore behavior is in accordance with our results for the FoT. Gray seals in the EEC performed their trips mostly in specific areas along the coast, as the EEC is known to be a major ground for flatfish (Carpentier et al., 2009;Riou et al., 2001;Selleslagh et al., 2009), which are observed in gray seal diet (Planque et al., 2021). In the ICS and Iroise Sea, gray seals made shorter trips.
These two sites are known as highly biologically productive regions in the Eastern North Atlantic (upwelling and area enclosed by specific currents, respectively, for the ICS and the Iroise Sea (Hily & Glémarec, 1999;Raine & McMahon, 1998). In the Iroise Sea, gray seals selected their foraging habitat around their haulout sites, in the kelp forest. In this area, they are known to forage on wrasse (Ridoux et al., 2007) found on rocky sediment around their haulout sites.
However, sediment selected by seals with foraging habitat modeling did not match with rocky substrate. This could be explained by the fact that a selection was made on dives deeper than 3 m, thus excluding those in very shallow waters around their haulout sites (Huon et al., 2015;Vincent et al., 2016). Gray seals in the ICS are located in and soles (Soleidae) (Gosch et al., 2019), which can be found over muddy and sandy area in shallow waters (Ellis et al., 2000), corresponding to habitat selected by gray seals in this area.
Distance from the last haulout, distance from shore, and/or bathymetry explained most of the deviance (>90%) in harbor seal foraging habitat selection. These three variables had a negative influence, but at different degrees depending on the site configuration. Harbor seals were very coastal and sedentary in the six study areas, which was also supported by their spatial patterns. This is consistent with previous findings on the species; in the Moray Firth (East of Scotland) for instance, seals forage within 30 Km of their haulout sites and dive at a maximum depth of 50 m (Bailey et al., 2014;Tollit et al., 1998). This was also highlighted on the other side of the Atlantic, in the Saint Laurent estuary (Lesage et al., 2004), seals were coastal (with distances shorter than 11 Km from the shore) and in shallow waters (<50 m deep). Tidal current and sediment types accounted less in the explained deviance than the other variables. Their influence was generally very low, with the exception of the tidal current for the FoT. Harbor seals feed on diverse fish species and their diet vary locally (Hall et al., 1998 (Spitz et al., 2015). Seal diet in BdM was not available so we could not compare seals' foraging habitat selection and known prey habitat in this region. However, on the East coast of Scotland, harbor seals selected habitat over sandy areas in front of the Tay river mouth, corresponding to the habitat of sandeelthe main harbor seal prey in the area (Wilson & Hammond, 2019).
Harbor seals in Inner Hebrides essentially forage on demersal fish (e.g., whitings Merlangius merlangus) and pelagic fish (e.g., herrings Clupea harengus) (Wilson & Hammond, 2019); therefore, their benthic component may be less important than in other sites and that could be concordant with their habitat selection along the sounds of the fjords. Sediments and current features selected by harbor seals might correspond to the habitat features of their prey.
Distance from the last haulout was the most important factor influencing habitat selection for both species and all colonies. This is in accordance with the theory of central place foraging where animals will minimize distance-dependent travel costs (Orians & Pearson, 1979) and shows the importance of including not only potential quality of habitat, but also elements related to species biology. This has been demonstrated in many central place foraging species, both herbivores (Gerwing et al., 2013;Shrader et al., 2012) and carnivores (Rainho & Palmeirim, 2011). Results of habitat selection and spatial usage for both species highlighted the importance of considering local-scale effects within a metapopulation. Indeed, despite the fact that the colonies are a few hundred kilometers apart, the colonies were located in contrasting environments with different habitat characteristics and associated prey availability. The predatorprey relationship varies spatially in association to the underlying physical conditions (Santora et al., 2014). Gray seal and harbor seal are two species with a high behavioral plasticity and the importance of behavioral plasticity in environmental adaptation has been ably demonstrated for multiple species, as for example in squirrels (Hefty & Stewart, 2019). All animals experience dynamic physiological and environmental demands and must adjust their activity patterns to minimize associated cost and increases fitness. Due to differences in habitat availability and habitat selection among colonies, our results on spatial usage did not highlight any potential effect of density dependence on foraging ranges. Various factors other than density dependence can also impact spatial patterns. Animals should adjust their foraging movements for habitat exploitation, depending on environmental heterogeneity and hierarchical distribution of resources (Pinaud & Weimerskirch, 2007). This was already demonstrated in seabird species, where the influence of the shape of the coast and the availability of habitat further away from the central place forces birds to forage further offshore without there necessarily being a density effect (Wakefield et al., 2017). In our study, this may have been the case of gray seals in the FoT foraging on sandbanks in the middle of the North Sea. Phenotypic plasticity was also highlighted to influence spatial usage for fur seals (Baylis et al., 2018).
The number of harbor seals has declined at some colonies in the North Sea during the last decade (Thompson et al., 2019). The decline of sandeel numbers (main prey of harbor seals around these colonies) and interspecific competition with gray seals were suggested as one of the potential causes in local harbor seal declines (Wilson & Hammond, 2019). In our study, gray seals and harbor seals were tracked in two areas where both species haul-out (FoT and EEC), and spatial partitioning between seal species was highlighted in habitat selection and spatial patterns in these areas. Gray seals made longer trips than harbor seals (15 Km for median maximum extents in gray seals in both areas, versus 5 Km and 4 Km respectively for harbor seals in the EEC and FoT). In both cases, harbor seals tended to forage in inshore areas, while at least some of the gray seals foraged further offshore, which is consistent with previous studies that also found differences in the use of marine environment between these two species Sharples et al., 2012). However, Planque et al. (2021) found trophic overlap between both species in the EEC.
In the North Atlantic, gray seals and harbor seals are managed at local scales, and in the absence of genetic information on population structuring, haulout groups are often considered as "Seal Management Units" (Russell et al., 2013). Both species are considered as generalist, using a variety of habitats and prey. This study highlights the importance of studying foraging habitat selection at local scale and considering the variability between colonies, as physical habitat features and seals' prey resources vary between regions. As marine top predators, both seal species are listed in the Annex II of the European habitat directive requiring establishment of protected areas to maintain favorable conservation status. At a local scale, our predictive maps of foraging habitat selection could be used by managers to implement specific areas of conservation to maintain a good ecological state of their habitat and prey resources potentially at risk due to anthropic activities. Such a foraging habitat selection analysis could be applied and/or adapted for other central place foraging species, in both marine and terrestrial ecosystems.

ACK N OWLED G M ENTS
We would like to thank the Sea Mammal Research Unit for providing telemetry data and raw numbers for gray seals and harbor Geographic Global Exploration Fund.

CO N FLI C T O F I NTE R E S T
The authors have no conflict of interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
The data that support the findings of this study are openly avail-