Geology Reports⌕ Search

Geology topics

Research about Montana, North Dakota, South Dakota, Wyoming

Source-linked reports with geographic coverage including Montana, North Dakota, South Dakota, Wyoming.

18 recordsLinked to original sources

A three-dimensional geologic framework model of the northern Great Plains region of Montana, North Dakota, South Dakota, and Wyoming, USA

This report presents a new three-dimensional geologic framework model (GFM) of the northern Great Plains region, encompassing parts of Montana, North Dakota, South Dakota, and Wyoming. The model provides a regionally consistent, geographic information system (GIS)-ready representation of Phanerozoic sedimentary strata, major fault systems, and Precambrian basement geometry across two sedimentary basins and adjacent uplifts. More than 300,000 geologic and geophysical data inputs were synthesized to model 41 stratigraphic horizons and 47 faults, yielding an internally coherent, sealed-volume interpretation of the subsurface. The modeling workflow developed for this study demonstrates an efficient and scalable approach for constructing basin-to regional-scale GFMs in geologically complex and data-variable settings. Although model fidelity varies with data density and quality, the resulting geometry is broadly consistent with 1:500,000-scale geologic mapping and highlights areas where additional geologic study is most needed. The three-dimensional GFM provides a foundational framework to support groundwater, energy, and mineral resource assessments, and offers a transferable methodology for potential future U.S. Geological Survey efforts to build large-area subsurface models in underexplored regions of the United States.

Montana, North Dakota, South Dakota, Wyoming↗

Methods for peak-flow frequency analysis for streamgages in or near Montana, North Dakota, South Dakota, and Wyoming

