USGS Science⌕ Search

SEARCH · USGS Science

Results for “Population Ecology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74Linked to original sources
USGS2015

Spatial structure of morphological and neutral genetic variation in Brook Trout

Brook Trout Salvelinus fontinalis exhibit exceptional levels of life history variation, remarkable genetic variability, and fine-scale population structure. In many cases, neighboring populations may be highly differentiated from one another to an extent that is comparable with species-level distinctions in other taxa. Although genetic samples have been collected from hundreds of populations and tens of thousands of individuals, little is known about whether differentiation at neutral markers reflects phenotypic differences among Brook Trout populations. We compared differentiation in morphology and neutral molecular markers among populations from four geographically proximate locations (all within 24 km) to examine how genetic diversity covaries with morphology. We found significant differences among and/or within streams for all three morphological axes examined and identified the source stream of many individuals based on morphology (52.3% classification efficiency). Although molecular and morphological differentiation among streams ranged considerably (mean pairwise F ST : 0.023–0.264; pairwise P ST : 0.000–0.339), the two measures were not significantly correlated. While in some cases morphological characters appear to have diverged to a greater extent than expected by neutral genetic drift, many traits were conserved to a greater extent than were neutral genetic markers. Thus, while Brook Trout exhibit fine-scale spatial patterns in both morphology and neutral genetic diversity, these types of biological variabilities are being structured by different ecological and evolutionary processes. The relative influences of genetic drift versus selection and phenotypic plasticity in shaping morphology appear to vary among populations occupying nearby streams.

Maryland↗
USGS2009

A comparison of pectoral fin contact between two different wild dolphin populations

Contact behaviour involving the pectoral fin has been documented in a number of dolphin species, and various explanations about its function have been offered. Pectoral fin contact can take a variety of forms, and involves a number of body parts and movements, likely differing depending upon social or ecological context. For this study, we compare the pectoral fin contact behaviour of two species of wild dolphins: Indo-Pacific bottlenose dolphins (Tursiops aduncus) from around Mikura Island, Japan, and Atlantic spotted dolphins (Stenella frontalis) from The Bahamas. The two study populations exhibit surprising similarity in the ways in which pectoral fin contacts are used, despite differences in species and environmental conditions at the two sites. Differences in contact rates for calves between the two sites suggest that calf-focused aggression from adult dolphins is more prevalent at Mikura than in The Bahamas. Our results suggest that pectoral fin contact behaviour seems to be driven primarily by social pressures, and may be similar in function to allogrooming described in primates. ?? 2008 Elsevier B.V.

Behavioural Processes↗
USGS2004

Combined impacts of Black-crowned Night-Heron predation/disturbance and various management activities on Roseate Tern productivity in 2003, and testing of a video surveillance system for recording the diurnal and nocturnal behavior of terns and night-herons at Falkner Island, Stewart B. McKinney National Wildlife Refuge, Connecticut, in 2004: Report to U.S. Fish and Wildlife Service, Stewart B. McKinney National Wildlife Refuge, Westbrook, Connecticut and U.S. Fish and Wildlife Service, Region 5 Regional Office, Hadley, Massachusetts

