Search USGS⌕ Search

SEARCH · Search USGS

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 613 records · Page 34Linked to original sources

Robustness of closed capture-recapture methods to violations of the closure assumption

Closed-population capture–recapture methods have been used extensively in animal ecology, both by themselves and within the context of Pollock’s robust design and multistate models, to estimate various parameters of population and community dynamics. The defining assumption of geographic and demographic closure (i.e., no births, deaths, immigration, or emigration) for the duration of sampling is restrictive and is likely to be violated in many field situations. I evaluated several types of violations of the closure assumption and found that completely random movement in and out of a study area does not introduce bias to estimators from closed-population methods, although it decreases precision. In addition, if capture probabilities vary only with time, the closed-population Lincoln-Petersen estimator is unbiased for the size of the superpopulation when there are only births/immigration or only deaths/emigration. However, for other cases of nonrandom movement, closed-population estimators were biased when movement was Markovian (dependent on the presence/absence of the animal in the previous time period), when an animal was allowed one entry to and one exit from the study area, or when there was trap response or heterogeneity among animals in capture probability. In addition, the probability that an animal is present and available for capture (e.g., breeding propensity) can be estimated using Pollock’s robust design only when movement occurs at a broader temporal scale than that of sampling.

Ecology↗

Comparative demography of new world populations of thrushes (Turdus spp.): Comment

Survival and fecundity are fundamental to the study of evolutionary ecology, as they are two of the key variables influencing the constrained optimization process we call natural selection. Likewise, population managers require accurate estimates of these parameters; along with dispersal, they govern population dynamics ( Lack 1954 ) and thus are essential for predicting population change and the effects of management actions. It can be frustrating, then, when survival is difficult to estimate for certain populations, as in some avian species that exhibit a substantial degree of dispersal.

Ecology↗

Pinyon and juniper encroachment into sagebrush ecosystems impacts distribution and survival of greater sage-grouse

In sagebrush ( Artemisia spp.) ecosystems, encroachment of pinyon ( Pinus spp.) and juniper ( Juniperus spp.; hereafter, “pinyon-juniper”) trees has increased dramatically since European settlement. Understanding the impacts of this encroachment on behavioral decisions, distributions, and population dynamics of greater sage-grouse (Centrocercus urophasianus) and other sagebrush obligate species could help benefit sagebrush ecosystem management actions. We employed a novel two-stage Bayesian model that linked avoidance across different levels of pinyon-juniper cover to sage-grouse survival. Our analysis relied on extensive telemetry data collected across 6 yr and seven subpopulations within the Bi-State Distinct Population Segment (DPS), on the border of Nevada and California. The first model stage indicated avoidance behavior for all canopy cover classes on average, but individual grouse exhibited a high degree of heterogeneity in avoidance behavior of the lowest cover class (e.g., scattered isolated trees). The second stage modeled survival as a function of estimated avoidance parameters and indicated increased survival rates for individuals that exhibited avoidance of the lowest cover class. A post hoc frailty analysis revealed the greatest increase in hazard (i.e., mortality risk) occurred in areas with scattered isolated trees consisting of relatively high primary plant productivity. Collectively, these results provide clear evidence that local sage-grouse distributions and demographic rates are influenced by pinyon-juniper, especially in habitats with higher primary productivity but relatively low and seemingly benign tree cover. Such areas may function as ecological traps that convey attractive resources but adversely affect population vital rates. To increase sage-grouse survival, our model predictions support reducing actual pinyon-juniper cover as low as 1.5%, which is lower than the published target of 4.0%. These results may represent effects of pinyon-juniper cover in areas with similar ecological conditions to those of the Bi-State DPS, where populations occur at relatively high elevations and pinyon-juniper is abundant and widespread.

Rangeland Ecology and Management↗

Monitoring trends in bat populations of the United States and territories: Problems and prospects

