USGS Science⌕ Search

SEARCH · USGS Science

Results for “Input”

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

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

At least 1,153 records · Page 64Linked to original sources
USGS2013

Simulation of salinity intrusion along the Georgia and South Carolina coasts using climate-change scenarios

Potential changes in climate could alter interactions between environmental and societal systems and adversely affect the availability of water resources in many coastal communities. Changes in streamflow patterns in conjunction with sea-level rise may change the salinity-intrusion dynamics of coastal rivers. Several municipal water-supply intakes are located along the Georgia and South Carolina coast that are proximal to the present day saltwater-freshwater interface of tidal rivers. Increases in the extent of salinity intrusion resulting from climate change could threaten the availability of freshwater supplies in the vicinity of these intakes. To effectively manage these supplies, water-resource managers need estimates of potential changes in the frequency, duration, and magnitude of salinity intrusion near their water-supply intakes that may occur as a result of climate change. This study examines potential effects of climate change, including altered streamflow and sea-level rise, on the dynamics of saltwater intrusion near municipal water-supply intakes in two coastal areas. One area consists of the Atlantic Intracoastal Waterway (AIW) and the Waccamaw River near Myrtle Beach along the Grand Strand of the South Carolina Coast, and the second area is on or near the lower Savannah River near Savannah, Georgia. The study evaluated how future sea-level rise and a reduction in streamflows can potentially affect salinity intrusion and threaten municipal water supplies and the biodiversity of freshwater tidal marshes in these two areas. Salinity intrusion occurs as a result of the interaction between three principal forces—streamflow, mean coastal water levels, and tidal range. To analyze and simulate salinity dynamics at critical coastal gaging stations near four municipal water-supply intakes, various data-mining techniques, including artificial neural network (ANN) models, were used to evaluate hourly streamflow, salinity, and coastal water-level data collected over a period exceeding 10 years. The ANN models were trained (calibrated) to learn the specific interactions that cause salinity intrusions, and resulting models were able to accurately simulate historical salinity dynamics in both study areas. Changes in sea level and streamflow quantity and timing can be simulated by the salinity intrusion models to evaluate various climate-change scenarios. The salinity intrusion models for the study areas are deployed in a decision support system to facilitate the use of the models for management decisions by coastal water-resource managers. The report describes the use of the salinity-intrusion models decision support system to evaluate salinity-intrusion dynamics for various climate-change scenarios, including incremental increases in sea level in combination with incremental decreases in streamflow. Operation of municipal water-treatment plants is problematic when the specific-conductance values for source water are greater than 1,000 to 2,000 microsiemens per centimeter (µS/cm). High specific-conductance values contribute to taste problems that require treatment. Data from a gage downstream from a municipal water intake indicate specific conductance exceeded 1,000 µS/cm about 5.4 percent of the time over the 14-year period from August 1995 to August 2008. Simulations of specific conductance at this gaging station that incorporates sea-level rises resulted in a doubling of the exceedances to 11.0 percent for a 1-foot increase and 17.6 percent for a 2-foot increase. The frequency of intrusion of water with specific conductance values of 1,000 µS/cm was less sensitive to incremental reductions in streamflow than to incremental increases in sea level. Simulations of conditions associated with a 10-percent reduction in streamflow, in combination with a 1-foot rise in sea level, increased the percentage of time specific conductance exceeded 1,000 µS/cm at this site from 11.0 to 13.3 percent, and a 20-percent reduction in streamflow increased the percentage of time to 16.6 percent. Precipitation and temperature data from a global circulation model were used, after scale adjustments, as input to a watershed model of the Yadkin-Pee Dee River basin, which flows into the Waccamaw River and Atlantic Intracoastal Waterway study area in South Carolina. The simulated streamflow for historical conditions and projected climate change in the future was used as input for the ANN model in decision support system. Results of simulations incorporating climate-change projections for alterations in streamflow indicate an increase in the frequency of salinity-intrusion events and a shift in the seasonal occurrence of the intrusion events from the summer to the fall.

South Carolina↗
USGS2013

Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

The Ogallala aquifer is an important water resource for the Rosebud Sioux Tribe in Gregory and Tripp Counties in south-central South Dakota and is used for irrigation, public supply, domestic, and stock water supplies. To better understand groundwater flow in the Ogallala aquifer, conceptual and numerical models of groundwater flow were developed for the aquifer. A conceptual model of the Ogallala aquifer was used to analyze groundwater flow and develop a numerical model to simulate groundwater flow in the aquifer. The MODFLOW–NWT model was used to simulate transient groundwater conditions for water years 1985–2009. The model was calibrated using statistical parameter estimation techniques. Potential future scenarios were simulated using the input parameters from the calibrated model for simulations of potential future drought and future increased pumping. Transient simulations were completed with the numerical model. A 200-year transient initialization period was used to establish starting conditions for the subsequent 25-year simulation of water years 1985–2009. The 25-year simulation was discretized into three seasonal stress periods per year and used to simulate transient conditions. A single-layer model was used to simulate flow and mass balance in the Ogallala aquifer with a grid of 133 rows and 282 columns and a uniform spacing of 500 meters (1,640 feet). Regional inflow and outflow were simulated along the western and southern boundaries using specified-head cells. All other boundaries were simulated using no-flow cells. Recharge to the aquifer occurs through precipitation on the outcrop area. Model calibration was accomplished using the Parameter Estimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. This program was designed to estimate parameter values that are statistically the most likely set of values to result in the smallest differences between simulated and observed values, within a given set of constraints. The potentiometric surface of the aquifer calculated during the 200-year initialization period established initial conditions for the transient simulation. Water levels for 38 observation wells were used to calibrate the 25-year simulation. Simulated hydraulic heads for the transient simulation were within plus or minus 20 feet of observed values for 95 percent of observation wells, and the mean absolute difference was 5.1 feet. Calibrated hydraulic conductivity ranged from 0.9 to 227 feet per day (ft/d). The annual recharge rates for the transient simulation (water years 1985–2009) ranged from 0.60 to 6.96 inches, with a mean of 3.68 inches for the Ogallala aquifer. This represents a mean recharge rate of 280.5 ft 3 /s for the model area. Discharge from the aquifer occurs through evapotranspiration, discharge to streams through river leakage and flow from springs and seeps, and well withdrawals. Water is withdrawn from wells for irrigation, public supply, domestic, and stock uses. Simulated mean discharge rates for water years 1985–2009 were about 185 cubic feet per second (ft 3 /s) for evapotranspiration, 66.7 ft 3 /s for discharge to streams, and 5.48 ft 3 /s for well withdrawals. Simulated annual evapotranspiration rates ranged from about 128 to 254 ft 3 /s, and outflow to streams ranged from 52.2 to 79.9 ft 3 /s. A sensitivity analysis was used to examine the response of the calibrated model to changes in model parameters for horizontal hydraulic conductivity, recharge, evapotranspiration, and spring and riverbed conductance. The model was most sensitive to recharge and maximum potential evapotranspiration and least sensitive to riverbed and spring conductances. Two potential future scenarios were simulated: a potential drought scenario and a potential increased pumping scenario. To simulate a potential drought scenario, a synthetic drought record was created, the mean of which was equal to 60 percent of the mean estimated recharge rate for the 25-year simulation period. Compared with the results of the calibrated model (non-drought simulation), the simulation representing a potential drought scenario resulted in water-level decreases of as much as 30 feet for the Ogallala aquifer. To simulate the effects of potential future increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 25-year simulation period. Compared with the results of the calibrated model, the simulation representing an increased pumping scenario resulted in water-level decreases of as much as 26 feet for the Ogallala aquifer. Groundwater budgets for the potential future scenario simulations were compared with the transient simulation representing water years 1985–2009. The simulation representing a potential drought scenario resulted in lower aquifer recharge from precipitation and decreased discharge from streams, springs, seeps, and evapotranspiration. The simulation representing a potential increased pumping scenario was similar to results from the transient simulation, with a slight increase in well withdrawals and a slight decrease in discharge from river leakage and evapotranspiration. This numerical model is suitable as a tool that could be used to better understand the flow system of the Ogallala aquifer, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the study area. The model also is useful to help assess the response of the aquifer to additional stresses, including potential drought conditions and increased well withdrawals.

