Open Access
How to translate text using browser tools
1 December 2015 Nesting Pair Density and Abundance of Ferruginous Hawks (Buteo regalis) and Golden Eagles (Aquila chrysaetos) from Aerial Surveys in Wyoming
Lucretia E. Olson
Author Affiliations +
Abstract

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.

i0892-1016-49-4-400-t03.tif

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).

Figure 1. 

Study area in Wyoming, U.S.A. Detail shows our study area, comprising all townships within the distribution of Ferruginous Hawks in Wyoming, color-coded by ecoregion. Black squares indicate townships randomly selected for sampling, and outlined gray regions are mountainous areas not included in our study.

i0892-1016-49-4-400-f01.tif

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.

Figure 2. 

Diagram of the aerial survey flight lines in a survey township in Wyoming. Each transect is 9.65 km long and 600 m apart.

i0892-1016-49-4-400-f02.tif

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):

i0892-1016-49-4-400-e01.tif

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:

i0892-1016-49-4-400-e02.tif

where n is the number of occupied nests detected per township, L is township length, W is township width, and 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).

Figure 3. 

Occupied Ferruginous Hawk and Golden Eagle nests detected during surveys in 2010 and 2011 in shrubland and grassland habitat in Wyoming, U.S.A.

i0892-1016-49-4-400-f03.tif

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.

i0892-1016-49-4-400-t01.tif

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.

i0892-1016-49-4-400-t02.tif

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.

Literature Cited

1.

H. Akaike 1973. Information theory as an extension of the maximum likelihood principle. Pages 267–281 in B.N. Petrovand F. Caski [Eds.], Second international symposium on information theory. Akademiai Kiado, Japan. Google Scholar

2.

R. Alpizar-Jaraand K.H. Pollock 1996. A combination line transect and capture-recapture sampling model for multiple observers in aerial surveys. Environmental and Ecological Statistics 3:311–327. Google Scholar

3.

D.E. Andersen 2007. Survey techniques. Pages 89–100 in D. Birdand K. Bildstein [Eds.], Raptor research and management techniques. Hancock House Publishers, Blaine, WA U.S.A. Google Scholar

4.

D.E. Andersen O.J. Rongstadand W.R. Mytton 1985. Line transect analysis of raptor abundance along roads. Wildlife Society Bulletin 13:533–539. Google Scholar

5.

L.W. Ayersand S.H. Anderson 1999. An aerial sightability model for estimating Ferruginous Hawk population size. Journal of Wildlife Management 63:85–97. Google Scholar

6.

L.W. Ayers R.J. Oakleafand T. Filipi 2009. Distribution and abundance of breeding Ferruginous Hawk pairs in southcentral Wyoming, Completion report. Pages 72–85 in A. Orabona [Ed.], Threatened, endangered, and nongame bird and mammal investigations. Wyoming Game and Fish Department, Laramie, WY U.S.A. Google Scholar

7.

M.J. Bechard R.L. Knight D.G. Smithand R.E. Fitzner 1990. Nest sites and habitats of sympatric hawks (Buteo spp.) in Washington. Journal of Field Ornithology 61:159–170. Google Scholar

8.

T.L. Booms P.F. Schempf B.J. McCaffery M.S. Lindbergand M.R. Fuller 2010. Detection probability of cliff-nesting raptors during helicopter and fixed-wing aircraft surveys in western Alaska. Journal of Raptor Research 44:175–187. Google Scholar

9.

D.L. Borchers W. Zucchiniand R.M. Fewster 1998. Mark-recapture models for line transect surveys. Biometrics 54:1207–1220. Google Scholar

10.

D.L. Borchers J.L. Laake C. Southwelland C.G.M. Paxton 2006. Accommodating unmodeled heterogeneity in double‐observer distance sampling surveys. Biometrics 62:372–378. Google Scholar

11.

L.A. Brennanand W.P. Kuvlesky 2005. North American grassland birds: an unfolding conservation crisis? Journal of Wildlife Management 69:1–13. Google Scholar

12.

S.T. Buckland D.R. Anderson K.P. Burnhamand J.L. Laake 1993. Distance sampling: estimating abundance of biological populations. Chapman and Hall, London, U.K. Google Scholar

13.

S.T. Buckland D.L. Borchers A. Johnston P.A. Henrysand T.A. Marques 2007. Line transect methods for plant surveys. Biometrics 63:989–998. Google Scholar

14.

K.P. Burnhamand D.R. Anderson 2002. Model selection and multimodel inference: a practical information-theoretic approach. Springer-Verlag, New York, NY U.S.A. Google Scholar

15.

F.J. Camenzind 1969. Nesting ecology and behavior of the Golden Eagle Aquila chrysaetos L.Brigham Young University Science Bulletin, Biological Series 10:4–15. Google Scholar

16.