Bats are ecologically and economically important mammals. The life histories of bats (particularly their low reproductive rates and the need for some species to gather in large aggregations at limited numbers of roosting sites) make their populations vulnerable to declines. Many of the species of bats in the United States (U.S.) and territories are categorized as endangered or threatened, have been candidates for such categories, or are considered species of concern. The importance and vulnerability of bat populations makes monitoring trends in their populations a goal for their future management. However, scientifically rigorous monitoring of bat populations requires well-planned, statistically defensible efforts. This volume reports findings of an expert workshop held to examine the topic of monitoring populations of bats. The workshop participants included leading experts in sampling and analysis of wildlife populations, as well as experts in the biology and conservation of bats. Findings are reported in this volume under two sections. Part I of the report presents contributed papers that provide overviews of past and current efforts at monitoring trends in populations of bats in the U.S. and territories. These papers consider current techniques and problems, and summarize what is known about the status and trends in populations of selected groups of bats. The contributed papers in Part I also include a description of the monitoring program developed for bat populations in the United Kingdom, a critique of monitoring programs in wildlife in general with recommendations for survey and sampling strategies, and a compilation and analysis of existing data on trends in bats of the U.S. and territories. Efforts directed at monitoring bat populations are piecemeal and have shortcomings. In Part II of the report, the workshop participants provide critical analyses of these problems and develop recommendations for improving methods, defining objectives and priorities, gaining mandates, and enhancing information exchange to facilitate future efforts for monitoring trends in U.S. bat populations.

Information and Technology Report↗

Large-scale variation in wave attenuation of oyster reef living shorelines and the influence of inundation duration

One of the paramount goals of oyster reef living shorelines is to achieve sustained and adaptive coastal protection, which requires meeting ecological (i.e., develop a self-sustaining oyster population) and engineering (i.e., provide coastal defense) targets. In a large-scale comparison along the Atlantic and Gulf coasts of the United States, the efficacy of various designs of oyster reef living shorelines at providing wave attenuation was evaluated accounting for the ecological limitations of oysters with regard to inundation duration. A critical threshold for intertidal oyster reef establishment is 50% inundation duration. Living shorelines that spent less than one-half of the time (<50%) inundated were not considered suitable habitat for oysters, however, were effective at wave attenuation (68% reduction in wave height). Reefs that experienced >50% inundation were considered suitable habitat for oysters, but wave attenuation was similar to controls (no reef; ~5% reduction in wave height). Many of the oyster reef living shoreline approaches therefore failed to optimize the ecological and engineering goals. In both inundation regimes, wave transmission decreased with an increasing freeboard (difference between reef crest elevation and water level), supporting its importance in the wave attenuation capacity of oyster reef living shorelines. However, given that the reef crest elevation (and thus freeboard) should be determined by the inundation duration requirements of oysters, research needs to be refocused on understanding the implications of other reef parameters (e.g., width) for optimizing wave attenuation. A broader understanding of the reef characteristics and seascape contexts that result in effective coastal defense by oyster reefs is needed to inform appropriate design and implementation of oyster-based living shorelines globally.

Alabama, Florida, Louisiana, New Jersey, Virginia↗

Abundance of introduced species at home predicts abundance away in herbaceous communities

Many ecosystems worldwide are dominated by introduced plant species, leading to loss of biodiversity and ecosystem function. A common but rarely tested assumption is that these plants are more abundant in introduced vs. native communities, because ecological or evolutionary-based shifts in populations underlie invasion success. Here, data for 26 herbaceous species at 39 sites, within eight countries, revealed that species abundances were similar at native (home) and introduced (away) sites - grass species were generally abundant home and away, while forbs were low in abundance, but more abundant at home. Sites with six or more of these species had similar community abundance hierarchies, suggesting that suites of introduced species are assembling similarly on different continents. Overall, we found that substantial changes to populations are not necessarily a pre-condition for invasion success and that increases in species abundance are unusual. Instead, abundance at home predicts abundance away, a potentially useful additional criterion for biosecurity programmes.

Ecology Letters↗

Range-wide population trend analysis for greater sage-grouse (Centrocercus urophasianus)—Updated 1960–2021

