Dens are places for cavernicolous animals to hibernate, reproduce, and avoid predators and harsh weather conditions, and thus they have a vital impact on their survival. M. himalayana is the main large cavernicolous rodent on the Qinghai-Tibet Plateau. The analysis of den traits and their ecological functions can reveal mechanisms by which marmots have adapted to their environment, which is important for further understanding the ecological significance of this species. From July to August 2019 (warm season), we used unmanned aerial vehicles to fly at low altitudes and slow speeds to locate 131 marmot burrows (45 on shaded slopes, 51 on sunny slopes, and 35 on flat areas) in the northeastern Qinghai-Tibet Plateau. We then measured the physical characteristics (den density, entrance size, first tunnel length, volume, orientation and plant characteristics near the den entrance) of these dens on site. The physical parameters of the M. himalayana dens showed that they function to protect the marmots from natural enemies and bad weather, provide good drainage, and maintain a stable microclimate around the entrance. This is a result of the marmot’s adaptation to the harsh environment (cold and humidity) of the Qinghai-Tibet Plateau.

1. The description and analysis of animal behaviour over long periods of time is one of the most important challenges in ecology. However, most of these studies are limited due to the time and cost required by human observers. The collection of data via video recordings allows observation periods to be extended. However, their evaluation by human observers is very time-consuming. Progress in automated evaluation, using suitable deep learning methods, seems to be a forwardlooking approach to analyse even large amounts of video data in an adequate time frame. 2. In this study we present amulti-step convolutional neural network system for detecting animal behaviour states, which works with high accuracy. An important aspect of our approach is the introduction of model averaging and post-processing rules to make the system robust to outliers. 3. Our trained system achieves an in-domain classification accuracy of >0.92, which is improved to >0.96 by a postprocessing step. In addition, the whole system performs even well in an out-of-domain classification task with two unknown types, achieving an average accuracy of 0.93. We provide our system at https://github.com/Klimroth/Video-Action-Classifier-for-African-Ungulates-in-Zoos/tree/main/mrcnn_based so that interested users can train their own models to classify images and conduct behavioural studies of wildlife. 4. The use of a multi-step convolutional neural network for fast and accurate classification of wildlife behaviour facilitates the evaluation of large amounts of image data in ecological studies and reduces the effort of manual analysis of images to a high degree. Our system also shows that post-processing rules are a suitable way to make species-specific adjustments and substantially increase the accuracy of the description of single behavioural phases (number, duration). The results in the out-of-domain classification strongly suggest that our system is robust and achieves a high degree of accuracy even for new species, so that other settings (e.g. field studies) can be considered.

Artiodactyl prey species of Chile, especially guanacos (Lama guanicoe) are reported to be very susceptible to predation by pack hunting feral dogs. It has been previously suggested that guanacos and endemic South American deer may have evolved in the absence of pack-hunting cursorial predators. However, the paleoecology of canid presence in southern South America and Chile is unclear. Here, we review the literature on South American and Chilean canids, their distributions, ecologies and hunting behaviour. We consider both wild and domestic canids, including Canis familiaris breeds. We establish two known antipredator defense behaviours of guanacos: predator inspection of ambush predators, e.g. Puma concolor, and rushing at and kicking smaller cursorial predators, e.g. Lycalopex culpaeus. We propose that since the late Pleistocene extinction of hypercarnivorous group-hunting canids east of the Andes, there were no native species creating group-hunting predation pressures on guanacos. Endemic deer of Chile may have never experienced group hunting selection pressure from native predators. Even hunting dogs (or other canids) used by indigenous groups in the far north and extreme south of Chile (and presumably the center as well) appear to have been used primarily within ambush hunting strategies. This may account for the susceptibility of guanacos and other prey species to feral dog attacks. We detail seven separate hypotheses that require further investigation in order to assess how best to respond to the threat posed by feral dogs to the conservation of native deer and camelids in Chile and other parts of South America.