The U.S. Geological Survey, in cooperation with the Montana Department of Natural Resources and Conservation, North Dakota Department of Water Resources, South Dakota Department of Transportation, and the Wyoming Water Development Office, has developed standard methods of peak-flow frequency analysis for studies in Montana, North Dakota, South Dakota, and Wyoming. These methods describe the implementation of national flood frequency guidelines described in Bulletin 17C ( https://doi.org/10.3133/tm4B5 ) for the four States and deviations from Bulletin 17C standard procedures to accommodate unusual hydrologic conditions. A U.S. Geological Survey data release accompanying this report ( https://doi.org/10.5066/P1WHRK8H ) provides example peak-flow frequency analyses for selected streamgages in the study area. The methods described in this report can be used to publish similar data releases for other streamgages in the study area.

Montana, North Dakota, South Dakota, Wyoming↗

Biophysical drivers for predicting the distribution and abundance of invasive yellow sweetclover in the Northern Great Plains

Context Yellow sweetclover ( Melilotus officinalis ; YSC) is an invasive biennial legume that bloomed across the Northern Great Plains in 2018–2019 in response to above-average precipitation. YSC can increase nitrogen (N) levels and potentially cause substantial changes in the composition of native plant species communities. There is little knowledge of the spatiotemporal variability and conditions causing substantial widespread blooms of YSC across western South Dakota (SD). Objectives We aimed to develop a generalized prediction model to predict the relative abundance of YSC in suitable habitats across rangelands of western South Dakota for 2019. Our research questions are: (1) What is the spatial extent of YSC across western South Dakota? (2) Which model can accurately predict the habitat and percent cover of YSC? and (3) What significant biophysical drivers affect its presence across western South Dakota? Methods We trained machine learning models with in situ data (2016–2021), Sentinel 2A-derived surface reflectance and indices (10 m, 20 m) and site-specific variables of climate, topography, and edaphic factors to optimize model performance. Results We identified moisture proxies (Shortwave Infrared reflectance and variability in Tasseled Cap Wetness) as the important predictors to explain the YSC presence. Land Surface Water Index and variability in summer temperature were the top predictors in explaining the YSC abundance. We demonstrated how machine learning algorithms could help generate valuable information on the spatial distribution of this invasive plant. We delineated major YSC hotspots in Butte, Pennington, and Corson Counties of South Dakota. The floodplains of major rivers, including White and Bad Rivers, and areas around Badlands National Park also showed a higher occurrence probability and cover percentage. Conclusions These prediction maps could aid land managers in devising management strategies for the regions that are prone to YSC outbreaks. The management workflow can also serve as a prototype for mapping other invasive plant species in similar regions.

Montana, North Dakota, South Dakota, Wyoming↗

Predictions and drivers of sub-reach-scale annual streamflow permanence for the upper Missouri River basin: 1989-2018

The presence of year-round surface water in streams (i.e., streamflow permanence) is an important factor for identifying aquatic habitat availability, determining the regulatory status of streams, managing land use change, allocating water resources, and designing scientific studies. However, accurate, high resolution, and dynamic prediction of streamflow permanence that accounts for year-to-year variability at a regional extent is a major gap in modeling capability. Herein, we expand and adapt the U.S. Geological Survey (USGS) PRObability of Streamflow PERmanence (PROSPER) model from its original implementation in the Pacific Northwest (PROSPER PNW ) to the upper Missouri River basin (PROSPER UM ), a geographical region that includes mountain and prairie ecosystems of the northern United States. PROSPER UM is an empirical model used to estimate the probability that a stream channel has year-round flow in response to climatic conditions (monthly and annual) and static physiographic predictor variables of the upstream basin. The structure and approach of PROSPER UM are generally consistent with the PROSPER PNW model but include improved spatial resolution (10 m) and a longer modeling period. Average model accuracy was 81 %. Drainage area, upstream proportion as wetlands, and upstream proportion as developed land cover were the most important predictor variables. The PROSPER UM model identifies decreases in streamflow permanence during climatically drier years, although there is variability in the magnitude across basins highlighting geographically varying sensitivity to drought. Variability in the response of perennial streams to drought conditions among basins in the study area was also observed.

Montana, North Dakota, South Dakota, Wyoming↗

Nesting, brood rearing, and summer habitat selection by translocated greater sage‐grouse in North Dakota, USA

Human enterprise has led to large‐scale changes in landscapes and altered wildlife population distribution and abundance, necessitating efficient and effective conservation strategies for impacted species. Greater sage‐grouse ( Centrocercus urophasianus ; hereafter sage‐grouse) are a widespread sagebrush ( Artemisia spp.) obligate species that has experienced population declines since the mid‐1900s resulting from habitat loss and expansion of anthropogenic features into sagebrush ecosystems. Habitat loss is especially evident in North Dakota, USA, on the northeastern fringe of sage‐grouse’ distribution, where a remnant population remains despite recent development of energy‐related infrastructure. Resource managers in this region have determined a need to augment sage‐grouse populations using translocation techniques that can be important management tools for countering species decline from range contraction. Although translocations are a common tool for wildlife management, very little research has evaluated habitat following translocation, to track individual behaviors such as habitat selection and fidelity to the release site, which can help inform habitat requirements to guide selection of future release sites. We provide an example where locations from previously released radio‐marked sage‐grouse are used in a resource selection function framework to evaluate habitat selection following translocation and identify areas of seasonal habitat to inform habitat management and potential restoration needs. We also evaluated possible changes in seasonal habitat since the late 1980s using spatial data provided by the Rangeland Analysis Platform coupled with resource selection modeling results. Our results serve as critical baseline information for habitat used by translocated individuals across life stages in this study area, and will inform future evaluations of population performance and potential for long‐term recovery.

Montana, North Dakota, South Dakota, Wyoming↗

Riverine turtles select habitats maintained by natural discharge regimes in an unimpounded large river

Turtle populations are imperiled worldwide, but limited ecological information from unaltered systems hampers science‐based management and conservation of some species, especially riverine turtles such as the spiny softshell ( Apalone spinifera ). We therefore investigated movements and spatial habitat selection of 54 A . spinifera in 633 river kilometres (rkm) of the least‐altered river in the conterminous United States—the Yellowstone River in Montana—from 2005 to 2009. Movement rates and home ranges were smaller than in fragmented, altered river systems because nesting and overwintering habitats were common and in close proximity. Habitat selection also differed. A . spinifera in the Yellowstone River overwintered in unaltered bluff pools and summered in complex reaches with side channels, islands, and diverse habitats. However, those in the highly altered Missouri River used deep alluvial pools for overwintering and flooded, inundated, or backwatered tributary mouths in spring and summer. Importantly, selected habitats in both rivers were functionally similar, including complex river reaches (with multiple channels, islands, and diverse habitats) and natural pool types. Unfortunately, these are the very habitats that are limited in rivers affected by dams, bank stabilization, and channelization. Therefore, preservation of natural and diverse riverine habitats—and the fluvial dynamics that maintain them—may enhance conservation of A . spinifera in large rivers.

Montana, North Dakota, South Dakota, Wyoming↗

Construction and calibration of a groundwater-flow model to assess groundwater availability in the uppermost principal aquifer systems of the Williston Basin, United States and Canada

The U.S. Geological Survey developed a groundwater-flow model for the uppermost principal aquifer systems in the Williston Basin in parts of Montana, North Dakota, and South Dakota in the United States and parts of Manitoba and Saskatchewan in Canada as part of a detailed assessment of the groundwater availability in the area. The assessment was done because of the potential for increased demands and stresses on groundwater associated with large-scale energy development in the area. As part of this assessment, a three-dimensional groundwater-flow model was developed as a tool that can be used to simulate how the groundwater-flow system responds to changes in hydrologic stresses at a regional scale. The three-dimensional groundwater-flow model was developed using the U.S. Geological Survey’s numerical finite-difference groundwater model with the Newton-Rhapson solver, MODFLOW–NWT, to represent the glacial, lower Tertiary, and Upper Cretaceous aquifer systems for steady-state (mean) hydrological conditions for 1981‒2005 and for transient (temporally varying) conditions using a combination of a steady-state period for pre-1960 and transient periods for 1961‒2005. The numerical model framework was constructed based on existing and interpreted hydrogeologic and geospatial data and consisted of eight layers. Two layers were used to represent the glacial aquifer system in the model; layer 1 represented the upper one-half and layer 2 represented the lower one-half of the glacial aquifer system. Three layers were used to represent the lower Tertiary aquifer system in the model; layer 3 represented the upper Fort Union aquifer, layer 4 represented the middle Fort Union hydrogeologic unit, and layer 5 represented the lower Fort Union aquifer. Three layers were used to represent the Upper Cretaceous aquifer system in the model; layer 6 represented the upper Hell Creek hydrogeologic unit, layer 7 represented the lower Hell Creek aquifer, and layer 8 represented the Fox Hills aquifer. The numerical model was constructed using a uniform grid with square cells that are about 1 mile (1,600 meters) on each side with a total of about 657,000 active cells. Model calibration was completed by linking Parameter ESTimation (PEST) software with MODFLOW–NWT. The PEST software uses statistical parameter estimation techniques to identify an optimum set of input parameters by adjusting individual model input parameters and assessing the differences, or residuals, between observed (measured or estimated) data and simulated values. Steady-state model calibration consisted of attempting to match mean simulated values to measured or estimated values of (1) hydraulic head, (2) hydraulic head differences between model layers, (3) stream infiltration, and (4) discharge to streams. Calibration of the transient model consisted of attempting to match simulated and measured temporally distributed values of hydraulic head changes, stream base flow, and groundwater discharge to artesian flowing wells. Hydraulic properties estimated through model calibration included hydraulic conductivity, vertical hydraulic conductivity, aquifer storage, and riverbed hydraulic conductivity in addition to groundwater recharge and well skin. The ability of the numerical model to accurately simulate groundwater flow in the Williston Basin was assessed primarily by its ability to match calibration targets for hydraulic head, stream base flow, and flowing well discharge. The steady-state model also was used to assess the simulated potentiometric surfaces in the upper Fort Union aquifer, the lower Fort Union aquifer, and the Fox Hills aquifer. Additionally, a previously estimated regional groundwater-flow budget was compared with the simulated steady-state groundwater-flow budget for the Williston Basin. The simulated potentiometric surfaces typically compared well with the estimated potentiometric surfaces based on measured hydraulic head data and indicated localized groundwater-flow gradients that were topographically controlled in outcrop areas and more generalized regional gradients where the aquifers were confined. The differences between the measured and simulated (residuals) hydraulic head values for 11,109 wells were assessed, which indicated that the steady-state model generally underestimated hydraulic head in the model area. This underestimation is indicated by a positive mean residual of 11.2 feet for all model layers. Layer 7, which represents the lower Hell Creek aquifer, is the only layer for which the steady-state model overestimated hydraulic head. Simulated groundwater-level changes for the transient model matched within plus or minus 2.5 feet of the measured values for more than 60 percent of all measurements and to within plus or minus 17.5 feet for 95 percent of all measurements; however, the transient model underestimated groundwater-level changes for all model layers. A comparison between simulated and estimated base flows for the steady-state and transient models indicated that both models overestimated base flow in streams and underestimated annual fluctuations in base flow. The estimated and simulated groundwater budgets indicate the model area received a substantial amount of recharge from precipitation and stream infiltration. The steady-state model indicated that reservoir seepage was a larger component of recharge in the Williston Basin than was previously estimated. Irrigation recharge and groundwater inflow from outside the Williston Basin accounted for a relatively small part of total groundwater recharge when compared with recharge from precipitation, stream infiltration, and reservoir seepage. Most of the estimated and simulated groundwater discharge in the Williston Basin was to streams and reservoirs. Simulated groundwater withdrawal, discharge to reservoirs, and groundwater outflow in the Williston Basin accounted for a smaller part of total groundwater discharge. The transient model was used to simulate discharge to 571 flowing artesian wells within the model area. Of the 571 established flowing artesian wells simulated by the model, 271 wells did not flow at any time during the simulation because hydraulic head was always below the land-surface altitude. As hydraulic head declined throughout the simulation, 68 of these wells responded by ceasing to flow by the end of 2005. Total mean simulated discharge for the 571 flowing artesian wells was 55.1 cubic feet per second (ft 3 /s), and the mean simulated flowing well discharge for individual wells was 0.118 ft 3 /s. Simulated discharge to individual flowing artesian wells increased from 0.039 to 0.177 ft 3 /s between 1961 and 1975 and decreased to 0.102 ft 3 /s by 2005. The mean residual for 34 flowing wells with measured discharge was 0.014 ft 3 /s, which indicates the transient model overestimated discharge to flowing artesian wells in the model area. Model limitations arise from aspects of the conceptual model and from simplifications inherent in the construction and calibration of a regional-scale numerical groundwater-flow model. Simplifying assumptions in defining hydraulic parameters in space and hydrologic stresses and time-varying observational data in time can limit the capabilities of this tool to simulate how the groundwater-flow system responds to changes in hydrologic stresses, particularly at the local scale; nevertheless, the steady-state model adequately simulated flow in the uppermost principal aquifer systems in the Williston Basin based on the comparison between the simulated and estimated groundwater-flow budget, the comparison between simulated and estimated potentiometric surfaces, and the results of the calibration process.

Montana, North Dakota, South Dakota, Wyoming↗

Decreased runoff response to precipitation, Little Missouri River Basin, northern Great Plains, USA

High variability in precipitation and streamflow in the semiarid northern Great Plains causes large uncertainty in water availability. This uncertainty is compounded by potential effects of future climate change. We examined historical variability in annual and growing season precipitation, temperature, and streamflow within the Little Missouri River Basin and identified differences in the runoff response to precipitation for the period 1976-2012 compared to 1939-1975 ( n = 37 years in both cases). Computed mean values for the second half of the record showed little change (<5%) in annual or growing season precipitation, but average annual runoff at the basin outlet decreased by 22%, with 66% of the reduction in flow occurring during the growing season. Our results show a statistically significant ( p < 0.10) 27% decrease in the annual runoff response to precipitation (runoff ratio). Surface-water withdrawals for various uses appear to account for <12% of the reduction in average annual flow volume, and we found no published or reported evidence of substantial flow reduction caused by groundwater pumping in this basin. Results of our analysis suggest that increases in monthly average maximum and minimum temperatures, including >1°C increases in January through March, are the dominant driver of the observed decrease in runoff response to precipitation in the Little Missouri River Basin.

Montana, North Dakota, South Dakota, Wyoming↗

Management of plains cottonwood at Theodore Roosevelt National Park, North Dakota

Establishment of cottonwood trees is driven by flood-induced channel migration, which provides the new surfaces necessary for successful germination and survival. Along the Little Missouri River the largest floods typically result from snowmelt in March or April. Seed release occurs in early summer, and seedlings usually germinate in moist, open locations on point bars at relatively low elevations above the channel. Subsequent channel migration allows seedlings to mature by protecting them from scour in floods and ice jams. Management actions that decrease channel movement will reduce cottonwood reproduction. Growth and survival of cottonwood trees are strongly decreased by extreme low flows. As a result, management activities that decrease low flows could strongly reduce growth or kill trees. Surface-flow diversions are less damaging to trees if carried out during the spring when flows are relatively high. Herbicide application by helicopter to control leafy spurge appears to have inadvertently damaged or killed about 25% of the cottonwood forest along the Little Missouri River in the South Unit. Area of adult trees sprayed has been reduced since 2007 to limit this damage. It is not known whether spraying of cottonwood seedlings on unforested point bars is reducing cottonwood reproduction in the South Unit. Warmer temperatures since 1976 have reduced flood peaks and the ice jamming that magnifies those peaks; as a result channel movement, cottonwood establishment and cottonwood growth have decreased. Increasing temperatures associated with global climate change could continue this trend.

Montana, North Dakota, South Dakota, Wyoming↗

Bakken, Three Forks largest continuous US oil accumulation

The recent reversal of the declining trend of US oil production is largely due to production from unconventional or "continuous" low-permeability reservoirs by use of multistage hydraulic fracturing of horizontal legs of exploration wells. The US currently produces about 7.4 million bo/d, and the increasing trend in domestic production has led to speculation that the US could become energy independent in oil in the near future. The US still imports an additional 11 million bo/d to meet consumption requirements).1 Critical to the discussion of energy independence are estimates of resources contained in low-porosity and low-permeability (or "tight") continuous oil and gas accumulations, which have formed a critical component of national energy policies in recent years.

