Resource Selection by the California Condor (Gymnogyps californianus) Relative to Terrestrial-Based Habitats and Meteorological Conditions James W. Rivers 1 *, J. Matthew Johnson 2¤ , Susan M. Haig 2 , Carl J. Schwarz 3 , John W. Glendening 4 , L. Joseph Burnett 5 , Daniel George 6 , Jesse Grantham 7 1 Department of Forest Ecosystems and Society, Oregon State University, Corvallis, Oregon, United States of America, 2 U.S. Geological Survey Forest and Rangeland Ecosystem Science Center, Corvallis, Oregon, United States of America, 3 Department of Statistics and Actuarial Science, Simon Fraser University, Burnaby, British Columbia, Canada, 4 Salinas, California, United States of America, 5 Ventana Wildlife Society, Salinas, California, United States of America, 6 Pinnacles National Park, National Park Service, Paicines, California, United States of America, 7 California Condor Recovery Program, U.S. Fish and Wildlife Service, Ventura, California, United States of America Abstract Condors and vultures are distinct from most other terrestrial birds because they use extensive soaring flight for their daily movements. Therefore, assessing resource selection by these avian scavengers requires quantifying the availability of terrestrial-based habitats, as well as meteorological variables that influence atmospheric conditions necessary for soaring. In this study, we undertook the first quantitative assessment of habitat- and meteorological-based resource selection in the endangered California condor (Gymnogyps californianus) within its California range and across the annual cycle. We found that condor use of terrestrial areas did not change markedly within the annual cycle, and that condor use was greatest for habitats where food resources and potential predators could be detected and where terrain was amenable for taking off from the ground in flight (e.g., sparse habitats, coastal areas). Condors originating from different release sites differed in their use of habitat, but this was likely due in part to variation in habitats surrounding release sites. Meteorological conditions were linked to condor use of ecological subregions, with thermal height, thermal velocity, and wind speed having both positive (selection) and negative (avoidance) effects on condor use in different areas. We found little evidence of systematic effects between individual characteristics (i.e., sex, age, breeding status) or components of the species management program (i.e., release site, rearing method) relative to meteorological conditions. Our findings indicate that habitat type and meteorological conditions can interact in complex ways to influence condor resource selection across landscapes, which is noteworthy given the extent of anthropogenic stressors that may impact condor populations (e.g., lead poisoning, wind energy development). Additional studies will be valuable to assess small-scale condor movements in light of these stressors to help minimize their risk to this critically endangered species. Citation: Rivers JW, Johnson JM, Haig SM, Schwarz CJ, Glendening JW, et al. (2014) Resource Selection by the California Condor (Gymnogyps californianus) Relative to Terrestrial-Based Habitats and Meteorological Conditions. PLoS ONE 9(2): e88430. doi:10.1371/journal.pone.0088430 Editor: Fabio S. Nascimento, Universidade de Sa ˜o Paulo, Faculdade de Filosofia Cie ˆncias e Letras de Ribeira ˜o Preto, Brazil Received June 14, 2013; Accepted January 8, 2014; Published February 11, 2014 This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication. Funding: Funding came from the National Fish and Wildlife Foundation, Wildlife Encounters, Inc., William R. Hearst III, The Oregon Zoo Foundation, Pacific Gas and Electric, and J. Glendening. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors note that they received funding from Wildlife Encounters, Inc. but this does not alter their adherence to all the PLOS ONE policies on sharing data and materials. * E-mail: [email protected]¤ Current address: U.S. Forest Service, Plumas National Forest, Quincy, California, United States of America Introduction All animals require resources critical to their survival, and determining how and why an organism selects among available resources is fundamental to understanding its ecological niche. Resource selection is considered to occur across sequential spatial scales, with the broadest level of selection being that of a geographic range (i.e., first-order selection), followed by selection of individual home ranges (i.e., second-order selection), selection of coarse-scale habitats within the home range (i.e., third-order selection), and selection of microhabitats within coarse-scale habitats (i.e., fourth-order selection; [9,24,33]).Traditionally, studies of vertebrate resource selection have focused on quantify- ing use of terrestrial resources because most terrestrial-based species have their needs met by resources that are located on or near ground level. Nevertheless, a number of animals have resource requirements that extend beyond terrestrial habitats, especially for organisms that use extensive, long-ranging flight [14]. Animals that fly above the earth’s surface do so within the convective ‘‘boundary layer’’ of the atmosphere, and meteorolog- ical conditions that occur within this stratum can strongly influence space use and movements for some taxa. Animal flight can be facilitated by two types of vertical air movement within the boundary layer: thermal lift and orographic lift [5,6,32,42,44,57]. Thermal lift occurs when solar radiation heats the earth’s surface and creates convective thermals of vertically rising air within the boundary layer [18]. Warm, rising air provides lift used by large birds to move vertically via soaring flight within convective PLOS ONE | www.plosone.org 1 February 2014 | Volume 9 | Issue 2 | e88430
13
Embed
Resource Selection by the California Condor (Gymnogyps ...people.forestry.oregonstate.edu/jim-rivers/sites... · Resource Selection by the California Condor (Gymnogyps californianus)
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Resource Selection by the California Condor (Gymnogypscalifornianus) Relative to Terrestrial-Based Habitats andMeteorological ConditionsJames W. Rivers1*, J. Matthew Johnson2¤, Susan M. Haig2, Carl J. Schwarz3, John W. Glendening4,
L. Joseph Burnett5, Daniel George6, Jesse Grantham7
1 Department of Forest Ecosystems and Society, Oregon State University, Corvallis, Oregon, United States of America, 2 U.S. Geological Survey Forest and Rangeland
Ecosystem Science Center, Corvallis, Oregon, United States of America, 3 Department of Statistics and Actuarial Science, Simon Fraser University, Burnaby, British
Columbia, Canada, 4 Salinas, California, United States of America, 5 Ventana Wildlife Society, Salinas, California, United States of America, 6 Pinnacles National Park,
National Park Service, Paicines, California, United States of America, 7 California Condor Recovery Program, U.S. Fish and Wildlife Service, Ventura, California, United States
of America
Abstract
Condors and vultures are distinct from most other terrestrial birds because they use extensive soaring flight for their dailymovements. Therefore, assessing resource selection by these avian scavengers requires quantifying the availability ofterrestrial-based habitats, as well as meteorological variables that influence atmospheric conditions necessary for soaring. Inthis study, we undertook the first quantitative assessment of habitat- and meteorological-based resource selection in theendangered California condor (Gymnogyps californianus) within its California range and across the annual cycle. We foundthat condor use of terrestrial areas did not change markedly within the annual cycle, and that condor use was greatest forhabitats where food resources and potential predators could be detected and where terrain was amenable for taking offfrom the ground in flight (e.g., sparse habitats, coastal areas). Condors originating from different release sites differed intheir use of habitat, but this was likely due in part to variation in habitats surrounding release sites. Meteorologicalconditions were linked to condor use of ecological subregions, with thermal height, thermal velocity, and wind speedhaving both positive (selection) and negative (avoidance) effects on condor use in different areas. We found little evidenceof systematic effects between individual characteristics (i.e., sex, age, breeding status) or components of the speciesmanagement program (i.e., release site, rearing method) relative to meteorological conditions. Our findings indicate thathabitat type and meteorological conditions can interact in complex ways to influence condor resource selection acrosslandscapes, which is noteworthy given the extent of anthropogenic stressors that may impact condor populations (e.g., leadpoisoning, wind energy development). Additional studies will be valuable to assess small-scale condor movements in lightof these stressors to help minimize their risk to this critically endangered species.
Citation: Rivers JW, Johnson JM, Haig SM, Schwarz CJ, Glendening JW, et al. (2014) Resource Selection by the California Condor (Gymnogyps californianus)Relative to Terrestrial-Based Habitats and Meteorological Conditions. PLoS ONE 9(2): e88430. doi:10.1371/journal.pone.0088430
Editor: Fabio S. Nascimento, Universidade de Sao Paulo, Faculdade de Filosofia Ciencias e Letras de Ribeirao Preto, Brazil
Received June 14, 2013; Accepted January 8, 2014; Published February 11, 2014
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone forany lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: Funding came from the National Fish and Wildlife Foundation, Wildlife Encounters, Inc., William R. Hearst III, The Oregon Zoo Foundation, Pacific Gasand Electric, and J. Glendening. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors note that they received funding from Wildlife Encounters, Inc. but this does not alter their adherence to all the PLOS ONEpolicies on sharing data and materials.
throughout the course of the study. GPS transmitters were
programmed to collect location data at hourly intervals each day
from 0500–2000 h PST with a resolution of 618 m (based on
manufacturer specifications). Although we attempted to assign
transmitters to sex and age classes in a balanced manner,
management needs required the non-random assignment of
transmitters in some cases. The number of individuals that were
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 2 February 2014 | Volume 9 | Issue 2 | e88430
fitted with GPS transmitters began with two individuals in 2003
and increased to 50 individuals in 2010, reflecting the growth of
the free-flying condor population in California. Because lead
exposure is a serious threat to free-ranging condors [17,47] we
regularly re-captured condors to collect blood and evaluate lead
levels; individuals with high levels were retained in captivity for
varying lengths of time. This created gaps in location data for
some individuals and led us to restrict our analysis of resource use
to individuals that had a minimum of 100 locations per month to
ensure adequate sampling. We used month as the temporal scale
for our analysis because this time period provided an opportunity
to assess potential changes in resource use across the annual cycle
and enabled us to incorporate meteorological data, which were
available in monthly summaries from January 2007–December
2009 (see below).
Figure 1. Habitat classifications of landcover in California. The range of primary concern according to the 1984 California Condor RecoveryPlan [60] is shown by a solid black outline with condor release sites considered in this study (crosses). Note that Bitter Creek NWR and HopperMountain NWR are combined for analysis because they are both part of the Hopper Mountain National Wildlife Refuge Complex.doi:10.1371/journal.pone.0088430.g001
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 3 February 2014 | Volume 9 | Issue 2 | e88430
Utilization Distributions and Home Range DelineationWe used location data to calculate utilization distributions (UDs)
on a monthly basis, using individual condors as our sampling unit.
UDs quantify relative space use of an individual by predicting its
probability of occurrence as a function of relocation points within
an area of interest, such as a monthly home range [63]. Thus, an
animal’s use of space can be envisioned as a three dimensional
plot, with the height of the UD at any particular location within an
area of interest (e.g., monthly home range) being a function of the
likelihood an animal will use that location, with a greater UD
height indicating a greater likelihood of use [35]. To calculate
UDs, we used 99% fixed kernel density analysis [22,54]. Initially,
we evaluated the reference (href) and the least squares cross-
validation (hlscv) smoothing parameters [25,64] for use with our
dataset. However, both smoothing parameters failed during initial
analysis (see [48]) so we calculated utilization distributions by
initially estimating href with the Home Range Tools in ArcGIS
[49], and then used Hawth’s Analysis Tools for ArcGIS [4] to
calculate 99% fixed-kernel monthly home ranges using a grid cell
size of 100 m. We used an ad hoc smoothing parameter (hadhoc) to
choose the smallest increment of href that resulted in a contiguous
99% kernel polygon (i.e., 0.3*href = hadhoc) as this minimizes
overestimation of the outward boundary of the utilization
distribution [3,12]. We deemed this a reasonable approach
because [1] field observations found that using the href overesti-
mated condor space use in habitats immediately surrounding high
use areas that contained many overlapping location points, and [2]
condors often concentrated their perching and roosting at the
same distinct locations throughout the annual cycle, leading to
location data that contained a large number of overlapping
individual location points.
Terrestrial-based Habitat and Meteorological DataWe used two independent datasets to quantify the terrestrial-
based habitat and meteorological conditions available to free-
ranging condors. In the first, we obtained Geographic Information
Systems (GIS) landcover data from the U.S. Department of
Agriculture that spanned the geographic range of condors in
California and covered 47 ecological subregions (hereafter,
ecoregions) as delineated by [11]. Next, we took 244 distinct
landcover classifications from [13] for the state of California and
reclassified them into 12 distinct habitat types (Document S1) to
eliminate redundancy in landcover classifications and to express
landcover in units that were ecologically relevant to condors. We
then overlaid habitat data on ecoregion data so that we could
assess resource selection on two nested spatial scales (habitat was
nested within ecoregion). We used this approach because wind
resource data were available at a coarser spatial scale than habitat
data, requiring analysis on the ecoregion scale (see below).
In the second dataset, we obtained output from the North
American Mesoscale (NAM) weather prognostication model run
by NOAA’s National Centers for Environmental Prediction
(NCEP). This model forecasts atmospheric properties over the
United States at a resolution of 12 km [23]. This model solves the
atmospheric primitive differential equations, in particular incor-
porating surface factors such as solar heating, soil temperature, soil
moisture, and vegetative type which vary over the 12 km NAM
grid. Because this output did not directly provide the parameters
most relevant to soaring flight, we used it to estimate three types of
wind resources: (1) maximum thermal height, equivalent to the
boundary layer height under convective conditions [18], [2]
convective thermal velocity, which depends upon the surface
heating rate and boundary layer depth [18], and (3) wind speed,
which was averaged vertically over the boundary layer depth.
Values were computed near mid-day (21Z) and then averaged into
monthly values at each point in the 12 km NAM grid.
Meteorological data were only available for the years 2007–
2009, so we used data from these years to calculate an average
monthly value for each 12612 km cell for each of the three wind
Quantifying Condor Resource SelectionWe used UDs to predict the probability of occurrence within
each habitat type within each monthly home range, which served
as proportional, habitat-specific measures of resource use. We then
compared habitat use measures to the proportion that each habitat
was available for monthly home range estimates by calculating
resource selection ratios [i.e., ln(rf)] as follows:
Ybird,year,month,habitat~
ln rfbird,year,month,habitat
� �~
lnProportion Usebird,year,month,habitat
Proportion Availablebird,year,month,habitat
� �
Values of ln(rf).0 indicate a habitat whose use exceeded its
availability (i.e., the habitat is selected), values of ln(rf) = 0 indicate
a habitat used in proportion to its availability, and values of
ln(rf),0 indicate a habitat whose use was less than its availability
(i.e., the habitat is avoided). We used a ln transformation to make
the analysis ‘‘symmetric’’ around 1 (e.g., a selection ratio of K is
the same distance away from 1 as a selection ratio of 2).
We also used UDs to predict the probability of occurrence at the
ecoregion level within the monthly home range as a measure of
use, compared it to the area of the ecoregion within the monthly
home range as a measure of availability, and related the ln(rf) to
the availability of meteorological parameters in each ecoregion.
Values for meteorological parameters were found by determining
which 12612 km cells overlapped each ecoregion within the
monthly home range and then calculating an average over all
12612 km cells within an ecoregion for each meteorological
parameter. We calculated a selection ratio each month for each
bird’s use of each ecoregion relative to the area of each ecoregion
that fell within that month’s home range:
Ybird,year,month,ecoregion~
ln rfbird,year,month,ecoregion
� �~
lnProportion Usebird,year,month,ecoregion
Proportion Availablebird,year,month,ecoregion
� �
Statistical AnalysisWe performed two levels of statistical analysis. In the first, we
restricted our examination to assess the relative use of a habitat
relative to its availability for each bird-year-month combination
and did not consider meteorological parameters. After calculating
a selection ratio, ln(rf), for each habitat, we then used a repeated-
measures mixed linear modeling approach to model resource use
for each habitat relative to intrinsic characteristics of individuals
(i.e., age, sex, breeding status [two levels each]) and factors related
to the recovery program (i.e., rearing method [two levels], release
location [three levels]). For this analysis, we classified individuals as
either immature (i.e., juvenile [0–2 years], sub-adult [3–5 years])
or adult [$6 years]) because breeding does not occur regularly
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 4 February 2014 | Volume 9 | Issue 2 | e88430
until 6 years of age [56]. We classified adults as breeders if they
frequented a nesting site and were found to have laid an egg. We
classified the rearing method of individuals as either raised in the
wild or reared in captivity. We used the PROC MIXED
procedure in SAS/STAT version 9.3 for Windows to account
for the repeated nature of the measurements, with ln(rf) as the
repeated measure; age, sex, breeding status, release site, and
rearing method as fixed categorical effects; and individual bird was
included as a random factor. Initially, we found that an
autoregressive covariance structure outperformed a compound
symmetrical covariance structure so it was retained for all
subsequent models. Because our data were unbalanced, we
calculated least-squares marginal means (LSMEANS) for effect
sizes, with a Tukey-Kramer adjustment for all multiple compar-
isons, and used the Kenward-Rogers method to calculate degrees
of freedom for contrasts and estimates. We note that although
spatial autocorrelation arises through the use of a fixed-kernel
procedure to construct the UD, we used individual condors as our
sampling unit. This allowed us to ignore spatial autocorrelation
between individual locations because previous work has demon-
strated that individual model coefficients are unbiased even when
autocorrelation is present [1,16,27,29,36].
Figure 2. Coarse-scale condor habitat use. Mean selection ratios for each of 12 habitat types assessed across the annual cycle. Gray shaded arearepresents 95% confidence intervals. A = agriculture, B = coastal dune, C = coastal rock, D = deciduous forest, E = evergreen forest, F = grassland,G = modified land, H = savanna, I = shrubland, J = sparse vegetation, K = unsuitable habitat, L = wetland.doi:10.1371/journal.pone.0088430.g002
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 5 February 2014 | Volume 9 | Issue 2 | e88430
In the second analysis, we again computed a relative use of each
ecoregion to its availability for each bird-year-month-ecoregion
combination. We then fit a mixed-effects linear model for each
ecoregion (again using an autoregressive covariance structure to
account for repeated measurements on the same bird over time)
where the ln(rf) was a function of bird characteristics, month (as a
categorical variable), time, and the meteorological parameters. We
used an information-theoretic approach [8] because bird charac-
teristics were not of interest for this part of our analysis. We
created an initial model set containing the effects of bird-specific
Significant P-values (i.e., ,0.01) are highlighted in bold text.doi:10.1371/journal.pone.0088430.t001
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 6 February 2014 | Volume 9 | Issue 2 | e88430
ecoregions there was no evidence of a sex, age, breeding status, or
rearing effect (Table 3). We found models that included at least one
of the meteorological parameters often had the largest weight among
the model set for many ecoregions (Document S5). Most of the model
weight favored models with meteorological parameters, but meteo-
rological parameters were not equally important in all ecoregions
(Table 4, Fig. 3). Finally, changes in meteorological parameters
across the months were linked with changes in usage in some
ecoregions (Fig. 4, Documents S6, S7). For example, in ecoregion 18,
the model averaged coefficient of thermal height (km) of 1.080 (see
Table 4) indicates that for every 1 km of change in the thermal
height, the ln(rf) increased by 1.080, or exp(1.080) = 2.95 times.
Discussion
Resource selection relative to terrestrial-based habitatsWe found that free-ranging individual condors varied signifi-
cantly in the terrestrial habitats they used within monthly home
ranges when we assessed habitat use separately from meteorolog-
ical data. Selection ratios were greatest for coastal dune, deciduous
forest (including oak woodland), and sparse vegetation habitats;
lowest for grassland and savanna habitats; and condors avoided
evergreen forest and shrubland. Previous authors indicated that
grassland and oak (Quercus spp.) savanna comprises historically
important foraging areas for condors [10,26,56], so it was
somewhat surprising that our analysis revealed that condors did
not exhibit strong selection of either grassland or savanna habitats
within their home range. This may be because selection of
grassland and savanna habitats occurs at broader scales than we
examined (i.e., first- or second-order selection; see [24]), or
because previous authors merely described use of habitats and did
not report selection relative to availability. Selection for sparse
vegetation and coastal habitats likely reflects several important
features of condor foraging ecology: habitats where food resources
(i.e., animal carcasses) and potential predators can be detected,
and habitats that have terrain that is amenable for taking off from
Figure 3. Meteorological parameters influenced condor space use. Estimated mean marginal effects and associated 95% CIs for each ofthree monthly-averaged meteorological parameters by ecoregion. A = thermal height (km), B = thermal velocity (m/s), C = wind speed (m/s).doi:10.1371/journal.pone.0088430.g003
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 7 February 2014 | Volume 9 | Issue 2 | e88430
the ground in flight [56]. Strong selection for coastal habitats is
noteworthy for two reasons. First, condor diets during recent times
(i.e., 1993 to 2001) were reported to have shifted away from
marine mammals and towards terrestrial animals, especially
domestic cattle [10]. However, since 1999, condors have been
observed foraging on coastal marine mammals (e.g., gray whale
[Eschrichtius robustus] and California sea lion [Zalophus californianus])
in central California and concern has increased regarding whether
environmental contaminants found in marine mammal tissues
threaten condor populations in this region [7]. Thus, use of coastal
habitats by the condors in our dataset, coupled with foraging
observations of condors on marine mammal resources, suggests a
return to feeding in coastal areas in recent times. Second, coastal
habitats comprise a very small portion of the total landcover in
California (,0.1%) but the high use and selection of these areas by
condors in central California indicates it is an especially important
habitat for this species, probably due to a combination of foraging
resource and onshore winds that facilitate soaring flight.
We found little evidence that intrinsic characteristics of
individuals and factors related to the recovery program, aside
from release site, influenced the use of terrestrial-based habitats.
Although we did detect a significant effect of month, habitat
selection ratios were rather consistent throughout the year for most
habitats, indicating that the month effects are relatively small when
compared to the type of habitat. In contrast, we did find a
significant and consistent effect of release site on habitat use by
condors. Release sites are spatially distinct and exhibit coarse-scale
differences in habitat types (Fig. 1), so these effects are likely
influenced by differences in habitat availability surrounding the
release sites, and they may not be indicative of consistent
differences in habitat selection by condors among the different
release sites.
Figure 4. Condor selection ratios by ecoregion. Mean monthly ln(rf) values for each ecoregion (filled points) and the effect of themeteorological parameters and the range of monthly average values for each meteorological value (vertical bars). Thus, the vertical bars are not errorestimates but instead represent the range of mean selection ratio across the annual cycle for each ecoregion. A = thermal height (km), B = thermalvelocity (m/s), C = wind speed (m/s).doi:10.1371/journal.pone.0088430.g004
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 8 February 2014 | Volume 9 | Issue 2 | e88430
Resource selection relative to meteorological conditionsOur results provide the first quantitative evidence that resource
selection by condors is linked to meteorological parameters in a
manner that is thought to facilitate soaring flight. These effects
were present despite the relatively coarse scale at which weather
data were available and an analysis based on monthly averages of
meteorological parameters, both of which are expected to reduce
our ability to detect significant effects. Thus, our results indicate
that meteorological conditions can have a particularly strong effect
on condor use of some ecoregions. Our data indicate that the
thermal characteristics we measured (i.e. thermal height and
velocity) generally had a stronger influence on selection by condors
among ecoregions than wind speed averaged across the depth of
the boundary layer. Selection for thermal characteristics is
expected for condors, which are similar to vultures and other
large birds in their use of thermals for large-scale movements
[6,41–43,57]. Thermals are exploited by soaring birds because
they allow individuals to minimize energetic output during large-
scale movements [43,53], such as searching for carrion. Despite
the recognized importance of thermals for movement, California
condors are likely to use orthographic lift to facilitate soaring flight
as has been shown in a Peruvian coastal population of the closely-
related Andean condor (Vultur gryphus [37,44]. Orthographic lift is
likely to be used most often in coastal areas (e.g., Big Sur region)
where on-shore wind conditions create lift and facilitate foraging of
marine mammal carrion [7], in addition to mountaintop areas that
experience updrafts that are strong enough to support the weight
of condors (ca. 8.2 kg [56]). Thus, consideration of the diverse
habitats in which we detected condors indicates that all of the
atmospheric properties we measured are likely to be important for
condors throughout their California range.
Our results also indicate that wind resources were linked to
differences in the extent of condor use at the ecoregion level. The
magnitude of model averaged effects for the three meteorological
parameters were often similar among ecoregions, and for many
ecoregions effects did not differ from zero. However, meteorolog-
ical parameters did have a significant influence on condor use for
several geographically distinct ecoregions (e.g., 101, 116, 125, 128;
see Fig. 4, Document S4), and these effects included selection for
and avoidance of ecoregions relative to the atmospheric parameter
being examined. We also detected a link between meteorological
parameters and individual characteristics as well as factors related to
the recovery program; however, the strength of these relationships
varied substantially relative to habitat type and ecoregion, with
month and release site having the strongest effect. Thus, a complex
picture is emerging regarding how individual condors vary in their
use of space relative to meteorological variables, how these
relationships can change across landscape scales, and how they
differ relative to individual and program-related characteristics. It is
worth noting that measures of model-averaged effect size (Figure 3)
Table 2. Maximum and minimum average values for meteorological parameters observed across the annual cycle by ecoregion.
PLOS ONE | www.plosone.org 9 February 2014 | Volume 9 | Issue 2 | e88430
and selection ratios and the degree of selection exhibited by condors
(Figure 4) can vary independently of each other. For example,
consider ecoregion 117 in panels Figure 3C and Figure 4C. In
Figure 3C, the value for ecoregion 117 indicates that there was a
significant effect size for condor selection for wind speed in this
location (i.e., there is less use when wind speed increases). In
Figure 4, the value for ecoregion 117 indicates that the mean range
of selection ratios across the annual cycle was rather broad and
covered zero. Taken together, this indicates that condors appeared
to avoid this ecoregion during times of high winds yet find it suitable
otherwise. These results highlight how meteorological variables are
dynamic throughout the annual cycle within given ecoregions and
indicate such changes can have strong influence on ecoregion use by
condors. Furthermore, they suggest research that considers
decisions made by an individual on small spatial and temporal
scales may be especially useful for furthering our understanding of
resource selection in this species [19,39].
Implications of resource selection for condorconservation
The California condor is noteworthy because it is one of the
most endangered birds in the world and, as demonstrated in this
study, its use of space is influenced by meteorological conditions.
Our analysis indicates that condors use a wide range of terrestrial
habitats in California, and their movements in some areas are
influenced in part by meteorological conditions. These results
therefore add significantly to previous data and observations on
condor movement and spatial ecology [38,48] by demonstrating
that condors are not restricted in their use of any single habitat but
instead use all available coarse-scale habitat types in California.
Lead poisoning is considered to be the most serious threat to
recovery of the condor population at the current time [17,47];
however, the spatial extent of lead availability on the landscape
and its potential to poison condors is currently unknown. In
addition to the lead issue, concern has also increased in recent
years regarding potential impacts of wind energy developments on
condors in California, where the greatest number of free-flying
individuals currently reside [61]. This has arisen in part because of
recent legislation that requires California to increase the amount of
electricity generated from renewable energy resources to 33% of
total sales of electricity by the end of 2020 [46]. Our results are
informative in both of these contexts because they demonstrate
that condors alter their use of habitat types across the annual cycle
and relative to local meteorological conditions, and that condors
Table 3. Summary of main effects tests for resource selection ratios [i.e., ln(rf)] for 25 ecoregions within the geographic range ofthe California condor in California.
Effect
Month Release site Sex Sex*Month Age class Breeder Rearing method
See Document S3 for the names that correspond to ecoregion codes. Significant P-values (i.e., ,0.01) are highlighted in bold text.doi:10.1371/journal.pone.0088430.t003
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 10 February 2014 | Volume 9 | Issue 2 | e88430
use all available habitats at least to some degree during the course
of the year. Thus, additional studies that focus on documenting
condor lead exposure across the landscape will be especially
valuable to understand the spatio-temporal dynamics of lead
exposure and condor poisoning events. In addition, study of finer
scale movements of condors in areas being evaluated for wind
development would inform siting considerations, development of
curtailment measures, design of deterrence devices, and other
measures to reduce the risk of collision of condors with
anthropogenic structures linked to wind energy developments.
Finally, movement studies that go beyond space use and
incorporate a temporal scale will be especially valuable for
understanding decisions condors make when selecting resources.
Such studies have already shown great promise in enhancing our
understanding of the fine-scale decisions made by obligate
scavengers during large-scale movements [6,53]), and should
provide important data that can help with the conservation and
management of this critically endangered species.
Supporting Information
Document S1 Table of 244 distinct landcover classifica-tions taken from [13] for California and reclassified into12 distinct habitat types.
(PDF)
Document S2 Plots of marginal mean ln(rf) by habitatfor sex effects, age and age class effects, release site,rearing method, and breeding status.
(PDF)
Document S3 Ecoregions as delineated by [11] and usedto assess California Condor resource selection relativeto meteorological parameters.
(PDF)
Document S4 Map of 25 ecoregions in which CaliforniaCondors were observed in reasonable numbers toquantify meteorological parameters in the study.
(PDF)
Document S5 Summary of the number of ecoregionmodels and their parameters that had DAIC#3 fromcandidate models.
(PDF)
Document S6 Plots for three meteorological parametersand raw ln(rf) values plotted against months in theannual cycle for each of the 25 California ecoregionsexamined in the study.
(PDF)
Table 4. Model weighted estimates of meteorological parameters for 25 ecoregions in which condors were observed.
Estimates whose 95% confidence intervals did not overlap with zero are highlighted in bold text.doi:10.1371/journal.pone.0088430.t004
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 11 February 2014 | Volume 9 | Issue 2 | e88430
Document S7 Plots for three meteorological parametersand mean ln(rf) values plotted against months in theannual cycle for each of the 25 California ecoregionsexamined in the study.(PDF)
Acknowledgments
Any use of trade, product, or firm names is for descriptive purposes only
and does not imply endorsement by the U.S. Government. We thank
Bureau of Land Management, California Department of Fish and Game;
National Parks Service, Pinnacles National Park; U.S. Fish and Wildlife
Service, Condor Recovery Program; U.S. Fish and Wildlife Service,
Ventura Field Office; U.S. Geological Survey, Forest and Rangeland
Ecosystem Science Center; and the Ventana Wildlife Society for logistical
support. We thank H. Beyer, J. Brandt, P. Haggarty, D. Howard, K.
Jenkins, J. Kern, B. Kertson, B. Massey, J. Marzluff, S. Murphy, A.
Rodgers, N. Snyder, D. Wiens, and S. Wilbur for logistical support and
feedback, and R. Anthony and two anonymous reviewers for helpful
comments on a previous version of the manuscript.
Author Contributions
Conceived and designed the experiments: JMJ SMH. Performed the
25. Kernohan BJ, Gitzen RA, Millspaugh JJ (2001) Analysis of animal space use andmovements. In: Millspaugh JJ, Marzluff JM, editors. Radio tracking and animal
populations San Diego: Academic Press. pp. 125–166.
26. Koford CB (1953) The California Condor. New York: National AudubonSociety. 154 p.
27. Kutner MH, Nachtsheim CJ, Neter J, Li W (2004) Applied linear statistical
models. Boston: McGraw-Hill. 1396 p.28. Kuvlesky JR WP, Brennan LA, Morrison ML, Boydston KK, Ballard BM, et al.
(2007) Wind energy development and wildlife conservation: challenges andopportunities. J Wildl Manage 71: 2487–2498.
29. Liang KY, Zeger SL (1986) Longitudinal data analysis using generalized linear
models. Biometrika 73: 13–22
30. de Lucas M, Ferrer M, Janss GFE (2012) Using wind tunnels to predict birdmortality in wind farms: the case of griffon vultures. PLoS ONE 11: e48092.
31. Madders M, Whitfield DP (2006) Upland raptors and the assessment of wind
Migration path annotation: cross-continental study of migration-flight response
to environmental conditions. Ecol Appl 21: 2258–2268.33. Manly BFJ, McDonald LL, Thomas DL, McDonald TL, Erickson WP (2002)
Resource selection by animals: statistical design and analysis for field studies.Dordrecht: Kluwer Academic Publishers. 221 p.
34. Martinez-Abrain A, Tavecchia G, Regan HM, Jimenez J, Surroca M, et al.
(2012) Effects of wind farms and food scarcity on a large scavenging bird speciesfollowing an epidemic of bovine spongiform encephalopathy. J Appl Ecol 49:
109–117.
35. Marzluff JM, Millspaugh JJ, Hurvitz P, Handcock MS (2004) Relating resourcesto a probabilistic measure of space use: forest fragments and Steller’s jays.
Ecology 85: 1411–1427.
36. McCullagh P, Nelder JA (1989) Generalized linear models. London: Chapmanand Hall. 532 p.
37. McGahan J (1973) Gliding flight of the Andean condor in nature. J Exp Biol 58:
225–237.
38. Meretsky VJ, Snyder NFR (1992) Range use and movements of Californiacondors. Condor 94: 313–335.
39. Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, et al. (2008)A movement ecology paradigm for unifying organismal movement research.
Proc Natl Acad Sci USA 105: 19052–19059.
40. Pennycuick CJ (1969) The mechanics of bird migration. Ibis 111: 525–556.
41. Pennycuick CJ (1972) Soaring behavior and performance of some East Africanbirds, observed from a motor glider. Ibis 114: 178–218.
42. Pennycuick JC (1998) Field observations of thermals and thermal streets, and the
theory of cross-country soaring flight. J Avian Biol 29: 33–43.43. Pennycuick CJ (2008) Modeling the flying bird. London: Elsevier. 480 p.
gryphs and turkey vultures Cathartes aura around the Paracas Peninsula, Peru. Ibis126: 253–256.
45. Piorkowski MD, Farnsworth AJ, Fry M, Rohrbaugh RW, Fitzpatrick JW, et al.
(2012) Research priorities for wind energy and migratory wildlife. J WildlManage 76: 451–456.
46. Renewable Energy Action Team (2010) Best management practices and
guidance manual: desert renewable energy projects. California EnergyCommission. 265 pp.
47. Rideout BA, Stalis I, Papendick R, Pessier A, Puschner B, et al. (2012) Patterns
of mortality in free-ranging California condors (Gymnogyps californianus). J WildlDis 48: 95–112.
48. Rivers JW, Johnson M, Haig SM, Schwarz CJ, Burnett LJ, et al. (2014) An
analysis of monthly home range size in the critically endangered CaliforniaCondor. Bird Conserv Int 24: in press.
49. Rodgers AR, Carr AP, Smith L, Kie JG (2005) HRT: home range tools for
ArcGIS. Available: http://www.blueskytelemetry.com/downloads.asp. Accessed1 November 2011.
50. Ruxton GD, Houston DC (2004) Obligate vertebrate scavengers must be largesoaring fliers. J Theor Biol 228: 431–436.
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 12 February 2014 | Volume 9 | Issue 2 | e88430
51. Shamoun-Baranes J, Leshem Y, Yom-Tov Y, Liechti O (2003) Differential use
of thermal convection by soaring birds over central Israel. Condor 105: 208–218.
52. Shannon HD, Young GS, Yates MA, Fuller MR, Seegar WS (2002) American
white pelican soaring flight times and altitudes relative to changes in thermaldepth and intensity. Condor 104: 679–683.
53. Shepard ELC, Lambertucci SA, Vallmitjana D, Wilson RP (2011) Energybeyond food: foraging theory informs time spent in thermals by a large soaring
bird. PLoS ONE 6:e27375.
54. Silverman BW (1986) Density estimation for statistics and data analysis. London:Chapman and Hall. 176 p.
55. Smallwood KS, Thelander C (2008) Bird mortality in the Altamont Pass WindResource Area, California. J Wildl Manage 72: 215–223.
56. Snyder NF, Schmitt NJ (2002) California Condor (Gymnogyps californianus). TheBirds of North America Online. Available: http://bna.birds.cornell.edu/bna/
species/610. Accessed 1 November 2011.
57. Spaar R, Bruderer B (1996) Soaring migration of steppe eagles Aquila nipalensis insouthern Israel: flight behavior under various wind and thermal conditions.
J Avian Biol 27: 289–301.
58. Thomas PJ, Labrosse AK, Pomeroy AC, Otter KA (2011) Effects of weather on
avian migration in proposed ridgeline wind energy sites. J Wildl Manage 75:
805–815.
59. Tucker VA (1988) Gliding birds: descending flight of the White-backed Vulture,
Gyps africanus. J Exp Biol 140: 325–344.
60. U.S. Fish and Wildlife Service (1996) California Condor Recovery Plan.
Portland: U.S. Fish and Wildlife Service. 76 p.
61. U.S. Fish and Wildlife Service (2013) http://www.fws.gov/hoppermountain/
cacorecoveryprogram/CurrentStatus.html. Accessed 1 December 2013.
62. Walters JR, Derrickson SR, Fry DM, Haig SH, Marzluff JM, et al. (2010) Status
of the California condor (Gymnogyps californianus) and efforts to achieve its
recovery. Auk 127: 969–1001.
63. White GC, Garrott RA (1990) Analysis of wildlife radiotracking data. San Diego:
Academic Press. 383 p.
64. Worton BJ (1989) Kernel methods for estimating the utilization distribution in
home-range studies. Ecology 70: 164–168.
65. Young GS (1988) Convection in the atmospheric boundary layer. Earth Sci Rev
25: 179–198.
Resource Selection in the California Condor
PLOS ONE | www.plosone.org 13 February 2014 | Volume 9 | Issue 2 | e88430