The Influence of Sex and Season on Conspecific Spatial Overlap in a Large, Actively-Foraging Colubrid Snake

Home Range of the Caspian Whipsnake Dolichophis caspius (Gmelin, 1789) in a Threatened Peri-Urban Population

Semi-natural environments within cities can provide habitats for vulnerable reptile species. Better understanding of their habitat use and home range sizes is important for their conservation. We investigated the spatial ecology of Caspian whipsnakes (Dolicophis caspius) in a peri-urban habitat in Budapest, Hungary. We used radiotelemetry to track five adult snakes and analyzed their microhabitat preferences, home range sizes and daily movements. The Caspian whipsnakes intensively utilized areas covered with woody vegetation, with a high density of hibernacula. The tracked snakes used an area of 40.15 ha during the activity period from spring to autumn, but for the winter, they withdrew to a central area of 1.75 ha, abundant in hibernacula. During the activity period the individual home range sizes varied between 6.1 and 15.5 ha, estimated using the minimum convex polygon (MCP); however, for the entire datasets of the individuals, the adaptive kernel method gave the highest mean (13.8 ha), while the LoCoH-R yielded the smallest home ranges (5.19 ha). We found that the average daily displacement for the different individuals ranged between 12.6 and 36.6 m during their main activity season. In the study area, the whipsnakes currently have enough space for foraging, but the restricted spatial distribution of hibernacula, which is mainly available in the central dry rocky forest and partly in the shrubby areas, can limit the extent of the suitable habitat. Human activities and anthropogenic disturbances, especially around hibernacula, may exert pressure on this peri-urban snake population.

Accounting for geographic variation in species-habitat associations during habitat suitability modeling

Range-wide species conservation efforts are facilitated by spatially explicit estimates of habitat suitability. However, species-environment relationships often vary geographically and models assuming geographically constant relationships may result in misleading inferences. We present the first range-wide habitat suitability model (HSM) for the federally threatened eastern indigo snake (Drymarchon couperi) as a case study illustrating an approach to account for known latitudinal variation in habitat associations. Specifically, we modeled habitat suitability using interactive relationships between minimum winter temperature and several a priori environmental covariates and compared our results to those from models assuming geographically constant relationships. We found that multi-scale models including interactive effects with winter temperature outperformed single-scale models and models not including interactive effects with winter temperature. Our top-ranked model had suitable range-wide predictive performance and identified numerous large (i.e., ≥1000 ha) potential habitat patches throughout the indigo snake range. Predictive performance was greatest in southern Georgia and northern Florida likely reflecting more restrictive indigo snake habitat associations in these regions. This study illustrates how modeling interactive effects between temperature and environmental covariates can improve the performance of HSMs across geographically varying environmental gradients.

Multiscale assessment of functional connectivity: Landscape genetics of eastern indigo snakes in an anthropogenically fragmented landscape in central Florida

Landscape features can strongly influence gene flow and the strength and direction of these effects may vary across spatial scales. However, few studies have evaluated methodological approaches for selecting spatial scales in landscape genetics analyses, in part because of computational challenges associated with optimizing landscape resistance surfaces (LRS). We used the federally threatened eastern indigo snake (Drymarchon couperi) in central Florida as a case study with which to compare the importance of landscape features and their scales of effect in influencing gene flow. We used genetic algorithms (ResistanceGA) to empirically optimize LRS using categorical land cover surfaces, multiscale resource selection surfaces (RSS), and four combinations of landscape covariates measured at multiple spatial scales (multisurface multiscale LRS). We compared LRS where scale was selected using pseudo- and full optimization. Multisurface multiscale LRS received more empirical support than LRS optimized from categorical land cover surfaces or RSS. Multiscale LRS with scale selected using full optimization generally outperformed those with scale selected using pseudo-optimization. Multiscale LRS with large spatial scales (1200-1800 m) received the most empirical support. Our results highlight the importance of considering landscape features across multiple spatial scales in landscape genetic analyses, particularly broad scales relative to species movement potential. Different effects of scale on home range-level movements and dispersal could explain weak associations between habitat suitability and gene flow in other studies. Our results also demonstrate the importance of large tracts of undeveloped upland habitat with heterogenous vegetation communities and low urbanization for promoting indigo snake connectivity.

Dynamics of animal joint space use: a novel application of a time series approach

BACKGROUND