Greater sage-grouse ( Centrocercus urophasianus ) are at the center of state and national land use policies largely because of their unique life-history traits as an ecological indicator for health of sagebrush ecosystems. This updated population trend analysis provides state and federal land and wildlife managers with best-available science to help guide current management and conservation plans aimed at benefitting sage-grouse populations. This analysis relied on previously published population trend modeling methodology from Coates and others (2021) and includes the addition of three analytical updates: (1) identification of population nadirs (lowest points within cycles) at the lek (breeding ground) and neighborhood cluster (group of leks) spatial scales, (2) truncation of prior distributions on rate of change in apparent abundance values to more realistic boundaries for leks with missing data, and (3) addition of 2 years of population lek count data (2020 and 2021) to the current dataset (1953–2021). Bayesian state-space models estimated 2.9 percent average annual decline in sage-grouse populations across their geographical range, which varied among subpopulations at the largest scale of analysis, termed climate clusters (2.2–4.6). Cumulative declines were 42.5, 65.6, and 80.1 percent range-wide across short (19 years), medium (35 years), and long (55 years) temporal periods, respectively. These results indicate that range-wide populations continued to decline during 2020 and 2021, although two climate clusters (eastern area and Bi-State area) have shown growth in population abundance in recent years, indicating they have surpassed a recent population abundance nadir.

Data Report↗

Using stable isotopes to investigate individual diet specialization in California sea otters (Enhydra lutris nereis)

Differences in diet composition among conspecifics (dietary specialization) have been documented across a broad range of taxonomic groups and habitats, and such variation at the individual level is increasingly recognized as an important component of diversity in trophic interactions. Accurate identification of individual dietary specialization, however, requires longitudinal dietary records that are labor-intensive and cost-prohibitive to obtain for many species. Here we explore the use of stable isotopes (δ 13 C and δ 15 N) as a promising technique for detecting and quantifying patterns of individual dietary specialization. Southern sea otters ( Enhydra lutris nereis ) offer a unique opportunity for testing this approach because (1) they consume a wide variety of prey that span multiple trophic levels, habitats, and ecologically defined functional groups; and (2) individual diet specialization can be validated with existing observational data. We analyzed the isotopic composition of sea otter vibrissae ( n = 31) in order to characterize inter- and intra-individual variation in sea otter diets at Monterey Bay, California, USA. At the population level, sea otters showed substantial variation in both δ 13 C and δ 15 N values, occupying nearly all of the “isotopic space” created by the diversity of isotopic signatures of potential prey taxa. Most of the variation in sea otter vibrissae was accounted for by differences between individuals, with much less contributed by within-individual variation. A majority of sea otters (∼80%) showed relatively little temporal variability in isotopic composition, suggesting that the proportional composition of most individuals' diets is relatively constant over time; a few individuals (∼20%) exhibited a high degree of intra-vibrissa isotopic variability, suggesting seasonal shifts in diet composition. These results and our interpretation of them were supported by long-term observational data on the diets of radio-tagged sea otters from the same population ( n = 23). Our results demonstrate that stable isotopes can provide an efficient tool for measuring individual- and population-level dietary breadth and may be useful for studying populations where longitudinal data on individuals would otherwise be impossible to acquire. This will be critical for examining the causes and consequences of dietary variation within and among consumer populations, thereby improving our understanding of these important ecological and evolutionary processes at the community level.

Ecology↗

Using spatial capture–recapture to elucidate population processes and space-use in herpetological studies

The cryptic behavior and ecology of herpetofauna make estimating the impacts of environmental change on demography difficult; yet, the ability to measure demographic relationships is essential for elucidating mechanisms leading to the population declines reported for herpetofauna worldwide. Recently developed spatial capture–recapture (SCR) methods are well suited to standard herpetofauna monitoring approaches. Individually identifying animals and their locations allows accurate estimates of population densities and survival. Spatial capture–recapture methods also allow estimation of parameters describing space-use and movement, which generally are expensive or difficult to obtain using other methods. In this paper, we discuss the basic components of SCR models, the available software for conducting analyses, and the experimental designs based on common herpetological survey methods. We then apply SCR models to Red-backed Salamander ( Plethodon cinereus ), to determine differences in density, survival, dispersal, and space-use between adult male and female salamanders. By highlighting the capabilities of SCR, and its advantages compared to traditional methods, we hope to give herpetologists the resource they need to apply SCR in their own systems.

Journal of Herpetology↗

Range-wide population trend analysis for greater sage-grouse (Centrocercus urophasianus)—Updated 1960–2022

