Range dynamics of small mammals along an elevational gradient over an 80‐year interval

One expected response to observed global warming is an upslope shift of species elevational ranges. Here, we document changes in the elevational distributions of the small mammals within the Ruby Mountains in northeastern Nevada over an 80‐year interval. We quantified range shifts by comparing distributional records from recent comprehensive field surveys (2006–2008) to earlier surveys (1927–1929) conducted at identical and nearby locations. Collector field notes from the historical surveys provided detailed trapping records and locality information, and museum specimens enabled confirmation of species' identifications. To ensure that observed shifts in range did not result from sampling bias, we employed a binomial likelihood model (introduced here) using likelihood ratios to calculate confidence intervals around observed range limits. Climate data indicate increases in both precipitation and summer maximum temperature between sampling periods. Increases in winter minimum temperatures were only evident at mid to high elevations. Consistent with predictions of change associated with climate warming, we document upslope range shifts for only two mesic‐adapted species. In contrast, no xeric‐adapted species expanded their ranges upslope. Rather, they showed either static distributions over time or downslope contraction or expansion. We attribute these unexpected findings to widespread land‐use driven habitat change at lower elevations. Failure to account for land‐use induced changes in both baseline assessments and in predicting shifts in species distributions may provide misleading objectives for conservation policies and management practices.     

Temperature, but not productivity or geometry, predicts elevational diversity gradients in ants across spatial grains

Aim  This research aims to understand the factors that shape elevational diversity gradients and how those factors vary with spatial grain. Specifically, we test the predictions of the species–productivity hypothesis, species–temperature hypothesis, the metabolic theory of ecology and the mid-domain effects null model. We also examine how the effects of productivity and temperature on richness depend on spatial grain.
Location  Deciduous forests along an elevational gradient in Great Smoky Mountains National Park, USA.
Methods  We sampled 22 leaf litter ant assemblages at three spatial grains, from 1-m² quadrats to 50 × 50 m plots using Winkler samplers.
Results  Across spatial grains, warmer sites had more species than did cooler sites, and primary productivity did not predict ant species richness. We found some support for the predictions of the metabolic theory of ecology, but no support for the mid-domain effects null model. Thus, our data are best explained by some version of a species–temperature hypothesis.
Main conclusions  Our results suggest that temperature indirectly affects ant species diversity across spatial grains, perhaps by limiting access to resources. Warmer sites support more species because they support more individuals, thereby reducing the probability of local extinction. Many of our results from this elevational gradient agree with studies at more global scales, suggesting that some mechanisms shaping ant diversity gradients are common across scales.     

Dispersal by cattle of salt-marsh and dune species into salt-marsh and dune communities

We investigated the spatial cohort structure of Trillium populations in old-growth cove forests in Great Smoky Mountains National Park (TN, USA). We mapped the locations of all Trillium erectum L., Trillium grandiflorum (Michaux) Salisbury, and Trillium vaseyi Harbison occurring within two 10 × 10 m sample plots at each of three old-growth sites—Anthony Creek, Cove Mountain, and Kalanu Prong. The height and life stage of each individual were recorded and a randomly selected subset was excavated for age determination. Our results suggest that Trillium populations in cove forests of the southern Appalachians display a high degree of spatial aggregation and are relatively stable, spatially, over long time periods (i.e., decades). Individual patches (aggregations of plants) within populations were typically multi-aged and no clear spatial cohort structure was observed. Surprisingly, more isolated plants (distal from large aggregations) were among the oldest plants in the population, rather than recent colonists dispersing away from parent populations. Individual species were less mingled than expected given that they share a common dispersal agent (ants). This study provides a double-baseline for Trillium population structure in old, primary forests with low browse pressure.     

Degrees of compositional shift in tree communities vary along a gradient of temperature change rates over one decade: Application of an individual‐based temporal beta‐diversity concept

