USGS Science⌕ Search

Geology topics

Research about Montana

Source-linked reports with geographic coverage including Montana.

At least 235 records · Page 13Linked to original sources

Distributions of small nongame fishes in the lower Yellowstone River

The Yellowstone River is the longest unimpounded river in the conterminous United States. It has a relatively natural flow regime, which helps maintain diverse habitats and fish assemblages uncommon in large rivers elsewhere. The lower Yellowstone River was thought to support a diverse nongame fish assemblage including several species of special concern. However, comprehensive data on the small nongame fish assemblage of the lower Yellowstone River is lacking. Therefore, we sampled the Yellowstone River downstream of its confluence with the Clark’s Fork using fyke nets and otter trawls to assess distributions and abundances of small nongame fishes. We captured 42 species (24 native and 18 nonnative) in the lower Yellowstone River with fyke nets. Native species constituted over 99% of the catch. Emerald shiners Notropis atherinoides , western silvery minnows Hybognathus argyritis , flathead chubs Platygobio gracilis , sand shiners Notropis stramineus , and longnose dace Rhinichthys cataractae composed nearly 94% of fyke net catch and were caught in every segment of the study area. We captured 24 species by otter trawling downstream of the Tongue River. Sturgeon chubs Macrhybopsis gelida , channel catfish Ictalurus punctatus , flathead chubs, stonecats Noturus flavus , and sicklefin chubs Macrhybopsis meeki composed 89% of the otter trawl catch. The upstream distributional limit of sturgeon chubs in the Yellowstone River was the Tongue River; few sicklefin chubs were captured above Intake Diversion Dam. This study not only provides biologists with baseline data for future monitoring efforts on the Yellowstone River but serves as a benchmark for management and conservation efforts in large rivers elsewhere as the Yellowstone River represents one of the best references for a naturally functioning Great Plains river.

Montana↗

Energetic costs of mange in wolves estimated from infrared thermography

Parasites, by definition, extract energy from their hosts and thus affect trophic and food web dynamics even when the parasite may have limited effects on host population size. We studied the energetic costs of mange ( Sarcoptes scabiei ) in wolves ( Canis lupus ) using thermal cameras to estimate heat losses associated with compromised insulation during the winter. We combined the field data of known, naturally infected wolves with data set on captive wolves with shaved patches of fur as a positive control to simulate mange-induced hair loss. We predict that during the winter in Montana, more severe mange infection increases heat loss by around 5.2 to 12 MJ per night (1240 to 2850 kcal, or a 65% to 78% increase) for small and large wolves, respectively accounting for wind effects. To maintain body temperature would require a significant proportion of a healthy wolf's total daily energy demands (18-22 MJ/day). We also predict how these thermal costs may increase in colder climates by comparing our predictions in Bozeman, Montana to those from a place with lower ambient temperatures (Fairbanks, Alaska). Contrary to our expectations, the 14°C differential between these regions was not as important as the potential differences in wind speed. These large increases in energetic demands can be mitigated by either increasing consumption rates or decreasing other energy demands. Data from GPS-collared wolves indicated that healthy wolves move, on average, 17 km per day, which was reduced by 1.5, 1.8 and 6.5 km for light, medium, and severe hair loss. In addition, the wolf with the most hair loss was less active at night and more active during the day, which is the converse of the movement patterns of healthy wolves. At the individual level mange infections create significant energy demands and altered behavioral patterns, this may have cascading effects on prey consumption rates, food web dynamics, predator-prey interactions, and scavenger communities.

Montana↗

Water-quality, bed-sediment, and biological data (October 2013 through September 2014) and statistical summaries of data for streams in the Clark Fork Basin, Montana

Water, bed sediment, and biota were sampled in streams from Butte to near Missoula, Montana, as part of a monitoring program in the upper Clark Fork Basin of western Montana. The sampling program was led by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, to characterize aquatic resources in the Clark Fork Basin, with emphasis on trace elements associated with historic mining and smelting activities. Sampling sites were located on the Clark Fork and selected tributaries. Water samples were collected periodically at 20 sites from October 2013 through September 2014. Bed-sediment and biota samples were collected once at 14 sites during August 2014. This report presents the analytical results and qualityassurance data for water-quality, bed-sediment, and biota samples collected at sites from October 2013 through September 2014. Water-quality data include concentrations of selected major ions, trace elements, and suspended sediment. At 12 sites, dissolved organic carbon and turbidity samples were collected. In addition, nitrogen (nitrate plus nitrite) samples were collected at two sites. Daily values of mean suspended-sediment concentration and suspended-sediment discharge were determined for four sites. Seasonal daily values of turbidity were determined for four sites. Bed-sediment data include trace-ele­ment concentrations in the fine-grained fraction. Biological data include trace-element concentrations in wholebody tissue of aquatic benthic insects. Statistical summaries of water-quality, bed-sediment, and biological data for sites in the upper Clark Fork Basin are provided for the period of record.

