Quantifying intraspecific variation in host resistance and tolerance to a lethal pathogen.

Testing for intraspecific variation for host tolerance or resistance in wild populations is important for informing conservation decisions about captive breeding, translocation, and disease treatment. Here, we test the importance of tolerance and resistance in multiple populations of boreal toads (Anaxyrus boreas boreas) against Batrachochytrium dendrobatidis (Bd), the amphibian fungal pathogen responsible for the greatest host biodiversity loss due to disease. Boreal toads have severely declined in Colorado (CO) due to Bd, but toad populations challenged with Bd in western Wyoming (WY) appear to be less affected. We used a common garden infection experiment to expose post-metamorphic toads sourced from four populations (2 in CO and 2 in WY) to Bd and monitored changes in mass, pathogen burden and survival for 8 weeks. We used a multi-state modelling approach to estimate weekly survival and transition probabilities between infected and cleared states, reflecting a dynamic infection process that traditional approaches fail to capture. We found that WY boreal toads are more tolerant to Bd infection with higher survival probabilities than those in CO when infected with identical pathogen burdens. WY toads also appeared more resistant to Bd with a higher probability of infection clearance and an average of 5 days longer to reach peak infection burdens. Our results demonstrate strong intraspecific differences in tolerance and resistance that likely contribute to why population declines vary regionally across this species. Our multi-state framework allowed us to gain inference on typically hidden disease processes when testing for host tolerance or resistance. Our findings demonstrate that describing an entire host species as 'tolerant' or 'resistant' (or lack thereof) is unwise without testing for intraspecific variation.

Addressing Detection Uncertainty in Bombus affinis (Hymenoptera: Apidae) Surveys Can Improve Inferences Made From Monitoring.

The U.S. Fish and Wildlife Service developed national guidelines to track species recovery of the endangered rusty patched bumble bee [Bombus affinis Cresson (Hymenoptera: Apidae)] and to investigate changes in species occupancy across space and time. As with other native bee monitoring efforts, managers have specifically acknowledged the need to address species detection uncertainty and determine the sampling effort required to infer species absence within sites. We used single-season, single-species occupancy models fit to field data collected in four states to estimate imperfect detection of B. affinis and to determine the survey effort required to achieve high confidence of species detection. Our analysis revealed a precipitous, seasonal, decline in B. affinis detection probability throughout the July through September sampling window in 2021. We estimated that six, 30-min surveys conducted in early July are required to achieve a 95% cumulative detection probability, whereas >10 surveys would be required in early August to achieve the same level of confidence. Our analysis also showed B. affinis was less likely to be detected during hot and humid days and at patches of reduced habitat quality. Bombus affinis was frequently observed on Monarda fistulosa (Lamiales: Lamiaceae), followed by [Pycnanthemum virginianum Rob. and Fernald (Lamiales: Lamiaceae)], Eutrochium maculatum Lamont (Asterales: Asteraceae), and Veronicastrum virginicum Farw. (Lamiales: Plantaginaceae). Although our research is focused on B. affinis, it is relevant for monitoring other bumble bees of conservation concern, such as B. occidentalis Greene (Hymenoptera: Apidae) and B. terricola Kirby (Hymenoptera: Apidae) for which monitoring efforts have been recently initiated and occupancy is a variable of conservation interest.

Range-wide sources of variation in reproductive rates of northern spotted owls.

We conducted a range-wide investigation of the dynamics of site-level reproductive rate of northern spotted owls using survey data from 11 study areas across the subspecies geographic range collected during 1993-2018. Our analytical approach accounted for imperfect detection of owl pairs and misclassification of successful reproduction (i.e., at least one young fledged) and contributed further insights into northern spotted owl population ecology and dynamics. Both nondetection and state misclassification were important, especially because factors affecting these sources of error also affected focal ecological parameters. Annual probabilities of site occupancy were greatest at sites with successful reproduction in the previous year and lowest for sites not occupied by a pair in the previous year. Site-specific occupancy transition probabilities declined over time and were negatively affected by barred owl presence. Overall, the site-specific probability of successful reproduction showed substantial year-to-year fluctuations and was similar for occupied sites that did or did not experience successful reproduction the previous year. Site-specific probabilities for successful reproduction were very small for sites that were unoccupied the previous year. Barred owl presence negatively affected the probability of successful reproduction by northern spotted owls in Washington and California, as predicted, but the effect in Oregon was mixed. The proportions of sites occupied by northern spotted owl pairs showed steep, near-monotonic declines over the study period, with all study areas showing the lowest observed levels of occupancy to date. If trends continue it is likely that northern spotted owls will become extirpated throughout large portions of their range in the coming decades.