Diffusible iodine-based contrast-enhanced Computed-Tomography (diceCT) visualizes soft-tissue from microCT (µCT) scans of specimens to uncover internal features and natural history information without incurring physical damage via dissection. Unlike hard-tissue imaging, diceCT datasets are currently limited to a few individual specimens and taxonomically underrepresented. To initiate best practices for diceCT in a non-model group, we outline a guide for staining and high-throughput µCT scanning in snakes. We scanned the entire body and one region of interest (i.e., head) for 23 specimens representing 23 species from the clades Aniliidae, Dipsadinae, Colubrinae, Elapidae, Lamprophiidae and Viperidae. We generated 82 scans that include 1.25% Lugols iodine stained (soft tissue) and unstained (skeletal) data for each specimen. We found that duration of optimal staining time increased linearly with body size; head radius was the best indicator. Post-reconstruction of scans, optimal staining was evident by evenly distributed grayscale values and clear differentiation among soft-tissue anatomy. Under and over stained specimens produced poor contrast among soft-tissues, which was often exacerbated by user bias during “digital dissections” (i.e., segmentation). Regardless, all scans produced usable data from which we assessed a range of downstream analytical applications within ecology and evolution (e.g., predator-prey interactions, life history, and morphological evolution). Ethanol de-staining reversed the known effects of iodine on the exterior appearance of physical specimens, but required substantially more time than reported for other de-staining methods. We discuss the feasibility of implementing diceCT techniques for a new user, including approximate financial and temporal commitments, required facilities, and potential effects of staining on specimens. We present the first high-throughput workflow for full-body skeletal and diceCT scanning in snakes, which can be generalized to any elongate vertebrates, and increases publicly available diceCT scans for reptiles by an order of magnitude.

Identifying factors that influence sea otter (Enhydra lutris) population density can provide insight into why it varies spatially and temporally and when a recovering population has reached an equilibrium density because of food resources (i.e., carrying capacity K). Although food availability is widely recognized as an important extrinsic factor affecting sea otter density, how do we determine when a population has reached K? The goal of this study was to estimate K for Simpson Bay, Alaska by measuring the abundance of edible bivalves, the primary prey for sea otters for over 40 years. We then compared prey abundance and estimated replacement rate (i.e., the mean age of bivalves predated by sea otters) to estimated annual prey consumption based on the mean population density for the past 18 years. On average, 110 adult sea otters (5.2 km−2) have occupied Simpson Bay annually since 2001 consuming an estimated 176,660 kg of bivalves. The total mass (standing stock) of the major bivalves (predominately butter clams and stained macomas) was 785,730 kg, so adult sea otters consumed about 22% annually. Based on these observations and calculations, the estimated annual number of sea otters occupying Simpson Bay appears to be at or near K based on the replacement rate of food resources. However, other intrinsic (e.g., male territoriality and emigration) and extrinsic (e.g., predation, disease, human-related mortality) factors may influence equilibrium density, which varies spatially and temporally, resulting in a mosaic of subpopulations with different densities, rates of growth and discontinuous distributions. Understanding the balance among these factors may be one of the most challenging ecological questions for sea otter conservation and management as populations recover from their range-wide decimation during the Maritime Fur Trade in the late 18th and 19th centuries

The muskrat (Ondatra zibethicus) is an iconic species in Canada, valued for both its fur and its integral role in wetland ecosystems, and widely regarded for its perseverance. However, the resilience of this semi-aquatic mammal seems to be in question now as increasing evidence points to widespread population declines. Recent analyses of harvest data across North America suggest a reduction in their numbers, but this has not been widely corroborated by population surveys. In this study we replicated historic muskrat house count surveys at two large Great Lakes coastal wetlands and present confirmation that declines in muskrat harvest correspond to actual declines in muskrat abundance. At the Point Pelee National Park marsh and the Matchedash Bay-Gray Marsh wetland we found that mean muskrat house counts declined by 93% and 91% respectively between historic surveys 40-50 years ago and contemporary surveys over the past five years. The factors responsible for these dramatic declines remain unclear but there may be a relationship with changes in the habitat quality of these wetlands that have occurred over the same time frame. Not only is the loss of muskrats an issue for the resulting loss of the wetland ecosystem services they provide, but it may be an indication of broader marsh ecosystem degradation. As such, a scarcity of muskrats should be considered a red flag for the state of biodiversity in our wetlands. Continued surveys and ongoing research are needed to shed more light on the current status of muskrat populations and their marsh habitats across their native range. Keywords: Fur harvest; Muskrat; Ondatra; Population decline; Typha; Wetlands