Montana↗

Integrating resource selection into spatial capture-recapture models for large carnivores

Wildlife managers need reliable methods to estimate large carnivore densities and population trends; yet large carnivores are elusive, difficult to detect, and occur at low densities making traditional approaches intractable. Recent advances in spatial capture-recapture (SCR) models have provided new approaches for monitoring trends in wildlife abundance and these methods are particularly applicable to large carnivores. We applied SCR models in a Bayesian framework to estimate mountain lion densities in the Bitterroot Mountains of west central Montana. We incorporate an existing resource selection function (RSF) as a density covariate to account for heterogeneity in habitat use across the study area and include data collected from harvested lions. We identify individuals through DNA samples collected by (1) biopsy darting mountain lions detected in systematic surveys of the study area, (2) opportunistically collecting hair and scat samples, and (3) sampling all harvested mountain lions. We included 80 DNA samples collected from 62 individuals in the analysis. Including information on predicted habitat use as a covariate on the distribution of activity centers reduced the median estimated density by 44%, the standard deviation by 7%, and the width of 95% credible intervals by 10% as compared to standard SCR models. Within the two management units of interest, we estimated a median mountain lion density of 4.5 mountain lions/100 km 2 (95% CI = 2.9, 7.7) and 5.2 mountain lions/100 km 2 (95% CI = 3.4, 9.1). Including harvested individuals (dead recovery) did not create a significant bias in the detection process by introducing individuals that could not be detected after removal. However, the dead recovery component of the model did have a substantial effect on results by increasing sample size. The ability to account for heterogeneity in habitat use provides a useful extension to SCR models, and will enhance the ability of wildlife managers to reliably and economically estimate density of wildlife populations, particularly large carnivores.

Montana↗

Spawning and rearing behavior of bull trout in a headwater lake ecosystem

Numerous life histories have been documented for bull trout Salvelinus confluentus . Lacustrine-adfluvial bull trout populations that occupy small, headwater lake ecosystems and migrate short distances to natal tributaries to spawn are likely common; however, much of the research on potamodromous bull trout has focused on describing the spawning and rearing characteristics of bull trout populations that occupy large rivers and lakes and make long distance spawning migrations to natal headwater streams. This study describes the spawning and rearing characteristics of lacustrine-adfluvial bull trout in the Quartz Lake drainage, Glacier National Park, USA, a small headwater lake ecosystem. Many spawning and rearing characteristics of bull trout in the Quartz Lake drainage are similar to potamodromous bull trout that migrate long distances. For example, subadult bull trout distribution was positively associated with slow-water habitat unit types and maximum wetted width, and negatively associated with increased stream gradient. Bull trout spawning also occurred when water temperatures were between 5 and 9 °C, and redds were generally located in stream segments with low stream gradient and abundant gravel and cobble substrates. However, this study also elucidated characteristics of bull trout biology that are not well documented in the literature, but may be relatively widespread and have important implications regarding general characteristics of bull trout ecology, use of available habitat by bull trout, and persistence of lacustrine-adfluvial bull trout in small headwater lake ecosystems.

Montana↗

Identifying priority chronic wasting disease surveillance areas for mule deer in Montana