Greater sage-grouse ( Centrocercus urophasianus ) are at the center of state and national land-use policies largely because of their unique life-history traits as an ecological indicator for health of sagebrush ecosystems. This updated population trend analysis provides state and federal land and wildlife managers with best-available science to help guide current management and conservation plans aimed at benefitting sage-grouse populations. This analysis relied on previously published population trend modeling methodology from Coates and others (2021, 2022a) and incorporated population lek count data through 2022. Bayesian state-space models estimated 2.9 percent average annual decline in sage-grouse populations across their geographical range, which varied among subpopulations at the largest scale of analysis, termed climate clusters (2.2–4.7). Cumulative declines were 40.9, 65.0, and 79.6 percent range-wide across short (19 years), medium (35 years), and long (55 years) temporal periods, respectively. These results indicate that the most recent nadir for range-wide populations occurred during 2021. However, growth during 2022 was modest, making 2021 a tentative final nadir at this point.

Data Report↗

Range-wide population trend analysis for greater sage-grouse (Centrocercus urophasianus)—Updated 1960–2023

Greater sage-grouse ( Centrocercus urophasianus ) are at the center of state and national land-use policies largely because of their unique life-history traits as an ecological indicator for health of sagebrush ecosystems. This updated population trend analysis provides state and federal land and wildlife managers with best-available science to help guide management and conservation plans aimed at benefitting sage-grouse populations. This analysis relied on previously published population trend modeling methodology from Coates and others (2021, 2022a) and incorporates population lek count data for 1960–2023. Included in this update are changes in terminology. Specifically, we now use the terms Period 1 (previously Long), Period 2 (previously Medium/Long), Period 3 (previously Medium), Period 4 (previously Short/Medium), Period 5 (previously Short), and Period 6 (previously Recent) to identify specific trends. State-space models estimated 2.8-percent average annual decline in sage-grouse populations between 1966 and 2021 (Period 1, six population oscillations) across their geographical range. Average annual decline among climate clusters for the same number of oscillations ranged between 2.1 and 3.1 percent. Cumulative declines were 41.1, 64.5, and 78.4 percent range-wide during Period 5 (19 years), Period 3 (35 years), and Period 1 (55 years), respectively. Population growth during 2022 and 2023 continue to point to 2021 as the most recent range-wide nadir.

Data Report↗

Forecasting sea otter recolonization: Insights from isotopic analysis of modern and zooarchaeological populations

Retrospective datasets offer essential context for conservation by revealing species’ ecological roles before industrial-era human impacts. We analysed isotopic compositions of pre-industrial and modern sea otters ( Enhydra lutris ) to reconstruct pre-extirpation ecology and offer insights for management. Our study focuses on southeast Alaska (SEAK), where sea otters are recolonizing, and northern Oregon, where translocations are being considered. We measured bulk bone collagen δ 13 C and δ 15 N values and essential amino acid δ 13 C values of extirpated sea otters from archaeological contexts, and bulk isotopic values from vibrissae of modern SEAK sea otters. We compare these results with published isotopic data of potential prey and additional archaeological datasets. In SEAK, our data show pre-industrial sea otter populations consumed infaunal bivalves and used soft-sediment (33%) and kelp forest habitats (67%), with sub-regional variation. We anticipate current populations will expand into this historical niche, and conflict with regional traditional/subsistence bivalve fisheries will persist. In northern Oregon, isotopic data from extirpated sea otters indicate past consumption of low trophic level invertebrates and a stronger reliance on kelp forests (88%) rather than soft-sediment habitats, highlighting the importance of kelp forests for future translocations. Our work exemplifies the value of historical ecology in informing conservation strategies for recovering species.

Alaska, Oregon↗

Informed breeding dispersal following stochastic changes to patch quality in a pond-breeding amphibian