Predation is the most common cause of nest failure in birds. While nest predation is relatively well studied in general, our knowledge is unevenly distributed across globe and taxa, with for example limited information on shorebirds breeding in sub-tropics. Importantly, we know fairly little about the timing of predation within a day and season. Here, we followed 499 nests of red-wattled lapwings (Vanellus indicus), a ground-nesting shorebird, to estimate a nest predation rate, and continuously monitored 231 of these nests for a sum of 2951 days to reveal how timing of predation changes over the day and season in a sub-tropical desert. We found that 324 nests hatched, 77 nests were predated, 38 failed for other reasons and 60 had unknown fate. Daily predation rate was 0.97% (95%CrI: 0.77% – 1.2%), which for a 30-day long incubation period translates into ~25% chance of nest being predated. Such predation rate is low compared to most other species. Predation events were distributed evenly across day and night, with a tendency for increased predation around sunrise. Predation rate and events were distributed evenly also across the season, although night predation was more common later in the season, perhaps because predators reduce their activity during daylight to avoid extreme heat. Indeed, nests were never predated upon when mid-day ground temperatures exceeded 45°C. Whether the activity pattern of predators indeed changes across the breeding season and whether the described predation patterns hold for other populations, species and geographical regions awaits future investigations.

Intralocus conflict has been well documented between the sexes, but much less is known about the potential for this genetic conflict in other polymorphisms, such as alternative reproductive tactics. Here we investigate two of three criteria necessary for demonstrating intralocus tactical conflict: the ARTs have different phenotypic optima for a shared trait, and one or both ARTs are not at their phenotypic optima for a shared trait. We address these two criteria in Xiphophorus multilineatus, a live-bearing freshwater fish that has two male alternative reproductive tactics, a behaviorally fixed courter male and a behaviorally plastic sneaker male that switches between courtship and force-copulatory behavior. We used measures of reproductive success from a wild population to estimate selection gradients on three tactically dimorphic traits involved in sexual selection: body size, body shape, and sword length. We present evidence that both body size and sword length are experiencing tactically antagonistic selection, providing evidence for both criteria. Additionally, selection on body shape in sneaker males appears to be buffered due to behavioral plasticity. Our study provides novel insight from a wild population into the role that intralocus tactical conflict can play in constraining ARTs from reaching their respective phenotypic optima despite tactical dimorphism.

At least nine felid species can co-occur in Southeast Asia, thus providing an unusual opportunity to investigate poorly known guild structure and the factors controlling it. Using camera-trap data, we quantified space use, temporal activity, and multi-dimensional niche overlap of tiger, clouded leopard, marbled cat, leopard cat, and Asiatic golden cat in Htamanthi Wildlife Sanctuary of Myanmar. We hypothesized that the spatio-temporal behaviour of smaller cats can reflect avoidance of the larger cats, which are both potential competitors and predators, and similar-sized guild members would partition their niches in space or time to reduce competition for resources. Our approach involved single-species occupancy modelling to identify site covariates, pairwise spatial overlap using Bayesian inference, and activity overlap with Kernel density estimation and multivariate analyses to test hypotheses. We found tiger and marbled cats were primarily diurnal, clouded leopard and leopard cat were nocturnal and golden cat exhibited cathemeral activity. We observed a complex pattern of guild assembly and potential competition involving strong niche displacement between the golden cat and marbled cat, but high overlap between the relatively similarly-sized pairing of clouded leopard and golden cat, and the markedly differently-sized tiger -- golden cat pairing. No significant evidence of mesopredator release was observed and the felid assembly in Northern Myanmar appeared to be partitioned mainly on a spatial, rather than temporal, dimension. Nonetheless, the temporal association between the three mesopredators was inversely related to the similarity in their body sizes. The insights into this felid guild revealed that the largest niche differences in the use of space and time occurred between the three smaller species, most evidently between the Asiatic golden cat-marbled cat pairing, followed by marbled cat - leopard cat pairing. This study offers new insight into carnivore guild assembly and, adds substantially to knowledge of five of the least known felids of conservation concern.