Chronic wasting disease (CWD) is a fatal prion disease that affects a variety of ungulate species including mule deer ( Odocoileus hemionus ). As of 2014, no CWD cases had been reported in free-ranging ungulates in Montana. However, nearby cases in Canada, Wyoming, and the Dakotas indicated that the disease was encroaching on Montana's borders. Mule deer are native and common throughout Montana, and they represent a significant portion of the total hunter-harvested cervids in the state. The arrival of CWD in Montana may have significant ecosystem and socioeconomic impacts as well as potential consequences for wildlife management. We used 18,879 mule deer locations from 892 individual deer collected during 1975–2011 and modeled habitat selection for 7 herds in 5 of the 7 wildlife management regions in Montana. We estimated resource selection functions (RSF) in a Bayesian framework to predict summer and winter habitat preferences for mule deer. We estimated deer abundance from flyover counts for each region, and used the RSF predictions as weights to distribute the deer across the region. We then calculated the distance to the nearest known infected herds. We predicted areas of high risk of CWD infection in mule deer as areas with densities above the median density estimate and within the lowest quartile of distances to known infected herds. We identified these areas, the southeast corner of Montana and the north-central border near Alberta and Saskatchewan, as priority areas for CWD surveillance and management efforts.

Montana↗

Consequences of actively managing a small Bull Trout population in a fragmented landscape

Habitat fragmentation, which affects many native salmonid species, is one of the major factors contributing to the declines in distribution and abundance of Bull Trout Salvelinus confluentus . Increasingly, managers are considering options to maintain and enhance the persistence of isolated local populations through active management strategies. Understanding the ecological consequences of such actions is a necessary step in conservation planning. We used an individual-based model to evaluate the consequences of an ongoing management program aimed at mitigating the anthropogenic fragmentation of the lower Clark Fork River in Montana. Under this program juvenile Bull Trout are trapped and transported from small, headwater source populations to Lake Pend Oreille, Idaho, for rearing, and adults are subsequently recaptured in their upstream migration and returned to the natal population for spawning. We examined one of these populations and integrated empirical estimates of demographic parameters to simulate different management scenarios where moderate ( n = 4) and high ( n = 8) numbers of age-2, age-3, or age-4 Bull Trout were removed for transport with variable return rates under both demographic stochasticity and environmental perturbations. Our results indicated the risks from removal with no returns increased substantially when removal totals and age of Bull Trout removed from the simulated population increased. Specifically, removing eight age-3 or age-4 individuals resulted in 26% and 62% reductions in average adult population size, respectively, across simulations. We found the risks of transport were not likely alleviated with low (3%) or moderate (6%) return rates, and there were considerable risks of declines for the source population even when return rates were extremely high (>12%). Our simulations indicated little risk of declines for the source population with removals of age-2 Bull Trout, and any risks were alleviated with low return rates. However, we found higher return rates were particularly beneficial in the presence of large, density-independent perturbations.

Montana↗

Time scales of porphyry Cu deposit formation: insights from titanium diffusion in quartz

Porphyry dikes and hydrothermal veins from the porphyry Cu-Mo deposit at Butte, Montana, contain multiple generations of quartz that are distinct in scanning electron microscope-cathodoluminescence (SEM-CL) images and in Ti concentrations. A comparison of microprobe trace element profiles and maps to SEM-CL images shows that the concentration of Ti in quartz correlates positively with CL brightness but Al, K, and Fe do not. After calibrating CL brightness in relation to Ti concentration, we use the brightness gradient between different quartz generations as a proxy for Ti gradients that we model to determine time scales of quartz formation and cooling. Model results indicate that time scales of porphyry magma residence are ~1,000s of years and time scales from porphyry quartz phenocryst rim formation to porphyry dike injection and cooling are ~10s of years. Time scales for the formation and cooling of various generations of hydrothermal vein quartz range from 10s to 10,000s of years. These time scales are considerably shorter than the ~0.6 m.y. overall time frame for each porphyry-style mineralization pulse determined from isotopic studies at Butte, Montana. Simple heat conduction models provide a temporal reference point to compare chemical diffusion time scales, and we find that they support short dike and vein formation time scales. We interpret these relatively short time scales to indicate that the Butte porphyry deposit formed by short-lived episodes of hydrofracturing, dike injection, and vein formation, each with discrete thermal pulses, which repeated over the ~3 m.y. generation of the deposit.

Montana↗

Population genetic structure and disease in montane boreal toads: More heterozygous individuals are more likely to be infected with amphibian chytrid