Montana, North Dakota, South Dakota, Wyoming↗

Simulated responses of soil organic carbon stock to tillage management scenarios in the Northwest Great Plains

Background Tillage practices greatly affect carbon (C) stocks in agricultural soils. Quantification of the impacts of tillage on C stocks at a regional scale has been challenging because of the spatial heterogeneity of soil, climate, and management conditions. We evaluated the effects of tillage management on the dynamics of soil organic carbon (SOC) in croplands of the Northwest Great Plains ecoregion of the United States using the General Ensemble biogeochemical Modeling System (GEMS). Tillage management scenarios included actual tillage management (ATM), conventional tillage (CT), and no-till (NT). Results Model simulations show that the average amount of C (kg C ha -1 yr -1 ) released from croplands between 1972 and 2000 was 246 with ATM, 261 with CT, and 210 with NT. The reduction in the rate of C emissions with conversion of CT to NT at the ecoregion scale is much smaller than those reported at plot scale and simulated for other regions. Results indicate that the response of SOC to tillage practices depends significantly on baseline SOC levels: the conversion of CT to NT had less influence on SOC stocks in soils having lower baseline SOC levels but would lead to higher potentials to mitigate C release from soils having higher baseline SOC levels. Conclusion For assessing the potential of agricultural soils to mitigate C emissions with conservation tillage practices, it is critical to consider both the crop rotations being used at a local scale and the composition of all cropping systems at a regional scale.

