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.

Hidden Markov models for extended batch data.

Batch marking provides an important and efficient way to estimate the survival probabilities and population sizes of wild animals. It is particularly useful when dealing with animals that are difficult to mark individually. For the first time, we provide the likelihood for extended batch-marking experiments. It is often the case that samples contain individuals that remain unmarked, due to time and other constraints, and this information has not previously been analyzed. We provide ways of modeling such information, including an open N-mixture approach. We demonstrate that models for both marked and unmarked individuals are hidden Markov models; this provides a unified approach, and is the key to developing methods for fast likelihood computation and maximization. Likelihoods for marked and unmarked individuals can easily be combined using integrated population modeling. This allows the simultaneous estimation of population size and immigration, in addition to survival, as well as efficient estimation of standard errors and methods of model selection and evaluation, using standard likelihood techniques. Alternative methods for estimating population size are presented and compared. An illustration is provided by a weather-loach data set, previously analyzed by means of a complex procedure of constructing a pseudo likelihood, the formation of estimating equations, the use of sandwich estimates of variance, and piecemeal estimation of population size. Simulation provides general validation of the hidden Markov model methods developed and demonstrates their excellent performance and efficiency. This is especially notable due to the large numbers of hidden states that may be typically required.

Adjusting for undercoverage of access-points in creel surveys with fewer overflights.

Creel surveys are used in recreational fisheries to estimate angling effort, catch, and harvest. Aerial-access creel surveys rely on two components: (1) a ground component in which fishing parties returning from their trips are interviewed at some access-points of the fishery; (2) an aerial component in which the number of fishing parties is counted. A common practice is to sample fewer aerial survey days than ground survey days. This is thought by practitioners to reduce the cost of the survey, but there is a lack of sound statistical methodology for this case. In this article, we propose various estimation methods to handle this situation and evaluate their asymptotic properties from a design-based perspective. We also propose formulas for the optimal allocation of the effort between the ground and the aerial portion of the survey, for given costs and budget. A simulation study investigates the performance of the estimators. Finally, we apply our methods to data from an annual Kootenay Lake survey (Canada).

Resource selection by the California condor (Gymnogyps californianus) relative to terrestrial-based habitats and meteorological conditions.

Condors and vultures are distinct from most other terrestrial birds because they use extensive soaring flight for their daily movements. Therefore, assessing resource selection by these avian scavengers requires quantifying the availability of terrestrial-based habitats, as well as meteorological variables that influence atmospheric conditions necessary for soaring. In this study, we undertook the first quantitative assessment of habitat- and meteorological-based resource selection in the endangered California condor (Gymnogyps californianus) within its California range and across the annual cycle. We found that condor use of terrestrial areas did not change markedly within the annual cycle, and that condor use was greatest for habitats where food resources and potential predators could be detected and where terrain was amenable for taking off from the ground in flight (e.g., sparse habitats, coastal areas). Condors originating from different release sites differed in their use of habitat, but this was likely due in part to variation in habitats surrounding release sites. Meteorological conditions were linked to condor use of ecological subregions, with thermal height, thermal velocity, and wind speed having both positive (selection) and negative (avoidance) effects on condor use in different areas. We found little evidence of systematic effects between individual characteristics (i.e., sex, age, breeding status) or components of the species management program (i.e., release site, rearing method) relative to meteorological conditions. Our findings indicate that habitat type and meteorological conditions can interact in complex ways to influence condor resource selection across landscapes, which is noteworthy given the extent of anthropogenic stressors that may impact condor populations (e.g., lead poisoning, wind energy development). Additional studies will be valuable to assess small-scale condor movements in light of these stressors to help minimize their risk to this critically endangered species.

A comparison of abundance estimates from extended batch-marking and Jolly-Seber-type experiments.