Amphibians are more threatened than any other vertebrate group, with 41 % of species classified as threatened. The causes of most declines are not well understood, though many declines have been linked to disease. Additionally, amphibians are physiologically constrained to moist habitats and considered poor dispersers; thus, they may suffer genetic consequences of population isolation. To understand threats to the persistence of boreal toads ( Bufo boreas ) in Glacier National Park, USA, we genotyped 551 individuals at 11 microsatellite loci and used Bayesian clustering methods to describe population genetic structure and identify barriers to gene flow. We found evidence of two primary genetic groups that differed substantially in elevation and two secondary groups within the high elevation group. There was also evidence of further substructure within the southern high elevation group, suggesting mountain ridges are barriers to gene flow at local scales. Overall, genetic variation was high, but allelic richness declined with increasing elevation, reflecting greater isolation or smaller effective population sizes of high altitude populations. We tested for Batrachochytrium dendrobatidis (Bd), the fungal pathogen which causes chytridiomycosis, and we found that 35 of 199 toads were positive for Bd. Unexpectedly, more heterozygous individuals were more likely to be infected. This suggests that dispersal facilitates the spread of disease because heterozygosity may be highest where dispersal and gene flow are greatest.

Montana↗

Genomics and introgression: Discovery and mapping of thousands of species-diagnostic SNPs using RAD sequencing

Invasive hybridization and introgression pose a serious threat to the persistence of many native species. Understanding the effects of hybridization on native populations (e.g., fitness consequences) requires numerous species-diagnostic loci distributed genome-wide. Here we used RAD sequencing to discover thousands of single-nucleotide polymorphisms (SNPs) that are diagnostic between rainbow trout (RBT, Oncorhynchus mykiss ), the world’s most widely introduced fish, and native westslope cutthroat trout (WCT, O. clarkii lewisi ) in the northern Rocky Mountains, USA. We advanced previous work that identified 4,914 species-diagnostic loci by using longer sequence reads (100 bp vs . 60 bp) and a larger set of individuals ( n = 84). We sequenced RAD libraries for individuals from diverse sampling sources, including native populations of WCT and hatchery broodstocks of WCT and RBT. We also took advantage of a newly released reference genome assembly for RBT to align our RAD loci. In total, we discovered 16,788 putatively diagnostic SNPs, 10,267 of which we mapped to anchored chromosome locations on the RBT genome. A small portion of previously discovered putative diagnostic loci (325 of 4,914) were no longer diagnostic (i.e., fixed between species) based on our wider survey of non-hybridized RBT and WCT individuals. Our study suggests that RAD loci mapped to a draft genome assembly could provide the marker density required to identify genes and chromosomal regions influencing selection in admixed populations of conservation concern and evolutionary interest.

Montana↗

Long term changes in trout numbers following channel reconstruction, instream wood placement, and livestock removal from a spring creek in the Blackfoot Basin, Montana

To restore habitat for wild trout, Kleinschmidt Creek, a low-gradient, groundwater-dominated stream in the Blackfoot Basin, Montana, was reconstructed using natural channel design principles. Reconstruction increased stream sinuosity from a ratio of 1.1 to 1.6, decreased mean channel width from 14.5 to 2.8 m, and increased sediment transport capacity to reduce accumulations of fine instream sediment. To further improve trout habitat, coarse woody debris (CWD) was variably placed within the new channel and livestock were excluded to promote the vegetative recovery of the riparian area. To evaluate the response of wild trout (92% Brown Trout Salmo trutta ) to channel restoration, the abundance (number of trout per linear meter) and biomass (g/linear m) of age 1+ trout were monitored for 15 years (1998–2012) in a reach with low density CWD (1.3 stems/100 m) and compared with regional (reference) trends. Posttreatment (2002–2012) trout numbers in the low-density CWD reach were also compared with those in a reach with high-density CWD (18.2 stems/100 m). Long-term trends for the reference reaches showed a significant negative trend in trout abundance and no significant trend for biomass. Long-term trends for the low-density CWD reach showed a significant positive trend in abundance, as well as a significant trend in biomass. Trout abundance and biomass increased over the posttreatment period in the low-density CWD reach. However, in the high density CWD reach, while posttreatment abundance increased significantly, there was no significant trend in biomass. These results demonstrated that channel restoration increased wild trout populations in a deep, narrow, vegetated stream and that instream wood provided primarily short-term benefits during the early phase of habitat recovery.

Montana↗

Quantifying stream thermal regimes at management-pertinent scales: combining thermal infrared and stationary stream temperature data in a novel modeling framework.