Falkner Island (FICT), a unit of the U.S. Fish and Wildlife Service's Stewart B. McKinney National Wildlife Refuge (SBMNWR) since 1985, is located in Long Island Sound 5 km south of Guilford, CT. For more than three decades it has been the site of the only large breeding colony in Connecticut of the federally endangered Northwest Atlantic population of Roseate Terns (Sterna dougallii) and the state's largest colony of Common Terns (S. hirundo). Both species have been studied at this site since 1978 as part of the Falkner Island Tern Project (FITP), and since 1987 also as part of a regional Cooperative Roseate Tern Metapopulation Dynamics and Ecology Project (CRTMP), both coordinated by Dr. Jeffrey A. Spendelow of the U.S. Geological Survey's Patuxent Wildlife Research Center (USGS-PWRC). From 1997-2002 the Roseate Tern breeding population at this site declined by more than 50% from about 150 to about 70 nesting pairs, mostly as a result of the nocturnal predation and disturbance of tern chicks and eggs by Black-crowned Night-Herons (Nycticorax nycticorax). Here we report the results of research done with the goal of improving management of nocturnal predators and developing new practices/structures to reduce losses of tern eggs and chicks so as to prevent the abandonment of this site by Roseate Terns. Notification of release of the USGS 'Quick Response Funds' (QRF) that were to be used to support the part of this study entitled 'Nocturnal behavior/interactions of endangered Roseate Terns and Black-crowned Night-Herons', and final approval of the Study Plan for this research did not occur until after the breeding season in 2003 was well underway. As a result, some work will need to be completed during the 2004 field season. There are two major objectives of this study. The first is to collect basic information (a) on the nocturnal behavior and interactions of Roseate (and Common) Terns with predatory Black-crowned Night-Herons, and (b) on how the behavior of the night-herons and the behavior and productivity of the terns are affected by the activities of management personnel (or researchers) that may be moving around the island at night. The second objective is to develop new structures (or to modify existing ones) that can be used by adult Roseate Terns for nesting sites and/or by chicks as hiding sites so that tern eggs and chicks will be less susceptible to night-heron predation. Work on the development and comparative use of various man-made structures began prior to the start of the 2003 nesting season; work on the evaluation and testing of a video surveillance system to study the nocturnal behavior of the terns and night-herons began in June 2003. Reported here are (1) the results of work on the placement and use by Roseate Terns of several types of structures in different habitats, and (2) a description of the set-up, testing, and evaluation of a multi-camera video system to record tern and predator behavior during both day and night. Also given here is some additional information on other factors that affected the productivity of the Roseate Terns at FICT collected as part of the CRTMP which served as the underlying foundation for the additional work done with support of the 2003 Quick Response Funds.

Book↗
USGS2009

Temporal and maternal effects on reproductive ecology of the giant gartersnake ( Thamnophis gigas )

We used mixed-effects models to examine relationships of reproductive characteristics of the giant gartersnake ( Thamnophis gigas ) to improve population modeling and conservation planning for this species. Neonates from larger litters had lower mass, and mass of neonates also was affected by random variation among mothers. Length of mother did not affect relative mass of litters; however, our data suggest that longer mothers expended less reproductive effort per offspring than shorter mothers. We detected random variation in length of neonates among mothers, but these lengths were not related to length of mother or size of litter. Mean size of litter varied among years, but little evidence existed for a relationship between size of litter or mass of litter and length of mother. Sex ratios of neonates did not differ from 1:1.

Southwestern Naturalist↗
USGS2019

Public access for pheasant hunters: Understanding an emerging need

Ring-necked pheasant ( Phasianus colchicus ; i.e., pheasant) hunting participation is declining across North America, reflecting a larger downward trend in American hunting participation and threatening benefits to grassland conservation and rural economies. To stabilize and expand the pheasant hunting population, we must first identify factors that influence pheasant hunter participation. We used an extensive in-person hunter survey to test the hypothesis that hunter demographics interact with social-ecological traits of hunting locations to affect hunter decisions, outcomes, and perceptions. We built a series of Bayesian mixed effects models to parse variation in demographics, perceptions, and hunt outcomes of pheasant hunters interviewed at public access hunting sites across 3 regions in Nebraska, USA, that varied in pheasant abundance and proximity to urban population centers. Among pheasant hunters in Nebraska, access to private lands was negatively related to the human population density of a pheasant hunter's home ZIP code and the distance a hunter had traveled to reach a hunting location. Pheasant hunters interviewed closer to metropolitan areas tended to be more urban and travel shorter distances, and their parties were more likely to include youth but less likely to include dogs. Hunter satisfaction was positively associated with seeing and harvesting pheasants and hunting with youth. Whereas youth participation and the number of pheasants seen varied by study region, hunter satisfaction did not differ across regions, suggesting that hunters may calibrate their expectations and build their parties based on where they plan to hunt. The variation in hunter demographics across hunting locations and disconnects between social and ecological correlates of hunter satisfaction suggests that diverse pheasant hunting constituencies will be best served by diverse pheasant hunting opportunities.

Nebraska↗
USGS2003

Managing a subsidized predator population: Reducing common raven predation on desert tortoises

