Estimating the relative use of anthropogenic resources by Herring Gull (Larus argentatus) in the Bay of Fundy, Canada

Gulls (Larus spp.) are described as generalist, opportunistic feeders that show great flexibility in habitat use. Despite an apparent advantage in changing landscapes, many Larus populations have declined in eastern North America since the 1990s. The main hypothesis explaining gull declines at a broad scale is a decrease in total food availability, especially anthropogenically derived fisheries discards and human refuse as industries and cities have improved their management practices. However, it is difficult to quantify the total proportion of gull diet subsidized by humans to test this hypothesis because many common prey items can be traced to both anthropogenic and nonanthropogenic sources. Our aim was to estimate the proportion of diet derived from anthropogenic food sources for Herring Gull (L. argentatus) during the breeding season at the two largest colonies in the Bay of Fundy, Canada, which are located 36 km apart. GPS loggers were deployed to quantify spatiotemporal movement patterns, and whole blood and feather samples were collected for δC and δN stable-isotope analysis to estimate diet composition during the incubation and chick-rearing stages of the breeding season. Results indicate that there is spatial segregation in the foraging areas used by gulls from the two colonies. All gulls relied on a variety of anthropogenic food sources, with some individuals selecting heavily on fisheries (i.e., active town wharfs, fish packaging plants, aquaculture pens) and mink (Neovison vison) fur farms. Landfills were not a significant source of food during the breeding season. Our study provides valuable information about the relative reliance of gulls on anthropogenic food subsidies, providing insight into how changing industry practices may affect patterns in nesting and foraging by gulls in the region. Estimation de l'utilisation relative de ressources anthropiques par le Goéland argenté (Larus argentatus) dans la baie de Fundy, Canada RÉSUMÉ. Les goélands (Larus spp.) sont décrits comme des mangeurs généralistes et opportunistes qui font preuve d'une grande souplesse dans l'utilisation de leur habitat. Malgré un certain avantage procuré par la modification des paysages, de nombreuses populations de Larus ont diminué dans l'est de l'Amérique du Nord depuis les années 1990. La principale hypothèse derrière cette diminution à grande échelle serait la baisse de la disponibilité générale de nourriture, en particulier de rejets de pêche d'origine anthropique et de déchets d'origine humaine, les industries et les villes ayant amélioré leurs pratiques de gestion. Cependant, il est difficile de quantifier la proportion totale du régime alimentaire des goélands fournie par l'homme pour tester cette hypothèse, car de nombreuses proies communes peuvent provenir à la fois de sources anthropiques ou non anthropiques. Notre objectif était de calculer la proportion du régime alimentaire provenant de sources anthropiques chez le Goéland argenté (L. argentatus) pendant la saison de reproduction dans les deux plus grandes colonies de la baie de Fundy, au Canada, qui sont situées à 36 km l'une de l'autre. Nous avons déployé des enregistreurs GPS pour quantifier l'ensemble des déplacements spatiotemporels, et pris des échantillons de sang entier et de plumes pour faire l'analyse des isotopes stables δ13C et δ15N afin d'estimer la composition du régime alimentaire pendant les étapes d'incubation et d'élevage des poussins. Nos résultats indiquent qu'il y a une ségrégation spatiale dans les zones d'alimentation utilisées par les goélands des deux colonies. Tous les goélands dépendaient d'une variété de ressources alimentaires anthropiques, certains individus se nourrissant principalement des produits d'activités de pêcheries (c.-à-d. aux quais municipaux actifs, aux usines de transformation des poissons, aux bassins d'aquaculture) et aux élevages de visons (Neovison vison). Les dépotoirs n'étaient pas une source importante de nourriture pendant la saison de reproduction. Notre étude fait état d'informations précieuses sur la dépendance relative des goélands aux ressources alimentaires d'origine anthropique, et révèle comment l'évolution des pratiques industrielles peut affecter les habitudes de nidification et de recherche de nourriture des goélands dans la région.


