Spatial capture-recapture design and modelling for the study of small mammals

Analyzing the spatiotemporal habitat colonization dynamics of the common vole (Microtus arvalis Pallas) in Castilla y León, northwest Spain, using a hotspot-based approach†

BACKGROUND

The common vole is an agricultural pest and population outbreaks have caused significant crop damage in the agricultural areas of the west-central part of the Iberian Peninsula (Spain) during the last few decades. Thus, monitoring is imperative to gain a comprehensive insight on its spatiotemporal habitat colonization dynamics. This work was performed on a long-term database with the primary objective of establishing an operational framework for understanding the spatial and temporal distribution of the common vole, all of it based on the Getis-Ord statistics.

RESULTS

The temporal evolution of the vole abundance index (VAI) was consistent for the three studied habitats, that is, crops, reservoirs and paths. Furthermore, the majority of common vole abundance peaks coincide with summer periods, especially in annual herbaceous crops. The spatial distribution of vole abundance exhibited a non-random pattern, characterized by spatial clustering. Particularly, the areas with higher significance of this clustering were located at the so-called 'Tierra de Campos' county, covering northern Valladolid, southern Palencia and north-eastern Zamora provinces. Periods of major incidence and concentration were 2013-2014, 2016-2017 and 2019. Common vole temporal occupancy patterns demonstrate that colonization occurred simultaneously in alfalfa reservoirs and dispersion paths, both preceding settlement in annual herbaceous crop plots.

CONCLUSION

The geographic information system (GIS)-hotspots-based methodology proposed here can be valuable for stakeholders involved in integrated management of the common vole, serving as a detector of pest-prone areas in both space and time. These hotspots are useful for predicting future surveillance areas that accurately reflect pest colonization patterns. It was found that common vole abundance along dispersal paths acts as a source of dissemination, preceding the risk of colonization in annual herbaceous crop plots. © 2024 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Molecular detection of Cryptosporidium parvum in wild rodents (Phyllotis darwini) inhabiting protected and rural transitional areas in north-central Chile

Wild rodents often harbor Cryptosporidium species that can be transmitted to multiple mammal hosts. In Chile, little is known about Cryptosporidium in wild rodents, and available studies have been focused on morphological findings with no molecular-based evidence. A longitudinal survey was conducted between 2021 and 2022 to investigate the occurrence of Cryptosporidium spp. in populations of the Darwin's leaf-eared mouse (Phyllotis darwini) living in protected and rural transitional areas in north-central Chile, using staining and molecular methods. A total of 247 fecal samples were collected and examined by the modified Ziehl-Neelsen (ZN) staining test, 54 of which were positive for Cryptosporidium-like oocysts. Molecular analyses were carried out by PCR of the partial 18S ribosomal RNA and 60 kDa glycoprotein (gp60) genes. Cryptosporidium infection was confirmed in 34 samples (13.7 %) based on the PCR amplification, and individual (i.e., sex, and body mass index) and ecological variables (i.e., type of site and season) were not statistically significant (p > 0.05). Using the nucleotide sequencing of the partial 18S rRNA gene, Cryptosporidium parvum was identified in nine isolates. Also, C. parvum subgenotype family IIa was determined in seven samples by the partial gp60 gene, including the subtype IIaA17G4R1 in two samples. This is the first molecular evidence of Cryptosporidium parvum IIa in Phyllotis darwini in Chile. These results indicate potential cross-species transmition between wild rodents and domestic-wild animals in north-central Chile. More research is needed to understand better the role of wild rodents in the transmission of Cryptosporidium spp. in Chile.

Exploring patterns of female house mouse spatial organisation among outbreaking and stable populations

The size and distribution of home ranges reflect how individuals within a population use, defend, and share space and resources, and may thus be an important predictor of population-level dynamics. Eruptive species, such as the house mouse in Australian grain-growing regions, are an ideal species in which to investigate variations in space use and home range overlap between stable and outbreaking populations. In this study, we use spatially explicit capture-recapture models to explore if space use and home range overlap among female mice could serve as indicators of changes in population density leading into summer. Additionally, we assess the sensitivity of space use and home range estimates to reduced recapture rates. Our analysis did not reveal variations in the spring spatial organisation of female mice based on existing capture-mark-recapture data. However, our study highlights the need to balance monitoring efforts within regions, emphasising the importance of exploring studies that can improve spatial recaptures by optimising trapping efforts. This is particularly important in Australian agricultural systems, where varying farm management practices may drive differences in population dynamics.