S.S. Chapman S.A. Bryce J.M. Omernik D.G. Despain J. Zumbergeand M. Conrad 2004. Ecoregions of Wyoming (color poster with map, descriptive text, summary tables, and photographs). U.S. Geological Survey, Reston, VA U.S.A. Google Scholar

17.

C.J. Conwayand J.C. Simon 2003. Comparison of detection probability associated with Burrowing Owl survey methods. Journal of Wildlife Management 67:501–511. Google Scholar

18.

H.E. Copeland A. Pocewiczand J.M. Kiesecker 2011. Geography of energy development in western North America: potential impacts on terrestrial ecosystems. Pages 7–22 in D.E. Naugle [Ed.], Energy development and wildlife conservation in western North America. Island Press/Center for Resource Economics, Washington, DC U.S.A. Google Scholar

19.

K.W. Davies C.S. Boyd J.L. Beck J.D. Bates T.J. Svejcarand M.A. Gregg 2011. Saving the sagebrush sea: an ecosystem conservation plan for big sagebrush plant communities. Biological Conservation 144:2573–2584. Google Scholar

20.

J.A. Dechant M.L. Sondreal D.H. Johnson L.D. Igl C.M. Goldade A.L. Zimmermanand B.R. Euliss 2002. Effects of management practices on grassland birds: Ferruginous Hawk. Northern Prairie Wildlife Research Center, Jamestown, ND U.S.A. Google Scholar

21.

J.P. DeLong 2004. Effects of management practices on grassland birds: Golden Eagle. Northern Prairie Wildlife Research Center, Jamestown, ND U.S.A. Google Scholar

22.

B.C. Fedyand K.E. Doherty 2011. Population cycles are highly correlated over long time series and large spatial scales in two unrelated species: Greater Sage-Grouse and cottontail rabbits. Oecologia 165:915–924. Google Scholar

23.

R.M. Fewsterand A.R. Pople 2008. A comparison of mark–recapture distance-sampling methods applied to aerial surveys of eastern grey kangaroos. Wildlife Research 35:320–330. Google Scholar

24.

M.R. Fullerand J.A. Mosher 1987. Raptor survey techniques. Pages 37–65 in B.A. Giron Pendleton B.A. Millsap K.W. Clineand D.M. Bird [Eds.], Raptor management techniques manual. National Wildlife Federation, Washington, DC U.S.A. Google Scholar

25.

D.S. Gilmerand R.E. Stewart 1983. Ferruginous Hawk populations and habitat use in North Dakota. Journal of Wildlife Management 47:146–157. Google Scholar

26.

R.E. Good R.M. Nielson H.H. Sawyerand L.L. McDonald 2007. A population estimate for Golden Eagles in the western United States. Journal of Wildlife Management 71:395–402. Google Scholar

27.

M. Grenierand T. Filipi 2009. Estimate of abundance for the white-tailed prairie dog in Wyoming. Pages 149–172 in A. Orabona [Ed.], Threatened, endangered, and nongame bird and mammal investigations. Wyoming Game and Fish Department, Lander, WY U.S.A. Google Scholar

28.

J.E. Gross L.C. Stoddartand F.H. Wagner 1974. Demographic analysis of a northern Utah jackrabbit population. Wildlife Society Monographs 40:3–68. Google Scholar

29.

J.M. Hoekstra T.M. Boucher T.H. Rickettsand C. Roberts 2005. Confronting a biome crisis: global disparities of habitat loss and protection. Ecology Letters 8:23–29. Google Scholar

30.

C.S. Houstonand M.J. Bechard 1984. Decline of the Ferruginous Hawk in Saskatchewan. American Birds 38:166–170. Google Scholar

31.

R.P. Howardand M.L. Wolfe 1976. Range improvement practices and Ferruginous Hawks. Journal of Range Management 29:33–37. Google Scholar

32.

W.G. Hunt 1998. Raptor floaters at Moffat’s equilibrium. Oikos 82:191–197. Google Scholar

33.

D.A. Keinath M.D. Andersenand G.P. Beauvais 2010. Range and modeled distribution of Wyoming’s species of greatest conservation need. Report prepared by the Wyoming Natural Diversity Database, Laramie, WY, for the Wyoming Game and Fish Department, Cheyenne, WY and the U.S. Geological Survey. Fort Collins, CO U.S.A. Google Scholar

34.

M.L. Kisslingand E.O. Garton 2006. Estimating detection probability and density from point-count surveys: a combination of distance and double-observer sampling. Auk 123:735–752. Google Scholar

35.

S.T. Knick D.S. Dobkin J.T. Rotenberry M.A. Schroeder W.M. Vander Haegenand C. Van Riper, III 2003. Teetering on the edge or too late? Conservation and research issues for avifauna of sagebrush habitats. Condor 105:611–634. Google Scholar

36.

