Carry-over effects of seasonal migration on reproductive success through breeding site retention in a partially migratory bird.

Understanding the maintenance and dynamics of phenotypic polymorphisms requires unpicking key ecological mechanisms shaping the fitness costs and benefits of expressing alternative phenotypes, generating selection. Seasonal migration versus year-round residence expressed in partially migratory populations represents one common polymorphism that can experience strong selection through differential reproductive success. Yet, key hypothesised pathways that could generate such selection remain to be empirically tested. One hypothesis is that migratory tactics affect subsequent reproductive success through carry-over effects on breeding site retention and resulting breeding dispersal. By remaining in breeding areas all year round, residents could retain their preferred breeding site between years, and consequently have higher reproductive success. Conversely, migrants that escape harsh non-breeding season conditions could return in better condition, with high resource holding potential, and outcompete residents to retain their site. Such effects could further depend on migration timing and vary between years. Yet, such pathways have not been quantified, precluding empirical parameterisation of partial migration theory. We used 4 years of breeding and non-breeding season data from partially migratory European shags (Gulosus aristotelis) to test whether the three most frequent migratory tactics in this population (full resident, early migrant departing soon after breeding, and late migrant departing in late autumn) differed in their breeding site retention; whether site retention predicted reproductive success; and hence whether effects of migratory tactic on reproductive success were explicable through site retention. Overall, residents were much more likely to retain their breeding site between years than both early and late migrants, and site retention was associated with increased reproductive success. Yet, these effects varied somewhat among years: late migrants were always least likely to retain their site but had variable relative reproductive success. Path analyses revealed that effects of migratory tactic on reproductive success were only partly attributable to breeding site retention. These results indicate that multiple mechanisms underlie reproductive selection on migratory tactics, potentially contributing to maintaining behavioural polymorphisms. Yet, the clear associations between migratory tactics and local breeding dispersal reveal that these movements can be strongly interlinked across seasons, shaping overall spatioseasonal dynamics in partially migratory systems.

Hierarchical Variation in Phenotypic Flexibility across Timescales and Associated Survival Selection Shape the Dynamics of Partial Seasonal Migration.

AbstractPopulation responses to environmental variation ultimately depend on within-individual and among-individual variation in labile phenotypic traits that affect fitness and resulting episodes of selection. Yet complex patterns of individual phenotypic variation arising within and between time periods, as well as associated variation in selection, have not been fully conceptualized or quantified. We highlight how structured patterns of phenotypic variation in dichotomous threshold traits can theoretically arise and experience varying forms of selection, shaping overall phenotypic dynamics. We then fit novel multistate models to 10 years of band-resighting data from European shags to quantify phenotypic variation and selection in a key threshold trait underlying spatioseasonal population dynamics: seasonal migration versus residence. First, we demonstrate substantial among-individual variation alongside substantial between-year individual repeatability in within-year phenotypic variation ("flexibility"), with weak sexual dimorphism. Second, we demonstrate that between-year individual variation in within-year phenotypes ("supraflexibility") is structured and directional, consistent with the threshold trait model. Third, we demonstrate strong survival selection on within-year phenotypes-and hence on flexibility-that varies across years and sexes, including episodes of disruptive selection representing costs of flexibility. By quantitatively combining these results, we show how supraflexibility and survival selection on migratory flexibility jointly shape population-wide phenotypic dynamics of seasonal movement.

Foraging habitat and site selection do not affect feeding rates in European shags.

During feeding trips, central-place foragers make decisions on whether to feed at a single site, move to other sites and/or exploit different habitats. However, for many marine species, the lack of fine-resolution data on foraging behaviour and success has hampered our ability to test whether individuals follow predictions of the optimal foraging hypothesis. Here, we tested how benthic foraging habitat usage, time spent at feeding sites and probability of change of feeding sites affected feeding rates in European shags (Gulosus aristotelis) using time-depth-acceleration data loggers in 24 chick-rearing males. Foraging habitat (rocky or sandy) was identified from characteristic differences in dive patterns and body angle. Increase in body mass was estimated from changes in wing stroke frequency during flights. Bout feeding rate (increase in body mass per unit time of dive bout) did not differ between rocky and sandy habitats, or in relation to the order of dive bouts during trips. Bout feeding rates did not affect the duration of flight to the next feeding site or whether the bird switched habitat. However, the likelihood of a change in habitat increased with the number of dive bouts within a trip. Our findings that shags did not actively move further or switch habitats after they fed at sites of lower quality are in contrast to the predictions of optimal foraging theory. Instead, it would appear that birds feed probabilistically in habitats where prey capture rates vary as a result of differences in prey density and conspecific competition or facilitation.

