Geology Reports⌕ Search

SEARCH · Geology Reports

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,045 records · Page 58Linked to original sources

Simulation of streamflows and basin-wide hydrologic variables over several climate-change scenarios, Methow River basin, Washington

The purpose of this project was to demonstrate the capabilities of an existing watershed model and downscaling procedures to provide simulated hydrological data over various greenhouse gas emission scenarios for use in the Methow River framework prototype. An existing watershed model was used to simulate daily time series of streamflow and basin-wide hydrologic variables for baseline conditions (1990–2000), and then for all combinations of three greenhouse gas emission scenarios and five general circulation models for future conditions (2008–2095). Input data for 18 precipitation and 17 temperature model input sites were generated using statistical techniques to downscale general circulation model data. The simulated results were averaged using an 11-year moving window to characterize the central year of the window to provide simulated data for water years 2008–2095. Simulation results indicate that substantial changes of monthly mean streamflows will occur. For all greenhouse gas emission scenarios, the future streamflows are greater in the winter than baseline conditions because a greater percentage of future precipitation is projected to fall as rain rather than as snow. The simulated future spring streamflows are less than baseline conditions because the spring snowpacks are smaller; therefore, flow contributions from snowmelt are less. A database was developed to automate model execution and to provide users with Internet access to voluminous data products ranging from summary figures to model output timeseries. Database-enabled Internet tools were developed to allow users to create interactive graphs of output results based on their analysis needs. For example, users were able to create graphs by selecting time intervals, greenhouse gas emission scenarios, general circulation models, and specific hydrologic variables.

Washington↗

Nutrient and suspended-sediment trends, loads, and yields and development of an indicator of streamwater quality at nontidal sites in the Chesapeake Bay watershed, 1985-2010

The U.S. Geological Survey (USGS) updates information on loads of, and trends in, nutrients and sediment annually to help the Chesapeake Bay Program (CBP) investigators assess progress toward improving water-quality conditions in the Chesapeake Bay and its watershed. CBP scientists and managers have worked since 1983 to improve water quality in the bay. In 2010, the U.S. Environmental Protection Agency (USEPA) established a Total Maximum Daily Load (TMDL) for the Chesapeake Bay. The TMDL specifies nutrient and sediment load allocations that need to be achieved in the watershed to improve dissolved oxygen, water-clarity, and chlorophyll conditions in the bay. The USEPA, USGS, and state and local jurisdictions in the watershed operate a CBP nontidal water-quality monitoring network and associated database that are used to update load and trend information to help assess progress toward reducing nutrient and sediment inputs to the bay. Data collected from the CBP nontidal network were used to estimate loads and trends for two time periods: a long-term period (1985-2010) at 31 "primary" sites (with storm sampling) and a 10-year period (2001-10) at 33 primary sites and 16 "secondary" sites (without storm sampling). In addition, loads at 64 primary sites were estimated for the period 2006 to 2010. Results indicate improving flow-adjusted trends for nitrogen and phosphorus for 1985 to 2010 at most of the sites in the network. For nitrogen, 21 of the 31 sites showed downward (improving) trends, whereas 2 sites showed upward (degrading) trends, and 8 sites showed no trends. The results for phosphorus were similar: 22 sites showed improving trends, 4 sites showed degrading trends, and 5 sites indicated no trends. For sediment, no trend was found at 40 percent of the sites, with 10 sites showing improving trends and 8 sites showing degrading trends. The USGS, working with CBP partners, developed a new water-quality indicator that combines the results of the 10-year trend analysis with results from a greater number of sites (64 primary sites) where loads and yields of total nitrogen and phosphorus and sediment could be calculated. The new indicator shows fewer significant trends for the 10-year time period than for the long-term time period (1985-2010). For 2001-10, total nitrogen trends were downward (improving) at 14 sites and upward (degrading) at 2 sites; no trend was found at 17 sites. For total phosphorus, 12 sites showed improving trends, 4 sites showed degrading trends, and 17 sites showed no trend. For total sediment, most sites (21) did not exhibit a significant trend; 3 sites showed improving trends, and 10 sites showed degrading trends. Few significant trends were seen at the 16 secondary sites: improving trends for total nitrogen at 4 sites, improving trends for total phosphorus at 2 sites, and a degrading trend for sediment at 1 site. Total streamflow to the Chesapeake Bay was 20 percent higher in 2010 than in 2009 and is considered to be within the normal range of flow, whereas annual streamflow at 28 sites was greater in 2010 than in 2009. No trends in daily streamflow were detected at the 31 long-term sites. Combined loads for the farthest downstream nontidal monitoring sites (called "River Input Monitoring sites") increased 33 percent for total nitrogen, 120 percent for total phosphorus, and 330 percent for total sediment from 2009 to 2010. The large increase in phosphorus and sediment loads in 2010 was caused in large part by two large storm events that occurred during the spring in the Potomac River Basin. Yields (load per watershed area) of total nitrogen in the Chesapeake Bay watershed decreased from north to south (New York to Virginia). No spatial patterns were discernible for total phosphorus or sediment.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Estimates of gains and losses from unmeasured sources and sinks for streamflow and dissolved-solids load in selected reaches of the Arkansas River, southeastern Colorado, 2009-2010