South Dakota↗
USGS2013

Modeled future peak streamflows in four coastal Maine rivers

To safely and economically design bridges and culverts, it is necessary to compute the magnitude of peak streamflows that have specified annual exceedance probabilities (AEPs). Annual precipitation and air temperature in the northeastern United States are, in general, projected to increase during the 21st century. It is therefore important for engineers and resource managers to understand how peak flows may change in the future. This report, prepared in cooperation with the Maine Department of Transportation (MaineDOT), presents modeled changes in peak flows at four basins in coastal Maine on the basis of projected changes in air temperature and precipitation. To estimate future peak streamflows at the four basins in this study, historical values for climate (temperature and precipitation) in the basins were adjusted by different amounts and input to a hydrologic model of each study basin. To encompass the projected changes in climate in coastal Maine by the end of the 21st century, air temperatures were adjusted by four different amounts, from -3.6 degrees Fahrenheit (ºF) (-2 degrees Celsius (ºC)) to +10.8 ºF (+6 ºC) of observed temperatures. Precipitation was adjusted by three different percentage values from -15 percent to +30 percent of observed precipitation. The resulting 20 combinations of temperature and precipitation changes (includes the no-change scenarios) were input to Precipitation-Runoff Modeling System (PRMS) watershed models, and annual daily maximum peak flows were calculated for each combination. Modeled peak flows from the adjusted changes in temperature and precipitation were compared to unadjusted (historical) modeled peak flows. Annual daily maximum peak flows increase or decrease, depending on whether temperature or precipitation is adjusted; increases in air temperature (with no change in precipitation) lead to decreases in peak flows, whereas increases in precipitation (with no change in temperature) lead to increases in peak flows. As the magnitude of air temperatures increase in the four basins, peak flows decrease by larger amounts. If precipitation is held constant (no change from historical values), 17 to 26 percent decreases in peak flow occur at the four basins when temperature is increased by 7.2°F. If temperature is held constant, 26 to 38 percent increases in peak flow result from a 15-percent increase in precipitation. The largest decreases in peak flows at the four basins result from 15-percent decreases in precipitation combined with temperature increases of 10.8°F. The largest increases in peak flows generally result from 30-percent increases in precipitation combined with 3.6 °F decreases in temperatures. In many cases when temperature and precipitation both increase, small increases or decreases in annual daily maximum peak flows result. For likely changes projected for the northeastern United States for the middle of the 21st century (temperature increase of 3.6 °F and precipitation increases of 0 to 15 percent), peak-flow changes at the four coastal Maine basins in this study are modeled to be evenly distributed between increases and decreases of less than 25 percent. Peak flows with 50-percent and 1-percent AEPs (equivalent to 2-year and 100-year recurrence interval peak flows, respectively) were calculated for the four basins in the study using the PRMS-modeled annual daily maximum peak flows. Modeled peak flows with 50-percent and 1-percent AEPs with adjusted temperatures and precipitation were compared to unadjusted (historical) modeled values. Changes in peak flows with 50-percent AEPs are similar to changes in annual daily maximum peak flow; changes in peak flows with 1-percent AEPs are similar in pattern to changes in annual daily maximum peak flow, but some of the changes associated with increasing precipitation are much larger than changes in annual daily maximum peak flow. Substantial decreases in maximum annual winter snowpack water equivalent are modeled to occur with increasing air temperatures at the four basins in the study. (Snowpack is the snow on the ground that accumulates during a winter, and water equivalent is the amount of water in a snowpack if it were melted.) The decrease in modeled peak flows with increasing air temperature, given no change in precipitation amount, is likely caused by these decreases in winter snowpack and resulting decreases in snowmelt runoff. This Scientific Investigations Report, prepared in cooperation with the Maine Department of Transportation, presents a summary of modeled changes in peak flows at four basins in coastal Maine on the basis of projected changes in air temperature and precipitation. The full Fact Sheet (Hodgkins and Dudley, 2013) is available at http://pubs.usgs.gov/fs/2013/3021/.

Maine↗
USGS2016

Design and testing of a process-based groundwater vulnerability assessment (P-GWAVA) system for predicting concentrations of agrichemicals in groundwater across the United States

Efforts to assess the likelihood of groundwater contamination from surface-derived compounds have spanned more than three decades. Relatively few of these assessments, however, have involved the use of process-based simulations of contaminant transport and fate in the subsurface, or compared the predictions from such models with measured data—especially over regional to national scales. To address this need, a process-based groundwater vulnerability assessment (P-GWAVA) system was constructed to use transport-and-fate simulations to predict the concentration of any surface-derived compound at a specified depth in the vadose zone anywhere in the conterminous United States. The system was then used to simulate the concentrations of selected agrichemicals in the vadose zone beneath agricultural areas in multiple locations across the conterminous United States. The simulated concentrations were compared with measured concentrations of the compounds detected in shallow groundwater (that is, groundwater drawn from within a depth of 6.3 ± 0.5 meters [mean ± 95 percent confidence interval] below the water table) in more than 1,400 locations across the United States. The results from these comparisons were used to select the simulation approaches that led to the closest agreement between the simulated and the measured concentrations. The P-GWAVA system uses computer simulations that account for a broader range of the hydrologic, physical, biological and chemical phenomena known to control the transport and fate of solutes in the subsurface than has been accounted for by any other vulnerability assessment over regional to national scales. Such phenomena include preferential transport and the influences of temperature, soil properties, and depth on the partitioning, transport, and transformation of pesticides in the subsurface. Published methods and detailed soil property data are used to estimate a wide range of model input parameters for each site, including surface albedo, surface crust permeability, soil water content, Brooks-Corey parameters, saturated hydraulic conductivity, macroporosity and sizes of microbial populations, as well as solute partition coefficients, reaction rates, and meso-micropore diffusion rates. To ensure geographic consistency among the predictions, the only site-specific input data that are used are those that are available for all of the 48 conterminous states.

Scientific Investigations Report↗
USGS2015

Conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers, Pine Ridge Indian Reservation area, South Dakota, water years 1980-2009