Is bedding material a more effective thermal insulator than trap cover for small mammal trapping? A field experiment

Live trapping is a key technique for conducting ecological studies on small mammals. All-metal live traps are popular in monitoring schemes owing to their tested performance, lightweight design (aluminium) and foldability. However, capture represents a stressful situation for small mammals, particularly during cold seasons, when individuals are susceptible to cold weather starvation resulting from low temperature and insufficient food to maintain body temperature. Metal live traps provide limited protection against cold temperatures, and it is often recommended to use covers to buffer external temperature fluctuations and prevent entry of moisture. Here, we compared the insulative performance of a PVC cover designed for Sherman traps and of bedding material, using data loggers to record temperature and humidity inside traps. We conducted different experiments simulating field conditions (traps at night with a heat source inside) and different treatments (cover, bedding material) to test the thermal insulation capacity of three models of widely used commercial traps: Longworth, Sherman, and Heslinga. Our findings indicated that Longworth and Sherman traps were better insulated against ambient air temperature fluctuations than Heslinga traps (+2.0 °C warmer on average). Bedding material was paramount in reducing relative humidity and increasing thermal insulation capacity of traps (+3.1 °C), an effect that was strengthened when a PVC cover was additionally used (+4.2 °C). The covered traps prevented the direct entrance of rain and dew (reducing damp bedding), provided camouflage (reducing thefts), and improved thermal and humidity conditions of traps (potentially increasing survival of captive small mammals). Our results suggest that using covers and bedding materials can improve thermal and humidity conditions within live traps, thus reducing the metabolic costs of thermoregulation and increasing survival chances for trapped small mammals during cold seasons.

The first comprehensive population size estimations for the highly endangered largest diving beetle Dytiscus latissimus in Europe

Dytiscus latissimus (Coleoptera Dytiscidae) is an endangered diving beetle throughout its range. It is one of the two species of Dytiscidae listed in Annex II of the Habitats Directive, IUCN red list and in many national level legislations and therefore strictly protected. The conservation of endangered species first of all requires an assessment of their population size. Until now, a method has not been developed for estimating the size of D. latissimus populations. The article summarizes the results of two studies carried out independently in Germany and Latvia. Both studies were carried out in one water body used recapture method but with a different spatial placement of traps, which, according to our data, is an important factor in population estimation. We evaluated Jolly-Seber and Schnabel approaches of estimating aquatic beetle's populations and found that confidence intervals obtained by different methods in our research do not differ significantly, but combination of both models provide the most accurate estimates of population dynamics. As part of the study, we concluded that the populations of Dytiscus latissimus are relatively closed, so we accept that the Schnabel estimate shows more accurate data. By fixing the places of capture of each individual, it was found that females live mainly locally, and males actively move within the water body. This aspect indicates the advantage of the spatial placement of traps compared to the use of transects. The results of our study show a significantly higher number of both captured and recaptured males Such a sex ratio may indicate both a greater activity of males and differences in the sex ratio in the population. The study confirmed that environmental changes, such as the water level in a water body, can also significantly affect the result of a population assessment. In the frame of D. latissimus monitoring, to obtain an objective estimation of the species population size we recommend using four traps for each 100 m of water body shoreline with 4-8 censuses, dependently on the recapture rate.

Assessing small-mammal trapping design using spatially explicit capture recapture (SECR) modeling on long-term monitoring data