Individual heterogeneity influences the effects of translocation on urban dispersal of an invasive reptile.

BACKGROUND: Invasive reptiles pose a serious threat to global biodiversity, but early detection of individuals in an incipient population is often hindered by their cryptic nature, sporadic movements, and variation among individuals. Little is known about the mechanisms that affect the movement of these species, which limits our understanding of their dispersal. Our aim was to determine whether translocation or small-scale landscape features affect movement patterns of brown treesnakes (Boiga irregularis), a destructive invasive predator on the island of Guam. METHODS: We conducted a field experiment to compare the movements of resident (control) snakes to those of snakes translocated from forests and urban areas into new urban habitats. We developed a Bayesian hierarchical model to analyze snake movement mechanisms and account for attributes unique to invasive reptiles by incorporating multiple behavioral states and individual heterogeneity in movement parameters. RESULTS: We did not observe strong differences in mechanistic movement parameters (turning angle or step length) among experimental treatment groups. We found some evidence that translocated snakes from both forests and urban areas made longer movements than resident snakes, but variation among individuals within treatment groups weakened this effect. Snakes translocated from forests moved more frequently from pavement than those translocated from urban areas. Snakes translocated from urban areas moved less frequently from buildings than resident snakes. Resident snakes had high individual heterogeneity in movement probability. CONCLUSIONS: Our approach to modeling movement improved our understanding of invasive reptile dispersal by allowing us to examine the mechanisms that influence their movement. We also demonstrated the importance of accounting for individual heterogeneity in population-level analyses, especially when management goals involve eradication of an invasive species.

Patch utilization and flower visitations by wild bees in a honey bee-dominated, grassland landscape.

Understanding habitat needs and patch utilization of wild and managed bees has been identified as a national research priority in the United States. We used occupancy models to investigate patterns of bee use across 1030 transects spanning a gradient of floral resource abundance and richness and distance from apiaries in the Prairie Pothole Region (PPR) of the United States. Estimates of transect use by honey bees were nearly 1.0 during our 3.5-month sampling period, suggesting honey bees were nearly ubiquitous across transects. Wild bees more frequently used transects with higher flower richness and more abundant flowers; however, the effect size of the native flower abundance covariate ( ß ^ native  = 3.90 ± 0.65 [1SE]) was four times greater than the non-native flower covariate ( ß ^ n o n - n a t i v e  = 0.99 ± 0.17). We found some evidence that wild bee use was lower at transects near commercial apiaries, but the effect size was imprecise ( ß ^ distance  = 1.4 ± 0.81). Honey bees were more frequently detected during sampling events with more non-native flowers and higher species richness but showed an uncertain relationship with native flower abundance. Of the 4039 honey bee and flower interactions, 85% occurred on non-native flowers, while only 43% of the 738 wild bee observations occurred on non-native flowers. Our study suggests wild bees and honey bees routinely use the same resource patches in the PPR but often visit different flowering plants. The greatest potential for resource overlap between honey bees and wild bees appears to be for non-native flowers in the PPR. Our results are valuable to natural resource managers tasked with supporting habitat for managed and wild pollinators in agroecosystems.

Using a multistate occupancy approach to determine molecular diagnostic accuracy and factors affecting avian haemosporidian infections.

The use of a sensitive and accurate parasite detection methodology is crucial in studies exploring prevalence of parasites in host populations or communities, and uncertainty in identifying parasite genera and/or lineages may limit the understanding of host-parasite interactions. Here, we used a multistate occupancy approach that accounts for imperfect detection to assess whether sex and breeding season influenced the prevalence of a specific Haemoproteus lineage (TARUF02) in a white-lined tanager population. Likewise, we explored whether the probability of detecting the target parasite in an infected bird using PCR and sequencing analyses may be influenced by season and host sex. We found little evidence that sex influenced the probability of an individual host being infected by a haemosporidian parasite. Conversely, we found that the probability of infection by Haemoproteus TARUF02 was ~30% higher during the breeding season, reflecting a higher prevalence of this parasite in this season. The probability that PCR detects DNA of haemosporidian parasite was higher for female birds, suggesting that they are more prone to be parasitized with parasitemia levels that are more successfully detected by molecular analysis. Sequencing successfully determined the Haemoproteus TARUF02 lineage in 60% of samples collected during the breeding season and 84% of samples collected during the non-breeding season. Understanding the ecology of hosts and aspects of their physiology that may influence the parasite infection is essential to better understanding of hemoparasite infections and how parasites influence their native hosts, through decreasing reproductive success, lifespan, and/or survival.

