Climate Change Habitat Model Forecasts for Eight Owl Species in the Southwestern US.

The high-resolution forecasting of vegetation type shifts may prove essential in anticipating and mitigating the impacts of future climate change on bird populations. Here, we used the US Forest Service Ecological Response Unit (ERU) classification to develop and assess vegetation-based breeding habitat profiles for eight owl species occurring in the foothills and mountains of the Southwestern US. Shifts in mapped habitat were forecast using an ecosystem vulnerability model based on the pre-1990 climate envelopes of ERUs and the Intergovernmental Panel on Climate Change's (IPCC) A1B moderate-emission scenario for the future climate. For five of the eight owl species, the regional breeding habitat extent was projected to decline by at least 60% by 2090. Three species, the boreal owl (Aegolius funereus; at the trailing edge of its distribution), flammulated owl (Psiloscops flammeolus), and northern pygmy-owl (Glaucidium gnoma), were projected to experience the steepest habitat loss rates of 85%, 85%, and 76%, respectively. Projected vegetation shifts overlaid with well-documented flammulated owl breeding populations showed the complete or near complete loss of habitat by 2090 in areas of montane forest currently supporting dense aggregations of owl territories. Generalist or lower-elevation owl species were predicted to be less impacted, while, for the whiskered screech-owl (Megascops trichopsis), the contraction of the current habitat was nearly offset by a projected northward expansion. In general, the results of this study suggest high exposure to climate change impacts for the upper-elevation forest owls of semi-arid Southwestern North America. Long-distance migration and low natal philopatry may prove important to some montane owl populations in adapting to the regional loss of habitat.

Comparative spatially explicit approach for testing effects of soil chemicals on terrestrial wildlife bioindicator demographics.

Wildlife species are often used as bioindicators to evaluate the extent and severity of environmental contamination and the effectiveness of remediation practices. A common approach for investigating population- or community-level impacts on bioindicators compares demographic parameter estimates (e.g., population size or density) between sites that were subjected to different levels of contamination. However, the traditional analytical method used in such studies is nonspatial capture-recapture, which results in conclusions about potential relationships between demographics and contaminants being inferred indirectly. Here, we extend this comparative approach to the spatially explicit framework, allowing direct estimation of said relationships and comparisons between study areas, by applying spatial capture-recapture (SCR) models to bioindicator (deer mice [Peromyscus spp.]) detection data from two study areas that were subjected to different industrial activities and remediation practices. Bioindicator density differed by 178% between the neighboring study areas, and the area with the highest soil concentrations of polychlorinated biphenyls, chromium, and zinc had the highest bioindicator density. Under the traditional nonspatial approach, we might have concluded that soil chemical levels had negligible influences on demographics. However, by modeling density as a spatial function of select chemical concentrations using SCR models, we found strong support for a positive relationship between density and soil chromium concentrations in one study area (ß = 0.82), which was not masked by or associated with habitat-related metrics. To obtain reliable inferences about potential effects of environmental contamination on bioindicator demographics, we contend that a comparative spatially explicit approach using SCR ought to become standard.

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate.

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

Comparing western (Megascops kennicottii) and whiskered (M. trichopsis) screech-owl microbiomes in southern Arizona using a novel 16S rRNA sequencing method.

Microbiomes are essential to a host's physiology and health. Despite the overall importance of microbiomes to animal health, they remain understudied in wildlife. Microbiomes function as physical barriers to invading pathogens, and changes in the diversity or composition of microbes within a host may disrupt this barrier. In order to use microbiomes in wildlife ecology, knowledge of the natural variation within and among species is essential. We compare the diversity and composition of two avian species that share the same habitat and niche in our study area, the western screech-owl (Megascops kennicottii) and the whiskered screech-owl (M. trichopsis). We used a targeted 16S sequencing method to improve the taxonomic resolution of microbiomes. We found similar measures of alpha diversity between species and sample types (cloacal samples vs. fecal samples). However, there were significant differences in bacterial species richness among nestlings from different nest boxes, and the composition differed between the two bird species and among nestlings from different nest boxes. Western screech-owls had more variation in alpha diversity and composition and had fewer bacterial species in their core microbiome than whiskered screech-owls. Siblings are likely to yield similar findings for microbiomes; thus, sampling nestlings from different nests may be most informative for monitoring population-level changes.

Does Age, Residency, or Feeding Guild Coupled with a Drought Index Predict Avian Health during Fall Migration?

Birds are good indicators of environmental change and are often studied for responses to climate. Many studies focus on breeding birds, while fewer look at the migration period, which is a critical time for many birds. Birds are more susceptible to unusual climatic events during their migration due to the metabolic stress of long-distance movements. In the fall of 2020, an unusual cold weather event coupled with drought and wildfire smoke led to a large avian mortality event in New Mexico. Later analysis pointed to the mortality being largely due to starvation. This was the impetus for our research. We used 11 years of fall bird banding data from two locations, along with local drought indices, to determine what predicts avian health during the migration period. We used fat score data from over 15,000 individual birds to assess whether drought indices, age, diet, or residency influenced avian health using multiple logistic regression. We found that the probability of positive fat scores decreased as drought severity increased for younger, insectivorous, migratory birds. Insectivores had a higher probability of receiving a fat score greater than zero relative to local drought conditions, which is important, since many North American insectivores are in steep decline. Migratory birds showed a greater response than year-round residents, and older birds showed a lower but significant response compared to hatch-year birds. Our results suggest that migratory insectivores in the southwestern United States may be less resilient to drought-related climate change.

