Texas AgriLife Research golden-cheeked warbler surveys 979.845-1851 [email protected]irnr.tamu.edu The golden-cheeked warbler is a migratory songbird that breeds exclusively in central Texas. It was listed as endangered by the U.S. Fish and Wildlife Service in 1990 due to habitat loss and fragmentation. The Balcones Canyonlands National Wildlife Refuge and the Balcones Canyonlands Preserve in and near Austin were created, in part, to protect and manage habitat for the golden-cheeked warbler. These properties, along with the Fort Hood Military Reservation near Killeen, have been the site of many studies for more than 15 years to increase understanding of the golden-cheeked warbler’s biology and ecology. However, these properties are only a small proportion of the warbler’s range, as habitat can be found in parts of nearly 40 counties in central Texas. Questions still exist about its habitat use and abundance throughout much of that area. There have been only a few attempts over the years to estimate the population size of the golden-cheeked warbler and the extent of its habitat throughout its breeding range. Those estimates (some now over 20 years old) were often based on a limited amount of data and, depending on the researchers’ data and analyses, varied from roughly 9,000 to 54,000 individuals; other researchers have suggested that sufficient habitat exists to support over 228,000 individuals. Given this limited, outdated information, the Texas A&M Institute of Renewable Natural Resources, a unit of Texas AgriLife Research, initiated a study in 2008 to determine the distribution and abundance of golden-cheeked warblers in central Texas. Compared to other studies on the golden-cheeked warbler, the AgriLife Research study was the first to survey for the presence of warblers and their habitat over such a large area to estimate the warbler’s population size. Statistical models were developed using data from these surveys along with satellite imagery depicting potential warbler habitat. The imagery and models suggest approximately 262,000 male warblers occur within 4 million acres of potential habitat. The research team released a report in 2010 detailing the sampling methods, statistical analyses, and results of its study. This report was soon followed by three related publications in scientific journals that expanded upon and clarified sections of the original report. Drafts of these publications were reviewed by well-respected scientists before being submitted to the journals’ own peer-review processes. The population size suggested in the AgriLife Research study is many times higher than previous estimates, and it quickly drew the attention of proponents and critics. The U.S. Fish and Wildlife Service requested a peer review of the report and associated manuscripts because the importance of the study put it in the category of “Highly Influential Scientific Information”; the Service wanted to determine whether the best available science had been used and interpreted in a reasonable way. The peer review was facilitated by The Wildlife Society, which is the scientific society for professional wildlife biologists and managers. The reviewers concluded that the overall study design, analyses, and inferences are supported by sound scientific data and analysis. They identified some of the limitations of the approach but concluded that the models used in the study were “defensible and allow for reliable conclusions” and “appropriate given the scale and scope of the project.” The AgriLife Research scientists agree with the peer reviewers that additional research is needed to validate and refine the approach and to further improve the reliability of the results. The results of this study have generated a variety of new questions about golden-cheeked warblers. The models developed in the study are a representation of reality but should be refined over time as new data is collected and analyzed. However, the models cannot tell us how many golden-cheeked warblers existed in the past, or how the population may have increased or decreased over time. The models also do not apply to the warbler’s winter range in southern Mexico and Central America, where knowledge of the species distribution and habitat use is limited. The report and related publications do not advocate changes in the listing status of the golden-cheeked warbler nor do they imply that conservation measures to protect the species habitat are no longer needed. Rather, this study is one of many necessary steps in the evolving knowledge of the golden- cheeked warbler and is not intended as the final word on the matter. Science is an iterative process—new ideas develop as new data and analyses become available. Research on the golden-cheeked warbler is no exception, and focusing on a different approach to understanding the species (e.g., range- wide surveys in addition to site-specific studies) can help advance that process. Reports and publications by the Texas A&M Institute of Renewable Natural Resources can be found at http://irnr.tamu. edu/publications/.
31
Embed
Texas AgriLife Research golden-cheeked warbler surveystwri.tamu.edu/docs/ramses/gcw_survey_informational_packet.pdf · Texas AgriLife Research golden-cheeked warbler surveys 979.845-1851
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
Texas AgriLife Research golden-cheeked warbler surveys
The golden-cheeked warbler is a migratory songbird that breeds exclusively in central Texas. It was listed as endangered by the U.S. Fish and Wildlife Service in 1990 due to habitat loss and fragmentation. The Balcones Canyonlands National Wildlife Refuge and the Balcones Canyonlands Preserve in and near Austin were created, in part, to protect and manage habitat for the golden-cheeked warbler. These properties, along with the Fort Hood Military Reservation near Killeen, have been the site of many studies for more than 15 years to increase understanding of the golden-cheeked warbler’s biology and ecology. However, these properties are only a small proportion of the warbler’s range, as habitat can be found in parts of nearly 40 counties in central Texas. Questions still exist about its habitat use and abundance throughout much of that area.
There have been only a few attempts over the years to estimate the population size of the golden-cheeked warbler and the extent of its habitat throughout its breeding range. Those estimates (some now over 20 years old) were often based on a limited amount of data and, depending on the researchers’ data and analyses, varied from roughly 9,000 to 54,000 individuals; other researchers have suggested that sufficient habitat exists to support over 228,000 individuals. Given this limited, outdated information, the Texas A&M Institute of Renewable Natural Resources, a unit of Texas AgriLife Research, initiated a study in 2008 to determine the distribution and abundance of golden-cheeked warblers in central Texas.
Compared to other studies on the golden-cheeked warbler, the AgriLife Research study was the first to survey for the presence of warblers and their habitat over such a large area to estimate the warbler’s population size. Statistical models were developed using data from these surveys along with satellite imagery depicting potential warbler habitat. The imagery and models suggest approximately 262,000 male warblers occur within 4 million acres of potential habitat.
The research team released a report in 2010 detailing the sampling methods, statistical analyses, and results of its study. This report was soon followed by three related publications in scientific journals that expanded upon and clarified sections of the original report. Drafts of these publications were reviewed by well-respected scientists before being submitted to the journals’ own peer-review processes. The population size suggested in the AgriLife Research study is many times higher than previous estimates, and it quickly drew the attention of proponents and critics. The U.S. Fish and Wildlife Service requested a peer review of the report
and associated manuscripts because the importance of the study put it in the category of “Highly Influential Scientific Information”; the Service wanted to determine whether the best available science had been used and interpreted in a reasonable way. The peer review was facilitated by The Wildlife Society, which is the scientific society for professional wildlife biologists and managers. The reviewers concluded that the overall study design, analyses, and inferences are supported by sound scientific data and analysis. They identified some of the limitations of the approach but concluded that the models used in the study were “defensible and allow for reliable conclusions” and “appropriate given the scale and scope of the project.” The AgriLife Research scientists agree with the peer reviewers that additional research is needed to validate and refine the approach and to further improve the reliability of the results.
The results of this study have generated a variety of new questions about golden-cheeked warblers. The models developed in the study are a representation of reality but should be refined over time as new data is collected and analyzed. However, the models cannot tell us how many golden-cheeked warblers existed in the past, or how the population may have increased or decreased over time. The models also do not apply to the warbler’s winter range in southern Mexico and Central America, where knowledge of the species distribution and habitat use is limited.
The report and related publications do not advocate changes in the listing status of the golden-cheeked warbler nor do they imply that conservation measures to protect the species habitat are no longer needed. Rather, this study is one of many necessary steps in the evolving knowledge of the golden-cheeked warbler and is not intended as the final word on the matter. Science is an iterative process—new ideas develop as new data and analyses become available. Research on the golden-cheeked warbler is no exception, and focusing on a different approach to understanding the species (e.g., range-wide surveys in addition to site-specific studies) can help advance that process.
Reports and publications by the Texas A&M Institute of Renewable Natural Resources can be found at http://irnr.tamu.edu/publications/.
Estimating Breeding Season Abundance ofGolden-Cheeked Warblers in Texas, USA
HEATHER A. MATHEWSON,1 Texas A&M Institute of Renewable Natural Resources, 1500 Research Parkway, Suite 110, Texas A&MUniversity, College Station, TX 77843, USA
JULIE E. GROCE, Texas A&M Institute of Renewable Natural Resources, 2632 Broadway, Suite 301 S., Texas A&M University, San Antonio,TX 78215, USA
TIFFANY M. MCFARLAND, Texas A&M Institute of Renewable Natural Resources, 1500 Research Parkway, Suite 110, Texas A&M University,College Station, TX 77843, USA
MICHAEL L. MORRISON, Department of Wildlife and Fisheries Sciences, 210 Nagle Hall, Texas A&M University, CollegeStation, TX 77843, USA
J. CAL NEWNAM, Texas Department of Transportation, P.O. Box 15426, Austin, TX 78761, USA
R. TODD SNELGROVE, Texas A&M Institute of Renewable Natural Resources, 1500 Research Parkway, Suite 110, Texas A&M University,College Station, TX 77843, USA
BRET A. COLLIER, Texas A&M Institute of Renewable Natural Resources, 1500 Research Parkway, Suite 110, Texas A&M University, CollegeStation, TX 77843, USA
R. NEAL WILKINS, Texas A&M Institute of Renewable Natural Resources, 1500 Research Parkway, Suite 110, Texas A&M University, CollegeStation, TX 77843, USA
ABSTRACT Population abundance estimates using predictive models are important for describing habitatuse and responses to population-level impacts, evaluating conservation status of a species, and for establishingmonitoring programs. The golden-cheeked warbler (Setophaga chrysoparia) is a neotropical migratory birdthat was listed as federally endangered in 1990 because of threats related to loss and fragmentation of itswoodland habitat. Since listing, abundance estimates for the species have mainly relied on localizedpopulation studies on public lands and qualitative-based methods. Our goal was to estimate breedingpopulation size of male warblers using a predictive model based on metrics for patches of woodland habitatthroughout the species’ breeding range. We first conducted occupancy surveys to determine range-widedistribution.We then conducted standard point-count surveys on a subset of the initial sampling locations toestimate density of males. Mean observed patch-specific density was 0.23 males/ha (95% CI ¼ 0.197–0.252,n ¼ 301). We modeled the relationship between patch-specific density of males and woodland patchcharacteristics (size and landscape composition) and predicted patch occupancy. The probability of patchoccupancy, derived from a model that used patch size and landscape composition as predictor variables whileaddressing effects of spatial relatedness, best predicted patch-specific density. We predicted patch-specificdensities as a function of occupancy probability and estimated abundance of male warblers across 63,616woodland patches accounting for 1.678 million ha of potential warbler habitat. Using a Monte Carlosimulation, our approach yielded a range-wide male warbler population estimate of 263,339 (95% CI:223,927–302,620). Our results provide the first abundance estimate using habitat and count data from asampling design focused on range-wide inference. Managers can use the resulting model as a tool to supportconservation planning and guide recovery efforts. � 2012 The Wildlife Society.
KEY WORDS abundance, density, endangered species, golden-cheeked warbler, point count, population estimate,Setophaga chrysoparia, Texas.
Abundance estimates are of particular importance for evalu-ating conservation status and determining recovery goals,establishing monitoring programs, describing habitat usepatterns, and assessing population-level impacts driven byanthropogenic and natural factors (Campbell et al. 2002,Scott et al. 2005, Fitzgerald et al. 2009, Sirami et al. 2010).Population size estimation is a challenge for most species, butapproaches integrating remotely sensed data with predictive
models can assist in predicting abundance at large spatialscales (Thompson 2002a, Fitzgerald et al. 2009). The gold-en-cheeked warbler (Setophaga chrysoparia) is a neotropicalmigratory songbird that breeds only in central Texas andwinters in the highlands of southern Mexico and CentralAmerica (Pulich 1976, Groce et al. 2010). In 1990, theUnited States Fish and Wildlife Service (USFWS) listedthe golden-cheeked warbler (hereafter warbler) as endan-gered and cited habitat loss and fragmentation as primarythreats (USFWS 1990, 1992). Warbler occurrence, density,and recruitment rates appear to decrease as the size ofhabitat patches and the amount of habitat in the surrounding
Received: 25 March 2011; Accepted: 28 November 2011
The Journal of Wildlife Management; DOI: 10.1002/jwmg.352
Mathewson et al. �Warbler Population Abundance 1
landscape decline (DeBoer and Diamond 2006, Magnesset al. 2006, Peak 2007, Butcher et al. 2010, Collier et al.2010).Previously, approximations of golden-cheeked warbler
abundance within the breeding range were based on esti-mates of warbler density collected at a few (<20) study sitesand extrapolated based on projected extents of breedinghabitat (Pulich 1976, Wahl et al. 1990, Rowell et al.1995, Rappole et al. 2003; Table 1). Population estimateshave ranged between 9,644 and 32,032 individuals (Pulich1976, Wahl et al. 1990), whereas estimates of carryingcapacity have ranged from 64,520 to 228,426 individuals(Rowell et al. 1995, Rappole et al. 2003). Variation inpopulation estimates results from the methods used to esti-mate the extent of habitat, and assumptions regarding whatcharacteristics are representative of warbler habitat (Table 1).Our objective was to develop a range-wide estimate of
abundance for male warblers. We relied on a range-widemodel predicting patch level occupancy (Collier et al. 2012)to serve as our sampling frame. We then used point countsurveys to estimate patch-specific density of male warblers.We evaluated predictive relationships between patch-specificdensity and remotely sensedmetrics of habitat patches. Usingthese relationships, we combined our density estimates withthe range-wide occupancy model to yield estimates of malewarbler abundance across the species’ breeding range.
STUDY AREA
The golden-cheeked warbler breeding range is confined tocentral Texas, USA, on the eastern half of the EdwardsPlateau and the southern half of the Cross Timbers ecor-egions (Hatch et al. 1990).We conducted our research acrossthe breeding range on public and private properties in 35counties (USFWS 1992; Fig. 1). Our sampling units werepatches of potential warbler habitat, characterized as oak-juniper woodlands (Collier et al. 2010, 2012), and we defineda patch as a relatively homogenous unit of vegetation distinctfrom its surroundings (Kotliar and Wiens 1990). To guideus in determining survey locations in 2008, we used ahabitat classification developed by SWCA EnvironmentalConsultants (SWCA 2007) that delineated potential warblerhabitat as those patches having >50% canopy closure by amixture of mature or second-growth Ashe juniper (Juniperusashei) and deciduous hardwood trees based on 2004 NationalAgricultural Imagery Program color infrared digital imagery(1-m resolution). SWCA (2007) excluded patches <4 haunless immediately adjacent (unspecified distance) to otherpatches of potential habitat. SWCA’s restricted definition ofhabitat yielded 7,865 patches (approx. 552,000 ha) of po-tential habitat across the breeding range. Based on our fieldsurveys during 2008, as well as previous work within thissystem (Butcher et al. 2010, Collier et al. 2010), we found theSWCA (2007) classification narrowly defined warbler habi-tat; thus, we developed a broader classification of warblerhabitat patches for selecting survey sites in 2009 and forevaluating relationships between patch-specific densities andpatch metrics.
We delineated potential habitat using an unsupervisedclassification of woodlands from 2007 and 2008 cloud-freeLandsat 5 imagery collected during late spring (Morrisonet al. 2010). Using the ISODATA clustering algorithm inLeica ERDAS Imagine 9.3 (Intergraph Corporation,Norcross, GA), we grouped imagery spectral response pat-terns into 20 statistically different clusters. Using high reso-lution aerial photography from the National AgricultureImagery Program (NAIP; 2004 and 2006) and data collectedin the field, we identified those clusters that corresponded towoodland land cover, grouping all other land cover types asnon-habitat (e.g., wetlands, cropland, urban areas, water,barren, impervious surfaces, grassland). Several woodlandclusters were based on characteristics such as slope, aspect,shadows, and composition (e.g., evergreen or deciduous);however, because classification accuracy of these woodlandswas low and our intent was to sample woodland cover, weaggregated these clusters into 1 class. To further refine ourwoodland classification, we used the 2001 National LandCover Data set (NLCD; Homer et al. 2007) to eliminate anymisclassified non-woodland land cover types from our initialclassification. We also removed those pixels classified aswoodland but with canopy cover <30% using the 2001NLCD canopy cover layer. Using road data from theTexas Strategic Mapping Program (STRATMAP), we de-fined breaks between patches by removing pixels that inter-sected paved or public roads. We thus delineated 63,616patches (mean patch size ¼ 26.39 ha, range ¼ 2.8–26,967 ha) or approximately 1.678 million ha of potentialwarbler habitat (Fig. 1). We used ArcGIS 9.3.1(Environmental Systems Research Institute, Inc.,Redlands, CA) and calculated patch size (ha), landscapecomposition (% woodland within a 400-m radius of a givenpixel; Magness et al. 2006), patch core area, and edge-to-arearatio for each of the 63,616 patches of mixed woodlandidentified within the study area. We calculated the landscapecomposition value for patches as the mean value for all pixelswithin the patch.We calculated patch core habitat area as theinternal portion of the patch after buffering the exterior edgeof the patch internally by 30 m. The distribution of patchesused to determine 2008 survey points was biased towardlarge, contiguous patches, but combined with survey pointsfrom 2009 the distribution of surveyed patches included awider range of patch sizes, shapes, and amount of edgehabitat. Thus, we re-projected point count survey locationsfrom 2008 sampling to those patches identified in our habitatdelineation for analysis.
METHODS
Occupancy SurveysWe used a 2-phase sampling approach (Conroy et al. 2008)whereby we first conducted occupancy (i.e., detection or non-detection) surveys at the patch scale to determine areas thatlikely supported warblers (Collier et al. 2012). For occupancysurveys, we used a double-observer, removal approach(MacKenzie et al. 2006) in which 2 observers simultaneouslyyet independently surveyed a habitat patch up to a maximum
of 6 surveys (3 double-observer surveys) to determine warblerpresence (Collier et al. 2010). We surveyed 434 patches forwarbler presence in 2008 and 2009 and if we detected awarbler during occupancy surveys, we returned to the patchto conduct abundance surveys (see below).To select patches for occupancy surveys, we stratified
patches based on administrative units defined by USFWS(USFWS 1992) and used a probability proportional to size(PPS) sampling design (Thompson 2002a), in which werandomly selected habitat patches proportional to the distri-bution of patch sizes across the range. Because habitatpatches often spanned multiple public and private land-owners, we attempted to locate and contact all propertyowners of randomly selected patches. However, our accessto the entirety of each patch was sometimes limited becauseof inability to locate landowners, lack of response, unwill-ingness to participate by landowners, or other logisticalconstraints. Because landowners tend to base their partici-pation on multiple factors that are often not associated withland management practices that would influence warblerabundance (Hilty and Merenlender 2003, DeBoer andDiamond 2006, Sorice et al. 2011), we assumed that suchaccess restrictions did not bias our sampling. Furthermore,we assumed that properties for which we were unable toacquire access were missing from our sample completely atrandom (Stevens and Jensen 2007). We did not stratify byownership type so we focused considerable effort on acquir-
ing access to private property. Our approach further assumedthat variability in habitat conditions was not related consis-tently to designation as private versus public lands as fewpublic or private properties are managed specifically forwarbler habitat in Texas (Groce et al. 2010).
Abundance Surveys
Within each patch selected for abundance surveys, we buff-ered each survey patch to reduce the proportion of survey areathat would fall outside of a habitat patch, such that pointcount centers in patches >40 ha were located �100 m fromthe edge and points in patches�40 ha were�25 m from theedge. Under those constraints, we randomly established themaximum number of point count stations while ensuring a400-mminimum distance between point count centers usingArcGIS 9.3. Small patches (<10 ha) that were inadvertentlyskipped because they violated the 400-m spacing constraintwere processed a second time without the spacing constraintto ensure adequate coverage of point survey locations withinsmall patches. Given our spacing and random placement, weassumed abundance counts at each point were independentwithin patches. Because our sampling unit was the patch, wedid not incorporate any within-patch habitat variability thatcould influence local differences in male density. We as-sumed that our random distribution of point count stationscaptured any within-patch habitat heterogeneity in localvegetation or densities.
Figure 1. Distribution of mixed woodland patches, and locations of abundance surveys conducted during 2008 and 2009 on private or public lands in the 35county breeding range of the golden-cheeked warbler in Texas, USA.
4 The Journal of Wildlife Management
Patches selected for point count surveys were conditionalon positive detections of warblers during occupancy surveys,but we conducted point count surveys independent of theoccupancy sampling process. As such, a zero patch-levelabundance estimate was a recordable result when we detectedno warblers during point count surveys. During mid-Marchto mid-May 2008 and 2009, we conducted abundance sur-veys for warblers using 100-m fixed-radius point countsfollowing methods detailed by Laake et al. (2011).Surveys began at sunrise and ended no later than 13:00.When surveyors detected a bird during occupancy surveys,they immediately initiated point count surveys at the pre-determined point count locations. However, if sufficient timewas not available to conduct the surveys within the same day,we returned within a week of the occupancy survey. We useda dependent double-observer sampling approach (Cook andJacobson 1979, Nichols et al. 2000). At each survey point, werandomly assigned primary or secondary observer status.During the 5-minute survey, the primary observer commu-nicated visual and auditory detections of male warblers to thesecondary observer, noting the direction and classifying theapproximate distance into a distance bin (0–50 m or 50–100 m). Concurrently, the secondary observer recordeddetections of any individuals missed by the primary observer(Laake et al. 2011). The combined data from the primary andsecondary observers represents a 2-sample capture history.To maximize the spatial distribution of sample locations, wevisited each patch on 1 occasion for abundance surveys(Thompson 2002b).We did not conduct point counts duringinclement weather or periods of high wind. We followed thestandard assumptions for point count surveys: 1) populationclosure, 2) observers correctly identified birds, 3) no doublecounting, and 4) observers correctly estimated distances tobirds (Buckland 2006, Johnson 2008).
Analysis
Tomodel male warbler abundance across the breeding range,we first estimated observed patch-specific density using ourpoint count survey data. Next, using an information-theo-retic approach to select the best fitting model given the data(Burnham and Anderson 2002); we regressed patch-specificpredictor variables against the observed patch-specific den-sity estimates and developed a predictive equation relatingbiological metrics to density at the patch scale. Using thepredictive equation, we predicted densities (as well as lowerand upper 95% CL) for each patch of potential habitat fromour delineation across the range and converted those topatch-specific abundance. Finally, we used a Monte Carlosimulation to randomly identify occupied patches based ontheir predicted occupancy probability (Collier et al. 2012)and summed the resultant abundance estimates for eachrandomly selected patch to generate a population abundanceestimate with lower and upper 95% confidence limits for each1,000 realizations.We combined point count data collected during 2008
and 2009 because we independently selected point countlocations for each study year and we assumed minimal annualvariation of patch-level warbler densities. We examined our
assumption using the patches randomly selected for abun-dance surveys in both study years. We estimated observedpatch-specific density using the count of male warblers stan-dardized to the total surveyed area (i.e., 3.14 ha for 100-mfixed-radius sample area) within each patch. We assumedthat the effect of 100-m radius point count locations thatincluded non-woodland vegetation based on our patch de-lineation was minimal. We used counts of abundance un-corrected for detection probability because a preliminaryevaluation revealed estimated double observer detection ratesfor warblers using our survey design were high (probability ofdetection ¼ 0.97; Laake et al. 2011). Furthermore, our un-corrected density estimates will be comparable to previousresearch on warbler densities and abundance that have notincluded detection corrections in their estimates (Table 1).We acknowledge that by implicitly assuming a detectionprobability of 1.0 we likely underestimated density; thus,our estimates may be conservative (Pollock et al. 2002,Thompson 2002b). We assumed that density was constantacross patches and that the random distribution of oursampling points captured any within-patch vegetationheterogeneity.To evaluate the relationships between observed density of
male warblers and our predictor variables, we used general-ized linear modeling (GLM) with a negative binomial dis-tribution of the raw count data (McCullagh and Nelder1989, White and Bennets 1996). We included total areasurveyed at the patch scale as an offset term to standardizeour count data to a density estimate for comparison acrosspatches. We used a negative binomial distribution of ourcount data instead of the Poisson distribution commonlyused for count data because the negative binomial provides aflexible approach to addressing overdispersion (White andBennets 1996).We assembled 8 competing models (Table 2) for predicting
male density that were based on patch-specific metrics andprobability of patch occupancy derived from the model de-veloped during a concurrent study and described in Collieret al. (2012). We examined correlations among our predictorvariables before determining our final set of candidate mod-els. Core area correlated with patch size (r ¼ 0.997) andedge-to-area ratio correlated with landscape composition(r ¼ �0.816), thus we used patch size and landscape com-position as predictors of patch density. We modeled thelinear relationship between male patch-specific densityand patch size because golden-cheeked warblers are consid-ered to be an area-sensitive species (Butcher et al. 2010) inthat demographic parameters are positively associated withpatch size (Donovan and Flather 2002). Additionally, weexamined a quadratic trend in patch size assuming that onceterritory density reached a threshold, density remained con-stant regardless of increases in patch size. We used landscapecomposition as an alternative predictor because landscape-scale factors might drive mechanisms associated with settle-ment decisions of territorial males, such as conspecificattraction (Campomizzi et al. 2008, Farrell 2011). Wealso considered an interaction between patch-size and land-scape composition. We did not consider distance between
Mathewson et al. �Warbler Population Abundance 5
patches because our landscape composition metric capturedthis variability and previous research on warbler occurrenceat the point-scale indicated that these variables were notas predictive as is landscape composition (DeBoer andDiamond 2006, Magness et al. 2006).We also used the patch occupancy prediction from Collier
et al. (2012) as a predictor variable of abundance. The geo-additive semiparametric approach detailed in Collier et al.(2012) predicted occupancy probability (ci) as a function ofcovariate data (patch size and landscape composition) andspatial location within the warbler’s range expressed as
logitðciÞ ¼ b0 þ blXl þX20
k¼1ukðLocationi � kkÞ þ "i
where the bl Xl represents a vector of l predictor variables(patch size, landscape composition, patch size-landscapecomposition interaction, and patch-specific UniversalTransverseMercator [UTM] coordinates) entering the mod-el linearly and ðLocation� kkÞ represents the effect of spatiallocation for each surveyed habitat patch (Gimenez et al.2009). Collier et al. (2012) used a temporal covariate rep-resenting sample survey date for detection modeling becausedate of survey has been shown to adequately predict detectionrates of warblers at the patch scale (Collier et al. 2010).We used an information-theoretic approach (Burnham and
Anderson 2002) to evaluate competing models. We rankedmodels using Akaike’s Information Criterion (AICc) cor-rected for small sample sizes. We calculated Akaike weights(wi) to indicate relative likelihood of each model as the bestapproximating model given the data for our set of candidatemodels. We considered models �4 DAIC units to be com-petitive. We judged model adequacy by performing a 10-foldcross-validation. We performed a likelihood ratio test of asaturated model that included occupancy and an interactionbetween patch size and landscape composition compared to aconstant model to examine over-dispersion (Venables andRipley 2002).We used the best fitting model given our model set and
predicted patch-specific density of male warblers for eachhabitat patch that we delineated across the range
(n ¼ 63,616 patches). We converted predicted density topredicted patch-specific abundance and ran Monte Carlosimulations using the equation
N ¼Xn
i¼1ðyi � area� kiÞ
where yi is the predicted patch-specific density estimate andarea (ha) is the size of each patch. The covariate k representeda random assignment of occupancy (0 or 1) based on randomdraws from a Bernoulli distribution using the predicted patchoccupancy rates from Collier et al. (2012). We generated adistribution of abundance estimates and associated upper andlower bounds using 1,000 replicates of the above simulationwhere the summed abundance values for those patches wherek ¼ 1 represented the predicted population.When surveying avian species, counting a single bird mul-
tiple times is possible, which would bias point survey data(Buckland 2006). Given our short point count period(5 min), we concluded that double counting was unlikely;however, we were interested in exploring the potentialimpacts of being wrong in reaching such a conclusion.Thus, we estimated minimum density and population sizebased on 2 possible scenarios of double counting. First, giventhat double counting could not occur if observers detected 0or 1 warbler at a point during a survey, and that the likelihoodof observers having over counted when reporting 2 birds isunlikely given our survey design, we re-examined our data bylimiting total counts at a point to a maximum of 2 (Table 3).Second, to exclude potential bias of double counting warblersbetween the 2 observers, we used the number of birdscounted only by the primary observer, representing the min-imum count that would be obtained from a single-visitstandard point count (Nichols et al. 2000).
RESULTS
Abundance Counts and Density EstimationWe conducted 1,057 point count surveys in 301 mixedwoodland patches (2008: n ¼ 151 patches; 2009: n ¼ 150)across the warbler’s breeding range (Fig. 1). Survey points per
Table 2. Model selection statistics formodels using negative binomial regression explaining patch-specific density of golden-cheekedwarblers in Texas, in 2008and 2009.
a Occ ¼ estimated occupancy from semiparametric model of Collier et al. (2012), Size ¼ size of habitat patch (ha), Size2 ¼ size of habitat patch quadratictrend, LS ¼ landscape composition.
b K ¼ no. of parameters.c AICc ¼ Akaike’s Information Criterion corrected for small sample sizes.d DAICc ¼ AICc relative to the most parsimonious model.e wi ¼ AICc model weight.
6 The Journal of Wildlife Management
patch ranged from 1 to 35 (mean ¼ 3.5, SD ¼ 4.11). Wedetected 796 warblers during both years of our surveys. Forpatches surveyed for abundance estimation, patch size rangedfrom 2.8 to 26,967 ha (mean ¼ 743.7, 95% CI ¼ 497.2–990.2, median ¼ 141.2, n ¼ 301). Mean patch size forpatches primarily in private ownership ranged from 3.2 to26,970 ha (mean ¼ 740.1, 95% CI ¼ 388.3–1,093.5,median ¼ 147.3, n ¼ 178) and for public properties rangedfrom 2.8 to 11,880 ha (mean ¼ 747.8, 95% CI ¼ 420.3–1,075, median ¼ 119.2, n ¼ 123). Mean landscape compo-sition (% woodland vegetation within 400-m radius) rangedfrom 15% to 92% (mean ¼ 65.4, 95% CI ¼ 63.7–67.0,median ¼ 67.1, n ¼ 301). Mean landscape compositionfor patches primarily in private ownership ranged from15% to 91% (mean ¼ 64, 95%CI ¼ 62–66, median ¼ 65.7,n ¼ 178) and for public properties ranged from 26% to 92%(mean ¼ 67.2, 95% CI ¼ 64.3–70.1, median ¼ 72.8,n ¼ 123). For the sub-sample of patches surveyed in bothyears (n ¼ 34), density estimates were similar (2008:0.25 males/ha, 95% CI ¼ 0.18–0.32; 2009: 0.23 males/ha,95% CI ¼ 0.15–0.29), supporting our assumption of mini-mal annual variation between our survey years.Thirty-four percent (n ¼ 301) of patches in which we
conducted abundance point counts had no warbler detec-tions. Mean observed patch-specific density of male warblersfor both years, was 0.23 males/ha (95% CI ¼ 0.197–0.252,n ¼ 301 patches). Density varied between USFWS recoveryregions (USFWS 1992); mean observed density in the northregion was 0.15 males/ha (95% CI ¼ 0.10–0.17, n ¼ 86patches), the central region equaled 0.19 males/ha (95%CI ¼ 0.15–23, n ¼ 106), and the south region equaled0.32 males/ha (95% CI ¼ 0.26–37, n ¼ 109).
Model Results and Abundance Estimation
The goodness-of-fit test indicated that the saturated modelfit the data (x24 ¼ 52.1, P < 0.001). Based on our modelselection results, patch-specific occupancy probability bestpredicted male warbler density (Table 2). The standardizedmodel parameter estimates for the fitted model indicated thatfor every 0.25 increase in predicted patch occupancy proba-bility (i.e., 1-unit increase for scaled data), the density ofmales increased by 61% (b ¼ 0.478, 95% CI: 0.326–0.636).Using our predicted patch-specific density estimates as a
function of predicted patch occupancy probability and basedon 1,000 simulated realizations of population distribution,we estimated median singing male warbler abundance acrossthe range as 262,013 singing males (95% CI: 223,164–
301,081; Fig. 2). Under the assumption that the maximumnumber of male warblers detected was 2, predicted densityand estimated population size declined 10.9% (Table 4).Removal of the secondary observer’s data resulted in a13.4% decrease in density and population size estimates(Table 4). Cross-validation indicated a predictive error of15% but upon removal of 1 extreme outlier the error droppedto 5%. We summarized estimates of golden-cheeked warblerhabitat area and abundance from the literature (Table 1).Methods used to estimate the amount of habitat varied butgenerally used Landsat imagery. In all cases, however,authors gathered data at only a few locations and thenextrapolated across the estimated habitat range.
DISCUSSION
We estimated male golden-cheeked warbler population sizeacross the breeding range in central Texas, USA using arelationship that accounts for predicted occupancy of a habi-tat patch, which in turn reflects the geographic spatial loca-tion, patch size, and landscape composition (Collier et al.2012). The relationship with occupancy suggests that abun-dance decreased from south to north, and with decreasingpatch size and woodland cover. We did not find a relation-ship between abundance and these environmental covariatesper se, but instead the relationship relied on spatial locationwithin the warbler range. Indeed, estimates from our rawcount data (observed patch-specific density), when groupedby USFWS Recovery Regions, indicated that density wasgreater for regions in the south than in the north. Thispattern was not unexpected because the southern portionof the range contains the greatest percentage of availablewarbler habitat (Groce et al. 2010), and the habitat typicallyoccurs in less fragmented, larger patches that are consideredmore suitable habitat for warblers (DeBoer and Diamond2006, Baccus et al. 2007; Fig. 3). Our study suggested thatmore warblers exist than previously estimated (Table 1), orthat the carrying capacity of available habitat is greater, andbelow we discuss how differences in study design, analysis,and assumptions contributed to discrepancies in populationsize estimates.A positive relationship between occupancy probability and
local abundance is well formulated in biology (Gaston et al.2000, He and Gaston 2003). Indeed, some degree of thisrelationship is inevitable at the extreme limits of populationsize: local abundance necessarily increases when habitat issaturated, and at low levels of abundance, occurrence patternsare not maintained (Gaston 1999). Density-dependent hab-itat selection is one mechanism that might explain the abun-dance–occupancy relationship in golden-cheeked warblerssince the species demonstrates conspecific attraction(Campomizzi et al. 2008, Farrell 2011) but other mecha-nisms also may contribute (see review Gaston et al. 2000).Given that the occupancy estimates are a function of patchsize and landscape composition (Collier et al. 2012), theabundance–occupancy relationship is likely a response tovariations in habitat conditions if these metrics predict habi-tat quality (i.e., survival or reproductive success). Similarly,lower predicted occupancy and abundance in the northern
Table 3. Number of survey points and total number of singingmale golden-cheeked warblers counted during point-count surveys in central Texas in2008 and 2009.
portion of the warbler’s range might be a consequence ofinterrupted dispersal dynamics because of landscapes withsmaller or isolated habitat patches (Donovan and Flather2002).The relationship between warbler density and patch-scale
metrics that we used to predict abundance across the species’range was consistent with patch-scale metrics previouslyshown to affect warbler density at local scales (Magnesset al. 2006, Baccus et al. 2007). Species–habitat relationshipsare complex and studies of this broad-scale application nec-essarily rely on assumptions concerning the validity of se-lected habitat metrics to define the relationship (Boyce andMcDonald 1999, Fitzgerald et al. 2009). Warbler densitiesare likely to vary within patches in response to habitatcharacteristics such as tree species composition, density ofjunipers or oaks, and relative age of woodland (see review inGroce et al. 2010), variables that are important to definehabitat at finer scales (Fitzgerald et al. 2009). We did notattempt to incorporate these relatively fine-scale variables,because our intent was to analyze patch-scale relationshipsacross a broad geographic range.We depended upon random
sampling to account for the within-patch variation in warblerdensities, including any variation that might result fromwithin-patch differences in habitat. With future use ofmore advanced remote sensing technologies, our modelcould incorporate additional landscape-scale and local-scalemetrics. As with any large-scale study, precision and accuracyat finer scales are inherently reduced and further refinementof our model is needed for application at local scales.
Patch-Specific Density
Our study differs from previous approaches that estimatedpopulation size of this species in that we estimated densityfrom random samples across the species’ breeding range andwe based our method of inference to potential habitat on apredictive species–habitat relationship at the patch-scale(Table 1). Our density estimates of males were from ran-domly sampled point count locations across the 35-countybreeding range, spanning a range of public and privateproperties. Even though inferential expansion to a broadgeographic extent assumes that sampling is representativeof the conditions across that area (Gutzwiller and Barrow
Table 4. Estimates and 95% confidence intervals of golden-cheeked warbler density and predicted warbler population size in central Texas in 2008 and 2009based on 3 data analysis options: observed estimates with no correction for detection probability, limiting individual counts to a maximum of 2 detections, andestimates for data from observer 1 only.
Method Patch-specific density (95% CI) Predicted population size (95% CI)
Uncorrected 0.224 (0.197–0.252) 262,013 (223,164–301,081)Maximum 2 count 0.205 (0.182–0.228) 234,862 (202,031–267,601)Observer 1, not corrected for detection 0.194 (0.169–0.219) 226,871 (191,658–262,091)
Figure 2. Predicted distribution for the range-wide population abundance of male golden-cheeked warbler estimates and associated 95% lower (LCL) andupper (UCL) confidence limits in Texas, USA, 2008–2009, based on predicted patch-specific density estimates as a function of predicted patch occupancyprobability.
8 The Journal of Wildlife Management
2001) no previous study reporting population sizes of gold-en-cheeked warblers has used density estimates from ran-domly located study sites. Previous calculations of warblerabundance used density estimates primarily from study sitesin the north or central portions of the breeding range (e.g.,Pulich 1976, Rappole et al. 2003). Furthermore, some pop-ulation size estimates were based on a mean or mediandensity estimate from a few locations (Wahl et al. 1990,Rowell et al. 1995, Rappole et al. 2003), implying a constantdensity and thus failing to address variation in densitiesamong habitat patches and across the range (Royle andNichols 2003). Although Wahl et al. (1990) estimatedsite-specific densities at 17 sites in 11 counties, they didnot incorporate variation in abundance in their populationsize estimate because they used the median density estimate,assuming constant density across the range. Pulich (1976)addressed regional variation in density by deriving estimatesfrom 3 study sites across the range; however, for projection topotential habitat he assumed a constant density within 3habitat assessment ranks (Table 1). We addressed regionaland among-patch heterogeneity in density by defining habi-tat patches as our sampling unit and projecting patch-specific
densities to our potential habitat delineation. Although ourapproach requires the assumption of constant density withinpatches, we suggest that this is an improvement from previ-ous estimates of golden-cheeked warbler abundance.We used standard methods for avian point count surveys
and we assumed a detection probability of 1.0. It is widelyrecognized that assuming a constant detection probabilityresults in biased estimates because failing to account fordetection typically produces underestimates of abundance(Pollock et al. 2002, Thompson 2002b, Johnson 2008).We accepted the potential for a negative bias because evi-dence that the use of a dependent double-observer closecapture model resulted in a high probability (>0.97) ofdetecting a male warbler by at least 1 of the 2 observers(Laake et al. 2011). Multiple factors influence detectability,including skills of observers, habitat, ambient noise, weather,distance, and temporal effects (Johnson 2008). Correctingfor detectability requires additional assumptions concerningthese sources of heterogeneity in detection.Johnson (2008) warned against making additional assump-
tions concerning detection heterogeneity until further inves-tigations identify these sources because it could result in the
Figure 3. Predicted patch-specific occupancy probability for male golden-cheeked warblers in Texas, USA, 2008–2009. Figure from Collier et al. (2012).
Mathewson et al. �Warbler Population Abundance 9
application of an inadequate correction factor and, subse-quently, a false confidence in accuracy. Several factors incombination likely affected detectability of warblers differ-ently across the range in this study and by not systematicallycorrecting for detection we eliminated any concerns regard-ing which sources of variability should be controlled statisti-cally (Diefenbach et al. 2007, Johnson 2008). However, wecollected our data in a manner conducive to future evaluationof some these effects on detection heterogeneity. For exam-ple, the assumption of constant detection could be relaxed byapplying a distance-based or double-observer approach toour data (Nichols et al. 2000, Laake et al. 2011). We furtheraccepted the potential for negative bias from our uncorrecteddensity estimates to produce more conservative estimatesgiven the endangered status of the warbler.We attempted to reduce variation in detectability and to
minimize bias associated with standard avian point countsthrough our study design (Pollock et al. 2002, Johnson2008). We used 5-minute point count periods to minimizedouble counting of birds within the survey area or by count-ing birds that have moved through the area (Buckland 2006).We used a dependent double-observer method because it isrobust to violations of the closure assumption (Moore et al.2004). Although we consistently use well-trained observers,another benefit of the double-observer approach is that itreduces the prospect of species misidentification (Mooreet al. 2004). We assumed a closed population because weconducted single-occasion point counts in each patch. Giventhe distribution of our survey patches, individuals did notlikely move among our survey patches. By conducting single-occasion point counts, we were unable to detect and controlfor any temporal variation, but this was a trade-off to increaseour survey efforts across a greater spatial extent (Thompsonet al. 2002). To minimize temporal variation, we limited oursurveys to the peak of the warbler breeding season, but werecognize that detection probability decreases across thisperiod and we likely negative biased counts as the seasonprogressed (Collier et al. 2010).We report an estimate of the population size of male
warblers during the breeding season, making no assumptionsregarding breeding status of males. In other words, ourcounts included males that were territorial and paired, terri-torial and unpaired, or non-territorial floaters. Previous stud-ies have reported population sizes of total breedingindividuals that required assumptions associated with popu-lation-level pairing success to obtain an estimate of femalepopulation size (Table 1). The bias associated with assump-tions of male breeding status, and consequently pairingsuccess, depends on the study design, survey or monitoringeffort, or adjustments to male-based detections. Becausemany songbird populations, including the golden-cheekedwarbler, consist of territorial but unpaired males (Newton1992), assuming that each male detected during a surveyrepresents a breeding pair overestimates the female propor-tion of the population. Nevertheless, even with informationgarnered from monitoring protocols concerning pairing suc-cess, these studies tend to overlook non-territorial males thatcan constitute a substantial proportion of a population. For
many songbirds, particularly those in fragmented habitats,male breeding status varies with habitat condition and de-mographic factors, such as densities or age structure (Newton1992, Jette et al. 1998, Bayne and Hobson 2001). Althoughstudies that monitor territories can capture variability in malebreeding status at a local scale, these methods are impracticalfor broad-scale application (Buckland 2006). Until furtherassessment of how male breeding success varies across therange, we are unable to determine the degree of precision andbias in reported population sizes adjusted by pairing successestimates from a few locations. In our study, we eliminatedthe need for uninformed adjustments to calculate femalepopulation size and instead we present an approach thatcan incorporate future consideration of male breeding statusacross the range.
Other Habitat Delineations
Population size estimates necessarily depend on the habitatdelineation used to project a local density estimate to therange of potential habitat. Thus, the underlying question ofall attempts to estimate population size is how habitat isdefined and what habitat characteristics can be estimatedusing remotely sensed techniques. Much of the variation inestimated population sizes for the golden-cheeked warblerare created because of differences in data sources used toidentify potential habitat (e.g., Landsat imagery, NLCD),each having differing resolutions and classification accuracies(Table 1). In addition, studies use different characterizationsof the landscape to define habitat, oftentimes assuming aminimum threshold for warbler occurrence or other restric-tions based on evaluations of habitat quality (Table 1). Theseapproaches are inherently conservative, producing negativelybiased habitat and population size estimates. Given theuncertainties in identifying warbler habitat in terms of habi-tat quality, we made no assumption regarding habitat qualityin our broadly defined habitat delineation and provided anestimate of available potential habitat. Indeed, our totalhabitat estimate is similar to other recent studies whenassumptions regarding habitat quality are removed(Diamond and True 1998, Diamond 2007). There may beconcern that overestimating habitat would inflate our popu-lation size estimates, but the additional woodland included inour liberal estimation of habitat was primarily small patchesand edge habitat (Fig. 3; Morrison et al. 2010, Collier et al.2012). Occupancy probability was �0.10 for 59% of thepatches in our habitat delineation (Collier et al. 2012) result-ing in low predicted density estimates for a large portion ofthe habitat; thus, including habitat often assumed as lowerquality contributed little to our total population sizeestimates.Previous population size estimates are from data acquired
on primarily public properties, whereas over half of thehabitat patches in our study were private properties acrossthe warbler’s range. Given the differences in land uses ob-served on our private properties and that there was nosignificant difference in our habitat metrics between privateand public lands, we assumed that warbler density andhabitat conditions did not vary from properties for which
10 The Journal of Wildlife Management
we did not acquire access. Similarly, habitat patches withinwhich landowners actively managed for golden-cheekedwarblers did not bias our sample of public properties. Ofpublic lands, TPWD state parks dominated our sample andfew of these implement management plans for warblers(Groce et al. 2010). Only 5% of our survey patches wereon properties managed specifically for warblers (e.g.,Balcones Canyonland Preserve, Travis County).Our study is an example of how researchers and land
managers can use predictive models to estimate populationabundance. Application of this model is at the range-widescale but implementation of our approach at the regionalscale is possible by incorporating additional information torefine abundance relationships. As remote-sensing technol-ogy evolves, we can examine additional habitat variables thatmight inform within-patch heterogeneity. Furthermore, ourframework provides the opportunity to include detectionprobability to move beyond an index of warbler abundance(Johnson 2008).
MANAGEMENT IMPLICATIONS
Our results indicate that warblers were more abundant inlarger woodland patches of which most occur in the southernportion of the breeding range. Having this knowledge willhelp direct limited resources to the most effective areas forconservation of the species and will help guide managementdecisions (Fitzgerald et al. 2009). The population abundanceestimate we generated can also inform recovery planning forthe warbler. A better understanding of the influence of patchsize and landscape composition on the occupancy and abun-dance of golden-cheeked warblers will provide a tool forrecovery planning. If recovery planning is to involve thedesignation of focal areas for prioritizing conservationefforts, then the model developed here could provide avaluable planning aid in the process. As future habitat con-servation plans and biological assessments are developed, theresults of this work could inform impact assessments forestimating incidental take as well as supporting the develop-ment of more accurate metrics for crediting and mitigationprograms. By providing a framework for more reliably pre-dicting abundance and changes in abundance, our resultsprovide the basis for assessment of this species’ status. Withfurther refinement of the model to account for differences inwithin-patch habitat quality, our results form the basis foradvancing a system for monitoring of the warbler’s abun-dance patterns, ultimately supporting decisions leading torecovery.
ACKNOWLEDGMENTS
The Texas Department of Transportation provided supportfor our work. Additionally, we acknowledge Texas Parks andWildlife Department (TPWD) for support. We appreciatethe private landowners who graciously allowed us access totheir property, as well as TPWD, the U. S. Fish andWildlifeService, The Nature Conservancy, City of Austin, TravisCounty, and the Lower Colorado River Authority. Wethank J. Dunk, P. Hamel, R. Peak, W. Thogmartin, and
F. Thompson III for helpful discussions regarding our re-search. We also thank J. Sedinger, B. Block, C. Farquhar,D. Wolfe, and numerous anonymous referees for insightfulreviews on previous versions of this manuscript. We grate-fully acknowledge A. Snelgrove, K. Skow, B. Stevener, andA. Dube for logistical support, as well as the many techni-cians and graduate students from Texas A&M University.B. A. Collier acknowledges partial support from Award No.KUS-C1-016-04 given by the King Abdullah University ofScience and Technology (KAUST).
LITERATURE CITEDBaccus, J. T., M. E. Tolle, and J. D. Cornelius. 2007. Response of golden-cheeked warblers (Dendroica chrysoparia) to wildfires at Fort Hood, Texas.Texas Ornithological Society. Occasional Publication 7:1–37.
Bayne, E. M., and K. A. Hobson. 2001. Effects of habitat fragmentation onpairing success of ovenbirds: importance of male age and floater behavior.Auk 118:380–388.
Benson, R. H. 1990. Habitat area requirements of the golden-cheekedwarbler on the Edwards Plateau. Texas Parks and WildfireDepartment, Austin, USA.
Boyce, M. S., and L. L. McDonald. 1999. Relating populations to habitatsusing resource selection functions. Trends in Ecology & Evolution14:268–272.
Buckland, S. T. 2006. Point transect surveys for songbirds: robust method-ologies. Auk 123:345–357.
Burnham, K. P., and D. R. Anderson. 2002. Model selection and multi-model inference: a pratical information-theoretic approach. Springer, NewYork, New York, USA.
Butcher, J. A., M. L. Morrison, D. Ransom, Jr., R. D. Slack, and R. N.Wilkins. 2010. Evidence of a minimum patch size threshold of reproduc-tive success in an endangered songbird. Journal of Wildlife Management74:133–139.
Campbell, S. P., J. A. Clark, L. H. Crampton, A. D. Guerry, L. T. Hatch,P. R. Hosseini, J. J. Lawler, and R. J. O’Connor. 2002. An assessment ofmonitoring efforts in endangered species recovery plans. EcologicalApplications 12:674–681.
Campomizzi, A. J., J. A. Butcher, S. L. Farrell, A. G. Snelgrove, B. A.Collier, K. J. Gutzwiller, M. L. Morrison, and R. N. Wilkins. 2008.Conspecific attractions is a missing component in wildlife habitat model-ing. Journal of Wildlife Management 72:331–336.
Collier, B. A., J. E. Groce, M. L.Morrison, C. Newnam, A. J. Campomizzi,S. L. Farrell, H. A. Mathewson, R. T. Snelgrove, R. J. Carroll, and R. N.Wilkins. 2012. Predicted patch occupancy in fragmented landscapes at therangewide scale for an endangered species: an example of an Americanwarbler. Diversity and Distributions 18:158–167.
Collier, B. A., M. L. Morrison, S. L. Farrell, A. J. Campomizzi, J. A.Butcher, K. B. Hays, D. I. MacKenzie, and R. N. Wilkins. 2010.Monitoring endangered species occupying private lands: case study usingthe golden-cheeked warbler. Journal of Wildlife Management 74:140–147.
Conroy, M. J., J. P. Runge, R. J. Barker, M. R. Schofield, and C. J.Fonnesbeck. 2008. Efficient estimation of abundance for patchily distrib-uted populations via two-phase, adaptive sampling. Ecology 89:3362–3370.
Cook, R. D., and J. O. Jacobson. 1979. A design for estimating visibility biasin aerial surveys. Biometrics 35:735–742.
DeBoer, T. S., and D. D. Diamond. 2006. Predicting presence-absenceof the endangered golden-cheeked warbler (Dendroica chrysoparia).Southwestern Naturalist 51:181–190.
Diamond, D. D. 2007. Range-wide modeling of golden-cheeked warblerhabitat. Project Final Report to Texas Parks & Wildlife, Austin, Texas,USA.
Diamond, D. D., and C. D. True. 1998. Golden-cheeked warbler nestinghabitat area, habitat distribution, and change. Project Final Report toUSFWS, Region 2, Albuquerque, New Mexico, USA.
Diefenbach, D. R., M. R. Marshall, J. A. Mattice, and D. W. Brauning.2007. Incorporating availability for detection in estimates of bird abun-dance. Auk 124:96–106.
Mathewson et al. �Warbler Population Abundance 11
Original Article
Utilization of a Species Occupancy Model forManagement and Conservation
TIFFANY M. MCFARLAND,1 Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
HEATHER A. MATHEWSON, Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
JULIE E. GROCE, Institute of Renewable Natural Resources, Texas A&M University, San Antonio, TX 78215, USA
MICHAEL L. MORRISON, Department of Wildlife and Fisheries Sciences, Texas A&M University, College Station, TX 77843, USA
J. CAL NEWNAM, Texas Department of Transportation, P.O. Box 15426, Austin, TX 78761, USA
R. TODD SNELGROVE, Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
KEVIN L. SKOW, Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
BRET A. COLLIER, Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
R. NEAL WILKINS, Institute of Renewable Natural Resources, Texas A&M University, College Station, TX 77843, USA
ABSTRACT Conserving habitat is increasingly challenging as human populations grow. Remote-sensingtechnology has provided a means to delineate species’ habitat on large spatial scales. However, by combininghabitat delineations with predictions of species’ occurrence, habitat models can provide additional utilityapplications for conservation by allowing us to forecast how changing environmental and landscape con-ditions affect species’ occurrence and distribution.We demonstrate how a spatially explicit habitat occupancymodel for the golden-cheeked warbler (Setophaga chrysoparia) can be used as an impact assessment andconservation planning tool. We used predictions of patch-level occupancy rates and created several scenariosthat simulated the removal or protection of warbler habitats. Resulting changes to habitat structure andavailability were used to assess the resulting impacts of removal or protection on the occurrence probability forremaining habitat patches. By recalculating occupancy based on changes to habitat, our approach provides theability to assess and compare impacts of location and orientation of development so that the least harmfuloption relative to predicted occurrence can be chosen. Potential applications of our modeling approach aremany because our methods provide a useful tool for identifying potential impacts and assisting withmitigation efforts focused on the conservation and management of a species. � 2012 The Wildlife Society.
KEY WORDS conservation planning, golden-cheeked warbler, habitat, landscape management, mitigation, Setophagachrysoparia, species distribution.
Implementation of successful management strategies requireknowledge of animal distribution relative to environmentalconditions and predictions of how animal populations willrespond to direct and indirect environmental impacts(Kantrud and Stewart 1984, Wiens and Rotenberry 1985,Debinski and Brussard 1994, Colwell and Dodd 1995).Spatially explicit habitat models that predict the probabilityof species occurrence on broad spatial scales are useful toolsfor assessing potential impacts of land use or environmentalchanges, or for forecasting the outcome of managementactions (Sagarin et al. 2006, Brotons et al. 2007). For speciesof conservation concern, the implications of habitat loss areoften significant and loss can have long-term implications onpopulation sustainability.The golden-cheeked warbler (Setophaga chrysoparia; here-
after, ‘‘warbler’’) is a Neotropical migratory songbird thatbreeds exclusively in central Texas, USA, in mixed wood-
lands of mature Ashe juniper (Juniperus ashei), oak (Quercusspp.), and other deciduous species (Pulich 1976, Wahl et al.1990). In 1990, the U.S. Fish andWildlife Service (USFWS)listed the warbler as endangered, citing habitat loss andfragmentation as the primary threats to the species(USFWS 1990). Increasing infrastructure development(e.g., roads, utilities) within the warbler’s limited breedingrange creates networks of corridors that can deterioratewoodland habitat due to fragmentation (Peak 2007) andpotentially inhibit exchange and dispersal of individualsbetween patches (Gobeil and Villard 2002, Goodwin andFahrig 2002, Haynes and Cronin 2006, Richard andArmstrong 2010). However, the impacts of environmentalor anthropogenic activities that could fragment potentialhabitat have typically been limited by our knowledge ofthe likelihood of warbler occurrence within those habitats.Any loss of habitat, either from natural phenomena (i.e., fire)or human activities, has the potential to impact the warblerdue to the species’ limited range.Several warbler habitat classifications created in the past
20 years have provided estimates of range-wide habitat
Received: 10 August 2011; Accepted: 14 November 2011
McFarland et al. � Utilization of a Species Occupancy Model 1
distribution (Wahl et al. 1990, Rowell et al. 1995, Rappoleet al. 2003, DeBoer and Diamond 2006, Loomis-Austin,Inc. 2008) based on remotely sensed metrics and researcher-defined descriptions of warbler habitat. Several previoussurveys also evaluated warbler occurrence and habitat asso-ciations, but these studies focused on a few public propertieswith known populations of warblers (Anders and Dearborn2004, Peak 2007, Reidy et al. 2008), reducing the ability ofthe authors to extrapolate avian–habitat relationships acrossthe breeding range. Two studies classified habitat qualityacross the range. DeBoer and Diamond (2006) used sampledata from across the range to model habitat quality.However, their sample size was small (n ¼ 50 patches),and several habitat metrics identified as relevant cannot beassessed remotely and, therefore, cannot be modeled ormapped at the range-wide scale. Diamond (2007) modeledhabitat quality across the range, but estimates of quality werebased on researchers’ evaluation of quality, rather than usingwarbler occurrence data collected in the field to drive thedevelopment of a habitat model.Here we demonstrate how a spatially explicit habitat model
can be utilized as a tool for conservation planning and impactassessment. We used a patch-level habitat occupancy modelfor the warbler built on survey data collected on public andprivate properties across the warbler’s breeding range inTexas (Collier et al. 2012). These patch occupancy predic-tions provided high accuracy in predicting warbler occur-rence at the patch scale and were based on patch size andlandscape composition (hereafter, ‘‘woodland composition’’)metrics combined with a north to south spatial gradient.Unlike previous models, which only delineate potential hab-itat, Collier et al. (2012) can categorize a habitat patch basedon the probability of occupancy and, thus, can be used topredict the impacts of habitat loss or modification on warbleroccupancy, rather than simply calculating the amount ofhabitat change. Herein, we illustrate the use of such a modelby first showing the changes in predicted warbler occupancybased on simulated fragmentation or removal of habitatpatches, and then showing how these effects can provideguidance for warbler conservation and management.Although we were able to sample across the entirety ofthe warbler’s breeding range, occupancy models such asours could also be developed for subsets of the range ofmost species.
METHODS
Occupancy ModelCollier et al. (2012) created a spatially explicit occupancymodel for the warbler across its range in Texas. Collier et al.defined potential warbler habitat as mixed woodland patchesdelineated on 2007 and 2008 LANDSAT 5 imagery (Fig. 1;Collier et al. 2010). Collier et al. eliminated nonwoodlandland-cover types (e.g., wetlands, cropland, urban areas,water) using the 2001 National Land Cover Data set(Homer et al. 2007), as well as pixels misclassified as wood-land (canopy cover<30%). Collier et al. identified pixels thatintersected paved or public roads using the Texas strategic
mapping program, and removed those pixels from woodlandclassification, which defined breaks between some patches.As such, this delineation resulted in 63,616 distinct patchesof potential habitat (mean patch size ¼ 26.39 ha, range ¼2.8–26,967 ha), equaling approximately 1.678 million ha,which were the sample units for the model (Fig. 1;Collier et al. 2010, 2012). Collier et al. used ArcGIS9.3.1 to calculate patch size (ha) and woodland composition(percent of woodland within a 400-m radius of each pixel inthe patch, averaged over the entire patch) for each patch ofmixed woodland identified within the study area (modifiedfrom Magness et al. 2006). The model included these 2metrics, plus an interaction between them, as well as spatiallocation (latitude and longitude) to assign an occurrenceprobability to each of the 63,616 patches of potential habitat.Predictions from Collier et al. (2012) can be visualizedas a map depicting each habitat patch with an associatedprobability of occupancy, ranging from 0 to 1. We used thesepatch-specific estimates of warbler occurrence probabilityfromCollier et al. (2012) as the foundation for our evaluationof natural and human-induced impacts on the golden-cheeked warbler across its breeding range in Texas.
Modeling Impacts
We used 2 hypothetical (Scenarios 1 and 2) and 2 real-life(Scenarios 3 and 4) scenarios of development activities andnatural impacts to demonstrate the utility of the modelfor predicting impacts of habitat loss and fragmentation.Scenario 1 demonstrates how land managers can use themodel predictions to determine where to protect warblerhabitat. Scenario 2 details potential impacts on warblerhabitat due to residential development outside of an urbanarea. Scenario 3 evaluates the impact of a large wildfire that
Figure 1. The 35-county breeding range of the golden-cheeked warbler(USFWS 1992) and the mixed-woodland habitat patches that create thepotential habitat layer used in the model (Collier et al. 2012). Distinctpatches of this cover type were used as the sample units.
2 Wildlife Society Bulletin
occurred in April 2011 within the northern range of thewarbler, while Scenario 4 demonstrates how model predic-tions can support selection of the best route for a powertransmission corridor that minimizes impacts on habitatpatches.We used ArcGIS 9.3.1 to remove habitat from our delin-
eation of habitat specific to Scenarios 2–4 and to evaluate thearea of patches occurring within given boundaries. ForScenario 2, we deleted habitat patches that mimicked actualsizes and shapes of other developments in the region but donot represent actual locations of planned developments. ForScenarios 3 and 4 we obtained location information from anApril 2011 wildfire and for a transmission corridor proposedby a utility company in Texas (Lower Colorado RiverAuthority 2010). After removing habitat, we recalculatedpatch area and woodland composition for each scenario usingthe methods detailed in Collier et al. (2012). We thenrecalculated patch-specific occupancy probabilities for allpatches using the new patch area and woodland compositionvalues to evaluate the indirect impacts of habitat loss. Wedid not consider patch fragments �2 ha as potential habitatbased on the approximate mean territory size of warblers(2.9 ha; Pulich 1976) and assigned them an occupancyprobability of zero.
MODEL UTILITY
Model BasicsThe underlying occupancy model calculates predicted occu-pancy based on patch size, patch-specific woodland compo-sition, and space; thus, some general trends in predictedpatch occupancy emerged (Collier et al. 2012). First, asexpected, larger patches generally have higher predictedoccupancy rates than smaller patches. Second, because themodel accounts for woodland composition, an increase ordecrease in the amount of woodland in one habitat patch willaffect the predicted occupancy of all other patches locatedwithin 400 m of the altered patch; therefore, patches that aremore isolated from other patches of woodland habitat typi-cally had lower occupancy rates (Fig. 2). Third, patch occu-pancy predictions and projections from Collier et al. (2012)incorporated spatial location across the range; thus, patchesin the southern part of the range generally have higheroccupancy probabilities than those in the north. Based onthe characteristics included in the model, we can use themodel to recalculate occupancy probability when habitatpatches are lost or fragmented. Fragmentation of patchesresults not only in habitat loss (i.e., the habitat removedwithin the corridors that fragment the patch), but alsoreduces the size of patches and the woodland compositionof surrounding patches, causing a decrease in probability ofoccupancy of the patch fragments (as illustrated below).
ScenariosScenario 1: Prioritizing regions for conservation.—Habitat
delineation alone provides a minimal approach to determin-ing levels of conservation importance among patchesand properties because it ignores potential differences inoccupancy probabilities among patches within property
boundaries. Without assigning occupancy probabilities tohabitat, a habitat distribution map would assume equalwarbler occupancy to 2 properties (property A and B inFig. 3a) that have approximately equal amounts of potentialhabitat (250.6 ha and 254.1 ha, respectively). Using thewarbler occupancy model, property A has higher overallwarbler occupancy probability, with 220.2 ha (88% of total
Figure 2. Predicted patch occupancy from Collier et al. (2012) for 3 patchesof similar size (4.9–5.0 ha) and in the same geographic location (KerrCounty, TX, USA), but with different landscape composition due to theproximity of other patches. Patches of the same size that are less isolatedhave a higher landscape composition and, therefore, a higher probability ofoccupancy.
Figure 3. Demonstration of predicted occupancy wherein (a) 2 properties,A and B, contain approximately equal amount of potential golden-cheekedwarbler habitat (250.6 ha and 254.1 ha, respectively) designated by greenpolygons, but (b) predicted patch occupancy from Collier et al. (2012)differed among patches in each property because of variation in patch-specific sizes and landscape composition. Warm colors in (b) representhigher predicted patch occupancy probabilities.
McFarland et al. � Utilization of a Species Occupancy Model 3
area) above 93% occupancy probability, than does property B(Fig. 3b). Property B only has 75.7 ha (30% of total area)above 70% probability of occupancy, and 176.8 ha above 50%probability of occupancy.Scenario 2: Urban expansion.—Urban expansion is one of
the major threats to woodland and, thus, warbler habitat(USFWS 1990); therefore, our next example evaluatedimpacts of potential urban development by simulating thefragmentation of a large habitat patch near San Antonio,Texas. We chose this patch due to its size, high occupancyprobability, and relative proximity to a major city. We sim-ulated removal of habitat by residential development byreplicating the shape and size of neighborhoods, shoppingcenters, and connecting roads in the area and placing themover potential warbler habitat, delineated by our habitatclassification. We determined the amount of habitat lostand the changes in probability of warbler occurrence usingthe occupancy model. Before the simulation, the total habitatwithin the patch was 5,104.8 ha with a predicted occupancyprobability of 0.99 (Fig. 4). We simulated 1,093.2 haof developed area, a 21.4% reduction in total habitatpostdevelopment, which fragmented the original patchinto 20 individual patches. The largest remaining patchwas 1,139.6 ha and still had an occupancy probability of0.97. Fifteen patches (1,095.1 ha) had an occupancy proba-bility <0.9, of which 8 patches had an occupancy probability<0.8 (290.3 ha).Scenario 3: Wildfire.—In April 2011, a major wildfire
started west of Possum Kingdom Lake in StephensCounty, Texas, and ultimately burned approximately51,000 ha of land in Stephens, Palo Pinto, and Eastland
counties (Fig. 5). Stephens and Palo Pinto counties comprise2 of the 3 counties that make up the USFWS RecoveryRegion 1 for the warbler (Fig. 5). Because wildfires burn atdifferent intensities across the burn area, we cannot assumeall area within the fire scar was lost warbler habitat. However,we evaluated the impact of this fire on warbler habitat interms of maximum potential habitat lost by running thesimulation as if all habitat within the fire scar was lost,and addressing the potential effect on nearby patches andresulting patch fragments based on changes in probability ofoccupancy.Prior to the wildfire, Recovery Region 1 contained approx-
imately 114,550 ha of potential habitat, with 38,180 ha(33.3%) of the potential habitat having >50% probabilityof occupancy (Table 1). The fire caused a potential loss of19,994 ha of habitat (17.5% of total habitat in the region);however, nearly 40% of the habitat with occupancy proba-bility between 0.7 and 0.9 was potentially lost or affected bythe fire. The region did not contain habitat with an occu-pancy probability >0.9 (Table 1).Scenario 4: Transmission line corridor.—In September 2010,
our habitat model (Collier et al. 2012) was used to assess theroute of the Competitive Renewable Energy Zone (CREZ)transmission line from northern Schleicher County to thewest side of Kendall County, Texas (Wilkins 2010, directtestimony, State Office of Administrative Hearings, SOAHDocket no. 473-10-5546, 28 Sep 2010, Austin, TX). Weevaluated the impact of the planned 48.8-m-wide corridortransmission line and 2 proposed route alternatives, referred
Figure 4. Demonstration of the effect of fragmentation of a patch of golden-cheeked warbler habitat. In this hypothetical example, we selected an areacontaining a large (5,129-ha) patch of potential habitat in Bexar County,Texas, USA (a), with a predicted probability of warbler occupancy of 0.99,based on the Collier et al. (2012) occupancy model. Warm colors representhigher predicted estimates of warbler occupancy and numbers are the patch-specific probability of occupancy. After 3,923 ha of potential habitat wereremoved from the area, the largest remaining patch was only 1,085 ha (b),but still had a 0.955 predicted occupancy. The other patches created hada predicted occupancy ranging from 0 (for the smallest patches) to 0.909.
Figure 5. Extent and occupancy probability of golden-cheeked warblerhabitat patches potentially impacted by the Possum Kingdom (PK),Stephens, Palo Pinto, and Eastland counties (TX, USA) complex fire inApril 2011. The fire potentially eliminated 17.5% of the potential habitatavailable within Recovery Region 1.
4 Wildlife Society Bulletin
to here as the northern route and southern route, with theintention of demonstrating which route option would have alower impact on warbler habitat. The southern route followsInterstate 10 for much of its length while the northern routecuts a shorter, more direct path (Fig. 6). For this analysis,we removed all potential warbler habitat within the 48.8-m-wide transmission corridors for the length of both routes andcalculated total habitat lost. Additionally, if patch fragmentsresulting from the development were �2 ha, we assumedthat warblers would not occupy habitat patches of this size.Thus, we included the size of these patch fragments in theestimate of total habitat lost for each transmission lineoption. We then assessed the indirect impacts of the routesby recalculating patch occupancy probabilities for all patchesthat occurred within, or were intersected by, the transmissioncorridors.The northern route traversed 18.5 km of potential habitat
and would eliminate 87.5 ha. The southern route traversed27.0 km of potential habitat with a total potential habitatloss of 125.1 ha. Fragments smaller than 2 ha accounted for
142 patches along the northern route, totaling 53 ha ofpotential habitat, and these <2-ha patches accounted for176 patches along the southern route, totaling 54 ha ofpotential habitat. Combining these values with the amountof habitat that would be directly lost within the corridormeant that approximately 141 ha of potential habitat wouldbe lost in the preferred route, and approximately 180 hawould be lost in the southern route, initially indicatingthat the preferred route would impact less of the warblerhabitat.Looking at only the patches that intersected each route, our
assessment of habitat lost based on recalculated occupancyprobabilities showed that habitat categorized as >0.9 occu-pancy probability was reduced by a greater amount in thenorthern route (487 ha) than the southern route (100 ha;Table 2; Fig. 7) due to a decrease in occupancy probabilityof the patches that remained after being fragmented bythe corridors. However, the comparison changes dependingon the occupancy probabilities of interest. For instance, thehectares lost in all probability categories>0.6 actually show agreater amount of loss in the southern route (385 ha) thanthe northern route (164 ha), but the hectares lost in proba-bility categories >0.7 are approximately equal between the2 routes.
DISCUSSION
Occupancy models allow land managers to identify not onlythe amount of habitat lost from proposed or actual activitiesor natural events, but also the probability that the habitatlost would be occupied by the species of interest. Further,the potential changes in probability of occupancy of nearbypatches and any remaining fragments of patches can becalculated and mapped. As we demonstrated here usingthe golden-cheeked warbler, the ability to adapt the habitatmodel to realized and potential changes in habitat and thenmodel the associated changes in occupancy probability with-in surrounding patches creates endless conservation applica-tions. In our scenarios, we demonstrated how estimating onlytotal habitat can result in misleading or incorrect conclusionsregarding conservation value of different areas or assessmentof fragmentation and other impacts on warbler habitat.
Table 1. Summary of total golden-cheeked warbler habitat potentially lost within U.S. Fish andWildlife Service Recovery Region 1 resulting from a wildfire inApril 2011 west of PossumKingdomLake, Stephens County, Texas, USA. The amount of habitat falling within the fire scar accounted for 17.5% of the habitatin the region; however, nearly 40% of the habitat with occupancy probability >0.7 was potentially lost.
Figure 6. Two of many proposed routes for the CREZ transmission linefrom Schleicher County to Kendall County in central Texas, USA. Thenorthern route is shown in blue and the southern route that parallelsInterstate 10 for the majority of its length is shown in red.
McFarland et al. � Utilization of a Species Occupancy Model 5
Additionally, assuming that the relationship between occur-rence and density of warblers is similar to other species(Royle and Nichols 2003, Royle 2004), which is indicatedby the findings of Mathewson et al. (2012), a decrease inpredicted occupancy may have more impact than just adecrease in the probability that the patch will be occupiedby a single bird. Implementation of the warbler occupancymodel can guide land managers in selecting advantageouslocations for conservation or management activities, evalu-ating least-cost configurations, locations, and pathwaysof developments or corridors, and assessing the impact ofnatural or human-induced habitat loss.We presented a scenario in which assessments of occupancy
probability can guide land managers in determining advan-tageous locations for conservation. In Scenario 1 we dem-onstrated how a planner or manager could determine whereto prioritize conservation efforts. Without a calculatedoccupancy probability for each patch, the 2 areas inour scenario would appear to be of equal conservation
importance; however, one property contained habitat patcheswith higher occupancy probability estimates. It is advanta-geous for land managers to seek optimal locations for con-servation or habitat management practices. Many songbirds,including the golden-cheeked warbler, demonstrate high sitefidelity on the breeding grounds in that breeding adults aremore likely to return to areas or adjacent to areas that theybreed at in previous years (Holmes and Sherry 1992, Jetteet al. 1998). Furthermore, golden-cheeked warblers demon-strate conspecific attraction, wherein breeding birds use thepresence of conspecifics to determine locations to settle forbreeding (Campomizzi et al. 2008, Farrell 2011). Theseecological phenomena of songbirds suggests that occupiedareas are more likely to harbor birds in subsequent years;thus, conservation of habitat is likely more beneficial in areaswith higher occupancy probabilities.We also demonstrated how the model can assess the impact
of habitat loss and subsequent fragmentation of habitat,either after the loss of habitat has occurred (Scenario 3) orpreemptively for proposed developments (Scenarios 2 and 4).In Scenario 3, we show that habitat loss due to a wildfireappeared to eliminate only 17.5% of the available habitat inRecovery Region 1, but it potentially eliminated nearly 40%of the habitat with a probability of occupancy >0.7.Therefore, the impact of this fire on the region may havea great impact on the recovery of the species, especiallybecause recovery of the species defined by USFWS is depen-dent upon recovery in all individual recovery regions.In Scenario 2, we demonstrate how fragmentation of a large
patch can reduce the occupancy probabilities of resultingpatch fragments. Although fragmentation of habitat is sel-dom desirable for the warbler, fragmentation of large patchesmay not always decrease the occupancy probability of result-ing patches. For instance, if habitat fragmentation occurs in alarge area of habitat and results in large patch fragments(�160 ha), the predicted occupancies of the new patchesmay not decrease significantly from the original occupancyvalue, supporting the idea of a patch-size threshold (Butcheret al. 2010, Collier et al. 2010). We demonstrated this in
Table 2. Count of patches and total hectares within 10 golden-cheeked warbler occupancy probability classes both pre- and postimpact for 2 alternateCompetitive Renewable Energy Zone routes in Texas, USA. Although the north route only directly eliminates 87.4 ha of habitat, 487.6 ha having a predictedoccupancy >0.9 are either directly lost or result in a lowered predicted occupancy due to fragmentation. However, looking at the change in patches withoccupancy probability >0.7, or even patches with occupancy probability >0.5, the 2 routes are similar.
Occupancyprobability
North South
No. of patches Sum (ha) Change(ha)
No. of patches Sum (ha) Change(ha)Preimpact Postimpact Preimpact Postimpact Preimpact Postimpact Preimpact Postimpact
Figure 7. Change in hectares in each golden-cheeked warbler occupancyprobability category for the preferred route and an alternate southern route inTexas, USA. More habitat is lost in the highest occupancy-probability cate-gory in the preferred route, indicating that the preferred route will have moreof an impact on potential warbler habitat than the longer southern route.
6 Wildlife Society Bulletin
Scenario 2, wherein the predicted patch occupancy of thelargest patch was not significantly reduced even thoughconsiderable habitat was removed. Although complete con-servation of the largest patches may often be unreasonabledue to financial, social, and logistical limitations, our findingssuggest that the conservation of suitably sized and strategi-cally placed patches may be more efficient for maintainingviable habitat (Zuidema et al. 1996).In Scenario 4, the proposed CREZ routes illustrate how
including information on occupancy probabilities whenprioritizing location options for development provides great-er insight into impacts to warbler habitat beyond simplecalculations of the total amount of habitat that is removedfor development. Based only on total habitat area, the north-ern route for the CREZ line would remove less total habitatthan the southern route; however, total habitat area fails toaccount for differences in the sizes and distribution of habitatpatches. In our example, the southern route followed theexisting interstate corridor and habitat along this route wasalready fragmented into smaller sized patches, resulting inlower occupancy probabilities of these patches. Althoughmore total habitat would be lost along the southern route,less high-occupancy habitat (>0.9) would be removed com-pared to habitat along the northern route. The northernroute would bisect multiple large and contiguous patchesof warbler habitat, increasing fragmentation of these patchesand removing more high-occupancy habitat patches.However, both routes eliminated approximately the sameamount of habitat with an occupancy probability >0.7. Thisexercise demonstrated the fallacy of evaluating impacts onhabitat when considering only total habitat loss estimations.Similarly, the model could be used for assessing impacts ofdevelopment of residential neighborhoods, shopping centers,infrastructure like roads and transmission lines, clearing forranching, or natural disasters in which habitat would be lost.Two main trends resulting from our model are demon-
strated by our scenarios. First, concurring with previous work(Coldren 1998, DeBoer and Diamond 2006, Baccus et al.2007, Butcher et al. 2010), our model showed the importanceof large patches for maintaining high rates of warbler occu-pancy. Large patches harbor more potential habitat and alsotend to have a higher predicted occupancy. Thus, there is lessimpact on predicted occupancy when habitat is removedalong the edges of patches as opposed to bisecting patches.Second, the model showed that isolated patches had a lowerprobability of occupancy, indicating the importance of main-taining connectivity of habitat (through close proximityof patches) across the landscape, particularly in areas wherehabitat patches are small.Although we were able to calculate changes in predicted
occupancy on any spatial scale, as currently formulated, ourmodel is most appropriately applied to broader spatial extents(e.g., network of patches) as opposed to a patch-by-patch orwithin-patch basis, and for determining the potential changein habitat area and occupancy between multiple managementscenarios. As for any model, however, additional field sam-pling will allow for improvements in the underlying database, as will field testing to continually improve model
predictions, especially as applied to increasingly smallerspatial extents (Vaughan and Ormerod 2003).Given the availability of 1) a range-wide map of a species’
habitat, 2) remotely sensed, species-relevant habitat data,and 3) occupancy data from randomly sampled units acrossthe habitat, our modeling technique can be applied to anyspecies. Unlike static habitat maps, models like these havethe dynamic ability to project the probability of occupancythroughout the range and also adapt to changes in thehabitat. The dynamic nature of these models providesnumerous applications for conservation of habitat andassessment of actual or potential impacts.
ACKNOWLEDGMENTS
We thank the Texas Department of Transportation forsupporting this work, and the many private landownerswho graciously allowed us access to their property, as wellas Texas Parks and Wildlife, the U.S. Fish and WildlifeService, The Nature Conservancy, City of Austin, TravisCounty, and the Lower Colorado River Authority for prop-erty access. We gratefully acknowledge A. Snelgrove, B.Stevener, and A. Dube for logistical support, as well asthe many technicians and graduate students from TexasA&M University. We thank Dr. Y. Wiersma ofMemorial University for reviewing this manuscript.
LITERATURE CITEDAnders, A. D., and D. C. Dearborn. 2004. Population trends of theendangered golden-cheeked warbler at Fort Hood, Texas, from 1992–2001. The Southwestern Naturalist 49:39–47.
Baccus, J. T., M. E. Tolle, and J. D. Cornelius. 2007. Response of golden-cheeked warblers (Dendroica chrysoparia) to wildfires at Fort Hood, Texas.Texas Ornithological Society, Occasional Publication 7:1–37.
Brotons, L., S. Herrando, and M. Pla. 2007. Updating bird species distri-bution at large spatial scales: applications of habitat modeling to data fromlong-term monitoring programs. Diversity and Distributions 13:276–288.
Butcher, J. A., M. L. Morrison, D. Ransom, Jr., R. D. Slack, and R. N.Wilkins. 2010. Evidence of a minimum patch size threshold of repro-ductive success in an endangered songbird. The Journal of WildlifeManagement 74:133–139.
Campomizzi, A. J., J. A. Butcher, S. L. Farrell, A. G. Snelgrove, B. A.Collier, K. J. Gutzwiller, M. L. Morrison, and R. N. Wilkins. 2008.Conspecific attraction is a missing component in wildlife habitat model-ing. The Journal of Wildlife Management 72:331–336.
Coldren, C. L. 1998. The effects of habitat fragmentation on the golden-cheeked warbler. Dissertation, Texas A&M University, College Station,USA.
Collier, B. A., J. E. Groce, M. L. Morrison, J. C. Newnam, A. J.Campomizzi, S. L. Farrell, H. A. Mathewson, R. T. Snelgrove, R. J.Carroll, and R. N. Wilkins. 2012. Predicting patch occupancy in frag-mented landscapes at the range-wide scale for an endangered species: anexample of an American warbler. Diversity and Distributions 18:158–167.
Collier, B. A., M. L. Morrison, S. L. Farrell, A. J. Campomizzi, J. A.Butcher, K. B. Hays, D. I. MacKenzie, and R. N. Wilkins. 2010.Monitoring endangered species occupying private lands: case study usingthe golden-cheeked warbler. The Journal of Wildlife Management 74:1–12.
Colwell, M. A., and S. L. Dodd. 1995. Water bird communities and habitatrelationships in coastal pastures of Northern California. ConservationBiology 9:827–834.
Debinski, D. M., and P. F. Brussard. 1994. Using biodiversity data to assessspecies–habitat relationships in Glacier National Park, Montana.Ecological Applications 4:833–843.
McFarland et al. � Utilization of a Species Occupancy Model 7
DeBoer, T. S., and D. D. Diamond. 2006. Predicting presence–absenceof the endangered golden-cheeked warbler (Dendroica chrysoparia). TheSouthwestern Naturalist 51:181–190.
Diamond, D. D. 2007. Range-wide modeling of golden-cheeked warblerhabitat. Technical report, Texas Parks and Wildlife Department, Austin,Texas, USA.
Farrell, S. L. 2011. Use of social information for habitat selectionin songbirds. Dissertation, Texas A&M University, College Station,USA.
Gobeil, J.-F., and M.-A. Villard. 2002. Permeability of three boreal forestlandscape types to bird movements as determined from experimentaltranslocations. Oikos 98:447–458.
Goodwin, B. J., and L. Fahrig. 2002.How does landscape structure influencelandscape connectivity? Oikos 99:552–570.
Haynes, K. J., and J. T. Cronin. 2006. Inter-patch movement and edgeeffects: the role of behavioral responses to the landscape matrix. Oikos113:43–54.
Holmes, R. T., and T.W. Sherry. 1992. Site fidelity of migratory warblers intemperate breeding and Neotropical wintering areas: implications forpopulation dynamics, habitat selection, and conservation. Pages 563–575 in J. M. Hagan and D. W. Johnston, editors. Ecology and conserva-tion of Neotropical migrant land birds. Smithsonian Institution Press,Washington, D.C., USA.
Homer, C., J. Dewitz, J. Fry, M. Coan, N. Hossain, C. Larson, N. Herold,A. McKerrow, J. N. VanDriel, and J. Wickham. 2007. Completion of the2001 national land cover database for the conterminous United States.Photogrammetric Engineering and Remote Sensing 73:337–341.
Jette, L. A., T. J. Hayden, and J. D. Cornelius. 1998. Demographics of thegolden-cheeked warbler (Dendroica chrysoparia) on FortHood, Texas. U.S.ArmyCorps of Engineers, U.S. ArmyConstruction Engineering ResearchLaboratory Technical Report 98/52, Fort Hood, Texas, USA.
Kantrud, H. A., and R. E. Stewart. 1984. Ecological distribution and crudedensity of breeding birds on prairie wetlands. The Journal of WildlifeManagement 48:426–437.
Lower Colorado River Authority. 2010. Lower Colorado River Authorityand Competitive Renewable Energy Zones. <www.lcra.org/CREZ>Accessed 1 Sep 2010.
Loomis-Austin, Inc. 2008. Mapping potential golden-cheeked warblerbreeding habitat using remotely sensed forest canopy cover data.County of Hays, San Marcos, Texas, USA.
Magness, D. R., R. N. Wilkins, and S. J. Hejl. 2006. Quantitative relation-ships among golden-cheeked warbler occurrence and landscape size,composition, and structure. Wildlife Society Bulletin 34:473–479.
Mathewson, H. A., J. E. Groce, T. M. McFarland, M. L. Morrison, J. C.Newnam, R. T. Snelgrove, B. A. Collier, and R. N. Wilkins. 2012.Estimating breeding season abundance of golden-cheeked warblers inTexas. Journal of Wildlife Management 76:in press.
Peak, R. G. 2007. Forest edges negatively affect golden-cheeked warblernest survival. The Condor 109:628–637.
Pulich, W. M. 1976. The golden-cheeked warbler: a bioecological study.Texas Parks and Wildlife Department, Austin, USA.
Rappole, J. H., D. I. King, and J. Diez. 2003. Winter- vs. breeding-habitatlimitation for an endangered avian migrant. Ecological Applications13:735–742.
Reidy, J. L., M. M. Stake, and F. R. Thompson, III. 2008. Golden-cheekedwarbler nest mortality and predators in urban and rural landscapes. TheCondor 110:458–466.
Richard, Y., and D. Armstrong. 2010. The importance of integratinglandscape ecology in habitat models: isolation-driven occurrence of northisland robins in a fragmented landscape. Landscape Ecology 25:1363–1374.
Rowell, G. A., D. P. Keddy-Hector, D. D. Diamond, J. Lloyd, B.McKinney, R. C. Maggio, and T. L. Cook. 1995. Remote sensing andGIS of golden-cheeked warbler breeding habitat and vegetation types onthe Edwards Plateau. Texas Parks and Wildlife Department, Austin,USA.
Royle, J. A. 2004. Modeling abundance index data from anuran callingsurveys. Conservation Biology 18:1378–1385.
Royle, J. A., and J. D. Nichols. 2003. Estimating abundance from repeatedpresence–absence data or point counts. Ecology 84:777–790.
Sagarin, R. D., S. D. Gaines, and B. Gaylord. 2006. Moving beyondassumptions to understand abundance distributions across the ranges ofspecies. Trends in Ecology & Evolution 21:524–530.
United States Fish and Wildlife Service [USFWS]. 1990. Endangered andthreatened wildlife and plants; final rule to list the golden-cheeked warbleras endangered. Federal Register 55:53153–53160.
United States Fish and Wildlife Service [USFWS]. 1992. Golden-cheekedwarbler (Dendroica chrysoparia) recovery plan. U.S. Fish and WildlifeService, Austin, Texas, USA.
Vaughan, I. P., and S. J. Ormerod. 2003. Improving the quality of distribu-tion models for conservation by addressing shortcomings in the fieldcollection of training data. Conservation Biology 17:1601–1611.
Wahl, R., D. D. Diamond, and D. Shaw. 1990. The golden-cheekedwarbler: a status review. Texas Parks and Wildlife Department, Austin,USA.
Wiens, J. A., and J. T. Rotenberry. 1985. Response of breeding passerinebirds to rangeland alteration in a North American shrub steppe locality.Journal of Applied Ecology 22:655–668.
Zuidema, P. A., J. A. Sayer, and W. Dijkman. 1996. Forest fragmentationand biodiversity: the case for intermediate-sized conservation areas.Environmental Conservation 23:290–297.
Associate Editor: White.
8 Wildlife Society Bulletin
BIODIVERSITYRESEARCH
Predicting patch occupancy in fragmentedlandscapes at the rangewide scale for anendangered species: an example of anAmerican warbler
Bret A. Collier1*, Julie E. Groce2, Michael L. Morrison3, John C. Newnam4,
Andrew J. Campomizzi3, Shannon L. Farrell3, Heather A. Mathewson1,
Robert T. Snelgrove1, Raymond J. Carroll5 and Robert N. Wilkins1
INTRODUCTION
Species inhabiting human-dominated environments often exist
in locations where habitat loss and fragmentation have reduced
patch contiguity, patch size, and increased edge and isolation
effects (Marzluff, 2001; Bolger, 2002). Such changes in
structural features at the local scale also influence dynamics
in surrounding areas (Forman, 1995; Saab, 1999). Moreover,
1Institute of Renewable Natural Resources,
Texas A and M University, College Station,
TX 77843, USA, 2Institute of Renewable
Natural Resources, Texas A and M University,
San Antonio, TX 77843, USA, 3Department of
Wildlife and Fisheries Sciences, Texas A and
M University, College Station, TX 77843,
USA, 4Texas Department of Transportation,
PO Box 15426, Austin, TX 78761, USA,5Department of Statistics, Texas A and M
Donovan, T. M., and C. H. Flather. 2002. Relationships amongNorth American songbird trends, habitat fragmentation, and landscapeoccupancy. Ecological Applications 12:364–374.
Emlen, J. T. 1971. Population densities of birds derived from transectcounts. Auk 88:323–342.
Farrell, S. L. 2011. Use of social information for habitat selection insongbirds. Dissertation, Texas A&M University, College Station,USA.
Fitzgerald, J. A.,W. E. Thogmartin, R. Dettmers, T. Jones, C. Rustay, J. M.Ruth, F. R. Thompson, III, and T. Will. 2009. Application of models toconservation planning for terrestrial birds in North America. Pages 593–624 in J. J. Millspaugh and F. R. Thompson, III, editors. Models forplanning wildlife conservation in large landscapes. Academic Press, SanDiego, California, USA.
Gaston, K. J. 1999. Implications of interspecific and intraspecific abundance-occupancy relationships. Oikos 86:195–207.
Gaston, K. J., T. M. Blackburn, J. J. D. Greenwood, R. D. Gregory, R. M.Quinn, and J. H. Lawton. 2000. Abundance-occupancy relationships.Journal of Applied Ecology 37:39–59.
Gimenez, O., A. Gregoire, and T. Lenormand. 2009. Estimating andvisualizing fitness surfaces using mark-recapture data. Evolution63:3097–3105.
Groce, J. E., H. A. Mathewson, M. L. Morrison, and R. N. Wilkins. 2010.Scientific evaluation for the 5-year state review of the golden-cheekedwarbler. Texas A&M Institute of Renewable Natural Resources, CollegeStation, USA.
Gutzwiller, K. J., and W. C. Barrow. 2001. Bird-landscape relations in theChihuahuan desert: coping with uncertainties about predictive models.Ecological Applications 11:1517–1532.
Hatch, S. L., K. N. Gandhi, and L. E. Brown. 1990. Checklist of thevascular plants of Texas. MP-1655. Texas A&M University, TexasAgriculture Experiment Station, College Station, USA.
He, F., and K. J. Gaston. 2003. Occupancy, spatial variance, and theabundance of species. American Naturalist 162:366–375.
Hilty, J., and A. M. Merenlender. 2003. Studying biodiversity on privatelands. Conservation Biology 17:132–137.
Homer, C., J. Dewitz, J. Fry, M. Coan, N. Hossain, C. Larson, N.Herold, A. McKerrow, J. N. VanDriel, and J. Wickham. 2007.Completion of the 2001 national land cover database for the conterminousUnited States. Photogrammetric Engineering and Remote Sensing73:337–341.
Jette, L. A., T. J. Hayden, and J. D. Cornelius. 1998. Demographics of thegolden-cheeked warbler (Dendroica chrysoparia) on FortHood, Texas. U.S.ArmyCorps of Engineers, U.S. ArmyConstruction Engineering ResearchLaboratory Technical Report 98/52, Fort Hood, Texas, USA.
Johnson, D. H. 2008. In defense of indices: the case of bird surveys. Journalof Wildlife Management 72:857–868.
Kotliar, N. B., and J. A.Wiens. 1990.Multiple scales of patchiness and patchstructure: a hierarchical framework for the study of heterogeneity. Oikos59:253–260.
Laake, J. L., B. A. Collier, M. L. Morrison, and R. N.Wilkins. 2011. Point-based mark-recapture distance sampling. Journal of Agricultural,Biological, and Environmental Statistics 16:389–408.
MacKenzie, D. I., J. D. Nichols, J. A. Royle, K. H. Pollock, L. L. Bailey, andJ. E. Hines. 2006. Occupancy estimation and modeling: inferring patternsand dynamics of species occurrence. Elselvier Publishing Academic Press,Burlington, Massachusetts, USA.
Magness, D. R., R. N. Wilkins, and S. J. Hejl. 2006. Quantitative relation-ships among golden-cheeked warbler occurrence and landscape size, com-position, and structure. Wildlife Society Bulletin 34:473–479.
McCullagh, P., and J. A. Nelder. 1989. Generalized linear models. Secondedition. Chapman & Hall/CRC Press, Boca Raton, Florida, USA.
Moore, J. E., D.M. Scheiman, and R. K. Swihart. 2004. Field comparison ofremoval and modified double-observer modeling for estimating detect-ability and abundance of birds. Auk 121:865–876.
Morrison, M. L., R. N. Wilkins, B. A. Collier, J. E. Groce, H. A.Mathewson, T. M. McFarland, A. G. Snelgrove, R. T. Snelgrove, andK. Skow. 2010. Golden-cheeked warbler population distribution andabundance. Texas A&M Institute of Renewable Natural Resources,College Station, USA.
Newton, I. 1992. Experiments on the limitation of bird numbers by territo-rial behavior. Biological Reviews of the Cambridge Philosophical Society67:129–173.
Nichols, J. D., J. E. Hines, J. R. Sauer, F. W. Fallon, J. E. Fallon, and P. J.Heglund. 2000. A double-observer approach for estimating detectionprobability and abundance from point counts. Auk 117:393–408.
Peak, R. G. 2003. Population trends of the golden-cheeked warbler on FortHood, Texas 1992–2003. Chapter 5 in Endangered species monitoringand management at Fort Hood, Texas: 2003 annual report. The NatureConservancy, Fort Hood Project, Fort Hood, Texas, USA.
Peak, R. G. 2007. Forest edges negatively affect golden-cheeked warblernest survival. Condor 109:628–637.
Pollock, K. H., J. D. Nichols, T. R. Simons, G. L. Farnsworth, L. L. Bailey,and J. R. Sauer. 2002. Large scale wildlife monitoring studies: statisticalmethods for design and analysis. Environmetrics 13:105–119.
Pulich, W. M. 1976. The golden-cheeked warbler: a bioecological study.Texas Parks and Wildlife Department, Austin, Texas, USA.
Rappole, J. H., D. I. King, and J. Diez. 2003. Winter- vs. breeding-habitatlimitation for an endangered avian migrant. Ecological Applications13:735–742.
Rowell, G. A., D. P. Keddy-Hector, D. D. Diamond, J. Lloyd, B.McKinney, R. C. Maggio, and T. L. Cook. 1995. Remote sensing andGIS of golden-cheeked warbler breeding habitat and vegetation types onthe Edwards Plateau. Texas Parks and Wildlife Department, Austin,Texas, USA.
Royle, J. A., and J. D. Nichols. 2003. Estimating abundance from repeatedpresence-absence data or point counts. Ecology 84:777–790.
Scott, J. M., D. D. Goble, J. A. Wiens, D. S. Wilcove, M. Bean, and T.Male. 2005. Recovery of imperiled species under the Endangered SpeciesAct: the need for a new approach. Frontiers in Ecology and theEnvironment 3:383–389.
Shaw, D. M. 1989. Applications of GIS and remote sensing for thecharacterization of habitat for threatened and endangered species.Dissertation, University of North Texas, Denton, USA.
Sirami, C., L. Brotons, and J.-L. Martin. 2010. Erratum to: do bird spatialdistribution patterns reflect population trends in changing landscapes?Landscape Ecology 25:159–160.
Sorice, M. G., W. Haider, J. R. Conner, and R. B. Ditton. 2011. Incentivestructure of and private landowner participation in an endangered speciesconservation program. Conservation Biology 25:587–596.
Stevens, D. L., and S. F. Jensen. 2007. Sample design, execution, andanalysis for wet land assessment. Wetlands 27:515–523.
SWCA Environmental Consultants. 2007. Preliminary deliverable golden-cheeked warbler status review. Texas Department of Transportation, SanAntonio, USA.
Thompson, F. R. III, D. E. Burhans, and B. Root. 2002. Effects of pointcount protocol on bird abundance and variability estimates and power todetect population trends. Journal of Field Ornithology 73:141–150.
Thompson, S. K. 2002a. Sampling. Second edition. John Wiley and Sons,New York, New York, USA.
Thompson, W. L. 2002b. Towards reliable bird surveys: accounting forindividuals present but not detected. Auk 119:18–25.
U.S. Fish and Wildlife Service. 1990. Endangered and threatened wildlifeand plants; final rule to list the golden-cheeked warbler as endangered.Federal Register 55:53153–53160.
U.S. Fish and Wildlife Service. 1992. Golden-cheeked warbler (Dendroicachrysoparia) recovery plan. U.S. Fish and Wildlife Service, Albuquerque,New Mexico, USA.
U.S. Fish and Wildlife Service. 2004. Biological Opinion 2-12-05-F-021. U.S. Department of Agriculture, Natural Resources ConservationService, Temple, Texas, USA.
Venables,W. N., and B. D. Ripley. 2002. Modern Applied Statistics with S.Fourth edition. Springer, New York, New York, USA.
Wahl, R., D. D. Diamond, and D. Shaw. 1990. The golden-cheekedwarbler: a status review. Texas Parks and Wildlife Department, Austin,USA.
White, G. C., and R. E. Bennets. 1996. Analysis of frequency count datausing the negative binomial distribution. Ecology 77:2549–2557.