Temporal patterns in communities have gained widespread attention recently, to the extent that temporal changes in community composition are now termed “temporal beta‐diversity.” Previous studies of beta‐diversity have made use of two classes of dissimilarity indices: incidence‐based (e.g., Sørensen and Jaccard dissimilarity) and abundance‐based (e.g., Bray–Curtis and Ružička dissimilarity). However, in the context of temporal beta‐diversity, the persistence of identical individuals and turnover among other individuals within the same species over time have not been considered, despite the fact that both will affect compositional changes in communities. To address this issue, I propose new index concepts for beta‐diversity and the relative speed of compositional shifts in relation to individual turnover based on individual identity information. Individual‐based beta‐diversity indices are novel dissimilarity indices that consider individual identity information to quantitatively evaluate temporal change in individual turnover and community composition. I applied these new indices to individually tracked tree monitoring data in deciduous and evergreen broad‐leaved forests across the Japanese archipelago with the objective of quantifying the effect of climate change trends (i.e., rates of change in both annual mean temperature and annual precipitation) on individual turnover and compositional shifts at each site. A new index explored the relative contributions of mortality and recruitment processes to temporal changes in community composition. Clear patterns emerged showing that an increase in the temperature change rate facilitated the relative contribution of mortality components. The relative speed of compositional shift increased with increasing temperature change rates in deciduous forests but decreased with increasing warming rates in evergreen forests. These new concepts provide a way to identify novel and high‐resolution temporal patterns in communities.     

Black Bear Movement and Food Conditioning in an Exurban Landscape

Conflicts between humans and wildlife have become increasingly important challenges for resource managers along the urban-wildland interface. Food conditioning (i.e., reliance by an animal on anthropogenic foods) of American black bears (Ursus americanus) is related to conflict behavior (i.e., being bold or aggressive toward humans, consuming human food or garbage, causing property damage) and often occurs in communities adjacent to Great Smoky Mountains National Park (GRSM or Park), USA. The goal of our study was to evaluate black bear space use in GRSM and in exurban areas on surrounding private lands and to identify factors associated with food conditioning and conflict behavior. We radio-collared 53 bears (29 males, 24 females) from 2015 to 2017 to compare space use characteristics and used carbon isotopic signatures (δ¹³C) from bear hair to assess food conditioning. We then performed an integrated step selection function (iSSF) analysis to characterize and compare movement and resource use as related to food conditioning. Based on the stable isotope analyses, 24 bears were classified as food conditioned (FC; 16 males and 8 females) and 37 were not food conditioned (NFC; 14 males and 23 females). Annual 95% kernel density estimate (KDE) home ranges and 50% KDE core area estimates of female and male bears did not differ by level of food conditioning (i.e., mean δ¹³C), but 95% and 50% home ranges of FC females were smaller than NFC females when data from 2015, a year of food scarcity and abnormally large home ranges, were excluded. The mean proportion of exurban development (e.g., roads, buildings, openings) within 95% KDE and 50% KDE home ranges of females increased with mean δ¹³C (i.e., greater food conditioning). The iSSF models indicated that FC bears were more likely to use forest openings associated with higher levels of development than NFC bears. We used those models to demonstrate how landscape modifications can reduce bear use of exurban areas, particularly for NFC bears. Our stable isotope, movement, and resource use data indicate that conflict behaviors displayed by many bears within GRSM were learned in areas outside Park boundaries. © 2020 The Wildlife Society.     

A new cost‐distance model for human accessibility and an evaluation of accessibility bias in permanent vegetation plots in Great Smoky Mountains National Park,USA

Question: Can a new cost‐distance model help us to evaluate the potential for accessibility bias in ecological observations? How much accessibility bias is present in the vegetation monitoring plots accumulated over the last three decades in Great Smoky Mountains National Park? Location: Great Smoky Mountains National Park, North Carolina and Tennessee, USA. Methods: Distance, slope, stream crossings, and vegetation density were incorporated into a least‐cost model of energetic expenditure for human access to locations. Results: Estimated round‐trip energy costs for the park ranged from 0 to 1.62 × 10⁵ J kg^?1. The estimated round‐trip energetic expenditure for the surveys ranged from 53 to 1.51 × 10⁵ J kg^?1. Their distribution was more accessible than the random expectation. Ten (17%) of the vegetation types in the park are significantly under‐sampled relative to their area, and 16 (29%) are over‐sampled. Plots in 18 of the 40 vegetation types exhibited a significant positive correlation with accessibility. Conclusions: The least‐cost model is an improvement over previous attempts to quantify accessibility. The bias in plot locations suggests using a least‐cost model to test for bias in cases in which human accessibility is confounded with other sources of ecosystem variation.     

Recovery of Sasajiscymnus tsugae,released against hemlock woolly adelgid,Adelges tsugae,in the southern Appalachians