Individual Nest Site Preferences Do Not Explain Upslope Population Shifts of a Secondary Cavity-Nesting Species.

Geographic ranges of plants and animals are shifting due to environmental change. While some species are shifting towards the poles and upslope in elevation, the processes leading to these patterns are not well known. We analyzed 22 years of western bluebird (Sialia mexicana) data from a large nest box network in northern New Mexico at elevations between 1860 m and 2750 m. This population has shifted to higher elevations over time, but whether this is due to changes in nesting behavior and preference for higher elevation within the population or driven by immigration is unclear. We banded adults and nestlings from nest boxes and examined nesting location and elevation for individual birds captured two or more times. Most recaptured birds nested at the same nest boxes in subsequent years, and the number of birds that moved upslope did not significantly differ from the number that moved downslope. Fledglings moved greater distances and elevations than adults, but these movements were not upslope specific. Female fledglings showed greater changes in elevation and distance compared to male fledglings, but again, movements were not consistently upslope. The upslope shift in this population may be due to birds immigrating into the population and not from changes in individual nesting behavior.

Sex ratio of Western Bluebirds Sialia mexicana is mediated by phenology and clutch size.

Mothers may produce more of one sex to maximize their fitness if there are differences in the cost of producing each sex or there are differences in their relative reproductive value. Breeding date and clutch size are known to influence offspring sex ratios in birds through sex differences in dispersal, social behaviours, differential mortality, and available food resources. We tested if breeding date, clutch size and drought conditions influenced offspring sex ratios in a sexually size-monomorphic species, the Western Bluebird, by interrogating a 21-year dataset. After controlling for differential mortality, we found that hatch dates late in the breeding season were associated with the production of more females, suggesting that the value of producing males declines as the breeding season progresses. When clutch size was taken into account, small clutches yielded significantly more females late in the breeding season compared to the early and middle parts of the breeding season that produced significantly more males. Large clutches early in the season tended to produce more females, although this was not significant. Drought severity was not correlated with sex ratio adjustment. We propose and discuss several explanations for these patterns, including male offspring, but not female offspring, acting as helpers, increased female nestling provisioning late in the breeding season, differences in food abundance, and egg-laying order. Future work will help to uncover the mechanisms leading to these patterns. Identifying patterns and mechanisms of sex ratio skew from long-term datasets is important for informing predictions regarding life-history trade-offs in wildlife populations.

Investigating effects of soil chemicals on density of small mammal bioindicators using spatial capture-recapture models.

Monitoring the ecological impacts of environmental pollution and the effectiveness of remediation efforts requires identifying relationships between contaminants and the disruption of biological processes in populations, communities, or ecosystems. Wildlife are useful bioindicators, but traditional comparative experimental approaches rely on a staunch and typically unverifiable assumption that, in the absence of contaminants, reference and contaminated sites would support the same densities of bioindicators, thereby inferring direct causation from indirect data. We demonstrate the utility of spatial capture-recapture (SCR) models for overcoming these issues, testing if community density of common small mammal bioindicators was directly influenced by soil chemical concentrations. By modeling density as an inhomogeneous Poisson point process, we found evidence for an inverse spatial relationship between Peromyscus density and soil mercury concentrations, but not other chemicals, such as polychlorinated biphenyls, at a site formerly occupied by a nuclear reactor. Although the coefficient point estimate supported Peromyscus density being lower where mercury concentrations were higher (ß = -0.44), the 95% confidence interval overlapped zero, suggesting no effect was also compatible with our data. Estimated density from the most parsimonious model (2.88 mice/ha; 95% CI = 1.63-5.08), which did not support a density-chemical relationship, was within the range of reported densities for Peromyscus that did not inhabit contaminated sites elsewhere. Environmental pollution remains a global threat to biodiversity and ecosystem and human health, and our study provides an illustrative example of the utility of SCR models for investigating the effects that chemicals may have on wildlife bioindicator populations and communities.

Long-term phenology of two North American secondary cavity-nesters in response to changing climate conditions.

Wildlife populations can respond to changes in climate conditions by either adapting or moving to areas with preferred climate regimes. We studied nesting responses of two bird species, western bluebird (Sialia mexicana) and ash-throated flycatcher (Myiarchus cinerascens), to changing climate conditions (i.e., rising temperatures and increased drought stress) over 21 years in northern New Mexico. We used data from 1649 nests to assess whether the two species responded to changing climate conditions through phenological shifts in breeding time or shifts in nesting elevation. We also examined changes in reproductive output (i.e., clutch size). Our data show that western bluebirds significantly increased nesting elevation over a 19-year period by approximately 5 m per year. Mean spring temperature was the best predictor of western bluebird nesting elevation. Higher nesting elevations were not correlated with hatch dates or clutch sizes in western bluebirds, suggesting that nesting at higher elevations does not affect breeding time or reproductive output. We did not observe significant changes in nesting elevation or breeding dates in ash-throated flycatchers. Nesting higher in elevation may allow western bluebirds to cope with the increased temperatures and droughts. However, this climate niche conservatism may pose a risk for the conservation of the species if climate change and habitat loss continue to occur. The lack of significant changes detected in nesting elevation, breeding dates, and reproductive output in ash-throated flycatchers suggests a higher tolerance for changing environmental conditions in this species. This is consistent with the population increases reported for flycatchers in areas experiencing dramatic climate changes.