# Can attraction to and competition for high-quality habitats shape breeding propensity?

Paul Acker, Michael Schaub, Aurélien Besnard, Jean-Yves Monnat,
Emmanuelle Cam

► To cite this version:

Paul Acker, Michael Schaub, Aurélien Besnard, Jean-Yves Monnat, Emmanuelle Cam. Can attraction to and competition for high-quality habitats shape breeding propensity?. Journal of Animal Ecology, 2022, 91 (5), pp.933-945. 10.1111/1365-2656.13676. hal-04071559

HAL Id: hal-04071559

https://hal.univ-brest.fr/hal-04071559v1

Submitted on 8 Jul 2024

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

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

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

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

# Can attraction to and competition for high-quality habitats shape breeding propensity?

Paul Acker<fnref n="1" />,²,³ ID | Michael Schaub⁴ ID | Aurélien Besnard³ ID | Jean-Yves Monnat⁵ |
Emmanuelle Cam²,⁶ ID

**Correspondence**

Paul Acker

Email: paul.acker@ntnu.no

**Funding information**

LabEx TULIP, Grant/Award Number: ANR-10-LABX-41; ANR-11-IDEX-0002-02

Handling Editor: Thierry Boulinier

## Abstract

1. In many animal species, sexually mature individuals may skip breeding opportunities despite a likely negative impact on fitness. In spatio-temporally heterogeneous environments, habitat selection theory predicts that individuals select habitats where fitness prospects are maximized. Individuals are attracted to high-quality habitat patches where they compete for high-quality breeding sites. Since failures in contests to secure a site may prevent individuals from breeding, we hypothesized that attraction to and competition for high-quality habitats could shape breeding propensity.

2. Under this hypothesis, we predicted the two following associations between breeding propensity and two key population features. (1) When mean habitat quality in the population increases in multiple patches such that availability of high-quality sites increases across the population, the resulting decrease in competition should positively affect breeding propensity. (2) When the number of individuals increases in the population, the resulting increase in competitors should negatively affect breeding propensity (negative density dependence).

3. Using long-term data from Kittiwakes *Rissa tridactyla*, we checked the prerequisite of prediction (1), that availability of high-quality sites is positively associated with current mean habitat quality in the population (represented by breeding success). We then applied integrated population modelling to quantify annual fluctuations in population mean breeding success, breeding propensity and number of individuals by breeding status (pre-breeders, breeders, skippers and immigrants), and tested our predictions.

4. Our results showed that breeding propensity acts as an important driver of population growth. As expected, breeding propensity was positively associated with preceding mean habitat quality in the population, and negatively with the number of competitors. These relationships varied depending on breeding status, which likely reflects status dependence in competitive ability.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2022 The Authors. *Journal of Animal Ecology* published by John Wiley & Sons Ltd on behalf of British Ecological Society.

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

<footnote n="1">¹Centre for Biodiversity Dynamics, Institutt for Biologi, NTNU, Trondheim, Norway; ²Laboratoire EDB (UMR 5174), Université Paul Sabatier - CNRS - IRD, Toulouse, France; ³CEFE, Univ Montpellier, CNRS, EPHE-PSL University, IRD, Montpellier, France; ⁴Swiss Ornithological Institute, Sempach, Switzerland; ⁵Pennarun d'An Traon, Goulien, France and ⁶Univ Brest, CNRS, IRD, Ifremer, LEMAR, Plouzané, France</footnote>
5. These findings highlight the importance of competition for high-quality breeding sites in shaping breeding propensity. Thereby, we draw attention towards alternative and complementary explanations to more standard considerations regarding the energetic cost of reproduction, and point to possible side effects of habitat selection behaviours on individual life histories and population dynamics.

### KEYWORDS

breeding decision, breeding habitat selection, density dependence, immigration, integrated population model, intraspecific competition, recruitment, reproductive skipping

# 1 | INTRODUCTION

At each reproductive occasion, sexually mature individuals experience various external and internal constraints that may alter their breeding propensity. This can have major impacts on individual fitness, population growth and demographic structure (Lee et al., 2016; Stearns, 1992). Non-breeders lose their current reproductive value and are at risk of dying before the next breeding opportunity, yet they often represent a non-negligible part of the population in long-lived species. These can be individuals that have not recruited yet ('pre-breeders') or have already bred previously ('skippers'), as found in a wide range of taxa, spanning fish (e.g. Rideout & Tomkiewicz, 2011), reptiles (e.g. Shine & Brown, 2008), amphibians (e.g. Cayuela

et al., 2014), birds (e.g. Bruinzeel, 2007) and mammals (e.g. Desprez et al., 2018). Identifying the factors that lead to non-breeding is thus critical to understand key eco-evolutionary processes underlying population dynamics.

Individuals may not breed simply because they do not fulfil essential requirements, that is they did not accumulate sufficient energy reserves (Giudici et al., 2010; Meijer & Drent, 1999) or failed to acquire a mate or a *breeding site* (Box 1; Danchin & Cam, 2002; Bruinzeel, 2007). Meeting these breeding requirements is costly, and subsequent breeding activities are costly too. According to life-history theory, the proportion of finite resources allocated to current reproduction is traded off against the proportion allocated to survival and/or future reproduction (Stearns, 1992). Consequently,

## BOX 1 Glossary

**Breeding habitat selection:** The choice made by an individual to occupy a given breeding habitat. This choice typically involves the use of cues allowing an organism to assess habitat quality (e.g. conspecific breeding success). Such a choice may not be attained, for example if competitive inferiority prevents the individual from acquiring a breeding site.

**Breeding patch:** The space containing a contiguous set of breeding sites. Breeding patches can be considered at various spatial scales: for example, in Kittiwakes, a patch can be a subpart of a cliff wall, an entire cliff wall, a cove consisting of several cliff walls or a colony consisting of contiguous coves.

**Breeding propensity:** The tendency of individuals to breed at a given occasion. In a population, breeding propensity is typically measured using breeding probability, independent of the patch where individuals will breed (since individuals may disperse between patches). For immigrants, it is represented by the proportion of individuals in the local population (the immigration rate), because the source population is unknown.

**Breeding site:** Space that is used by an individual (or a pair) to reproduce (e.g. in the Kittiwake, where a pair builds a nest to lay eggs and rear chicks). It is also termed 'breeding territory' in species where individuals defend a delimited location against intruders.

**Habitat quality:** The expected fitness prospects offered to an individual by a given habitat (i.e. a breeding site or a breeding patch or the full set of patches in the population), and that varies according to spatio-temporally heterogeneous factors (e.g. climate, vegetation, predation, food availability and parasitism). In temporally autocorrelated environments, it is best approximated by the preceding mean fitness of individuals in the habitat.

**Partial correlation:** Value of the correlation between two variables when other covariates are held constant in the sample (i.e. controlling for the confounding effect of the other covariates). If two processes corresponding to non-mutually exclusive hypotheses are responsible for a relationship (e.g. the number of competitors and population habitat quality both influence breeding propensity), partial correlations allow a process to be detected while the other is also operating. This approach is not designed to discount the hypothesis corresponding to the process that is controlled for.
when individuals face high costs of current reproduction, non-breeding could allow maximizing long-term fitness prospects and be selectively advantageous (Desprez et al., 2018; Erikstad et al., 1998). This can explain why unfavourable environmental conditions experienced by a population (e.g. lower overall food availability) are associated with lower subsequent *breeding propensity* (Cayuela et al., 2014; Hoy et al., 2016; Rideout & Tomkiewicz, 2011; Shine & Brown, 2008).

Yet environmental conditions are typically varying not only across time, but also across space. In species moving actively, mechanisms of *breeding habitat selection* (Box 1) that allow individuals to assess *habitat quality* (Box 1) and occupy the best possible habitats are expected to have evolved, on condition that environments are temporally autocorrelated (Doligez et al., 2003; Fretwell & Lucas, 1969; Johnson, 2007). Individuals looking for a breeding site should be attracted to *breeding patches* (Box 1) that they perceive as high-quality ones, and therefore by high-quality breeding sites that are likely already occupied or targeted by others, generating competition (Acker et al., 2017; Fretwell & Lucas, 1969; Lima & Zollner, 1996; Pulliam & Danielson, 1991). Consequently, non-breeding could result from failure in the contest for dominance on a high-quality site if individuals do not have enough time and energy to secure another —potentially lower quality— site or mate while competing for a high-quality one (Bruinzeel, 2007; Danchin & Cam, 2002; Kokko et al., 2004). Non-breeding could also result from queueing behaviour. Indeed, waiting for vacancy of a high-quality site may offer better long-term fitness prospects than breeding on a lower quality site (Ens et al., 1995; Zack & Stutchbury, 1992). Territorial competition for high-quality breeding sites can therefore be hypothesized to contribute to shaping breeding propensity and hence contribute to population dynamics (Kokko & Sutherland, 1998; Newton, 1992), but empirical evidence is lacking.

