Geology Reports⌕ Search

Geology topics

Research about Idaho, Montana, Wyoming

Source-linked reports with geographic coverage including Idaho, Montana, Wyoming.

At least 19 recordsLinked to original sources

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗

Predicted genetic consequences of alternative population control strategies for North American plains bison in Yellowstone National Park

North American plains bison ( Bison bison bison ) in Yellowstone National Park are descendant from a population low of 23 indigenous individuals from the late 1800s, and 21 individuals introduced from outside herds in the early 1900s to strengthen the possibility of recovery. Within the park, the herd has rebounded and now averages around 5,000 animals. The herd is managed as a closed population within the park and adjacent areas in Montana, USA, which involves annual removals of animals to control numbers. We used genetic simulation to first isolate the effects of population size on retention of genetic diversity, and second to predict the genetic consequences of different population control strategies that vary in the number of bison removed annually, and their relatedness, while considering varying degrees of male dominance and environmental catastrophes. We found through simulation of a stable population that a herd size >2,500 had a high probability (>0.9) of retaining genetic diversity at levels compatible with long-term conservation (>90% of the initial genetic variation over 200 years), but a herd of 3,500 or more had a high likelihood (>0.9) of retaining levels of genetic variation (95%) suggested for maintaining evolutionary potential. When simulating the removals based on plausible management strategies available to the National Park Service, we found that maintaining a herd size of 3,500 or more, limiting removals to less than 40% of the population at a time, keeping a balanced sex ratio, and preferentially removing related individuals would likely retain more than 95% of existing genetic variation. With a herd size of 3,500 animals or more, population catastrophes and male dominance only influenced genetic retention in the most extreme cases, such as a catastrophic population decline of over 95%, a prolonged decline from 90% to 10% over 5 years, or extreme male dominance where 1% of males produced all offspring each year.

Idaho, Montana, Wyoming↗

Yellowstone Volcano Observatory 2024 annual report

The Yellowstone Volcano Observatory (YVO) monitors volcanic and hydrothermal activity associated with the Yellowstone magmatic system, carries out research into magmatic processes occurring beneath Yellowstone Caldera, and issues timely warnings and guidance related to potential future geologic hazards. YVO is a collaborative consortium that includes the U.S. Geological Survey (USGS), Yellowstone National Park, University of Utah, University of Wyoming, Montana State University, EarthScope Consortium, Wyoming State Geological Survey, Montana Bureau of Mines and Geology, and Idaho Geological Survey. The USGS component of YVO also has the operational responsibility for monitoring volcanic activity in the Intermountain West of the United States, including Arizona, New Mexico, Utah, and Colorado. This report summarizes the activities and findings of YVO during the year 2024, focusing on the Yellowstone volcanic system.

Idaho, Montana, Wyoming↗

Genetic structure in a previously extirpated population of gray wolves following reintroduction and natural recolonization

Genetic structuring in wildlife populations is driven by barriers that restrict gene flow as well as the history of population demography. Mechanisms driving genetic structuring can be nuanced in group-living species, such as gray wolves ( Canis lupus ). Behavioral factors, such as social affiliation and resistance, natal habitat imprinting, and trade-offs between dispersal from natal packs and territorial biding, affect habitat selection of wolves despite landscape barriers providing little resistance to their extensive dispersal capabilities. Wolves were previously extirpated from Idaho, USA, and current populations are the result of both reintroductions in 1995 and 1996 and natural dispersal from Canada. In this context we examined genetic structure of wolves in Idaho using 101 individuals genotyped at 18 nuclear DNA microsatellite loci and a subset of 38 individuals genotyped at 1019 single nucleotide polymorphism markers. We hypothesized panmictic (i.e., random mating) genetic structure in Idaho due to the long-distance dispersal abilities of gray wolves. Contrary to our hypothesis, we found three genetic clusters of gray wolves in Idaho, primarily supported by SNP markers. Microsatellite data suggested similar patterns, but permutation tests indicated these differences were not statistically significant. The extent of differentiation and evidence of gene flow, however, suggests that the three genetic clusters are not wholly isolated from one another. The distinctions between clusters spatially align with areas of reintroduction into central Idaho and Yellowstone National Park, as well ongoing natural recolonization from adjacent populations in Canada and Montana. Wolves at the periphery of analysis areas showed more admixture than those in the core, consistent with territoriality and mating behaviors contributing to genetic structuring. We demonstrate how management history, including reintroduction efforts, and animal behavior may interact and contribute to patterns of genetic structure in wild populations.