Eastern hemlock in the Great Smoky Mountains National Park is currently threatened by the hemlock woolly adelgid, Adelges tsugae Annand (Hemiptera: Adelgidae). As part of a management plan against this invasive insect pest, about 350,000 adults of the predatory beetle Sasajiscymnus tsugae (Sasaji and McClure) (Coleoptera: Coccinellidae) were released at ca. 150 sites in the Park from 2002 to 2007. Of these adult release sites, 33 were sampled in 2008 and 2009 using beat-sheet sampling for 4 man-hours. Sasajiscymnus tsugae adults (n=78) and/or larvae (n=145) were recovered from seven sites (21.2% of the release sites sampled). Recovery of S. tsugae was significantly associated with older release sites, with the most beetles recovered from 2002 release sites. These results indicate that S. tsugae may require more time (i.e., 5–7 years) than anticipated for population densities to reach readily detectable levels in some areas.     

Temporal patterns of diversity: assessing the biotic and abiotic controls on ant assemblages

In this study, we use 12 months of data from 11 ant assemblages to test whether seasonal variation in ant diversity is governed by either the structuring influences of interspecific competition or environmental conditions. Because the importance of competition might vary along environmental gradients, we also test whether the signature of competition depends on elevation. We find little evidence that competition structures the seasonal patterns of activity in the ant assemblages considered, but find support for the effects of temperature on seasonal patterns of diversity, especially at low-elevation sites. Although, in general, both competition and the environment interact to structure ant assemblages, our results suggest that environmental conditions are the primary force structuring the seasonal activity of the ant assemblages studied here.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 191–201.     

Intraspecific trait variability and community assembly in hawkmoths (Lepidoptera: Sphingidae) across an elevational gradient in the eastern Himalayas,India