The major histocompatibility complex (MHC) is a multiple-copy immune gene family in vertebrates. Its genes are highly variable and code for antigen-presenting molecules. Characterization of MHC genes in different species and investigating the mechanisms that shape MHC diversity is an important goal in understanding the evolution of biological diversity. Here we developed a next generation sequencing (NGS) protocol to genotype the MHC class I genes of 326 Godlewski’s buntings (Emberiza godlewskii) sampled in the Western mountain area of Beijing from 2014 to 2016. A total of 184 functional alleles were identified, including both non-classical and classical alleles. Classical alleles could be clustered into nine supertypes. Compared with other passerine birds, the individual diversity of MHC class I genes in Godlewski’s buntings is intermediate. Ten amino acid sites in the antigen-binding domain showed signatures of positive selection and eight of them exhibit high amino acid polymorphism. These findings indicate the action of balancing selection and provide a framework for subsequent investigation of selection acting on MHC genes in Godlewski’s buntings.

Hybridization is a common and important stage in species formation in plants and animals. The evolutionary consequences of hybridization depend not only on reproductive compatibility between sympatric species, but also on factors like vulnerability to each other’s predators and parasites. We examine infection patterns of the blood parasite Haemoproteus lophortyx, a causative agent of avian malaria, at a site in the contact zone between California quail (Callipepla californica) and Gambel’s quail (C. gambelii). We tested whether species identity, sex, and year predicted infection status and intensity. While we found no effect of sex on the status or intensity of infection, we found differences in infection status and intensity across species and between years. The prevalence of infection in California and hybrid quail was lower than in Gambel’s quail. Once infected, however, California and hybrid quail had higher infection intensities than Gambel’s quail. California and hybrid quail exhibited no significant differences in prevalence or intensity of infection. These findings suggest that infection by H. lophortyx has the potential to influence species barrier dynamics in this system, however, more work is necessary to determine the exact evolutionary consequences of this blood parasite.

Understanding what variables affect ungulate neonate survival is imperative to successful conservation and management of the species. Predation is commonly cited as a cause-specific source of mortality and ecological covariates often influence neonate survival. However, variation in survival estimates related to capture methodology has been documented with opportunistically captured neonates generally displaying greater survival than those captured via aid of vaginal implant transmitters (VITs), likely because of increased left truncation observed in the opportunistically captured datasets. Our goal was to assess if 3- and 6-month survival estimates varied by capture method while simultaneously assessing if capture method affected model selection and interpretation of ecological covariates for white-tailed deer neonates captured from three study sites in North Dakota and South Dakota, USA. We found survival varied by capture method for 3-month neonate survival with opportunistically captured neonates displaying up to 26% greater survival than their counterparts captured via VITs; however, this relationship was not present for 6-month survival. We also found model selection and subsequent interpretation of ecological covariates varied when analyzing datasets comprised of neonates captured via VITs, neonates captured opportunistically, and all neonates combined regardless of capture method. When interpreting results from our VIT only analysis for 3-month survival, we found survival varied by three time intervals and was lowest in the first two weeks of life. Capture method did not affect 6-month survival which was most influenced by total precipitation occurring during 3 – 8 weeks of a neonate’s life and percent canopy cover found at a neonate’s capture site. Our results support previous research that capture method must be accounted for when deriving survival estimates for ungulate neonates as it can impact derived estimates and subsequent interpretation of results.