Montana, North Dakota, South Dakota, Wyoming↗

Ground Water Atlas of the United States: Segment 8, Montana, North Dakota, South Dakota, Wyoming

The States of Montana, North Dakota, South Dakota, and Wyoming compose the 392,764-square-mile area of Segment 8, which is in the north-central part of the continental United States. The area varies topographically from the high rugged mountain ranges of the Rocky Mountains in western Montana and Wyoming to the gently undulating surface of the Central Lowland in eastern North Dakota and South Dakota (fig. 1). The Black Hills in southwestern South Dakota and northeastern Wyoming interrupt the uniformity of the intervening Great Plains. Segment 8 spans the Continental Divide, which is the drainage divide that separates streams that generally flow westward from those that generally flow eastward. The area of Segment 8 is drained by the following major rivers or river systems: the Green River drains southward to join the Colorado River, which ultimately discharges to the Gulf of California; the Clark Fork and the Kootenai Rivers drain generally westward by way of the Columbia River to discharge to the Pacific Ocean; the Missouri River system and the North Platte River drain eastward and southeastward to the Mississippi River, which discharges to the Gulf of Mexico; and the Red River of the North and the Souris River drain northward through Lake Winnipeg to ultimately discharge to Hudson Bay in Canada. These rivers and their tributaries are an important source of water for public-supply, domestic and commercial, agricultural, and industrial uses. Much of the surface water has long been appropriated for agricultural use, primarily irrigation, and for compliance with downstream water pacts. Reservoirs store some of the surface water for flood control, irrigation, power generation, and recreational purposes. Surface water is not always available when and where it is needed, and ground water is the only other source of supply. Ground water is obtained primarily from wells completed in unconsolidated-deposit aquifers that consist mostly of sand and gravel, and from wells completed in semi-consolidated- and consolidated-rock aquifers, chiefly sandstone and limestone. Some wells withdraw water from volcanic rocks, igneous and metamorphic rocks, or fractured fine-grained sedimentary rocks, such as shale; however, wells completed in these types of rocks generally yield only small volumes of water. Most wells in the four-State area of Segment 8 are on privately owned land (fig. 2). Agriculture, primarily irrigation, is one of the largest uses of ground water. The irrigation generally is on lowlands close to streams (fig. 3). Lowlands within a few miles of major streams usually are irrigated with surface water that is diverted by gravity flow from the main stream or a reservoir and transported through a canal system. Surface water also is pumped to irrigate land that gravity systems cannot supply. In addition, ground water is pumped from large-capacity wells to supplement surface water during times of drought or during seasons of the year when surface water is in short supply. Ground water is the only source of water for irrigation in much of the segment. The thickness and permeability of aquifers in the area of Segment 8 vary considerably, as do yields of wells completed in the aquifers. Ground-water levels and artesian pressures (hydraulic head) have declined significantly in some places as a result of excessive withdrawals by wells. State governments have taken steps to control the declines by enacting programs that either limit the number of additional wells that can be completed in a particular aquifer or prevent further ground-water development altogether. The demand for water is directly related to the distribution of people. In 1990, Montana had a population of 799,065; North Dakota, 638,800; South Dakota, 696,004; and Wyoming, 453,588. The more densely populated areas are on lowlands near major streams. Many of the mountain, desert, and upland areas lack major population centers, particularly in Montana and Wyoming, where use of much of the land is controlled by the Federal Government and withdrawal of ground water is restricted. Average annual precipitation (1951-80) in Segment 8 ranges from less than 8 inches in parts of Montana and Wyoming to more than 40 inches in some of the mountainous areas (fig. 4). Most storms move eastward through Segment 8 and are particularly common during the winter months. Moisture that evaporates from the Pacific Ocean is absorbed by eastward- moving air. As the moisture-laden air masses move eastward, they rise and cool as they encounter mountain ranges and lose some of their moisture to condensation. Consequently, the western sides of mountain ranges receive the most precipitation, much of it as snow during the winter months. In contrast, the eastern sides of some of the higher mountain ranges are in rain shadows and receive little precipitation. East of the Continental Divide, precipitation that falls during many summer storms results from northward-moving, moisture-laden air masses from the Gulf of Mexico. These air masses move northward when the polar front recedes; accordingly, a major part of the annual precipitation falls on the plains during the growing season. Average annual precipitation minus the total of average annual runoff plus evapotranspiration (the combination of evaporation and transpiration by plants) is the amount of water potentially available for recharge to the aquifers. Average annual runoff (1951-80) in the area of Segment 8 varies greatly, and the distribution of runoff (fig. 5) generally parallels that of precipitation. In arid and semiarid areas of the segment, most precipitation replenishes soil moisture, evaporates, or is transpired by vegetation, and only a small part of the precipitation is left to maintain streamflow or recharge aquifers. In wetter areas of the segment, much of the precipitation runs off the land surface directly to perennial streams. Because a smaller percentage of precipitation in wet areas usually is lost to evapotranspiration than in dry areas, more water is, therefore, available to recharge aquifers where more precipitation falls. Precipitation that falls as snow generally does not become runoff until spring thaws begin. Runoff is affected in some areas by reservoirs that have been constructed on major streams to mitigate flooding and to store water for irrigation, electrical power generation, and recreation. Water stored in reservoirs during times when runoff is great is subsequently released during drier periods to maintain downstream flow.