Human communities often are an inadvertent source of food, water, and other resources to native species of wildlife. Because these resources are more stable and predictable than those in a natural environment, animals that subsist on them are able to increase in numbers and expand their range, much to the detriment of their competitors and species they prey upon. In the Mojave Desert, common ravens (Corvus corax) have benefited from human-provided resources to increase in population size precipitously in recent years. This trend has caused concern because ravens prey on juvenile desert tortoises (Gopherus agassizii), a federally threatened species. In this paper, I discuss management strategies to reduce raven predation on desert tortoises. The recommendations fall into three categories: (1) managing raven populations by reducing access to anthropogenic resources; (2) removing offending ravens or other birds in specially targeted tortoise management zones; and (3) continuing research on raven ecology, raven behavior, and methods of reducing raven predation on tortoises. I also recommend approaching the problem within an adaptive management framework: management efforts should first be employed as scientific experiments - with replicates and controls - to yield an unbiased assessment of their effectiveness. Furthermore, these strategies should be implemented in concert with actions that reduce other causes of desert tortoise mortality to aid the long-term recovery of their populations. Overall, the approaches outlined in this paper are widely applicable to the management of subsidized predators, particularly where they present a threat to a declining species of prey.

Environmental Management↗
USGS2007

Occupancy estimation and modeling with multiple states and state uncertainty

The distribution of a species over space is of central interest in ecology, but species occurrence does not provide all of the information needed to characterize either the well-being of a population or the suitability of occupied habitat. Recent methodological development has focused on drawing inferences about species occurrence in the face of imperfect detection. Here we extend those methods by characterizing occupied locations by some additional state variable ( e. g., as producing young or not). Our modeling approach deals with both detection probabilities,1 and uncertainty in state classification. We then use the approach with occupancy and reproductive rate data from California Spotted Owls ( Strix occidentalis occidentalis ) collected in the central Sierra Nevada during the breeding season of 2004 to illustrate the utility of the modeling approach. Estimates of owl reproductive rate were larger than naive estimates, indicating the importance of appropriately accounting for uncertainty in detection and state classification.

Ecology↗
USGS2013

Linking anti-predator behaviour to prey demography reveals limited risk effects of an actively hunting large carnivore

Ecological theory predicts that the diffuse risk cues generated by wide-ranging, active predators should induce prey behavioural responses but not major, population- or community-level consequences. We evaluated the non-consumptive effects (NCEs) of an active predator, the grey wolf (Canis lupus), by simultaneously tracking wolves and the behaviour, body fat, and pregnancy of elk (Cervus elaphus), their primary prey in the Greater Yellowstone Ecosystem. When wolves approached within 1 km, elk increased their rates of movement, displacement and vigilance. Even in high-risk areas, however, these encounters occurred only once every 9 days. Ultimately, despite 20-fold variation in the frequency of encounters between wolves and individual elk, the risk of predation was not associated with elk body fat or pregnancy. Our findings suggest that the ecological consequences of actively hunting large carnivores, such as the wolf, are more likely transmitted by consumptive effects on prey survival than NCEs on prey behaviour.

Montana;Wyoming;Idaho↗
USGS2022

Scat as a source of DNA for population monitoring

Sampling fecal droppings (scat) to genetically identify individual animals is an established method for monitoring mammal populations and could be highly useful for monitoring reptile populations. Whereas existing protocols for obtaining DNA from reptile scat focus on analyses of whole, fresh scat deposited during animal handling, the collection of scat naturally deposited by reptiles in situ, as required for non-invasive population monitoring, requires protocols to extract highly degraded DNA. Using surface swabs from such scats can reduce PCR inhibition, ecological impacts of removing scat, and zoonotic risks. We report on three related but independently designed studies of DNA analyses from scat swabs of herbivorous reptiles under natural desert conditions: two free-ranging desert tortoise species (Agassiz's desert tortoise, Gopherus agassizii , California, US, and Morafka's desert tortoise, G. morafkai , Arizona, US) and the common chuckwalla (Sauromalus atar ) (Arizona, US, and Sonora, MX). We analyzed samples from both tortoise species with the same set of 16 microsatellites and chuckwalla samples with four mtDNA markers; studies also varied in swab preservation medium and DNA extraction method. Microsatellite amplification success per sample, defined as ≥9 loci with amplification, varied between studies, with 15% for Agassiz's desert tortoise and 42% Morafka's desert tortoise. For chuckwallas, we successfully amplified and sequenced 50% of samples. We recovered fragments up to 400 bp for tortoises and 980 bp for chuckwallas from scat swab samples. This study demonstrates that genotypes can successfully be obtained from swabs of herbivorous reptile scat collected in the field under natural environmental conditions and emphasizes that repeat amplifications are necessary for genetic identification of individuals from non-invasive samples.