Increased parental effort fails to buffer the cascading effects of warmer seas on common guillemot demographic rates.

Climate warming can reduce food resources for animal populations. In species exhibiting parental care, parental effort is a 'barometer' of changes in environmental conditions. A key issue is the extent to which variation in parental effort can buffer demographic rates against environmental change. Seabirds breed in large, dense colonies and globally are major predators of small fish that are often sensitive to ocean warming. We explored the causes and consequences of annual variation in parental effort as indicated by standardised checks of the proportions of chicks attended by both, one or neither parent, in a population of common guillemots Uria aalge over four decades during which there was marked variation in marine climate and chick diet. We predicted that, for parental effort to be an effective buffer, there would be a link between environmental conditions and parental effort, but not between parental effort and demographic rates. Environmental conditions influenced multiple aspects of the prey delivered by parents to their chicks with prey species, length and energy density all influenced by spring sea surface temperature (sSST) in the current and/or previous year. Overall, the mean annual daily energy intake of chicks declined significantly when sSST in the current year was higher. In accordance with our first prediction, we found that parental effort increased with sSST in the current and previous year. However, the increase was insufficient to maintain chick daily energy intake. In contrast to our second prediction, we found that increased parental effort had major demographic consequences such that growth rate and fledging success of chicks, and body mass and overwinter survival of breeding adults all decreased significantly. Common guillemot parents were unable to compensate effectively for temperature-mediated variation in feeding conditions through behavioural flexibility, resulting in immediate consequences for breeding population size because of lower adult survival and potentially longer-term impacts on recruitment because of lower productivity. These findings highlight that a critical issue for species' responses to future climate change will be the extent to which behavioural buffering can offer resilience to deteriorating environmental conditions.

Raising offspring increases ageing: Differences in senescence among three populations of a long-lived seabird, the Atlantic puffin.

Actuarial senescence, the decline of survival with age, is well documented in the wild. Rates of senescence vary widely between taxa, to some extent also between sexes, with the fastest life histories showing the highest rates of senescence. Few studies have investigated differences in senescence among populations of the same species, although such variation is expected from population-level differences in environmental conditions, leading to differences in vital rates and thus life histories. We predict that, within species, populations differing in productivity (suggesting different paces of life) should experience different rates of senescence, but with little or no sexual difference in senescence within populations of monogamous, monomorphic species where the sexes share breeding duties. We compared rates of actuarial senescence among three contrasting populations of the Atlantic puffin Fratercula arctica. The dataset comprised 31 years (1990-2020) of parallel capture-mark-recapture data from three breeding colonies, Isle of May (North Sea), Røst (Norwegian Sea) and Hornøya (Barents Sea), showing contrasting productivities (i.e. annual breeding success) and population trends. We used time elapsed since first capture as a proxy for bird age, and productivity and the winter North Atlantic Oscillation Index (wNAO) as proxies for the environmental conditions experienced by the populations within and outside the breeding season, respectively. In accordance with our predictions, we found that senescence rates differed among the study populations, with no evidence for sexual differences. There was no evidence for an effect of wNAO, but the population with the lowest productivity, Røst, showed the lowest rate of senescence. As a consequence, the negative effect of senescence on the population growth rate (λ) was up to 3-5 times smaller on Røst (Δλ = -0.009) than on the two other colonies. Our findings suggest that environmentally induced differences in senescence rates among populations of a species should be accounted for when predicting effects of climate variation and change on species persistence. There is thus a need for more detailed information on how both actuarial and reproductive senescence influence vital rates of populations of the same species, calling for large-scale comparative studies.

Site-dependent regulation of breeding success: Evidence for the buffer effect in the common guillemot, a colonially breeding seabird.