The Ogallala and Arikaree aquifers are the largest sources of groundwater on the Pine Ridge Indian Reservation and are used extensively for irrigation and public and domestic water supplies. To assess the potential for decreased water levels and discharge to streams in the Pine Ridge Indian Reservation, conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers in southwestern South Dakota were developed by the U.S. Geological Survey in cooperation with the Oglala Sioux Tribe. The study area includes most of the Pine Ridge Reservation in Jackson and Shannon Counties and Indian trust lands in Bennett County in southwestern South Dakota. The High Plains aquifer, which includes the Ogallala and Arikaree aquifers, generally is less developed in South Dakota compared with other areas underlain by this aquifer; therefore, water levels in the High Plains aquifer in South Dakota generally fluctuated by less than 5 feet (ft) from 1980 to 1999. Despite minimal water-level changes in the High Plains aquifer in South Dakota, extensive withdrawals of groundwater for irrigation have caused water-level declines in many areas and increased concerns about the long-term sustainability of the aquifer; therefore, continued or increased withdrawals from the aquifer or prolonged drought may have the potential to affect water levels within the aquifer and discharge to important streams in the area. The Ogallala and Arikaree aquifers generally consist of poorly consolidated claystones, siltstones, sandstones, and shale deposited in fluvial and lacustrine environments. Saturated thicknesses ranged from 10 to 314 ft for the Ogllala aquifer and from 10 to 862 ft for the Arikaree aquifer. Previous hydraulic conductivity estimates ranged from less than 1 to 180 feet per day (ft/d) for the Ogallala aquifer and from less than 1 to 13 ft/d for the Arikaree aquifer. Recharge to the Ogallala and Arikaree aquifers is from precipitation on the outcrop areas, and discharge occurs through evapotranspiration, discharge to streams, and well withdrawals. Evapotranspiration generally occurs in topographically low areas along streams, and maximum evapotranspiration occurs when the water level is at the land surface. The generalized groundwater-flow direction is to the northeast with local flow towards streams. Precipitation for water years 1980–2009 ranged from about 11 to 39 inches per year (in/yr) and averaged about 19 in/yr. Estimated mean recharge for water years 1980–2009 was about 17.3 percent of precipitation for the Ogallala aquifer and 7.9 percent of precipitation for the Arikaree aquifer. The estimated mean maximum evapotranspiration for water years 1980–2009 was about 35 in/yr. Estimated mean base flow for gaged streams was about 0.06 cubic foot per second (ft 3 /s) per square mile of drainage area. Estimated mean total water use for water years 1980–2009 was 5.4 ft 3 /s from the Ogallala aquifer and 7.1 ft 3 /s from the Arikaree aquifer. A two-layer numerical groundwater-flow model was constructed using MODFLOW–NWT with a uniformly spaced grid consisting of 166 rows and 288 columns with cells 1,640 ft on a side. The numerical model of the Ogallala and Arikaree aquifers was used to simulate steady-state and transient conditions for water years 1980–2009. Model calibration was accomplished using the Parameter ESTimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. Aquifer boundaries were no-flow on the northern and western sides and constant-head on the southern and eastern sides. The mean arithmetic difference was 1.4 ft between the 731 simulated and observed hydraulic heads in the Ogallala aquifer and 9.8 ft between the 2,754 simulated and observed hydraulic heads in the Arikaree aquifer. Simulated mean discharge from the Ogallala and Arikaree aquifers to selected stream reaches was 92.1 ft 3 /s compared to estimated discharge of 88.7 ft 3 /s. Calibrated recharge for the transient simulation averaged 3.3 in/yr for the Ogallala aquifer and 1.1 in/yr for the Arikaree aquifer. The mean maximum potential evapotranspiration rate was 35.4 in/yr. Streambed conductance for perennial stream reaches averaged 530 feet squared per day. Horizontal hydraulic conductivity averaged 27 ft/d for the Ogallala aquifer and 1.0 ft/d for the Arikaree aquifer. The vertical hydraulic conductivity averaged 1.4 ft/d for the Ogallala aquifer and 0.004 ft/d for the Arikaree aquifer. Specific yield for the Ogallala aquifer was 0.15 (dimensionless) and averaged 0.02 for the Arikaree aquifer. Specific storage for the Arikaree aquifer was 1.7x10 -6 per foot. Simulated steady-state model inflow and outflow was 459 ft 3 /s. The percentages of inflows were 17 percent from constant-head boundaries, 9 percent from streams, and 74 percent from recharge. Percentages of outflow were 8 percent to constant-head boundaries, 1 percent to wells, 31 percent to streams, and 59 percent to evapotranspiration. Simulated net inflow from the Ogallala aquifer to the Arikaree aquifer ranged from about 22 ft 3 /s in dry years to about 37 ft 3 /s in wet years. Two hypothetical future stress scenarios were simulated using input from the 30-year calibrated simulation of water years 1980–2009. The first hypothetical scenario represented an increase in groundwater withdrawals from 50 hypothetical production wells completed in the Arikaree aquifer. At the end of the 30-year hypothetical increased pumping simulation, water levels declined as much as 66 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 26 percent (2.6 ft 3 /s) of increased withdrawals, and decreased evapotranspiration accounted for about 53 (5.3 ft 3 /s) percent of increased withdrawals. The second hypothetical scenario represented a 30-year period of decreased recharge (drought) by decreasing recharge 0.2 inch (24 ft 3 /s) for each water year. At the end of the hypothetical drought simulation, water levels declined as much as 10.9 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 23 percent (5.5 ft 3 /s) of decreased recharge, and decreased evapotranspiration accounted for about 72 percent (17.3 ft 3 /s) of decreased recharge. The numerical model is a tool that could be used to better understand the flow system of the Ogallala and Arikaree aquifers, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the Pine Ridge Reservation area in Bennett, Jackson, and Shannon Counties. The model also is useful to help assess the response of the aquifer to additional stress, including potential increased well withdrawals and potential drought conditions.

South Dakota↗
USGS2015

Climate change and prairie pothole wetlands: mitigating water-level and hydroperiod effects through upland management

Prairie pothole wetlands offer crucial habitat for North America’s waterfowl populations. The wetlands also support an abundance of other species and provide ecological services valued by society. The hydrology of prairie pothole wetlands is dependent on atmospheric interactions. Therefore, changes to the region’s climate can have profound effects on wetland hydrology. The relevant literature related to climate change and upland management effects on prairie pothole wetland water levels and hydroperiods was reviewed. Climate change is widely expected to affect water levels and hydroperiods of prairie pothole wetlands, as well as the biota and ecological services that the wetlands support. In general, hydrologic model projections that incorporate future climate change scenarios forecast lower water levels in prairie pothole wetlands and longer periods spent in a dry condition, despite potential increases in precipitation. However, the extreme natural variability in climate and hydrology of prairie pothole wetlands necessitates caution when interpreting model results. Recent changes in weather patterns throughout much of the Prairie Pothole Region have been in increased precipitation that results in increased water inputs to wetlands above losses associated with warmer temperatures. However, observed precipitation increases are within the range of natural climate variability and therefore, may not persist. Identifying management techniques with the potential to affect water inputs to prairie pothole wetlands would provide increased options for managers when dealing with the uncertainties associated with a changing climate. Several grassland management techniques (for example, grazing and burning) have the potential to affect water levels and hydroperiods of prairie pothole by affecting infiltration, evapotranspiration, and snow deposition.

Prairie Pothole Region↗
USGS2015

Estimating natural recharge in San Gorgonio Pass watersheds, California, 1913–2012

A daily precipitation-runoff model was developed to estimate spatially and temporally distributed recharge for groundwater basins in the San Gorgonio Pass area, southern California. The recharge estimates are needed to define transient boundary conditions for a groundwater-flow model being developed to evaluate the effects of pumping and climate on the long-term availability of groundwater. The area defined for estimating recharge is referred to as the San Gorgonio Pass watershed model (SGPWM) and includes three watersheds: San Timoteo Creek, Potrero Creek, and San Gorgonio River. The SGPWM was developed by using the U.S. Geological Survey INFILtration version 3.0 (INFILv3) model code used in previous studies of recharge in the southern California region, including the San Gorgonio Pass area. The SGPWM uses a 150-meter gridded discretization of the area of interest in order to account for spatial variability in climate and watershed characteristics. The high degree of spatial variability in climate and watershed characteristics in the San Gorgonio Pass area is caused, in part, by the high relief and rugged topography of the area. Daily climate data developed from a network of monitoring sites and published average monthly precipitation maps were used to develop the climate inputs for the SGPWM. Geographic Information System (GIS) data defining land surface altitude, vegetation, soils, surficial geology, and land cover were used to define input parameters representing the physical characteristics of the land surface, root zone, and shallow subsurface underlying the root zone. Model parameterization was based on a previous INFILv3 model developed for an area including the upper parts of the San Timoteo Creek and Potrero Creek drainages and the western part of the San Gorgonio River watershed. The previous INFILv3 model was calibrated by using available streamflow records from the model area. The SGPWM uses an updated INFILv3 version to represent shallow groundwater flow better beneath the root zone that contributes to lateral, downslope seepage rather than deep recharge. The SGPWM calibration was tested by using available streamflow records in the San Gorgonio Pass region. The SGPWM was used to simulate a 100-year water budget, including recharge and runoff, for water years 1913 through 2012. Results indicated that most recharge came from episodic infiltration of surface-water runoff in the larger stream channels. Results also indicated periods of great variability in recharge and runoff in response to variability in precipitation. More recharge was simulated for the area of the groundwater basin underlying the more permeable alluvial fill of the valley floor compared to recharge in the neighboring upland areas of the less permeable mountain blocks. The greater recharge was in response to the episodic streamflow that discharged from the mountain block areas and quickly infiltrated the permeable alluvial fill of the groundwater basin. Although precipitation at the higher altitudes of the mountain block was more than double precipitation at the lower altitudes of the valley floor, recharge for inter-channel areas of the mountain block was limited by the lower permeability bedrock underlying the thin soil cover, and most of the recharge in the mountain block was limited to the main stream channels underlain by alluvial fill.