Mating systems and patterns in reproductive success of fishes play an important role in ecology and evolution. While information on the reproductive ecology of many anadromous salmonids (Oncorhynchus spp.) is well-detailed, there is less information for non-anadromous species including the Yellowstone Cutthroat Trout (O. clarkii bouvieri), a species of recreational angling importance and conservation concern. Here, we used data from a parentage-based tagging study to describe the mating system of Yellowstone Cutthroat Trout from a spawning tributary of the South Fork Snake River, Idaho, and identify predictors of relative reproductive success. We detected evidence of monogamy, polygyny, and polyandry and showed that reproductive success was best explained by arrival time at the spawning ground and total length. Specifically, the largest adults arrived earliest in the season and produced a disproportionate number of offspring. Lastly, we estimated the effective number of breeders (Nb) and effective population size (Ne) and showed that while Nb was lower than Ne, both are sufficiently high to suggest Yellowstone Cutthroat Trout in Burns Creek represent a genetically stable and diverse population.

International socioeconomic relationships form the background that underlies the history of invasive species. Species with economic value, such as the North American bullfrog (Lithobates catesbeianus), are more likely to become internationally distributed and to be subsequently introduced to non-native areas and, consequently, become more difficult to control. Using population genetics methods, we investigated the invasion pathways, the connectivity among clusters in different countries and the native population of origin of globally introduced bullfrog populations. Throughout the analysis of seven microsatellite nuclear loci, one fragment of the mitochondrial cytochrome b locus, and historical information, four main lineages were identified and analyzed along with previous findings. This species’ capability to colonize several countries from few starting lineages highlights the necessity to control new propagule pressure to ensure successful management programs, as high inbreeding and bottleneck effect seem not to diminish the invasive success of this species. There is a consensus between markers that most areas of South America belong to the same genetic population while populations in Asia have a more complex history of introduction.

Identifying phenotypic characteristics of evolutionarily fit individuals provides important insight into the evolutionary processes that cause range shifts with climate warming. Female beluga whales (Delphinapterus leucas) from the Canadian high Arctic (BB) residing in the core region of the species’ geographic range are 14% larger than their conspecifics at the southern periphery in Hudson Bay (HB). We investigated the causal mechanism for this north (core)-south (periphery) difference as it relates to fitness by combining morphometric data with ovarian corpora counted in female reproductive tracts. We found evidence for reproductive senescence in older HB females from the southern peripheral population but not for BB whales. Female beluga whale fitness in the more-northern BB increased faster with age (48% partial variation explained) versus a more gradual slope (25%) in HB. In contrast, body length in HB female beluga accounted for five times more of the total variation in fitness compared to BB whales. We speculate that female HB beluga fitness was more strongly linked with body length due to higher density, as larger body size provides survival advantages during seasonal food limitations. Understanding the evolutionary mechanism of how fitness changes will assist conservation efforts in anticipating and mitigating future challenges to peripheral populations.