Few studies have evaluated the optimal sampling design for tracking small mammal population trends, especially for rare or difficult to detect species. Spatially explicit capture-recapture (SECR) models present an advancement over non-spatial models by accounting for individual movement when estimating density. The salt marsh harvest mouse (SMHM; Reithrodontomys raviventris) is a federal and California state listed endangered species endemic to the San Francisco Bay-Delta estuary, California, USA; where a population in a subembayment has been continually monitored over an 18-year period using mark-recapture methods. We analyzed capture data within a SECR modeling framework that allowed us to account for differences in detection and movement between sexes. We compared the full dataset to subsampling scenarios to evaluate how the grid size (area) of the trap design, trap density (spacing), and number of consecutive trapping occasions (duration) influenced density estimates. To validate the subsampling methods, we ran Monte Carlo simulations based on the true parameter estimates for each specific year. We found that reducing the area of the trapping design by more than 36% resulted in the inability of the SECR model to replicate density estimates within the SE of the original density estimates. However, when trapping occasions were reduced from 4 to 3-nights the density estimates were indistinguishable from the full dataset. Furthermore, reducing trap density by 50% also resulted in density estimates comparable to the full dataset and was a substantially better model than reducing the trap area by 50%. Overall, our results indicated that moderate reductions in the number of trapping occasions or trap density could yield similar density estimates when using a SECR approach. This approach allows the optimization of field trapping efforts and designs by reducing field efforts while maintaining the same population estimate compared to the full dataset. Using a SECR approach may help other wildlife programs identify sampling efficiencies without sacrificing data integrity for long term monitoring of population densities.

Biodiversity Monitoring of a Riparian Wetland in a Mixed-Use Watershed in the Central Appalachians, USA, before Restoration

Wetland mitigation efforts have increased in numbers over the past two decades to combat wetland loss in the United States. Data regarding wetland function such as biodiversity are required to be collected 5–10 years after a project is complete; however, pre-restoration data that can inform the effectiveness of mitigation are often not collected. We conducted pre-restoration surveys on various taxa along or within Ruby Run, a tributary of Deckers Creek in north-central West Virginia, USA, from 2016 to 2020 to determine the baseline relative abundance and diversity within the stream and the associated riparian zone. In five years, we observed 237 species (154 plant, 58 bird, 13 fish, 6 small mammal, and 6 anuran) and 25 families of macroinvertebrates. Seasonal fluctuations in diversity were present, but mean diversity was relatively consistent among years across taxa, except in anurans, where there was a decrease each year. Wetland mitigation efforts should continue to be monitored for success using multiple taxa, because land use change can affect taxa in different ways, resulting in well-rounded assessments that can improve wetland management practices.

Crested porcupine (Hystrix cristata) abundance estimation using Bayesian methods: first data from a highly agricultural environment in central Italy

Wildlife abundance estimation is one of the key components in conservation biology. Bayesian frameworks are widely used to adjust the potential biases derived by data collected in the field, as they can increase the precision of model parameter as a consequence of the combination of previous pieces of knowledge (priors) combined with data collected in the field to produce an a-posteriori distribution. Capture-recapture is one of the most common techniques used to assess animal abundance. However, the implementation with camera traps requires that animals present unique phenotypic traits for individual-based recognition. The crested porcupine Hystrix cristata is a semi-fossorial rodent with a continuous, but patchily distribution across Italy. Despite the species does not present evident individual-specific phenotypic traits, the information gathered using presence-only data obtained from camera traps, opportunistic observations, and road-killing events could be used to provide a rough estimate of the species abundance within an area. The main purpose of the present research was hence to provide the first preliminary estimate of the abundance of the crested porcupine in central Italy using presence-only data obtained from the above different monitoring methods. The results obtained estimated an average minimum number of 1803 individuals (SD = 26.89, CI 95% = 1750–1855) within an area covering about 17,111 km2. Since the porcupine is considered as “potentially problematic” because of damages to croplands and riverbanks, assessing its abundance is even more important to delineate adequate conservation and management actions to limit the potential trade-off effects over human activities.

A review of spatial capture-recapture: Ecological insights, limitations, and prospects

