Simulations inform design of regional occupancy‐based monitoring for a sparsely distributed, territorial species

Abstract Sparsely distributed species attract conservation concern, but insufficient information on population trends challenges conservation and funding prioritization. Occupancy‐based monitoring is attractive for these species, but appropriate sampling design and inference depend on particulars of the study system. We employed spatially explicit simulations to identify minimum levels of sampling effort for a regional occupancy monitoring study design, using white‐headed woodpeckers (Picoides albolvartus), a sparsely distributed, territorial species threatened by habitat decline and degradation, as a case study. We compared the original design with commonly proposed alternatives with varying targets of inference (i.e., species range, space use, or abundance) and spatial extent of sampling. Sampling effort needed to achieve adequate power to observe a long‐term population trend (≥80% chance to observe a 2% yearly decline over 20 years) with the previously used study design consisted of annually monitoring ≥120 transects using a single‐survey approach or ≥90 transects surveyed twice per year using a repeat‐survey approach. Designs that shifted inference toward finer‐resolution trends in abundance and extended the spatial extent of sampling by shortening transects, employing a single‐survey approach to monitoring, and incorporating a panel design (33% of units surveyed per year) improved power and reduced error in estimating abundance trends. In contrast, efforts to monitor coarse‐scale trends in species range or space use with repeat surveys provided extremely limited statistical power. Synthesis and applications. Sampling resolutions that approximate home range size, spatially extensive sampling, and designs that target inference of abundance trends rather than range dynamics are probably best suited and most feasible for broad‐scale occupancy‐based monitoring of sparsely distributed territorial animal species.