S.T. Knick S.E. Hanser M. Leu C.L. Aldridgeand M.J. Wisdom 2011. Introduction: an ecoregional assessment of the Wyoming Basins. Pages 1–9 in S.E. Hanser M. Leu S.T. Knickand C.L. Aldridge [Eds.], Sagebrush ecosystem conservation and management: ecoregional assessment tools and models for the Wyoming Basin. Allen Press, Lawrence, KS U.S.A. Google Scholar

37.

D.H. Knight 1996. Mountains and plains: the ecology of Wyoming landscapes. Yale Univ. Press, New Haven, CT U.S.A. Google Scholar

38.

M.N. Kochert 1972. Population status and chemical contamination in Golden Eagles in southwestern Idaho. M.S. thesis, Univ. Idaho, Moscow, ID U.S.A. Google Scholar

39.

M.N. Kochertand K. Steenhof 2002. Golden Eagles in the US and Canada: status, trends, and conservation challenges. Journal of Raptor Research 36:32–40. Google Scholar

40.

M.N. Kochert K. Steenhof C.L. McIntyreand E.H. Craig 2002. Golden Eagle (Aquila chrysaetos). in A. Poole [Ed.], The birds of North America online, No. 684. Cornell Lab of Ornithology, Ithaca, NY U.S.A.  http://bna.birds.cornell.edu/bna/species/684 (last accessed 22 May 2014). Google Scholar

41.

J. Laake M.J. Dawsonand J. Hone 2008. Visibility bias in aerial survey: mark-recapture, line-transect or both? Wildlife Research 35:299–309. Google Scholar

42.

R.N. Lehman P.L. Kennedyand J.A. Savidge 2007. The state of the art in raptor electrocution research: a global review. Biological Conservation 136:159–174. Google Scholar

43.

J.T. Lokemoenand H.F. Duebbert 1976. Ferruginous Hawk nesting ecology and raptor populations in northern South Dakota. Condor 78:464–470. Google Scholar

44.

P.A. MacLaren S.H. Andersonand D.E. Runde 1988. Food habits and nest characteristics of breeding raptors in southwestern Wyoming. Western North American Naturalist 48:548–553. Google Scholar

45.

B.F. Manly 2007. Randomization, bootstrap and Monte Carlo methods in biology, Third Ed. Chapman and Hall, Boca Raton, FL U.S.A. Google Scholar

46.

B.A. Millsap T.G. Grubb R.K. Murphy T. Swemand J.W. Watson 2015. Conservation significance of alternative nests of Golden Eagles. Global Ecology and Conservation 3:234–241. Google Scholar

47.

B.A. Millsapand M.N. LeFranc Jr 1988. Road transect counts for raptors: how reliable are they? Journal of Raptor Research 22:8–16. Google Scholar

48.

J.D. Nichols J.E. Hines J.R. Sauer F.W. Fallon J.E. Fallonand P.J. Heglund 2000. A double-observer approach for estimating detection probability and abundance from point counts. Auk 117:393–408. Google Scholar

49.

R.F. Noss E.T. LaRoe, IIIand J.M. Scott 1995. Endangered ecosystems of the United States: a preliminary assessment of loss and degradation. Biological Report 28. U.S.D.I. National Biological Service, Washington, DC U.S.A. Google Scholar

50.

R.J. Oakleaf L.E. Olson J.R. Squiresand Z.P. Wallace 2014. The status of Golden Eagles (Aquila chrysaetos) in Wyoming: a preliminary review. Pages 357–396 in A.C. Orabonaand N. Cudworth [Eds.], Threatened, endangered, and nongame bird and mammal investigations. Annual completion report. Wyoming Game and Fish Department, Lander, WY U.S.A. Google Scholar

51.

R.R. Olendorff 1993. Status, biology, and management of Ferruginous Hawks: a review. Raptor Research and Technical Assistance Center, Spec. Rep. U.S.D.I. Bureau of Land Management, Boise, ID U.SA. Google Scholar

52.

A.C. Orabona C. Rudd M. Grenier Z. Walker S. Patlaand R. Oakleaf 2012. Atlas of birds, mammals, amphibians, and reptiles in Wyoming. Wyoming Game and Fish Department, Nongame Program, Lander, WY U.S.A. Google Scholar

53.

R.L. Phillips T.P. McEneaneyand A.E. Beske 1984. Population densities of breeding Golden Eagles in Wyoming. Wildlife Society Bulletin 12:269–273. Google Scholar

54.

K.H. Pollock H.D. Marsh I.R. Lawlerand M.W. Alldredge 2006. Estimating animal abundance in heterogeneous environments: an application to aerial surveys for dugongs. Journal of Wildlife Management 70:255–262. Google Scholar

55.