California↗
USGS2016

Modern (1992–2011) and projected (2012–99) peak snowpack and May–July runoff for the Fort Peck Lake and Lake Sakakawea watersheds in the Upper Missouri River Basin

Mountain snowpack is an important contributor to runoff in the Upper Missouri River Basin; for example, high amounts of winter and spring precipitation in the mountains and plains in 2010–11 were associated with the peak runoff of record in 2011 in the Upper Missouri River Basin. To project trends in peak mountain snowpack and runoff in the upcoming decades, multiple linear regression models of peak mountain snowpack and total May–July runoff were developed for the Fort Peck Lake (above Fort Peck Dam) and lower Lake Sakakawea watersheds (between Fort Peck and Garrison Dams) in the Upper Missouri River Basin. Input to regression models included seasonal estimates of precipitation, air temperature, and total reference evapotranspiration stratified by elevation. Calibration was based on records from 107 weather stations from 1991 to 2011. Regressed annual peak mountain snowpack was used as input to the transfer function of May–July runoff. Peak snowpack and May–July runoff were projected for 2012–99 on the basis of air temperature and precipitation from the Community Climate System Model (CCSM) output. Two estimates of projected peak snowpack and May–July runoff for 2012–99 were computed: one estimate was based on output from the CCSM, version 3.0 (CCSM3), and the second estimate was based on output from the CCSM, version 4.0 (CCSM4). The significance of projected trends was based on the Kendall’s tau nonparametric test. Annual peak snowpack was projected to have a downward trend for the Fort Peck Lake watershed and an upward trend for the lower Lake Sakakawea watershed. Projections of May–July runoff had a significant downward trend for the Fort Peck Lake, lower Lake Sakakawea, and Lake Sakakawea (combination of Fort Peck Lake and lower Lake Sakakawea) watersheds. Downward trends in projected May–July runoff indicated that power production at Fort Peck Dam might be affected particularly in the later part of the simulation (2061–99); however, confidence in projected May–July runoff for the later part of the simulation was less certain because bias-corrected air temperatures from CCSM3 and CCSM4 commonly fell outside of the observed range used for calibration. Projected May–July runoff combined for the Fort Peck Lake and lower Lake Sakakawea watersheds were on the order of magnitude of the 2011 flood for 1 simulation year for each of the CCSM-based simulations. High peak snowpack and precipitation in April, May, and June in the plains was associated with large May–July runoff events; therefore, high precipitation at lower elevations in the Fort Peck Lake and lower Lake Sakakawea watersheds was a factor in the simulation of extreme runoff events at the magnitude of the 2011 flood.

Montana, North Dakota, Wyoming↗
USGS2016

Decision analysis to support development of the Glen Canyon Dam long-term experimental and management plan

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, National Park Service, and Argonne National Laboratory, completed a decision analysis to use in the evaluation of alternatives in the Environmental Impact Statement concerning the long-term management of water releases from Glen Canyon Dam and associated management activities. Two primary decision analysis methods, multicriteria decision analysis and the expected value of information, were used to evaluate the alternative strategies against the resource goals and to evaluate the influence of uncertainty. A total of 18 performance metrics associated with 8 out of 12 resource goals (fundamental objectives) were developed by the Bureau of Reclamation and National Park Service in partnership with subject-matter teams composed of Federal, State, tribal, and private experts. A total of 19 long-term strategies associated with 7 alternatives were developed by the Bureau of Reclamation, National Park Service, Argonne National Laboratory, U.S. Geological Survey, and Cooperating Agencies. The 19 long-term strategies were evaluated against the 18 performance metrics using a series of coupled simulation models, taking into account the effects of several important sources of uncertainty. A total of 27 Federal, State, tribal, and nongovernmental agencies were invited by the Assistant Secretary of Interior to participate in a swing-weighting exercise to understand the range of perspectives about how to place relative value on the resource goals and performance metrics; 14 of the 27 chose to participate. The results of the swing-weighting exercise were combined with the evaluation of the alternatives to complete a multicriteria decision analysis. The effects of uncertainty on the ranking of long-term strategies were evaluated through calculation of the value of information. The alternatives and their long-term strategies differed across performance metrics, producing unavoidable tradeoffs; thus, there was no long-term strategy that was dominated by another across all performance metrics. When the performance of each alternative was weighted across performance metrics, three alternatives (B, D, and G) were top-ranked depending on the set of weights proposed: Alternative B was favored by those stakeholders that placed a high value on hydropower; Alternative G was favored by those stakeholders that placed a high value on the restoration of natural processes, like beachbuilding and natural vegetation; and Alternative D was favored by the remaining stakeholders. Surprisingly, these rankings were not sensitive to the critical uncertainties that were evaluated; that is, the choice of a preferred long-term strategy was sensitive to the value-based judgment about how to place relative weight on the resource goals but was not sensitive to the uncertainties in the system dynamics that were evaluated in this analysis. The one area of uncertainty that did slightly affect the ranking of alternatives was the long-term pattern of hydrological input; because of this sensitivity, some attention to the possible effects of climate change is warranted. The results of the decision analysis are meant to serve as only one of many sources of information that can be used to evaluate the alternatives proposed in the Environmental Impact Statement. These results only focus on those resource goals for which quantitative performance metrics could be formulated and evaluated; there are other important aspects of the resource goals that also need to be considered. Not all the stakeholders who were invited to participate in the decision analysis chose to do so; thus, the Bureau of Reclamation, National Park Service, and U.S. Department of Interior may want to consider other input.

Arizona↗
USGS2016

Estimating spatially and temporally varying recharge and runoff from precipitation and urban irrigation in the Los Angeles Basin, California