Under this hypothesis, breeding propensity would vary over time with the intensity of competition for high-quality sites in the population —which depends on the availability of disputed resources (high-quality sites) and the number of competitors (individuals already occupying a site or looking for a site). Two predictions can be made regarding how breeding propensity is associated with two key population factors linked to competition for breeding sites: mean breeding habitat quality (hereafter 'population habitat quality') and number of conspecific individuals. (1) If increased population habitat quality occurs through increased habitat quality across multiple patches (i.e. decreased spatial heterogeneity and decreased variation in attractivity among patches), this would imply higher availability of high-quality sites (whether they are occupied or not). Number of individuals being equal, the resulting competition decrease should be associated with increased breeding propensity. (2) If the number of individuals in the population increases, this would imply a higher number of competitors. Spatial heterogeneity of habitat quality being equal, the resulting competition increase should be associated with decreased breeding propensity.

Under prediction (1), we expect a positive relationship between mean reproductive success of breeders in the population (hereafter 'population breeding success', representing population habitat quality) and subsequent breeding propensity. Such a relationship would be detected while controlling for the confounding effect of the number of individuals in the population. This relationship could also result from spatially homogeneous temporal variation in environmental conditions —and thus in habitat quality— affecting the energetic cost of reproduction (e.g. food availability or weather conditions; Cayuela et al., 2018; Hoy et al., 2016; Nur & Sydeman, 1999). But this would contrast with situations where temporal variation in population habitat quality is spatially heterogeneous and where the prerequisite to prediction (1) is fulfilled: a tight negative relationship between population habitat quality and the degree of spatial heterogeneity of habitat quality. Under prediction (2), we expect a negative relationship between the numbers of individuals in the population and subsequent breeding propensity, that is negative density dependence in breeding propensity. Here again, such a relationship could also result from competition for food independent of competition for breeding sites. But if so, one would also expect competition for food to underlie a concomitant negative correlation between the number of individuals and population breeding success (e.g. Layton-Matthews et al., 2019). By controlling for the confounding effect of population breeding success when testing for a negative relationship between the number of individuals and breeding propensity, one will thus detect the distinctive effect of competition for high-quality breeding sites.

We tested our predictions in a population of black-legged kitting wakes *Rissa tridactyla*, using 28 years of monitoring data of all active nests (~1000 each year) and capture-resighting histories of >12,000 individuals. In this system, habitat selection for breeding at time *t* involves attraction to and intense competition for sites located in high-quality patches, and individuals identify such patches via the reproductive success of conspecifics at the end of the previous breeding season *t* − 1 (Acker et al., 2017; Cadiou et al., 1994; Danchin et al., 1998; Appendix S1.2). We first quantified the relationship between the degree of spatial heterogeneity and the mean breeding success in the population to check the prerequisite to prediction (1) that population habitat quality is positively associated with the availability of high-quality sites in the population. We then designed an integrated population model ('IPM') to jointly quantify fluctuations in population-wide numbers of pre-breeders, breeders, skippers and immigrants ('breeding status'), their breeding propensity and population breeding success. We used IPM estimates to quantify the relationships between breeding propensity and numbers of competitors or population habitat quality, which allowed us to test predictions (1) and (2). More precisely, we assessed whether breeding propensity at *t* was positively correlated with population breeding success at *t* − 1 and negatively correlated with the numbers of breeders or non-breeders at *t* − 1 (in each case, controlling for the other covariates using *partial correlations*, Box 1).
## 2 | MATERIALS AND METHODS

### 2.1 | Population monitoring

The data were collected in the Cap Sizun Kittiwake population (48°03′N, 4°39′W; Brittany, France), where thousands of chicks have been individually marked with colour rings since 1979 (Appendix S1.1). Our analyses are based on data from 1985 to 2012. Monitoring was carried out throughout each breeding season by visiting all colonies weekly from first arrivals to the fledging period (January–June), and then daily until bird departures (July–August). During visits, the content of every nest site was recorded to determine breeding success, and the location and behaviour of ringed birds were recorded to determine breeding status (Cam et al., 1998). All fieldwork was licensed by the Centre de Recherches sur la Biologie des populations d'Oiseaux ('CRBPO', Muséum National d'Histoire Naturelle, Paris, France), and carried out in accordance with standard animal care protocols approved by the CRBPO

Resighting probability is virtually equal to one once an individual is recruited to the breeding population (age 3 at the earliest; Cam et al., 1998; see also Section 3). Whether they breed or not, the intensive resighting effort allows individuals that attend the breeding cliffs to be detected. Known-age individuals are considered 'pre-breeders' before their first breeding attempt in the population (recruitment), 'breeders' when they completed nest building in the current year (Cullen, 1957) or 'skippers' when they bred in the past but did not complete nest building in the current year. Pre-breeders not always show up at the colonies in the breeding season, and those that attend the colonies may enter territorial contests for nest sites mainly at the end of the season. Skippers attend the breeding colonies, and their behaviour ranges from aterritorial floating to consistent territory holding throughout the season, including territorial contests for occupied and non-occupied sites.

The breeding success of each nest was assessed using the number of chicks that reached at least fledging age (35 days or more). Breeding population was counted using the annual number of breeders, derived as twice the number of nests. Pairs very rarely build two nests successively; for marked individuals, successive nests were assigned to a unique pair.

### 2.2 | Spatio-temporal variation in habitat quality

The breeding habitat consists of multiple patches, which can be considered at various spatial scales: geographically distinct colonies (2–5, distant from each other by 0.5–12 km) composed of contiguous coves (5–18) including cliff walls (20–44; separated from each other by rocky ridges or coastal segments without nesting birds), themselves divisible into smaller heterogeneous patches (Acker et al., 2017; Bled et al., 2011; Danchin et al., 1998; Appendix S1.1,2). There is substantial among- and within-patch heterogeneity in habitat quality, and this spatial heterogeneity is dynamic across years (Acker et al., 2017; Danchin et al., 1998; Appendix S1.1,2). At every

spatial scale, patches are connected through large natal and breeding dispersal flows, which are typically directed from current low-quality patches to high-quality ones (Appendix S1.2).

During the study period, spatio-temporal variation in breeding success is believed to have been mostly caused by predation on eggs by corvids and predation on young chicks by herring gulls in one colony, which eventually led to massive dispersal within the population and desertion of entire colonies (Acker et al., 2017; Cam et al., 2004; Danchin et al., 1998; Appendix S1.1,2). Ticks (*Ixodes uriae*) have also been suggested as a potential driver of variation in breeding success (Danchin et al., 1998). Food availability is unlikely to have caused the large spatial heterogeneity in breeding habitat quality that we observed, since Kittiwakes feed on non-defendable resources of which the availability varies at much larger spatial scales than within-population foraging destinations (Christensen-Dalsgaard et al., 2018; Oro & Furness, 2002; Suryan et al., 2002).