Idaho, Montana, Wyoming↗

Luminescence dating of hydrothermal explosions in the Yellowstone Plateau volcanic field

Hydrothermal explosions are a significant geological hazard in some active volcanic systems; however, the timing and triggering mechanisms of these explosions are poorly constrained. This study applies luminescence dating techniques to hydrothermal explosion deposits in the Yellowstone Plateau volcanic field to constrain explosion chronologies and evaluate potential triggering mechanisms. We tested four luminescence dating techniques: K-feldspar post-infrared infrared stimulated luminescence (pIRIR225), quartz blue light optically stimulated luminescence (BLOSL), quartz blue thermoluminescence (BTL), and quartz red thermoluminescence (RTL). The pIRIR225 and RTL protocols produce consistent age estimates that agree with independent radiocarbon ages and with the timing of the Pinedale deglaciation. This study focuses on two craters, Mary Bay, along the northern shore of Yellowstone Lake, and Pocket Basin in Lower Geyser Basin. The mean pIRIR225 ages from Mary Bay deposits (11.99 ± 0.68 ka) agree with previous radiocarbon constraints. The mean pIRIR225 results from Pocket Basin deposits (13.44 ± 1.06 ka) suggest a history of explosion following Pinedale deglaciation, followed by recent hydrothermal alteration. Luminescence dating techniques are a promising tool for reconstructing the timing of hydrothermal explosions in the Late Pleistocene and Holocene, helping to constrain recurrence intervals of the largest hydrothermal systems, informing risk, and improving hazard assessments.

Idaho, Montana, Wyoming↗

Dismembering the Lava Creek Tuff, Yellowstone volcano

The youngest large explosive eruption of Yellowstone volcano generated the 0.63 Ma Lava Creek Tuff (LCT) ignimbrite and a voluminous, widespread fall deposit accompanying formation of Yellowstone Caldera. The tuff is conventionally divided into members A (LCT-A) and B (LCT-B), defined by a reversal in density/welding profiles and the presence of phenocrystic amphibole in LCT-A. Remapping within the northeastern caldera sector has, however, documented multiple ignimbrites that yield LCT eruption ages but are distinct from documented members A and B. Thus, the framework by which the LCT was originally subdivided is re-examined here using field observations, petrologic analyses, and sanidine geochemistry from LCT-A and LCT-B ignimbrite (as previously mapped) at reference localities around the caldera. Mineral modal abundances in both LCT-A and LCT-B outcrops are dominated by quartz, followed by sanidine, with sizes and abundances that generally increase up-section, as previously documented. Amphibole is found by us, however, only in western and southern sectors, in deposits mapped as both LCT-A and LCT-B, and only as a trace phase. Further, a locality reported as LCT-A with amphibole present yields a sanidine 40 Ar/ 39 Ar date of 0.775 ± 0.002 Ma, precluding any LCT correlation. Ignimbrite samples from mapped LCT-A and LCT-B exhibit variable crystal sizes, sanidine cathodoluminescence image populations and Ba concentrations, with no clear distinction between the mapped members. Collectively this research found no evidence for a definable LCT-A versus LCT-B, meaning that the traditional A and B subdivision does not accurately represent the eruptive architecture of the Lava Creek Tuff and warrants revision.

Idaho, Montana, Wyoming↗

Bioclimatic, demographic, and anthropogenic correlates of grizzly bear activity patterns in the Greater Yellowstone Ecosystem