A daily precipitation-runoff model, referred to as the Los Angeles Basin watershed model (LABWM), was used to estimate recharge and runoff for a 5,047 square kilometer study area that included the greater Los Angeles area and all surface-water drainages potentially contributing recharge to a 1,450 square kilometer groundwater-study area underlying the greater Los Angeles area, referred to as the Los Angeles groundwater-study area. The recharge estimates for the Los Angeles groundwater-study area included spatially distributed recharge in response to the infiltration of precipitation, runoff, and urban irrigation, as well as mountain-front recharge from surface-water drainages bordering the groundwater-study area. The recharge and runoff estimates incorporated a new method for estimating urban irrigation, consisting of residential and commercial landscape watering, based on land use and the percentage of pervious land area. The LABWM used a 201.17-meter gridded discretization of the study area to represent spatially distributed climate and watershed characteristics affecting the surface and shallow sub-surface hydrology for the Los Angeles groundwater study area. Climate data from a local network of 201 monitoring sites and published maps of 30-year-average monthly precipitation and maximum and minimum air temperature were used to develop the climate inputs for the LABWM. Published maps of land use, land cover, soils, vegetation, and surficial geology were used to represent the physical characteristics of the LABWM area. The LABWM was calibrated to available streamflow records at six streamflow-gaging stations. Model results for a 100-year target-simulation period, from water years 1915 through 2014, were used to quantify and evaluate the spatial and temporal variability of water-budget components, including evapotranspiration (ET), recharge, and runoff. The largest outflow of water from the LABWM was ET; the 100-year average ET rate of 362 millimeters per year (mm/yr) accounted for 66 percent of the combined water inflow of 551 mm/yr, including 488 mm/yr from precipitation and 63 mm/yr from urban irrigation. The simulated ET rate varied from a minimum of 0 mm/yr for impervious areas to high values of more than 1,000 mm/yr for many areas, including the south-facing slopes of the San Gabriel Mountains, stream channels underlain by permeable soils and thick root zones, and pervious locations receiving inflows both from urban irrigation and surface water. Runoff was the next largest outflow, averaging 145 mm/yr for the 100-year period, or 26 percent of the combined precipitation and urban-irrigation inflow. Recharge averaged 45 mm/yr, or about 8 percent of the combined inflow from precipitation and urban irrigation. Simulation results indicated that recharge in response to urban irrigation was an important component of spatially distributed recharge, contributing an average of 56 percent of the total recharge to the eight LABWM subdomains containing the Los Angeles groundwater study area. The 100‑year average recharge rate for the eight subdomains was 41 mm/yr, or 8,473 hectare-meters per year (ha-m/yr), with urban irrigation included in the simulation compared to a recharge rate of 18 mm/yr, or 3,741 ha-m/yr, with urban irrigation excluded. In contrast to recharge, the effect of urban irrigation on runoff was slight; runoff was 72,667 ha-m/yr with urban irrigation included compared to 72,618 ha-m/yr with urban irrigation excluded, an increase of only 48 ha-m/yr (about 0.1 percent). Simulation results also indicated that potential recharge from hilly drainages outside of, but bordering and tributary to, the lower-lying area of the Los Angeles groundwater study area, in this study referred to as mountain-front recharge, could provide an important contribution to the total recharge for the groundwater basins. The time-averaged recharge rate was similar to the combined direct and mountain-front recharge components estimated in a previous study and used as input for a calibrated groundwater model. The annual (water year) recharge estimates simulated in this study, however, indicated much greater year-to-year variability, which was dependent on year-to-year variability in the magnitude and distribution of daily precipitation, compared to the previous estimates.

California↗
USGS2016

Generalized sediment budgets of the Lower Missouri River, 1968–2014

Sediment budgets of the Lower Missouri River were developed in a study led by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. The scope of the study included the development of a long-term (post-impoundment, 1968–2014) average annual sediment budget and selected annual, monthly, and daily sediment budgets for a reach and period that adequate data were available. Included in the analyses were 31 main-stem and tributary stations of the Lower Missouri River and two Mississippi River stations—the Mississippi River below Grafton, Illinois, and the Mississippi River at St. Louis, Missouri. Long-term average annual suspended-sediment loads of Missouri River main-stem stations ranged from 0.33 million tons at the Missouri River at Yankton, South Dakota, station to 71.2 million tons at Missouri River at Hermann, Mo., station. Gaged tributary gains accounted for 9–36 percent of the local reach budgets and cumulative gaged tributary contributions accounted for 84 percent of the long-term average suspended-sediment load of the Missouri River at Hermann, Mo., station. Although the sediment budgets for seven defined main-stem reaches generally were incomplete—missing bedload, reach storage, and ungaged tributary contributions—the budget residuals (net result of sediment inputs and outputs) for six of the seven reaches ranged from -7.0 to 1.7 million tons, or from -9.2 to 4.0 percent of the reach output suspended-sediment load, and were within the 10 percent reported measurement error of annual suspended-sediment loads for large rivers. The remaining reach, downstream from Gavin’s Point Dam, extended from Yankton, S. Dak., to Sioux City, Iowa, and had a budget residual of -9.8 million tons, which was -88 percent of the suspended-sediment load at Sioux City. The Lower Missouri River reach from Omaha, Nebraska, to Nebraska City, Nebr., had periods of concurrent sediment data for each primary budget component with which to analyze and determine a suspended-sediment budget for selected annual, monthly, and daily time increments. The temporal changes in the cumulative annual budget residuals were poorly correlated with the comparatively steady 1968–2011 annual stage trends at the Missouri River at Nebraska City, Nebr., station. An accurate total sediment budget is developed by having concurrent data available for all primary suspended and bedload components for a reach of interest throughout a period. Such a complete budget, with concurrent record for suspended-sediment load and bedload components, is unavailable for any reach and period in the Lower Missouri River. The primary data gaps are in bedload data, and also in suspended-sediment gains and losses including ungaged tributary inputs and sediment storage. Bedload data gaps in the Missouri River Basin are much more prevalent than suspended-sediment data gaps, and the first step in the development of reach bedload budgets is the establishment of a standardized bedload monitoring program at main-stem stations. The temporal changes in flow-adjusted suspended-sediment concentrations analyzed at main-stem Missouri River stations indicated an overall downward change in concentrations between 1968 and 2014. Temporary declines in flow-adjusted suspended-sediment concentrations during and following large floods were evident but generally returned to near pre-flood values within about 6 months. Data uncertainties associated with the development of a sediment budget include uncertainties associated with the collection of suspended-sediment and bedload data and the computation of suspended-sediment loads. These uncertainties vary depending on the frequency of data collection, the variability of conditions being represented by the discrete samples, and the statistical approach to suspended-sediment load computations. The coefficients of variation of suspended-sediment loads of Missouri River tributary stations for 1968–2014 were greater, 75.0 percent, than the main-stem stations, 47.1 percent. The lower coefficient of variation at main-stem stations compared to tributaries, primarily is the result of the lower variability in streamflow and sediment discharge identified at main-stem stations. To obtain similar accuracy between suspended-sediment loads at main-stem and tributary stations, a longer period of record is required of the tributary stations. During 1968–2014, however, the Missouri River main-stem station record was much more complete (87 percent) than the tributary station record (28 percent).

Missouri River Basin↗
USGS2017

Geomorphic response of the North Fork Stillaguamish River to the State Route 530 landslide near Oso, Washington