Little attention has been paid to the use of multi-sample batch-marking studies, as it is generally assumed that an individual's capture history is necessary for fully efficient estimates. However, recently, Huggins et al. (2010) present a pseudo-likelihood for a multi-sample batch-marking study where they used estimating equations to solve for survival and capture probabilities and then derived abundance estimates using a Horvitz-Thompson-type estimator. We have developed and maximized the likelihood for batch-marking studies. We use data simulated from a Jolly-Seber-type study and convert this to what would have been obtained from an extended batch-marking study. We compare our abundance estimates obtained from the Crosbie-Manly-Arnason-Schwarz (CMAS) model with those of the extended batch-marking model to determine the efficiency of collecting and analyzing batch-marking data. We found that estimates of abundance were similar for all three estimators: CMAS, Huggins, and our likelihood. Gains are made when using unique identifiers and employing the CMAS model in terms of precision; however, the likelihood typically had lower mean square error than the pseudo-likelihood method of Huggins et al. (2010). When faced with designing a batch-marking study, researchers can be confident in obtaining unbiased abundance estimators. Furthermore, they can design studies in order to reduce mean square error by manipulating capture probabilities and sample size.

Predicting grizzly bear density in western North America.

Conservation of grizzly bears (Ursus arctos) is often controversial and the disagreement often is focused on the estimates of density used to calculate allowable kill. Many recent estimates of grizzly bear density are now available but field-based estimates will never be available for more than a small portion of hunted populations. Current methods of predicting density in areas of management interest are subjective and untested. Objective methods have been proposed, but these statistical models are so dependent on results from individual study areas that the models do not generalize well. We built regression models to relate grizzly bear density to ultimate measures of ecosystem productivity and mortality for interior and coastal ecosystems in North America. We used 90 measures of grizzly bear density in interior ecosystems, of which 14 were currently known to be unoccupied by grizzly bears. In coastal areas, we used 17 measures of density including 2 unoccupied areas. Our best model for coastal areas included a negative relationship with tree cover and positive relationships with the proportion of salmon in the diet and topographic ruggedness, which was correlated with precipitation. Our best interior model included 3 variables that indexed terrestrial productivity, 1 describing vegetation cover, 2 indices of human use of the landscape and, an index of topographic ruggedness. We used our models to predict current population sizes across Canada and present these as alternatives to current population estimates. Our models predict fewer grizzly bears in British Columbia but more bears in Canada than in the latest status review. These predictions can be used to assess population status, set limits for total human-caused mortality, and for conservation planning, but because our predictions are static, they cannot be used to assess population trend.

Smoothing population size estimates for time-stratified mark-recapture experiments using Bayesian P-splines.

Petersen-type mark-recapture experiments are often used to estimate the number of fish or other animals in a population moving along a set migration route. A first sample of individuals is captured at one location, marked, and returned to the population. A second sample is then captured farther along the route, and inferences are derived from the numbers of marked and unmarked fish found in this second sample. Data from such experiments are often stratified by time (day or week) to allow for possible changes in the capture probabilities, and previous methods of analysis fail to take advantage of the temporal relationships in the stratified data. We present a Bayesian, semiparametric method that explicitly models the expected number of fish in each stratum as a smooth function of time. Results from the analysis of historical data from the migration of young Atlantic salmon (Salmo salar) along the Conne River, Newfoundland, and from a simulation study indicate that the new method provides more precise estimates of the population size and more accurate estimates of uncertainty than the currently available methods.

The Jolly-Seber model with tag loss.

Tag loss in mark-recapture experiments is a violation of one of the Jolly-Seber model assumptions. It causes bias in parameter estimates and has only been dealt with in an ad hoc manner. We develop methodology to estimate tag retention and abundance in double-tagging mark-recapture experiments. We apply this methodology to walleyes (Stizostedion vitreum) in Mille Lacs, Minnesota.

Capture-recapture studies using radio telemetry with premature radio-tag failure.

Radio tags, because of their high detectability, are often used in capture-recapture studies. A key assumption is that radio tags do not cease functioning during the study. Radio-tag failure before the end of a study can lead to underestimates of survival rates. We develop a model to incorporate secondary radio-tag failure data. This model was applied to chinook smolts (Oncorhynchus tshawytscha) on the Columbia River, Washington. Estimates of fish survival from this model were much larger than those from the standard Cormack-Jolly-Seber analysis.