Plasticity of diel activity rhythms may be a key element for adaptations of wildlife populations to changing environmental conditions. In the last decades, grizzly bears Ursus arctos in the Greater Yellowstone Ecosystem (GYE) have experienced notable environmental fluctuations, including changes in availability of food sources and severe droughts. Although substantial research has been conducted on grizzly bear diets, space use, and demographic parameters, studies on factors that may influence their diel activity patterns are lacking. We investigated diel activity of grizzly bears in the GYE as a function of anthropogenic landscape modification, maximum daily ambient temperature, drought severity, and bear density. Specifically, we used accelerometry readings of 169 bears (39 females, 130 males) from 2009 to 2022 to compute three complementary activity measures, hourly intensity of activity, daily active minutes, and active bout length, each used as a response variable within a Bayesian modeling framework. Grizzly bears generally exhibited bimodal diel activity, with crepuscular peaks and slight variations across seasons. Females with young (i.e. cubs or yearlings) were an exception, with more pronounced diurnal activity patterns, possibly as a strategy to avoid infanticide by dominant males. Landscape modification and maximum ambient temperature were the factors most strongly associated with activity patterns of grizzly bears, with greater nocturnality observed in lone females and males as these factors increased. Females with young were comparatively less affected. The GYE is changing because of increasing land development, human recreation pressures, and effects of climate change. Given their greater diurnal activity compared with other cohorts, female grizzly bears with dependent offspring may be more constrained in their ability to modify activity patterns. Our findings add to a growing body of research emphasizing the importance of the temporal dimension of wildlife behavior as a critical factor in assessing species adaptability and vulnerability in a changing world.

Idaho, Montana, Wyoming↗

Exploring the exploitation of migratory moths by radio-marked grizzly bears in Wyoming

Grizzly bear ( Ursus arctos ) feeding on migratory army cutworm moths ( Euxoa auxiliaris ) was first documented by the Interagency Grizzly Bear Study Team (IGBST) during the early 1980s in the southeastern portion of the Greater Yellowstone Ecosystem (GYE). Since those initial observations, use of this seasonally available food resource by grizzly bears has grown substantially. As of 2023, we have documented 4,754 observations of grizzly bears feeding or digging at high-elevation talus slopes. We used those records to identify 36 unique moth sites in the GYE, assessed geographic characteristics of these sites, and documented chronology and frequency of use by grizzly bears. We used occurrences of radio-collared grizzly bears to identify a sample of bears that foraged at moth sites and investigated characteristics of fidelity, duration, diel activity, and movement patterns. Grizzly bears exhibited high fidelity to specific moth sites within and across years. Bears showed approximately 50% reductions in movement metrics while using sites, albeit with increased activity patterns. Estimates of feeding metrics by females exceeded those of males in intensity and duration of both use and daily activity. Given increasing human impacts and interest in observing this feeding interaction, quantifying metrics that describe the spatial and temporal patterns of moth site use by grizzly bears could be beneficial for future management. Further analyses are needed to fully examine the relationship between caloric influences of moth use on grizzly bear demographic rates.

Idaho, Montana, Wyoming↗

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗

Unintended indirect effects limit elk productivity from supplemental feeding in the Greater Yellowstone Ecosystem

The widespread practice of supplemental feeding, a bottom-up forcing of resource availability, is intended to improve wildlife population health and survival. However, supplemental feeding could trigger indirect effects by altering predation rates and disease dynamics. We investigated the effects of feeding on three key elk ( Cervus canadensis ) population productivity metrics (calf:cow ratios, annual change in elk density, and harvestable surplus) across 13 regions in the Greater Yellowstone Ecosystem (GYE) over 26 years. Incorporating previous population size, climate, predator, and harvest data in a Bayesian regression framework revealed new insights about elk productivity metrics in the GYE. Supplemental feeding was associated with increased calf:cow ratios (4.9%) but was not substantially related to changes in elk density and harvestable surplus, which are both management targets. Notably, the feeding effect on calf:cow ratios appeared to be offset by increased wolf ( Canis lupus ) and grizzly bear ( Ursus arctos horribilis ) predation. We hypothesize that increased elk productivity resulting from supplemental feeding is primarily transferred to predator and pathogen trophic levels in this system with limited observed effects on elk abundance and harvestable surplus. Anthropogenic food resources may have unintended indirect consequences on other trophic levels that potentially limit the direct impacts of feeding.

Idaho, Montana, Wyoming↗

Hydrologic mechanisms for 2022 Yellowstone River flood and comparisons to recent historic floods