The past and future roles of competition and habitat in the range-wide occupancy dynamics of Northern Spotted Owls.

Slow ecological processes challenge conservation. Short-term variability can obscure the importance of slower processes that may ultimately determine the state of a system. Furthermore, management actions with slow responses can be hard to justify. One response to slow processes is to explicitly concentrate analysis on state dynamics. Here, we focus on identifying drivers of Northern Spotted Owl (Strix occidentalis caurina) territorial occupancy dynamics across 11 study areas spanning their geographic range and forecasting response to potential management actions. Competition with Barred Owls (Strix varia) has increased Spotted Owl territory extinction probabilities across all study areas and driven recent declines in Spotted Owl populations. Without management intervention, the Northern Spotted Owl subspecies will be extirpated from parts of its current range within decades. In the short term, Barred Owl removal can be effective. Over longer time spans, however, maintaining or improving habitat conditions can help promote the persistence of northern spotted owl populations. In most study areas, habitat effects on expected Northern Spotted Owl territorial occupancy are actually greater than the effects of competition from Barred Owls. This study suggests how intensive management actions (removal of a competitor) with rapid results can complement a slower management action (i.e., promoting forest succession).

Multistate occupancy modeling improves understanding of amphibian breeding dynamics in the Greater Yellowstone Area.

Discerning the determinants of species occurrence across landscapes is fundamental to their conservation and management. In spatially and climatologically complex landscapes, explaining the dynamics of occurrence can lead to improved understanding of short- vs. long-term trends and offer novel insight on local vs. regional change. We examined the changes in occupancy for two species of anurans with different life histories over a decade using hundreds of wetland sites in Yellowstone and Grand Teton National Parks. To account for the joint dynamics of wetland drying and amphibian breeding, we adopted a multistate occupancy model as a means to investigate mechanistic relationships of observed occurrence patterns with climatological drivers of wetland hydrologic variability. This approach allowed us to decompose occupancy dynamics into habitat changes caused by wetland drying and amphibian breeding activity, conditional on available water and previous breeding state. Over our 10-yr time series, we observed considerable variability in climate drivers and the proportion of dry wetlands. Boreal chorus frogs (Pseudacris maculata) were more responsive to changes in wetland inundation status than Columbia spotted frogs (Rana luteiventris), as indicated by higher breeding colonization probabilities under favorable (wet) conditions. Both species had high probabilities of breeding persistence in permanently inundated wetlands with prior breeding. Despite the absence of multi-year drought in our time series, mechanistic relationships described here offer insights on how future climate variation may result in reduced and/or shifted occurrence patterns for pond-breeding anurans in the Greater Yellowstone Area. Further, our modeling approach may prove valuable in evaluating determinants of occurrence for other species that are dependent on wetlands or other dynamic habitats.

Quantifying climate sensitivity and climate-driven change in North American amphibian communities.

Changing climate will impact species' ranges only when environmental variability directly impacts the demography of local populations. However, measurement of demographic responses to climate change has largely been limited to single species and locations. Here we show that amphibian communities are responsive to climatic variability, using >500,000 time-series observations for 81 species across 86 North American study areas. The effect of climate on local colonization and persistence probabilities varies among eco-regions and depends on local climate, species life-histories, and taxonomic classification. We found that local species richness is most sensitive to changes in water availability during breeding and changes in winter conditions. Based on the relationships we measure, recent changes in climate cannot explain why local species richness of North American amphibians has rapidly declined. However, changing climate does explain why some populations are declining faster than others. Our results provide important insights into how amphibians respond to climate and a general framework for measuring climate impacts on species richness.

Beyond the swab: ecosystem sampling to understand the persistence of an amphibian pathogen.