Density-dependent regulation can offer resilience to wild populations experiencing fluctuations in environmental conditions because, at lower population sizes, the average quality of habitats or resources is predicted to increase. Site-dependent regulation is a mechanism whereby individuals breed at the highest quality, most successful, sites, leaving poorer quality, less successful sites vacant. As population size increases, higher quality sites become limiting but when populations decline, lower quality sites are vacated first, offering resilience. This process is known as the 'buffer effect'. However, few studies have tested whether such regulation operates in populations experiencing changes in size and trend. We used data from a population of common guillemots Uria aalge, a colonially breeding seabird, to investigate the relationship between site occupancy probability, site quality and population size and trend. These data were collected at five sub-colonies spanning a 38-year period (1981-2018) comprising phases of population increase, decrease and recovery. We first tested whether site quality and population size in sub-colonies explained which sites were occupied for breeding, and if this was robust to changes in sub-colony trend. We then investigated whether disproportionate use of higher quality sites drove average site quality and breeding success across sub-colony sizes and trends. Finally, we tested whether individuals consistently occupied higher quality sites during periods of decline and recovery. Higher quality sites were disproportionality used when sub-colony size was smaller, resulting in higher average site quality and breeding success at lower population sizes. This relationship was unaffected by changes in sub-colony trend. However, contrary to the predictions of the buffer effect, new sites were established at a similar rate to historically occupied sites during sub-colony decline and recovery despite being of lower quality. Our results provide support for the buffer effect conferring resilience to populations, such that average breeding success was consistently higher at lower population size during all phases of population change. However, this process was tempered by the continued establishment of new, lower quality, sites which could act to slow population recovery after periods when colony size was low.

Opposing effects of spatiotemporal variation in resources and temporal variation in climate on density dependent population growth in seabirds.

Understanding how ecological processes combine to shape population dynamics is crucial in a rapidly changing world. Evidence has been emerging for how fundamental drivers of density dependence in mobile species are related to two differing types of environmental variation-temporal variation in climate, and spatiotemporal variation in food resources. However, to date, tests of these hypotheses have been largely restricted to mid-trophic species in terrestrial environments and thus their general applicability remains unknown. We tested if these same processes can be identified in marine upper trophic level species. We assembled a multi-decadal data set on population abundance of 10 species of colonial seabirds comprising a large component of the UK breeding seabird biomass, and covering diverse phylogenies, life histories and foraging behaviours. We tested for evidence of density dependence in population growth rates using discrete time state-space population models fit to long time-series of observations of abundance at seabird breeding colonies. We then assessed if the strength of density dependence in population growth rates was exacerbated by temporal variation in climate (sea temperature and swell height), and attenuated by spatiotemporal variation in prey resources (productivity and tidal fronts). The majority of species showed patterns consistent with temporal variation in climate acting to strengthen density dependent feedbacks to population growth. However, fewer species showed evidence for a weakening of density dependence with increasing spatiotemporal variation in prey resources. Our findings extend this emerging theory for how different sources of environmental variation may shape the dynamics and regulation of animal populations, demonstrating its role in upper trophic marine species. We show that environmental variation leaves a signal in long-term population dynamics of seabirds with potentially important consequences for their demography and trophic interactions.

Modelling the responses of partially migratory metapopulations to changing seasonal migration rates: From theory to data.

Among-individual and within-individual variation in expression of seasonal migration versus residence is widespread in nature and could substantially affect the dynamics of partially migratory metapopulations inhabiting seasonally and spatially structured environments. However, such variation has rarely been explicitly incorporated into metapopulation dynamic models for partially migratory systems. We, therefore, lack general frameworks that can identify how variable seasonal movements, and associated season- and location-specific vital rates, can control system persistence. We constructed a novel conceptual framework that captures full-annual-cycle dynamics and key dimensions of metapopulation structure for partially migratory species inhabiting seasonal environments. We conceptualize among-individual variation in seasonal migration as two variable vital rates: seasonal movement probability and associated movement survival probability. We conceptualize three levels of within-individual variation (i.e. plasticity), representing seasonal or annual variation in seasonal migration or lifelong fixed strategies. We formulate these concepts as a general matrix model, which is customizable for diverse life-histories and seasonal landscapes. To illustrate how variable seasonal migration can affect metapopulation growth rate, demographic structure and vital rate elasticities, we parameterize our general models for hypothetical short- and longer-lived species. Analyses illustrate that elasticities of seasonal movement probability and associated survival probability can sometimes equal or exceed those of vital rates typically understood to substantially influence metapopulation dynamics (i.e. seasonal survival probability or fecundity), that elasticities can vary non-linearly, and that metapopulation outcomes depend on the level of within-individual plasticity. We illustrate how our general framework can be applied to evaluate the consequences of variable and changing seasonal movement probability by parameterizing our models for a real partially migratory metapopulation of European shags Gulosus aristotelis assuming lifelong fixed strategies. Given observed conditions, metapopulation growth rate was most elastic to breeding season adult survival of the resident fraction in the dominant population. However, given doubled seasonal movement probability, variation in survival during movement would become the primary driver of metapopulation dynamics. Our general conceptual and matrix model frameworks, and illustrative analyses, thereby highlight complex ways in which structured variation in seasonal migration can influence dynamics of partially migratory metapopulations, and pave the way for diverse future theoretical and empirical advances.