In June 2022, a historic flood event occurred in the headwaters of the Yellowstone River Basin. The flood resulted in millions of dollars in damages and substantial interruptions to Yellowstone National Park. The 2022 flood event was substantially higher in magnitude than other high-peak flow events over the last 30 years. The high discharge was primarily due to the combination of hydrologic mechanisms initiated by rain-on-snow, including a high-elevation snowpack that peaked later than average. However, the contributions of each hydrologic driver, rain and snow, have not been quantified and could be important for understanding future flood events in the region. The contribution of snowmelt to the total terrestrial water input (TWI) varied throughout the area, yet was concentrated in the headwaters of the Yellowstone, Stillwater, and Boulder rivers, along with the headwaters of Rock Creek in Wyoming and Montana. The primary atmospheric contributions to the TWI during the 2022 event were precipitation from moisture transported from the Pacific Ocean that converged over the Greater Yellowstone Area (GYA) and snowmelt from residual snowpack in the northeast part of Yellowstone National Park.

Idaho, Montana, Wyoming↗

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗

The progression of basaltic–rhyolitic melt storage at Yellowstone Caldera

Yellowstone Caldera is one of the largest volcanic systems on Earth, hosting three major caldera-forming eruptions in the past two million years, interspersed with periods of less explosive, smaller-volume eruptions 1 . Caldera-forming eruptions at Yellowstone are sourced by rhyolitic melts stored within the mid- to upper crust. Seismic tomography studies have suggested that a broad region of rhyolitic melt extends beneath Yellowstone Caldera, with an estimated melt volume that is one to four times greater than the eruptive volume of the largest past caldera-forming eruption, and an estimated melt fraction of 6–28 per cent 2 , 3 , 4 , 5 . Seismic velocity is strongly influenced by temperature, pressure and melt; however, magnetotelluric data are primarily sensitive to the presence of melt, making these data ideal for constraining volcanic systems. Here we utilize magnetotelluric data to model the resistivity structure of Yellowstone Caldera’s crustal magma reservoir and constrain the region’s potential for producing major volcanic eruptions. We find that rhyolitic melts are stored in segregated regions beneath the caldera with low melt fractions, indicating that the reservoirs are not eruptible. Typically, these regions have melt volumes equivalent to small-volume post-caldera Yellowstone eruptions. The largest region of rhyolitic melt storage, concentrated beneath northeast Yellowstone Caldera, has a storage volume similar to the eruptive volume of Yellowstone’s smallest caldera-forming eruption. We identify regions of basalt migrating from the lower crust, merging with and supplying heat to the northeast region of rhyolitic melt storage. On the basis of our analysis, we suggest that the locus of future rhyolitic volcanism has shifted to northeast Yellowstone Caldera.

Idaho, Montana, Wyoming↗

June 2022 floods in the Upper Yellowstone River Basin

Extraordinary floods surged down the Yellowstone River and its tributaries in northwestern Wyoming and south-central Montana on June 13–15, 2022. During the flood, U.S. Geological Survey staff worked to maintain real-time data from streamgages by making field measurements of streamflow and repairing damaged equipment while communicating the latest streamflow information with the public and with local, State, and Federal agencies. After the flood, staff surveyed high-water marks, computed peak streamflow at streamgages unreachable during the flood, and updated flood-frequency estimates for streamgages in the Upper Yellowstone River Basin. Streamflows were the highest on record at 17 streamgages in the Upper Yellowstone River Basin. River stages were highest on record at most of those streamgages. The flood-related data and analyses are summarized in this fact sheet.

Idaho, Montana, Wyoming↗

Application of normalized radar backscatter and hyperspectral data to augment rangeland vegetation fractional classification