First described by Efford (2004), spatial capture-recapture (SCR) has become a popular tool in ecology. Like traditional capture-recapture, SCR methods account for imperfect detection when estimating ecological parameters. In addition, SCR methods use the information inherent in the spatial configuration of individual detections, thereby allowing spatially explicit estimation of population parameters, such as abundance, survival, and recruitment. Paired with advances in noninvasive survey methods, SCR has been applied to a wide range of species across different habitats, allowing for population- and landscape-level inferences with direct consequences for conservation and management. I conduct a literature review of SCR studies published since the first description of the method and provide an overview of their scope in terms of the ecological questions answered with this tool, taxonomic groups targeted, geography, spatio-temporal extent of analyses, and data collection methods. In addition, I review approaches for analytical implementation and provide an overview of parameters targeted by SCR studies and conclude with current limitations and future directions in SCR methods.

Linking Zoonosis Emergence to Farmland Invasion by Fluctuating Herbivores: Common Vole Populations and Tularemia Outbreaks in NW Spain

The expansion and intensification of agriculture are driving profound changes in ecosystems worldwide, favoring the (re)emergence of many human infectious diseases. Muroid rodents are a key host group for zoonotic infectious pathogens and frequently invade farming environments, promoting disease transmission and spillover. Understanding the role that fluctuating populations of farm dwelling rodents play in the epidemiology of zoonotic diseases is paramount to improve prevention schemes. Here, we review a decade of research on the colonization of farming environments in NW Spain by common voles (Microtus arvalis) and its public health impacts, specifically periodic tularemia outbreaks in humans. The spread of this colonizing rodent was analogous to an invasion process and was putatively triggered by the transformation and irrigation of agricultural habitats that created a novel terrestrial-aquatic interface. This irruptive rodent host is an effective amplifier for the Francisella tularensis bacterium during population outbreaks, and human tularemia episodes are tightly linked in time and space to periodic (cyclic) variations in vole abundance. Beyond the information accumulated to date, several key knowledge gaps about this pathogen-rodent epidemiological link remain unaddressed, namely (i) did colonizing vole introduce or amplified pre-existing F. tularensis? (ii) which features of the “Francisella—Microtus” relationship are crucial for the epidemiology of tularemia? (iii) how virulent and persistent F. tularensis infection is for voles under natural conditions? and (iv) where does the bacterium persist during inter-epizootics? Future research should focus on more integrated, community-based approaches in order to understand the details and dynamics of disease circulation in ecosystems colonized by highly fluctuating hosts.

Patterns of flea infestation in rodents and insectivores from intensified agro-ecosystems, Northwest Spain

BACKGROUND

Fleas frequently infest small mammals and play important vectoring roles in the epidemiology of (re)emerging zoonotic disease. Rodent outbreaks in intensified agro-ecosystems of North-West Spain have been recently linked to periodic zoonotic diseases spillover to local human populations. Obtaining qualitative and quantitative information about the composition and structure of the whole flea and small mammal host coexisting communities is paramount to understand disease transmission cycles and to elucidate the disease-vectoring role of flea species. The aims of this research were to: (i) characterise and quantify the flea community parasiting a small mammal guild in intensive farmlands in North-West Spain; (ii) determine and evaluate patterns of co-infection and the variables that may influence parasitological parameters.

METHODS

We conducted a large-scale survey stratified by season and habitat of fleas parasitizing the small mammal host guild. We report on the prevalence, mean intensity, and mean abundance of flea species parasitizing Microtus arvalis, Apodemus sylvaticus, Mus spretus and Crocidura russula. We also report on aggregation patterns (variance-to-mean ratio and discrepancy index) and co-infection of hosts by different flea species (Fager index) and used generalized linear mixed models to study flea parameter variation according to season, habitat and host sex.

RESULTS

Three flea species dominated the system: Ctenophthalmus apertus gilcolladoi, Leptopsylla taschenbergi and Nosopsyllus fasciatus. Results showed a high aggregation pattern of fleas in all hosts. All host species in the guild shared C. a. gilcolladoi and N. fasciatus, but L. taschenbergi mainly parasitized mice (M. spretus and A. sylvaticus). We found significant male-biased infestation patterns in mice, seasonal variations in flea abundances for all rodent hosts (M. arvalis, M. spretus and A. sylvaticus), and relatively lower infestation values for voles inhabiting alfalfas. Simultaneous co-infections occurred in a third of all hosts, and N. fasciatus was the most common flea co-infecting small mammal hosts.