Arizona, California↗
USGS2026

DNA retention in sea lamprey digestive tracts: Insights from controlled feeding experiments

The sea lamprey ( Petromyzon marinus ), a non-native species in the Laurentian Great Lakes, has significantly impacted native fish communities and commercial fisheries, requiring population suppression efforts. While traditional control methods such as lampricides and barriers have reduced sea lamprey population abundance, questions remain regarding sea lamprey dietary composition given the focus of current damage assessments on economically and ecologically important host species. Recent advances in molecular technology offer promising methods of sea lamprey dietary assessment. Specifically, DNA metabarcoding enables species-specific identification of taxonomically diverse prey items from gut and fecal samples, and has proven effective in many taxa, including hematophagous species such as Arctic lamprey ( Lethenteron camtschaticum ) and sea lamprey. However, studies on DNA retention within digestive tracts are limited, particularly given the potential effects of environmental and dietary factors among hematophagous species. We used controlled feeding experiments to understand the effects these factors may have on DNA retention and host detectability within sea lamprey digestive tracts. Additionally, we evaluated the utility of metabarcoding for identifying multiple host species from consecutive feedings. Results indicate that host DNA can be detected up to 30 days post-feeding, with detection probability decreasing with increasing time following feeding. Temperature effects were dependent upon fasting periods, and host-switching trials indicated multiple previous host species could be detected from a single lamprey. Findings provide valuable insights for refining dietary analysis protocols for wild-caught sea lamprey within native and introduced ranges.

Environmental DNA↗
USGS1990

Behavioral toxicology: Stimulating challenges for a growing discipline

Since the early 1970s, contaminants have been shown to affect virtually every aspect of behavior in terrestrial and aquatic organisms. Behavior inte- grates many cellular processes and is essential to the viability of the organism, the population and the community. Therefore, observations of behavior provide a unique toxicological perspective - one that links the biochemical and ecological conse- quences of environmental contamination.

Environmental Toxicology and Chemistry↗
USGS2017

Evaluating population expansion of black bears using spatial capture-recapture

The population of American black bears ( Ursus americanus ) in southern New York, USA has been growing and expanding in range since the 1990s. This has motivated a need to anticipate future patterns of range expansion. We conducted a non-invasive, genetic, spatial capture-recapture (SCR) study to estimate black bear density and identify spatial patterns of population density that are potentially associated with range expansion. We collected hair samples in a 2,519-km 2 study area in southern New York with barbed-wire hair snares and identified individuals and measured genetic diversity using 7 microsatellite loci and 1 sex-linked marker. We estimated a mean density of black bears in the region of 13.7 bears/100 km 2 , and detected a slight latitudinal gradient in density consistent with the documented range expansion. However, elevation and the amounts of forest, crop, and developed landcover types did not influence density, suggesting that bears are using a diversity of resources in this heterogeneous landscape outside their previously described distribution. These results provide the first robust baseline estimates for population density and distribution associated with different landcover types in the expanded bear range. Further, genetic diversity was comparable to that of non-expanding black bear populations in the eastern United States, and in combination with the latitudinal density gradient, suggest that the study area is not at the colonizing front of the range expansion. In addition, the diversity of landcover types used by bears in the study area implies a possible lack of constraints for further northern expansion of the black bear range. Our non-invasive, genetic, spatial capture-recapture approach has utility for studying populations of other species that may be expanding in range because SCR allows for the testing of explicit, spatial ecological hypotheses.

New York↗
USGS2018

Biomonitoring using invasive species in a large Lake: Dreissena distribution maps hypoxic zones