INTRODUCTION
Generally, migratory birds are of conservation interest because of the essential roles they play in ecosystems, their cultural appeal and symbolic importance, and because population declines have been observed on local to global scales (Cury et al. 2011, Rosenberg et al. 2019. Avian species have additionally served as indicators of ecosystem health in conservation and restoration projects worldwide (Weseloh et al. 2002, Frederick et al. 2009, Scopel et al. 2018, Velarde et al. 2019. Despite these values, gulls (Larus spp.) are something of an outlier, often labeled as a "nuisance" to human health and industry (Blokpoel and Tessier 1986, Cleary et al. 2006, Clark et al. 2013, Province of New Brunswick 2013 as well as a threat to other species of conservation (e.g., terns, salmon) and social or economic (e.g., eiders) interest (Donehower and Bird 2008, LaRue 2016, Scopel and Diamond 2017. Nuisance concerns about gulls have caused the focus of wildlife management to shift away from mitigation of human practices that are incompatible with the maintenance of ecological integrity to the control of gulls. Wildlife managers may consider mitigating both the causes and symptoms of apparent nuisance situations created by wildlife, but feasibility and perceived simplicity often favor strategies directed at the symptoms. This approach may lead to a commitment to the inefficient allocation of limited conservation resources because the causes of environmental imbalance are not addressed.
Long-term monitoring indicates that some of the most common North American Larus species have declined asynchronously throughout North America since the 1990s or earlier (Nisbet et al. 2013, Mittelhauser et al. 2016) despite their reportedly flexible foraging habits and habitat use Annett 1991, Nisbet et al. 2017). The family Laridae, which includes both gulls and terns, has declined in abundance in North America by 50.5% (90% confidence interval: 39.9-58.4%) between 1970 and 2017, with 16 of 22 of Laridae species having negative trends, including Herring Gull (L. argentatus) and Great Black-backed Gull (L. marinus; Rosenberg et al. 2019). Although rooftop nesting, which may be missed by traditional coastal seabird surveys, might account for some of this reported decline, past surveys that include urban areas generally report the proportion of rooftop nests to be < 15% of all nests (Dwyer et al. 1996, Raven and Coulson 1997, Nager and O'Hanlon 2016, Blight et al. 2019. Therefore, a potential increase in rooftop nesting is insufficient to explain the reported pattern of long-term continental population declines. Overall, these declines suggest that research should focus on the drivers and related pressures on ecosystems and ecosystem processes resulting from human activities rather than on the species-level interactions and impacts. Ecosystem-based management will require a shift to more integrated and comprehensive management approaches to mitigate anthropogenic contributing factors to threats. Despite the general complexity of predicting the species-specific outcomes of interacting biotic and abiotic factors, diet is understood to be an important driver in gull population dynamics (Anderson et al. 2016, Bennett et al. 2017. In particular, the population-level implications of the use of open landfills (Kadlec and Drury 1968, Coulson 2015, Ackerman et al. 2018, Seif et al. 2018 and fisheries discards (Annett and Pierotti 1999, Bennett et al. 2017, Foster et al. 2017 by gulls have long been debated, especially regarding the extent to which these food subsidies contributed to the increase in gull numbers in the early 20th century (Drury and Kadlec 1974). Nevertheless, major shifts in diet since the 1990s have been recorded in Larus spp. These shifts include a switch from wildcaught prey to novel anthropogenic food sources , Ramírez et al. 2015, Bond 2016) and from marine prey to terrestrial-sourced foods , Bond 2016, Ackerman et al. 2018, Maynard and Ronconi 2018 or lower trophic-level marine prey species (Farmer andLeonard 2011, Ronconi et al. 2014). Importantly, studies in the latter group do not systematically capture whether diet items are from anthropogenic or nonanthropogenic sources.
The number of rooftop-nesting gulls has grown in Europe (Spelt et al. 2019) and may be growing in North America , Kroc 2018, suggesting that gulls are increasingly interacting with humans. This growing overlap, especially at a time when gull numbers are declining overall, heightens the need to understand how prey selection and availability influence gull distribution and population dynamics. This overlap is particularly true of industries that may inadvertently supply prey to gulls or reduce natural prey abundance or availability, especially when industries have options to mitigate such practices. Anthropogenic sources of prey commonly include fishing, agricultural, and waste management industries. In the fishing industry, harvesting, processing, and packaging activities produce bycatch, bait, and waste that become available to gulls when discarded at sea. Likewise, gulls have been known to glean crop fields and access uncontained waste from animal husbandry operations. Agricultural activities such as field tilling and fertilizer application also enhance the availability of naturally occurring prey, including soil-dwelling invertebrates. Finally, the waste management industry, including landfills, transfer stations, and composting facilities, imperfectly contains its product during operations.
The effect of anthropogenic food availability on Larus spp. prey selection has been investigated using stable isotope analysis since Hobson's (1987) work, and modern studies now frequently incorporate GPS tracking with this and other traditional methods of diet description, such as pellet analysis (Steenweg et al. 2011, Ronconi et al. 2014, Enners et al. 2018. One standing challenge is to link prey items to their respective sources, especially in human-affected landscapes. For example, pellet or regurgitate sampling may be used to identify the bones of Atlantic mackerel (Scomber scombrus), but it cannot indicate whether the fish was caught in the wild, scavenged near a wharf, or taken from a fishing vessel at sea (Steenweg et al. 2011, O'Hanlon et al. 2017. Similarly, a fish's trophic level and association with a marine vs. freshwater environment may be inferred using stable isotope analysis, but the point source of those fish remain unknown. GPS tracking allows for detailed analysis of spatiotemporal patterns in travel outside the colony, but cannot be used to confirm which foraging bouts led to successful prey capture. Combining multiple methods to observe foraging behavior is therefore necessary to interpret data about the contribution of anthropogenic food subsidies to bird diet.
Our aim was to quantify the relative use of anthropogenic food sources by breeding Herring Gull in a rapidly changing ecosystem: the Bay of Fundy, Canada. Understanding the contribution of anthropogenic resources to the gull's diet can deepen the understanding of gull-human interactions within an ecosystem context and may help to uncover the role that human activities play in gull distribution and population dynamics.

Study area
The Gulf of Maine/Bay of Fundy system ( Fig. 1) has large, shallow areas uniquely influenced by strong tidal currents and coastal inputs, producing one of the world's richest and most diverse marine biota (Parker 2012). From 2005 to 2015, sea surface temperatures in the Gulf of Maine warmed faster than 99.9% of the global ocean (Pershing et al. 2015). This warming, combined with a highly intensive fishing industry, has caused widespread changes to the entire marine food web, from phytoplankton to top predators such as cod (Gadus morhua; Pershing et al. 2015). As highly visible top predators, marine-dependent seabirds such as terns (Sterna spp.) and Atlantic Puffin (Fratercula arctica) have been used as biological indicators to monitor bottom-up changes occurring in the Gulf of Maine Devlin 2003, Gaston et al. 2009). Because gulls are unique among seabirds in their success in accessing terrestrial resources, changes in measures such as reproductive success and adult survival are likely to lag in gulls compared to other marine-dependent indicator species. However, a rapid, extreme shift in food resources should be reflected quickly in studies of gull diet and foraging movements and, ultimately, may be detected as local or regional declining nest numbers if gulls do not compensate for the ecosystem changes. Kent Island,New Brunswick,Canada (44.581° N,66.756° W) was the largest Herring Gull colony in the Gulf of Maine/Bay of Fundy in the 1930s-1940s, with a population > 20,000 pairs, and probably served as a source for new colonization by gulls to the east and south during that period (Cannell and Maddox 1983). By 1983, the colony declined to 5000 nesting pairs, and then declined further to 1400 by 1989, possibly because of egg collection by local hunters (Hébert 1989). The latest estimates of approximately 6000 pairs in 2001 (Ronconi and Wong 2003) and approximately 3000 pairs in 2015 (Bennett et al. 2017) show that the population rebounded slightly after Hébert's (1989) study but may be declining again . By contrast, nesting habitat on Brier Island, Nova Scotia, Canada (44.275° N,66.345°W ) became available only after drainage of the island's bog in 1958. Herring Gull and Great Black-backed Gull colonized Brier Island thereafter, and the number of nesting pairs grew from 0 in 1970 to > 2400 by 2013 (Toms 2015). At a total colony size of 5100 breeding pairs, Brier Island was the largest Herring Gull colony in Atlantic Canada in 2016. Kent Island and Brier Island are located 36 km apart, so there could be high spatial overlap in the foraging areas around each colony Nisbet 1972, Thaxter et al. 2012).
Fieldwork was conducted during the breeding season (May-August) on Kent Island in 2015 and on Brier Island in 2014-2015 ( Fig. 1). Adults were captured during incubation on nests containing three-egg clutches using walk-in traps. Basic morphometric measurements were taken to determine the sex of sampled individuals according to the discriminant function developed by Robertson et al. (2016). An 80-μL sample of blood was drawn from the tarsal or brachial vein of captured individuals, and the distal 3 cm from the oldest feather among the innermost primaries, P1-P3, was clipped (the significance of these specific samples is described below; Steenweg et al. 2011). On a subset of birds sampled for stable isotopes, GPS loggers were deployed using a leg-loop harness made of Teflon ribbon following Mallory and Gilbert (2008). Two logger types were used, including batterypowered i-GotU (MobileAction Technology, 17-g logger with waterproofing + 2-g harness, N = 15, Brier Island only) and solarpowered HARRIER-L (Ecotone Telemetry, 15-g logger + 2-g harness, N = 26, both islands). Loggers were programmed to record speed and latitude-longitude coordinates at intervals of 5 or 15 min for i-GotU and Ecotone loggers, respectively. For birds equipped with i-GotU loggers, recapture was attempted for data retrieval approximately 2 wk after logger deployment. No recapture was attempted (or required) at nests where the bird carried a solar-powered logger; instead, a long-range antenna was placed at the edge of both colonies to receive GPS data remotely from tagged gulls. We returned to each colony briefly in 2016 to set up a long-range antenna to collect GPS data from any previously tagged gulls returning to the colonies.

Analysis of prey sources using GPS data
In total, 11 solar-powered and 5 battery-powered GPS loggers contributed data from Brier Island, Nova Scotia, and 11 solarpowered loggers contributed data from Kent Island, New Brunswick (Table 1; Appendix 1 Table). Results from a discriminant function of morphometric measurements  suggested that 14 of the tagged individuals were female, 11 were male, and 2 could not be determined. To annotate the GPS records with breeding stage categories, average hatch dates were calculated for both colonies based on nest monitoring data (2015) and observations made regarding nesting activity while trapping (e.g., average clutch size, frequency of hatching; 2014-2015). Average hatch dates were then used to establish incubation periods (beginning 6 wk prior to average colony hatch) and chick rearing periods (ending 6 wk after average colony hatch) in each year for each colony. We reviewed raw GPS files to address data collection gaps, duplicate records, and erroneous records (e.g., speed > 100 km/h or timestamp outside the deployment period). In years of logger deployment (Brier Island 2014-2015, Kent Island 2015), we excluded GPS records from the first two days of tracking data to avoid analyzing the abnormal sedentary behavior that was observed immediately following gull capture and tagging. Of the loggers we analyzed, no significant long-term effects of capture or tagging on gull movement or productivity were observed during daily nest monitoring visits (Shlepr 2017). Our failure to detect a logger effect is in accordance with previous gull tagging studies with harness attachment (Thaxter et al. 2016, Kavelaars et al. 2018. Because our aim was to analyze foraging destinations in particular, we developed a filtering protocol to exclude GPS records in the immediate vicinity of the nesting area (< 500 m from an individual's nest) as well as locations that were likely to have been recorded on the way to or from foraging or roosting destinations. To remove these in-transit points, we assumed that gulls would come to rest at (or near) foraging and roosting destinations, and we filtered the GPS data to retain only records where the GPS logger recorded a speed of < 2.0 km/h (Fig. 1). The 2 km/h cutoff is equivalent to gull walking speed or slower and is similar to the data-filtering strategy used in previous studies (Shamoun-Baranes et al. 2011, Maynard andRonconi 2018). We were initially concerned that this decision would bias our results against marine foragers, which are more likely than terrestrial or coastal foragers to fly steadily while foraging because the marine environment is patchy (Weimerskirch 2007) and gulls are visual hunters (Nisbet et al. 2017). However, marine points in our quality-controlled data set were uncommon, and the proportion of points located in the marine environment changed only minimally after filtering the data set for near-stationary records (14.2% of off-colony locations prior to filtering for speeds < 2 km/h; 9.2% of off-colony locations after filtering for speeds < 2 km/h).
Roosting destinations were effectively separated from foraging destinations via identification of known anthropogenic sources of prey using on-ground site observations (in 2015) and GPS results (in 2014-2016). To categorize locally relevant prey sources, polygons of potential point sources of anthropogenic prey items were manually delineated using latitude and longitude coordinates in Google Earth and then entered into R. These categories were: "fisheries" (i.e., active town wharfs, seafood processing and packaging plants, and salmon aquaculture pens), "mink" (i.e., active mink [Neovison vison] farms, where gulls could access mink feed that had fallen through animal cages, and animal composting facilities), and "other" (i.e., active waste transfer stations, solid waste parks, and municipal compost sites). All offcolony GPS records located outside of these "known anthropogenic" polygons defaulted to a category called "unassigned". Unassigned locations represented all nonanthropogenic sources of food, unknown anthropogenic sources of food, and nonforaging roosting locations. To improve explanatory potential for records in the unassigned category, unassigned locations were split into three generic geographic classifications: "terrestrial" (> 3 km inland), "coastal" (± 3 km of shoreline), and "marine" (> 3 km offshore).
In R, we delineated unique off-colony trips by assigning a trip ID to all GPS records between the time an individual left the vicinity of its nest (> 500 m from nest) to the time it reentered its nest area (< 500 m from nest). The mean distance from nest for each offcolony trip was calculated and then compared using analysis of variance (ANOVA). Additionally, the proportion of GPS records in known anthropogenic categories was averaged for all available years Brier Island, 2015) and compared at the colony level using ANOVA. Statistical analyses were conducted using the "aov" package in R (R Core Development Team 2019).

Analysis of prey selection using stable isotopes
Stable isotope blood and feather sample collection followed procedures described in Steenweg et al. (2011). Because wholebody turnover of blood takes approximately 4 wk in birds (Hobson and Clark 1992a), blood samples drawn at the time of trapping reflect the nest building, egg synthesis (females only), and incubation stages of the gull breeding cycle. Additionally, gulls in the Bay of Fundy begin molting their innermost primary feathers during the incubation period (Nisbet et al. 2017). The innermost primaries are synthesized during incubation and are the first to be replaced (i.e., prior to being dropped, the innermost is the oldest feather with a with molt score = 0). Because it was the most consistently available at capture, the tip of the oldest available primary feather was collected for isotope analysis and was assumed to best reflect the individual's diet during the incubation stage one year prior (2013 or 2014 for birds sampled in 2014 and 2015, respectively). Although we were not yet studying prey conditions in the Bay of Fundy in 2013, the oldest primary feather is a useful indicator of local Herring Gull diet because these birds tend to have high interannual site fidelity (Nisbet et al. 2017), which was confirmed at our study sites by the high return of gulls that were tagged (Brier Island: six of nine gulls that were present during chick-rearing in 2014 returned to the colony for at least two consecutive years; Kent Island: 10 of 11 gulls that were present during chick-rearing in 2015 returned to the colony in the following year; Table 1; Appendix 1 Table). The mean molt scores recorded for P1 (the innermost primary feather), P2, and P3 were 0.4, 0.1, and 0.0, respectively (N = 41). Finally, potential prey types were determined by local knowledge, our observations of gull foraging from boats (Bay of Fundy) and inland (Nova Scotia, New Brunswick;-2015, and a regionally focused scan of the literature. We sampled these prey types opportunistically throughout the 2015 field season within the reported mean foraging range of gulls by harvesting or obtaining samples from fishers and mink farmers. Terrestrial prey items, including earthworms (Lumbricus spp.), beetles (family Scarabaeidae), and slugs (class Gastropoda), were sampled and analyzed in the lab but were ultimately excluded from isotope mixing models when preliminary results showed that both the δ 13 C and δ 15 N values of this group were considerably lower than any of the gull tissue samples (mean ± standard deviation; δ 13 C: −27.9 ± 2.5‰; δ 15 N: 5.5 ± 1.7‰; see Results for more information).
Whole blood samples (N = 90; Table 1) were dried in an oven at 60°C for 48 h and then ground by hand using a mortar and pestle. Feathers (N = 73; Table 1) were cut, washed with a 2:1 chloroform methanol solution to remove surface impurities, and then placed in a 60°C oven for 1-2 h to dry. Prey source samples were dried in a 60°C oven, ground, and, if relevant, treated with hydrochloric acid (to dissolve carbonate shells) or a 2:1 chloroform methanol solution (to extract lipids). Samples were weighed and submitted to the Stable Isotopes in Nature Laboratory (Fredericton, New Brunswick, Canada) for chemical analysis. Continuous flowisotope ratio mass spectrometry was used for stable-isotope analysis of δ 13 C and δ 15 N. Values were calibrated to International Atomic Energy Agency primary standards: Vienna Pee Dee Belemnite for δ 13 C, and atmospheric nitrogen for δ 15 N. Experimental error was determined to be 0.1‰ for both isotopes based on repeat analysis of an in-house standard, muskellunge liver tissue (Jardine and Cunjak 2005).
Discrimination factors for whole blood (carbon: −0.1 ± 0.7‰; nitrogen: 2.8 ± 0.4‰) and primary feathers (carbon: 1.2 ± 0.7‰; nitrogen: 4.1 ± 0.4‰) were taken from Herring Gull studies by Steenweg et al. (2011) and Ronconi et al. (2014), who averaged the discrimination factors from experiments with Ring-billed Gull (Hobson and Clark 1992b) and other piscivorous seabirds (Bearhop et al. 2002, Cherel et al. 2005, Becker et al. 2007). Potential prey items were analyzed for δ 13 C and δ 15 N values and then grouped according to isotopic and biological similarity to reduce the number of prey categories considered in mixing models. All statistical analyses were run in R (R Core Development Team 2019). We used ANOVA to test whether year explained any of the variation in δ 13 C and δ 15 N of gull tissues from Brier Island, where we had two years of sampling (2014)(2015), and we also used ANOVA to test whether colony-level differences in δ 13 C and δ 15 N were significant. Likewise, year and colony were included in the Bayesian mixing model as random and fixed effects, respectively, using the package MixSIAR (Stock et al. 2018).

Estimating anthropogenic resource use with combined methods
Gulls can obtain some prey from wild and anthropogenic sources, so both aspects (prey source and type) must be known to quantify anthropogenic resource use as a proportion of the whole diet. Because no single method directly measures both aspects, we devised a way to objectively visualize information overlaps between our GPS (prey source) and stable-isotope (prey type) results. We restricted this analysis of anthropogenic resource use to the incubation period because the stable-isotope analysis used tissue samples that only reflect diet in the early stages of breeding.
First, by colony, we used the GPS results (Table 2) to calculate the mean percentage of tagged gull locations that fell into each prey source category (incubation only, all available years). This process provided a snapshot of where gulls tend to spend time foraging. Next, by colony, we used the stable-isotope results to isolate the percentage estimates of prey type by category (blood only, averaged years when available; alternate analysis using feather sample results can be found in Fig. S1 in Appendix 1). This process provided a snapshot of what prey types gulls tend to consume. Finally, to combine this information on where (GPS results) and what (stable-isotope results) gulls tend to eat, we averaged estimates from every ecologically viable combination of prey source categories (GPS results) and prey type categories (stable-isotope results). For example, we averaged the use scores from the logical combination of high trophic marine prey (stableisotope category) in known anthropogenic fisheries sites (GPS category), but we did not average the use scores from the illogical combination of high trophic marine prey (stable-isotope category) in unassigned terrestrial sites (GPS category), labeling illogical prey source-type combinations as such. Even though these wild and anthropogenic resource use values are the result of averaging pairs of percentage estimates, we acknowledge that we do not account for variation around means (standard deviation [SD] in the GPS analysis; 95% confidence interval [CI] in stableisotope mixing models) by displaying resource use values as unitless scores.
The scores produced by this exercise can be used objectively to weigh the importance of a specific wild or anthropogenic resource relative to other resources available to gulls. These scores theoretically range from 0 (0% mean estimate of a prey source category and 0% mean estimate of a prey type category) to 100 (100% mean estimate of a prey source category and 100% mean estimate of a prey type category). As such, higher resource use scores indicate that gulls consumed high levels of a prey type that is obtained from a habitat type where gulls spent a lot of time foraging, therefore providing greater confidence that a particular wild or anthropogenic resource is an important part of the gull diet. Conversely, a low resource use score could indicate: (1) a mismatch between the amount of time gulls spent at a habitat type and what they were found to have consumed, which means we have low confidence about the importance of a particular resource; or (2) high confidence that gulls do not depend on a particular resource because associated prey source and prey type estimates were low in both the GPS and stable-isotope analyses. To target foraging and roosting destinations in the analysis, GPS data were filtered to include only off-colony locations where speed was < 2 km/h. Potential anthropogenic prey sources were defined a priori, and GPS records falling within these polygons were classified as "known anthropogenic" sources. GPS records falling outside anthropogenic source polygons defaulted to an "unassigned" category and were classified as terrestrial (> 3 km inland), coastal (± 3 km of shoreline), or marine (> 3 km offshore). See Table 1

RESULTS
Of the 115,723 GPS fixes recorded during the incubation and chick-rearing periods, 51,425 (44%) were recorded during 4255 independent off-colony trips at speeds < 2 km/h; these were the data used in subsequent analyses. These data came from 16 individuals from Brier Island (N = 11 solar loggers, mean ± SD: 124 ± 81 d/bird tracked, 2014-2016; N = 5 battery loggers, mean ± SD: 19 ± 5 d/bird tracked, 2014 incubation only) and 11 individuals from Kent Island (N = 11 solar loggers, mean ± SD: 113 ± 30 d/bird tracked, 2015-2016; Table S1 in Appendix 1). For the stable-isotope analysis, we collected blood and feather samples from all tagged individuals and additional captures at both colonies ( Table 1).

Analysis of prey sources using GPS data
Gulls from Brier Island tended to travel further away from their colony (N = 16, mean ± SD: 38 ± 81 km) than gulls from Kent Island (N = 11, mean ± SD: 21 ± 24 km; ANOVA P < 0.001; Fig.  1). Brier Island gulls spent a range of 15-85% (mean ± SD: 64 ± 20%) of their off-colony time at rest in a known anthropogenic site, which is significantly more than Kent Island gulls, which spent a range of 3-82% (mean ± SD: 35 ± 24%) of their time in the same category (ANOVA P = 0.016; Table 2). For gulls of both colonies, most GPS points that were located outside a known anthropogenic site occurred in coastal (± 3 km from shoreline) habitat (Brier Island: 24 ± 16%, Kent Island: 58 ± 23%), which may reflect intertidal foraging or time spent roosting without foraging. Gulls overall spent little time foraging or resting in the marine environment (Brier Island: 7 ± 11%, Kent Island: 7 ± 12%; Fig. 1, Table 2).
As is reflected by the large standard deviations in the known anthropogenic site category (Table 2), individual use of known anthropogenic sites varied widely, especially among gulls from Kent Island (Fig. 2). The use of fisheries (including town wharfs, fish processing plants, and salmon aquaculture rings) varied widely among individuals and could not be explained by colony of origin (ANOVA P = 0.272; Fig. 2). Conversely, individuals appeared to use mink prey sources either heavily or not at all. Colony appeared to be the biggest determinant for this food source; 11 of 16 gulls from Brier Island spent at least 25% of their off-colony time at mink sites compared to 1 of 11 Kent Island gulls (ANOVA P < 0.001; Fig. 2). Only three gulls (all from Kent Island) ever visited an alternative anthropogenic site (including municipal and regional landfills), and that site type (i.e., other) constituted ≤ 4% of records for each of the three individuals.
Overall, visits to municipal or regional landfills constituted 0.3% of all GPS locations in the Bay of Fundy.

Analysis of prey selection using stable isotopes
Potential prey items were homogenized, analyzed for mean δ 13 C and δ 15 N and then grouped according to isotopic and biological similarity, resulting in five prey categories that served as the basis for subsequent MixSIAR analyses (Table 3). Isotopic separation was best along the δ 15 N gradient (category averages ranged from δ 15 N of 5.1‰ to 11.8‰) and overlapped more extensively along the δ 13 C gradient (category averages ranged from δ 13 C of −21.5‰ to −15.9‰; Table 3, Fig. 3). An increase in δ 15 N between samples is typically interpreted to reflect an animal's increase in trophic position in the food web, whereas an increase in δ 13 C between samples is typically indicative of an animal that lives in a more marine (vs. in-shore or terrestrial) environment.
Year also had a significant effect on δ 13 C (P < 0.001) for the Brier Island colony, where two breeding seasons of data were available for analysis (2014: −19.5 ± 0.3‰, 2015: −19.0 ± 0.4‰). We did not detect a significant effect of year on δ 15 N values derived from Brier Island (P = 0.365) blood samples (Fig. 3A).   Table  4). Mink and industrial discards, the two solely anthropogenic categories, were moderately important in both colonies, constituting mean contribution estimates up to 21% and 19%, respectively (Table 4). The low trophic marine prey category, i.e., prey that gulls catch wild, was also of moderate importance in both colonies, constituting a mean contribution estimate of up to 18% of gull diet at Brier Island and up to 12% of gull diet at Kent Island (Table 4). Isotope analysis cannot be used to determine what proportion of prey from the high trophic marine and crabs categories is anthropogenically sourced.  Table 3). Prey categories were adjusted by discrimination factors relative to bird tissue samples. Brier Island: N blood = 56, feathers = 59; Kent Island: N blood = 34, feathers = 14.

Estimating anthropogenic resource use with combined methods
For the Brier Island colony, resource use (prey source-prey type) scores ranged from 9 to 29, with 50% of them ranging from 17 to 23 (Fig. 4A). Resources with the highest scores include known anthropogenic mink farms-mink feed/carcasses (score = 29), unassigned coastal habitat-high trophic marine prey (28), and known anthropogenic fisheries-high trophic marine prey (score = 26). Resources with the lowest scores include unassigned marine habitat-low trophic marine prey (9), unassigned marine habitatindustrial discards (score = 12), and both logical unassigned terrestrial habitat pairings (scores = 13-14).
Resource use at the Kent Island colony ranged from 5 to 59, with 50% of resource use scores ranging from 7 to 19 (Fig. 4B). All logical unassigned coastal habitat-prey type pairings scored ≥ 35, including the resource with the highest score, unassigned coastal habitat-crab (score = 59). The second highest scoring resource was known anthropogenic fisheries-crab (score = 42). Resources with the lowest scores included industrial discards from known anthropogenic landfills (score = 5) and unassigned terrestrial (score = 5) or marine habitats (score = 6), as well as marine prey from unassigned marine sources (score = 7).

Fig. 4.
Estimation of the relative importance of various wild and anthropogenic food resources available to Herring Gull (Larus argentatus) captured on Brier Island, Nova Scotia (A) and Kent Island, New Brunswick (B), Canada, during the breeding seasons in May-July 2014-2016. Prey type: mean estimates (%) of the proportion of the diet that comprises each of five prey-type categories are the results of a separate stable-isotope analysis (see Table 4). Source: mean estimates (%) of GPS points that were assigned to each of six prey-source categories are the results of a separate GPS analysis (see Table 2). Each score in the table is the average between a prey-type category and a prey-source category; ecologically illogical prey source-prey type combinations are labeled. Scores can be used objectively to weigh the importance of specific wild (green) or anthropogenic (yellow) resources relative to other resources available to gulls. Higher resource use scores (darker cells) indicate that gulls consumed high levels of a prey type that is obtained from a habitat type where gulls spent a lot of time foraging, therefore giving confidence that a particular wild or anthropogenic resource is an important part of the gull diet. 3 Anthropogenic † Standard deviation. ‡ Indicates whether a particular prey item is expected to have been obtained by gulls from the wild, from an anthropogenic location, or when the origin could be either wild or anthropogenic.

Anthropogenic resource use
Together, the GPS and stable-isotope results suggest that prey from anthropogenic sources are likely to be an important part of the diet of herring gulls that nest in the Bay of Fundy, Canada. Gulls from both colonies consumed high levels of crab and high trophic marine prey categories (e.g., Atlantic herring [Clupea harengus], Atlantic mackerel), which are available and accessed by gulls at both known anthropogenic locations (i.e., fisheries). They also consumed either unidentified anthropogenic or nonanthropogenic prey sources.
Our results also provide strong evidence that landfills and terrestrial destinations in general, excluding mink sites, do not provide significant food subsidies to gulls in the Bay of Fundy during the breeding season. In the GPS analysis, the terrestrial categories (known anthropogenic other and unassigned terrestrial) had the lowest or second lowest use scores among all six prey source categories on both islands. In the stable-isotope mixing model results, the industrial discards category, which featured human refuse (compost), received a moderate use score relative to the other four prey type categories for Brier Island gulls and the lowest relative use score for Kent Island gulls. Taken together, the anthropogenic resource use estimate showed that where gulls scored higher in the industrial discards prey type category, the prey source is much likelier to be from known anthropogenic fisheries or unassigned coastal habitats than from an inland location such as a landfill. Finally, we reiterate that we ultimately excluded terrestrial arthropod prey types from the stable-isotope mixing models because their δ 13 C and δ 15 N values were extreme compared to all gull tissue samples.
On average, gulls from Brier Island travelled further from their nests than gulls from Kent Island to access locations where mink sites inadvertently supplied food in the form of fallen mink feed or improperly disposed mink carcasses. Approximately half of the gulls from Brier Island spent the majority of their off-colony time at sites where anthropogenic food is known to be available, particularly mink sources. Even though the mink prey category received the second highest use estimate by stable-isotope mixing models, it received the highest overall score in the anthropogenic resource use estimate (i.e., combined GPS and stable-isotope results). Gulls nesting on Kent Island relied less on anthropogenic food types overall, although approximately one-third of Kent Island birds' off-colony movements still fell into one of the anthropogenic source categories, predominantly known anthropogenic fisheries. The stable-isotope results were less definitive because all prey type categories had wide confidence intervals and no single category was estimated to be a majority of gulls' diet at either colony.
Blood and feather samples produced similar results in the stableisotope analysis, with the strongest pattern being clear colonylevel separation in δ 13 C rather than trophic level (δ 15 N). Because the innermost primary feathers reflect the diet at the time of feather growth, these samples reflect the individuals' diets during incubation (approximately May) of the year prior to sample collection, whereas blood, which turns over on a monthly basis, reflects the individual's diet during the incubation period of sample collection (Hobson and Clark 1992a). Therefore, the observed consistency in blood and feather model outputs over both tissue types and years suggests consistency in colony-level diets across at least two years.

Study limitations and future research
We have noted differences in the estimates produced by GPS and stable-isotope analyses of gulls' reliance on anthropogenic food sources. From the GPS perspective, this difference could indicate that a large portion of gulls' off-colony stopovers are nonforaging destinations (i.e., roost sites) or are foraging destinations where prey is infrequently or unpredictably available. Alternatively, the estimates produced by stable-isotope mixing models may be influenced by the way we combined potential prey items into prey Table 4. Mean percentage (and 95% confidence interval) contribution of each prey type category as estimated by Bayesian mixing models. Dietary estimates were derived from δ 13 C and δ 15 N values from whole blood and primary feathers sampled from nesting Herring Gull (Larus argentatus) at Brier Island, Nova Scotia (N blood = 56, feathers = 59) and Kent Island,New Brunswick (N blood = 34,feathers = 14), Canada. Average values are shown for Brier Island only where sampling took place over two years (see Table 1).  (Phillips et al. 2005) or by uncertainty about discrimination factors (Bond and Diamond 2011). Such limitations are especially troublesome when analyzing generalist predators such as the Herring Gull, for which the list of possible prey sources spans nearly the entire range of δ 13 C and δ 15 N values that can be found in the study area. Additionally, in eastern Canada, ground Atlantic herring is often used as the primary ingredient in both mink feed and aquaculture salmon feed, so these two prey items, in combination with wild-caught herring, are indistinguishable to isotope mixing models. Moreover, Atlantic herring and other high trophic-level fish are used as bait in local lobster fisheries, which are active in May and June, and herring gulls frequently scavenge for discarded bait around vessels at sea.
We circumvented these challenges of stable-isotope analysis by combining the results with results from the GPS analysis. Because the GPS analysis suggests that gulls spend much of their time at coastal or inland sites and little time offshore, we can infer that gulls in the Bay of Fundy are more likely to obtain Atlantic herring, and thus high trophic marine prey, through anthropogenic sources such as fish processing plants or mink (Brier Island) and salmon (Kent Island) feed. Moreover, the availability of herring as lobster bait in the marine environment likely declines after the lobster season's annual closure (May 31, Brier Island region; June 29, Kent Island region), making terrestrial sources of herring (i.e., mink feed, aquaculture salmon feed) the most likely sources during the gulls' chick-rearing stage (early June to mid-August).
It is possible that our GPS analysis is biased against at-sea foraging. We filtered GPS records to include only near-stationary speeds (< 2 km/h) after the preliminary observation that gull attendance in the marine environment truly was low relative to coastal and inland areas. However, possible specialists in marine foraging deserve more attention in future studies. Another possible biasing issue relates to our decision to analyze GPS data from Brier Island in 2014, including data from five gulls with an alternate logger type (battery instead of solar), even though Kent Island data collection did not begin until 2015. Including the additional year (2014) and logger type on Brier Island maximized what could be learned about interannual and intracolony variation in prey sources but may have biased the colony-level averages that were highlighted in our combined GPS-isotope analysis. We believe this bias to be insignificant because categorical averages and variation were consistent between years in both the GPS and stable-isotope analyses, even in the incubation stage of 2014, when the alternate logger type contributed data. Nevertheless, further study of individual variation should be investigated.

Evolutionary implications
Though the Brier Island and Kent Island colonies are located 36 km apart (well within the typical foraging range of these species), we found colony-level differences in both foraging areas and prey types obtained in each category. These results are consistent with results from gull tagging studies in coastal California (Shaffer et al. 2017) and Germany (Corman et al. 2016, Enners et al. 2018, which identified diet and foraging habitat differences between gulls from different colonies within the same region. Optimal foraging theory posits that this apparent niche partitioning is an adaptive response to intraspecific competition (Kotler andBrown 2007, Corman et al. 2016).
It is rarely known whether individual variation in gull foraging behaviors is the result of behavioral plasticity, local adaptation, or both (Maynard and Ronconi 2018). A single experimental translocation study of the Lesser Black-backed Gull in northern Europe and western Asia suggests that plasticity is the primary driver of observed geographic variation in diet (van Toor et al. 2017); however, more work is needed to understand the implications of individual variation for population dynamics and species conservation. If individuals are truly generalist opportunists, they may be able to compensate for the loss of one food type by switching to another (Kotler and Brown 2007). However, if gulls are locally adapted (i.e., specialized) or are limited by the types of prey available to them, rapid changes in a single prey type may have dramatic consequences to gull survival and reproduction. For example, Kent Island gulls have shown long-term variation in egg size, which tracks the availability of fisheries-related food sources over a 28-yr period (Bennett et al. 2017). In the near term, food shortages may cause individuals to skip breeding seasons until prey return to adequate abundance or availability so as to elevate adult condition and warrant investment in breeding efforts. Long term, individuals may emigrate to new breeding areas in search of more abundant or available prey. Gull responses on both time scales may have contributed to the decline in nest numbers that has been witnessed on Kent Island since the mid-20th century Wong 2003, Bennett et al. 2017).
Variation in diet at the individual and colony levels may ultimately explain differences in annual productivity and colony return rates if the energetic gains from these prey types also vary, either because different prey types vary in quality Annett 1991, Weiser andPowell 2010) or because different prey types demand alternative foraging strategies that may trade off with nest attendance (Hunt 1972, Shaffer et al. 2017 or other behavioral factors (O'Hanlon et al. 2017). Future research should therefore focus directly on the relationships between prey choice and productivity at the nest level, which may reveal further understanding of population dynamics within colonies.

Conservation implications
The standing hypothesis explaining the widespread, long-term declines in Larus species as observed in eastern North America is that gulls have been unable to adapt to the decrease in total prey availability due to improvements in waste management and the collapse of multiple major marine fisheries (MacKinnon and Kennedy 2014, Anderson et al. 2016. The literature has focused on a decreased volume of available prey, but changes in the overall quality of the diet, especially in relation to artificial prey types, may play an equally important role. Regardless, a critical information gap has been in understanding whether pivotal changes to diet are related to changing anthropogenic sources, natural sources, or both. If the total prey decrease is caused primarily by the elimination of human industrial practices that inadvertently supply food, then we might predict that earlier peak populations were "artificially" inflated and are predictably returning to sustainable, presubsidy breeding numbers. If, however, the reduction in total prey available to gulls comes partly or entirely from reduced natural food sources, then wildlife managers must be prepared to mitigate unsustainable declines in gull numbers as anthropogenic sources of prey are eliminated. Our results indicate that crab and high trophic marine prey from both wild and anthropogenic sources are key components of the diet of gulls at the two largest colonies in the Bay of Fundy. Wildlife managers should therefore exercise caution in assuming that future elimination of anthropogenic food subsidies will result in a gull population decline to sustainable, preindustrial colony sizes. Elimination of mink farms as a prey source may very well lead to a local decline in gull numbers (and is perhaps already being observed), but it is unclear whether gulls, especially those from Brier Island, have access to wild food types that they are not currently using. The future of gulls within foraging range of modern mink farm activity centers therefore remains to be seen.

Conclusion
GPS tracking and stable-isotope analysis are two complementary approaches to understanding gull foraging ecology, and each has strengths and biases. Along with site-level field observations pertaining to the natural history of the species, these methods allowed us to better explore sources of variation in gull diets, especially in relation to these birds' relative reliance on anthropogenic sources of food. Understanding how human activities influence gull behavior permits improvements in the design and implementation of wildlife management strategies that mitigate undesired anticipated effects on gulls and gull populations and on the ecosystems to which they belong. Ultimately, this information can be used to inform the adaptive management of gulls, which is especially vital in the face of global climate change and other causes of intensive landscape alteration.