Rangeland ecosystems in the western United States are vulnerable to climate change, fire, and anthropogenic disturbances, yet classification of rangeland areas remains difficult due to frequently sparse vegetation canopies that increase the influence of soils and senesced vegetation, the overall abundance of senesced vegetation, heterogeneity of life forms, and limited ground-based data. The Rangeland Condition Monitoring Assessment and Projection (RCMAP) project provides fractional vegetation cover maps across western North America using Landsat imagery and artificial intelligence from 1985 to 2023 at yearly time-steps. The objectives of this case study are to apply hyperspectral data from several new data streams, including Sentinel Synthetic Aperture Radar (SAR) and Earth Surface Mineral Dust Source Investigation (EMIT), to the RCMAP model . We run a series of five tests (Landsat-base model, base + SAR, base + EMIT, base + SAR + EMIT, and base + Landsat NEXT [LNEXT] synthesized from EMIT) over a difficult-to-classify region centered in southwest Montana, USA. Our testing results indicate a clear accuracy benefit of adding SAR and EMIT data to the RCMAP model, with a 7.5% and 29% relative increase in independent accuracy ( R 2 ), respectively. The ability of SAR data to observe vegetation height allows for more accurate classification of vegetation types, whereas EMIT’s continuous characterization of the spectral response boosts discriminatory power relative to multispectral data. Our spectral profile analysis reveals the enhanced classification power with EMIT is related to both the improved spectral resolution and representation of the entire domain as compared to legacy Landsat. One key finding is that legacy Landsat bands largely miss portions of the electromagnetic spectrum where separation among important rangeland targets exists, namely in the 900–1250 nm and 1500–1780 nm range. Synthesized LNEXT data include these gaps, but the reduced spectral resolution compared to EMIT results in an intermediate 18% increase in accuracy relative to the base run. Here, we show the promise of enhanced classification accuracy using EMIT data, and to a smaller extent, SAR.

Idaho, Montana, Wyoming↗

Valuing wildlife sightings at the species-wide and individual animal levels: An approach and application to bear viewing in Yellowstone National Park

Wildlife viewing is a popular recreation activity in parks and protected areas around the world, yet information on the nonmarket value that visitors derive from viewing specific species, and how that value is affected by small-scale population changes, is lacking. We devised an approach to fill the gap. First, we applied the travel cost method to obtain the value of a wildlife viewing trip. Next, we estimated a value per animal sighting, which we aggregated by the number of visitors who experienced and valued a sighting. Finally, we incorporated the probability of viewing an animal to determine an individual animal’s contribution to sighting value, which varies by species type, visibility, and seasonality. We examined bear viewing in Yellowstone National Park. For grizzly and black bears, respectively, per-sighting values are $16 and $14, aggregate annual sighting values are $6.9 and $9.7 million, and annual per-bear viewing values are $46,000 and $15,000.

Idaho, Montana, Wyoming↗

A unified approach to long-term population monitoring of grizzly bears in the Greater Yellowstone Ecosystem

Long-term wildlife research and monitoring programs strive to maintain consistent data collections and analytical methods. Incorporating new techniques is important but can render data sets incongruent and limit their potential to discern trends in demographic parameters. Integrated population models (IPMs) can address these limitations by combining data sources that may span different periods into a unified statistical framework while providing a holistic view of population dynamics. We developed an IPM in a Bayesian framework for grizzly bears ( Ursus arctos ) in the Greater Yellowstone Ecosystem. We coupled demographic data with multiple, independent population count data to link annual changes in abundance with vital rates over 4 decades (1983–2023). Abundance increased threefold from an estimated 270 individuals in 1984 to 1030 individuals in 2023. Parameter estimates indicated survival of bears ≥2 years of age was high, contributing to robust population growth during the 1980s (λ = 1.023 [50 % interquartile range = 0.993–1.082]) and 1990s (λ = 1.064 [1.023–1.103]). A slowing of population growth started around 2000 (2000s: λ = 1.030 [0.989–1.068]) and continued into the 2010s (λ = 1.021 [0.985–1.057]), due primarily to reductions in survival of bears <2 years of age. These findings corroborate previous research that identified density-dependent effects as a likely cause. The IPM framework provided greater certainty and understanding regarding the dynamic demographic characteristics of the population and serves as a powerful monitoring tool for this long-lived species. Implementation of the IPM allows timely dissemination of demographic data to help inform adaptive management strategies and policy decisions necessary for the continued management and conservation of this population. This robust and flexible monitoring system allows scientists to investigate the effects of a changing ecosystem on population dynamics, incorporate new data sources and statistical models, and respond to changes in monitoring needs for the population. We highlight the efficacy of the IPM in estimating and tracking demographic parameters for a long-lived species, while accommodating shifts in monitoring techniques and data collections typical of long-term wildlife conservation programs worldwide.

Idaho, Montana, Wyoming↗