Due to cultural eutrophication and global climate change, an exponential increase in the number and extent of hypoxic zones in marine and freshwater ecosystems has been observed in the last few decades. Hypoxia, or low dissolved oxygen (DO) concentrations, can produce strong negative ecological impacts and, therefore, is a management concern. We measured biomass and densities of Dreissena in Lake Erie, as well as bottom DO in 2014 using 19 high frequency data loggers distributed throughout the central basin to validate a three-dimensional hydrodynamic-ecological lake model. We found that a deep, offshore hypoxic zone was formed by early August, restricting the Dreissena population to shallow areas of the central basin. Deeper than 20 m, where bottom hypoxia routinely develops, only young of the year mussels were found in small numbers, indicating restricted recruitment and survival of young Dreissena . We suggest that monitoring Dreissena distribution can be an effective tool for mapping the extent and frequency of hypoxia in freshwater. In addition, our results suggest that an anticipated decrease in the spatial extent of hypoxia resulting from nutrient management has the potential to increase the spatial extent of profundal habitat in the central basin available for Dreissena expansion.

Lake Erie↗
USGS2023

Out of the frying pan and into the fire: Effects of volcanic heat and other stressors on the conservation of a critically endangered plant in Hawaiʻi

Loss of local biodiversity resulting from abrupt environmental change is a significant environmental problem throughout the world. Extinctions of plants are particularly important yet are often overlooked. Drawing from a case in Hawai‘i, a global hotspot for plant and other extinctions, we demonstrate an effort to better understand and determine priorities for the management of an endangered plant (‘Ihi makole or Portulaca sclerocarpa ) in the face of rapid and extreme environmental change. Volcanic heat emissions and biological invasions have anecdotally been suggested as possible threats to the species. We integrated P. sclerocarpa outplanting with efforts to collect geological and ecological data to gauge the role of elevated soil temperatures and invasive grasses in driving P. sclerocarpa mortality and population decline. We measured soil temperature, soil depth, surrounding cover and P. sclerocarpa survivorship over three decades. The abundance of wild P. sclerocarpa decreased by 99.7% from the 1990s to 2021. Only 51% of outplantings persisted through 3–4 years. Binomial regression and structural equation modelling revealed that, among the variables we analysed, high soil temperatures were most strongly associated with population decline. Finding the niche where soil temperatures are low enough to allow P. sclerocarpa survival but high enough to limit other agents of P. sclerocarpa mortality may be necessary to increase population growth of this species.

Hawaii↗
USGS2006

Carnivore re-colonisation: Reality, possibility and a non-equilibrium century for grizzly bears in the southern Yellowstone ecosystem

Most large native carnivores have experienced range contractions due to conflicts with humans, although neither rates of spatial collapse nor expansion have been well characterised. In North America, the grizzly bear (Ursus arctos) once ranged from Mexico northward to Alaska, however its range in the continental USA has been reduced by 95-98%. Under the U. S. Endangered Species Act, the Yellowstone grizzly bear population has re-colonised habitats outside Yellowstone National Park. We analysed historical and current records, including data on radio-collared bears, (1) to evaluate changes in grizzly bear distribution in the southern Greater Yellowstone Ecosystem (GYE) over a 100-year period, (2) to utilise historical rates of re-colonisation to project future expansion trends and (3) to evaluate the reality of future expansion based on human limitations and land use. Analysis of distribution in 20-year increments reflects range reduction from south to north (1900-1940) and expansion to the south (1940-2000). Expansion was exponential and the area occupied by grizzly bears doubled approximately every 20 years. A complementary analysis of bear occurrence in Grand Teton National Park also suggests an unprecedented period of rapid expansion during the last 20-30 years. The grizzly bear population currently has re-occupied about 50% of the southern GYE. Based on assumptions of continued protection and ecological stasis, our model suggests total occupancy in 25 years. Alternatively, extrapolation of linear expansion rates from the period prior to protection suggests total occupancy could take > 100 years. Analyses of historical trends can be useful as a restoration tool because they enable a framework and timeline to be constructed to pre-emptively address the social challenges affecting future carnivore recovery. ?? 2004 The Zoological Society of London.

Idaho, Montana, Wyoming↗
USGS2014

Annual variation of spawning Cutthroat Trout in a small Western USA stream: A case study with implications for the conservation of potamodromous trout life history diversity