1. We investigated some aspects of hawkmoth community assembly at 13 elevations along a 200‐ to 2770‐m transect in the eastern Himalayas, a little studied biodiversity hot spot of global importance. We measured the morphological traits of body mass, wing loading, and wing aspect ratio of 3,301 free‐ranging individuals of 76 species without having to collect or even constrain them. We used these trait measurements and T‐statistic metrics to assess the strength of intracommunity (“internal") and extra‐community (“external”) filters which determine the composition of communities vis‐a‐vis the regional pool of species.
2. The trait distribution of constituent species turned out to be nonrandom subsets of the community‐trait distribution, providing strong evidence for internal filtering in all elevational communities. The external filter metric was more ambiguous. However, the elevational dependence of many metrics including that of the internal filter provided evidence for external (i.e., environmental) filtering. On average, a species occupied as much as 50%–75% of the total community‐trait space, yet the T‐statistic metric for internal filter was sufficiently sensitive to detect a strong nonrandom structure in the trait distribution.
3. We suggest that the change in T‐statistic metrics along the environmental gradient may provide more clues to the process of community assembly than previously envisaged. A large, smoothly varying and well‐sampled environmental span would make it easier to discern them. Developing T‐statistics for combined analysis of multiple traits will perhaps provide a more accurate picture of internal/filtering and niche complementarity. Moths are a hyperdiverse taxon and a very important component of many ecosystems. Our technique for accurately measuring body and wing dimensions of free‐ranging moths can generate trait database for a large number of individuals in a time‐ and resource‐efficient manner for a variety of community assembly studies using this important taxon.

     

Effects of experimental watering but not warming on herbivory vary across a gradient of precipitation

Climate change can affect biotic interactions, and the impacts of climate on biotic interactions may vary across climate gradients. Climate affects biotic interactions through multiple drivers, although few studies have investigated multiple climate drivers in experiments. We examined the effects of experimental watering, warming, and predator access on leaf water content and herbivory rates of woolly bear caterpillars (Arctia virginalis) on a native perennial plant, pacific silverweed (Argentina anserina ssp. pacifica), at two sites across a gradient of precipitation in coastal California. Based on theory, we predicted that watering should increase herbivory at the drier end of the gradient, predation should decrease herbivory, and watering and warming should have positive interacting effects on herbivory. Consistent with our predictions, we found that watering only increased herbivory under drier conditions. However, watering increased leaf water content at both wetter and drier sites. Warming increased herbivory irrespective of local climate and did not interact with watering. Predation did not affect herbivory rates. Given predictions that the study locales will become warmer and drier with climate change, our results suggest that the effects of future warming and drying on herbivory may counteract each other in drier regions of the range of Argentina anserina. Our findings suggest a useful role for range‐limit theory and the stress‐gradient hypothesis in predicting climate change effects on herbivory across stress gradients. Specifically, if climate change decreases stress, herbivory may increase, and vice versa for increasing stress. In addition, our work supports previous suggestions that multiple climate drivers are likely to have dampening effects on biotic interactions due to effects in different directions, though this is context‐dependent.     

The distributions of morphologically specialized hummingbirds coincide with floral trait matching across an Andean elevational gradient

Morphological trait matching between species affects resource partitioning in mutualistic systems. Yet, the determinants of spatial variation in trait matching remain largely unaddressed. Here, we generate a hypothesis that is based on the geographical distributions of species morphologies. To illustrate our hypothesis, as a study system we use hummingbirds in the tropical Andes. Hummingbirds with specialized morphologies (i.e., long or curved bills) may forage on flowers that are inaccessible to hummingbirds with generalized bill morphologies (i.e., small‐to‐medium‐sized bills with no curvature), yet the vast majority of hummingbirds have generalized bill morphologies. Thus, we propose that trait matching across space is determined by the distribution of morphological specialists. In the Andes, we observe the richness of specialized hummingbird morphotypes to peak at high and low elevations. Therefore, we hypothesize that trait matching should be most influential in predicting pairwise interactions at high and low elevations. We illustrate our hypothesis by field observations along an elevational gradient in Podocarpus National Park (Ecuador). Using Bayesian hierarchical modeling of interaction frequencies in combination with network analyzes, we found that hummingbirds at high and low elevations contributed to resource partitioning by foraging on morphologically close‐matching flowers. Moreover, at high and low elevations, hummingbirds with specialized morphologies showed a stronger tendency to visit close‐matching flowers than morphological non‐specialists did. In contrast, at mid‐elevations, hummingbirds were not attracted to morphologically close‐matching flowers. These results suggest that the spatial distribution of specialized morphotypes determines trait matching and the partitioning of interactions within hummingbird–plant communities. Abstract in Spanish is available with online material.     

Sharing detection heterogeneity information among species in community models of occupancy and abundance can strengthen inference

1. The estimation of abundance and distribution and factors governing patterns in these parameters is central to the field of ecology. The continued development of hierarchical models that best utilize available information to inform these processes is a key goal of quantitative ecologists. However, much remains to be learned about simultaneously modeling true abundance, presence, and trajectories of ecological communities.
2. Simultaneous modeling of the population dynamics of multiple species provides an interesting mechanism to examine patterns in community processes and, as we emphasize herein, to improve species‐specific estimates by leveraging detection information among species. Here, we demonstrate a simple but effective approach to share information about observation parameters among species in hierarchical community abundance and occupancy models, where we use shared random effects among species to account for spatiotemporal heterogeneity in detection probability.
3. We demonstrate the efficacy of our modeling approach using simulated abundance data, where we recover well our simulated parameters using N‐mixture models. Our approach substantially increases precision in estimates of abundance compared with models that do not share detection information among species. We then expand this model and apply it to repeated detection/non‐detection data collected on six species of tits (Paridae) breeding at 119 1 km² sampling sites across a P. montanus hybrid zone in northern Switzerland (2004–2020). We find strong impacts of forest cover and elevation on population persistence and colonization in all species. We also demonstrate evidence for interspecific competition on population persistence and colonization probabilities, where the presence of marsh tits reduces population persistence and colonization probability of sympatric willow tits, potentially decreasing gene flow among willow tit subspecies.
4. While conceptually simple, our results have important implications for the future modeling of population abundance, colonization, persistence, and trajectories in community frameworks. We suggest potential extensions of our modeling in this paper and discuss how leveraging data from multiple species can improve model performance and sharpen ecological inference.

     

Distinct responses and range shifts of lizard populations across an elevational gradient under climate change

Ongoing climate change has profoundly affected global biodiversity, but its impacts on populations across elevations remain understudied. Using mechanistic niche models incorporating species traits, we predicted ecophysiological responses (activity times, oxygen consumption and evaporative water loss) for lizard populations at high-elevation (<3600 m asl) and extra-high-elevation (≥3600 m asl) under recent (1970–2000) and future (2081–2100) climates. Compared with their high-elevation counterparts, lizards from extra-high-elevation are predicted to experience a greater increase in activity time and oxygen consumption. By integrating these ecophysiological responses into hybrid species distribution models (HSDMs), we were able to make the following predictions under two warming scenarios (SSP1-2.6, SSP5-8.5). By 2081–2100, we predict that lizards at both high- and extra-high-elevation will shift upslope; lizards at extra-high-elevation will gain more and lose less habitat than will their high-elevation congeners. We therefore advocate the conservation of high-elevation species in the context of climate change, especially for those populations living close to their lower elevational range limits. In addition, by comparing the results from HSDMs and traditional species distribution models, we highlight the importance of considering intraspecific variation and local adaptation in physiological traits along elevational gradients when forecasting species' future distributions under climate change.     

Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization,Ascaphus truei,resolved using two types of genetic markers

Reconstruction of historical relationships between geographic regions within a species’ range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques—microsatellite and genotype‐by‐sequencing (GBS)—to compare the within‐region genetic diversity of populations near the northern extent of the species’ range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.     

Weather variation affects the dispersal of grasshoppers beyond their elevational ranges

1. Understanding how abiotic conditions influence dispersal patterns of organisms is important for understanding the degree to which species can track and persist in the face of changing climate.
2. The goal of this study was to understand how weather conditions influence the dispersal pattern of multiple nonmigratory grasshopper species from lower elevation grassland habitats in which they complete their life‐cycles to higher elevations that extend beyond their range limits.
3. Using over a decade of weekly spring to late‐summer field survey data along an elevational gradient, we explored how abundance and richness of dispersing grasshoppers were influenced by temperature, precipitation, and wind speed and direction. We also examined how changes in population sizes at lower elevations might influence these patterns.
4. We observed that the abundance of dispersing grasshoppers along the gradient declined 4‐fold from the foothills to the subalpine and increased with warmer conditions and when wind flow patterns were mild or in the downslope direction. Thirty‐eight unique grasshopper species from lowland sites were detected as dispersers across the survey years, and warmer years and weak upslope wind conditions also increased the richness of these grasshoppers. The pattern of grasshoppers along the gradient was not sex biased. The positive effect of temperature on dispersal rates was likely explained by an increase in dispersal propensity rather than by an increase in the density of grasshoppers at low elevation sites.
5. The results of this study support the hypothesis that the dispersal patterns of organisms are influenced by changing climatic conditions themselves and as such, that this context‐dependent dispersal response should be considered when modeling and forecasting the ability of species to respond to climate change.

     

No long-term effect of black bear removal on elk calf recruitment in the southern Appalachians

In 2001 and 2002, 52 elk (Cervus canadensis; 21 males, 31 females), originally obtained from Elk Island National Park, Alberta, Canada, were transported and released into Cataloochee Valley in the northeastern portion of Great Smoky Mountains National Park (GRSM, Park), North Carolina, USA. The annual population growth rate (λ) was negative (0.996, 95% CI = 0.945–1.047) and predation by black bears (Ursus americanus) on elk calves was identified as an important determinant of population growth. From 2006 to 2008, 49 bears from the primary elk calving area (i.e., Cataloochee Valley) were trapped and translocated about 70 km to the southwestern portion of the Park just prior to elk calving. Per capita recruitment (i.e., the number of calves produced per adult female that survive to 1 year of age) increased from 0.306 prior to bear translocation (2001–2005) to 0.544 during years when bears were translocated (2006–2008) and λ increased to 1.118 (95% CI = 1.096–1.140). Our objective was to determine whether per capita calf recruitment rates after bear removal (2009–2019) at Cataloochee were similar to the higher rates estimated during bear removal (i.e., long-term response) or if they returned to rates before bear removal (i.e., short-term response), and how those rates compared with recruitment from portions of our study area where bears were not relocated. We documented 419 potential elk calving events and monitored 129 yearling and adult elk from 2001 to 2019. Known-fate models based on radio-telemetry and observational data supported calf recruitment returning to pre-2006 levels at Cataloochee (short-term response); recruitment of Cataloochee elk before and after bear relocation was lower (0.184) than during bear relocation (0.492). Recruitment rates of elk outside the removal area during the bear relocation period (0.478) were similar to before and after rates (0.420). In the Cataloochee Valley, cause-specific annual calf mortality rates due to predation by bears were 0.319 before, 0.120 during, and 0.306 after bear relocation. In contrast, the cause-specific annual mortality rate of calves in areas where bears were not relocated was 0.033 after the bear relocation period, with no bear predation on calves before or during bear relocation. The mean annual population growth rate for all monitored elk was 1.062 (95% CI = 0.979–1.140) after bear relocation based on the recruitment and survival data. Even though the effects of bear removal were temporary, the relocations were effective in achieving a short-term increase in elk recruitment, which was important for the reintroduction program given that the elk population was small and vulnerable to extirpation.