Variation and correlation in the timing of breeding of North Atlantic seabirds across multiple scales.

Timing of breeding, an important driver of fitness in many populations, is widely studied in the context of global change, yet despite considerable efforts to identify environmental drivers of seabird nesting phenology, for most populations we lack evidence of strong drivers. Here we adopt an alternative approach, examining the degree to which different populations positively covary in their annual phenology to infer whether phenological responses to environmental drivers are likely to be (a) shared across species at a range of spatial scales, (b) shared across populations of a species or (c) idiosyncratic to populations. We combined 51 long-term datasets on breeding phenology spanning 50 years from nine seabird species across 29 North Atlantic sites and examined the extent to which different populations share early versus late breeding seasons depending on a hierarchy of spatial scales comprising breeding site, small-scale region, large-scale region and the whole North Atlantic. In about a third of cases, we found laying dates of populations of different species sharing the same breeding site or small-scale breeding region were positively correlated, which is consistent with the hypothesis that they share phenological responses to the same environmental conditions. In comparison, we found no evidence for positive phenological covariation among populations across species aggregated at larger spatial scales. In general, we found little evidence for positive phenological covariation between populations of a single species, and in many instances the inter-year variation specific to a population was substantial, consistent with each population responding idiosyncratically to local environmental conditions. Black-legged kittiwake Rissa tridactyla was the exception, with populations exhibiting positive covariation in laying dates that decayed with the distance between breeding sites, suggesting that populations may be responding to a similar driver. Our approach sheds light on the potential factors that may drive phenology in our study species, thus furthering our understanding of the scales at which different seabirds interact with interannual variation in their environment. We also identify additional systems and phenological questions to which our inferential approach could be applied.

Episodes of opposing survival and reproductive selection cause strong fluctuating selection on seasonal migration versus residence.

Quantifying temporal variation in sex-specific selection on key ecologically relevant traits, and quantifying how such variation arises through synergistic or opposing components of survival and reproductive selection, is central to understanding eco-evolutionary dynamics, but rarely achieved. Seasonal migration versus residence is one key trait that directly shapes spatio-seasonal population dynamics in spatially and temporally varying environments, but temporal dynamics of sex-specific selection have not been fully quantified. We fitted multi-event capture-recapture models to year-round ring resightings and breeding success data from partially migratory European shags (Phalacrocorax aristotelis) to quantify temporal variation in annual sex-specific selection on seasonal migration versus residence arising through adult survival, reproduction and the combination of both (i.e. annual fitness). We demonstrate episodes of strong and strongly fluctuating selection through annual fitness that were broadly synchronized across females and males. These overall fluctuations arose because strong reproductive selection against migration in several years contrasted with strong survival selection against residence in years with extreme climatic events. These results indicate how substantial phenotypic and genetic variation in migration versus residence could be maintained, and highlight that biologically important fluctuations in selection may not be detected unless both survival selection and reproductive selection are appropriately quantified and combined.

Meeting Paris agreement objectives will temper seabird winter distribution shifts in the North Atlantic Ocean.

We explored the implications of reaching the Paris Agreement Objective of limiting global warming to <2°C for the future winter distribution of the North Atlantic seabird community. We predicted and quantified current and future winter habitats of five North Atlantic Ocean seabird species (Alle alle, Fratercula arctica, Uria aalge, Uria lomvia and Rissa tridactyla) using tracking data for ~1500 individuals through resource selection functions based on mechanistic modeling of seabird energy requirements, and a dynamic bioclimate envelope model of seabird prey. Future winter distributions were predicted to shift with climate change, especially when global warming exceed 2°C under a "no mitigation" scenario, modifying seabird wintering hotspots in the North Atlantic Ocean. Our findings suggest that meeting Paris agreement objectives will limit changes in seabird selected habitat location and size in the North Atlantic Ocean during the 21st century. We thereby provide key information for the design of adaptive marine-protected areas in a changing ocean.