Prediction (1) relies on the prerequisite that increased population habitat quality arises from decreased spatial heterogeneity of habitat quality among patches, implying increased availability of high-quality sites. To check this prerequisite, we evaluated spatial heterogeneity of habitat quality every year by inspecting the distribution of mean breeding success among patches (weighted by the number of sites occupied by breeders in each patch) and measuring its dispersion via the Gini coefficient (Appendix S1.3). We show our results at the cliff scale (we found similar patterns at the cove and colony scale, but no smaller scale was investigated; Appendix S1.3). Specifically, low Gini coefficients (low heterogeneity) corresponded to distributions packed around the mean, and high Gini coefficients (high heterogeneity) corresponded to years when a large proportion of patches had very low breeding success and only a few had high breeding success (Appendix S1.3). We found a strong negative correlation between the degree of spatial heterogeneity of habitat quality (measured by the Gini coefficient) and population breeding success (Figure 1; Pearson's *r* = -0.79). Such a pattern demonstrates that our prerequisite is fulfilled.

### 2.3 | Modelling population dynamics

To test the population-level predictions from our hypothesis linking competition for high-quality habitats to breeding propensity, we required robust quantification of numbers of individuals, breeding propensity in each breeding status and population breeding success. We developed an IPM (Schaub & Abadi, 2011; Schaub & Kéry, 2021) to model population dynamics from the joint analyses of population counts, individual resightings and breeding success observations. Such a model allows the estimation of key demographic parameters while fully propagating uncertainty across the different types of observations. Our IPM notably allowed estimating the numbers of immigrants, unmarked skippers and pre-breeders that cannot be directly counted in the field or directly estimated from a single dataset.

The core of the IPM is a matrix population model (Caswell, 2001) depicting changes in the number of individuals in each state in year
t as a function of demographic rates and the number of individuals in each state in year t - 1. We designed the life cycle (Figure 2) and population matrix (pre-breeding census; Appendix S2.1) using prior knowledge of the population (Cam et al., 1998, 2002; Link et al., 2002). We defined nine life-history states: yearlings, pre-breeders of age 2, 3, 4, 5 and 6, first-time breeders (locals and immigrants), experienced breeders and skippers (Figure 2). Demographic rates were modelled as time dependent. To account for demographic stochasticity, the number of individuals in each state was modelled using a Poisson or binomial distribution (Appendix S2.1).

We considered two age classes for survival probability ($\phi_0$ at age 0-1 and $\phi_2$ at age $\ge$ 2; Figure 2; Link et al., 2002). Given the

<figure><img src="image_4.png" /><figcaption>FIGURE 1 Relationship between the degree of spatial heterogeneity of habitat quality (measured by the Gini coefficient, at the cliff scale) and population breeding success (mean number of fledglings per nest). Grey background: 95% confidence interval of regression line</figcaption></figure>

complexity of the IPM, we made simplifying assumptions regarding heterogeneity in survival from age 2 to achieve reasonable computing times. We assumed equal survival of immigrants and locals, which is necessary because immigrants are not individually monitored. If this assumption does not hold, the estimated number of immigrants could be negatively (or positively) biased when immigrants have lower (or higher) survival than locals. We have no a priori hypothesis concerning this point, but there is no indication that major differences are likely (Appendix S2.2). Even if our estimates were systematically biased, temporal variation in the number of immigrants relative to the mean should still be correctly inferred, and derived relationships between immigration and demographic features would be properly assessed.

We considered four age classes for recruitment rate ($\rho_i$, i.e. probability that a pre-breeder at age $i-1$ is a breeder at age $i$, $i \in [[3, 6]]$, conditional on survival from $i-1$ to $i$; Figure 2; Link et al., 2002). We assumed no recruitment before age 3 and after age 7 (recruitment rate at age 7 was fixed to 1), since no individual was ever recorded breeding at age 1, and we ignored the very few cases of recruitment at age 2 (~0.05% of individuals) and between ages 8 and 14 (~0.4%; Appendix S2.1). We considered status-dependent breeding rates ($\Psi_b$ and $\Psi_s$, i.e. probability of breeding at t, respectively for individuals that bred and skipped breeding at t-1; Figure 2; Cam et al., 1998).

We considered experience-dependent per capita breeding success rates (i.e. half the number of fledglings produced in the nest) of first-time breeders ($\pi_f$; Figure 2) and experienced breeders whatever their number of previous breeding attempts ($\pi_e$; Figure 2; Link et al., 2002). We assumed equal breeding success for immigrants and local first-time breeders. This assumption has been shown to have a negligible impact on estimates of immigrant numbers and demographic rates in the common tern (Sterna hirundo), which has a very similar life cycle to the Kittiwake (Szostek et al., 2014). This is expected because the population growth rate is not very sensitive to variation in fecundity parameters in long-lived species (Caswell, 2001).

<figure><img src="image_5.png" /><figcaption>FIGURE 2 Kittiwake life cycle graph underlying the integrated population model. Life-history states (black circles): yearlings (Y), pre-breeders of age $i$ ($\Pi_i$), first-time breeders (F), experienced breeders (E) and skippers (S). Black arrows: state transitions; subscripts: transition rates. Demographic parameters: survival at age 0 and 1 ($\phi_0$) and from age 2 ($\phi_2$), recruitment rate at age $i$ ($\rho_i$), breeding rate of former breeders ($\psi_b$) and former skippers ($\psi_s$), per capita breeding success of first-time breeders ($\pi_f$) and experienced breeders ($\pi_e$). In grey is the annual pulse of immigrants (I) into first-time breeders</figcaption></figure>
## 2.4 | IPM datasets

We analysed three datasets with the IPM: population count data, capture-recapture data and reproduction data. Population count data consisted of annual numbers of breeders, ranging 1316-2402 with large fluctuations (Figure 3a). Capture-recapture data consisted of capture-resighting histories of ringed birds (12,091 individuals, of which 642 were marked in 1979-1984, before the period modelled here) indicating age and breeding status (pre-breeder, breeder or skipper) at resighting. Reproduction data consisted of annual numbers of fledglings and corresponding numbers of nests belonging to pairs of either (a) first-time breeders (both mates ringed, 1962 breeding attempts), (b) experienced breeders (both mates ringed, 8785 breeding attempts) or (c) a first-time breeder mated with an experienced breeder (both ringed) or at least one unringed mate (25,366 breeding attempts).

## 2.5 | Likelihood

The IPM likelihood is the product of likelihoods of three models for the three datasets, assuming independence between datasets. In practice, this assumption of independence is not completely fulfilled, but simulations have shown that its violation has a very limited effect on parameter estimates (Schaub & Fletcher, 2015; Weegman et al., 2021). The likelihood given the population count data was formulated using a state-space model (Appendix S2.3). The state process was defined by the matrix population model in which fluctuations in class-specific population sizes are described. We assumed a log-normal distribution for the observation process with time-independent standard deviation. The likelihood given the individual capture-resighting histories was formulated using a state-space formulation of a multistate capture-recapture model (Appendix S2.4). We assumed different time-varying resighting rates for yearlings and pre-breeders, equal and temporally constant resighting rates for breeders and skippers, and no error in state assignment at resighting (Cam et al., 1998, 2002). The likelihood given the reproduction data was formulated using three Poisson regressions for fledgling numbers as a function of the number of nests and experience-dependent breeding success (Appendix S2.5). The three regressions were for pairs of (a) first-time breeders, (b) experienced breeders and (c) individuals of unknown or different categories of experience. For the latter, we ignored pair characteristics and assumed that their breeding success rate was the population breeding success, that is mean breeding success rate of inexperienced and experienced breeders weighted by their proportions among breeders.

## 2.6 | Inference and model assessment

To estimate model parameters, the joint likelihood was analysed in the Bayesian framework. We specified vague prior distributions with reasonable bounds for all parameters (Appendix S2.6). We used the

uniform distribution over [-5 to 1,000] as prior for the number of immigrants; the inclusion of negative values enables testing whether there is immigration at all (Schaub & Fletcher, 2015). We performed Markov Chain Monte Carlo simulation for posterior sampling using JAGS 3.4.0 (Plummer, 2003; model code and full details of sampling are in Appendices S2.7 and S3.1 respectively). While the capture-recapture data had already been analysed with similar multistate model structures (e.g. Cam et al., 1998, 2002; Link et al., 2002; yielding similar estimates, see Section 3), the additional analysis of population counts and reproduction data allowed estimating parameters that had not yet been estimated. We used posterior predictive checks (Gelman & Hill, 2006) to evaluate the fit of our IPM to the population count data and the reproduction data (Appendix S4). Overall, these checks indicated a good fit (Appendix S4).

## 2.7 | Derived quantities

We derived the posterior distribution of key quantities from model parameters, synthesizing compound biological effects of interest while accounting for their uncertainty (Appendix S5). Specifically, to characterize population dynamics with respect to breeding propensity, we derived the breeding population growth rate as the number of breeders in year *t* divided by the number of breeders in year *t* - 1 (Appendix S5.4). To characterize population composition, we derived the among-breeder proportions of former breeders (individuals that bred at *t* - 1), former skippers (individuals that skipped breeding at *t* - 1), local first-time breeders and immigrants. To synthesize the breeding propensity of all pre-breeders, we derived the age-independent 'integrative recruitment rate', that is the proportion of first-time breeders at *t* among the individuals of all age classes (3 to 6) alive and available for recruitment in the current year *t* (i.e. that have never bred before). Breeding propensity of immigrants was represented by the immigration rate, that is the proportion of immigrants among breeders in the current year (note that similar results were obtained using the absolute number of immigrants; Appendix S5.4). We also derived the number of non-breeders (i.e. pre-breeders plus skippers) present at the breeding colonies by correcting the number of non-breeders in the population by their resighting rate (Appendix S5.3).

## 2.8 | Correlates of demographic features and test of hypothesis regarding breeding propensity

Before testing our predictions, and to place our working hypothesis in the general demographic context, we assessed the contribution of breeding propensity at *t* of each breeding status at *t* - 1 to population dynamics. Specifically, we assessed the contributions of demographic rates to population fluctuations using estimates from the IPM. We derived posterior distributions of partial correlations between breeding population growth rate and survival rate, breeding rates of former breeders and former skippers, integrative local recruitment rate and
immigration rate (while controlling for each other's effects; Appendix S5.4; Szostek et al., 2014). Because information was insufficient in the first year to properly estimate the number of individuals that could not be counted in the field, we considered all parameters from the second year onwards to estimate partial correlations.

We then tested the population-level predictions of our hypothesis linking competition for high-quality habitats to breeding propensity using estimates from the IPM. We derived posterior distributions of partial correlations to test for relationships between breeding propensity (at *t*), and population habitat quality as well as abundance of competitors (at *t* − 1; Appendix S5.4). We assessed the relationship between breeding propensity at *t* in each breeding status (i.e. breeding rates of former breeders and former skippers, integrative recruitment rate and immigration rate) and population breeding success, number of breeders or number of present non-breeders at *t* − 1 (while controlling for each other's effects). We assessed the evidence for a partial correlation by computing the proportion of its posterior distribution that had the same sign as its posterior mean ('P'). Values of *P* close to 1 indicate strong evidence for a correlation with a given sign, while values close to 0.5 indicate no clear evidence (i.e. similar evidence for a negative or positive correlation).

## 3 | RESULTS

The estimates of breeding population size from the IPM closely matched the population count data (Figure 3a). Detailed posterior summaries of IPM parameters and derived quantities are given in Appendices S3 and S5 respectively. Hereafter, estimates are reported as the posterior mean with 95% credible interval ('95%CRI') in brackets.

### 3.1 | General demographic context

At the scale of the study period, breeding population size was stationary or nearly so (average growth rate: 1.001 [0.999, 1.004]; Appendix S5.1), despite large annual fluctuations (Figure 3a). The estimates of the demographic rates were consistent with those reported in previous studies not using an IPM (e.g. Cam et al., 1998, 2002; Link et al., 2002). Mean breeding success across years was 0.16 [0.14, 0.19] fledglings per capita for first-time breeders and 0.36 [0.33, 0.40] fledglings per capita for experienced breeders, resulting in population breeding success of 0.65 [0.64, 0.66] fledglings per nest, with large annual fluctuations indicating pronounced temporal variability in population habitat quality (Appendix S3.2; Figure 4). Mean local survival probability was 0.65 [0.59, 0.71] at ages 0 and 1, and 0.81 [0.78, 0.83] afterwards. Mean resighting probability, indicative of presence at the breeding grounds, was 0.05 [0.04, 0.07] for yearlings, 0.81 [0.78, 0.84] for older pre-breeders and 0.998 [0.997, 0.999] for recruited individuals. Mean recruitment rate at ages 3, 4, 5 and 6 was 0.13 [0.08, 0.18], 0.41 [0.34, 0.47], 0.53 [0.48, 0.59] and 0.67 [0.58, 0.76] respectively. The resulting mean integrative

recruitment rate (i.e. breeding propensity of pre-breeders) was 0.34 [0.33, 0.35]. Mean breeding rate was 0.90 [0.87, 0.92] for former breeders and 0.69 [0.62, 0.75] for former skippers. These breeding propensities were clearly lower than 1, and highly variable across years (Appendix S3.2; Figure 4), indicating the demographic importance of breeding propensity.

### 3.2 | Population composition

Among local individuals (i.e. locally born or already established in the population), there was a prominent proportion of breeders

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

**FIGURE 3** Population dynamics over 1985–2012. (a) Estimates of the numbers of pre-breeders (orange triangles), skippers (green diamonds) and breeders (blue dots), and nest count data (red circles). (b) Breakdown of the numbers of breeders into immigrants (grey downward triangles), local first-time breeders (orange upward triangles), former skippers (green diamonds) and former breeders (blue dots). Points: posterior means; shaded areas: 95%CRIs. In 1989, 2003, 2006–2008 and 2010–2012, 95%CRIs of the number of immigrants included negative values, suggesting that immigration may have been absent in these years
(1985–2012 mean: 62.1% [61.3, 62.8]), a moderate proportion of pre-breeders (30.2% [29.4, 31.0]) and a small proportion of skippers (7.7% [7.2, 8.3]), with large fluctuations (Figure 3a). There was a high turnover among breeders, with a mean of 30% of current breeders that had not bred in the population in the previous year (Figure 3b). Across years, the breeding population was composed on average of 7.6% [7.2, 8.0] local first-time breeders, 7.0% [6.5, 7.5] former skippers and 14.0% [12.9, 15.0] immigrants, versus 71.4% [70.4, 72.3] former breeders (Figure 3b). These results highlight how status-dependent breeding propensity shaped the highly dynamic compositions of the breeding and non-breeding segments of the population.

## 3.3 | Contribution of breeding propensity to population dynamics

The partial correlation with breeding population growth rate was 0.59 [0.29, 0.87] for immigration rate, 0.56 [0.36, 0.75] for breeding propensity of former breeders, 0.32 [0.05, 0.59] for breeding propensity of former skippers and 0.08 [-0.20, 0.36] for the integrative recruitment rate. For comparison, this partial correlation was 0.47 [0.27, 0.67] for local survival probability from age 2, that is the rate responsible for permanent disappearance of individuals from the breeding population. These values indicate that breeding propensity in all statuses except pre-breeders had non-negligible effects on temporal variation in breeding population growth, and these effects were especially high (and higher than the effect of survival) for immigrants and former breeders.

## 3.4 | Correlates of breeding propensity

We found positive associations between population habitat quality and breeding propensity in former pre-breeders and breeders, but not in former skippers and immigrants (Table 1; Figure 4). There was evidence of positive partial correlations between population breeding success in year *t* − 1 and both the breeding rate of former breeders in year *t* and the integrative recruitment rate (Table 1; Figure 4a,d). We also found negative associations between the numbers of competitors and breeding propensity (independently of population habitat quality) in all breeding statuses, with some status-dependent modulation in strength (Table 1; Figure 4). This is shown by negative partial correlations between the number of breeders at *t* − 1 and the breeding rates of former breeders and skippers, the integrative recruitment rate and the immigration rate at *t* — which was of lower magnitude for the latter two (Table 1, Figure 4b,c,e,f). There was also evidence of a negative partial correlation between the immigration rate at *t* and the number of non-breeders (pre-breeders plus skippers) present at *t* − 1 (Table 1; Figure 4g).

## 4 | DISCUSSION

Reproduction in animals is often contingent on acquisition or retention, or simply access to a breeding site (except when fertilization is external, or in non-territorial species). Accordingly, we hypothesized that habitat selection processes, including attraction to and competition for high-quality breeding sites, could influence whether individuals will breed or not. Given the finite availability of high-quality sites, we predicted that larger numbers of competitors generate lower breeding propensity due to competition for breeding sites, regardless of temporal variation in population habitat quality. We also predicted that higher habitat quality across a population, if realized through greater availability of high-quality sites across patches, relaxes competition in each high-quality patch, generating higher breeding propensity. Our integrated population model applied to long-termittiwake monitoring data allowed us to evidence relationships that match these predictions, in addition to demonstrating the critical role of breeding propensity for population growth and composition. Overall, in complement to standard energy-cost views on the achievement of reproductive careers, our study sheds light on the importance of competition for high-quality sites in shaping breeding propensity, individual life histories and population dynamics.

Previous studies have used the occurrence of delayed or skipped breeding and territorial behaviour in heterogeneous habitats to hypothesize that intense competition for high-quality breeding sites can drive non-breeding (Zack & Stutchbury, 1992), implying density dependence in breeding propensity (Kokko & Sutherland, 1998). This is corroborated by many studies that have shown that non-breeding is associated with subordination in territorial contests for high-quality breeding sites, from behavioural observations to experiments in taxa spanning fish, reptiles, birds, mammals and arthropods (e.g. Baird & Timanus, 1998; Gołąb et al., 2013; Kokko et al., 2004; Newton, 1992; Piper et al., 2000; Stiver et al., 2006; Wauters & Lens, 1995). Other studies have matched age dependence in recruitment patterns with expectations of adaptive queuing for high-quality sites (Ens et al., 1995; van de Pol et al., 2007). While no empirical studies investigating the role of competition for high-quality sites in breeding propensity have previously demonstrated negative density dependence in breeding propensity, such a mechanism has been suggested by several studies that found high recruitment rates subsequent to high adult mortality (e.g. Porter & Coulson, 1987; Pradel et al., 1997; Sæther et al., 2002; Votier et al., 2008). Our study unifies and generalizes these previous findings by providing evidence of negative density dependence in breeding propensities (Table 1; Figure 4), and highlights the importance of competition in shaping breeding propensity at the population level.

In general, negative density dependence of breeding propensity can be mediated through competition for other resources than breeding space, namely food resources. This is an inherent part of competition for breeding sites when food resources are spatially heterogeneous and defended in the breeding territory (e.g. Aho
<figure><img src="image_7.png" /><figcaption>FIGURE 4 Associations between breeding propensity in year *t* (y-axis) and key population features in year *t* − 1 (x-axis). Different breeding propensities were considered depending on the individual's status in year *t* − 1: breeding rate for former breeders (a, b) and skippers (c), integrative recruitment rate for pre-breeders (d, e) and the immigration rate for immigrants (f, g). The key population features considered were population breeding success (mean number of fledglings per nest (a, e)) and number of competitors (number of breeders, b, c, d, f, or number of present non-breeders, g). Relationships presented here are those with strong evidence for a positive or negative partial correlation (Table 1; see Appendix S5.4 for other relationships). These are partial residual plots representing partial correlations controlling for the remaining population features (e.g. in panel (a) the partial correlation between the integrative recruitment rate and number of breeders controls for population breeding success and the number of present non-breeders); residuals were centred on the variable mean to rescale variation within the original range. Points: posterior means of rescaled residuals; segments: 95%CRIs. Solid line: posterior mean of regression line; grey background: 95%CRI</figcaption></figure>

et al., 1999; Ens et al., 1995; Wauters & Lens, 1995). Alternatively, the
link between competition for food resources and for breeding hab-
itat is loose or inexistent when food resources are limited but non-
defendable (as in central place foragers such as Kittiwakes) or spatially

homogeneous across the breeding habitat. Nonetheless, if food is a
limiting resource for which individuals compete independently of the
breeding site, density should also be negatively associated with popu-
lation breeding success (e.g. Arcese & Smith, 1988; Layton-Matthews
<table><caption>TABLE 1 Summary of the associations between status-specific breeding propensity (rows) and key population features: population breeding success or numbers of competitors (columns)</caption><thead><tr><th colspan="2">Breeding propensity (year t)</th><th colspan="3">Population feature (year t - 1)</th></tr><tr><th>Former status (year t - 1)</th><th>Parameter (year t)</th><th>Population breeding success</th><th>Number of breeders</th><th>Number of present non-breeders</th></tr></thead><tbody><tr><td>Breeder</td><td>Breeding rate of former breeders</td><td>0.38 [0.21, 0.55] (1.00)</td><td>-0.46 [-0.64, -0.27] (1.00)</td><td>0.00 [-0.18, 0.18] (0.51)</td></tr><tr><td>Skipper</td><td>Breeding rate of former skippers</td><td>-0.09 [-0.39, 0.21] (0.72)</td><td>-0.40 [-0.64, -0.15] (1.00)</td><td>-0.11 [-0.37, 0.15] (0.80)</td></tr><tr><td>Pre-breeder</td><td>Integrative recruitment rate</td><td>0.34 [0.20, 0.48] (1.00)</td><td>-0.18 [-0.37, 0.01] (0.97)</td><td>-0.04 [-0.18, 0.10] (0.72)</td></tr><tr><td>Immigrant</td><td>Immigration rate</td><td>0.07 [-0.21, 0.35] (0.70)</td><td>-0.26 [-0.52, 0.01] (0.96)</td><td>-0.35 [-0.59, -0.09] (0.99)</td></tr></tbody></table>

Notes: coefficients are partial correlations controlling for the confounding effect of the remaining population features (e.g. the partial correlation between immigration rate and population breeding success controls for the number of breeders and number of present non-breeders). Estimates are posterior means with 95%CRI between brackets, and the proportion of the posterior distribution that had the same sign as the posterior mean (P) between parentheses. Relationships with strong evidence of correlation with a given sign (P > 0.95) are highlighted in bold.

et al., 2019; Wauters et al., 2004). Here we controlled for the effect of population breeding success when estimating the correlation between the number of competitors and breeding propensity (Table 1; Figure 4), which is why the observed relationship is to be explained by competition for breeding habitats, not for food.

The energetic requirements of reproduction and food intake remain a major potential determinant of breeding propensity. And indeed, it has been shown that improved environmental conditions implying lower energetic demand or simply increased food availability are associated with both increased breeding propensity and breeding success (e.g. Hoy et al., 2016; Nur & Sydeman, 1999). Following this view, previous studies documenting positive relationships between population habitat quality and subsequent breeding propensity have referred to physiological condition or perceived chances to overcome reproductive costs (e.g. Cayuela et al., 2018; Frederiksen & Bregnballe, 2001). However, competition for breeding sites on its own can also generate a positive association between population breeding success and breeding propensity, as found in our study (Table 1; Figure 4). Increased mean population breeding success can reflect increased habitat quality in multiple patches across the population, which results in decreased spatial heterogeneity of habitat quality across space (Figure 1; Appendix S1) and decreased competition for high-quality sites. Competition for high-quality breeding sites would appear to better explain temporal variation in breeding propensity than energetic requirements in systems where temporal variation in habitat quality is spatially heterogeneous rather than homogeneous. Where possible in the future, the relative importance and joint contribution of these two explanations could be addressed by analyses that explicitly distinguish between these two forms of variation (e.g. using unambiguous measures of food availability or experiments relying on supplementary feeding).

Spatial heterogeneity of the environment, attraction to high-quality habitats and competition for breeding space are commonplace in animal taxa and are the basis of theory on spatial distribution of individuals (Fretwell & Lucas, 1969; Pulliam & Danielson, 1991). Our hypothesis linking habitat selection mechanisms and breeding propensity should thus be of general relevance, but its importance should depend on key factors underlying competition intensity. For example, the use of information on habitat quality emanating from

conspecifics (e.g. their breeding success) is a common habitat selection mechanism that necessarily makes individuals aggregate and covet the same sites, enhancing competition (Danchin et al., 1998; Doligez et al., 2003; Schmidt et al., 2010). However, such a process depends on predictability and spatial heterogeneity in habitat quality (Acker et al., 2018; Doligez et al., 2003): the more predictable (i.e. temporally autocorrelated) and heterogeneous the habitats (i.e. stronger site-dependent differences in fitness prospects), the higher the competition for high-quality breeding sites. The strength of competition will also depend on the degree to which the availability of high-quality sites is limited. The limitation as perceived by individuals will be conditioned by the type of information used to assess habitat quality and the overall strategy for habitat search (Acker et al., 2017; Lima & Zollner, 1996; Piper, 2011; Rushing et al., 2021; Schmidt et al., 2010), as well as by any physical limitation in the number of breeding sites (Kokko & Sutherland, 1998). Further, competition has led to the evolution of territorial behaviour characterized by costly defence and active contests for exclusive space suitable for breeding (Adams, 2001; Stamps, 1994). By modulating the benefits of occupying a high-quality site versus a low-quality one through related costs of site acquisition and retention in face of competitors, key features of territorial behaviour like territory size and reducibility or risk of injury should modulate the influence of competition for high-quality sites on breeding propensity (Kokko & Sutherland, 1998; López-Sepulcre & Kokko, 2005).

Our study system provides a valuable example of the behavioural and environmental characteristics leading to strong competition for high-quality sites and of its consequences for breeding propensity. Several studies have shown that kittiwakes breed in spatio-temporally heterogeneous but predictable habitats, use public information to target high-quality habitats at all spatial scales, devote substantial time and energy to acquiring and defending breeding sites and show positive associations between breeding propensity and competitive behaviour claiming territorial dominance (Acker et al., 2017; Aubry et al., 2009; Boulinier et al., 2008; Cadiou et al., 1994; Cam et al., 2002; Danchin et al., 1998; Appendices S1.2 and S5.5). Kittiwakes feed on non-defendable resources that vary at regional scales, and although food availability can affect their reproductive success (Frederiksen et al., 2005; Golet et al., 2004; Suryan et al., 2002), previous studies
have not found evidence of effects of food availability on breeding
propensity (Golet et al., 2004; Oro & Furness, 2002) or of density
dependence mediated by food limitation (Frederiksen et al., 2005).
Studies of breeding propensity in other taxa would be valuable to fur-
ther clarify the role that competition for high-quality breeding sites
could play in shaping breeding propensity.

In our study population, competitive asymmetries among in-
dividuals in different breeding status likely modulate how com-
petition for high-quality sites influences breeding propensity. In
general, the most competitive individuals are assumed to occupy
the highest quality habitats and force others to settle in lower qual-
ity habitats (Fretwell & Lucas, 1969; Pulliam & Danielson, 1991) or
to skip breeding (Ens et al., 1995; Piper et al., 2000). Our results sug-
gest that breeders and skippers might benefit from a lower number
of breeding competitors in the population to a greater extent than
pre-breeders and immigrants, and that immigrants might be the
only status affected by competition with non-breeders (Table 1;
Figure 4). This probably reflects the lack of behavioural maturity
of pre-breeders compared to experienced individuals, which would
lead to inferiority of many pre-breeders under any competitive in-
tensity (Aubry et al., 2009; Cam et al., 2002), and to an even greater
inferiority of immigrants due to their lack of knowledge and famil-
iarity with the local competitive context (e.g. Germain et al., 2017).
Our results suggest that decreased competition for high-quality
sites when population breeding success increased benefited pre-
breeders and breeders, while this was not clear for skippers and
immigrants (Table 1; Figure 4). This could be because skippers and
immigrants tend to target less attractive sites (located in patches of
lower quality; e.g. Bruinzeel, 2007) where their chances of acquir-
ing a site are not (or less) impacted by variation in the availability of
high-quality sites across the population.

Overall, the process of competition for high-quality breeding sites emphasized here may explain some major variations in individual life histories. Through despotism exercised by some individuals that manage to breed in high-quality habitats, less competitive ones are forced to poorer reproductive careers (e.g. Bruinzeel, 2007; van de Pol et al., 2007). In our study population, out-competed Kittiwakes could skip breeding opportunities (Cadiou et al., 1994; Cam et al., 2002) or access lower quality breeding sites (Aubry et al., 2009), where they are likely to fail and then disperse to avoid failing again (Acker et al., 2017), re-enter competition to obtain a new site and repeat this cycle ('the spiral of failure'; Cam et al., 2004, 2013). However, our results suggest that higher population habitat quality or lower density may soften competition by offering better breeding opportunities or enhanced access to good opportunities. The same mechanisms should also affect population dynamics. The positive association between previous population breeding success and breeding propensity should accentuate the impacts of temporal variation in habitat quality on population growth (Brown et al., 2000; Danchin et al., 1998). Yet, given the negative association between breeding propensity and the number of competitors, the impact of habitat quality is likely to be counteracted by the variation in competition intensity arising from breeding density changes. These results

open valuable future opportunities to evaluate the relative impor-
tance of competition for high-quality breeding sites in amplifying or
buffering population dynamics via breeding propensity.

ACKNOWLEDGEMENTS

The authors thank everyone involved in data collection, especially Jacques Nisser and Morgane Huteau. This work was supported by the LabEx TULIP (ANR-10-LABX-41; ANR-11-IDEX-0002-02). Analyses were performed using CALMIP supercomputing facilities (allocation P1427). The Institut Universitaire Européen de la Mer provided financial support for this long-term study. They are grateful to the Petzl and Courant companies for sponsoring working-at-height equipment. They thank the Conseil Départemental du Finistère for allowing them to work in the Pointe du Raz and the nature reserve of Goulien Cap Sizun (managed by Bretagne Vivante).

CONFLICT OF INTEREST

The authors have no conflict of interest to declare.

AUTHORS' CONTRIBUTIONS

P.A. and E.C. formulated the ideas; J.-Y.M. conceived the monitoring study; E.C. and J.-Y.M. collected the data; P.A. and M.S. designed the modelling methodology; P.A. analysed the data and led the writing of the manuscript, assisted by E.C. and A.B. All authors contributed critically to the drafts and gave final approval for publication.

DATA AVAILABILITY STATEMENT

Data are available from the Zenodo Repository https://doi.org/10.5281/zenodo.6009808 (Acker et al., 2022).

ORCID

Paul Acker https://orcid.org/0000-0002-3815-772X

Michael Schaub https://orcid.org/0000-0003-2745-0557

Aurélien Besnard https://orcid.org/0000-0002-2289-9761

Emmanuelle Cam https://orcid.org/0000-0001-7324-6958

REFERENCES

Acker, P., Besnard, A., Monnat, J.-Y., & Cam, E. (2017). Breeding habitat selection across spatial scales: Is grass always greener on the other side? *Ecology*, *98*(10), 2684–2697. https://doi.org/10.1002/ecy.1962

Acker, P., Francesiaz, C., Béchet, A., Sadoul, N., Lessells, C. M., Pijl, A. S., & Besnard, A. (2018). Insights on dispersal and recruitment paradigms: Sex-and age-dependent variations in a nomadic breeder. *Oecologia*, *186*(1), 85–97. https://doi.org/10.1007/s0044 2-017-3972-7

Acker, P., Schaub, M., Besnard, A., Monnat, J.-Y., & Cam, E. (2022). Data from: Can attraction to and competition for high-quality habitats shape breeding propensity? *Zenodo*. https://doi.org/10.5281/ zenodo.6009808

Adams, E. S. (2001). Approaches to the study of territory size and shape. *Annual Review of Ecology and Systematics*, 32(1), 277–303. https://doi.org/10.1146/annurev.ecolsys.32.081501.114034

Aho, T., Kuitunen, M., Suhonen, J., Jäntti, A., & Hakkari, T. (1999). Reproductive success of Eurasian treecreepers, *Certhia familiaris*, lower in territories with wood ants. *Ecology*, 80(3), 998–1007.
https://doi.org/10.1890/0012-9658(1999)080[0998:RSOETC]2.0.CO;2

Arcese, P., & Smith, J. N. M. (1988). Effects of population density and supplemental food on reproduction in song sparrows. *Journal of Animal Ecology*, 57(1), 119-136. https://doi.org/10.2307/4768

Aubry, L. M., Cam, E., & Monnat, J.-Y. (2009). Habitat selection, age-specific recruitment and reproductive success in a long-lived seabird. In D. L. Thomson, E. G. Cooch, & M. J. Conroy (Eds.), *Modeling demographic processes in marked populations* (pp. 365-392). Springer. https://doi.org/10.1007/978-0-387-78151-8_16

Baird, T. A., & Timanus, D. K. (1998). Social inhibition of territorial behaviour in yearling male collared lizards, *Crotaphytus collaris*. *Animal Behaviour*, 56(4), 989-994. https://doi.org/10.1006/anbe.1999.1295

Bled, F., Royle, J. A., & Cam, E. (2011). Assessing hypotheses about nesting site occupancy dynamics. *Ecology*, 92(4), 938-951. https://doi.org/10.1890/10-0392.1

Boulinier, T., McCoy, K. D., Yoccoz, N. G., Gasparini, J., & Tveraa, T. (2008). Public information affects breeding dispersal in a colonial bird: Kittiwakes cue on neighbours. *Biology Letters*, 4(5), 538-540. https://doi.org/10.1098/rsbl.2008.0291

Brown, C. R., Brown, M. B., & Danchin, E. (2000). Breeding habitat selection in cliff swallows: The effect of conspecific reproductive success on colony choice. *Journal of Animal Ecology*, 69(1), 133-142. https://doi.org/10.1046/j.1365-2656.2000.00382.x

Bruinzeel, L. W. (2007). Intermittent breeding as a cost of site fidelity. *Behavioral Ecology and Sociobiology*, 61(4), 551-556. https://doi.org/10.1007/s00265-006-0283-0

Cadiou, B., Monnat, J. Y., & Danchin, E. (1994). Prospecting in the kittiwake, *Rissa tridactyla*: Different behavioural patterns and the role of squatting in recruitment. *Animal Behaviour*, 47(4), 847-856. https://doi.org/10.1006/anbe.1994.1116

Cam, E., Cadiou, B., Hines, J. E., & Monnat, J.-Y. (2002). Influence of behavioural tactics on recruitment and reproductive trajectory in the kittiwake. *Journal of Applied Statistics*, 29(1-4), 163-185. https://doi.org/10.1080/02664760120108502

Cam, E., Gimenez, O., Alpizar-Jara, R., Aubry, L. M., Authier, M., Cooch, E. G., Koons, D. N., Monnat, J.-Y., Nichols, J. D., Rotella, J. J., Royle, J. A., & Pradel, R. (2013). Looking for a needle in a haystack: Inference about individual fitness components in a heterogeneous population. *Oikos*, 122(5), 739-753. https://doi.org/10.1111/j.1600-0706.2012.20532.x

Cam, E., Hines, J. E., Monnat, J.-Y., Nichols, J. D., & Danchin, E. (1998). Are adult nonbreeders prudent parents? *The kittiwake model*. *Ecology*, 79(8), 2917-2930. https://doi.org/10.2307/176526

Cam, E., Monnat, J., & Royle, A. (2004). Dispersal and individual quality in a long-lived species. *Oikos*, 106(2), 386-398. https://doi.org/10.1111/j.0030-1299.2003.13097.x

Caswell, H. (2001). *Matrix population models* (2nd ed.). John Wiley & Sons.

Cayuela, H., Besnard, A., Bonnaire, E., Perret, H., Rivoalen, J., Miaud, C., & Joly, P. (2014). To breed or not to breed: Past reproductive status and environmental cues drive current breeding decisions in a long-lived amphibian. *Oecologia*, 176(1), 107-116. https://doi.org/10.1007/s00442-014-3003-x

Christensen-Dalsgaard, S., May, R., & Lorentsen, S. H. (2018). Taking a trip to the shelf: Behavioral decisions are mediated by the proximity to foraging habitats in the black-legged kittiwake. *Ecology and Evolution*, 8(2), 866-878. https://doi.org/10.1002/ece3.3700

Cullen, E. (1957). Adaptations in the kittiwake to cliff-nesting. *Ibis*, 99(2), 275-302. https://doi.org/10.1111/j.1474-919X.1957.tb01950.x

Danchin, E., Boulinier, T., & Massot, M. (1998). Conspecific reproductive success and breeding habitat selection: Implications for the study of coloniality. *Ecology*, 79(7), 2415-2428. https://doi.org/10.2307/176832

Danchin, E., & Cam, E. (2002). Can non-breeding be a cost of breeding dispersal? *Behavioral Ecology and Sociobiology*, 51(2), 153-163. https://doi.org/10.1007/s00265-001-0423-5

Desprez, M., Gimenez, O., McMahon, C. R., Hindell, M. A., & Harcourt, R. G. (2018). Optimizing lifetime reproductive output: Intermittent breeding as a tactic for females in a long-lived, multiparous mammal. *Journal of Animal Ecology*, 87(1), 199-211. https://doi.org/10.1111/1365-2656.12775

Doligez, B., Cadet, C., Danchin, E., & Boulinier, T. (2003). When to use public information for breeding habitat selection? The role of environmental predictability and density dependence. *Animal Behaviour*, 66(5), 973-988. https://doi.org/10.1006/anbe.2002.2270

Ens, B. J., Weissing, F. J., & Drent, R. H. (1995). The despotic distribution and deferred maturity: Two sides of the same coin. *The American Naturalist*, 146(4), 625-650. https://doi.org/10.1086/285818

Erikstad, K. E., Fauchald, P., Tveraa, T., & Steen, H. (1998). On the cost of reproduction in long-lived birds: The influence of environmental variability. *Ecology*, 79(5), 1781-1788. https://doi.org/10.2307/176796

Frederiksen, M., & Bregnballe, T. (2001). Conspecific reproductive success affects age of recruitment in a great cormorant, *Phalacrocorax carbo sinensis*, colony. *Proceedings of the Royal Society of London. Series B: Biological Sciences*, 268(1475), 1519-1526. https://doi.org/10.1098/rspb.2001.1701

Frederiksen, M., Wright, P. J., Harris, M. P., Mavor, R. A., Heubeck, M., & Wanless, S. (2005). Regional patterns of kittiwake Rissa tridactyla breeding success are related to variability in sandeel recruitment. *Marine Ecology Progress Series*, 300, 201-211. https://doi.org/10.3354/meps300201

Fretwell, S. D., & Lucas, H. L. (1969). On territorial behavior and other factors influencing habitat distribution in birds I. Theoretical development. *Acta Biotheoretica*, 19(1), 16-36. https://doi.org/10.1007/BF01601953

Gelman, A., & Hill, J. (2006). *Data analysis using regression and multilevel/hierarchical models*. Cambridge University Press.

Germain, M., Pärt, T., & Doligez, B. (2017). Lower settlement following a forced displacement experiment: Nonbreeding as a dispersal cost in a wild bird? *Animal Behaviour*, 133, 109-121. https://doi.org/10.1016/j.anbehav.2017.09.001

Giudici, A., Navarro, J., Juste, C., & González-Solis, J. (2010). Physiological ecology of breeders and sabbaticals in a pelagic seabird. *Journal of Experimental Marine Biology and Ecology*, 389(1), 13-17. https://doi.org/10.1016/j.jembe.2010.04.002

Golet, G. H., Schmutz, J. A., Irons, D. B., & Estes, J. A. (2004). Determinants of reproductive costs in the long-lived black-legged kittiwake: A multiyear experiment. *Ecological Monographs*, 74(2), 353-372. https://doi.org/10.1890/02-4029

Hoy, S. R., Millon, A., Petty, S. J., Whitfield, D. P., & Lambin, X. (2016). Food availability and predation risk, rather than intrinsic attributes, are the main factors shaping the reproductive decisions of a long-lived predator. *Journal of Animal Ecology*, 85(4), 892-902. https://doi.org/10.1111/1365-2656.12517

Johnson, M. D. J. (2007). Measuring habitat quality: A review. *The Condor*, 109(3), 489-504. https://doi.org/10.1650/8347.1

Kokko, H., Harris, M. P., & Wanless, S. (2004). Competition for breeding sites and site-dependent population regulation in a highly colonial seabird, the common guillemot *Uria aalge*. *Journal of Animal Ecology*, 73(2), 367-376. https://doi.org/10.1111/j.0021-8790.2004.00813.x

Kokko, H., & Sutherland, W. J. (1998). Optimal floating and queuing strategies: Consequences for density dependence and habitat loss. *The American Naturalist*, 152(3), 354-366. https://doi.org/10.1086/286174

Layton-Matthews, K., Loonen, M. J. J. E., Hansen, B. B., Coste, C. F. D., Sæther, B.-E., & Grøtan, V. (2019). Density-dependent population dynamics of a high Arctic capital breeder, the barnacle goose. *Journal of Animal Ecology*, 88(8), 1191-1201. https://doi.org/10.1111/1365-2656.13001

Lee, A. M., Reid, J. M., & Beissinger, S. R. (2016). Modelling effects of nonbreeders on population growth estimates. *Journal of Animal Ecology*, 86(1), 75-87. https://doi.org/10.1111/1365-2656.12592
Lima, S. L., & Zollner, P. A. (1996). Towards a behavioral ecology of eco-
logical landscapes. *Trends in Ecology & Evolution*, *11*(3), 131–135.
https://doi.org/10.1016/0169-5347(96)81094-9

Link, W. A., Cooch, E. G., & Cam, E. (2002). Model-based estimation of individual fitness. *Journal of Applied Statistics*, *29*(1-4), 207-224.
https://doi.org/10.1080/02664760120108700a

López-Sepulcre, A., & Kokko, H. (2005). Territorial defense, territory size, and population regulation. *The American Naturalist*, *166*(3), 317-329. https://doi.org/10.1086/432560

Meijer, T., & Drent, R. (1999). Re-examination of the capital and income dichotomy in breeding birds. *Ibis*, 141(3), 399-414. https://doi.org/10.1111/j.1474-919X.1999.tb04409.x

Newton, I. (1992). Experiments on the limitation of bird numbers by territorial behaviour. *Biological Reviews*, 67(2), 129-173. https://doi.org/10.1111/j.1469-185X.1992.tb01017.x

Nur, N., & Sydeman, W. J. (1999). Survival, breeding probability and reproductive success in relation to population dynamics of Brandt's cormorants *Phalacrocorax penicillatus*. *Bird Study*, 46(sup1), S92-S103. https://doi.org/10.1080/00063659909477236

Oro, D., & Furness, R. W. (2002). Influences of food availability and pre-
dation on survival of kittiwakes. *Ecology*, *83*(9), 2516–2528. https://
doi.org/10.1890/0012-9658(2002)083[2516:IOFAAP]2.0.CO;2

Piper, W. H. (2011). Making habitat selection more "familiar": A review. *Behavioral Ecology and Sociobiology*, *65*(7), 1329–1351. https://doi.org/10.1007/s00265-011-1195-1

Piper, W. H., Tischler, K. B., & Klich, M. (2000). Territory acquisition in loons: The importance of take-over. *Animal Behaviour*, *59*(2), 385-394. https://doi.org/10.1006/anbe.1999.1295

Plummer, M. (2003). JAGS: A program for analysis of Bayesian graphical
models using Gibbs sampling. *Proceedings of the 3rd International
Workshop on Distributed Statistical Computing*, 124, 125. Retrieved
from https://www.r-project.org/conferences/DSC-2003/Drafts/
Plummer.pdf

Porter, J. M., & Coulson, J. C. (1987). Long-term changes in recruitment to the breeding group, and the quality of recruits at a kittiwake *Rissa tridactyla* colony. *Journal of Animal Ecology*, *56*(2), 675-689. https://doi.org/10.2307/5076

Pradel, R., Johnson, A. R., Viallefont, A., Nager, R. G., & Cézilly, F. (1997). Local recruitment in the greater flamingo: A new approach using capture-mark-recapture data. *Ecology*, *78*(5), 1431-1445. https://doi.org/10.1890/0012-9658(1997)078[1431:LRITGF]2.0.CO;2

Pulliam, H. R., & Danielson, B. J. (1991). Sources, sinks, and habi-
tat selection: A landscape perspective on population dynamics.
*The American Naturalist*, 137, S50-S66. https://doi.org/10.1086/
285139

Rodeout, R. M., & Tomkiewicz, J. (2011). Skipped spawning in fishes: more common than you might think. *Marine and Coastal Fisheries*, *3*(1), 176-189. https://doi.org/10.1080/19425120.2011.556943

Rushing, C. S., Ryder, T. B., Valente, J. J., Sillett, T. S., & Marra, P. P. (2021). Empirical tests of habitat selection theory reveal that conspecific density and patch quality, but not habitat amount, drive long-distance immigration in a wild bird. *Ecology Letters*, *24*(6), 1167-1177. https://doi.org/10.1111/ele.13729

Sæther, B.-E., Engen, S., & Matthysen, E. (2002). Demographic character-
istics and population dynamical patterns of solitary birds. *Science*,
295(5562), 2070-2073. https://doi.org/10.1126/science.1068766

Schaub, M., & Abadi, F. (2011). Integrated population models: A novel analysis framework for deeper insights into population dynamics. *Journal of Ornithology*, *152*(1), 227-237. https://doi.org/10.1007/s10336-010-0632-7

Schaub, M., & Fletcher, D. (2015). Estimating immigration using a Bayesian integrated population model: Choice of parametrization and priors. *Environmental and Ecological Statistics*, *22*, 535-549. https://doi.org/10.1007/s10651-015-0309-8

Schaub, M., & Kéry, M. (2021). *Integrated population models. Theory and ecological applications using R and JAGS* (1st ed.). Academic Press.

Schmidt, K. A., Dall, S. R., & Van Gils, J. A. (2010). The ecology of in-
formation: An overview on the ecological significance of mak-
ing informed decisions. *Oikos*, *119*(2), 304–316. https://doi.
org/10.1111/j.1600-0706.2009.17573.x

Shine, R., & Brown, G. P. (2008). Adapting to the unpredictable: Reproductive biology of vertebrates in the Australian wet-dry tropics. *Philosophical Transactions of the Royal Society B: Biological Sciences*, *363*(1490), 363–373. https://doi.org/10.1098/rstb.2007.2144

Stamps, J. (1994). Territorial behavior: Testing the assumptions. *Advances in the Study of Behavior*, 23, 173-232 http://10.0.3.248/S0065-3454(08)60354-X

Stearns, S. C. (1992). *The evolution of life histories*. Oxford University Press.

Stiver, K. A., Fitzpatrick, J., Desjardins, J. K., & Balshine, S. (2006). Sex differences in rates of territory joining and inheritance in a cooperatively breeding cichlid fish. *Animal Behaviour*, *71*(2), 449-456. https://doi.org/10.1016/j.anbehav.2005.06.011

Suryan, R. M., Irons, D. B., Kaufman, M., Benson, J., Jodice, P. G., Roby, D. D., & Brown, E. D. (2002). Short-term fluctuations in forage fish availability and the effect on prey selection and brood-rearing in the black-legged kittiwake *Rissa tridactyla*. *Marine Ecology Progress Series*, *236*, 273-287. https://doi.org/10.3354/meps236273

Szostek, K. L., Schaub, M., & Becker, P. H. (2014). Immigrants are attracted by local pre-breeders and recruits in a seabird colony. *Journal of Animal Ecology*, *83*(5), 1015-1024. https://doi.org/10.1111/1365-2656.12206

van de Pol, M., Pen, I., Heg, D., & Weissing, F. J. (2007). Variation in habitat choice and delayed reproduction: Adaptive queuing strategies or individual quality differences? *The American Naturalist*, *170*(4), 530-541. https://doi.org/10.1086/521237

Votier, S. C., Birkhead, T. R., Oro, D., Trinder, M., Grantham, M. J., Clark, J. A., McCleery, R. H., & Hatchwell, B. J. (2008). Recruitment and survival of immature seabirds in relation to oil spills and climate variability. *Journal of Animal Ecology*, *77*(5), 974-983. https://doi.org/10.1111/j.1365-2656.2008.01421.x

Wauters, L. A., & Lens, L. (1995). Effects of food availability and density on red squirrel (*Sciurus vulgaris*) reproduction. *Ecology*, *76*(8), 2460-2469. https://doi.org/10.2307/2265820

Wauters, L. A., Matthysen, E., Adriaensen, F., & Tosi, G. (2004). Within-sex density dependence and population dynamics of red squirrels *Sciurus vulgaris*. *Journal of Animal Ecology*, *73*(1), 11-25. https://doi.org/10.1111/j.1365-2656.2004.00792.x

Weegman, M. D., Arnold, T. W., Clark, R. G., & Schaub, M. (2021). Partial and complete dependency among data sets has minimal consequence on estimates from integrated population models. *Ecological Applications*, *31*(3), e2258. https://doi.org/10.1002/eap.2258

Zack, S., & Stutchbury, B. J. (1992). Delayed breeding in avian social systems: The role of territory quality and "floater" tactics. *Behaviour*, *123*(3-4), 194-219. https://doi.org/10.1163/156853992X00020

SUPPORTING INFORMATION

Additional supporting information may be found in the online
version of the article at the publisher's website.

How to cite this article: Acker, P., Schaub, M., Besnard, A., Monnat, J-Y, & Cam, E. (2022). Can attraction to and competition for high-quality habitats shape breeding propensity?. *Journal of Animal Ecology*, *91*, 933-945. https://doi.org/10.1111/1365-2656.13676