We conduct a phylogeographic and population genetic study of the Asiatic Toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate of the region. A total of 292 individuals from 94 locations were genotyped for two mitochondrial DNA loci (cytb, ND2 gene) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, DBI-2). We performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: West (China mainland) and Northeast (Northeast China, Russia and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical—robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Cooperative breeding, which is commonly characterized by non-breeding individuals that assist others with reproduction, is common in avian species. However, few accounts have been reported in Charadriiformes, particularly island-nesting species. We present observations of cooperative breeding behaviors in Hawaiian Stilts during the 2012-2020 nesting seasons on the Hawaiian islands of O‘ahu and Moloka‘i. We describe three different behaviors that indicate cooperative breeding: (1) nest sharing; (2) helper at the nest; (3) cooperative chick rearing. Our observations suggest an ideal opportunity to examine the evolution of cooperative breeding behaviors in the order Charadriiformes.

The relative contributions of adaptation and drift to morphological diversification of the crania of echolocating mammals was investigated using two horseshoe bat species, Rhinolophus simulator and R. cf. simulator as test cases. We used 3D geometric morphometrics to compare the shapes of skulls of the two lineages collected at various localities in southern Africa. Shape variation was predominantly attributed to selective forces; the between population variance (B) was not proportional to the within population variance (W). Modularity was evident in the crania of R. simulator but absent in the crania of R. cf. simulator and the mandibles of both species. The skulls of the two lineages thus appeared to be under different selection pressures, despite the overlap in their distributions. Selection acted mainly on the nasal dome region of R. cf. simulator whereas selection acted more on the cranium and mandibles than on the nasal domes of R. simulator. Probably the relatively higher echolocation frequencies used by R. cf. simulator, the shape of the nasal dome, which acts as a frequency dependent acoustic horn, is more crucial than in R. simulator, allowing maximization of the intensity of the emitted calls and resulting in comparable detection distances. In contrast, selection pressure is probably more pronounced on the mandibles and cranium of R. simulator to compensate for the loss in bite force because of its elongated rostrum. The predominance of selection probably reflects the stringent association between environment and the optimal functioning of phenotypic characters associated with echolocation and feeding in bats.

The main strategy for animal diversity conservation is to increase the territory size but little consideration is given to habitat characteristics requirement, which lead to a decrease in effectiveness for protected areas. Marginal of protected areas are considered to have higher species richness due to the edge effect. Strategy in these sites are still adopts to increase territory size or pay no attention to needs of specific habitat characteristics that is an important topic for the planner and manager. In this study, camera traps was used to estimate composition, diversity and habitat characteristics of mammals in a non-protected area near Huangshan Mountains in Anhui Province, China. We ran 49 liner models with the relative abundance index and 13 habitat characteristic factors of 11 mammals. To answer the question of habitat characteristics or territory size: which is more important to composition and diversity of mammals in non-protect area? We hypothesized that: (1) Non-protected areas have more mammal species than protected areas with the edge effect. (2) Non-protected areas have more species associated with habitat characteristics. We predicted that the habitat characteristics should be firstly considered, territory size secondly in non-protected areas, would provide a last refuge for mammals. Cameras were operated from June 2017 to October 2019, for a total of 29 months, 2,212 independent photos, 9,485 trap-days, recorded 18 species of mammals more than any other protected areas confirmed first hypothesis 1. The model analysis results showed that, habitat characteristics of mammals were different and showed a significant correlation, supported hypothesis 2. In addition, most species are related to vegetation characteristics except to primates (Macaca. thibetana) and rodent (Leopoldamys edwardsi) confirmed our prediction. We suggested conservation policies in non-protected areas: Habitat characteristics should be concerned at first and then increasing protected areas to provide the last refuge for species conservation.

Increasingly animal behaviour studies are enhanced through the use of accelerometry. To allow translation of raw accelerometer data to animal behaviours requires the development of classifiers. Here, we present the “rabc” package to assist researchers with the interactive development of such animal-behaviour classifiers based on datasets consisting out of accelerometer data with their corresponding animal behaviours. Using an accelerometer and a corresponding behavioural dataset collected on white stork (Ciconia ciconia), we illustrate the workflow of this package, including raw data visualization, feature calculation, feature selection, feature visualization, extreme gradient boost model training, validation, and, finally, a demonstration of the behaviour classification results.

Stable isotope analysis is a universally recognized and efficient method of indicating trophic relationships that is widely applied in research. However, variations in natural isotopic abundance may lead to inaccuracies due to the effects of complex environmental conditions. This research compared the carbon and nitrogen isotopic niches of fish communities between diverse biotopes around the Yellow River estuary and adjacent sea areas, with the aim of revealing distinctions in stable isotopic niche metrics, trophic positions, and feeding preferences. Stable isotopic niche results indicated that the communities of estuarine habitants were compatible in most study biotopes, and may provide a corridor for energy and material transportation between Laizhou Bay and the open water. Local biocoenosis was embodied in the wider isotopic niche corresponding to frequent environmental changes and abiotic gradients. This implied that they used various food sources to adapt to the fickle environment, including marine-terrestrial boundaries and the estuary. Our analysis of the food source contribution indicated that allochthonous sources were considered major energy sources in estuarine areas directly affected by Yellow River-diluted water, while autochthonous benthic and pelagic producers dominated carbon input into the food web in Laizhou Bay and the open water. A significant variation in the fish δ15N characteristic was found within estuarine adjacent regions, so, together with the results from previous studies, we deemed the local high concentration of dissolved inorganic nitrogen as the original trigger of the abnormal δ15N characteristic in fishes via a transport process along food chains. These results provide a new perspective on the natural distinction of carbon and nitrogen isotopic niches. The detailed data reported here enhance our understanding of variations in fish communities in estuarine ecosystems.

1. Understanding the animal-habitat relationship at local scale is crucial in ecology, particularly to develop strategies for wildlife management and conservation. As this relationship is governed by environmental features and intra and inter-specific interactions, habitat selection of a population may vary locally between its core and edges. 2. This is particularly true for central place foragers, such as grey and harbour seals, whose trends in numbers vary among different regions in the Northeast Atlantic. Here, we aimed at studying how foraging habitat selection may vary locally with the influence of population trends and physical habitat features 3. Using GPS/GSM tags deployed in grey and harbour seal colonies of contrasting sizes, we investigate spatial patterns and foraging habitat selection by comparing trip characteristics and home range similarities, and fitting GAMM to the seal distribution and environmental data respectively. 4. We show that grey seal foraging habitat selection and spatial patterns differed markedly between regions. Grey seals may select environmental characteristics for their foraging habitat accounting for local differences in prey consumed. Spatial patterns were different might depend on local seal density and regional productivity, located from inshore to offshore areas for the limit ranges and core population respectively. Our results on foraging habitat selection reflected the coastal and sedentary behaviour of harbour seals. We found no difference in spatial patterns between colonies, except for the Inner Hebrides where seals foraged further, potentially reflecting density dependence pressure, as the number in this colony is higher. 5. These results suggest that local conditions might have a strong influence on population spatial ecology, highlighting as well the relevance of studying foraging habitat selection based on foraging behaviour at fine geographical scale, particularly if species are managed within regional units.

Aim: In recent studies, three species of bat, Hypsignathus monstrosus, Myonycteris torquata and Epomops franqueti were identified as the most likely candidates to be reservoir for Ebolavirus. To help in epidemic preparedness and surveillance, species distribution modeling techniques are useful for predicting where these species are likely to occur in DR Congo. Methods: MaxEnt software was used to model the current and future distribution of the three species in DR Congo based on the occurrence data collected from global biodiversity information facility and environmental covariates collected from Worldclim and USGS. The future distribution were obtained based on two scenario (RCP4.5 and RCP 8.5) of the HadGEM2-CC an IPPC5 climate projections from global climate models (GCMs). Results: The suitable habitat of H. monstrosus, M. torquata and E. franqueti are essentially located from latitude -5° to 5° and longitude 17° to 30°. Overall, their favourable living areas are located in the territories along the entire northern border and the entire northern part of the eastern border of the DR Congo with varying degrees of importance depending on the territories and species. Three bioclimatic variables have shown to play major role in their spread, the precipitation of the driest quarter, the precipitation of coldest quarter and the temperature annual range. In the future, the suitable area of these species will be decreasing and being essentially located in the Kivu provinces. The change in the ecological niche of these species will differ with respect to future climate scenario. Climate of RCP 8.5 has shown to induce major decrease of their suitable habitat in DR Congo. Main conclusions: The MaxEnt model is potentially useful for forecasting the future adaptive distribution of the three bat species under climate change, and it provides important guidance for comprehensive management of the Ebolavirus risk.

Buck Creek is a spring-fed, cool-water tributary of the West Fork White River, Indiana. The Muncie Bureau of Water Quality sampled fishes and monitored water temperature in Buck Creek annually from 1986-2018. For this study, we utilized long-term fish data from the Bureau of Water Quality to evaluate spatial and temporal changes in the fish assemblages of Buck Creek in Delaware County, Indiana, USA. Non-metric multidimensional scaling (NMDS) was used to describe changes in the fish assemblages over space and time. Linear mixed effects models were used to evaluate the relationship between environmental factors and the fish assemblages. The spatial NMDS results were separated in distinct groups of upstream and downstream assemblages. This was characterized by a shift of headwater specialists shifting to large-river species. The temporal NMDS results were separated into distinct annual assemblages. This was characterized by a drop in pollution-tolerant species and an increase in intolerant species. Our findings indicate that the fish assemblages have improved in Buck Creek over space and time.