Environmental heterogeneity promotes individual specialisation in habitat selection in a widely distributed seabird.

Individual specialisations in behaviour are predicted to arise where divergence benefits fitness. Such specialisations are more likely in heterogeneous environments where there is both greater ecological opportunity and competition-driven frequency dependent selection. Such an effect could explain observed differences in rates of individual specialisation in habitat selection, as it offers individuals an opportunity to select for habitat types that maximise resource gain while minimising competition; however, this mechanism has not been tested before. Here, we use habitat selection functions to quantify individual specialisations while foraging by black-legged kittiwakes Rissa tridactyla, a marine top predator, at 15 colonies around the United Kingdom and Ireland, along a gradient of environmental heterogeneity. We find support for the hypothesis that individual specialisations in habitat selection while foraging are more prevalent in heterogeneous environments. This trend was significant across multiple dynamic habitat variables that change over short time-scales and did not arise through site fidelity, which highlights the importance of environmental processes in facilitating behavioural adaptation by predators. Individual differences may drive evolutionary processes, and therefore these results suggest that there is broad scope for the degree of environmental heterogeneity to determine current and future population, species and community dynamics.

No evidence for fitness signatures consistent with increasing trophic mismatch over 30 years in a population of European shag Phalacrocorax aristotelis.

As temperatures rise, timing of reproduction is changing at different rates across trophic levels, potentially resulting in asynchrony between consumers and their resources. The match-mismatch hypothesis (MMH) suggests that trophic asynchrony will have negative impacts on average productivity of consumers. It is also thought to lead to selection on timing of breeding, as the most asynchronous individuals will show the greatest reductions in fitness. Using a 30-year individual-level dataset of breeding phenology and success from a population of European shags on the Isle of May, Scotland, we tested a series of predictions consistent with the hypothesis that fitness impacts of trophic asynchrony are increasing. These predictions quantified changes in average annual breeding success and strength of selection on timing of breeding, over time and in relation to rising sea surface temperature (SST) and diet composition. Annual average (population) breeding success was negatively correlated with average lay date yet showed no trend over time, or in relation to increasing SST or the proportion of principal prey in the diet, as would be expected if trophic mismatch was increasing. At the individual level, we found evidence for stabilising selection and directional selection for earlier breeding, although the earliest birds were not the most productive. However, selection for earlier laying did not strengthen over time, or in relation to SST or slope of the seasonal shift in diet from principal to secondary prey. We found that the optimum lay date advanced by almost 4 weeks during the study, and that the population mean lay date tracked this shift. Our results indicate that average performance correlates with absolute timing of breeding of the population, and there is selection for earlier laying at the individual level. However, we found no fitness signatures of a change in the impact of climate-induced trophic mismatch, and evidence that shags are tracking long-term shifts in optimum timing. This suggests that if asynchrony is present in this system, breeding success is not impacted. Our approach highlights the advantages of examining variation at both population and individual levels when assessing evidence for fitness impacts of trophic asynchrony.

Strong survival selection on seasonal migration versus residence induced by extreme climatic events.