Accurately quantifying stream thermal regimes can be challenging because stream temperatures are often spatially and temporally heterogeneous. In this study, we present a novel modeling framework that combines stream temperature data sets that are continuous in either space or time. Specifically, we merged the fine spatial resolution of thermal infrared (TIR) imagery with hourly data from 10 stationary temperature loggers in a 100 km portion of the Big Hole River, MT, USA. This combination allowed us to estimate summer thermal conditions at a relatively fine spatial resolution (every 100 m of stream length) over a large extent of stream (100 km of stream) during during the warmest part of the summer. Rigorous evaluation, including internal validation, external validation with spatially continuous instream temperature measurements collected from a Langrangian frame of reference, and sensitivity analyses, suggests the model was capable of accurately estimating longitudinal patterns in summer stream temperatures for this system Results revealed considerable spatial and temporal heterogeneity in summer stream temperatures and highlighted the value of assessing thermal regimes at relatively fine spatial and temporal scales. Preserving spatial and temporal variability and structure in abiotic stream data provides a critical foundation for understanding the dynamic, multiscale habitat needs of mobile stream organisms. Similarly, enhanced understanding of spatial and temporal variation in dynamic water quality attributes, including temporal sequence and spatial arrangement, can guide strategic placement of monitoring equipment that will subsequently capture variation in environmental conditions directly pertinent to research and management objectives.

Montana↗

Terrain parameters of glide snow avalanches and a simple spatial glide snow avalanche model

Glide snow avalanches are dangerous and difficult to predict. Despite substantial recent research there is still inadequate understanding regarding the controls of glide snow avalanche release. Glide snow avalanches often occur in similar terrain or the same locations annually, and repeat observations and prior work suggest that specific topography may be critical. Thus, to gain a better understanding of the terrain component of these types of avalanches we examined terrain parameters associated with the specific area of glide snow avalanche release in comparison to avalanche starting zones where no glide snow avalanches were observed (i.e. non-glide snow avalanche terrain). Glide snow avalanche occurrences visible from the Going-to-the-Sun Road corridor in Glacier National Park, Montana from 2003 to 2013 are investigated using a database of all avalanche occurrences derived of daily observations each year from 1 April to 1 June. This yielded 192 glide snow avalanches in 53 distinct avalanche paths. Each avalanche was digitized in a GIS using satellite, oblique, and aerial imagery as reference. A set of 117 non-glide snow avalanche starting zones were also selected in this manner. These were start zones with avalanche activity potential, but without glide avalanches observed. Topographical parameters such as area, slope, aspect, curvature, potential incoming solar radiation, distance from ridge, and elevation were then derived for the entire dataset utilizing tools with a GIS and a 10 m DEM. Ground class and a glide factor were calculated using a four level classification index with in-situ observations and a land surface type layer in a GIS. A total of 21 terrain variables were examined using a univariate analysis between areas where glide snow avalanches occurred and areas where glide snow avalanches were never observed, despite crack formation. Only two variables were not significantly different. The significantly different variables were then used to train a classification tree to distinguish between glide and non-glide snow avalanche terrain. A 10-fold cross validated tree resulted in four decision nodes to classify the data. The nodes split on glide factor, maximum slope angle, seasonal sum of incoming solar radiation, and maximum curvature to distinguish between glide snow avalanche and non-glide snow avalanche terrain with an unweighted average accuracy (RPC) of 0.95 and probability of detection of events (POD) of 0.99. Finally, the results of the cross-validated tree were used in a GIS to examine other areas, not used in the training dataset of the classification tree, of potential glide snow avalanche release within Glacier National Park. Using this understanding of the role of topographic parameters on glide snow avalanche activity, a spatial terrain based model was developed to identify other areas with high glide snow avalanche potential outside of the immediate observation area. This simple spatial model correctly classified 78 percent of actual glide snow avalanche terrain (pixel count) of a small test area of four independent observed glide snow avalanches.

Montana↗

Meteorological variables to aid forecasting deep slab avalanches on persistent weak layers