Little is known about the variability in the spatial and temporal distribution of spawning potamodromous trout despite decades of research directed at salmonid spawning ecology and the increased awareness that conserving life history diversity should be a focus of management. We monitored a population of fluvial–resident Bonneville Cutthroat Trout Oncorhynchus clarkii utah in a tributary to the Logan River, Utah, from 2006 to 2012 to gain insight into the distribution and timing of spawning and what factors may influence these spawning activities. We monitored Bonneville Cutthroat Trout using redd surveys with multiple observers and georeferenced redd locations. We documented an extended spawning period that lasted from late April to mid‐July. The onset, median, and end of spawning was best predicted by the mean maximum water temperature during the first 13 weeks of the year ( F = 130. 4, df = 5, R 2 = 0.96, P < 0.0001) with spawning beginning and ending earlier in years that had warmer water temperatures prior to spawning. The distribution of redds was clumped each year and the relative density of redds was greater in a reach dominated by dams constructed by beavers Castor canadensis . Both dam failure and construction appeared to be responsible for creating new spawning habitat that was quickly occupied, demonstrating rapid temporal response to local habitat changes. Bonneville Cutthroat Trout appeared to establish and defend a redd for up to 2 d, and spawning most often occurred between similar‐sized individuals. Spawning surveys for potamodromous trout are an underutilized tool that could be used to better understand the distribution and timing of spawning as well as determine the size and trends of the reproducing portion of populations of management concern. Without efforts to document the diversity of this important aspect of potamodromous trout life history, prioritization of conservation will be problematic.

Utah↗
USGS2015

Evaluating multi-level models to test occupancy state responses of Plethodontid salamanders

Plethodontid salamanders are diverse and widely distributed taxa and play critical roles in ecosystem processes. Due to salamander use of structurally complex habitats, and because only a portion of a population is available for sampling, evaluation of sampling designs and estimators is critical to provide strong inference about Plethodontid ecology and responses to conservation and management activities. We conducted a simulation study to evaluate the effectiveness of multi-scale and hierarchical single-scale occupancy models in the context of a Before-After Control-Impact (BACI) experimental design with multiple levels of sampling. Also, we fit the hierarchical single-scale model to empirical data collected for Oregon slender and Ensatina salamanders across two years on 66 forest stands in the Cascade Range, Oregon, USA. All models were fit within a Bayesian framework. Estimator precision in both models improved with increasing numbers of primary and secondary sampling units, underscoring the potential gains accrued when adding secondary sampling units. Both models showed evidence of estimator bias at low detection probabilities and low sample sizes; this problem was particularly acute for the multi-scale model. Our results suggested that sufficient sample sizes at both the primary and secondary sampling levels could ameliorate this issue. Empirical data indicated Oregon slender salamander occupancy was associated strongly with the amount of coarse woody debris (posterior mean = 0.74; SD = 0.24); Ensatina occupancy was not associated with amount of coarse woody debris (posterior mean = -0.01; SD = 0.29) . Our simulation results indicate that either model is suitable for use in an experimental study of Plethodontid salamanders provided that sample sizes are sufficiently large. However, hierarchical single-scale and multi-scale models describe different processes and estimate different parameters. As a result, we recommend careful consideration of study questions and objectives prior to sampling data and fitting models.

PLoS ONE↗
USGS2019

Where has turtle ecology been, and where is it going?

Over 9000 articles have been published on turtles and tortoises (excluding sea turtles) since 1950 according to the Web of Science, including over 8000 contained in a personal bibliography that we analyze in this paper. Research had a slow start from 1900 to 1950, with mostly anecdotal additions to our knowledge until the contributions of F. Cagle and A. Carr took turtle research to new levels as the cofathers of turtle ecology in the middle of the last century. Books written in 1939, 1952, and 1972 that compiled existing literature on turtles in the United States and Canada set the stage for growing interest in turtles. The first global compilation of turtle species did not become available until 1961. Publication frequency increased in the 1960s and especially the 1970s as interest in turtles grew, and a wave of turtle biologists emerged from doctoral degree programs. We briefly review the contributions of scientists who published extensively on turtle ecology in those and later decades up to the present. We also review advances in our knowledge of various topics, including the global distribution of turtle research efforts; changes in our perceptions of turtle species diversity over time; turtle community ecology; sex ratios, sex-determination, and climate change; overwintering behavior; sexual size dimorphism and sexual dichromatism; analyses of population genetics; turtles and vocalization; and the emergence of turtle conservation biology efforts. We conclude with a discussion of future opportunities and challenges for working with turtles.

Herpetologica↗