On March 22, 2014, the State Route 530 Landslide near Oso, Washington mobilized 8 million cubic meters of unconsolidated Pleistocene material, creating a valley‑spanning deposit that fully impounded the North Fork Stillaguamish River. The river overtopped the 8-meter high debris impoundment within 25 hours and began steadily incising a new channel through the center of the deposit. Repeat topographic surveys, sediment transport measurements, bedload transport models, and observations of downstream channel change were used to document the establishment of that new channel through the landslide and assess the potential for downstream aggradation or channel change that might increase downstream flood hazards. Efficient erosion of the landslide deposit, associated with the steep knickzone formed by the downstream edge of the deposit, resulted in the re-establishment of a 20–40 meters wide, deeply inset channel through the entire deposit by May 2014, 2 months after the landslide. The mean water-surface elevation of the channel through the landslide decreased 7 meters during that 2-month period, and was about 1 meter above the pre-landslide profile in July 2014. The 2014–15 flood season, which included flows near the 0.5 annual exceedance probability discharge (2-year flood), widened the channel tens of meters, and further lowered the water-surface profile 0.5 meter. The planform position evolved slowly as a result of 5–20-meter high banks predominantly composed of clay-rich, cohesive lacustrine material. Erosion of the landslide deposit delivered a total of 820 thousand metric tons of sediment to the North Fork Stillaguamish River over the 18 months following the landslide. The sediment delivery from the deposit was predominantly fine grained: 77 percent (by mass) of the eroded material was silt or clay (less than 0.063 millimeter [mm]), 19 percent sand (0.063–2 mm), and 4 percent pebbles and cobbles (greater than 2 mm). Over the 18 months following the landslide, the bedload at a site 5 kilometers downstream of the landslide was estimated to be 310±65 thousand metric tons, and the suspended load at that same site was estimated to be 990±110 thousand metric tons. These loads represent the combined input from the landslide and ambient upstream sources; over the study interval, landslide sediment made up about 20–40 percent of the bedload, and 65–85 percent of the suspended-sediment load at this site. At a site 70 kilometers downstream of the landslide, near the mouth of the main‑stem Stillaguamish River, suspended sediment loads were estimated to be about 1,440 thousand metric tons, of which about 600 thousand metric tons, or 30 percent, likely was derived from the landslide. The mass of landslide sediment in suspension at the mouth of the river, and the timing of arrival of that sediment, indicates that about 70 percent of the landslide sediment eroded during the study period was quickly transported through the entire basin, exiting into Puget Sound within weeks of initial entrainment. Empirical bedload transport equations, in conjunction with surficial grain-size data and output from a one‑dimensional hydraulic model, were used to estimate spatial trends in bedload transport capacity, highlighting areas where reach-scale conditions would be most likely to promote deposition of coarse landslide sediment. Transport capacities decreased sharply over a reach about 5 kilometers downstream of the landslide and remained relatively low over the next 10 kilometers downstream. However, the magnitude of calculated transport capacities are large relative to the coarse sediment input from the landslide, suggesting that substantial deposition of landslide sediment was not likely to occur. These assessments were corroborated by observations of channel change, which indicated that the downstream channel response to the landslide was modest and short-lived. The most pronounced downstream effects included a wedge of aggradation just downstream of the landslide, about 1 meter high and extending a kilometer downstream, and a 0.3-meter pulse of aggradation observed 5 kilometers downstream of the landslide. In both locations, peak aggradation and channel response occurred within about a month of the landslide, and both sites had largely recovered to pre-landslide conditions by July 2014. No substantial channel change clearly linked to the landslide was observed after July 2014 except for a modest fining of surficial gravel size distributions and continued recovery and incision of the reach just downstream of the landslide. The muted downstream response of the North Fork Stillaguamish River to the State Route 530 Landslide primarily can be attributed to the cohesive, silt- and clay-rich material that bounded most of the new channel. Although the river efficiently incised a new channel through the deposit, subsequent rates of lateral erosion were slowed by the tall, cohesive banks, limiting the total volume of sediment delivery. Once entrained, however, most landslide material was rapidly transported downstream in suspension with little geomorphic effect. Landslide material coarse enough to travel as bedload was predominantly sand and fine gravel, and sediment transport models and observations of downstream change indicated that the rate of coarse sediment delivery from the landslide did not exceed the rivers ability to transport that material. The generally muted downstream response to sediment delivery from the State Route 530 Landslide, as well as the mechanics of that delivery and response, were generally consistent with observations made following the intentional removal of constructed dams. The rate and efficiency of erosion from the landslide decreased over the period of analysis, as the new channel approached a quasi-equilibrium form. In the absence of additional hillslope activity, rates of erosion from the landslide are likely to be small compared to those over the first 18 months after the landslide. The modest channel response to the highest rates of sediment delivery, and rapid recovery thereafter, indicate that the river should be able to convey the continued supply of landslide-derived sediment effectively with little effect on the downstream morphology and flood risks.

Washington↗
USGS2017

Estimation of the groundwater resources of the bedrock aquifers at the Kettle Moraine Springs State Fish Hatchery, Sheboygan County, Wisconsin

Groundwater resources information was needed to understand regional aquifer systems and water available to wells and springs for rearing important Lake Michigan fish species at the Kettle Moraine Springs State Fish Hatchery in Sheboygan County, Wisconsin. As a basis for estimating the groundwater resources available, an existing groundwater-flow model was refined, and new groundwater-flow models were developed for the Kettle Moraine Springs State Fish Hatchery area using the U.S. Geological Survey (USGS) finite-difference code MODFLOW. This report describes the origin and construction of these groundwater-flow models and their use in testing conceptual models and simulating the hydrogeologic system. The study area is in the Eastern Ridges and Lowlands geographical province of Wisconsin, and the hatchery property is situated on the southeastern edge of the Kettle Moraine, a north-south trending topographic high of glacial origin. The bedrock units underlying the study area consist of Cambrian, Ordovician, and Silurian units of carbonate and siliciclastic lithology. In the Sheboygan County area, the sedimentary bedrock sequence reaches a thickness of as much as about 1,600 feet (ft). Two aquifer systems are present at the Kettle Moraine Springs State Fish Hatchery. A shallow system is made up of Silurian bedrock, consisting chiefly of dolomite, overlain by unconsolidated Quaternary-age glacial deposits. The glacial deposits of this aquifer system are the typical source of water to local springs, including the springs that have historically supplied the hatchery. The shallow aquifer system, therefore, consists of the unconsolidated glacial aquifer and the underlying bedrock Silurian aquifer. Most residential wells in the area draw from the Silurian aquifer. A deeper confined aquifer system is made up of Cambrian- and Ordovician-age bedrock units including sandstone formations. Because of its depth, very few wells are completed in the Cambrian-Ordovician aquifer system (COAS) near the Kettle Moraine Springs State Fish Hatchery. Three groundwater-flow models were used to estimate the water resources available to the hatchery from bedrock aquifers under selected scenarios of well placement and seasonal water requirements and subject to constraints on the effects of pumping on neighboring wells, local springs, and creeks. Model input data (recharge, water withdrawal, and boundary conditions) for these models were compiled from a number of data and information sources. The first model, named the “KMS model,” (KMS stands for Kettle Moraine Springs) is an inset model derived from a published USGS regional Lake Michigan Basin model and was constructed to simulate groundwater pumping from the semiconfined Silurian aquifer. The second model, named the “Pumping Test model,” was constructed to evaluate an aquifer pumping test conducted in the COAS as part of this project. The Pumping Test model was also used to simulate the local effects of 20 years of groundwater pumping from this deep bedrock aquifer for future hatchery operations. The third model, named the “LMB modified model,” is a version of the published Lake Michigan Basin (LMB) model that was modified with aquifer parameters refined in an area around the hatchery (approximately a 5-mile radius circle, corresponding to the area stressed by the aquifer pumping test). This LMB modified model was applied to evaluate regional effects of pumping from the confined COAS. The available Silurian aquifer groundwater resource was estimated using the KMS model with three scenarios—named “AllConstraints,” “Constraints2,” and “Constraints3”—that specified local water-level and flow constraints such as drawdown at nearby household wells, water levels inside pumping well boreholes, and flow in local streams and springs. Each scenario utilized the MODFLOW Groundwater Management Process (GWM) to select three locations from six candidate locations that provided the greatest combined flow while satisfying the constraints. The three constraint scenarios provided estimates of 430 gallons per minute (gal/min), 480 gal/min, and 520 gal/min pumping from three wells—AllConstraints, Constraints2, and Constraints3, respectively. The same three wells were selected for the scenarios that estimated 480 gal/min and 520 gal/min; the scenario that estimated 430 gal/min shared two of these same wells, but the third selected well was different. The available COAS groundwater resource was estimated by two scenarios with each conducted over a period of 20 years with the Pumping Test model and the LMB modified model . The Pumping Test model was used to simulate local effects of pumping, and the LMB modified model was used to simulate regional effects of pumping. The scenarios simulate a range of total and seasonal pumping rates potentially linked to site activities. Scenario 1 simulates two wells completed in the Cambrian-Ordovician aquifer system, each pumping for 8 months at 300 gal/min, followed by pumping for 4 months at 600 gal/min. The average yearly pumping rate of Scenario 1 is 800 gal/min. Scenario 2 simulates three wells completed in the Cambrian-Ordovician aquifer system pumping for 8 months at 200 gal/min, followed by pumping for 4 months at 500 gal/min. The average yearly pumping rate of Scenario 2 is 900 gal/min. The Pumping Test model simulations confirmed that drawdown in the boreholes of the pumping wells at the selected 2-well or 3-well rates will meet the desired condition that the pumping water level remains at least 100 ft above the highest Cambrian-Ordovician unit open to the well. The LMB modified model was used to evaluate the regional drawdown of the pumping from the confined COAS under the same 2-well and 3-well scenarios. At the nearest known existing COAS well, Campbellsport production well #4, the simulated drawdown for Scenario 1 after 20 years of cyclical pumping with two pumping wells averaging a total of 800 gal/min is 16.9 ft, whereas the simulated drawdown for Scenario 2 after 20 years of pumping with three pumping wells averaging a total of 900 gal/min is 19.0 ft. The total deep aquifer thickness at the Campbellsport location is on the order of 620 ft, meaning that the simulated drawdown for either scenario is about 3 percent of the confined aquifer thickness. The models developed as part of this project are archived in the project data release. The archive includes the model input and output files as well as MODFLOW source code and executables. ( Haserodt and others, 2017 ).