Deep slab avalanches are particularly challenging to forecast. These avalanches are difficult to trigger, yet when they release they tend to propagate far and can result in large and destructive avalanches. We utilized a 44-year record of avalanche control and meteorological data from Bridger Bowl ski area in southwest Montana to test the usefulness of meteorological variables for predicting seasons and days with deep slab avalanches. We defined deep slab avalanches as those that failed on persistent weak layers deeper than 0.9 m, and that occurred after February 1st. Previous studies often used meteorological variables from days prior to avalanches, but we also considered meteorological variables over the early months of the season. We used classification trees and random forests for our analyses. Our results showed seasons with either dry or wet deep slabs on persistent weak layers typically had less precipitation from November through January than seasons without deep slabs on persistent weak layers. Days with deep slab avalanches on persistent weak layers often had warmer minimum 24-hour air temperatures, and more precipitation over the prior seven days, than days without deep slabs on persistent weak layers. Days with deep wet slab avalanches on persistent weak layers were typically preceded by three days of above freezing air temperatures. Seasonal and daily meteorological variables were found useful to aid forecasting dry and wet deep slab avalanches on persistent weak layers, and should be used in combination with continuous observation of the snowpack and avalanche activity.

Montana↗

Climate-induced range contraction of a rare alpine aquatic invertebrate

Climate warming poses a serious threat to alpine-restricted species worldwide, yet few studies have empirically documented climate-induced changes in distributions. The rare stonefly, Zapada glacier (Baumann and Gaufin), endemic to alpine streams of Glacier National Park (GNP), Montana, was recently petitioned for listing under the US Endangered Species Act because of climate-change-induced glacier loss, yet little was known about its current status and distribution. We resampled streams throughout the historical distribution of Z. glacier to investigate trends in occurrence associated with changes in temperature and glacial extent. The current geographic distribution of the species was assessed using morphological characteristics of adults and DNA barcoding of nymphs. Bayesian phylogenetic analysis of mtDNA data revealed 8 distinct clades of the genus corresponding with 7 known species from GNP, and one potentially cryptic species. Climate model simulations indicate that average summer air temperature increased (0.67–1.00°C) during the study period (1960–2012), and glacial surface area decreased by ∼35% from 1966 to 2005. We detected Z. glacier in only 1 of the 6 historically occupied streams and at 2 new locations in GNP. These results suggest that an extremely restricted historical distribution of Z. glacier in GNP has been further reduced over the past several decades by an upstream retreat to higher, cooler sites as water temperatures increased and glacial masses decreased. More research is urgently needed to determine the status, distribution, and vulnerability of Z. glacier and other alpine stream invertebrates threatened by climate change in mountainous ecosystems.

Montana↗

Bird species turnover is related to changing predation risk along a vegetation gradient

Turnover in animal species along vegetation gradients is often assumed to reflect adaptive habitat preferences that are narrower than the full gradient. Specifically, animals may decline in abundance where their reproductive success is low, and these poor-quality locations differ among species. Yet habitat use does not always appear adaptive. The crucial tests of how abundances and demographic costs of animals vary along experimentally manipulated vegetation gradients are lacking. We examined habitat use and nest predation rates for 16 bird species that exhibited turnover with shifts in deciduous and coniferous vegetation. For most bird species, decreasing abundance was associated with increasing predation rates along both natural and experimentally modified vegetation gradients. This landscape-scale approach strongly supports the idea that vegetation-mediated effects of predation are associated with animal distributions and species turnover.

Montana↗

Modeling multi-scale resource selection for bear rubs in northwestern Montana

Both black ( Ursus americanus ) and grizzly bears ( U. arctos ) are known to rub on trees and other objects, producing a network of repeatedly used and identifiable rub sites. In 2012, we used a resource selection function to evaluate hypothesized relationships between locations of 887 bear rubs in northwestern Montana, USA, and elevation, slope angle, density of open roads and distance from areas of heightened plant-productivity likely containing forage for bears. Slope and density of open roads were negatively correlated with rub presence. No other covariates were supported as explanatory variables. We also hypothesized that bear rubs would be more strongly associated with closed roads and developed trails than with game trails. The frequencies of bear rubs on 30 paired segments of developed tracks and game trails were not different. Our results suggest bear rubs may be associated with bear travel routes, and support their use as “random” sampling devices for non-invasive spatial capture–recapture population monitoring.

Montana↗

Geochemical and petrographic data for intrusions peripheral to the Big Timber Stock, Crazy Mountains, Montana

The Paleocene Fort Union Formation hosts a compositionally diverse array of Eocene plugs, dikes, and sills arrayed around the Eocene Big Timber stock in the Crazy Mountains of south-central Montana. The geochemistry and petrography of the sills have not previously been characterized or interpreted. The purpose of this report is (1) to present available geochemical and petrographic data for several dozen samples of these rocks and (2) to provide a basic interpretive synthesis of these data.

Montana↗