K.H. Pollock J.D. Nichols T.R. Simons G.L. Farnsworth L.L. Baileyand J.R. Sauer 2002. Large scale wildlife monitoring studies: statistical methods for design and analysis. Environmetrics 13:105–119. Google Scholar

56.

P.X. Quangand E.F. Becker 1997. Combining line transect and double count sampling techniques for aerial surveys. Journal of Agricultural, Biological, and Environmental Statistics 2:230–242. Google Scholar

57.

R Core Team. 2013. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.  http://www.R-project.org/(last accessed 20 May 2013). Google Scholar

58.

S.D. Roth Jrand J.M. Marzluff 1989. Nest placement and productivity of Ferruginous Hawks in western Kansas. Transactions of the Kansas Academy of Science 92:132–148. Google Scholar

59.

J.K. Schmutz 1984. Ferruginous and Swainson’s Hawk abundance and distribution in relation to land use in southeastern Alberta. Journal of Wildlife Management 48:1180–1187. Google Scholar

60.

J.K. Schmutz S.H. Brechtel K.D. DeSmet D.G. Hjertaas G.L. Holroyd C.S. Houstonand R.W. Nero 1992. Recovery plan for the Ferruginous Hawk in Canada. Recovery of National Endangered Wildlife (RENEW). Ottawa, Canada. Google Scholar

61.

J.K. Schmutz D.T.T. Flockhart C.S. Houstonand P.D. McLoughlin 2008. Demography of Ferruginous Hawks breeding in western Canada. Journal of Wildlife Management 72:1352–1360. Google Scholar

62.

J.K. Schmutzand D.J. Hungle 1989. Populations of Ferruginous and Swainson’s Hawks increase in synchrony with ground squirrels. Canadian Journal of Zoology 67:2596–2601. Google Scholar

63.

D.G. Smithand J.R. Murphy 1973. Breeding ecology of raptors in the eastern Great Basin of Utah. Brigham Young University Science Bulletin-Biological Series 18:1–76. Google Scholar

64.

K. Steenhofand I. Newton 2007. Assessing nesting success and productivity. Pages 181–190 in D.M. Birdand K.L. Bildstein [Eds.], Raptor research and management techniques. Raptor Research Foundation. Hancock House Publishers, Blaine, WA U.S.A. Google Scholar

65.

G.W. Suterand J.L. Joness 1981. Criteria for Golden Eagle, Ferruginous Hawk, and Prairie Falcon nest site protection. Raptor Research 15:12–18. Google Scholar

66.

L. Thomas S.T. Buckland E.A. Rexstad J.L. Laake S. Strindberg S.L. Hedley J.R.B. Bishop T.A. Marquesand K.P. Burnham 2010. Distance software: design and analysis of distance sampling surveys for estimating population size. Journal of Applied Ecology 47:5–14. Google Scholar

67.

R.W. Thorington J.L. Koprowski M.A. Steeleand J.F. Whatton 2012. Squirrels of the world. The Johns Hopkins Univ. Press, Baltimore, MD U.S.A. Google Scholar

68.

United States Energy Information Administration. 2014. Number of producing gas wells. U.S. Department of Energy, Washington, DC U.S.A.  http://www.eia.gov (last accessed March 2014). Google Scholar

69.

B. Van Horne 1983. Density as a misleading indicator of habitat quality. Journal of Wildlife Management 47:893–901. Google Scholar

70.

C.M. Whiteand T.L. Thurow 1985. Reproduction of Ferruginous Hawks exposed to controlled disturbance. Condor 87:14–22. Google Scholar

71.

Wyoming Oil and Gas Conservation Commission. 2014. Well header data 052013 DB5.  http://wogcc.state.wy.us/(last accessed 30 May 2013). Google Scholar

72.

N.D. Woffindenand J.R. Murphy 1977. Population dynamics of the Ferruginous Hawk during a prey decline. Western North American Naturalist 37:411–425. Google Scholar

73.

N.D. Woffindenand J.R. Murphy 1983. Ferruginous Hawk nest site selection. Journal of Wildlife Management 47:216–219. Google Scholar

74.

N.D. Woffindenand J.R. Murphy 1989. Decline of a Ferruginous Hawk population: a 20-year summary. Journal of Wildlife Management 53:1127–1132. Google Scholar
© 2015 The Raptor Research Foundation, Inc.
Lucretia E. Olson "Nesting Pair Density and Abundance of Ferruginous Hawks (Buteo regalis) and Golden Eagles (Aquila chrysaetos) from Aerial Surveys in Wyoming," Journal of Raptor Research 49(4), 400-412, (1 December 2015). https://doi.org/10.3356/rapt-49-04-400-412.1
Received: 18 June 2014; Accepted: 1 February 2015; Published: 1 December 2015
KEYWORDS
Aquila chrysaetos
Buteo regalis
distance sampling
double-observer sampling
ferruginous hawk
golden eagle
grassland Ecosystems
Back to Top