The Arkansas River is an important municipal water supply and is the primary supply for about 400,000 acres of irrigated land in southeastern Colorado. The suitability of this water for domestic, agricultural, and industrial use is affected by high salinity in parts of the Arkansas River. There is a need to quantify mass loading of dissolved solids (DS) in the Arkansas River. In 2009, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group and the Colorado Water Conservation Board, began a study to estimate gains and losses from unmeasured sources and sinks for streamflow and DS load in selected reaches of the Arkansas River in southeastern Colorado. Two study reaches were selected for investigation—Canon City to just upstream from Pueblo Reservoir (UARB) and Avondale to Las Animas (LARB). The results from the water-budget analyses indicated that potential areas of unmeasured sources and sinks of streamflow were identifiable in the two study reaches. In the UARB, a substantial volume of water in the subreach from Ark at Canon City to the seasonal gaging station 5 miles downstream (Ark nr Canon City) was unaccounted for by the methodology used in this analysis. The daily gain from unmeasured sources in this subreach was estimated to be about 100 cubic feet per second (ft 3 /s) or about 20 ft 3 /s per river mile. Water-budget estimates for the remaining 18 miles of the UARB study reach indicated that gains or losses from unmeasured sources or sinks were within the measurement error as defined for this report. In the LARB, gains and losses from unmeasured sources and sinks were identified in some of the subreaches but the magnitude of the flux generally was small. Unmeasured sources ranging from less than 2 to 3 ft 3 /s per mile were identified in the river subreaches from Ark at Catlin Dam downstream to Ark at Swink. A streamflow loss was indicated along the subreach from Ark at Nepesta to Ark at Catlin Dam, particularly in 2010. The mechanism and spatial extent of this sink was not identified, and further investigation would be required to better quantify the loss. The results from the analyses of unmeasured sources of DS load indicated that potential source areas were identifiable in the study areas. It might be expected that unmeasured DS load flux would be identified along the same reaches where unmeasured streamflow flux was identified. To that extent, some of the observed results from the analysis of daily DS loading did mirror the streamflow results. In some subreaches of the Arkansas River, however, unmeasured sources and sinks of DS load did not appear to be directly associated with unmeasured sources and sinks of streamflow. In the UARB from Ark at Canon City to Ark nr Canon City, unmeasured gains in DS load were estimated to range from 11 to 22 tons per day per mile in 2009 and from about 8 to 13 tons per day per mile in 2010; streamflow from unmeasured sources was estimated to be about 20 ft 3 /s per mile along this same reach. Downstream from this short reach, DS load to the river from unmeasured sources was estimated to range from 5.4 to 7.6 tons per day per mile in 2010 for Ark nr Canon City to Ark at Portland and from 11 to 16 tons per day per mile in 2009 for Ark at Portland to Ark nr Portland. Unmeasured gains in streamflow were not identified in either of these subreaches. Several small tributaries with DS concentrations ranging from 3,000 mg/L to as high as 6,000 mg/L enter the river along these subreaches. These inputs may indicate a potential source of groundwater that could affect DS loading in the river. Further investigation would be needed to identify the unmeasured source or sources of DS load to determine the nature and extent of unmeasured inputs. In the LARB, gains in DS load from unmeasured sources were identified for the subreach from Ark nr Avondale to Ark at Nepesta, although no substantial amounts of streamflow from unmeasured sources were identified for this subreach. In 2009, the estimated gain in DS load from unmeasured sources for this subreach was 4.7 tons per day per mile. An increase in DS load from unmeasured sources also was identified along the subreach of the river from Ark at Catlin to Swink; the DS load from unmeasured sources was estimated to range from 10 to 28 tons per day per mile. The only loss of DS load was identified for the subreach from Nepesta to Catlin Dam in 2010. The mechanism and spatial extent of the losses were not identified, and further investigation would be required to better understand the results.

Colorado↗

Evaluation of the groundwater-flow model for the Ohio River alluvial aquifer near Carrollton, Kentucky, updated to conditions in September 2010

The Ohio River alluvial aquifer near Carrollton, Ky., is an important water resource for the cities of Carrollton and Ghent, as well as for several industries in the area. The groundwater of the aquifer is the primary source of drinking water in the region and a highly valued natural resource that attracts various water-dependent industries because of its quantity and quality. This report evaluates the performance of a numerical model of the groundwater-flow system in the Ohio River alluvial aquifer near Carrollton, Ky., published by the U.S. Geological Survey in 1999. The original model simulated conditions in November 1995 and was updated to simulate groundwater conditions estimated for September 2010. The files from the calibrated steady-state model of November 1995 conditions were imported into MODFLOW-2005 to update the model to conditions in September 2010. The model input files modified as part of this update were the well and recharge files. The design of the updated model and other input files are the same as the original model. The ability of the updated model to match hydrologic conditions for September 2010 was evaluated by comparing water levels measured in wells to those computed by the model. Water-level measurements were available for 48 wells in September 2010. Overall, the updated model underestimated the water levels at 36 of the 48 measured wells. The average difference between measured water levels and model-computed water levels was 3.4 feet and the maximum difference was 10.9 feet. The root-mean-square error of the simulation was 4.45 for all 48 measured water levels. The updated steady-state model could be improved by introducing more accurate and site-specific estimates of selected field parameters, refined model geometry, and additional numerical methods. Collection of field data to better estimate hydraulic parameters, together with continued review of available data and information from area well operators, could provide the model with revised estimates of conductance values for the riverbed and valley wall, hydraulic conductivities for the model layer, and target water levels for future simulations. Additional model layers, a redesigned model grid, and revised boundary conditions could provide a better framework for more accurate simulations. Additional numerical methods would identify possible parameter estimates and determine parameter sensitivities.

Kentucky↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