The unidirectional movement of animals between breeding patches (i.e. breeding dispersal) has profound implications for the ecological and evolutionary dynamics of spatially structured populations. In spatiotemporally variable environments, individuals are expected to adjust their dispersal decisions according to information gathered on the environmental and/or social cues that reflect the fitness prospects in a given breeding patch (i.e. informed dispersal). A paucity of empirical work limited our understanding of the ability of animals to depart from low-quality breeding patches and settle in high-quality breeding patches. We examined the capacity of individuals to respond to stochastic changes in habitat quality via informed breeding dispersal in a pond-breeding amphibian. We conducted a 5-year (2015–2019) capture–recapture study of boreal toads Anaxyrus boreas boreas ( n = 1,100) that breed in beaver ponds in western Wyoming, USA. During early spring of 2017, an extreme flooding event destroyed several beaver dams and resulted in the loss of breeding habitat. We used multi-state models to investigate how temporal changes in pond characteristics influenced breeding dispersal, and determine whether movement decisions were in accordance with prospects for reproductive fitness. Boreal toads more often departed from low-quality breeding ponds (without successful metamorphosis) and settled in high-quality breeding ponds (with successful metamorphosis). Movement decisions were context-dependent and associated with pond characteristics altered by beaver dam destruction. Individuals were more likely to depart from shallow ponds with high vegetation cover and settle in deep ponds with low vegetation cover. The probability of metamorphosis was related to the same environmental cues, suggesting that boreal toads assess the fitness prospects of a breeding patch and adjust movement decisions accordingly (i.e. informed breeding dispersal). We demonstrated that stochastic variability in environmental conditions and habitat quality can underpin dispersal behaviour in amphibians. Our study highlighted the mechanistic linkages between habitat change, movement behaviour and prospects for reproductive performance, which is critical for understanding how wild animals respond to rapid environmental change.

Wyoming↗

Long-term spotlight surveys of American alligators in Mississippi, USA

Accurate population estimates and assessments of trajectory are an essential part of harvest management for game species and conservation action plans for protected species. Long-term monitoring can lead to ecological understanding by identifying biotic and abiotic drivers of population dynamics. Spotlight surveys are a widely used method to monitor abundance and size-class structure of crocodilian populations. The American Alligator (Alligator mississippiensis) has recovered from significant population reductions in the southeastern United States. The Mississippi Department of Wildlife, Fisheries, and Parks (MDWFP) has conducted alligator spotlight surveys since 1971 to monitor populations. We analyzed this long-term alligator survey dataset to assess possible trends in counts as a proxy for potential population changes. We tested for a positive trend in count data over 46 y and evaluated covariates that could influence counts to assist future survey protocols. Alligator counts during 1971– 2016 increased across survey routes in Mississippi. This observed positive response may represent an increase of the alligator population in Mississippi as a result of conservation benefits accrued from improved wetland conditions and species-specific management policies. Evaluation of survey covariates indicated recent rainfall and increasing wind velocity had negative effects on alligator counts while increasing water temperature had a positive effect. Implementing robust survey techniques will improve the reliability of alligator monitoring data and their application to the management of alligator populations. Further, these improved approaches may be useful to other conservation and management agencies as well as for other crocodilian species.

Mississippi↗

Assessing hypotheses about nesting site occupancy dynamics

Hypotheses about habitat selection developed in the evolutionary ecology framework assume that individuals, under some conditions, select breeding habitat based on expected fitness in different habitat. The relationship between habitat quality and fitness may be reflected by breeding success of individuals, which may in turn be used to assess habitat quality. Habitat quality may also be assessed via local density: if high‐quality sites are preferentially used, high density may reflect high‐quality habitat. Here we assessed whether site occupancy dynamics vary with site surrogates for habitat quality. We modeled nest site use probability in a seabird subcolony (the Black‐legged Kittiwake, Rissa tridactyla ) over a 20‐year period. We estimated site persistence (an occupied site remains occupied from time t to t + 1) and colonization through two subprocesses: first colonization (site creation at the timescale of the study) and recolonization (a site is colonized again after being deserted). Our model explicitly incorporated site‐specific and neighboring breeding success and conspecific density in the neighborhood. Our results provided evidence that reproductively “successful” sites have a higher persistence probability than “unsuccessful” ones. Analyses of site fidelity in marked birds and of survival probability showed that high site persistence predominantly reflects site fidelity, not immediate colonization by new owners after emigration or death of previous owners. There is a negative quadratic relationship between local density and persistence probability. First colonization probability decreases with density, whereas recolonization probability is constant. This highlights the importance of distinguishing initial colonization and recolonization to understand site occupancy. All dynamics varied positively with neighboring breeding success. We found evidence of a positive interaction between site‐specific and neighboring breeding success. We addressed local population dynamics using a site occupancy approach integrating hypotheses developed in behavioral ecology to account for individual decisions. This allows development of models of population and metapopulation dynamics that explicitly incorporate ecological and evolutionary processes.

Brittany↗

