
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70295
ECE370295
ECE-2024-05-00886.R1
Conservation Ecology
Life History Ecology
Population Ecology
Zoology
Research Article
Research Article
Moderate evidence for the sex‐dependent effect of poisoning on adult survival in a long‐lived raptor species
Zsinka et al.
Zsinka Bernadett https://orcid.org/0000-0002-8532-6688
1
Pásztory‐Kovács Szilvia 1 kovacs.szilvia@univet.hu

Kövér Szilvia https://orcid.org/0000-0002-5762-1022
1
Vili Nóra https://orcid.org/0000-0001-7241-231X
1
Horváth Márton https://orcid.org/0000-0002-9503-4128
2
1 Department of Zoology University of Veterinary Medicine Budapest Budapest Hungary
2 MME BirdLife Hungary Budapest Hungary
* Correspondence
Szilvia Pásztory‐Kovács, Department of Zoology, University of Veterinary Medicine Budapest, Rottenbiller utca 50, H‐1077 Budapest, Hungary.
Email: kovacs.szilvia@univet.hu

18 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7029511 8 2024
06 5 2024
29 8 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Survival rate is usually the greatest contributor to population growth in long‐lived species, and its accurate estimation along with the evaluation of the factors influencing it is therefore essential for effective conservation. Here, we studied the survival of breeding eastern imperial eagles Aquila heliaca in Hungary between 2011 and 2022 and investigated the effect of poisoning, the leading known anthropogenic cause of mortality. We used the Cormack‐Jolly‐Seber mark‐recapture model to estimate annual apparent survival and encounter probabilities based on the capture histories of 208 males and 411 females. We obtained encounter data from the DNA profiles of shed feathers collected at the nest sites, which we also supplemented with presences inferred from parentage analysis. The most supported model estimated a constant 91.6% annual survival over the study period, but models including the effect of sex and poisoning rate on survival had similar support. Sex difference in survival was less than 1% on average, but the survival of males decreased more with poisoning rate than the survival of females. However, due to smaller encounter probabilities, the estimates for males were less precise compared to females. Males may be more at risk from poisoning than females not only due to their more active foraging behaviour during incubation and chick‐rearing but also due to their smaller body size. Apart from providing direct practical information for the conservation management of imperial eagles, our results also highlight the importance of long‐term studies for estimating population parameters of birds of prey.

Survival estimation of adult imperial eagles revealed that poisoning activity may have a stronger negative effect on the survival of the smaller‐bodied, more actively foraging males than on females. The genetic‐based mark‐recapture analysis was conducted on data of over 600 birds breeding in Hungary between 2011 and 2022.

bird of prey
conservation genetics
eastern imperial eagle
mark‐recapture
non‐invasive sampling
Doctoral School of Veterinary Science, University of Veterinary Medicine BudapestLIFE programme 10.13039/100013288 LIFE10NAT/HU/019 LIFE15 NAT/HU/000902 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:18.09.2024
Zsinka, B. , Pásztory‐Kovács, S. , Kövér, S. , Vili, N. , & Horváth, M. (2024). Moderate evidence for the sex‐dependent effect of poisoning on adult survival in a long‐lived raptor species. Ecology and Evolution, 14 , e70295. 10.1002/ece3.70295
==== Body
pmc1 INTRODUCTION

Survival is one of the key life‐history traits affecting the population growth rate (Begon et al., 1996) and, consequently, the viability of a population (Sibly & Hone, 2002). Therefore, the accurate estimation of age‐, life‐stage‐ or sex‐specific survival rates is of great importance in understanding the dynamics of a population (Hernández‐Matías et al., 2011; Heuck et al., 2020). In addition, assessing the changes in these demographic parameters in response to certain human pressures, such as habitat alteration, electrocution or persecution, can facilitate the effective conservation of endangered species (Demerdzhiev et al., 2015; López‐López et al., 2011; Selwood et al., 2015; Sergio et al., 2021; Whitfield et al., 2004).

In the case of raptors (Accipitriformes, Falconiformes, Cathartiformes and Strigiformes), survival rate is a critical determinant of population growth: since raptors are generally long‐lived species with late maturity and low fecundity rates, their population growth rates are more sensitive to survival than to reproductive parameters (Katzner et al., 2006; Newton et al., 2016; Sæther & Bakke, 2000; Sergio et al., 2011; Stahl & Oli, 2006). Therefore, studies on raptor survival are especially important, considering that approximately 20% of raptor species are threatened, and more than half of all raptor species have a decreasing global population (McClure et al., 2018).

One of the major threats to raptors is poisoning, which can be intentional (directly aimed at controlling predator species) or accidental (resulting from the misuse of pesticides/rodenticides or from secondary poisoning through consumed poisoned prey) (McClure et al., 2018). Even though poisoning is known to have detrimental effects on populations (Köhler & Triebskorn, 2013), direct evidence linking poisoning to increased mortality and population decline remains scarce (Hernández & Margalida, 2009; Mateo‐Tomás et al., 2020; Ortega et al., 2009; Tenan et al., 2012; Whitfield et al., 2004).