Understanding the ecosystem-level persistence of pathogens is essential for predicting and measuring host-pathogen dynamics. However, this process is often masked, in part due to a reliance on host-based pathogen detection methods. The amphibian pathogens Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) are pathogens of global conservation concern. Despite having free-living life stages, little is known about the distribution and persistence of these pathogens outside of their amphibian hosts. We combine historic amphibian monitoring data with contemporary host- and environment-based pathogen detection data to obtain estimates of Bd occurrence independent of amphibian host distributions. We also evaluate differences in filter- and swab-based detection probability and assess inferential differences arising from using different decision criteria used to classify samples as positive or negative. Water filtration-based detection probabilities were lower than those from swabs but were > 10%, and swab-based detection probabilities varied seasonally, declining in the early fall. The decision criterion used to classify samples as positive or negative was important; using a more liberal criterion yielded higher estimates of Bd occurrence than when a conservative criterion was used. Different covariates were important when using the liberal or conservative criterion in modeling Bd detection. We found evidence of long-term Bd persistence for several years after an amphibian host species of conservation concern, the boreal toad (Anaxyrus boreas boreas), was last detected. Our work provides evidence of long-term Bd persistence in the ecosystem, and underscores the importance of environmental samples for understanding and mitigating disease-related threats to amphibian biodiversity.

Host-pathogen metapopulation dynamics suggest high elevation refugia for boreal toads.

Emerging infectious diseases are an increasingly common threat to wildlife. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), is an emerging infectious disease that has been linked to amphibian declines around the world. Few studies exist that explore amphibian-Bd dynamics at the landscape scale, limiting our ability to identify which factors are associated with variation in population susceptibility and to develop effective in situ disease management. Declines of boreal toads (Anaxyrus boreas boreas) in the southern Rocky Mountains are largely attributed to chytridiomycosis but variation exists in local extinction of boreal toads across this metapopulation. Using a large-scale historic data set, we explored several potential factors influencing disease dynamics in the boreal toad-Bd system: geographic isolation of populations, amphibian community richness, elevational differences, and habitat permanence. We found evidence that boreal toad extinction risk was lowest at high elevations where temperatures may be suboptimal for Bd growth and where small boreal toad populations may be below the threshold needed for efficient pathogen transmission. In addition, boreal toads were more likely to recolonize high elevation sites after local extinction, again suggesting that high elevations may provide refuge from disease for boreal toads. We illustrate a modeling framework that will be useful to natural resource managers striving to make decisions in amphibian-Bd systems. Our data suggest that in the southern Rocky Mountains high elevation sites should be prioritized for conservation initiatives like reintroductions.

ATLANTIC-CAMTRAPS: a dataset of medium and large terrestrial mammal communities in the Atlantic Forest of South America.

Our understanding of mammal ecology has always been hindered by the difficulties of observing species in closed tropical forests. Camera trapping has become a major advance for monitoring terrestrial mammals in biodiversity rich ecosystems. Here we compiled one of the largest datasets of inventories of terrestrial mammal communities for the Neotropical region based on camera trapping studies. The dataset comprises 170 surveys of medium to large terrestrial mammals using camera traps conducted in 144 areas by 74 studies, covering six vegetation types of tropical and subtropical Atlantic Forest of South America (Brazil and Argentina), and present data on species composition and richness. The complete dataset comprises 53,438 independent records of 83 species of mammals, includes 10 species of marsupials, 15 rodents, 20 carnivores, eight ungulates and six armadillos. Species richness averaged 13 species (±6.07 SD) per site. Only six species occurred in more than 50% of the sites: the domestic dog Canis familiaris, crab-eating fox Cerdocyon thous, tayra Eira barbara, south American coati Nasua nasua, crab-eating raccoon Procyon cancrivorus and the nine-banded armadillo Dasypus novemcinctus. The information contained in this dataset can be used to understand macroecological patterns of biodiversity, community, and population structure, but also to evaluate the ecological consequences of fragmentation, defaunation, and trophic interactions.

A Model to Inform Management Actions as a Response to Chytridiomycosis-Associated Decline.

Decision-analytic models provide forecasts of how systems of interest will respond to management. These models can be parameterized using empirical data, but sometimes require information elicited from experts. When evaluating the effects of disease in species translocation programs, expert judgment is likely to play a role because complete empirical information will rarely be available. We illustrate development of a decision-analytic model built to inform decision-making regarding translocations and other management actions for the boreal toad (Anaxyrus boreas boreas), a species with declines linked to chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd). Using the model, we explored the management implications of major uncertainties in this system, including whether there is a genetic basis for resistance to pathogenic infection by Bd, how translocation can best be implemented, and the effectiveness of efforts to reduce the spread of Bd. Our modeling exercise suggested that while selection for resistance to pathogenic infection by Bd could increase numbers of sites occupied by toads, and translocations could increase the rate of toad recovery, efforts to reduce the spread of Bd may have little effect. We emphasize the need to continue developing and parameterizing models necessary to assess management actions for combating chytridiomycosis-associated declines.