Elucidating the full eco-evolutionary consequences of climate change requires quantifying the impact of extreme climatic events (ECEs) on selective landscapes of key phenotypic traits that mediate responses to changing environments. Episodes of strong ECE-induced selection could directly alter population composition, and potentially drive micro-evolution. However, to date, few studies have quantified ECE-induced selection on key traits, meaning that immediate and longer-term eco-evolutionary implications cannot yet be considered. One widely expressed trait that allows individuals to respond to changing seasonal environments, and directly shapes spatio-seasonal population dynamics, is seasonal migration versus residence. Many populations show considerable among-individual phenotypic variation, resulting in 'partial migration'. However, variation in the magnitude of direct survival selection on migration versus residence has not been rigorously quantified, and empirical evidence of whether seasonal ECEs induce, intensify, weaken or reverse such selection is lacking. We designed full annual cycle multi-state capture-recapture models that allow estimation of seasonal survival probabilities of migrants and residents from spatio-temporally heterogeneous individual resightings. We fitted these models to 9 years of geographically extensive year-round resighting data from partially migratory European shags Phalacrocorax aristotelis. We thereby quantified seasonal and annual survival selection on migration versus residence across benign and historically extreme non-breeding season (winter) conditions, and tested whether selection differed between females and males. We show that two of four observed ECEs, defined as severe winter storms causing overall low survival, were associated with very strong seasonal survival selection against residence. These episodes dwarfed the weak selection or neutrality evident otherwise, and hence caused selection through overall annual survival. The ECE that caused highest overall mortality and strongest selection also caused sex-biased mortality, but there was little overall evidence of sex-biased selection on migration versus residence. Our results imply that seasonal ECEs and associated mortality can substantially shape the landscape of survival selection on migration versus residence. Such ECE-induced phenotypic selection will directly alter migrant and resident frequencies, and thereby alter immediate spatio-seasonal population dynamics. Given underlying additive genetic variation, such ECEs could potentially cause micro-evolutionary changes in seasonal migration, and thereby cause complex eco-evolutionary population responses to changing seasonal environments.

Among-individual and within-individual variation in seasonal migration covaries with subsequent reproductive success in a partially migratory bird.

Within-individual and among-individual variation in expression of key environmentally sensitive traits, and associated variation in fitness components occurring within and between years, determine the extents of phenotypic plasticity and selection and shape population responses to changing environments. Reversible seasonal migration is one key trait that directly mediates spatial escape from seasonally deteriorating environments, causing spatio-seasonal population dynamics. Yet, within-individual and among-individual variation in seasonal migration versus residence, and dynamic associations with subsequent reproductive success, have not been fully quantified. We used novel capture-mark-recapture mixture models to assign individual European shags (Phalacrocorax aristotelis) to 'resident', 'early migrant', or 'late migrant' strategies in two consecutive years, using year-round local resightings. We demonstrate substantial among-individual variation in strategy within years, and directional within-individual change between years. Furthermore, subsequent reproductive success varied substantially among strategies, and relationships differed between years; residents and late migrants had highest success in the 2 years, respectively, matching the years in which these strategies were most frequently expressed. These results imply that migratory strategies can experience fluctuating reproductive selection, and that flexible expression of migration can be partially aligned with reproductive outcomes. Plastic seasonal migration could then potentially contribute to adaptive population responses to currently changing forms of environmental seasonality.

Interactions between Environmental Contaminants and Gastrointestinal Parasites: Novel Insights from an Integrative Approach in a Marine Predator.

Environmental contaminants and parasites are ubiquitous stressors that can affect animal physiology and derive from similar dietary sources (co-exposure). To unravel their interactions in wildlife, it is thus essential to quantify their concurring drivers. Here, the relationship between blood contaminant residues (11 trace elements and 17 perfluoroalkyl substances) and nonlethally quantified gastrointestinal parasite loads was tested while accounting for intrinsic (sex, age, and mass) and extrinsic factors (trophic ecology inferred from stable isotope analyses and biologging) in European shags Phalacrocorax aristotelis. Shags had high mercury (range 0.65-3.21 µg g-1 wet weight, ww) and extremely high perfluorooctanoic acid (PFOA) and perfluorononanoic acid (PFNA) residues (3.46-53 and 4.48-44 ng g-1 ww, respectively). Males had higher concentrations of arsenic, mercury, PFOA, and PFNA than females, while the opposite was true for selenium, perfluorododecanoic acid (PFDoA), and perfluooctane sulfonic acid (PFOS). Individual parasite loads (Contracaecum rudolphii) were higher in males than in females. Females targeted pelagic-feeding prey, while males relied on both pelagic- and benthic-feeding organisms. Parasite loads were not related to trophic ecology in either sex, suggesting no substantial dietary co-exposure with contaminants. In females, parasite loads increased strongly with decreasing selenium:mercury molar ratios. Females may be more susceptible to the interactive effects of contaminants and parasites on physiology, with potential fitness consequences.

Environmental heterogeneity decreases reproductive success via effects on foraging behaviour.