Sex differences in survival are generally small in raptors (Hernández‐Matías et al., 2011; Hunt et al., 2017; Morrison, 2003; Newton et al., 2016), with either males or females having lower survival rates, and the direction of difference can be inconsistent even within a species (Dobson, 1987; Newton et al., 1983). Sex differences in survival rates are considered to be driven mainly by differences in selective forces on life‐history traits resulting in different trade‐off optima between survival and reproduction (dimorphic parental behaviour, defence of territory and body size), while the genetic hypotheses (unguarded X, mother's curse) seem to be less important (Connallon et al., 2022; Maklakov & Lummaa, 2013; Sumasgutner et al., 2019). In raptors, females tend to be larger and spend more time at the nest during the breeding season, while the males are responsible for providing food for their mates and chicks and defending the territory (Newton, 1979; Warkentin et al., 2016). These sex differences may lead to different survival rates in various ways. For example, the larger size of females could promote not only a higher chance of survival during the winter (Payevsky, 2016) but also a larger risk of electrocution (Ferrer & Hiraldo, 1992). Furthermore, the more sedentary behaviour of females during the breeding season could either result in less exposure to anthropogenic mortality factors (Hernández‐Matías et al., 2011) or a higher risk of mortality if they are targeted at the nest (Ewing et al., 2023).

Estimating the survival of raptors can be challenging because they appear in generally low densities and are often difficult to capture (Newton et al., 2016). Hence, the population dynamics of most raptor species are still poorly understood (Buechley et al., 2019; McClure et al., 2022). Survival is often estimated using mark‐recapture models, which require marking individuals and monitoring their fate through live encounters or dead recoveries (Lebreton et al., 1992; Sandercock, 2020). Conventional marking techniques (e.g. rings, wing tags, transmitters) all require the capture of birds, which is often not feasible for breeding adults, especially for large‐sized species. In such cases, DNA profiling from non‐invasively collected samples (e.g. shed feathers from the nest site) can provide an alternative method for identifying and monitoring individuals (Gil‐Sánchez et al., 2021; Horváth et al., 2005; Kenward et al., 2007; Rudnick et al., 2005; Vili, Szabó, et al., 2013). Another genetic‐based identification used to estimate the turnover of breeding individuals utilises DNA samples from chicks (Kylmänen et al., 2023; Ponnikas et al., 2017). In this case, the turnover of breeders is estimated based on relatedness between chicks of the same nest from consecutive years. While the plucked feathers and blood obtained from chicks are considered more reliable sources of DNA than shed feathers (Bush et al., 2005), this technique requires a high‐resolution marker set for reliable relatedness estimation.

Here, we investigate the survival of breeding individuals in the eastern imperial eagle Aquila heliaca (hereafter ‘imperial eagle’). As a Vulnerable (BirdLife International, 2019), large‐sized, long‐lived raptor, the survival rate of breeding individuals is a particularly important demographic parameter for the species. The imperial eagle's distribution is scattered throughout the Palearctic region, with only a few thousand breeding pairs worldwide (BirdLife International, 2019). Its westernmost population is found in the Pannonian Region (Demerdzhiev, Horváth, et al., 2011; Horváth et al., 2002) and consists of ca. 356–381 nesting pairs as of 2019 (Karyakin, 2020), making it the largest unified population outside of Russia and Kazakhstan (Horváth et al., 2011). Most of these Pannonian population (287 pairs as of 2019) reside in Hungary (Horváth et al., 2020). Imperial eagles exhibit floater behaviour in their first years and usually start breeding in their third or fourth calendar year (Horváth, 2022). The breeding success of these immature birds is lower than that of adults (Horváth et al., 2014). Following 1–2 years of successful breeding, adult breeders display high territory and mate fidelity (Rudnick et al., 2005; Vili, Szabó, et al., 2013). Sexual dimorphism in eastern imperial eagles only involves differences in size, with females being the larger sex, and behavioural dimorphism during the breeding season. Similarly to its sister species, the Spanish imperial eagle Aquila adalberti (Margalida, González, Sánchez, Oria, & Prada, 2007), eastern imperial eagle males also generally spend less time around the nest than females since they rarely take part in incubation, and they play the major role in food provisioning during the chick‐rearing period (Dobrev, 2009; Horváth, 2022). The imperial eagle is threatened by various anthropogenic factors, such as habitat fragmentation and alteration (Demerdzhiev et al., 2022), electrocution and persecution, including the illegal poisoning and shooting of birds (Deák, Fatér, et al., 2020; Lazarova et al., 2020). Even though the most dangerous pesticides (e.g. carbofuran) were banned in the EU in 2008, poisoning cases due to these substances still occur in imperial eagles, as well as in several other raptors (Deák, Árvay, & Horváth, 2020; Kitowski et al., 2020). In Hungary, poisoning was the leading known cause of mortality between 2005 and 2019, representing 25%–35% of all detected mortality cases (Deák, Fatér, et al., 2020; Horváth et al., 2011). Most of these were the result of intentional poisoning incidents, where baits (such as chicken eggs, carcasses of smaller prey animals or parts of large animals) poisoned with legally banned pesticides or insecticides were deployed with the aim of eliminating avian or mammal predators (Deák, Árvay, & Horváth, 2020). The other cases resulted from accidental poisoning, which occurs from the misuse of chemicals aimed to control agricultural pests or rodents. These were usually the result of improperly installed bait stations and/or the use of banned pesticides, which raptors encounter by consuming poisoned prey animals or their carcasses (Deák, Árvay, & Horváth, 2020; Deák, Fatér, et al., 2020).

We used a mark‐recapture method based on genetic identification to estimate annual apparent survival and encounter probabilities for breeding imperial eagles in East Hungary between 2011 and 2022. We aimed to explore possible sex differences in survival and investigate the relationship between poisoning and the annual survival probabilities of breeding males and females. Due to the aforementioned behavioural dimorphism, we expected males to have lower survival rates than females, as they spend more time away from the nest during breeding and, consequently, were assumed to have a higher risk of encountering anthropogenic mortality factors. Additionally, we hypothesised that the effect of poisoning on survival will be sex‐dependent due to the seasonal patterns of this behavioural dimorphism and that of poisoning activity. Most cases of poisoning during the study years occurred in the first half of the year, with a peak in early spring (February–April), which coincides with the egg‐laying and incubation period (March–May) of the imperial eagle (Deák, Árvay, & Horváth, 2020; Dobrev, 2009; Horváth, 2022). Since this is the period when females rarely leave the nest and only males hunt for food, we expected males to be more exposed to poisoning than females. Besides, males are also significantly smaller in body size than females (even by 20%–25% in weight); therefore, the same amount of a poisonous chemical in a prey or bait could be more detrimental or even fatal for males. Since breeding imperial eagles are difficult to capture and mark with conventional methods, we used DNA profiles obtained from shed feathers for the individual identification of breeders. This method has been previously applied with success in two studies estimating survival in imperial eagles (Rudnick et al., 2005; Vili, Szabó, et al., 2013). In addition, we also utilised the long‐term genetic monogamy of the species to obtain additional presence data of breeding birds through parentage analysis (Rudnick et al., 2005).

2 MATERIALS AND METHODS

2.1 Sample collection

DNA samples for genetic profiling were collected between June and September each year between 2011 and 2022 in East Hungary in the frame of the national monitoring scheme for the imperial eagle in Hungary (Figure 1). The eastern part of the country holds about 95% of the national population of imperial eagles. On average, 67% of the nests here were sampled yearly (SD: 13%), with the highest coverage in 2017 (86%, 191 nests) and the lowest in 2012 (51%, 77 nests). Armpit feathers were plucked from nestlings during ringing, and breeding individuals were non‐invasively sampled by collecting their shed feathers around the 100 m radius of the nest. Plucked feathers were stored at −20°C in 2 mL microtubes filled with 96% ethanol, while shed feathers were stored in tagged plastic bags in dark, dry, and cool conditions to preserve DNA (Vili, Nemesházi, et al., 2013). We processed the majority (90%) of the samples within 1 year and all samples within 3 years of collection.

FIGURE 1 Distribution of the eastern imperial eagle in Hungary between 2011 and 2022 with sampled (n = 369, black) and not sampled (n = 73, grey) nests.

2.2 DNA extraction

We extracted the whole‐genome DNA using the Omega E.Z.N.A.® Tissue DNA Kit (Omega Bio‐tek Inc.) following the manufacturer's instructions but using an additional 20 μL of dithiotreitol (1 M) to aid in the digestion of keratin (Weigmann, 1968). We extracted DNA from the tip of the calamus in plucked nestling feathers and from the superior umbilicus region of shed feathers (Horváth et al., 2005).

2.3 Molecular sexing

We conducted molecular sexing on each feather by amplifying introns of the sex chromosome‐linked CHD1 gene, using the primers CHD‐i16F/CHD‐i16R (Suh et al., 2011). The PCR reaction included 0.065 μL DreamTaq polymerase (Fermentas), 1.7 μL 10X DreamTaq Green Buffer (Fermentas), 0.65 μL 25 mM MgCl2 (Thermo Scientific), 0.65 μL 2 mM dNTP mix (Thermo Scientific), 1–1 μL 10 pmol/μl forward and reverse primers, 8 μL H2O and 4 μL ca. 50 ng/μl concentration DNA. The PCR programme for molecular sexing constituted of an initial denaturation step at 95°C for 2 min, a touchdown section of nine cycles (denaturation: 95°C for 30s, annelation: temperature lowering by 1°C each cycle from 60 to 52°C and lasting 45 s, elongation: 72°C for 45 s), followed by 28 cycles of 95°C for 30s, 52°C for 45 s and 72°C for 45 s, ending with a final elongation of 7 min at 72°C. The PCR products were visualised through gel electrophoresis by UV illumination (2% agarose gel stained with EcoSafe (Pacific Image Electronics Co., Ltd) intercalator, 100 V, 45 min); the heterogametic females display two bands while the homogametic males only one.

2.4 Individual genotyping

We used nine microsatellite markers for individual identification. Out of these, two tetranucleotide loci (IEAAAG09 and IEAAAG11) (Busch et al., 2005) were optimised for the imperial eagle, five dinucleotide loci (Aa02, Aa35, Aa36, Aa39 and Aa43) (Martínez‐Cruz et al., 2002) were published for the Spanish imperial eagle Aquila adalberti and two dinucleotide loci (Hal04 and Hal10) (Hailer et al., 2005) for the white‐tailed eagle Haliaeetus albicilla. The 5′ end of the forward primers was modified with the following fluorescent dyes (Applied Biosystems™): 6‐FAM™ for Aa02, Aa39, Aa43, IEAAAG09, IEAAAG11; HEX™ for Aa35 and Aa36; and NED™ for Hal04 and Hal10. The 5′ end of the reverse primer was modified with a pigtail 5′GTTT sequence in the case of Aa02, Aa36 and Aa39.

We performed PCR reactions in a 10 μL volume, containing 2 μL 5xFIREPol® Master Mix (Solis BioDyne), which consists of dNTP mix, MgCl2 and Taq DNA‐polymerase, 0.5–0.5 μL 10 pmol/μl forward and reverse primers, 5 μL H2O and 2 μL ca. 50 ng/μl concentration DNA. Aa36 and Aa39, Aa35 and Aa43, Hal04 and Hal10, IEAAAG09 and IEAAAG11 were also amplifiable as duplexes. We used the PCR procedure described in Martínez‐Cruz et al. (2002) for all Aa loci and applied a modified version of it (touchdown scheme: 66–60°C, annealing at 60°C for 31 cycles) to the IEAAAG loci. For the Hal loci, we used the PCR profile described in Hailer et al. (2005), with the following modifications: 37 cycles, 45 s for both annealing and amplification. Fragment lengths were determined using capillary electrophoresis: PCR products were run on an ABI3130 sequencer (Applied Biosystems, using Gene Scan™ ‐500LIZ™ Size Standard), and we identified and scored alleles with OSIRIS v2.16 (Goor et al., 2011). We performed the fragment analysis similarly to the suggestions of Beja‐Pereria et al. (2009) by scoring each sample three times independently. Genotypes were assigned blind to the origin of the sample. We checked the deviation from the Hardy–Weinberg equilibrium with Genepop v1.2.2. (Rousset, 2008) and investigated the possible occurrence of null alleles and allelic dropouts using microchecker v.2.2.2 (Van Oosterhout et al., 2004). We calculated probabilities of identity (PI and PISIB) (Waits et al., 2001) and exclusion probabilities for parentage analysis (P1X and P2X) with GenAlEx v.6.503 (Peakall & Smouse, 2006).

2.5 Constructing capture histories

We constructed capture histories (yearly presence–absence data for each individual) from two types of presences: (i) direct presences, when a breeding bird was sampled and genetically profiled directly from its shed feathers and (ii) indirect presences, when the presence of a breeding bird was inferred through parentage analysis due to the lack of shed feathers in that specific year. We could only obtain indirect presences for an individual if both the profile of its mate and at least one of its chicks were known from that specific year and the individual in question had also been profiled in another year (preceding or following the year in focus). We performed parentage analyses for indirect presences manually.

2.6 Survival analysis

We conducted the survival analysis in the MARK v9.0 software (White & Burnham, 1999) using the RMark v3.0.0 R interface (Laake, 2013) by fitting the Cormack‐Jolly‐Seber (CJS) open population model (Amstrup et al., 2006) to the capture histories. The CJS model estimates annual apparent survival probabilities (ϕ) for each one‐year interval and encounter probabilities (p) for each year as parameters of a generalised linear model. As sampling took place in June each year, these apparent survival probabilities refer to one‐year intervals. The estimated survival probability is apparent by definition, as it is the joint probability of surviving and remaining in the sampling area. However, since imperial eagles display high territory fidelity between years (Rudnick et al., 2005; Vili, Szabó, et al., 2013), we considered the apparent survival as a good proxy for survival in this case. Therefore, we will refer to ‘apparent survival’ as ‘survival’ throughout the article.

We constructed a candidate model set based on previous knowledge of the factors potentially influencing survival and encounter probabilities in this species. Both parameters were assumed to depend on sex (sex) based on the difference in behaviour the sexes display during breeding. We expected males to have lower survival and encounter probabilities than females, because we assumed that their more active foraging behaviour during the breeding results in a higher risk of mortality and a lesser probability of finding their shed feathers at the nest site. We investigated annual variation in survival and encounter probabilities with time‐dependency (time) models. We also assumed that survival probability would show variation across the years not only because of annual variation, for example, in food supply or weather conditions, but also due to varying levels of poisoning in the last decade. The number of poisoning incidents was taken from the BirdCrime Database of MME BirdLife Hungary, which incorporates all detected cases by the conservation organisations of the country (Deák, Árvay, & Horváth, 2020). To investigate the relationship between survival and poisoning, we modelled survival as a function of annual poisoning rate (poison), which we calculated as follows: number of poisoned imperial eagles found in East Hungary in each interval (e.g. between 1 June 2011 and the 31 May 2012)/total number of known nesting individuals (twice the number of known nesting pairs) at the beginning of the interval, multiplied by 100 to be expressed as a percentage. Therefore, a 1% poisoning rate means that the observed number of poisoned imperial eagles equals 1% of the known nesting individuals. We presumed that the number of known nesting pairs accounts for at least 95% of the total number of nesting pairs in East Hungary (Horváth et al., 2011). We emphasise that the number of poisoned imperial eagles includes not only nesting individuals but also floaters as well as individuals of unknown age. That is because determining whether the poisoned bird was a breeder or a floater was not possible in a high number of cases. First, there was no available information on age for 34% of the poisoned birds. Second, since imperial eagles can start breeding as early as three calendar years (Horváth, 2022), poisoned birds with immature colouration (10%) could have been either floaters or breeders. Only birds in their first or second calendar year (30%) could have been safely excluded as non‐breeders. However, due to the high number of poisoned birds with unknown age or unsure status, we decided to include all poisoned individuals when calculating poisoning rates instead of discarding a significant number of detected cases. Therefore, the poisoning rate used here is only a proxy for the true poisoning rate of breeding individuals. Assuming that the observed poisoning rate of all birds is directly proportional to the poisoning rate of breeding birds, yearly changes in this observed poisoning rate reflect the yearly changes in the true poisoning rate of breeding birds. Since August 2013, a poison and carcass detection dog (PCDD) unit has also been used to help uncover poisoning incidents and proved to be more successful in both carcass and bait detection than human investigators (Deák, Árvay, & Horváth, 2020). Therefore, we corrected the number of poisoned carcasses found before August 2013 for the possibly undetected poisoning events to gain a more accurate estimation of poisoning rates in these early intervals. We did this by assuming the detection probability by the PCDD unit to be 1 and then dividing the number of poisoned imperial eagle carcasses found before August 2013 by the detection probability of human investigators. This probability was estimated to be 0.81 since out of the 42 imperial eagle carcasses detected between August 2013 and August 2020, 34 were recovered by human investigators (Deák, Árvay, & Horváth, 2020; Deák's personal communication). Lastly, we also modelled encounter probabilities as a function of sampling effort (effort), which we calculated for each year as the number of genotyped feathers (shed and chick feathers for parentage analysis included)/total number of known nesting individuals. Controlling for the number of nesting pairs in the case of both poison and effort was necessary due to the rapid expansion of the population throughout the study (Horváth, 2022).

We used the standard CJS model {ϕ(sex × time), p(sex × time)} as the general model, and the candidate model set included the variables mentioned above with both additive (+) and interaction terms (×). We also considered models with constant (.) survival or encounter probability. We constructed the sex models with males as the reference category.

Assumptions of the CJS model include the homogeneity of survival and encounter probabilities among groups of marked individuals (Amstrup et al., 2006). Therefore, we used goodness‐of‐fit tests to assess whether these conditions are met in our dataset. We investigated the goodness‐of‐fit of our general model to the capture histories with the χ2 tests of programme RELEASE v.3.0 (Burnham et al., 1987). We estimated the overdispersion parameter (c^) by dividing the overall χ2 of the component tests by the overall degrees of freedom. Then, we used this c^ estimate to adjust the AICc (Akaike Information Criterion corrected for small sample size) values of models, resulting in a QAICc‐based model selection. We also adjusted AICc values by matching the parameter counts (K) to the model structure and by considering confounded parameters, as recommended by Cooch and White (2014, Chapter 4). We ranked the models using their QAICc values and considered them equally parsimonious if the difference in their QAICc values (∆QAICc) was less than two (Burnham & Anderson, 2002). We used the Akaike weights (w i ) to evaluate the relative support of the competing models.

We used R v.4.3.1. (R Core Team, 2023) to run RMark and calculate descriptive statistics of genotype data. To create figures, we used the ggplot2 v.3.4.4. R Package (Wickham, 2016) and QGIS v.3.28.1. (QGIS Development Team, 2023). Data and code are provided as supplementary files (S1).

3 RESULTS

3.1 Individual identification

During the study, we genotyped 1730 shed feathers, which belonged to 619 breeding individuals (208 males and 411 females) with an average of 2.8 samples per individual. Out of the 15,570 genotyped loci, 1172 (7.53%) failed to amplify adequately, and we detected 106 genotyping errors (0.74%) by comparing multiple samples of the same individual. In total, we found 58 alleles at the nine microsatellite loci, ranging from four (IEAAAG09) to ten (Aa35), with an average of 6.4 alleles per locus. PI and PISIB values for the complete marker set were 9.5 × 10−9 and 5.7 × 10−4, respectively. Exclusion probabilities P1X and P2X were 0.997 and 0.963, respectively. We found deviations from the Hardy–Weinberg equilibrium (p < .05) in three loci (Aa35, Aa36 and Aa43), with the possibility of null alleles on locus Aa36 (Tables S1 and S2). As null alleles do not amplify during PCR, their presence on a locus leads to falsely identifying heterozygotes as homozygotes, which can result in false mismatches during the comparison of true parent‐offspring genotypes, that is, falsely excluding breeders as parents of a chick. That is because the chick with the null allele appears to be homozygous, while it is expected to be heterozygous based on the putative parents' genotypes. To avoid such false exclusions, we allowed for such mismatches on locus Aa36 during parentage analysis if all other loci indicated a match between the breeder and the chick.

3.2 Survival analysis

Capture histories of the 208 males and 411 females consisted of a total of 1712 presences (436 presences for males and 1276 presences for females), out of which 181 were indirect presences, that is, encounter data obtained via parentage analysis (111 and 70 for males and females, respectively). 3% of the males and 20% of the females encountered on the first occasion were also encountered on the last occasion of the 12 years study (one male and 17 females).

Based on the RELEASE goodness‐of‐fit test results, our dataset satisfied the condition of homogeneous encounter probabilities for both males and females. However, we found some violations of the assumption of homogeneous survival probabilities for females (p = .023). This could be mainly attributed to the results of females encountered in 2012, 2014 and 2017. In all three cases, individuals encountered for the first time in the given year were indicated to have lower survival probabilities than those that had been encountered before. Overall, we detected overdispersion to be low, with an estimate of c^ = 1.148, and we carried out the model selection procedure after adjusting the AICc values with this c^ estimate.

In each of the four most supported models (∆QAICc < 2, Σw i  = 0.746), the encounter probability was a function of sex × time (Table 1). Regarding survival, several different models received similar support, as detailed below.

TABLE 1 Model selection results for Cormack‐Jolly‐Seber models estimating annual apparent survival (ϕ) and encounter probability (p) as a function of sex, time and poisoning rate for breeding eastern imperial eagles in East Hungary, 2011–2022. Poisoning rate was calculated as the number of poisoned imperial eagles found/number of nesting individuals x 100.

Model	QAICc	ΔQAICc	w i	K	Deviance	
{ϕ(.), p(sex × time)}	3250.20	0.00	0.268	23	1042.56	
{ϕ(poison), p(sex × time)}	3250.90	0.70	0.188	24	1041.19	
{ϕ(sex × poison), p(sex × time)}	3250.92	0.72	0.187	26	1037.07	
{ϕ(sex), p(sex × time)}	3252.09	1.89	0.104	24	1042.39	
{ϕ(sex + poison), p(sex × time)}	3252.77	2.57	0.074	25	1041.00	
{ϕ(.), p(sex + time)}	3253.11	2.91	0.062	13	1065.97	
{ϕ(poison), p(sex + time)}	3254.15	3.96	0.037	14	1064.97	
{ϕ(sex), p(sex + time)}	3254.38	4.18	0.033	14	1065.19	
{ϕ(sex + poison), p(sex + time)}	3255.42	5.22	0.020	15	1064.19	
{ϕ(sex × poison), p(sex + time)}	3256.21	6.00	0.013	16	1062.94	
{ϕ(sex × time), p(sex × time)}	3273.28	23.07	0.000	42	1025.93	
Note: QAICc: Quasi‐Akaike Information Criterion corrected for small sample size and adjusted with c^ = 1.148, ΔQAICc: Difference in the QAICc values of the model and the best‐supported model, w i : Model weight, K: Model parameter count. Models receiving little support (w i  < 0.01) are not shown, except for the general model {ϕ(sex × time), p(sex × time)}.

The minimum QAICc model {ϕ(.), p(sex × time)} estimated survival as a constant probability over time with a 0.916 ± 0.008 SE probability of survival independent of sex. Models in which survival was a function of poison, sex × poison, and sex were equally parsimonious (∆QAICc < 2, Table 1).

The poisoning rate (number of poisoned imperial eagles found/number of nesting individuals x 100) showed two peaks during the studied period, in intervals 2011–2012 and 2018–2019 (Figure 2a). Based on model {ϕ(poison), p(sex × time)}, poisoning had a negative relationship with survival, resulting in a 0.929 ± 0.013 SE estimated survival for the interval with the lowest poisoning rate (0.46%, 2017–2018) and a 0.884 ± 0.031 SE survival probability for the interval with the highest poisoning rate (4.68%, 2011–2012) (Figure 2b, black). However, the CI for the coefficient of the poisoning effect overlapped zero (β = −0.127, 95%CI = −0.333 to +0.080).

FIGURE 2 Poisoning rates (number of poisoned eastern imperial eagles found/number of nesting eastern imperial eagles × 100) in East Hungary between 2011 and 2022. Source: BirdCrime Database of MME BirdLife Hungary (a). Estimated apparent survival probabilities (ϕ± 95% CI) of breeding eastern imperial eagles from models {ϕ(sex × poison), p(sex × time)} (blue: Male, red: Female) and {ϕ (poison), p(sex × time)} (black: All) (b).

Males were estimated to have a lower survival probability than females when survival was assumed to be sex‐dependent but constant in time (model {ϕ(sex), p(sex × time)}). However, this difference was small, and the CI for the coefficient of the sex effect overlapped zero (males: 0.909 ± 0.018 SE, females: 0.918 ± 0.009 SE, β = +0.106, 95%CI = −0.392–0.605). This model received somewhat less support (w i  = 0.104) than model {ϕ(sex × poison), p(sex × time)} (w i  = 0.187), indicating that the difference in the survival of males and females may be attributed to their survival being affected by poisoning to a different degree. Based on the estimates from this model, only male survival probability was affected significantly by poisoning (β = −0.507, 95%CI = −0.927 to −0.087) (Figure 2b, blue), implying that poisoning may have a greater effect on the survival of males compared to females.

Regarding the unconstrained time model {ϕ(sex × time), p(sex × time)}, some of the survival parameters were not estimable (characterised by large SEs or SE = 0.000) due to the high number of parameters compared to the amount of data available. Those that were estimable were mostly in agreement with the estimates from the ϕ(sex × poison) model, as the estimated female survivals showed no clear relationship with poisoning and the male estimates were lower in the 2 years with the highest poisoning rates. The only interval when the estimates of the two models were clearly different is 2017–2018, when the unconstrained model estimated low male survival despite a low poisoning rate (Table S2).

Estimated encounter probabilities were lower for males than females in all models where encounter probability was a function of sex. Estimates of encounter probabilities from the best‐supported model {ϕ(.), p(sex × time)} ranged from 0.103 ± 0.032 SE (in 2021) to 0.624 ± 0.064 SE (in 2017) for males and from 0.410 ± 0.040 SE (in 2021) to 0.765 ± 0.040 SE (in 2017) for females (Figure 3a). This result generally coincides with the percentage of nesting individuals identified each year, which were lower for males compared to females in all years and showed variation across the years, ranging from 6.50% (in 2021) to 30.0% (in 2017) for males and from 31.5% (in 2021) to 67.7% (in 2017) for females (Figure 3b). Values of yearly sampling effort (number of genotyped feathers/number of nesting individuals) showed an average of 0.379 ± 0.136 SD and ranged from 0.192 (in 2021) to 0.605 (in 2015) (Figure 3b). Models where encounter probability was constrained as a function of this effort covariate, despite having fewer parameters, received substantially less support than fully time‐dependent models (Σw time /Σw effort  = 0.987/0.007 = 141), which implies that the number of genotyped feathers was not an adequate predictor of encounter probability.

FIGURE 3 Estimates of yearly encounter probabilities (p) from model {ϕ(.), p(sex × time)} for breeding male (blue) and female (red) eastern imperial eagles in East Hungary, 2011–2022 (a). Sample sizes of the mark‐recapture analysis on breeding eastern imperial eagles in East Hungary, 2011–2022 (b): Number of nesting pairs (black), number of males identified (all presences: Dark blue; only direct presences from shed feathers: Light blue), number of females identified (all presences: Red; only direct presences from shed feathers: Orange) and number of feathers (shed and chick feathers) genotyped (grey bars).

4 DISCUSSION

We estimated annual apparent survival probabilities for breeding eastern imperial eagles in the East Hungarian population using a mark‐recapture method based on genetic identification. Although the estimated survival probabilities are only apparent by definition, as the CJS model does not separate actual survival and emigration from the sampling area (Amstrup et al., 2006), we presume that since breeding imperial eagles exhibit high territory fidelity (Rudnick et al., 2005; Vili, Szabó, et al., 2013), the estimated apparent survivals are good proxies to the actual survival probabilities. The previously reported high territory fidelity of breeding imperial eagles is also supported by our data, as over the 12 years of our study, only 9 males (4% of all males studied) and 44 females (9% of all females studied) were identified from more than one territory. Such encounters correspond to less than 4% of all the 1712 presences detected and, in most cases, birds moved only a short distance to the neighbouring territory, similarly to the observations of Vili, Szabó, et al. (2013).

The East Hungarian population we studied is the largest in the Western range of the imperial eagle's distribution, and therefore, investigating its demography plays a crucial role in preserving this species in Europe. Adult survival is often considered the most important demographic parameter for long‐lived species (Katzner et al., 2006; Sæther & Bakke, 2000; Stahl & Oli, 2006); thus, the results of this study are essential for a future viability analysis for this population (Lacy, 2019). Our estimates refer to the period 2011–2022, which was characterised by variable levels of poisoning and the intensive conservation of the species in the context of the Helicon LIFE (2012–2016, LIFE10NAT/HU/019) and PannonEagle LIFE (2017–2023, LIFE15 NAT/HU/000902) projects. We estimated >90% annual survival on average, a typical estimate for a large raptor such as the imperial eagle (Newton et al., 2016). However, our estimates were lower for the first years of the study when poisoning rates were higher. These results agreed with the pattern in the population growth curve of imperial eagles during this period: a plateau was detected in the number of nesting pairs between 2011 and 2014 instead of the exponential growth observed before and after this period (Horváth, 2022, also shown in Figure 3b), implying that the high poisoning rates in 2011–2013 had major effects on population growth. Moreover, our data also provided moderate evidence for the sex‐dependent effect of poisoning on adult survival. The estimated negative effect of poisoning on survival was stronger in males than in females. Therefore, poisoning may lead to a female‐biassed adult sex ratio (ASR) in this population. This result implies that due to the strict monogamy in this species (Rudnick et al., 2005), the number of males may be a limiting factor to the growth of this population. A female‐biassed ASR can also affect the evolution of reproductive traits like age‐to‐maturation (Ancona et al., 2020) and may result in increased competition and therefore elevated mortality in females, which in turn could have adverse effects on population viability (Kokko & Jennions, 2008; Székely et al., 2014). To our knowledge, no direct evidence of sex difference in age‐to‐maturation is available in the eastern imperial eagle, but observations of a higher number of immature males than females in the breeding populations suggest that males may reach maturity at an earlier age than females (Demerdzhiev, Gradev, et al., 2011; Demerdzhiev, Stoychev, et al., 2011, Horváth's personal communication). Similar observations were also made for the Spanish imperial eagle (Margalida et al., 2008; Margalida, González, Sánchez, Oria, & Prada, 2007). An earlier maturation of males could mitigate the negative effect of male‐biassed mortality on population growth. However, it is important to note that a higher proportion of immature males in the breeding population may not necessarily mean that males reach maturity sooner than females but could also be the outcome of male‐biassed adult mortality.

The annual survival of breeding imperial eagles was estimated to be 91.6% on average. This survival estimate is close to the one estimated for the Bulgarian population based on an integrated population model (92.4%, Demerdzhiev et al., 2015). However, previous turnover‐based survival estimates for the imperial eagle were considerably lower: 84% for the population in Kazakhstan (Rudnick et al., 2005) and a maximum of 72.3% annual survival rate for females in this East Hungarian population between 1997 and 2006, with no estimate for males (Vili, Szabó, et al., 2013). The 72.3% estimate of Vili, Szabó, et al. (2013) is exceptionally low even compared to estimates obtained for the imperial eagle's sister species, the Spanish imperial eagle (92.7%–95.3%, Ferrer & Calderón, 1990; 91,5%, Ortega et al., 2009) and for other Aquila species (Bonelli's eagle Aquila fasciata: 83.9%–96.1%, Real & Mañosa, 1997; 87.0%, Hernández‐Matías et al., 2011; golden eagle Aquila chrysaetos: 91.0%, Whitfield et al., 2004; 90.0%, Millsap et al., 2022; 90.5% Hunt et al., 2017; 93.0%, Crandall et al., 2019).