Design- and model-based recommendations for detecting and quantifying an amphibian pathogen in environmental samples.

Accurate pathogen detection is essential for developing management strategies to address emerging infectious diseases, an increasingly prominent threat to wildlife. Sampling for free-living pathogens outside of their hosts has benefits for inference and study efficiency, but is still uncommon. We used a laboratory experiment to evaluate the influences of pathogen concentration, water type, and qPCR inhibitors on the detection and quantification of Batrachochytrium dendrobatidis (Bd) using water filtration. We compared results pre- and post-inhibitor removal, and assessed inferential differences when single versus multiple samples were collected across space or time. We found that qPCR inhibition influenced both Bd detection and quantification in natural water samples, resulting in biased inferences about Bd occurrence and abundance. Biases in occurrence could be mitigated by collecting multiple samples in space or time, but biases in Bd quantification were persistent. Differences in Bd concentration resulted in variation in detection probability, indicating that occupancy modeling could be used to explore factors influencing heterogeneity in Bd abundance among samples, sites, or over time. Our work will influence the design of studies involving amphibian disease dynamics and studies utilizing environmental DNA (eDNA) to understand species distributions.

Gene flow from an adaptively divergent source causes rescue through genetic and demographic factors in two wild populations of Trinidadian guppies.

Genetic rescue, an increase in population growth owing to the infusion of new alleles, can aid the persistence of small populations. Its use as a management tool is limited by a lack of empirical data geared toward predicting effects of gene flow on local adaptation and demography. Experimental translocations provide an ideal opportunity to monitor the demographic consequences of gene flow. In this study we take advantage of two experimental introductions of Trinidadian guppies to test the effects of gene flow on downstream native populations. We individually marked guppies from the native populations to monitor population dynamics for 3 months before and 26 months after gene flow. We genotyped all individuals caught during the first 17 months at microsatellite loci to classify individuals by their genetic ancestry: native, immigrant, F1 hybrid, F2 hybrid, or backcross. Our study documents a combination of demographic and genetic rescue over multiple generations under fully natural conditions. Within both recipient populations, we found substantial and long-term increases in population size that could be attributed to high survival and recruitment caused by immigration and gene flow from the introduction sites. Our results suggest that low levels of gene flow, even from a divergent ecotype, can provide a substantial demographic boost to small populations, which may allow them to withstand environmental stochasticity.

Multi-scale occupancy approach to estimate Toxoplasma gondii prevalence and detection probability in tissues: an application and guide for field sampling.

Increasingly, birds are recognised as important hosts for the ubiquitous parasite Toxoplasma gondii, although little experimental evidence exists to determine which tissues should be tested to maximise the detection probability of T. gondii. Also, Arctic-nesting geese are suspected to be important sources of T. gondii in terrestrial Arctic ecosystems, but the parasite has not previously been reported in the tissues of these geese. Using a domestic goose model, we applied a multi-scale occupancy framework to demonstrate that the probability of detection of T. gondii was highest in the brain (0.689, 95% confidence interval=0.486, 0.839) and the heart (0.809, 95% confidence interval=0.693, 0.888). Inoculated geese had an estimated T. gondii infection probability of 0.849, (95% confidence interval=0.643, 0.946), highlighting uncertainty in the system, even under experimental conditions. Guided by these results, we tested the brains and hearts of wild Ross's Geese (Chen rossii, n=50) and Lesser Snow Geese (Chen caerulescens, n=50) from Karrak Lake, Nunavut, Canada. We detected 51 suspected positive tissue samples from 33 wild geese using real-time PCR with melt-curve analysis. The wild goose prevalence estimates generated by our multi-scale occupancy analysis were higher than the naïve estimates of prevalence, indicating that multiple PCR repetitions on the same organs and testing more than one organ could improve T. gondii detection. Genetic characterisation revealed Type III T. gondii alleles in six wild geese and Sarcocystis spp. in 25 samples. Our study demonstrates that Arctic nesting geese are capable of harbouring T. gondii in their tissues and could transport the parasite from their southern overwintering grounds into the Arctic region. We demonstrate how a multi-scale occupancy framework can be used in a domestic animal model to guide resource-limited sample collection and tissue analysis in wildlife. Secondly, we confirm the value of traditional occupancy in optimising T. gondii detection probability in tissue samples.