Wisconsin↗
USGS2018

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↗
USGS2018

Conceptual model to assess water use associated with the life cycle of unconventional oil and gas development

As the demand for energy increases in the United States, so does the demand for water used to produce many forms of that energy. Technological advances, limited access to conventional oil and gas accumulations, and the rise of oil and gas prices resulted in increased development of unconventional oil and gas (UOG) accumulations. Unconventional oil and gas is developed using a method that combines directional drilling and hydraulic fracturing techniques, allowing for greater oil and gas production from previously unrecoverable reservoirs. Quantification of the water resources required for UOG development and production is difficult because of disparate data sources, variable reporting requirements across boundaries (local, State, and national), and incomplete or proprietary datasets. A topical study was started in 2015 under the U.S. Geological Survey’s Water Availability and Use Science Program, as part of the directive in the Secure Water Act for the U.S. Geological Survey to conduct a National Water Census, to better understand the relation between production of UOG resources for energy and the amount of water needed to produce and sustain this type of energy development in the United States. The Water Availability and Use Science Program goal for this topical study is to develop and apply a statistical model to better estimate the water use associated with UOG development, regardless of the location and target geologic formation. As a first step, a conceptual model has been developed to characterize the life cycle of water use in areas of UOG development. Categories of water use and the way water-use data are collected might change over time; therefore, a generic approach was used in developing the conceptual model to allow for greater flexibility in adapting to future changes or newly available data. UOG development can be summarized into four stages: predrilling construction, drilling, hydraulic fracturing, and ongoing production. The water used in UOG production can be categorized further as direct, indirect, or ancillary water use. Direct water use is defined as the water used for drilling and hydraulic fracturing a well and for maintaining the well during ongoing production. Indirect water use is defined as the water used at or near a well pad. The water used for dust abatement also is considered an indirect use but may be applied away from the well pad. Ancillary water use is defined as the additional local or regional water use resulting from a change (for example, population) directly related to UOG development throughout the life cycle that is not used directly in the well or indirectly for any other purpose at the well pad. The conceptual model presented in this report consists of five elements: (1) input data, (2) processes, (3) decisions, (4) output data, and (5) outcomes. The input data requirements for estimating water use associated with UOG development are somewhat onerous, and obtaining suitable datasets can be challenging because local, State, and Federal agencies do not collect data similarly. The quality of a water-use assessment that uses the conceptual model presented in this report is dependent on the quality and quantity of data that are available for a UOG play. The conceptual model can be used for an assessment with sparse data; however, having sparse data likely will result in greater uncertainty in the water-use estimates. The conceptual model presented in this report is designed to be robust to characterize and simulate the data processing to estimate water use associated with UOG development. Although the results of an analysis that includes missing data have greater uncertainty, the analysis still can be insightful because it can establish a baseline estimate of UOG water use that may be refined further as more data become available. Analysis of models that include missing data also could aid in identifying the data most needed for future water-use estimates. Characterizing individual model limitations is important because the conceptual model can be used in future water-use studies to facilitate data compiling, data processing, estimating, and assessing UOG activities regardless of location.

Scientific Investigations Report↗
USGS2018

Spatial and temporal trends in selenium in the upper Blackfoot River watershed, southeastern Idaho, 2001–16

Phosphate mining in southeastern Idaho has been an important economic driver for the region and State for over 100 years, but weathering of mining waste rock has also released selenium into the Blackfoot River. This report analyzes and presents data from three separate but complementary studies monitoring selenium in streams in the region. The U.S. Geological Survey (USGS), in cooperation with the Bureau of Land Management, has been collecting streamflow and water-quality samples year-round on the Blackfoot River above reservoir near Henry, Idaho, (USGS streamgage 13063000) since 2001. Over the same period, the Idaho Department of Environmental Quality (IDEQ) has collected streamflow and water-quality samples from the Blackfoot River and tributaries during spring runoff. Data collected from 2001 to 2012 during these two studies were analyzed previously. This report extends the analysis using new data collected through 2016. This report also presents the results of a joint USGS and IDEQ seepage study conducted in June 2016 in the Blackfoot River near Dry Valley. Although limited in scope, this study explored the hypothesis that unaccounted selenium loading (loading in excess of tributary inputs) in this reach could be caused by groundwater inflow. USGS dissolved selenium concentration data from streamgage 13063000 on the Blackfoot River and IDEQ data from the mainstem and mining-affected tributaries are highest shortly after peak runoff and correlate with streamflow magnitude. Although earlier analyses indicated increasing selenium concentrations from 2001 to 2012, this study shows that runoff and baseflow dissolved selenium concentrations increased and then decreased during 2001–16. High median runoff concentrations from 2005 through 2011 are associated with high snowpack and streamflow. This result suggests that more snowmelt moving through selenium-bearing waste rock leads to increased instream concentrations. The time lag between peak runoff and then peak selenium concentrations suggests that selenium mobilization may occur as snowmelt percolates through waste rock rather than by faster surface runoff. However, variability in local snow accumulation and snowmelt conditions likely affects interannual variability in selenium concentrations in the mainstem Blackfoot River and tributaries. In contrast to runoff selenium concentrations, median baseflow (August to October) dissolved selenium concentrations were highest from 2009 to 2013. Aquatic plant senescence and release of selenium is an unlikely explanation for this trend because plants are still growing during this time of year. In addition, this trend is observed during and shortly after the observed period of high snowpack. Thus, increased baseflow selenium concentrations suggest that increased selenium loading to alluvial groundwater may occur during periods of high snowmelt and manifest in later years as higher instream concentrations during baseflows when the majority of streamflow is attributable to groundwater gains. Runoff-period streamflow and selenium loads were calculated for the tributaries and mainstem Blackfoot River. Selenium loads vary from year to year with mainstem loads greater than the total tributary contributions in some years and less than tributary contributions in other years. In general, East Mill Creek usually accounted for the largest proportion of the total Blackfoot River load, and unaccounted loads (loads in excess of tributary inputs) often occurred in the vicinity of Spring Creek and Dry Valley. The latter observation led the USGS and IDEQ to conduct a seepage study to further investigate groundwater and selenium loading to the Blackfoot River near Dry Valley. The seepage study results show consistent albeit small unaccounted increases in streamflow and dissolved selenium load in the Blackfoot River near Dry Valley. Field observation of a spring to the north of the river and independent groundwater monitoring data from Dry Valley to the south of the river suggest that alluvial groundwater may discharge to the river from both sides. However, the small unaccounted selenium load measured in the June 2016 study relative to loads measured during runoff suggest that groundwater loading in this reach may occur primarily during runoff. An improved understanding of alluvial groundwater extent, gradient, hydraulic conductivity, and quality would aid in interpreting unaccounted gains and losses in selenium loads in the Blackfoot River. Finally, State of Idaho selenium water-quality criteria have recently shifted to a hierarchical fish tissue and water concentration scheme. This report summarizes existing fish tissue and water-quality data in the mainstem and offers considerations for future selenium monitoring in the Blackfoot River.

Idaho↗
USGS2018

Conceptual and numerical models of dissolved solids in the Colorado River, Hoover Dam to Imperial Dam, and Parker Dam to Imperial Dam, Arizona, California, and Nevada