Montana, North Dakota, South Dakota, Wyoming↗

Geochemistry of ground-water in two sandstone aquifer systems in the Northern Great Plains in parts of Montana and Wyoming, North Dakota, and South Dakota

The Kootenai Formation in the Judith Basin, Montana, and the Lance Formation and Fox Hills Sandstone in the Powder River Basin, Wyoming, constitute two important sandstone aquifer systems in the Northern Great Plains region. Ground waters in each of these systems evolve from low dissolved-solids concentration, near-neutral pH, predominantly calcium and magnesium bicarbonate types in their recharge areas, to high dissolved-solids concentration, high pH, predominantly sodium-bicarbonate types in the basins. Oxidation potentials decrease as the waters flow downgradient under confined conditions. Calculation of the saturation states of aquifer minerals suggests several groups of mineral phases that could control ground- water chemistry. Mass transfer modeling indicates, however, that the observed behavior of major and minor dissolved species in both systems can satisfactorily be explained only by equilibration with cal cite, dolomite, or calcite and dolomite. The geochemistry of these systems is probably controlled by the incongruent dissolution of dolo mite to form calcite. This reaction appears to be driven by cation ex change and the dissolution of carbon dioxide. Plausible carbon diox ide sources include organic carbon oxidation and lignite coalification. Aluminosilicates influence major element chemistry primarily as sub trates for cation exchange, which, in combination with carbonate equilibria, buffer ground water pH at values of 8.5 to 8.9. Dissolved- iron concentrations are controlled by equilibration with amorphous ferric oxyhydroxides in oxidizing waters, with amorphous ferric oxyhydroxides and siderite in moderately reducing waters, and with siderite and amorphous ferrous sulfide in strongly reducing waters. Measured variations in dissolved carbonate isotopic composition compare favorably with carbon isotopic evolution, calculated by assuming dedolomitization. Recharge areas of the two systems are characterized by ground waters with high tritium and carbon-14 activities and relatively low dissolved-solids concentrations, with calcium and magnesium as the predominant cations. Recharge temperatures, calculated from dissolved-argon concentrations and 5 18 0 and SD isotopic measure ments, indicate that recharge is derived primarily from spring snow- melt rather than late spring and summer storms. Ground-water flow directions are generally parallel to trends of increasing dissolved solids concentrations and decreasing divalent to monovalent cation concentration ratios. However, these trends are sometimes obscured in areas of leakage or mixing. Further indication of leakage between aquifers is provided by abrupt changes in major element and isotopic chemistry, which are not characteristic of normally observed geo chemical evolution. Ground-water flow rates, calculated by adjusting measured carbon-14 activities for carbonate mass transfer, are com parable to values calculated from aquifer tests and potentiometric data. These carbon-14 flow rates average 1.6 meters per year for the Second Cat Creek sandstone of the Kootenai Formation, and 1.3 meters per year for the Lance-Fox Hills aquifer.