| INTRODUCTION
Population monitoring informs biological conservation by revealing population trends, which inform conservation status and funding priorities (Marsh & Trenham, 2008). Conservationists focus on species experiencing severe or consistent declines due to anthropogenic impacts that elevate extinction risk (Male, Bean, & Schwartz, 2005;Rodrigues, Pilgrim, Lamoreux, Hoffmann, & Brooks, 2006). Species of uncertain status due to insufficient data are difficult to target, even if life history or declining habitat warrant concern. Information for prioritizing conservation is particularly limited for sparsely distributed species (Roberts, Taylor, & Joppa, 2016). Imperfect detectability and difficulties with modeling also impose challenges for territorial animals (Efford & Dawson, 2012;Latif, Ellis, & Amundson, 2016). Low detectability and an extensive range may necessitate broad and sustained effort to characterize population status, despite typically limited funding (Joseph, Field, Wilcox, & Possingham, 2006).
Study design for monitoring occupancy depends on desired inference and species ecology. Relatively large sampling units potentially occupied by multiple individuals can efficiently inform species range estimates, whereas smaller units may be better for tracking finer-scale changes in local abundance (Clare et al., 2015;Efford & Dawson, 2012;Noon et al., 2012). With smaller units, the timing of replicate samples used to correct for detectability in relation to territorial movement further shapes potential inference (Efford & Dawson, 2012;Latif et al., 2016;Valente, Hutchinson, & Betts, 2017). Sampling continuously distributed populations of mobile individuals with indefinite home range boundaries is especially challenging; such populations are inherently heterogeneous in ways not quantified by commonly used models, potentially obscuring inference (Efford & Dawson, 2012). More complex models that correctly specify this heterogeneity typically require more sampling effort, which may be infeasible or compromise sampling extent needed to document broad-scale trends (Welsh, Lindenmayer, & Donnelly, 2013). Simulation approaches can help inform design of occupancy-based monitoring with such inherent and unavoidable misspecification of spatial heterogeneity (Ellis, Ivan, Tucker, & Schwartz, 2015;Ellis et al., 2014).
Desired inference should primarily determine monitoring approach, but pragmatic considerations also influence study design.
Biologists may size sampling units for study area coverage or to match the resolution of available environmental data (Zielinski, Baldwin, Truex, Tucker, & Flebbe, 2013;e.g., Steenweg et al., 2016; but see Linden, Fuller, Royle, & Hare, 2017). Additionally, biologists select statistical models that best leverage available data. For example, despite a fundamental relationship between detectability and abundance (Royle & Nichols, 2003), analysts may hold detectability constant for parsimony (e.g., Zielinski et al., 2013). Sampling is often then designed to achieve adequate statistical power for tracking occupancy trends without a priori specifying desired targets of inference (e.g., species range, space use, or abundance). Inference of process, however, is ultimately needed to inform conservation.
To help fill this information gap, regional occupancy-based monitoring was established to evaluate population and distributional trends (Mellen-McLean, Saab, Bresson, Wales, & VanNorman, 2015). Repeat detection-nondetection surveys along transects in potential habitat of Oregon and Washington ( Figure 2) informed occupancy trends corrected for imperfect detection. Surveyors applied a common protocol for birds of point-based surveys oriented along transects (see also Rota, Fletcher, Dorazio, & Betts, 2009;Valente et al., 2017). Available funding was substantial (~$800 thousand) but nevertheless limited monitoring to 6 years at 30 transects while also accommodating other F I G U R E 1 Photograph of a White-headed Woodpecker objectives . Growing agency interest in white-headed woodpeckers could motivate expanded and more focused monitoring of long-term trends, which we aimed to inform.
To address the questions raised by this case study, we used simulations to evaluate alternative approaches to regional monitoring while explicitly considering potential inference and species ecology.
Spatially explicit simulations correctly represented model misspecification typically inherent with occupancy-based monitoring of continuously distributed populations, improving estimates of statistical power for informing study design (Ellis et al., 2014(Ellis et al., , 2015. We assessed minimum effort needed for desirable statistical power given the historical study design, and we explored how alternate sampling allocations influenced statistical power. We considered sampling allocations alternately suited for inferring coarse-scale distributional changes or range dynamics versus finer-scale changes in abundance or space use (e.g., Valente et al., 2017). Candidate designs varied in how they favored spatially extensive versus more intensive survey allocation, the value of which depends on population heterogeneity (Rhodes & Jonzén, 2011). Alternatives considered here represent commonly used designs for broad-scale occupancy-based studies, thus providing general guidance for monitoring sparsely distributed, territorial animals.

| White-headed Woodpecker regional monitoring
Occupancy-based monitoring of WHWO across the inland north- F I G U R E 2 National forests of the eastern Cascade Mountains, Oregon and Washington, U.S.A. White-headed Woodpecker regional monitoring focused on potential habitat (gray), where largecone pine species (mainly ponderosa) dominate Pavlacky, Blakesley, White, Hanni, & Lukacs, 2012;Rota et al., 2009), and for WHWO provided opportunity for broadcasting calls followed by a period of listening (max 5 min total) before proceeding to the next point. To conserve time, surveyors immediately proceeded to the next point along a transect when white-headed woodpecker was first detected. Thus, they strictly recorded detection-nondetection data, restricting the focus of monitoring to occupancy.

| Population and sampling simulations
Following the initial 6-year effort, we simulated occupancy-based monitoring of white-headed woodpeckers to inform potential future efforts. Recognizing the need for greater sampling effort to meaningfully quantify trends, however, we simulated surveys of ≥60 transects over 20 years.
Simulated populations experienced deterministic trends based on an exponential model, where N t is population abundance in year t and λ N is the proportion change in abundance per year (for theoretical basis, see Gotelli, 2001).
We considered a range of trend scenarios of potential conservation concern, λ N = {1.0, 0.98, 0.95, or 0.9}, that is, 0%, 33%, 64%, or 88% decline over 20 years. Simulated trends represented effect sizes for analyzing power. Positive trends (λ N > 1.0) were less of a concern for informing prioritization of WHWO for conservation action and therefore not considered. Although real populations fluctuate stochastically, we lacked information for simulating specific levels of stochasticity, and deterministic trends provided clearer effect sizes for interpreting power estimates (see also Ellis et al., 2015;MacKenzie & Royle, 2005).
We intended the range of simple deterministic trends considered here to inform surveillance monitoring aimed at documenting unanticipated change rather than particular ecological scenarios (Hutto & Belote, 2013).
We simulated population monitoring so as to explicitly represent the process of sampling discrete detection-nondetection data for continuously distributed populations (cf. Efford & Dawson, 2012;Ellis et al., 2014;contra MacKenzie & Royle, 2005). We conducted simulations using the rSPACE package (Ellis et al., 2014(Ellis et al., , 2015 in R (R Core Team, 2015). Simulations entailed (1) randomly distributing N 1 home ranges across suitable habitat, (2) calculating the probability of encountering ≥1 white-headed woodpecker at a survey point, (3) generating detection-nondetection data based on these encounter probabilities, and (4) randomly removing N t × (1 − λ N ) individuals from the landscape and repeating steps 2-3 for each remaining year t = 2-20 (Appendix S1). We collected data for a region-wide 300 m point grid, and later derived transect monitoring scenarios. Surveyors rarely recorded detections >150 m away (7% of 2011-2016 detections), so we simulated 150-m fixed-radius point surveys.
We consolidated point data to represent transect monitoring scenarios varying in sampling effort and allocation. A transect detection represented ≥1 detection at any given point along a transect on a given day. One home range (1-km radius) could include multiple neighboring points, so point-level detections within transects were spatially correlated, whereas ≥2 km transect spacing avoided spatial correlation in transect detections. Additionally, with transects, we were able to explore a fundamental issue in monitoring design: the relative merits of sampling intensively (e.g., more points per transect or repeat surveys) versus extensively (more transects).
We considered monitoring scenarios to accomplish two objectives: (1) identify levels of sampling effort capable of providing desirable power (≥80% chance to observe a decline given λ N ≤ 0.98); (2) compare commonly considered sampling allocation strategies representing alternative targets of inference and spatial extents of sampling. We addressed objective 1 by varying sampling effort (n transect ≥ 60; i.e., n pointsurveys-per-year ≥ 1,200) with different trends and the historical sampling allocation of 10 points per transect surveyed twice every year. For objective 2, we focused on a long-term decline scenario (λ N = 0.98) and fixed sampling effort (n point-surveys-per-year = 1,200) while varying monitoring strategies. The historical allocation scheme represented an intended inference of relatively coarse-scale trends. Alternative schemes included surveying shorter transects (8-3 points per transect), which targeted inference of finer-scale trends by sampling smaller areas potentially occupied by fewer individuals. We also considered surveying transects only once per year, representing single-survey occupancy approaches whose estimates provide temporal snapshots of populations, useful for inferring changes in abundance (Hutto, 2016;Latif et al., 2016). Finally, we considered surveying <100% of transects per year (i.e., panel designs; Bailey, Hines, Nichols, & MacKenzie, 2007;Urquhart & Kincaid, 1999) or repeating surveys at <100% of transects each year. Having fixed sampling effort, these alternate schemes allowed monitoring of more transects, which extended spatial sampling.
For simplicity, simulations assumed no false-negative observer error (hereafter observer error), that is, white-headed woodpeckers were always detected if present during a survey. Thus, detectability was determined exclusively by territorial movement between repeat surveys within a year. This assumption was defensible because callbroadcast surveys reduce observer error and standardized surveys limit potentially confounding interannual variation in observer error . Additionally, we calibrated simulations with pilot data (Appendix S2). Encounter probabilities during a survey therefore reflected the number, location, size, and spacing of home ranges, informed by white-headed woodpecker ecology; population size reflected calibration with pilot data and assumed trends (Appendices S1 and S2). Spatial variability in detectability (i.e., encounter probabilities in simulations) emerged from variation in nesting habitat and random placement of home ranges within this habitat, which caused local abundance and proximity to centers of activity to vary among transects. Reflecting likely realities, detectability at occupied transects increased with increasing abundance and decreased with distance from home range centers (see Appendix S1). Analyses ignored this spatial heterogeneity, and thus informed study design while accounting for likely constraints on model complexity due to limited sampling effort. We initially considered smaller home ranges (600 m radius), but calibration to pilot data required compensatory adjustments to initial abundance, resulting in similar patterns in statistical power (Q. Latif, unpublished data). We restricted simulations to national forests, representing 77% (7.7 × 10 6 ha) of potential habitat within the region ( Figure 2).

| Data analysis
For scenarios yielding repeat-survey data, we estimated trends with two different occupancy models representing commonly considered ways of correcting for detectability (p; e.g., Linden et al., 2017;Steenweg et al., 2016;Zielinski et al., 2013) to estimate occupancy probability (ψ; Figure 3a,b). One model allowed detectability estimates to vary interannually (hereafter the yearly-p model; Figure 3a), whereas the other held detectability constant (hereafter the constantp model; Figure 3b; for model structures, see Appendix S3). Because individuals could move in or out of the surveyed area between repeat surveys within a year, these models quantified the probability of a transect intersecting ≥1 home range, hereafter true occupancy, which describes species range or space use (Efford & Dawson, 2012;MacKenzie & Royle, 2005). The yearly-p model allows detectability to change with changing abundance (Royle & Nichols, 2003) to better estimate true occupancy. The constant-p model misspecifies true occupancy, but is frequently considered and may be selected for parsimony in applied studies (e.g., Zielinski et al., 2013). Additionally, having controlled for observer error (e.g., if nonexistent as in simulations, or controlled via standardized surveys), the constant-p model coerces occupancy estimates to reflect any interannual changes, shifting the target of inference to abundance (Figure 3b).
For scenarios yielding single-survey data, we estimated trends using logistic regression (see structure in Appendix S3). Having excluded observer error in simulations, logistic regression models estimated probability of ≥1 individual's physical presence during a survey, hereafter probability of physical presence. Single-survey scenarios represented single-survey occupancy approaches, in which replicate surveys occur within a narrow enough timeframe for detectability to quantify observer error so that occupancy estimates quantify probability of physical presence (e.g., double-observer and removal designs; Nichols et al., 2008;Rota et al., 2009). By omitting observer error from simulations, however, replicate surveys were unnecessary to quantify physical presence. Probability of physical presence represents a temporal snapshot of a population unaffected by territorial movement F I G U R E 3 How model estimates reflect underlying processes under alternative monitoring approaches. Repeat-survey occupancy estimates fundamentally quantify true occupancy (a, b) but can index abundance trends if detectability is held constant (b, contra A). Territorial movement influences repeat-survey occupancy estimates (a, b), whereas single-survey occupancy estimates represent population snapshots not influenced by movement (c). False-negative observer error was not simulated but could influence estimates in reality expected to closely track abundance (Figure 3c; Hutto, 2016;Latif et al., 2016). Additionally, surveying transects only once allowed us to monitor twice as many transects, increasing sampling extent. In practice, auxiliary sampling (e.g., recording detection timing or deploying multiple observers) would account for observer error (Nichols et al., 2008;Rota et al., 2009), likely adding to uncertainty in trend estimates. Ignoring observer error in both single-and repeat-survey scenarios, however, made their comparison informative.
We quantified occupancy trends as proportion yearly change in odds occupancy, λ ψ = et al., 2006). We analyzed detection-nondetection data with fixed year effects and subsequently calculated least-squares trends in yearly occupancy estimates (see also Ellis et al., 2014). We quantified statistical power as percent simulations when 95% Bayesian credible intervals (BCIs) for the estimated trend for the study period (λ ψ where logit(ψ t ) = β 0 + log (λ ψ ) × t;

| Occupancy, abundance, and estimator behavior
Comparing true occupancy (proportion transects with encounter p ≥ .05) and abundance informed understanding of statistical power and estimator properties. True occupancy related positively with abundance but plateaued at higher abundances (Figure 4a,c). True occupancy declines lagged abundance declines (Figure 4b,d). With shorter transects, true occupancy corresponded better but still imperfectly with abundance (Figure 4c,d).
Occupancy estimates remained constant with no abundance trend and declined with declining abundance ( Figure 5). Detectability estimates declined with declining abundance either across scenarios (constant-p estimates) or through time (yearly-p estimates; Figure 6).
Yearly-p estimator precision was less than for constant-p estimates and declined with declining abundance (Figures 5 and 6).
Detectability estimates followed the behavior of encounter probabilities at occupied transects but were generally higher, that is, positively biased (Figure 6), making occupancy estimates negatively biased ( Figure 5). Logistic regression estimates deviated even more from true occupancy (Figure 5c,f,i,l), reflecting the differing target of inference (i.e., probability of physical presence; see Section 2 and Appendix S3).

| Objective 1: Sampling effort
With historical survey allocation, statistical power increased with increasing sampling effort and stronger population declines (Figure 7). With no abundance or occupancy trends, models estimated actual trends with minimal error and no apparent bias, but error and bias grew with increasing trend (Figure 8)

| Objective 2: Sampling allocation
Monitoring strategies that targeted inference of finer-scale trends in space use or abundance and extended sampling spatially generally F I G U R E 5 Yearly occupancy estimates from simulated regional white-headed woodpecker monitoring. Simulated trends were λ N = 1 (a-c), 0.98 (d-f), 0.95 (g-i), and 0.9 (j-l). Repeat-survey occupancy estimates assumed constant detectability (a, d, g, and j) or variable detectability among years (b, e, h, and k). Single-survey estimates assumed perfect detectability (c, f, i, and l). Thirty simulations of monitoring transects of 10 points each for 20 years are represented for each scenario (n = 150 and 300 transects for repeat-and single-survey scenarios, respectively). Black dots and blue vertical bars show yearly estimates and 95% BCIs jittered for display. Black lines connect estimates from consecutive years for individual simulations. Red dots show mean true occupancy for 30 simulations, that is, proportion of all possible transects occupied provided more power and less estimation error than the historical strategy. Power improved and estimation error decreased when monitoring shorter but more transects (Figures 9 and 10). Power was greatest and error (RMSE N ) least when monitoring the probability of physical presence with single-survey data (Figures 8c and 10).
In contrast, we found the least power and greatest error (RMSE ψ ) when attempting to monitor true occupancy with repeat surveys and the yearly-p model (Figure 8b). Interestingly, despite more explicitly targeting inference of occupancy, yearly-p trends estimates did not estimate occupancy trends with any less error (RMSE ψ ≤ 0.105) than logistic regression (RMSE ψ ≤ 0.055). Although it provided adequate power in many scenarios, the constant-p model tended to overestimate both occupancy and abundance declines (Figures 8a   and 9).
Panel designs with relatively small panels (33% of transects surveyed each year) also improved power and reduced error, although larger panels (50% of transects surveyed each year) did not provide notable gains. Conducting fewer repeat surveys in exchange for more transects also did not substantively affect power and tended to increase estimation error (Figure 11).
The design that maximized power and minimized error with constant-p and logistic regression models was a 33% panel design with 3 points per transects. Even with this design, the yearly-p model provided inadequate power (13%), although trend estimation error was less than with the historical design (RMSE ψ = 0.009; n sim = 100; n transect = 600 over 20 years; compare with Figure 7b).

| Sampling resolution and scale of inference
Our results further emphasize the benefits of sampling at resolutions (i.e., unit size, grid cell size) approximating the size of an individual home range documented by others (Efford & Dawson, 2012;Linden et al., 2017). Finer resolution sampling generates occupancy estimates that more closely track abundance. This estimator property should be desirable for practitioners who monitor occupancy in lieu of abundance primarily on pragmatic grounds.
Single-survey sampling can similarly benefit monitoring of territorial animals by providing temporal snapshots of populations unaffected by movement and therefore closely related to abundance (Hutto, 2016;Latif et al., 2016). We simulated an ideal world with no observer error wherein single-survey estimates were readily interpre-  Nichols et al., 2008). Lacking data on observer error, naïve occupancy could usefully index abundance if we are confident that observer error does not vary interannually and therefore cannot confound trend estimation (e.g., with standardizing bird surveys; Hutto, 2016).
Observer error can vary with local abundance (Royle & Nichols, 2003), potentially introducing noise not represented in our simulations. Larger sampling units potentially occupied by multiple individuals would be most prone to such variability, so aligning sampling resolution with home range size would be desirable even with a singlesurvey design.
A single-survey design would require consistently conducting surveys when individuals are readily detectable. With sensitivity of nest survival to temperature (Hollenbeck et al., 2011), climate change may alter nesting phenology, potentially influencing responsiveness to call broadcasts. Such changes could necessitate adjusting the timing of surveys, which could be informed by targeted repeat surveys, as are commonly implemented for birds (Latif, Fleming, Barrows, & Rotenberry, 2012;Rota et al., 2009).
Our results indicate challenges for monitoring to infer changes in species range (coarse-scale) or space use (finer scale) as in yearly-p scenarios here. By definition, occupancy only declines when abundance declines enough to result in local extirpation, so true occupancy declines could indicate strong need for conservation. More F I G U R E 7 Simulation-based power to observe white-headed woodpecker regional occupancy trends (percent simulations with 95% BCI <1). Scenarios varied in number of transects, monitoring approach (constant-p or yearly-p occupancy models, or logistic regression), and trend (λ N = exponential change in abundance). For all scenarios, transects consisted of 10 survey points surveyed twice (occupancy models) or once (logistic regression) per year. Constant-p assumed constant detectability, whereas yearly-p allowed variable detectability among years. Logistic regression allowed double the number of transects by analyzing singlesurvey data intensive sampling in areas or years with low abundance, however, may be needed to correctly identify occupancy declines. At finer scales, spatial heterogeneity in detectability arising from variability in local abundance and home ranges that lack definitive boundaries can limit accurate estimation of space use (Efford & Dawson, 2012).
Biases in occupancy and detectability estimates observed here likely primarily reflect these effects. Including habitat relationships with occupancy in analytical models (omitted from simulations) might help by accounting somewhat for spatial heterogeneity in the data, but effects on detectability of varying local abundance and proximity to home ranges at occupied transects would remain. Effectively estimating species distribution at any scale may require substantial spatial or temporal replication within sampling units (Pavlacky et al., 2012;Valente et al., 2017). Given likely demands on funding, such approaches may be feasibly implemented only infrequently (e.g., Cruickshank, Ozgul, Zumbach, & Schmidt, 2016). Alternatively, predictive models (e.g., Hollenbeck et al., 2011;Latif et al., 2015;Wightman, Saab, Forristal, Mellen-McLean, & Markus, 2010) could supplement trend monitoring by identifying changes in habitat.
Nested surveys (e.g., points along transects) can inform hierarchically structured models capable of estimating patterns or trends at multiple scales (Pavlacky et al., 2012;Rota et al., 2009;Royle & Kéry, 2007). Multiscale inference would require sufficient sampling at all scales of interest, however, which may be beyond resources available for many monitoring programs (Valente et al., 2017). Our initial attempts found inadequate sampling for meaningful multiscale inference (Q. Latif, unpublished data), so we abandoned such approaches here.
In our simulations, spatial heterogeneity emerged from uneven distribution of habitat and random variation in local abundance among occupied transects. The benefits observed here with shorter transects and single surveys could reflect advantages of spatially extended sampling. Panel designs, however, did not inherently change the target of inference, and so their results more definitively demonstrated potential advantages with spatially extended sampling. In contrast, ignoring heterogeneity inherent in continuously distributed territorial species may obscure advantages of panel designs (Bailey et al., 2007;Urquhart & Kincaid, 1999).
Not all spatial extensions to sampling were beneficial. Given a repeat-survey design, we gained nothing by reducing repeat surveys to monitor more transects. Such strategies require high detectability (MacKenzie & Royle, 2005) likely uncharacteristic of sparsely and continuously distributed territorial species (see above). The lack F I G U R E 8 Correspondence of estimated occupancy trends (λ ψ ) with true abundance (λ N ) and occupancy (λ ψ ) trends. Trends were estimated with repeat-survey constant-p (a) and yearly-p (b) occupancy models, and logistic regression analyzing single-survey data (c). Red and blue dots mark perfect correspondence with actual abundance and occupancy trends, respectively. Root mean squared error quantifies the overall estimation error with respect to abundance (RMSE N ) and occupancy (RMSE ψ ) trends of benefit with 50% panels may reflect site fidelity, fixed at 100% in our simulations. By monitoring different transects in successive years, the number and distribution of individuals along surveyed transects varied interannually, potentially obscuring trends. Whiteheaded woodpecker do exhibit site fidelity (Garrett et al., 1996), so panel design benefits could trade-off with benefits of sampling the same sets of individuals in successive years. In reality, however, population processes (e.g., dispersal and turnover) could also obscure trends. The extent of paneling needed to benefit power would therefore depend on levels of spatial versus temporal heterogeneity in population trends (Rhodes & Jonzén, 2011), the latter of which was omitted from simulations here. Additionally, panel designs could limit study of processes underlying occupancy dynamics, for example, colonization and persistence (Bailey et al., 2007).

| Study limitations
Our simulations did not include spatial or temporal stochasticity in population dynamics, individual movement between years, or behavioral interactions between neighbors all of which could modulate occupancy estimates or trends (Reynolds, Wiens, Joy, & Salafsky, 2005;Sauer, Fallon, & Johnson, 2003;Warren, Veech, Weckerly, O'Donnell, & Ott, 2013). By not accounting for these realities, power estimates may be liberal, and therefore probably best used to inform a lower bound for sample size. Accordingly, we recommend ≥120 or ≥90 transects with single-survey or repeat-survey monitoring, respectively, of white-headed woodpeckers across our study region.
Our treatment of site fidelity, however, is conservative. For simplicity, we simulated populations with 100% site fidelity and zero immigration or recruitment, and to avoid artifacts of these assumptions, analysis models assumed occupancy varied independently among years. In reality, models correctly specifying uncertainty arising from additional population processes (e.g., Royle & Kéry, 2007) could improve power to observe trends (although with likely increased data demands). Models that correctly specify habitat relationships with occupancy could also help. Given the potentially counteracting features of simulations, we expect power estimates were sufficiently informative to compare alternative study designs.
Simulations ignored spatial variation in home range size, which can confound interpretation of occupancy estimates and trends drawn from repeat surveys (Efford & Dawson, 2012

| Additional considerations and broader implications
Agency biologists often conduct repeat surveys to estimate detectability and thereby improve credibility of trend estimates. This strategy potentially implies an overly rigid allocation of effort between sampling to inform occupancy versus detectability. Repeat surveys of mobile species may unwittingly focus effort toward tracking distributional shifts, which can be harder to observe and not necessarily more relevant to conservation than changes in abundance. Additionally, practitioners often discount the potential for detectability to change with changing abundance (e.g., Ahumada, Hurtado, & Lizcano, 2013;van Strien, van Swaay, & Termaat, 2013;Zielinski et al., 2013), which may limit explicit inference of abundance trends versus range dynamics from occupancy-based trend estimates. Estimating detectability is only useful if doing so improves inference of underlying population processes or accounts for interannual variability in observer error.
The former requires considering which processes can be more readily inferred by accounting for detectability at the scale it is measured.
If instead biologists are solely concerned with controlling observer error, monitoring of population indices may be more cost-effective while providing equivalent or stronger inference of population change (Hutto, 2016;Johnson, 2008;Welsh et al., 2013).
Despite growing sophistication of occupancy models (Bailey, MacKenzie, & Nichols, 2014), heterogeneity arising from locally varying abundance and poorly defined home range boundaries (Efford & Dawson, 2012) will continue to challenge monitoring efforts, especially where funding constrains data and, consequently, model complexity. Simulations can help explore our capacity for inference with models necessarily misspecified due to limited data. General power formulas available for occupancy models ignore spatial heterogeneity Spatially explicit simulations therefore complement these tools for tailoring sampling designs to particular study systems. Information on regional trends should be combined with information on various population parameters measured at different scales to fully inform species conservation status (Nichols & Williams, 2006). For example, other studies currently underway examine forest management effects on white-headed woodpecker nest densities, nest survival, and habitat use . Statistical models can now integrate multiple sources of data to better inform parameter estimation (Dorazio, 2014;Nichols et al., 2008). Simulations that explicitly and distinctly describe population from observation processes could inform sampling design to support these approaches.

ACKNOWLEDGMENTS
Region 6 of the U.S. Forest Service primarily funded the initial regional monitoring effort, including this study. Rocky Mountain Research Station and Montana State University also supported this work. J.
Dudley helped with some data management. We thank field crews, specifically crew leaders: E. Johnson, S. Mellmann-Brown, and A.
Woodrow. We thank B. Bird, J. Tucker, C. Amundson, J. Hines, and two anonymous reviewers for feedback on initial drafts. T. Kogut provided the Figure 1 photo.

CONFLICT OF INTEREST
None declared.

AUTHOR CONTRIBUTIONS
Q. Latif implemented simulations and data analysis, and drafted the manuscript. M. Ellis, V. Saab, and K. Mellen-McLean provided substantive feedback during conceptual framing, simulation design and implementation, and manuscript preparation.

DATA ACCESSIBILITY
R scripts for initiating simulations using the rSPACE package and BUGS code defining data analysis models are provided online supporting information.