Raptors that inhabit sagebrush steppe and grassland ecosystems in the western United States may be threatened by continued loss and modification of their habitat due to energy development, conversion to agriculture, and human encroachment. Actions to protect these species are hampered by a lack of reliable data on such basic information as population size and density. We estimated density and abundance of nesting pairs of Ferruginous Hawks (Buteo regalis) and Golden Eagles (Aquila chrysaetos) in sagebrush steppe and grassland regions of Wyoming, based on aerial line transect surveys of randomly selected townships. In 2010 and 2011, we surveyed 99 townships and located 62 occupied Ferruginous Hawk nests and 36 occupied Golden Eagle nests. We used distance sampling to estimate a nesting pair density of 94.7 km2 per pair (95% CI: 69.9–139.8 km2) for Ferruginous Hawks, and 165.9 km2 per pair (95% CI: 126.8–230.8 km2) for Golden Eagles. Our estimates were similar to or lower than those from other studies in similar locations in previous years; thus, we recommend continued monitoring to determine trends in nesting pair density over time. Additionally, we performed double-observer surveys on a subset of transects with a helicopter as the second observation aircraft. We estimated probability of detecting occupied nests from fixed-wing plane versus helicopter, as well as time and expense of each survey mode. Although observers surveying from helicopters were 1.19 and 1.12 times more likely to detect Ferruginous Hawk and Golden Eagle occupied nests, respectively, the helicopter survey was 4.55 times costlier due to longer flight time and the higher hourly costs. Thus, when systematically surveying large areas, we found cost and time of the helicopter surveys outweighed the increase in nest detection.
Ferruginous Hawks (Buteo regalis) and, to a lesser extent, Golden Eagles (Aquila chrysaetos) rely on grassland and sagebrush steppe ecosystems for foraging and nesting. Large undisturbed areas of these ecosystems are, however, becoming increasingly rare (Noss et al. 1995, Knick et al. 2003, Brennan and Kuvlesky 2005). Anthropogenic activities, altered fire regimes, and invasive species have diminished the sagebrush steppe ecosystem to only 56% of its historic range (Davies et al. 2011), and Wyoming contains one-fourth of the remaining sagebrush steppe habitat in the United States (Knick et al. 2011). Grassland ecosystems are threatened globally, due in part to habitat modification and conversion (Brennan and Kuvlesky 2005, Hoekstra et al. 2005).
Ferruginous Hawks are grassland or desert shrub habitat specialists (Schmutz 1984) and have a distribution centered on Wyoming. Although Golden Eagles are widespread habitat generalists, the grassland and shrubland habitats in which they are sympatric with Ferruginous Hawks in Wyoming have recently experienced elevated disturbance from energy development (Copeland et al. 2011). For example, the number of producing natural gas wells in Wyoming increased from 2600 in 1990 to 22,171 in 2012, a 753% increase (U.S. Energy Information Administration 2014). A decrease in suitable, undisturbed habitat may preclude future nesting efforts by Ferruginous Hawks and Golden Eagles, both of which are sensitive to disturbance while nesting (Suter and Joness 1981, White and Thurow 1985). Both species are also susceptible to mortality from power-lines, nest destruction, and human persecution (Kochert and Steenhof 2002, Lehman et al. 2007, Millsap et al. 2015). Given these concerns, documenting current nesting pair density of these species in Wyoming is central to conservation planning.
Surveys to determine density and abundance of raptors are frequently hampered by imperfect detection or nonrandom sampling methods (Andersen et al. 1985, Andersen 2007). Raptor surveys are often conducted near or along roads due to easy access and ability to cover large areas (Fuller and Mosher 1987). Roadside surveys for raptors may be biased, however, due to influence of roads on probability of detection and/or probability of occurrence (Andersen et al. 1985, Millsap and LeFranc 1988). In contrast, density estimates based on aerial surveys of a randomly sampled population should be more accurate, provided methods are used to account for imperfect detection. Line-transect distance-sampling methods, which control for failure to detect all objects on a transect, produce estimates of true density that can be compared over space and time (Buckland et al. 1993).
A major assumption of distance sampling is that all objects that fall directly on the survey transect are detected; if this assumption is not met, density will be underestimated (Buckland et al. 1993). This assumption is, however, frequently violated when surveying from a fixed-wing plane because visibility directly under the plane is limited. Double-observer methods, in which a second observer independently surveys the same transects as the first observer, can be used to determine the probability of detection directly on the transect (Borchers et al. 1998, Good et al. 2007). The use of a helicopter to obtain a second observation should be particularly effective as the glass cockpit affords a wider range of vision directly on or near the transect. Helicopters cost much more per hour to operate, however, and fly at a slower speed than a fixed-wing plane. A study by Booms et al. (2010) found that estimates of detection probability from a fixed-wing plane were generally similar or greater than those from a helicopter when searching for nests of cliff-nesting raptors, although their study did not conduct a direct comparison of the two. Fixed-wing planes may be less suitable when counting eggs or young, or surveying rugged terrain (Booms et al. 2010). Thus, the decision to undertake helicopter surveys must consider tradeoffs among objectives of the study, safety, and cost.
Spatial extent and location of a study area can also influence density estimates. For example, small study areas may not adequately capture large-scale spatial heterogeneity, and areas determined solely by proposed development boundaries may be biased by project location. Our study encompasses a large spatial scale, which allows us to estimate density for the majority of shrubland and grassland ecosystems in Wyoming. Previous estimates of nesting pair density for Ferruginous Hawks and Golden Eagles were derived from relatively smaller study areas (Table 3; Kochert et al. 2002). We acknowledge that changes in nesting pair density may misrepresent population trend without additional demographic information (Van Horne 1983). Therefore, our estimate of density represents only an index of population health, given our lack of knowledge of vital population rates such as recruitment or mortality.
Table 3.
Nesting pair density estimates for Ferruginous Hawks over their distribution in western North America.
For this study, we conducted extensive aerial surveys across sagebrush steppe and grassland habitats in Wyoming for Ferruginous Hawks and Golden Eagles in 2010 and 2011. Our main goal was to provide a current estimate of nesting pair density for these raptors from which to evaluate future trends relative to human disturbance of sagebrush steppe and grassland ecosystems. Specifically, we aimed to: (1) document nesting pair density of Ferruginous Hawks and Golden Eagles using distance sampling and double-observer methods across sagebrush steppe and grassland habitat of Wyoming, and (2) estimate detection probability of nests occupied by Ferruginous Hawks and Golden Eagles from fixed-wing plane versus helicopter. We also determined which aircraft (fixed-wing versus helicopter) was most efficient in terms of survey cost and occupied nest detection.
Methods
Study Area.
Our study area included non-mountainous sagebrush steppe and grassland regions of Wyoming, U.S.A. Because Ferruginous Hawks were the initial primary focus of the study, we further limited our study area to townships within the previously modeled distribution of Ferruginous Hawks in Wyoming (Keinath et al. 2010; total area = 110,112 km2; Fig. 1). We divided our study area based on level III ecoregions of Chapman et al. (2004), which categorize areas based on ecosystem similarity and environmental resources. The eastern portion of our study area, made up of the High Plains (9207 km2, 99 townships) and Northwestern Great Plains (29,295 km2, 315 townships) ecoregions, is primarily mixed grass prairie; the western and central portion, the Wyoming Basin ecoregion (71,610 km2, 770 townships), is primarily Wyoming big sagebrush (Artemisia tridentata) or salt desert scrub (Atriplex spp.; Chapman et al. 2004). Elevation in our study area ranged from 940 masl to 2200 masl, with a mean of 1780 masl; precipitation averaged 15 cm to 40 cm per year (Knight 1996).
Survey Methods.
Our sample unit was a Public Land Survey System “township,” a square unit of land with an area approximately 93 km2. We chose this unit to maximize our chances of detecting Ferruginous Hawks and Golden Eagles because the area is large enough to potentially incorporate home ranges of several nesting pairs of either species (Dechant et al. 2002, DeLong 2004). Potential sample townships were chosen if their centroid was within the previously modeled distribution of Ferruginous Hawks in Wyoming (Keinath et al. 2010) or contained historical records of occupied Ferruginous Hawk nests. We stratified townships by geographical ecoregion (Chapman et al. 2004), with the number of townships selected for survey in each ecoregion proportional to its area. For other study objectives, we selected townships to include three oil and natural gas well densities: none (0 wells/township), low (1−30 wells/township), and high (≥31 wells/township), with equal numbers of townships in each stratum.
We used fixed-wing aircraft to search for Ferruginous Hawk and Golden Eagle nests in the townships in our sample between 14 April and 15 May 2010 and 2011. We flew over each township at an altitude of 60 m above ground level and an average approximate speed of 130 km/hr, following 16 parallel north–south transects the length of the township sample unit (9.54 km) with a spacing of 600 m to ensure complete coverage (Fig. 2). To account for possible nonindependence between transects, we treated the 16 north–south transects in each township as a single transect 153.6 km long (16 × 9.54 km). Surveys were conducted with two fixed-wing planes (Bellanca Scout and Piper PA 18), each carrying one pilot and one technician, both of whom contributed to nest detection simultaneously; all had previous experience detecting raptors. We used GPS units to record actual flight paths and mark locations of all observed stick nests. We calculated distance to each nest by plotting recorded nest locations in ArcGIS v. 10.0 (ESRI, Redlands, California, U.S.A.) and measuring perpendicular distance from the nest to the transect. We recorded species and status of all stick nests located; we defined an “occupied” Ferruginous Hawk or Golden Eagle nest as one that contained eggs, young, incubating, or defending adults (Steenhof and Newton 2007); we refer to “unoccupied” nests as any nest that does not meet these criteria, including nests that were dilapidated or in disrepair, as long as they could still be distinguished by observers as a stick nest.
On a subset of townships (n = 57) in 2010, we resampled approximately 32 km of the transect in a helicopter (Bell 47 Soloy), flying at 80 km/hr and following the same flight path as the fixed-wing plane for three lengths of the township, with start point randomly selected prior to the survey. The helicopter pilot and a single observer acted as a team to record nest status and location using the same methods as the fixed-wing plane; the same pilot and observer performed all helicopter flights. We used GPS flight tracks and nest locations recorded during both fixed-wing and helicopter surveys to determine which nests were detected during one or both surveys. We plotted the location of all nests using ArcGIS, and considered nests to be repeat detections if they were within 100 m of each other and located on the same nest substrate. Helicopter surveys took place an average of 8 d (SD = 9.8 d) before or after the fixed-wing plane survey; given this short lag-time, we expected violations of the assumption of demographic closure between survey occasions to be minimal. Nonetheless, we further accounted for the possibility that nest failure or establishment occurred between survey occasions by considering detections of unoccupied nest structures during one survey as successful repeat detections if the unoccupied nest was located ≤100 m from an occupied nest detected during the other survey. We assume independence between the helicopter survey and the fixed-wing because surveys did not take place at the same time and surveyors did not share nest locations.
Nesting Pair Density Estimation.
Distance sampling
We used distance sampling to estimate nesting pair density for each species over the entire study area and within each geographic ecoregion (Buckland et al. 1993, Pollock et al. 2002). We binned distance data into 50-m intervals to account for location inaccuracy from GPS points taken while in the air, and truncated the top 5% of distances to avoid a long-tailed distribution (Thomas et al. 2010). To determine whether data between years 2010 and 2011 could be pooled, we performed a separate distance analysis for each year, and compared detection probabilities from the two years using a two-sample t-test with the statistical software R (R Core Team 2013).
We used program DISTANCE (Thomas et al. 2010) to calculate separate detection functions for Ferruginous Hawk and Golden Eagle occupied nests detected from the fixed-wing plane. We calculated nesting pair density separately for each ecoregion to investigate regional differences in abundance; however, sample sizes were insufficient to calculate region-specific detection functions. Therefore, we followed Thomas et al. (2010) and used the “Multiple Covariates Distance Sampling” engine in program DISTANCE to estimate a single detection function for the entire study area, and used this detection function to estimate nesting pair density separately for each region. Half-normal and hazard-rate detection functions with cosine and simple polynomial expansions were considered to model the detection function that best fit our data. We ranked models using Akaike’s Information Criterion for small sample sizes (AICc; Akaike 1973) and chose the model with the lowest value as the best fitting (Burnham and Anderson 2002).
We calculated density using the following equation, based on Buckland et al. (1993):
where n is the number of occupied nests, ĝ(W1) is the probability of detecting an occupied nest on the transect line, L is the length of the surveyed transect, and µ is the effective strip width as determined by distance modeling. This equation compensates for nests missed at the transect line, as well as imperfect detection that increases with distance from the transect line (Good et al. 2007). To determine ĝ(W1), the probability of detecting an occupied nest on the transect line, we used a logistic regression model with data from transects that were surveyed by both fixed-wing plane and helicopter (Good et al. 2007). We treated occupied nests detected while in the helicopter as random samples of Ferruginous Hawk and Golden Eagle nests available for detection. For each occupied nest detected by helicopter, we determined whether fixed-wing aircraft had detected or missed the nest, and used this binary value as the dependent variable; we used the distance of each occupied nest from the transect as the predictor variable. We then used the resulting equation to predict probability of detection at distance 0 for the fixed-wing plane.
We estimated variance around nesting pair density using a bootstrap procedure. We sampled individual transects within ecoregions with replacement 1000 times and calculated nesting pair density from the resulting number of occupied nests for each iteration. We performed this procedure separately for transects within each ecoregion and then summed occupied nests from all ecoregions for a density estimate over all grassland and shrubland ecosystems in our study area. We then used the middle 95% of the distribution of these values as a 95% confidence interval (Good et al. 2007, Manly 2007). We calculated nesting pair abundance based on this density estimate and total area of each region.
Double-observer sampling
Before beginning our survey, we considered the possibility that we would not locate enough occupied nests using the fixed-wing plane to accurately estimate nesting pair density using distance sampling methods. Therefore, we used a separate approach based on double-observer methods similar to Nichols et al. (2000), but with temporally sequential rather than simultaneous observers (Conway and Simon 2003), to calculate a second set of detection probabilities for each species from the fixed-wing plane. This method estimates probability of detection for primary and secondary observers based on number of occupied nests each detects as well as the number one detects and the other misses. We considered our primary observer to be the team in the fixed-wing plane and our secondary observer to be the team in the helicopter; observers were temporally separated by an average of 8 d, which we assumed to be a sufficient lag-time to ensure independence between observers while maintaining closure in the population. We calculated probability of detection from the plane as the number of occupied nests detected by both survey methods divided by total number of occupied nests detected only by helicopter. We calculated probability of detection from the helicopter as the number of occupied nests detected by both survey methods divided by total number of occupied nests detected only by plane. We estimated variance following Equation 8 in Nichols et al. (2000). We then used these probabilities of detection to calculate nesting pair density using all occupied nests detected by the fixed-wing plane on all surveyed townships with the following equation:
where n is the number of occupied nests detected per township, L is township length, W is township width, and P̂is detection probability as estimated by the double-observer method. Variance for this nesting pair density estimate was estimated using the same bootstrap procedure detailed above.
Survey Efficiency.
We also used our double-observer data to determine relative ability of observation teams in helicopters versus fixed wing planes to detect occupied nests. Based on double-observer probabilities of detection, we calculated amount of time each survey mode took per occupied nest, cost per occupied nest, and relative efficiency based on their differing probabilities of detecting an occupied nest. We included only occupied nests located on the subset of transects that were surveyed by both fixed-wing plane and helicopter. We summed total distance included in the double-observer survey, and used this to determine how long each survey mode took based on its average speed. We divided this survey time by number of occupied nests detected by each survey mode to determine survey time per occupied nest.
Results
Ferruginous Hawks.
We surveyed 60 townships from 15 April to 14 May 2010 and 39 additional townships from 18 April to 7 May 2011. In total, we surveyed approximately 15,206 km of transects and found 62 occupied Ferruginous Hawk nests within the 99 surveyed townships. We located 46 occupied nests in 2010 and 16 occupied nests in 2011 from the fixed-wing plane. Of these, five were in the High Plains, 14 were in the Northwestern Great Plains, and 43 were in the Wyoming Basin (Fig. 3). Occupied Ferruginous Hawk nests were located on a variety of substrates including 17 on artificial nest platforms, 12 ground nests, 13 on rocks or rims, and 20 in trees (13 in narrowleaf cottonwood (Populus angustifolia), four in junipers (Juniperus spp.), and three in other deciduous trees).
Distance sampling
In 2010, probability of detection as determined by program DISTANCE was 0.66 (SE = 0.05), and in 2011 probability of detection was 0.68 (SE = 0.13). We found no difference in detection probability between 2010 and 2011 (P = 0.83, t = 0.215, df = 66); thus, we pooled data over both years for the rest of the analysis.
Detection probability of occupied nests from the fixed-wing plane on the transect line as estimated by logistic regression was 0.71 (95% CI: 0.25–0.95). Based on AICc, we selected a half-normal key function with a cosine expansion to model the detection function for Ferruginous Hawks. The global detection probability associated with the distance analysis of both years combined was 0.55 (95% CI: 0.45–0.66), and effective strip width on transects was 273.5 m. We estimated a total density of 94.7 km2 per nesting pair (95% CI: 69.9–139.8 km2) in our study area. From this density, we calculated a total of 1163 nesting pairs of Ferruginous Hawks (95% CI: 788–1575 pairs) in our study area during the 2010–2011 nesting seasons (Table 1).
Table 1.
Nesting pair density and abundance estimates of Ferruginous Hawks in our Wyoming study area by ecoregion, 2010–2011. The upper table provides estimates from distance analysis, lower table is results from the double-observer survey.
Double-observer sampling
On double-observer surveyed transects (n = 57), we detected a total of 30 occupied Ferruginous Hawk nests from fixed-wing plane and helicopter combined. Of the 30 occupied nests found, 18 were detected by both plane and helicopter, eight were detected by helicopter but not plane, and four were detected by plane but not helicopter. We estimated a detection probability of 0.69 (SE = 0.023) from the fixed-wing plane and 0.82 (SE = 0.021) from the helicopter. The joint probability of detection from both survey modes was 0.94 (SE = 0.041). Using the fixed-wing plane probability of detection and occupied nests from all surveyed transects (n = 62 nests), we estimated a density of 99.5 km2 per nesting pair (95% CI: 73.4–146.9) and an abundance of 1107 (95% CI: 750–1500) nesting pairs of Ferruginous Hawks in the study area (Table 1).
Golden Eagles.
For Golden Eagles, we located 36 occupied nests during surveys of all 99 townships; 20 in 2010 and 16 in 2011. We found four occupied nests in the High Plains region, 16 in the Northwestern Great Plains region, and 16 in the Wyoming Basin (Fig. 3). Fourteen of the 36 occupied nests were found on cliffs, 17 in cottonwood trees, one on a man-made structure, and four on rocks or rims.
Distance sampling
Due to small sample sizes, we could not reliably estimate a detection function for each year to determine whether years could be pooled. We decided to proceed on the assumption that detectability of Golden Eagles was constant across years because this was true for Ferruginous Hawks, and the same survey methods and observers were used for both species in both years. Thus, we pooled detection data for Golden Eagles over years.
Detection probability of occupied nests from the fixed-wing plane on the transect line as estimated by logistic regression was 0.74 (95% CI: 0.20–0.97). We selected a half-normal key function with a cosine expansion to model the detection function for Golden Eagles. The global detection probability associated with the combined distance analysis of both years was 0.59 (95% CI: 0.44–0.78), and effective strip width was 235.3 m. We estimated a mean density of one nesting pair per 165.9 km2 (95% CI: 126.8–230.8 km2) across all ecoregions, and a total of 664 nesting pairs of Golden Eagles (95% CI: 477–868 pairs) within our study area during the 2010–2011 nesting seasons (Table 2).
Table 2.
Nesting pair density and abundance estimates of Golden Eagles in our Wyoming study area by ecoregion, 2010–2011.a The upper table provides estimates from distance analysis, lower table is results from the double-observer survey.
Double-observer sampling
On double-observer surveyed transects, we detected 13 occupied Golden Eagle nests. Of occupied nests, six were detected by both fixed-wing plane and helicopter surveys, four were detected by helicopter but not plane, and three were detected by the plane but missed by the helicopter. Detection probability from the plane was 0.60 (SE = 0.061) and from the helicopter 0.67 (SE = 0.062); combined detection probability was 0.87 (SE = 0.12). Using the fixed-wing plane probability of detection and occupied nests from all surveyed transects (n = 36 nests), we estimated a density of 149.0 km2 per nesting pair (95% CI: 127.7–233.2 km2) and a nesting pair abundance of 739 (95% CI: 472–862) Golden Eagles in the study area (Table 2).
Survey Efficiency.
Our average flight speed during nest surveys was 130 km/hr in the plane and 80 km/hr in the helicopter. The cost per hour for the plane and the helicopter increased annually but the helicopter averaged 2.8 times more expensive per hour than the plane (helicopter: $650/hr in 2010, $700/hr in 2011; plane: $230/hr in 2010, $250/hr in 2011). Based on a flight distance of 32.3 km per township, after 57 townships surveyed, the helicopter took 1.63 times as long to complete the survey as the fixed wing plane. Thus, the helicopter cost 4.55 times more than the plane for the same survey, after multiplying survey time by cost. In terms of detection probability, the observers surveying in a helicopter were 1.19 times more likely than those in a survey plane to detect an occupied Ferruginous Hawk nest, and 1.12 times more likely to detect an occupied Golden Eagle nest. During the double-observer survey, the observers surveying by plane detected a total of 22 occupied Ferruginous Hawk nests and nine occupied Golden Eagle nests, whereas those in the helicopter detected 26 occupied Ferruginous Hawk nests and 10 occupied Golden Eagle nests. Based on 2011 rates, cost per occupied Ferruginous Hawk nest for the helicopter survey was 3.85 times greater than for the fixed-wing plane. The cost per occupied nest for Golden Eagles was 4.09 times greater for the helicopter than the survey plane.
Discussion
We used both line-transect distance sampling and double-observer sampling methods to determine nesting pair density of Ferruginous Hawks and Golden Eagles throughout the distribution of Ferruginous Hawks in Wyoming. Because we accounted for probability of detection in our model, our estimates can be used for direct comparison with nesting pair density estimates from other areas. Our stratified random sampling scheme allowed us to make inference to the entire distribution of Ferruginous Hawks in the state of Wyoming, and provide estimates of nesting pair density for both species that will be useful for assessing current status of nesting pairs and as a baseline for future monitoring.
Our total estimate of 94.7 km2 per Ferruginous Hawk nesting pair suggests that density of nesting Ferruginous Hawks in Wyoming is slightly lower than previous nesting pair density estimates across the range of the species (Table 3). Numerous studies have indicated Ferruginous Hawk populations may be in decline (Woffinden and Murphy 1989, Schmutz et al. 1992, Olendorff 1993), with potential causes given as increased disturbance and loss of habitat from energy development and agriculture (Woffinden and Murphy 1977, Houston and Bechard 1984). Energy development is extensive across Wyoming, and is particularly intense within the distribution of Ferruginous Hawks (Keinath et al. 2010, Wyoming Oil and Gas Conservation Commission 2014). Our estimates document abundance of nesting Ferruginous Hawks and Golden Eagles given current levels of energy-related disturbance, and provide a benchmark for comparison as development continues to increase in our study area.
For Golden Eagles, we found one nesting pair per 165.9 km2 over our total survey area. Our study area and regional estimates of nesting pair density are generally lower than those from previous studies (Kochert et al. 2002). For instance, Phillips et al. (1984) found an average of one nesting pair per 60 km2 over 12 survey areas in Wyoming, with a range of 34–89 km2 per survey area, while Camenzind (1969) found 100–119 km2 per pair in Utah, and Kochert (1972) found 66 km2 per pair in southwest Idaho. This difference may be due in part to the low abundance of leporids, which are the main prey of Golden Eagles (MacLaren et al. 1988), during our study. Abundance of jackrabbits (Lepus spp.) and cottontail rabbits (Sylvilagus spp.) in Wyoming has been shown to follow an approximate 7- to 8-yr cycle (Gross et al. 1974, Fedy and Doherty 2011), although in recent years the population does not appear to have reached the peak in numbers seen previously (Oakleaf et al. 2014). Based on Fedy and Doherty (2011), the last population peak was in 2006, which indicates a low in population abundance should have occurred during our study period in 2010 or 2011. Additionally, because our survey was limited to the known distribution of Ferruginous Hawks in Wyoming, the density estimate we provide for Golden Eagles includes only those nesting in sagebrush steppe and grassland communities; eagles nesting in mountainous terrain were not included in this population estimate (Kochert et al. 2002).
Our results represent only an index of population health, as actual estimates of population vital rates are not available. Nesting pair density may underestimate the actual population of adults capable of reproduction, as territory availability or breeding opportunity may force some adult raptors into nonbreeding “floater” status (Hunt 1998). Ayers et al. (2009) found several nonbreeding Ferruginous Hawks and 11 breeding pairs in a 783 km2 study area near Baggs, Wyoming, whereas Schmutz et al. (2008) found no evidence of a nonbreeding adult population of Ferruginous Hawks in a study in western Canada. Golden Eagles have been generally found to maintain floaters in most populations (Kochert et al. 2002). Anecdotal results from our study agree with these findings: during surveys in 2010 we recorded 101 Golden Eagles and 13 Ferruginous Hawks that were perched or soaring and seemingly not associated with a nest. Thus, our estimate of nesting pair density may be a conservative estimate of the entire adult Ferruginous Hawk or Golden Eagle population, depending on the unknown abundance of nonbreeding adults.
Nesting pair density estimates for Ferruginous Hawks were highest in the Wyoming Basin, while Golden Eagle densities were highest in the Northwestern Great Plains and High Plains. Estimated densities for the two species were similar in the Northwestern Great Plains and High Plains ecoregions, but Ferruginous Hawks were approximately three times as abundant as Golden Eagles in the Wyoming Basin. The Wyoming Basin is composed largely of sagebrush communities, while the Northwestern Great Plains is made up of mixed-grass prairie. Both areas have extensive oil and gas development, and are otherwise relatively sparsely populated by humans. As mitigation for creating habitat disturbance, energy companies and government agencies erect artificial nesting structures within oil and gas fields. The higher densities of Ferruginous Hawks we found in the Wyoming Basin region may therefore be related to an abundance of anthropogenic nesting structures. Identifying causes of regional differences in density for these species, however, was beyond the scope of this study.
Our estimates of occupied nest density for both species may also have been influenced by nest substrate. Ayers and Anderson (1999) found that nest substrate had the greatest effect on detection probability of Ferruginous Hawks from a fixed-wing plane. Differences in detectability due to nest substrate could have produced unmodeled heterogeneity in our sample; however, small sample sizes limited our ability to include nest substrate as a detection covariate. Ferruginous Hawks nest on a variety of substrates, including hill tops, small rock structures, lone trees, and artificial nest platforms (Bechard et al. 1990, Woffinden and Murphy 1983). We found nests on all of these substrates during our study, which may have contributed to uncertainty in our detection estimates.
Within the northern portion of the Wyoming Basin, known as the Bighorn Basin, we found only one occupied Golden Eagle nest and one occupied Ferruginous Hawk nest in seven sampled townships. The Bighorn Basin is composed largely of salt desert shrub, rather than sagebrush steppe (Chapman et al. 2004), and does not support many mammalian prey species common to the southern portion of the Wyoming Basin, including pygmy rabbits (Brachylagus idahoensis), and Uinta (Urocitellus armatus), Wyoming (Urocitellus elegans), golden-mantled (Callospermophilus lateralis), and thirteen-lined ground squirrels (Ictidomys tridecemlineatus; Knick et al. 2011, Orabona et al. 2012, Thorington et al. 2012). Although white-tailed prairie dogs (Cynomys leucurus) occur in the Bighorn Basin, their abundance is much lower than in the southern portion of the Wyoming Basin (Grenier and Filipi 2009). Although we sampled few townships in this area, the lack of occupied nests found may be related to the scarcity of these prey, which make up a large portion of the diets of both Ferruginous Hawks and Golden Eagles (MacLaren et al. 1988).
The use of a helicopter for a second observer on a subset of transects provided a valuable correction to our nesting pair density estimates using distance sampling. Our regression of fixed-wing plane detections and misses over distance suggested the fixed-wing plane was not adequately detecting all occupied nests directly on the transect, an omission which violates a key assumption of distance modeling. With corrected probability of detection at the transect, however, distance sampling density estimates were very similar to our second set of density estimates using double-observer detection probability. This agreement between results of the two methods provides confidence in our estimates of nesting pair density for both species. Correction for imperfect detection at the transect line can be accomplished with a variety of methods (Alpizar-Jara and Pollock 1996, Quang and Becker 1997, Borchers et al. 1998, Laake et al. 2008). This type of correction is frequently used to improve distance sampling estimates in marine mammal research (Borchers et al. 2006, Pollock et al. 2006) and increasingly in a variety of terrestrial taxa (Kissling and Garton 2006, Buckland et al. 2007, Fewster and Pople 2008). Based on our results, we recommend inclusion of one of these types of correction to model detection at the transect when surveying in fixed-wing planes.
The efficiency of the observers in the helicopter at detecting occupied nests was insufficient to justify the added cost or time of helicopter surveys for estimating density in our study. There was a small difference in detection of occupied Ferruginous Hawk nests between observers in the fixed-wing plane and those in the helicopter, but detection probabilities were almost equal for Golden Eagles. This could be due to the larger size and conspicuousness of Golden Eagles and their occupied nests in the sagebrush steppe and grassland landscape, as both survey modes had a high probability of detecting an occupied Golden Eagle nest. Occupied Ferruginous Hawk nests were more likely to be located by the helicopter, but not enough to make up for the disparity in cost and time between the two modes. Thus, if time to completion or cost are limiting aspects of a survey, a fixed-wing plane will better meet these needs with little loss of detection ability. The results of our comparison, however, are based on linear transects over relatively smooth terrain, which may reduce the relative advantage of using a helicopter. In rugged terrain or with an active search, helicopters may be superior in safety and efficiency due to their greater maneuverability and visibility. Additionally, the efficiency and cost of both survey modes will likely be influenced by several outside factors, including existing occupied nest density, experience of observers, and habitat type. Thus, the decision to use a fixed-wing plane or helicopter should be determined not only by cost or detection ability, but by considerations of survey objectives, visibility and maneuverability requirements, and safety concerns.
Acknowledgments
Primary research support was provided by the U.S.D.A. Rocky Mountain Research Station, Bureau of Land Management, and Wyoming Game and Fish Department. Additional funding sources include Wyoming Wildlife Heritage Foundation; PacifiCorp; Pathfinder Renewable Wind Energy, LLC; and the Wyoming Governor’s Office. Elizabeth Furtaw provided essential GIS support. Pilots Dave Savage (Savage Air Service, Inc.), Ken Overfield (Sky Aviation, Corp.), Dave Stinson (Sky Aviation, Corp.), and Neil Cadwell (Elkhorn Aviation, LLC) were pivotal in project completion. Bob Oakleaf served as primary observer in the helicopter, Zach Wallace and Matt Wilson as primary observers for fixed-wing aerial surveys of transects and leaders of the field crews, which included Roderick Lawrence, Jon Lipke, and Ray Buchheit.