Quantitative evidence for the effects of multiple drivers on continental-scale amphibian declines.

Since amphibian declines were first proposed as a global phenomenon over a quarter century ago, the conservation community has made little progress in halting or reversing these trends. The early search for a "smoking gun" was replaced with the expectation that declines are caused by multiple drivers. While field observations and experiments have identified factors leading to increased local extinction risk, evidence for effects of these drivers is lacking at large spatial scales. Here, we use observations of 389 time-series of 83 species and complexes from 61 study areas across North America to test the effects of 4 of the major hypothesized drivers of declines. While we find that local amphibian populations are being lost from metapopulations at an average rate of 3.79% per year, these declines are not related to any particular threat at the continental scale; likewise the effect of each stressor is variable at regional scales. This result - that exposure to threats varies spatially, and populations vary in their response - provides little generality in the development of conservation strategies. Greater emphasis on local solutions to this globally shared phenomenon is needed.

ESTIMATING TOXOPLASMA GONDII EXPOSURE IN ARCTIC FOXES (VULPES LAGOPUS) WHILE NAVIGATING THE IMPERFECT WORLD OF WILDLIFE SEROLOGY.

Although the protozoan parasite Toxoplasma gondii is ubiquitous in birds and mammals worldwide, the full suite of hosts and transmission routes is not completely understood, especially in the Arctic. Toxoplasma gondii occurrence in humans and wildlife can be high in Arctic regions, despite apparently limited opportunities for transmission of oocysts shed by felid definitive hosts. Arctic foxes (Vulpes lagopus) are under increasing anthropogenic and ecologic pressure, leading to population declines in parts of their range. Our understanding of T. gondii occurrence in arctic foxes is limited to only a few regions, but mortality events caused by this parasite have been reported. We investigated the exposure of arctic foxes to T. gondii in the Karrak Lake goose colony, Queen Maud Gulf Migratory Bird Sanctuary, Nunavut, Canada. Following an occupancy-modeling framework, we performed replicated antibody testing on serum samples by direct agglutination test (DAT), indirect fluorescent antibody test (IFAT), and an indirect enzyme-linked immunosorbent assay (ELISA) that can be used in multiple mammalian host species. As a metric of test performance, we then estimated the probability of detecting T. gondii antibodies for each of the tests. Occupancy estimates for T. gondii antibodies in arctic foxes under this framework were between 0.430 and 0.758. Detection probability was highest for IFAT (0.716) and lower for DAT (0.611) and ELISA (0.464), indicating that the test of choice for antibody detection in arctic foxes might be the IFAT. We document a new geographic record of T. gondii exposure in arctic foxes and demonstrate an emerging application of ecologic modeling techniques to account for imperfect performance of diagnostic tests in wildlife species.

Conserving tigers in working landscapes.

Tiger (Panthera tigris) conservation efforts in Asia are focused on protected areas embedded in human-dominated landscapes. A system of protected areas is an effective conservation strategy for many endangered species if the network is large enough to support stable metapopulations. The long-term conservation of tigers requires that the species be able to meet some of its life-history needs beyond the boundaries of small protected areas and within the working landscape, including multiple-use forests with logging and high human use. However, understanding of factors that promote or limit the occurrence of tigers in working landscapes is incomplete. We assessed the relative influence of protection status, prey occurrence, extent of grasslands, intensity of human use, and patch connectivity on tiger occurrence in the 5400 km(2) Central Terai Landscape of India, adjacent to Nepal. Two observer teams independently surveyed 1009 km of forest trails and water courses distributed across 60 166-km(2) cells. In each cell, the teams recorded detection of tiger signs along evenly spaced trail segments. We used occupancy models that permitted multiscale analysis of spatially correlated data to estimate cell-scale occupancy and segment-scale habitat use by tigers as a function of management and environmental covariates. Prey availability and habitat quality, rather than protected-area designation, influenced tiger occupancy. Tiger occupancy was low in some protected areas in India that were connected to extensive areas of tiger habitat in Nepal, which brings into question the efficacy of current protection and management strategies in both India and Nepal. At a finer spatial scale, tiger habitat use was high in trail segments associated with abundant prey and large grasslands, but it declined as human and livestock use increased. We speculate that riparian grasslands may provide tigers with critical refugia from human activity in the daytime and thereby promote tiger occurrence in some multiple-use forests. Restrictions on human-use in high-quality tiger habitat in multiple-use forests may complement existing protected areas and collectively promote the persistence of tiger populations in working landscapes.