Understanding protected area resilience: a multi-scale, social-ecological approach

Protected areas (PAs) remain central to the conservation of biodiversity. Classical PAs were conceived as areas that would be set aside to maintain a natural state with minimal human influence. However, global environmental change and growing cross-scale anthropogenic influences mean that PAs can no longer be thought of as ecological islands that function independently of the broader social-ecological system in which they are located. For PAs to be resilient (and to contribute to broader social-ecological resilience), they must be able to adapt to changing social and ecological conditions over time in a way that supports the long-term persistence of populations, communities, and ecosystems of conservation concern. We extend Ostrom's social-ecological systems framework to consider the long-term persistence of PAs, as a form of land use embedded in social-ecological systems, with important cross-scale feedbacks. Most notably, we highlight the cross-scale influences and feedbacks on PAs that exist from the local to the global scale, contextualizing PAs within multi-scale social-ecological functional landscapes. Such functional landscapes are integral to understand and manage individual PAs for long-term sustainability. We illustrate our conceptual contribution with three case studies that highlight cross-scale feedbacks and social-ecological interactions in the functioning of PAs and in relation to regional resilience. Our analysis suggests that while ecological, economic, and social processes are often directly relevant to PAs at finer scales, at broader scales, the dominant processes that shape and alter PA resilience are primarily social and economic.

Ecological Applications↗

Detecting spatial ontogenetic niche shifts in complex dendritic ecological networks

Ontogenetic niche shifts (ONS) are important drivers of population and community dynamics, but they can be difficult to identify for species with prolonged larval or juvenile stages, or for species that inhabit continuous habitats. Most studies of ONS focus on single transitions among discrete habitat patches at local scales. However, for species with long larval or juvenile periods, affinity for particular locations within connected habitat networks may differ among cohorts. The resulting spatial patterns of distribution can result from a combination of landscape-scale habitat structure, position of a habitat patch within a network, and local habitat characteristics—all of which may interact and change as individuals grow. We estimated such spatial ONS for spring salamanders ( Gyrinophilus porphyriticus ), which have a larval period that can last 4 years or more. Using mixture models to identify larval cohorts from size frequency data, we fit occupancy models for each age class using two measures of the branching structure of stream networks and three measures of stream network position. Larval salamander cohorts showed different preferences for the position of a site within the stream network, and the strength of these responses depended on the basin-wide spatial structure of the stream network. The isolation of a site had a stronger effect on occupancy in watersheds with more isolated headwater streams, while the catchment area, which is associated with gradients in stream habitat, had a stronger effect on occupancy in watersheds with more paired headwater streams. Our results show that considering the spatial structure of habitat networks can provide new insights on ONS in long-lived species.

Virginia↗

Influences of supplemental feeding on winter elk calf:cow ratios in the southern Greater Yellowstone Ecosystem

Several elk herds in the Greater Yellowstone Ecosystem are fed during winter to alleviate interactions with livestock, reduce damage to stored crops, and to manage for high elk numbers. The effects of supplemental feeding on ungulate population dynamics has rarely been examined, despite the fact that supplemental feeding is partially justified as necessary for maintaining or enhancing population growth rates. We used linear regression to assess how the presence of feedgrounds, snowpack, summer rainfall, indices of grizzly bear density and wolves per elk, elk population trend counts, brucellosis seroprevalence, and survey date were correlated with midwinter calf:cow ratios, a metric correlated with population growth, from 1983&ndash;2010 from 12 ecologically similar elk herd units (7 fed and 5 unfed) in Wyoming, USA. Our statistical approach allowed for rigorous tests of the hypotheses that supplemental feeding had positive effects on calf:cow ratios and reduced sensitivity of calf:cow ratios to bottom-up limitation relative to top-down limitation from native predators. Calf:cow ratios generally declined across all herd units over the study period and varied widely among units with feedgrounds. We found no evidence that the presence of feedgrounds had positive effects on midwinter calf:cow ratios in Wyoming. Further, fed elk showed stronger correlations with environmental factors, whereas calf:cow ratios for unfed elk showed stronger correlations with predator indices. Although we found no consistent association between winter feeding and higher calf:cow ratios, we did not assess late winter mortality and differences in human offtake between fed and unfed regions, which remain a priority for future research.

Wyoming↗