Environmental heterogeneity shapes the uneven distribution of resources available to foragers, and is ubiquitous in nature. Optimal foraging theory predicts that an animal's ability to exploit resource patches is key to foraging success. However, the potential fitness costs and benefits of foraging in a heterogeneous environment are difficult to measure empirically. Heterogeneity may provide higher-quality foraging opportunities, or alternatively could increase the cost of resource acquisition because of reduced patch density or increased competition. Here, we study the influence of physical environmental heterogeneity on behaviour and reproductive success of black-legged kittiwakes, Rissa tridactyla. From GPS tracking data at 15 colonies throughout their British and Irish range, we found that environments that were physically more heterogeneous were associated with longer trip duration, more time spent foraging while away from the colony, increased overlap of foraging areas between individuals and lower breeding success. These results suggest that there is greater competition between individuals for finite resources in more heterogeneous environments, which comes at a cost to reproduction. Resource hotspots are often considered beneficial, as individuals can learn to exploit them if sufficiently predictable. However, we demonstrate here that such fitness gains can be countered by greater competition in more heterogeneous environments.

Earlier colony arrival but no trend in hatching timing in two congeneric seabirds (Uria spp.) across the North Atlantic.

A global analysis recently showed that seabird breeding phenology (as the timing of egg-laying and hatching) does not, on average, respond to temperature changes or advance with time (Keogan et al. 2018 Nat. Clim. Change8, 313-318). This group, the most threatened of all birds, is therefore prone to spatio-temporal mismatches with their food resources. Yet, other aspects of the breeding phenology may also have a marked influence on breeding success, such as the arrival date of adults at the breeding site following winter migration. Here, we used a large tracking dataset of two congeneric seabirds breeding in 14 colonies across 18° latitudes, to show that arrival date at the colony was highly variable between colonies and species (ranging 80 days) and advanced 1.4 days/year while timing of egg-laying remained unchanged, resulting in an increasing pre-laying duration between 2009 and 2018. Thus, we demonstrate that potentially not all components of seabird breeding phenology are insensitive to changing environmental conditions.

The energetic cost of parasitism in a wild population.

Parasites have profound fitness effects on their hosts, yet these are often sub-lethal, making them difficult to understand and quantify. A principal sub-lethal mechanism that reduces fitness is parasite-induced increase in energetic costs of specific behaviours, potentially resulting in changes to time and energy budgets. However, quantifying the influence of parasites on these costs has not been undertaken in free-living animals. We used accelerometers to estimate energy expenditure on flying, diving and resting, in relation to a natural gradient of endo-parasite loads in a wild population of European shags Phalacrocorax aristotelis We found that flight costs were 10% higher in adult females with higher parasite loads and these individuals spent 44% less time flying than females with lower parasite loads. There was no evidence for an effect of parasite load on daily energy expenditure, suggesting the existence of an energy ceiling, with the increase in cost of flight compensated for by a reduction in flight duration. These behaviour specific costs of parasitism will have knock-on effects on reproductive success, if constraints on foraging behaviour detrimentally affect provisioning of young. The findings emphasize the importance of natural parasite loads in shaping the ecology and life-history of their hosts, which can have significant population level consequences.

The role of parasitism in the energy management of a free-ranging bird.

Parasites often prompt sub-lethal costs to their hosts by eliciting immune responses. These costs can be hard to quantify but are crucial to our understanding of the host's ecology. Energy is a fundamental currency to quantify these costs, as energetic trade-offs often exist between key fitness-related processes. Daily energy expenditure (DEE) comprises of resting metabolic rate (RMR) and energy available for activity, which are linked via the energy management strategy of an organism. Parasitism may play a role in the balance between self-maintenance and activity, as immune costs can be expressed in elevated RMR. Therefore, understanding energy use in the presence of parasitism enables mechanistic elucidation of potential parasite costs. Using a gradient of natural parasite load and proxies for RMR and DEE in a wild population of breeding European shags (Phalacrocorax aristotelis), we tested the effect of parasitism on maintenance costs as well as the relationship between proxies for RMR and DEE. We found a positive relationship between parasite load and our RMR proxy in females but not males, and no relationship between proxies for RMR and DEE. This provides evidence for increased maintenance costs in individuals with higher parasite loads and suggests the use of an allocation energy management strategy, whereby an increase to RMR creates restrictions on energy allocation to other activities. This is likely to have fitness consequences as energy allocated to immunity is traded off against reproduction. Our findings demonstrate that understanding energy management strategies alongside fitness drivers is central to understanding the mechanisms by which these drivers influence individual fitness.