In the following, we discuss whether poisoning could explain the previously estimated exceptionally low breeder survival rate for the Hungarian population or if some other factors might be responsible. In Hungary, the first two carcasses of illegally poisoned imperial eagles were recovered in 2005, and poisoning has been considered the leading anthropogenic mortality cause of imperial eagles ever since (Horváth et al., 2011). However, we suppose that poisoning was also present in earlier years, as suggested by the two reported cases of alleged poisoning in 1980–2000 (Horváth et al., 2011). Only one poisoning rate estimate is available for the time range of the study of Vili, Szabó, et al. (2013). This poisoning rate calculated for 2005–2006 is only 2.74%, which is very low for the estimated female survival rate of 72.3%. Our models predicted much higher survival rates for this relatively low poisoning rate: a 91.6% annual survival for females from the ϕ(sex × poison) model and a 90.7% annual survival on average for both sexes from the ϕ(poison) model.

Reasons behind these seemingly contradicting results can include (1) actual poisoning rates may have been higher than detected for the previous period, with a more substantial effect on the survival of females than the one reported here; (2) other factors than poisoning may have contributed to this very low survival rate; or (3) methodological differences between the two studies.

First, many poisoning incidents could have remained uncovered between 1997 and 2006 since no organised poison searches were conducted at the time, and carcasses were rarely tested for poison compounds (Deák's personal communication). Therefore, actual poisoning rates may have been much higher than detected. This is also implied by the high poisoning rates observed in the following years when the monitoring of poisoning was already in focus after the recovery of two poisoned imperial eagles in 2005 (Horváth et al., 2011). In addition, the survival of females may react to poisoning detectably only above a specific rate. One plausible explanation is that the frequent turnover of males at higher levels of poisoning could introduce a striking effect on the survival of females. This hypothesis is based on studies conducted on other monogamous bird species, which revealed that mate change can lead to a lower apparent survival (Culina et al., 2015; Jankowiak et al., 2018), even in long‐lived species (Leach et al., 2020).

Second, factors other than poisoning may have also caused the previously low survival of females. The frequency of electrocution, the second most important mortality factor of imperial eagles, remained more or less constant over the years (Deák, Fatér, et al., 2020) and thus cannot explain the differences in survival. Another possible factor could be the change in habitat. During the 1970s, imperial eagles in Hungary resided almost exclusively in the North Hungarian Mountains (Bagyura et al., 2002). These lower‐quality mountainous habitats served as a refuge for the species since lower human population density and closed vegetation in these areas ensured a lower risk of persecution (Horváth et al., 2011). However, the population has expanded significantly over the decades (Horváth, 2022), and a larger proportion now resides in lowland areas. These lowland areas support higher fledging success due to better foraging opportunities than the previously occupied mountainous areas (Horváth et al., 2014) and may also promote higher survival in breeders.

Lastly, methodological differences may have also contributed to the difference between our estimate and that of Vili, Szabó, et al. (2013). The latter was based on turnover rate (annual % of breeding individuals replaced by new breeding individuals), which does not account for imperfect detectability, while mark‐recapture does. Consequently, breeding dispersal could have introduced a negative bias in the survival estimate of Vili, Szabó, et al. (2013). However, considering the high territory fidelity exhibited by this species (Rudnick et al., 2005; Vili, Szabó, et al., 2013), dispersal only is unlikely to have caused the large difference in the estimates. Additionally, the estimate of Vili, Szabó, et al. (2013) was based on four times fewer individuals than our study and, therefore, carries more uncertainty.

Since Vili, Szabó, et al. (2013) only gave estimates for females, our study is the first to explore sex differences in survival in this population of imperial eagles. Models featuring the sex‐dependence of survival indicated lower male survival compared to females mainly attributable to the lower survival probabilities of males in years of high poisoning, especially in the first interval. There was only one interval where male survival, in contrast to the low poisoning rate, was estimated to be low by the unconstrained time model. This contradiction may be explained by the imperfect detection of poisoning events or by the increase of another mortality factor in this interval.

A possible mortality factor that may show yearly fluctuations similarly to the poisoning rate and could, therefore, interact with the effect of poisoning on breeder survival, is food supply. Due to their body size dimorphism, food shortages may affect males and females differently (Payevsky, 2016). Therefore, prey availability may be enhancing or counteracting the possible sex‐dependent effect of poisoning. The latter might explain why models with the sex‐dependent effect of poisoning were not better supported.

Another reason why the model selection may not have strongly supported the relationship between poisoning and survival is that in most years of the study, only a relatively low level of poisoning was observed. The negative effect of poisoning on survival may only be prominent at higher rates of poisoning. This is also suggested by the estimated survival probabilities (Figure 2b), as in years with low poisoning rates, the CIs of male and female estimates largely overlap, and the sex difference in survival only appears to be relevant in the first intervals, when poisoning rates were much higher.

In contrast to our results of lower male survival, a previous study on imperial eagles implied higher survival for males in Kazakhstan (Rudnick et al., 2005) and in the case of the Spanish imperial eagle, both female‐biassed mortality (Ferrer & Hiraldo, 1992) and no difference in the survival of the sexes have been reported (Ortega et al., 2009). In raptors, sex‐dependent survival can be expected because of behavioural dimorphism (Newton et al., 2016), through which different environmental factors, including human persecution, influence the sexes differently (Ewing et al., 2023; Hernández‐Matías et al., 2011). Our results suggest this may also be true for the imperial eagle. In spring, females spend most of their time at the nest incubating and males are the ones foraging for food (Dobrev, 2009; Horváth, 2022). This period of primarily male foraging coincides with the reported annual peak of poisoning activity from February to April (Deák, Árvay, & Horváth, 2020). This was also true in the first interval of our study, when the difference between male and female survival was estimated to be the highest and 64% of the poisoned imperial eagles were recovered between the end of February and the end of April. We hypothesise that this temporal pattern of poisoning activity, coupled with the mentioned behavioural dimorphism, could explain why males may be more at risk from poisoning than females. Moreover, although males could bring back poisoned food to the nest for the females, their quick death after exposure to the poison might prevent them from doing so. During our study, the most commonly used pesticide for the intentional killing of raptors in Hungary was carbofuran (Deák, Árvay, & Horváth, 2020), which can lead to death in minutes following its ingestion (Lehel et al., 2010; López‐Bao & Mateo‐Tomás, 2022). Therefore, foraging males have only a slight chance of bringing poisoned baits back to the nest as they are very likely to die on site, explaining why they may fall victim to poisoning more frequently than their mates. To further investigate whether breeding males are more threatened by poisoning than females, assessing the age and sex of poisoned imperial eagles would be required.

This previously mentioned behavioural dimorphism also resulted in a smaller sample size and smaller encounter probabilities for males compared to females, since the shed feathers of males could be recovered with substantially less success around the nest. Due to this, our survival estimates for males are less precise than for females. This was especially true for 2011–2012, which was not only the first interval of the study but also had a much higher poisoning rate than the average. Both factors contributed to the fact that less data were available for the estimates for this period. Thus, the uncertainty of these estimates was greater than that of later years.

To obtain more encounter data, especially on males, we also used parentage analysis to infer the presence of breeders. Supplementing our dataset with these indirect presences was made possible by the high territory and mate fidelity of the species. However, this method could have introduced heterogeneity to encounter probabilities since indirect presences could only be obtained of those birds that had been identified previously at least once and whose mate was also identified in the year in question. Additionally, indirect presences could have also introduced upward bias to our survival estimates. This would be the case if the probability of obtaining an indirect presence of an individual was positively correlated with its survival. This can happen in two ways: (i) since indirect presences can only be obtained of birds that are identified in at least one other year, there is a higher chance to obtain indirect presences of individuals that live longer; and (ii) since indirect presences are determined through parentage analysis, more data could be obtained on those birds that successfully produced offspring over several years, meaning that if breeding success is positively correlated with survival probability (see for example Espie et al., 2004), the probability of obtaining an indirect presence was higher for birds with better survival. In both cases, birds with better survival would be overrepresented in our dataset. This would be especially true for males, considering that about 25% of the male encounter data were indirect presences. Nevertheless, we found no evidence of indirect presences violating the assumptions of the CJS model based on the goodness‐of‐fit tests: overdispersion was estimated to be of an acceptable value (c^ = 1.148), which was even smaller than the one estimated for the dataset without indirect presences (c^ = 1.200). Thus, indirect presences could serve as a valuable tool in cases when encounter data are otherwise scarce, a situation typical for raptors (Newton et al., 2016).

RELEASE goodness‐of‐fit tests revealed a slight violation of the assumption of homogeneous survival probabilities in the case of females, both with and without indirect presences: females identified for the first time in 2012, 2014 and 2017 were less likely to be encountered later than females that have been encountered before. Assuming that birds that are identified for the first time in a specific year are mainly young (3–4 calendar years) birds entering the breeding stage, an explanation for this pattern could be the lower apparent survival of younger breeders. Reasons behind this could include a lower actual survival for younger individuals or a lower territory fidelity, leading to their dispersal to another breeding site outside of our study area after unsuccessful breeding attempts (Catlin et al., 2005; Jiménez‐Franco et al., 2013). There are two options to investigate if the survival of younger breeders is lower compared to adults. The first is to apply a time‐since‐marking model, but it requires that all the birds are young when identified for the first time, which is not a valid assumption in this case. The second option is the age‐dependent model; however, we know the age of only a few birds in this study. Continued observation of breeding individuals and incorporation of ringing and satellite‐tracking data of young breeders could facilitate answering this question in the future.

Our study provided the first estimates of sex‐dependent survival for imperial eagles in Hungary, a population of high conservation value due to its large size and peripheric placement in the distribution of the species. We also investigated the relationship between survival and poisoning rates in imperial eagles for the first time. In conclusion, the annual survival probability of breeding imperial eagles in East Hungary is high despite current anthropogenic threats in the country. The 26% increase in the survival estimate compared to the period of 1997–2006 implies that the conservation measures have contributed to the increase in survival of breeding birds. However, we cannot exclude methodological differences as explanations for the difference between previous and current estimates. Furthermore, while we found no definite evidence for the effect of poisoning and sex on adult survival, our results alert that male survival might be lower and more sensitive to poisoning than female survival. If that is the case, it would be important from both a conservational and evolutionary perspective, as sex‐biassed adult mortality is considered the main driver of biassed adult sex ratios in birds, which in turn influence parental roles and mating competition and can affect population viability (Ancona et al., 2020; Kokko & Jennions, 2008; Székely et al., 2014). Furthermore, in a species, such as the imperial eagle, where food provisioning during breeding is mostly dependent on the male and experienced adults have better hunting capacity (Margalida et al., 2008), male‐biassed mortality can also have a direct impact on productivity. While the scarcity of male data makes estimates of male survival less precise than that of females, our method of indirectly assessing the presence of breeders through parentage analysis greatly contributed to the data on males. To obtain a deeper understanding of the demography of the population, more precise estimates of breeding male survival would be beneficial in the future. Furthermore, the survival rates of the floater age groups and the young breeders are unknown; a future task is to estimate these from ringing and satellite‐tracking data. Moreover, incorporating these results into a population viability analysis could further facilitate our understanding of the effect of poisoning on survival and its consequences on population growth in this protected, long‐lived raptor.

AUTHOR CONTRIBUTIONS

Bernadett Zsinka: Conceptualization (equal); data curation (equal); formal analysis (equal); investigation (equal); methodology (equal); visualization (lead); writing – original draft (lead). Szilvia Pásztory‐Kovács: Conceptualization (equal); data curation (equal); formal analysis (equal); investigation (equal); methodology (equal); project administration (equal); resources (equal); supervision (equal); visualization (supporting); writing – original draft (supporting); writing – review and editing (equal). Szilvia Kövér: Conceptualization (equal); formal analysis (equal); methodology (equal); supervision (equal); visualization (supporting); writing – original draft (supporting); writing – review and editing (equal). Nóra Vili: Investigation (supporting); writing – review and editing (supporting). Márton Horváth: Conceptualization (equal); data curation (equal); funding acquisition (lead); investigation (equal); project administration (equal); resources (equal); supervision (equal); writing – review and editing (supporting).

CONFLICT OF INTEREST STATEMENT

The authors declare that they have no conflict of interest.

Supporting information

Data S1.

Tables S1–S2.

ACKNOWLEDGEMENTS

Population monitoring and sample collection were conducted by the members of the Hungarian Imperial Eagle Working Group, which incorporates employees and volunteers of MME BirdLife Hungary and the ranger service of the relevant national park directories. We are particularly grateful to Márton Árvay, Attila Bereczky, Imre Fatér, Tibor Juhász, Péter Őze, János Sasvári and Nándor Seres for their contribution to sample collection. We are also grateful to Gábor Deák for his help in fieldwork and for providing data on poisoning incidents. We thank Krisztián Szabó, Péter Kacz, Szilvia Jakab and Ádám Gyarmati for their help in laboratory work. We thank two anonymous reviewers for their valuable suggestions that improved the manuscript. The field collection of samples and laboratory analyses were funded partially by the European Commission's LIFE‐Nature programme in the frame of Helicon LIFE (2012‐2016, LIFE10NAT/HU/019) and PannonEagle LIFE (2017‐2023, LIFE15 NAT/HU/000902) projects. Bernadett Zsinka received a PhD fellowship from the Doctoral School of Veterinary Science, University of Veterinary Medicine Budapest, Hungary. Feather and DNA samples from birds were collected under the research permission of MME BirdLife Hungary: 14/2583‐5/2011 (2011‐2015) and OKTF‐KP/56‐25/2015 (2015‐2022) issued by National Environmental and Nature Conservation Inspectorate, and its modifications PE‐KTFO/1867‐9/2018, PE‐KTFO/1867‐10/2018, PE‐KTFO/1867‐11/2018 and PE‐KTFO/507‐16/2019 issued by the Government Office of Pest County.

DATA AVAILABILITY STATEMENT

Data and code are provided as supplementary files S1.
==== Refs
REFERENCES

Amstrup, S. C. , McDonald, T. L. , & Manly, B. F. J. (Eds.). (2006). Handbook of capture‐recapture analysis. Princeton University Press. 10.1515/9781400837717
Ancona, S. , Liker, A. , Carmona‐Isunza, M. C. , & Székely, T. (2020). Sex differences in age‐to‐maturation relate to sexual selection and adult sex ratios in birds. Evolution Letters, 4 (1 ), 44–53. 10.1002/evl3.156 32055410
Bagyura, J. , Szitta, T. , Haraszthy, L. , Firmánszky, G. , Viszló, L. , Kovács, A. , Demeter, I. , & Horváth, M. (2002). Population increase of Imperial eagle (Aquila heliaca) in Hungary between 1980 and 2000. Aquila, 107–108 , 133–144.
Begon, M. , Mortimer, M. , & Thompson, D. J. (1996). A unified study of animals and plants (3rd ed.). Blackwell Science Ltd.
Beja‐Pereria, A. , Oliveira, R. , Alves, P. C. , Schwartz, M. K. , & Luikart, G. (2009). Advancing ecological understandings through technological transformations in noninvasive genetics. Molecular Ecology Resources, 9 (5 ), 1279–1301. 10.1111/j.1755-0998.2009.02699.x 21564900
BirdLife International . (2019). Aquila heliaca (amended version of 2017 assessment). The IUCN Red List of Threatened Species 2019: e.T22696048A155464885. Accessed on 16 March 2024.
Buechley, E. R. , Santangeli, A. , Girardello, M. , Neate‐Clegg, M. H. C. , Oleyar, D. , McClure, C. J. W. , & Şekercioğlu, Ç. H. (2019). Global raptor research and conservation priorities: Tropical raptors fall prey to knowledge gaps. Diversity and Distributions, 25 (6 ), 856–869. 10.1111/ddi.12901
Burnham, K. P. , & Anderson, D. R. (2002). Model selection and multimodel inference: A practical information‐theoretic approach. Springer.
Burnham, K. P. , Anderson, D. R. , White, G. C. , Brownie, C. , & Pollock, K. H. (1987). Design and analysis methods for fish survival experiments based on Release‐recapture. American Fisheries Society.
Busch, J. D. , Katzner, T. E. , Bragin, E. , & Keim, P. (2005). Tetranucleotide microsatellites for aquila and haliaeetus eagles. Molecular Ecology Notes, 5 (1 ), 39–41. 10.1111/j.1471-8286.2004.00823.x
Bush, K. L. , Vinsky, M. D. , Aldridge, C. L. , & Paszkowski, C. A. (2005). A comparison of sample types varying in invasiveness for use in DNA sex determination in an endangered population of greater sage‐grouse (Centrocercus uropihasianus). Conservation Genetics, 6 (5 ), 867–870. 10.1007/s10592-005-9040-6
Catlin, D. H. , Rosenberg, D. K. , & Haley, K. L. (2005). The effects of nesting success and mate fidelity on breeding dispersal in burrowing owls. Canadian Journal of Zoology, 83 (12 ), 1574–1580. 10.1139/z05-154
Connallon, T. , Beasley, I. J. , McDonough, Y. , & Ruzicka, F. (2022). How much does the unguarded X contribute to sex differences in life span? Evolution Letters, 6 (4 ), 319–329. 10.1002/evl3.292 35937469
Cooch, E. G. , & White, G. C. (2014). Program MARK: A gentle introduction (13th ed.). Cornell University. Available online at. http://www.phidot.org/software/mark/docs/book/ Accessed on 17 October 2023.
Crandall, R. H. , Craighead, D. J. , Bedrosian, B. , & Slabe, V. A. (2019). Survival estimates and cause of mortality of Golden eagles in south‐Central Montana. Journal of Raptor Research, 53 (1 ), 38–45. 10.3356/JRR-18-22
Culina, A. , Lachish, S. , & Sheldon, B. C. (2015). Evidence of a link between survival and pair fidelity across multiple tit populations. Journal of Avian Biology, 46 (5 ), 507–515. 10.1111/jav.00661
Deák, G. , Árvay, M. , & Horváth, M. (2020). Using detection dogs to reveal illegal pesticide poisoning of raptors in Hungary. Journal of Vertebrate Biology, 69 (3 ), 20110. 10.25225/jvb.20110
Deák, G. , Fatér, I. , Juhász, T. , & Horváth, M. (2020). Parlagi sasok (Aquila heliaca) pusztulási és kézre kerülési okainak alakulása Magyarországon 2010 és 2019 között (Causes of the death, injuries and diseases of Eastern Imperial Eagles (Aquila heliaca) in Hungary between 2010 and 2019‐in Hungarian with E). Heliaca‐the Annual of the Raptor Conservation Group of MME BirdLife Hungary, 16 , 114–117.
Demerdzhiev, D. , Dobrev, D. , Popgeorgiev, G. , & Stoychev, S. (2022). Landscape alteration affects the demography of an endangered avian predator by reducing the habitat quality. Avian Research, 13 , 100030. 10.1016/j.avrs.2022.100030
Demerdzhiev, D. , Horváth, M. , Kovács, A. , Stoychev, S. , & Karyakin, I. (2011). Status and population trend of the eastern imperial eagle (Aquila heliaca) in Europe in the period 2000‐2010. Acta Zoologica Bulgarica, 63 (SUPPL.3 ), 5–14.
Demerdzhiev, D. , Stoychev, S. , Dobrev, D. , Spasov, S. , & Oppel, S. (2015). Studying the demographic drivers of an increasing Imperial eagle population to inform conservation management. Biodiversity and Conservation, 24 (3 ), 627–639. 10.1007/s10531-014-0841-0
Demerdzhiev, D. A. , Gradev, G. Z. , Stoychev, S. A. , Ivanov, I. I. , Petrov, T. H. , & Marin, S. A. (2011). Increase of the population of the eastern Imperial eagle (Aquila heliaca) in Bulgaria. Acta Zoologica Bulgarica, 63 (SUPPL.3 ), 41–54.
Demerdzhiev, D. A. , Stoychev, S. A. , Terziev, N. G. , & Angelov, I. D. (2011). Status of the eastern Imperial eagle (Aquila heliaca) in the European part of Turkey. Acta Zoologica Bulgarica, 63 (SUPPL.3 ), 87–93.
Dobrev, D. D. (2009). Some aspects of the nest etology of the eastern Imperial eagle (Aquila heliaca) (Aves: Accipitridae) in Sakar Mountain. Ecologia Balkanica, 1 , 59–66.
Dobson, A. P. (1987). A comparison of seasonal and annual mortality for both sexes of fifteen species of common British birds. Ornis Scandinavica, 18 (2 ), 122–128. 10.2307/3676846
Espie, R. H. M. , James, P. C. , Oliphant, L. W. , Warkentin, I. G. , & Lieske, D. J. (2004). Influence of nest‐site and individual quality on breeding performance in Merlins Falco columbarius. Ibis, 146 (4 ), 623–631. 10.1111/j.1474-919x.2004.00294.x
Ewing, S. R. , Thomas, C. E. , Butcher, N. , Denman, B. , Douglas, D. J. T. , Anderson, D. I. K. , Anderson, G. Q. A. , Bray, J. , Downing, S. , Dugan, R. , Etheridge, B. , Hayward, W. , Howie, F. , Roos, S. , Thomas, M. , Weston, J. , Smart, J. , & Wilson, J. D. (2023). Illegal killing associated with gamebird management accounts for up to three‐quarters of annual mortality in hen harriers Circus cyaneus. Biological Conservation, 283 , 110072. 10.1016/j.biocon.2023.110072
Ferrer, M. , & Calderón, J. (1990). The Spanish imperial eagle aquila adalberti C. L. Brehm 1861 in Doñana National Park (south West Spain): A study of population dynamics. Biological Conservation, 51 (2 ), 151–161. 10.1016/0006-3207(90)90109-3
Ferrer, M. , & Hiraldo, F. (1992). Man‐induced sex‐biased mortality in the Spanish imperial eagle. Biological Conservation, 60 (1 ), 57–60. 10.1016/0006-3207(92)90799-S
Gil‐Sánchez, J. M. , Bautista, J. , Godinho, R. , & Moleón, M. (2021). Detection of individual replacements in a long‐lived bird species, the Bonelli's eagle (Aquila fasciata), using three noninvasive methods. Journal of Raptor Research, 55 (4 ), 552–564. 10.3356/JRR-20-53
Goor, R. M. , Neall, L. F. , Hoffman, D. , & Sherry, S. T. (2011). A mathematical approach to the analysis of multiplex DNA profiles. Bulletin of Mathematical Biology, 73 (8 ), 1909–1931. 10.1007/s11538-010-9598-0 21103945
Hailer, F. , Gautschi, B. , & Helander, B. (2005). Development and multiplex PCR amplification of novel microsatellite markers in the White‐tailed Sea eagle, Haliaeetus albicilla (Aves: Falconiformes, Accipitridae). Molecular Ecology Notes, 5 (4 ), 938–940. 10.1111/j.1471-8286.2005.01122.x
Hernández, M. , & Margalida, A. (2009). Poison‐related mortality effects in the endangered Egyptian vulture (Neophron percnopterus) population in Spain. European Journal of Wildlife Research, 55 (4 ), 415–423. 10.1007/s10344-009-0255-6
Hernández‐Matías, A. , Real, J. , Pradel, R. , Ravayrol, A. , & Vincent‐Martin, N. (2011). Effects of age, territoriality and breeding on survival of Bonelli's eagle Aquila fasciata. Ibis, 153 (4 ), 846–857. 10.1111/j.1474-919X.2011.01158.x
Heuck, C. , Herrmann, C. , Wendt, J. , Krone, O. , Brandl, R. , & Albrecht, J. (2020). Sex‐ but not age‐biased wind turbine collision mortality in the White‐tailed eagle Haliaeetus albicilla. Journal of Ornithology, 161 (3 ), 753–757. 10.1007/s10336-020-01757-w
Horváth, M. (2022). Parlagi sas Aquila heliaca Savigny, 1823 (Imperial Eagle Aquila heliaca Savigny, 1823) [in Hungarian]. In L. Haraszthy  & J. Bagyura (Eds.), Magyarország ragadadozómadarai és baglyai. 1. kötet Vágómadár‐alakúak. (Raptors and owls of Hungary. Volume 1. Accipitriformes) (pp. 310–357). Magyar Madártani és Természetvédelmi Egyesület.
Horváth, M. , Demeter, I. , Fatér, I. , Firmánszky, G. , Kleszó, A. , Kovács, A. , Szitta, T. , Tóth, I. , Zalai, T. , & Bagyura, J. (2011). Population dynamics of the eastern Imperial eagle (Aquila heliaca) in Hungary between 2001 and 2009. Acta Zoologica Bulgarica, 63 (SUPPL.3 ), 61–70.
Horváth, M. , Fatér, I. , Juhász, T. , Deák, G. , & Bereczky, A. (2020). A parlagi sas (Aquila heliaca) fészkelőállományának alakulása Magyarországon 2017 és 2019 között (the nesting population of eastern Imperial eagles (Aquila heliaca) in Hungary between 2017 and 2019‐ in Hungarian with English summary). Heliaca‐the Annual of the Raptor Conservation Group of MME BirdLife Hungary, 16 , 42–45.
Horváth, M. , Haraszthy, L. , Bagyura, J. , & Kovács, A. (2002). Eastern Imperial eagle (Aquila heliaca) populations in Europe. Aquila, 107–108 , 193–204.
Horváth, M. , Martínez‐Cruz, B. , Negro, J. J. , Kalmár, L. , & Godoy, J. A. (2005). An overlooked DNA source for non‐invasive genetic analysis in birds. Journal of Avian Biology, 36 (1 ), 84–88. 10.1111/j.0908-8857.2005.03370.x
Horváth, M. , Szitta, T. , Bagyura, J. , Fatér, I. , Firmánszky, G. , & Moskát, C. (2014). Simultaneous effect of habitat and age on reproductive success of imperial eagles (Aquila heliaca) in Hungary. Ornis Hungarica, 22 (1 ), 57–68. 10.2478/orhu-2014-0011
Hunt, W. G. , Wiens, J. D. , Law, P. R. , Fuller, M. R. , Hunt, T. L. , Driscoll, D. E. , & Jackman, R. E. (2017). Quantifying the demographic cost of humanrelated mortality to a raptor population. PLoS One, 12 (2 ), 1–22. 10.1371/journal.pone.0172232
Jankowiak, Ł. , Cholewa, M. , & Wysocki, D. (2018). Survival costs of within‐ and between‐season mate change in the European blackbird Turdus merula. Journal of Avian Biology, 49 (6 ), 1–9. 10.1111/jav.01643
Jiménez‐Franco, M. V. , Martínez, J. E. , Pagán, I. , & Calvo, J. F. (2013). Factors determining territory fidelity in a migratory forest raptor, the booted eagle Hieraaetus pennatus. Journal of Ornithology, 154 (1 ), 311–318. 10.1007/s10336-012-0895-2
Karyakin, I. V. (2020). Breeding population structure of the eastern Imperial eagle. Raptors Conservation, 41 , 64–332. 10.19074/1814-8654-2020-41-64-332
Katzner, T. E. , Bragin, E. A. , & Milner‐Gulland, E. J. (2006). Modelling populations of long‐lived birds of prey for conservation: A study of imperial eagles (Aquila heliaca) in Kazakhstan. Biological Conservation, 132 (3 ), 322–335. 10.1016/j.biocon.2006.04.024
Kenward, R. , Katzner, T. , Wink, M. , Marcström, V. , Walls, S. , Karlbom, M. , Pfeffer, R. , Bragin, E. , Hodder, K. , & Levin, A. (2007). Rapid sustainability modeling for raptors by Radiotagging and DNA‐fingerprinting. The Journal of Wildlife Management, 71 (1 ), 238–245. 10.2193/2005-471
Kitowski, I. , Łopucki, R. , Stachniuk, A. , & Fornal, E. (2020). A pesticide banned in the European Union over a decade ago is still present in raptors in Poland. Environmental Conservation, 47 (4 ), 310–314. 10.1017/S037689292000034X
Köhler, H. R. , & Triebskorn, R. (2013). Wildlife ecotoxicology of pesticides: Can we track effects to the population level and beyond? Science, 341 (6147 ), 759–765. 10.1126/science.1237591 23950533
Kokko, H. , & Jennions, M. D. (2008). Parental investment, sexual selection and sex ratios. Journal of Evolutionary Biology, 21 (4 ), 919–948. 10.1111/j.1420-9101.2008.01540.x 18462318
Kylmänen, A. , Karabanina, E. , Ollila, T. , Ponnikas, S. , & Kvist, L. (2023). Turnover and natal dispersal in the finnish golden eagle (aquila chrysaetos) population. Diversity, 15 (4 ), 567. 10.3390/d15040567
Laake, J. L. (2013). RMark: An R Interface for analysis of capture‐recapture data with MARK. https://CRAN.R‐project.org/package=RMark
Lacy, R. C. (2019). Lessons from 30 years of population viability analysis of wildlife populations. Zoo Biology, 38 (1 ), 67–77. 10.1002/zoo.21468 30585658
Lazarova, I. , Dobrev, D. , Gradev, G. , Petrov, R. , Stoychev, S. , Klisurov, I. , & Demerdzhiev, D. (2020). Main mortality factors for the eastern Imperial eagle (Aquila heliaca Savigny, 1809) in Bulgaria. Ornis Hungarica, 28 (2 ), 120–134. 10.2478/orhu-2020-0021
Leach, A. G. , Riecke, T. V. , Sedinger, J. S. , Ward, D. H. , & Boyd, S. (2020). Mate fidelity improves survival and breeding propensity of a long‐lived bird. Journal of Animal Ecology, 89 (10 ), 2290–2299. 10.1111/1365-2656.13286 32596854
Lebreton, J. D. , Burnham, K. P. , Clobert, J. , & Anderson, D. R. (1992). Modeling survival and testing biological hypotheses using marked animals: A unified approach with case studies. Ecological Monographs, 62 (1 ), 67–118. 10.2307/2937171
Lehel, J. , Laczay, P. , Déri, J. , Darin, E. G. , & Budai, P. (2010). Model study on the clinical signs and residue concentrations of sublethal carbofuran poisoning in birds. Journal of Wildlife Diseases, 46 (4 ), 1274–1278. 10.7589/0090-3558-46.4.1274 20966279
López‐Bao, J. V. , & Mateo‐Tomás, P. (2022). Wipe out highly hazardous pesticides to deter wildlife poisoning: The case of carbofuran and aldicarb. Biological Conservation, 275 , 109747. 10.1016/j.biocon.2022.109747
López‐López, P. , Ferrer, M. , Madero, A. , Casado, E. , & McGrady, M. (2011). Solving man‐induced large‐scale conservation problems: The Spanish Imperial eagle and power lines. PLoS One, 6 (3 ), e17196. 10.1371/journal.pone.0017196 21399692
Maklakov, A. A. , & Lummaa, V. (2013). Evolution of sex differences in lifespan and aging: Causes and constraints. BioEssays, 35 (8 ), 717–724. 10.1002/bies.201300021 23733656
Margalida, A. , González, L. M. , Sánchez, R. , Oria, J. , & Prada, L. (2007). A long‐term large‐scale study of the breeding biology of the Spanish imperial eagle (Aquila adalberti).. Journal of Ornithology, 148 (3 ), 309–322. 10.1080/00063650709461462
Margalida, A. , González, L. M. , Sánchez, R. , Oria, J. , Prada, L. , Caldera, J. , Aranda, A. , & Molina, J. I. (2007). A long‐term large‐scale study of the breeding biology of the Spanish imperial eagle (Aquila adalberti). Journal of Ornithology, 148 (3 ), 309–322. 10.1007/s10336-007-0133-5
Margalida, A. , Mañosa, S. , González, L. M. , Ortega, E. , Sánchez, R. , & Oria, J. (2008). Breeding of non‐adults and effects of age on productivity in the Spanish Imperial eagle Aquila adalberti. Ardea, 96 (2 ), 173–180. 10.5253/078.096.0203
Martínez‐Cruz, B. , David, V. A. , Godoy, J. A. , Negro, J. J. , O'Brien, S. J. , & Johnson, W. E. (2002). Eighteen polymorphic microsatellite markers for the highly endangered Spanish imperial eagle (Aquila adalberti) and related species. Molecular Ecology Notes, 2 (3 ), 323–326. 10.1046/j.1471-8286.2002.00231.x
Mateo‐Tomás, P. , Olea, P. P. , Mínguez, E. , Mateo, R. , & Viñuela, J. (2020). Direct evidence of poison‐driven widespread population decline in a wild vertebrate. Proceedings of the National Academy of Sciences of the United States of America, 117 (28 ), 16418–16423. 10.1073/pnas.1922355117 32601195
McClure, C. J. W. , Vargas, F. H. , Amar, A. , Concepcion, C. B. , MacColl, C. , & Sumasgutner, P. (2022). Conservation letter: Monitoring raptor populations‐a call for increased global collaboration and survey standardization. Journal of Raptor Research, 57 (1 ), 106–113. 10.3356/JRR-22-68
McClure, C. J. W. , Westrip, J. R. S. , Johnson, J. A. , Schulwitz, S. E. , Virani, M. Z. , Davies, R. , Symes, A. , Wheatley, H. , Thorstrom, R. , Amar, A. , Buij, R. , Jones, V. R. , Williams, N. P. , Buechley, E. R. , & Butchart, S. H. M. (2018). State of the world's raptors: Distributions, threats, and conservation recommendations. Biological Conservation, 227 , 390–402. 10.1016/j.biocon.2018.08.012
Millsap, B. A. , Zimmerman, G. S. , Kendall, W. L. , Barnes, J. G. , Braham, M. A. , Bedrosian, B. E. , Bell, D. A. , Bloom, P. H. , Crandall, R. H. , Domenech, R. , Driscoll, D. , Duerr, A. E. , Gerhardt, R. , Gibbs, S. E. J. , Harmata, A. R. , Jacobson, K. , Katzner, T. E. , Knight, R. N. , Lockhart, J. M. , … Watson, J. W. (2022). Age‐specific survival rates, causes of death, and allowable take of golden eagles in the western United States. Ecological Applications, 32 (3 ), 1–22. 10.1002/eap.2544
Morrison, J. L. (2003). Age‐specific survival of Florida's crested caracaras. Journal of Field Ornithology, 74 (4 ), 321–330. 10.1648/0273-8570-74.4.321
Newton, I. (1979). Population Ecology of Raptors. T. and A.D. Poyser.
Newton, I. , Marquiss, M. , & Rothery, P. (1983). Age structure and survival in a Sparrowhawk population. The Journal of Animal Ecology, 52 (2 ), 591. 10.2307/4574
Newton, I. , McGrady, M. J. , & Oli, M. K. (2016). A review of survival estimates for raptors and owls. Ibis, 158 (2 ), 227–248. 10.1111/ibi.12355
Ortega, E. , Mañosa, S. , Margalida, A. , Sánchez, R. , Oria, J. , & González, L. M. (2009). A demographic description of the recovery of the vulnerable spanish imperial eagle aquila adalberti. Oryx, 43 (1 ), 113–121. 10.1017/S0030605307991048
Payevsky, V. A. (2016). Sex‐biased survival and philopatry in birds: Do they interact? Biology Bulletin, 43 (8 ), 804–818. 10.1134/S1062359016080136
Peakall, R. , & Smouse, P. E. (2006). GENALEX 6: Genetic analysis in excel. population genetic software for teaching and research. Molecular Ecology Notes, 6 (1 ), 288–295. 10.1111/j.1471-8286.2005.01155.x
Ponnikas, S. , Ollila, T. , & Kvist, L. (2017). Turnover and post‐bottleneck genetic structure in a recovering population of Peregrine falcons Falco peregrinus. Ibis, 159 (2 ), 311–323. 10.1111/ibi.12460
QGIS Development Team . (2023). QGIS geographic information system. Open Source Geospatial Foundation Project. http://qgis.org/
R Core Team . (2023). R: A Language and Environment for Statistical Computing. https://www.r‐project.org/
Real, J. , & Mañosa, S. (1997). Demography and conservation of western European Bonelli's eagle Hieraaetus fasciatus populations. Biological Conservation, 79 (1 ), 59–66. 10.1016/S0006-3207(96)00100-0
Rousset, F. (2008). GENEPOP'007: A complete re‐implementation of the GENEPOP software for windows and Linux. Molecular Ecology Resources, 8 (1 ), 103–106. 10.1111/j.1471-8286.2007.01931.x 21585727
Rudnick, J. A. , Katzner, T. E. , Bragin, E. A. , Rhodes, O. E. , & Dewoody, J. A. (2005). Using naturally shed feathers for individual identification, genetic parentage analyses, and population monitoring in an endangered eastern imperial eagle (Aquila heliaca) population from Kazakhstan. Molecular Ecology, 14 (10 ), 2959–2967. 10.1111/j.1365-294X.2005.02641.x 16101766
Sæther, B. E. , & Bakke, Ø. (2000). Avian life history variation and contribution of demographic traits to the population growth rate. Ecology, 81 (3 ), 642–653. 10.1890/0012-9658(2000)081[0642:alhvac]2.0.co;2
Sandercock, B. K. (2020). Mark‐recapture models for estimation of demographic parameters. In D. L. Murray  & B. K. Sandercock (Eds.), Population ecology in practice (pp. 157–190). John Wiley & Sons Ltd.
Selwood, K. E. , McGeoch, M. A. , & Mac Nally, R. (2015). The effects of climate change and land‐use change on demographic rates and population viability. Biological Reviews, 90 (3 ), 837–853. 10.1111/brv.12136 25155196
Sergio, F. , Tavecchia, G. , Blas, J. , López, L. , Tanferna, A. , & Hiraldo, F. (2011). Variation in age‐structured vital rates of a long‐lived raptor: Implications for population growth. Basic and Applied Ecology, 12 (2 ), 107–115. 10.1016/j.baae.2010.11.004
Sergio, F. , Tavecchia, G. , Blas, J. , Tanferna, A. , & Hiraldo, F. (2021). Demographic modeling to fine‐tune conservation targets: Importance of pre‐adults for the decline of an endangered raptor. Ecological Applications, 31 (3 ), 1–12. 10.1002/eap.2266
Sibly, R. M. , & Hone, J. (2002). Population growth rate and its determinants: An overview. Philosophical Transactions of the Royal Society, B: Biological Sciences, 357 (1425 ), 1153–1170. 10.1098/rstb.2002.1117
Stahl, J. T. , & Oli, M. K. (2006). Relative importance of avian life‐history variables to population growth rate. Ecological Modelling, 198 (1–2 ), 23–39. 10.1016/j.ecolmodel.2006.04.001
Suh, A. , Kriegs, J. O. , Brosius, J. , & Schmitz, J. (2011). Retroposon insertions and the chronology of avian sex chromosome evolution. Molecular Biology and Evolution, 28 (11 ), 2993–2997. 10.1093/molbev/msr147 21633113
Sumasgutner, P. , Koeslag, A. , & Amar, A. (2019). Senescence in the city: Exploring ageing patterns of a long‐lived raptor across an urban gradient. Journal of Avian Biology, 50 (12 ), 1–14. 10.1111/jav.02247
Székely, T. , Liker, A. , Freckleton, R. P. , Fichtel, C. , & Kappeler, P. M. (2014). Sex‐biased survival predicts adult sex ratio variation in wild birds. Proceedings of the Royal Society B: Biological Sciences, 281 (1788 ), 20140342. 10.1098/rspb.2014.0342
Tenan, S. , Adrover, J. , Muñoz Navarro, A. , Sergio, F. , & Tavecchia, G. (2012). Demographic consequences of poison‐related mortality in a threatened bird of prey. PLoS One, 7 (11 ), 1–11. 10.1371/journal.pone.0049187
van Oosterhout, C. , Hutchinson, W. F. , Wills, D. P. M. , & Shipley, P. (2004). MICRO‐CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes, 4 (3 ), 535–538. 10.1111/j.1471-8286.2004.00684.x
Vili, N. , Nemesházi, E. , Kovács, S. , Horváth, M. , Kalmár, L. , & Szabó, K. (2013). Factors affecting DNA quality in feathers used for non‐invasive sampling. Journal of Ornithology, 154 (2 ), 587–595. 10.1007/s10336-013-0932-9
Vili, N. , Szabó, K. , Kovács, S. , Kabai, P. , Kalmár, L. , & Horváth, M. (2013). High turnover rate revealed by non‐invasive genetic analyses in an expanding eastern Imperial eagle population. Acta Zoologica Academiae Scientiarum Hungaricae, 59 (3 ), 279–295.
Waits, L. P. , Luikart, G. , & Taberlet, P. (2001). Estimating the probability of identity among genotypes in natural populations: Cautions and guidelines. Molecular Ecology, 10 (1 ), 249–256. 10.1046/j.1365-294X.2001.01185.x 11251803
Warkentin, I. G. , Espie, R. H. M. , Lieske, D. J. , & James, P. C. (2016). Variation in selection pressure acting on body size by age and sex in a reverse sexual size dimorphic raptor. Ibis, 158 (3 ), 656–669. 10.1111/ibi.12369
Weigmann, H. (1968). Reduction of disulfide bonds in keratin with 1,4‐dithiothreitol. I. Kinetic investigation. Journal of Polymer Science, Part A‐1: Polymer Chemistry, 6 (8 ), 2237–2253. 10.1002/pol.1968.150060819
White, G. C. , & Burnham, K. P. (1999). Program MARK: Survival estimation from populations of marked animals. Bird Study, 46 , S120–S139. 10.1080/00063659909477239
Whitfield, D. P. , Fielding, A. H. , Mcleod, D. R. A. , & Haworth, P. F. (2004). Modelling the effects of persecution on the population dynamics of golden eagles in Scotland. Biological Conservation, 119 (3 ), 319–333. 10.1016/j.biocon.2003.11.015
Wickham, H. (2016). ggplot2: Elegant graphics for data analysis. New York: Springer–Verlag. https://CRAN.R‐project.org/package=ggplot2