CONCLUSIONS

The generalist N. fasciatus and C. a. gilcolladoi dominated the flea community, and a high percentage of co-infections with both species occurred within the small mammal guild. Nosopsyllus fasciatus may show higher competence of inter-specific transmission, and future research should unravel its role in the circulation of rodent-borne zoonoses.

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.

Mark my words: experts’ choice of marking methods used in capture-mark-recapture studies of small mammals

Crucial to the success of studies based on capture-mark-recapture (CMR) designs is the retention (permanency) and recognition (readability) of marks to identify individuals. Several marking methods for small mammals (< 60 g) are available, but their efficacy and use is not well known. We implemented a targeted survey of experts to gather their experiences and opinions regarding marking small mammals. Respondents (n = 114) stated their beliefs, perceptions, and current and future use, of marking methods, as well as factors influencing their choices, based on Likert and rank order scale questions. We compared responses based on where researchers’ studies occurred, their level of experience, and their subfield of mammalogy. Most respondents (73%) had > 5 years experience marking small mammals, with 60% each marking > 1,000 individuals. Respondents believed that ear-tagging was most preferable in terms of efficiency, impact to affected animals (survival, pain), and personal ethics, whereas passive integrated transponder (PIT)-tagging was the most preferable with regard to retention and recognition, and toe-clipping with respect to cost. Most respondents plan to use ear-tagging (78%) or PIT-tagging (70%) in the future. PIT- and genetic-tagging are expected to increase, and toe-clipping to decline, in the future. The factors influencing which marking method respondents used were ranked—in order of decreasing preference—as impact, retention, recognition, cost, efficiency, and ethics. There were few differences in the mean response or consensus among respondents, regardless of their experience, location, or subfield. Most respondents (66%) agreed that additional studies on the performance and impact of various marking methods are needed to assess their costs and benefits for CMR-based studies. Ultimately, choice of marking method will depend on the species, research question, available resources, and local legislation and permitting. Our study, however, illustrates that collective insights by experienced mammalogists may aid individual researchers in deciding on study designs and protocols, particularly early career scientists.

Numerical response of a mammalian specialist predator to multiple prey dynamics in Mediterranean farmlands

The study of rodent population cycles has greatly contributed, both theoretically and empirically, to our understanding of the circumstances under which predator-prey interactions destabilize populations. According to the specialist predator hypothesis, reciprocal interactions between voles and small predators that specialize on voles, such as weasels, can cause multiannual cycles. A fundamental feature of classical weasel-vole models is a long time-lag in the numerical response of the predator to variations in prey abundance: weasel abundance increases with that of voles and peaks approximately 1 yr later. We investigated the numerical response of the common weasel (Mustela nivalis) to fluctuating abundances of common voles (Microtus arvalis) in recently colonized agrosteppes of Castilla-y-Léon, northwestern Spain, at the southern limit of the species' range. Populations of both weasels and voles exhibited multiannual cycles with a 3-yr period. Weasels responded quickly and numerically to changes in common-vole abundance, with a time lag between prey and weasel abundance that did not exceed 4 months and occurred during the breeding season, reflecting the quick conversion of prey into predator offspring and/or immigration to sites with high vole populations. We found no evidence of a sustained, high weasel abundance following vole abundance peaks. Weasel population growth rates showed spatial synchrony across study sites approximately 60 km apart. Weasel dynamics were more synchronized with that of common voles than with other prey species (mice or shrews). However, asynchrony within, as well as among sites, in the abundance of voles and alternative prey suggests that weasel mobility could allow them to avoid starvation during low-vole phases, precluding the emergence of prolonged time lag in the numerical response to voles. Our observations are inconsistent with the specialist predator hypothesis as currently formulated, and suggest that weasels might follow rather than cause the vole cycles in northwestern Spain. The reliance of a specialized predator on a functional group of prey such as small rodents does not necessarily lead to a long delay in the numerical response by the predator, depending on the spatial and interspecific synchrony in prey dynamics.