Montana, North Dakota, South Dakota, Wyoming↗

Maps showing formation temperatures and configurations of the tops of the Minnelusa Formation and the Madison Limestone, Powder River basin, Wyoming, Montana, and adjacent areas

This report is part of a study to describe the hydrologic framework needed to evaluate the water resources of the Paleozoic age aquifers in the Northern Great Plains coal region (fig 1). Preliminary studies by the U.S. Geological Survey and State agencies in Wyoming, Montana, and South Dakota have indicated that these aquifers might provide significant percentage of the water requirements for coal development. Data in this report are needed to help evaluate the potential of the Paleozoic age aquifers as a source for water supplies. The results will also be used to provide information for the orderly development of the aquifers.

Montana, North Dakota, South Dakota, Wyoming↗

Plan of study for the northern Great Plains regional aquifer-system analysis in parts of Montana, North Dakota, South Dakota, and Wyoming

The Northern Great Plains, an area of about 250,000 square miles in parts of Montana, North Dakota, South Dakota, and Wyoming, is underlain by an accumulation of sediments eroded from the Black Hills and from mountains to the west. Principal aquifers are areally extensive beds of sandstone within these sedimentary rocks, some at great depths. Anticipated future water needs dictate that available ground-water supplies be evaluated for management of this natural resource. The U.S. Geological Survey has started (1978) a 4-year study of the Northern Great Plains aquifer system. The objective of this study is to define availability and quality of ground water and to predict the effects of using this resource. To achieve this objective, the ground-water system will be described in terms of spatial distribution, hydraulics, geology, and geochemistry. Once described, the ground-water system will be simulated by mathematical models that will be used to define responses of the system to various management alternatives and assumed development patterns.

