Estimating species density and abundance is challenging but important for establishing conservation and management strategies. Significant progress has been made toward estimating Golden Eagle (Aquila chrysaetos) abundance in the conterminous United States of America (USA) but much less is known about eagle abundance in Alaska. Here, we paired migration count and GPS-tracking data collected near Gunsight Mountain, Alaska, in a Bayesian framework to estimate the number of Golden Eagles in south-central Alaska. We estimated 1204 (95% credible interval: 866, 1526) potentially breeding (≥4 yr old) Golden Eagles annually moved through the Gunsight Mountain migration corridor and summered over an area of 150,325 km2 in south-central Alaska, equating to a density of 0.80 potentially breeding eagles/100 km2. By extrapolating across the species' nesting range in Alaska (1,180,489 km2) and incorporating published productivity and age-specific survival rates for eagles <4 yr old into our hierarchical model, we estimated 12,717 (95% credible interval: 9043, 16,349) Golden Eagles of all ages occur in Alaska, annually. We propose this as a conservative statewide population estimate because we used methods that likely underestimated population size. Even so, our estimate is three to five times larger than previous estimates and likely represents about one quarter of the USA's population.
Introduction
Estimating the abundance of animal populations has been a central undertaking in wildlife ecology (Burnham and Overton 1979, Buckland et al. 1993, Collier et al. 2013) and such estimates are important components of conservation and management strategies (Dzul et al. 2013, Millsap et al. 2013, Dennhardt et al. 2015). Abundance can be relatively straightforward to estimate for sessile organisms occurring within a restricted range (Dzul et al. 2013); however, estimating abundance of mobile organisms with wide geographic distributions is more difficult (Whitehead 2002, Dennhardt et al. 2015). Geographically broad sampling strategies designed to estimate abundance over large areas can be expensive (Millsap et al. 2013, Nielson et al. 2014), require large numbers of people (Miller et al. 2016), and can be logistically difficult or near impossible to implement across remote landscapes (Dunn et al. 2005). These factors have collectively limited researchers from rigorously estimating densities and population sizes of many wildlife species in Alaska (Alaska Department of Fish and Game [ADFG] 2015), including the Golden Eagle (Aquila chrysaetos).
Golden Eagle abundance estimation has received substantial attention in the USA (US Fish and Wildlife Service [USFWS] 2009, Dennhardt et al. 2015, USFWS 2016), though most efforts have focused on estimates in the conterminous United States. In 2009 and 2014, the total population size of Golden Eagles in the United States was estimated at 39,000 and 40,000 individuals, respectively (USFWS 2016). Abundance estimates for specific regions within the USA were derived using late summer aerial surveys combined with Breeding Bird Survey data for the western USA (USFWS 2016) and migration count data and flight simulations for the eastern USA (Dennhardt et al. 2015). To estimate the population size in Alaska, the USFWS (2016) used the difference between the average late-summer estimate derived from surveys since 2006 and the average winter aerial survey estimates obtained from 2 yr of winter surveys within the conterminous western USA to estimate the number of migrant eagles present during winter. This number was then proportionally split by land area to approximate the proportion of migrants originating from Alaska (24%) and Canada (76%; USFWS 2016). The resulting estimate ranged from 1000 to 4000 mid-winter Golden Eagles that presumably originated from Alaska (USFWS 2016). The larger estimate (4000) was considered a liberal estimate and the midpoint (2500) a conservative estimate. The USFWS (2009) had previously estimated Alaska's summer population to be 2400 individuals.
McIntyre and Lewis (2016) observed 1364 migrating Golden Eagles during 9 d (39.5 hr) of fall migration counts from a single location in eastern Alaska. However, Golden Eagle migration in Alaska extends over a period of months (McIntyre et al. 2008, ADFG and USFWS unpubl. data) and the migration corridor through eastern Alaska covers an area much larger than what could be seen from a single count site (McIntyre and Lewis 2016, Bedrosian et al. 2018). Hence, McIntyre and Lewis (2016) suggested it was highly improbable that they counted approximately half of the state's estimated 2500 Golden Eagles from one site in 9 d, and the authors highlighted the need for a method to more accurately estimate eagle abundance in Alaska. Therefore, our goal was to provide a more plausible estimate by pairing migration count and GPS-tracking data from birds known to summer in Alaska, along with national demographic estimates, to coarsely estimate the total numbers of Golden Eagles that likely occur in Alaska.
Methods
Study Area. We conducted spring migration counts and captured migrating Golden Eagles in south-central Alaska during 2014–2018 in the Gunsight Mountain area, near Eureka, AK, USA (61°51.95′N, 147°20.96′W; Fig. 1). This area is approximately 155 km northeast of Anchorage, Alaska, in the Tahneta Pass between the Chugach and Talkeetna Mountain ranges. The two mountain ranges act as a natural funnel, concentrating migrating raptors moving westward through the pass where they can be counted by observers on the ground and sometimes captured. From this location, birds move westward to arrive on their summering grounds. We ultimately defined our study area as that over which Golden Eagles with satellite transmitters (hereafter, remotely tracked eagles) spent their summers (15 May–31 August) in south-central Alaska during 2014–2018 (Fig. 1 inset).
Figure 1.
Spring Golden Eagle count (stars) and trapping site locations (circles) at Gunsight Mountain, Alaska, and the approximate area over which remotely tracked eagles summered (inset).
![img-z3-1_496.jpg](ContentImages/Journals/rapt/55/4/JRR-20-81/graphic/WebImages/img-z3-1_496.jpg)
Field Methods. We conducted daily raptor migration counts from 7 March–15 May, 2016–2018 at pullouts along the Glenn Highway near milepost 120 (March, 61°51.95′N, 147°20.96′W) and milepost 121 (April and May, 61°52.91′N, 147°20.19′W; Fig. 1). We moved to milepost 121 after March to better observe the primary migratory flight line as it shifted seasonally. Counts were conducted by two experienced migration observers with at least several seasons of raptor migration experience; other HawkWatch International staff and local volunteers provided additional survey effort when available. Observers worked collectively to count migrating raptors; we did not use independent or double-observer methods. Weather permitting, we began observations at 0900 H and ended at 1700 H. We did not conduct counts on the five days when inclement weather (e.g., snow, fog) significantly obstructed views of the main migratory flight lines. All migrating raptors were identified to species, aged when possible, and tallied into 1-hr time intervals (e.g., 1000–1100 AKST). Local raptors (i.e., not actively migrating) were differentiated from migrants based on behavior (e.g., territorial displays), direction of flight, and unique characteristics (e.g., missing feathers, visible transmitters) and were excluded from counts.
We deployed 45-g solar-powered GPS/Argos platform transmitter terminals (Microwave Telemetry, Inc., Columbia, MD, USA) on 53 Golden Eagles (n= 15 in 2014, n = 23 in 2015, and n = 15 in 2016). Transmitters from three eagles failed to provide summer range locations and were not included in these analyses. We captured eagles between 17 March and 6 April annually within 7 km of the count sites (Fig. 1) and considered these eagles part of the population being counted. We used remote-controlled net launchers over carrion bait and placed transmitters on all eagles captured without discrimination by age, sex, body condition, or other variables until no transmitters remained for that year. We attached transmitters to eagles using custom-fit backpack harnesses and programmed transmitters to report 4–14 GPS locations/d depending on season and year (Eisaguirre et al. 2019). We filtered the data for outliers by speed and visual inspection and removed any locations that indicated a movement rate of >20 m/sec or appeared unreasonable. This removed only a few observations from our data set. All data were received through the Argos satellite system and managed in Movebank ( https://www.movebank.org/); we used only GPS locations that had a reported horizontal error of approximately ±18 m (Microwave Telemetry, Inc.) for our analyses. We used plumage to age eagles at time of capture (Liguori 2004), and all 50 eagles used for analyses were ≥4 yr old. Because eagles have been documented to first attempt breeding at 4 yr of age (Katzner et al. 2020) and because not all adults breed each year, we defined the population we studied and quantified in south-central Alaska as potential breeders (≥4 yr old).
Study Area Abundance Estimation. We implemented a Bayesian hierarchical model that integrated migration counts and observations of availability and detectability to estimate Golden Eagle abundance.
Availability. We defined availability as the proportion of remotely tracked eagles that likely were within view of the count site when observers were present (i.e., available to be counted). We mapped the observers' viewshed from count locations based on distances to landmarks and the maximum distance we could detect eagles with the naked eye or binoculars (<11 km; Fig. 2A). The boundaries of the viewshed were irregular and largely dictated by ridgetops of the surrounding mountains, which obscured eagles passing behind them. We used annual spring migration tracks from remotely tracked eagles that reported ≥5 GPS fixes/d on migration to assess availability of eagles at the count site. Not all transmitters provided adequate data to assess availability each year, and some were therefore excluded from our availability analysis. To account for uncertainty between discrete GPS locations of each eagle, we fit a dynamic Brownian bridge movement model to individual tracks for each spring migration (Farmer et al. 2010) in R (version 3.3.0) with the package move (Kranstauber et al. 2012, R Core Team 2016,).
Figure 2.
Example migration tracks, based on dynamic Brownian bridge movement models for spring 2016, of two Golden Eagles classified as (A) available and (B) unavailable to be counted based on the viewshed of observers on the ground (dashed line) near Gunsight Mountain, Alaska. Solid line denotes the 95% movement isopleth. Count site is identified by the crosshairs, and black dots indicate GPS locations of eagles during passage.
![img-z5-1_496.jpg](ContentImages/Journals/rapt/55/4/JRR-20-81/graphic/WebImages/img-z5-1_496.jpg)
To fit the model, we identified migration departure and arrival dates qualitatively. We defined the start migration point as the first point indicative of a northbound movement with no subsequent return to the individual's apparent winter area and the end migration point as the first point on the summer home range. The movement model was fit to these tracks plus approximately two margin lengths on each end to permit reasonable Brownian motion variance estimates for the entire migration track. The window and margin lengths, used for computing breakpoints in the dynamic Brownian motion variance (Kraunstaber et al. 2012), were set to an average of 3 d and 1 d, respectively, and calculated for each individual independently. We used 95% isopleths of the utilization distributions predicted by the movement models for each individual as estimates of the eagles' realized migratory movement corridors. We then overlaid movement isopleths along with location fix timestamps onto observer viewsheds and visually inspected for overlap.
We defined an eagle as available to be counted if it met two criteria: (1) the movement isopleth at least partially overlapped the observer viewshed (Fig. 2A); and (2) an observer was on duty for at least 1 hr during the possible passage window. If a movement isopleth did not overlap the observer viewshed (Fig. 2B) or no observer was on duty during the time of passage (including the 5 d when inclement weather precluded counts), we defined that eagle as unavailable to be counted.
Detectability. We assessed our subset of available eagles for their likelihood of being detected by on-the-ground observers. We defined detectability as the proportion of available eagles likely detected by observers. We classified a remotely tracked eagle as detected if observers counted any eagle during the interval when the movement isopleth passed the count site, irrespective of the recorded flight direction and locale of the available isopleth. We classified an eagle as undetected if observers saw no migrating eagles during that isopleth passage interval. We assumed the probability of eagle detection at a given count location and time was constant regardless of number or identity of observers, use of optics, age of eagles, weather, or flight conditions. We also used the average of two published detection probability estimates (Berthiaume et al. 2009 and Nolte et al. 2016) to create an informative prior in our hierarchical model of abundance (see below).
Counts. We adjusted eagle counts to exclude birds <4 yr old (hereafter, pre-breeders) because we did not collect movement data for these age classes and they can exhibit substantially different movements than potential breeders (Katzner et al. 2020). In 2017, highly experienced observers assigned ages to the greatest proportion of eagles counted at Gunsight Mountain in any year (89% of eagles aged; 86% of which we classified as potential breeders). We therefore multiplied each year's count by 0.86 and rounded it to the nearest integer to estimate the total number of potentially breeding eagles that migrated past Gunsight Mountain annually.
We assessed the independence of marked eagles using a Spearman test function in R (coin package) with a Monte Carlo distribution approximation to assess independence between our two samples: (1) average counts from 5-d passage windows; and (2) average daily passage of remotely tracked eagles from 5-d passage windows based on dynamic Brownian bridge movement models. We did not remove marked eagles from the count and compared marked eagles to all eagles. We used raw counts not corrected for availability or detection. We chose a 5-d window a priori because 5-d windows are often used to smooth count data, given large, daily migration count fluctuations. We rejected the null hypothesis that the two samples were independent (Z= 2.74, P= 0.002) and concluded that our sample of remotely tracked eagles was congruent with our sample of eagles counted. By concluding this, we assumed that (1) our availability and detectability estimates from remotely tracked individuals were representative of the population of migrating eagles being sampled, and (2) remotely tracked eagles were a representative sample of potentially breeding eagles migrating past Gunsight Mountain. Had we failed to reject the null hypothesis, our interpretation of the representativeness of our sample would have been more equivocal.
Hierarchical model. We implemented a Bayesian hierarchical model that integrated the observations of availability y and detectability z as well as the (adjusted) counts C to estimate the probabilities of an eagle being available pa and detected pd and the abundance of potential breeders in the study area Nad. For year t and individual i, the data model took the form:
and the process model took the form:
where X is a design matrix indicating the year; µβ, µα, and µλ are global availability, detectability, and potential breeder abundance, respectively; and Normal+ corresponds to the truncated (at zero) normal distribution. To complete the model hierarchy, we specified µβ ∼ Normal (0,1.752), µα ∼ Normal (2,0.52), µλ ∼ Normal+ (1000,3002), and InverseGamma(1,1) for all variances. The prior on µα is strongly informative, centered on the average of two previous detection probabilities estimated for Golden Eagles (89% in Berthiaume et al. [2009] and 84% in Nolte et al. [2016]).
We estimated parameters in the model with a custom Markov chain Monte Carlo (MCMC) algorithm, composed of a combination of Gibbs and Metropolis-Hastings updates, in R. We obtained an MCMC sample from the posterior with three chains of 200,000 iterations (discarding the first 100,000 as burn in) and assessed convergence to the posterior with traceplots. We summarized and reported results from the hierarchical model with posterior means and 95% credible intervals throughout.
Summer range. We used kernel density estimates (KDEs) to estimate the summer range size for potentially breeding eagles that used the Gunsight migration corridor; KDE is an approach commonly used in wildlife biology to estimate individual animal home ranges (Worton 1989, Kernohan et al. 2001, Kie et al. 2010, Kie 2013). Burt (1943) defined an animal's home range as the space traversed by an individual in its normal activities of foraging, mating, and caring for young. For this analysis, we expanded this definition to the population level and defined the summering range of a population as the collective space used by marked individuals during the summering period. This approach has been used in other taxa to identify population-level seasonal home ranges (e.g., Nicholson et al. 2016, Viengkone et al. 2018). Based on remote-tracking data and our understanding of this population's behavior, we defined the summering period as 15 May to 31 August. We believe this period represents a conservative estimate of the time over which the population is established on its summer range and avoids including data from late arriving or early departing eagles. To avoid biasing the KDEs toward birds with multiple years of summer locations, we used only the first summer of data for each bird even though eagles returned to the same areas annually. Due to duty cycles' inconsistency among deployment years, the year a transmitter was deployed had an inherent effect on the number of GPS locations reported. We therefore randomly sampled the track of each bird to match the number of locations provided by the duty cycle programmed to collect the fewest locations for a summer (eight locations per day and 872 for the summer). If a transmitter provided fewer than this number, we used all the locations it reported. We combined all of the resulting points from each eagle into one data set to estimate the collective summer range of the population. We summarized the mean number of fixes per bird ± standard errors, with ranges shown in parentheses.
We generated KDEs of the population's summer range in R with the package adehabitatLT version 0.3.25 (Calenge 2006). We chose the 95% isopleth of the KDE (hereafter 95% KDE; bivariate normal kernel) as the most appropriate representation of the population's use of space across the summer range (Kernohan et al. 2001). Choice of the bandwidth, or smoothing parameter h, is important to estimating home ranges with KDEs (Silverman 1986, Worton 1989, 1995, Seaman and Powell 1996, Seaman et al. 1999, Kernohan et al. 2001, Kie 2013). We therefore used two methods to select a bandwidth. The first method was the reference bandwidth (href; Silverman 1986, Worton 1989), which can result in over smoothing of multi-modal data and a positive bias in home range size estimates (Wand and Jones 1993, Kie 2013). In the second method, we used an adaptation of the “rule-based ad hoc” approach (hereafter ad hoc method), proposed by Kie (2013), which involves sequential scaling of href and visual assessment of the KDE until the KDE fits an a priori definition of home range. When the KDE fits the definition, the h used, and the KDE generated is accepted as the most appropriate for the animal's track. We defined a priori the summer range of this population as that which assigned the smallest continuous polygon for the primary range and a secondary polygon for a range comprised of two individuals that summered further west (i.e., smallest area that resulted in two polygons in the 95% KDE, Fig. 3). The ad hoc method resulted in a range much larger than the reference method, which itself is known to positively bias range size (Wand and Jones 1993, Kie 2013). Though we present results of both, we considered the ad hoc method more appropriate for cautiously estimating statewide abundance (a positive bias in our range size causes a negative bias in our statewide population estimate). We calculated eagles/100 km2 by dividing the population size by range size and multiplying by 100.
Figure 3.
Golden Eagle summer range estimates based on the collective space used by remotely tracked eagles trapped at Gunsight Mountain, Alaska 2014–2016 (n=50). Estimates created using two different 95% kernel density estimator (KDE) smoothing parameters: (A) Range One href, (B) Range Two had hoc. Gray areas denote the 95% KDE summer range; black dots are eagle GPS locations.
![img-z8-1_496.jpg](ContentImages/Journals/rapt/55/4/JRR-20-81/graphic/WebImages/img-z8-1_496.jpg)
We made the following assumptions for our summer range analysis: (1) remotely tracked eagles composed a representative sample of the migratory summering population of potentially breeding eagles; (2) all eagles in the adjusted counts summered within the calculated summering range; and (3) adjusted eagle counts represented the entire population of potentially breeding migratory eagles summering within the calculated summer range.
Statewide Abundance Estimation. To arrive at a statewide estimate that includes eagles of all ages, we first quantified the area over which Golden Eagles nest in Alaska by overlaying the breeding range map from Kochert et al. (2002) onto a georeferenced map of Alaska in ArcGIS 10.7 (Esri, Redlands, CA, USA) and clipped the state map by the breeding range polygon. We then generated posterior realizations from our model, such that for the kth MCMC iteration, where Astate is the area of the statewide summer range in Alaska and Asample is the area of the summer range of the remotely tracked sample of eagles. This provided us an estimate of the number of potential breeders in Alaska. To arrive at posterior realizations of total abundance including all age cohorts, we assumed a stable age distribution and incorporated previously estimated parameters. We specified the statewide abundances for potential breeder, fledgling, first year, second year, and third year age groups as follows:
∼Poisson(
),
∼Poisson
∼Poisson
∼Poisson
∼Poisso
.
Here, ψ and ϕ are the proportion of territorial eagles and productivity, respectively, and each s is age-specific survival. We used USFWS (2016) national estimates (1.13:1 floater to breeder ratio, 0.54 eagles fledged per occupied territory, and annual age-specific survival rates of 0.70, 0.77, and 0.84 for first-, second-, and third-year birds, respectively; survival rates were incorporated as Beta distributions parameterized by the USFWS), and we specified ψ∼ Normal(1/2.13, 0.012), ϕ∼ Normal(0.55, 0.0872), s1∼ Beta(353, 159), s2 ∼ Beta(287, 85), and s3 ∼ Beta(243, 47).
Finally, state wide total abundance was . All R code was archived at http://doi.org/10.5281/zenodo.4474843.
Results
We counted 1163, 705, and 1020 migrating Golden Eagles in 2016, 2017, and 2018, respectively, at Gunsight Mountain, Alaska. Our hierarchical model estimated availability and detectability as 0.68 (95% credible interval: 0.39, 0.89) and 0.95 (0.80, 0.99), respectively, although it suggested availability likely varied by year (Fig. 4). The model estimated 1204 (866, 1526) potentially breeding Golden Eagles migrated through the Gunsight Mountain migration corridor and summered within our study area, annually (Fig. 5).
Figure 4.
Probability distributions of eagle availability by year and summarized across years (global) for eagles migrating past Gunsight Mountain, Alaska, 2016–2018.
![img-z9-9_496.jpg](ContentImages/Journals/rapt/55/4/JRR-20-81/graphic/WebImages/img-z9-9_496.jpg)
Figure 5.
Probability distributions of potentially breeding Golden Eagle abundance in our study area in south-central Alaska from 2016–2018 and summarized across years (global).
![img-z10-1_496.jpg](ContentImages/Journals/rapt/55/4/JRR-20-81/graphic/WebImages/img-z10-1_496.jpg)
We collected a total of 52,658 GPS fixes from 15 May to 31 August from the 50 Golden Eagles in our study (1053 ± 39.2 per individual; range: 175–1518). Of those, we included 41,339 fixes (829 ± 19.3 per individual; range: 175–872) in our analyses. We estimated the total summer range for our south-central Alaska population at 108,945 km2 (Range One) and 150,325 km2 (Range Two; Fig. 3A and 3B, respectively) with resulting mean densities of 1.11 (0.79, 1.40) and 0.80 (0.58, 1.02) potentially breeding eagles/100 km2, respectively. The summer range consisted primarily of the western Alaska mountain range, as well as parts of the Talkeetna, Chugach, Kenai, and to a lesser extent, the Kuskokwim mountain ranges. These mountains support large areas of partially forested river valleys, shrub- and tundra-covered mountain slopes, alpine meadows, and barren rocky highlands. The summer range also included large areas of presumably lower quality eagle habitat including glaciers, ocean, and the lowland Susitna River valley.
We estimated the current breeding range of Golden Eagles in Alaska covers 1,180,489 km2 (69% of Alaska's 1,717,856 km2). By applying the Range One (108,945 km2) and Two (150,325 km2) summer area estimates, the posterior realizations from the hierarchical model that included demographic estimates for pre-breeding eagles suggested a total statewide population size of 17,547 (12,476, 22,553) and 12,717 (9,043, 16,349) Golden Eagles of all age cohorts, respectively.
Discussion
We combined Golden Eagle migration counts and movement data from remotely tracked eagles placed in a Bayesian hierarchical framework to conservatively estimate the number of potentially breeding Golden Eagles in our study area in south-central Alaska. We then applied this estimate to the geographic extent of the Golden Eagle's breeding range in Alaska along with published productivity and age-specific survival rates to arrive at a statewide late summer population estimate of approximately 12,700 eagles of all age classes. We propose this as a coarse conservative statewide population estimate. Though coarse, this estimate represents a three to five-fold increase in the number of eagles in Alaska compared to current USFWS population estimates (2500–4000; USFWS 2016). Hence, our study contributes substantial support to a previous investigation by McIntyre and Lewis (2016) that suggested the number of Golden Eagles in Alaska is much higher than previously estimated. In fact, our estimate shows Alaska likely has the largest population of any regional Golden Eagle management unit in the USA and supports roughly one quarter of the USA population.
Throughout this study, we purposely chose approaches that should guard against overestimating abundance, and we recognize this approach may seem unconventional. However, we consider it the most prudent method given the extrapolations necessary for a statewide estimate, lack of precision in our estimates, the substantial logistical challenges that currently prohibit a more robust, unbiased estimate for Alaska, and the conservation implications of potentially overestimating abundance. We believe our estimate likely underestimates the true abundance of Golden Eagles for several reasons. First, we likely overestimated availability (which lowers adjusted counts) by considering an eagle with any portion of its movement isopleth overlapping the observers' viewshed as available even if only the periphery of the isopleth overlapped. Second, we assumed a remotely tracked eagle was detected if any eagle was counted during that eagle's passage window, even though most eagles passed at distances too far to see a transmitter. This likely artificially inflated our detection probability estimate and lowered adjusted eagle counts. It is noteworthy that our eagle detectability estimate (95%) was quite high despite the strong prior used in the model, which was informed by previous estimates (89% in Berthiaume et al. 2009 and 84% in Nolte et al. 2016). Third, we assumed a constant eagle detection rate regardless of observer identity, weather, and flight conditions even though these factors are known to affect detectability (Berthiaume et al. 2009, Nolte et al. 2016). Accounting for these factors likely would have only reduced detection probability relative to our estimate and hence, raised adjusted counts. Fourth, we assumed all eagles in our population migrate and were accounted for in our migration counts. However, an unknown but potentially consequential proportion of Golden Eagles, including some tagged in this study, forego migration and overwinter in Alaska in some years (Kessel 1988; T. Booms, S. Lewis, unpubl. data). Any such eagles were not included in our counts and would again cause our population estimate to be biased low. Fifth, we selected the larger of two summer home range estimates, both of which were derived from methods known to biased ranges high (which result in lower density and population estimates because of the larger denominator). Lastly, we assumed no eagles arriving from other migration routes summered in our study area, which may be unlikely.
Our sample size of 50 remotely tracked Golden Eagles represents an immense effort and is much larger than any previous study of breeding-age Golden Eagles in Alaska (and many other locations), but it is still small compared to the size of the regional and statewide populations we estimated. Small sample sizes and population-level inference are known limitations of remote-tracking work, typically due to financial and logistical constraints (Hebblewhite and Haydon 2010), and are reasons that we chose methods that, if anything, likely underestimate abundance. We attempted to assess how representative our sample was by comparing the timing of passage between marked and all eagles, and we found no statistical difference. Additionally, past studies attempting to assess the number of remotely tracked individuals needed to quantify group-level space use concluded that tens of individuals were needed to approach an asymptote in median area estimates (Gutowsky et al. 2015). Relatively smaller sample sizes (n = 30) were needed to calculate space use during more sedentary periods of the annual cycle (e.g., breeding season as done here; Thums et al. 2018, reviewed in Sequeira et al. 2019). We therefore concluded that our sample of 50 tagged eagles represented a sufficient sample from the overall population of eagles that uses the Gunsight Mountain migration corridor.
We chose to not exclude substantial areas of water (including ocean) or glaciers from the south-central summer range estimate because one tracked eagle spent long periods of time on the coast and, given the discrete-time nature of tracking data, there is no way of knowing for certain that eagles do not make at least short forays over water or ice. Further, areas of presumably low-quality habitat are likewise included in the estimated statewide breeding range of the species. Thus, our approach should at least partially account for the variable nature of habitat quality and hence, eagle abundance, in a similar fashion to the large-scale delineation of the nesting range in Alaska.
We purposely excluded pre-breeding Golden Eagles from our study area population and density estimates because (1) we did not have GPS data with which to quantify their summer movements, (2) pre-breeding eagles often exhibit substantially different summer movements than potential breeders (McIntyre et al. 2008, Katzner et al. 2020), and (3) many pre-breeding eagles move outside the study area and Alaska's breeding range (e.g., coastal plain of northern Alaska) during their pre-breeding years (McIntyre and Lewis 2018). It is therefore not surprising that the proportion of potentially breeding eagles in our study area as estimated by the age composition of the count data (86%) was higher than what the national estimates at equilibrium would have predicted (78%, USFWS 2016). Our methods would not have been appropriate to estimate densities for pre-breeding eagles because an unknown portion of them do not occur in our study area or the statewide breeding distribution. Instead, we used national demographic rates (USFWS 2016) to calculate how many younger eagles are produced, survive, and occur in each pre-breeding age class in Alaska. We decided to use national estimates instead of local ones because few such estimates are available for Alaska, and our goal was to use approaches consistent with USFWS (2016) whenever possible. In the longest and most comprehensive study of Golden Eagles in Alaska, McIntyre and Schmidt (2012) estimated average Golden Eagle productivity in Denali National Park (just north of our study area) at 0.62 fledglings per occupied nesting territory, which is higher than the national estimate we used (0.54 fledglings). Using the Alaskan estimate would have increased our statewide population estimate by 300 eagles to a total of approximately 13,000 eagles.
Our density estimate (0.80 potentially breeding eagles/100 km2) is at the low end of other North American eagle density estimates, lending support to our conclusion that it is both reasonable and conservative in nature. Millsap et al. (2013) calculated late-summer Golden Eagle densities using western United States aerial survey data and presented the following mean densities within select Bird Conservation Regions (BCR; converted from eagles/km2 to eagles/100 km2): Southern Rockies/ Colorado Plateau (BCR 16) = 0.8; Great Basin (BCR 9) = 0.9; Northern Rockies (BCR 10) = 1.5; and Badlands and Prairies (BCR 17) = 2.7. One important caveat to this comparison is that the Millsap et al. (2013) estimates were calculated from data collected after the breeding season and therefore include recently fledged young and other pre-breeding eagles. Our study area density estimate only includes potentially breeding birds (i.e., ≥4 yr old). If pre-breeding eagles had been included, our estimate would have been slightly higher.
Estimating availability and its potential interannual variation at raptor migration sites is important to estimating actual population size as done here. It is also critically important in evaluating broader population indices derived from migratory count data to assess population change over time (Bildstein et al. 2008, Nolte et al. 2016). Importantly, our study is among the first to assess interannual variation in availability of birds at a count site (Fig. 4). Both Farmer et al. (2010) and Dennhardt et al. (2015) quantified availability of migrating birds at count sites but neither estimated how these measures change over time, which is inherently important to ascribing cause to changes. Without such an understanding, it is difficult to differentiate between population-induced changes in counts versus changes in bird availability at count sites. Our data suggest that availability may be influencing changes in raw count data, at least for one migration site.
We recognize that our estimates are based on key assumptions and that we extrapolate beyond our study area to arrive at a coarse estimate of Golden Eagle abundance in Alaska. The size, remoteness, and diversity of habitats in Alaska represent substantial challenges to using more rigorous methods to estimate statewide abundance. The current lack of resources to overcome these challenges leaves few or no feasible alternatives with which to reasonably estimate Golden Eagle abundance for the management community. Though based on relatively simple count and movement data, we contend our approach is logical, reasonable, and based on the best available science. Additionally, even without extrapolating beyond our study area, our corrected migration count of approximately 1200 potentially breeding golden eagles that summer on about 9% of the Alaska's total land mass is almost half of the current USFWS conservative population estimate for all of Alaska. If the USFWS estimate is accurate, then the other half of Alaska's eagles would be distributed over the remaining 91% of Alaska at extraordinarily low densities—a conclusion we find unlikely. Finally, based on our collective firsthand experience surveying and working with Golden Eagles across Alaska (TLB, 15 yr; SBL, 10 yr; JME, 7 yr), we consider our estimate more plausible and reflective of what we see in the field than previous population estimates. The true Golden Eagle population size in Alaska remains uncertain but given the information we present here and our field experiences, it is likely larger, possibly much larger, than even our conservative estimate of 12,700 eagles.
Acknowledgments
Funding was provided by the Alaska Department of Fish and Game through the Federal Aid State Wildlife Grant Program, Anchorage Audubon Society, Eppley Foundation for Research, HawkWatch International, and the US Fish and Wildlife Service. We thank hawk migration counters R. Rubenthaler, C. Davis, D. Crowson, F. Nicoletti, F. Simone, and S. Nelson for their hard work collecting data, and C. McIntyre, J. Schmidt, G. Pendleton, B. Millsap, J. Pagel, G. Zimmerman, and D. Oleyar who contributed valuable perspectives and suggestions in the development of this manuscript. We also thank A. Duerr and an anonymous peer reviewer for providing additional helpful suggestions. J. Roach assisted with the GIS analysis. The findings and conclusions in this article are those of the authors and do not necessarily represent the views of the US Fish and Wildlife Service.