Conceptual and numerical models were developed to understand and simulate monthly flow-weighted dissolved-solids concentrations in the Colorado River at Imperial Dam. The ability to simulate dissolved-solids concentrations at this location will help the Bureau of Reclamation satisfy the binational agreement on the volume and salinity of Colorado River water delivered to Mexico. A robust spatial- and temporal-resolution dataset that consists of river discharge and dissolved-solids concentration and load information between January 1990 and September 2016 for 10 sites on canals, drains, tributaries, and the main stem of the Colorado River between Hoover and Imperial Dams was generated. Daily mean dissolved-solids concentrations were estimated and monthly mean dissolved-solids loads were computed for each site. Spatial and temporal load patterns, and historical and current controls on loads and concentrations, were analyzed in order to develop a conceptual model of dissolved-solids transport between Hoover and Imperial Dams. Two numerical models describing the relations between dissolved-solids concentrations and components controlling dissolved-solids concentrations and loads were developed, calibrated, and verified. Between January 1990 and September 2016, there was a 98.8-million-acre-feet loss of water and a 57.0-million-ton loss of dissolved-solids load from the Colorado River between Hoover and Imperial Dams. Between Hoover and Parker Dams, about 69.0 million acre-feet of water was lost and 51.1 million tons of dissolved solids were lost; between Parker and Imperial Dams, about 29.8 million acre-feet of water was lost and 5.9 million tons of dissolved solids were lost. Water was removed from the river at a relatively consistent rate over the 25-year study period through water transfers to California and Arizona, evapotranspiration from crop irrigation, transpiration processes of riparian vegetation, and evaporation from the river main stem. Dissolved solids were removed from the river between Hoover and Parker Dams at a relatively constant rate through water transfers to California and Arizona, and water pumped from the river for irrigation within the Mohave Valley. A small amount of dissolved solids are gained by the river from inflow from the Bill Williams River. Between Parker and Imperial Dams, however, dissolved solids were not removed from the river at a consistent rate over the study period. Dissolved solids were generally removed from the river from 1990 to 2012, then gained by the river from 2012 to 2015, and then removed from the river from 2015 through 2016. Dissolved solids are assumed to be removed from the river and accumulated within the floodplain sediments and aquifers during irrigation processes; some dissolved solids may also be removed from the river through uptake by crops and riparian vegetation. Dissolved solids accumulated on the landscape and in the floodplain aquifer during irrigation are transported to the river during periods when the hydraulic gradient between the floodplain aquifer and the river is increased, causing a gain in dissolved solids in the river. Dissolved-solids gains in the river occur during periods of relatively low river discharge, such as during the winter months and during drier climatic conditions. Two numerical models were developed and coefficients were estimated by using data from a May 2008-September 2016 calibration period. One model simulates concentrations at Imperial Dam based on the Colorado River system downstream from Parker Dam, and the other model simulates concentrations at Imperial Dam based on the Colorado River system downstream from Hoover Dam. Both models simulated monthly flow-weighted concentrations of dissolved solids for the Colorado River at Imperial Dam, which corresponded well with observed concentrations for the entire study period. The models are more sensitive to input variables of monthly discharge of the Colorado River below Parker Dam and monthly flow-weighted dissolved-solids concentrations of the Colorado River below Hoover Dam and Parker Dam than to the rate of change in concentration with respect to time and the combined discharge of the Colorado River Indian Reservation Main Canal and the Palo Verde Canal. The calibrated models can be used to run scenarios of future monthly flow-weighted dissolved-solids concentrations in the Colorado River at Imperial Dam. Although the models are expected to provide concentration estimates within 18 milligrams per liter (Parker Dam to Imperial Dam model) to 22 milligrams per liter (Hoover Dam to Imperial Dam model), 95 percent of the time, the error of future scenarios increases as uncertainty in the estimated future input variables increases.

Arizona, California, Nevada↗
USGS2019

Assessing water quality from highway runoff at selected sites in North Carolina with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2015, the U.S. Geological Survey (USGS) entered into a cooperative agreement with the North Carolina Department of Transportation (NCDOT) to develop a North Carolina-enhanced variation of the national Stochastic Empirical Loading and Dilution Model (SELDM) with available North Carolina-specific streamflow and water-quality data and to demonstrate use of the model by documenting selected simulation scenarios. The USGS developed the national SELDM in cooperation with the Federal Highway Administration (FHWA) to provide the tools and techniques necessary for performing stormwater-quality simulations. SELDM uses a stochastic mass-balance approach to estimate combinations of flows, concentrations, and loads of stormwater constituents from the site of interest (often a highway catchment; nonhighway areas, such as a large impervious area at a shopping center complex, also can be used) and the basin upstream from the stormwater outfall to assess the risk for adverse effects of runoff. SELDM also can be used to simulate the effectiveness of volume reduction, hydrograph extension, and water-quality concentration reductions by stormwater best management practices (BMPs), which are designed to help mitigate the effects of runoff on receiving water bodies. Some of the statistical inputs needed for the North Carolina-enhanced SELDM were either calculated or augmented using local or regional data from North Carolina. Streamflow statistics used by SELDM were determined for 266 streamgages across North Carolina on the basis of data available through the 2015 water year. Recession ratio statistics used for triangular hydrographs were also developed for 30 streamgages across the State. The NCDOT identified previous research reports on highway-runoff and BMP studies in North Carolina for review of potential data addition to the national FHWA Highway-Runoff Database (HRDB). Following USGS review of these data, a total of 25,087 event mean concentration values and 1,140 storm events for 39 highway-runoff sites and 195 analytes were uploaded to the national HRDB from six North Carolina highway-runoff research reports and a recent USGS bridge deck runoff study. Using data for 27 streamgages in North Carolina, a total of 57 water-quality transport curves were developed for seven constituents for use in simulating water-quality conditions in the upstream basin. Performance data for three BMPs (bioretention, grass strip or swale, and wetland channel) from NCDOT research data were incorporated into the North Carolina-enhanced SELDM for volume-reduction statistics, including the effectiveness of treating four water-quality constituents (total suspended solids, total nitrogen, total phosphorus, nitrate plus nitrite) and turbidity. Simulations using the North Carolina-enhanced SELDM are presented for two hypothetical upstream basins in the Piedmont ecoregion and one hypothetical highway site to demonstrate how simulations can be used to provide risk-based information about potential effects of stormwater runoff on downstream water quality and the potential for mitigating those risks by using BMPs. The first group of simulations explores the stochastic variability in dilution factors (the ratio of the highway runoff to the total downstream stormflow) for a hypothetical Piedmont rural creek having drainage areas ranging from 1 to 100 square miles. The second group of simulations examines dilution factors based on variations in precipitation, streamflow, and recession ratios for two hypothetical Piedmont upstream basins (rural and urban) where the drainage area was held constant at 25 square miles. These simulations indicate the sensitivity of results to variations in each of the three variables. The third group of simulations examines the effects of varied concentrations in the upstream basin on water-quality conditions downstream from the highway crossing. Variations in upstream water-quality conditions for three constituents (suspended sediment concentration, total nitrogen, and total phosphorus) are based on water-quality transport curves selected from among the 57 curves developed as part of this study to represent low-, medium-, and high-concentration statistics. Simulations completed for this third group also examine the potential effects of grass swale and bioretention BMP treatment on total nitrogen and total phosphorus concentrations in highway runoff. The BMP performance data from the NCDOT research reports were applied in this group of simulations. The stochastic mass-balance approach used in SELDM analyses and simulations provides a strong tool for engineers and water-resource managers to use in exploring a wide range of possible hydrologic and water-quality inputs and their effects on downstream water quality. The results of this study can not only aid engineers and managers in planning for potential adverse effects of runoff at site-specific locations, they can also help the USGS and other Federal and State agencies with oversight responsibilities in stormwater-quality issues to continue gathering data on potential water-quality effects in receiving streams.

North Carolina↗