Animal use is a dynamic phenomenon, emerging from the movements of animals responding to a changing environment. Interactions between animals are reflected in patterns of joint space use, which are also dynamic. High frequency sampling associated with GPS telemetry provides detailed data that capture space use through time. However, common analyses treat joint space use as static over relatively long periods, masking potentially important changes. Furthermore, linking temporal variation in interactions to covariates remains cumbersome. We propose a novel method for analyzing the dynamics of joint space use that permits straightforward incorporation of covariates. This method builds upon tools commonly used by researchers, including kernel density estimators, utilization distribution intersection metrics, and extensions of linear models.

METHODS

We treat the intersection of the utilization distributions of two individuals as a time series. The series is linked to covariates using copula-based marginal beta regression, an alternative to generalized linear models. This approach accommodates temporal autocorrelation and the bounded nature of the response variable. Parameters are easily estimated with maximum likelihood and trend and error structures can be modeled separately. We demonstrate the approach by analyzing simulated data from two hypothetical individuals with known utilization distributions, as well as field data from two coyotes (Canis latrans) responding to appearance of a carrion resource in southern Texas.

RESULTS

Our analysis of simulated data indicated reasonably precise estimates of joint space use can be achieved with commonly used GPS sampling rates (s.e.=0.029 at 150 locations per interval). Our analysis of field data identified an increase in spatial interactions between the coyotes that persisted for the duration of the study, beyond the expected duration of the carrion resource. Our analysis also identified a period of increased spatial interactions before appearance of the resource, which would not have been identified by previous methods.

CONCLUSIONS

We present a new approach to the analysis of joint space use through time, building upon tools commonly used by ecologists, that permits a new level of detail in the analysis of animal interactions. The results are easily interpretable and account for the nuances of bounded serial data in an elegant way.

Taxonomic and conservation implications of population genetic admixture, mito-nuclear discordance, and male-biased dispersal of a large endangered snake, Drymarchon couperi

Accurate species delimitation and description are necessary to guide effective conservation of imperiled species, and this synergy is maximized when multiple data sources are used to delimit species. We illustrate this point by examining Drymarchon couperi (Eastern Indigo Snake), a large, federally-protected species in North America that was recently divided into two species based on gene sequence data from three loci and heuristic morphological assessment. Here, we re-evaluate the two-species hypothesis for D. couperi by evaluating both population genetic and gene sequence data. Our analyses of 14 microsatellite markers revealed 6–8 genetic population clusters with significant admixture, particularly across the contact zone between the two hypothesized species. Phylogenetic analyses of gene sequence data with maximum-likelihood methods suggested discordance between mitochondrial and nuclear markers and provided phylogenetic support for one species rather than two. For these reasons, we place Drymarchon kolpobasileus into synonymy with D. couperi. We suggest inconsistent patterns between mitochondrial and nuclear DNA are driven by high dispersal of males relative to females. We advocate for species delimitation exercises that evaluate admixture and gene flow in addition to phylogenetic analyses, particularly when the latter reveal monophyletic lineages. This is particularly important for taxa, such as squamates, that exhibit strong sex-biased dispersal. Problems associated with over-delimitation of species richness can become particularly acute for threatened and endangered species, because of high costs to conservation when taxonomy demands protection of more individual species than are supported by accumulating data.

Space fit for a king: spatial ecology of king cobras (Ophiophagus hannah) in Sakaerat Biosphere Reserve, Northeastern Thailand

A species’ spatial ecology has direct implications for that species’ conservation. Far-ranging species may be more difficult to conserve because their movements increase their chances of encountering humans. The movements can take them out of protected areas, which is especially risky for species that are routinely persecuted. The king cobra (Ophiophagus hannah), a large venomous elapid, is subject to anthropogenic pressures, such as persecution and habitat loss. Here we present results from a study using radio telemetry to quantify movements and habitat use of nine king cobras in and around a protected area in Northeast Thailand. This study is the first investigation into the movements and habitat use of king cobras outside of the Western Ghats, India. On average, the tracked king cobra’s use areas of 493.42 ± 335.60 ha (95% fixed kernel), moving 183.24 ± 82.63 m per day. King cobras did not remain in intact forested area. Five of the individuals frequently used the human-dominated agricultural areas surrounding the protected area, appearing to make regular use of irrigation canals. Two adult males showed increases in movements during the breeding season. One male’s increased breeding season range caused him to venture beyond the protected area, shifting his habitat use from intact forests to scrub in human-dominated areas. King cobras’ large home range and willingness to use anthropogenic landscapes merits special consideration from conservation planners.