Montana, North Dakota, South Dakota, Wyoming↗

Stratigraphy and geologic history of the Montana group and equivalent rocks, Montana, Wyoming, and North and South Dakota

During Late Cretaceous time a broad north-trending epicontinental sea covered much of the western interior of North America and extended from the Gulf of Mexico to the Arctic Ocean. The sea was bounded on the west by a narrow, unstable, and constantly rising cordillera which extended from Central America to Alaska and which separated the sea from Pacific oceanic waters. The east margin of the sea was bounded by the low-lying stable platform of the central part of the United States. Rocks of the type Montana Group in Montana and equivalent rocks in adjacent States, which consist of eastward-pointing wedges of shallow-water marine and nonmarine strata that enclose westward-pointing wedges of fine-grained marine strata, were deposited in and marginal to this sea. These rocks range in age from middle Santonian to early Maestrichtian and represent a time span of about 14 million years. Twenty-nine distinctive ammonite zones, each with a time span of about half a million years, characterize the marine strata. Persistent beds of bentonite in the transgressive part of the Claggett and Bearpaw Shales of Montana and equivalent rocks elsewhere represent periods of explosive volcanism and perhaps concurrent subsidence along the west shore in the vicinity of the Elkhorn Mountains and the Deer Creek volcanic fields in Montana. Seaward retreat of st randlines, marked by deposition of the Telegraph Creek, Eagle, Judith River, and Fox Hills Formations in Montana and the Mesaverde Formation in Wyoming, may be attributed to uplift in near-coastal areas and to an increase in volcaniclastic rocks delivered to the sea. Rates of transgression and regression determined for the Montana Group in central Montana reveal that the strandline movement was more rapid during times of transgression. The regression of the Telegraph Creek and Eagle strandlines averaged about 50 miles per million years compared with a rate of about 95 miles per million years for the advance of the strand-line during Claggett time. The Judith River regression averaged about 60 miles per million years compared with movement of the strandline during the Bearpaw advance of about 70 miles per million years. The final retreat of marine waters from Montana, marked by the Fox Hills regression, was about 35 miles per million years at first, but near the end of the regression it accelerated to a rate of about 500 miles per million years. Rates of sedimentation range from less than 50 feet per million years in the eastern parts of North and South Dakota to at least 2,500 feet in western Wyoming. The low rates in the Dakotas correspond well with modern rates in the open ocean, and the rates in western Wyoming approach the rate of present coastal sedimentation.

Montana, North Dakota, South Dakota, Wyoming↗