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.

SourceResearch connectionFrom yearThrough year

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

At least 559 records · Page 31Linked to original sources

USGS2010

Sedimentation, sediment quality, and upstream channel stability, John Redmond Reservoir, east-central Kansas, 1964-2009

A combination of available bathymetric-survey information, bottom-sediment coring, and historical streamgage information was used to investigate sedimentation, sediment quality, and upstream channel stability for John Redmond Reservoir, east-central Kansas. Ongoing sedimentation is reducing the ability of the reservoir to serve several purposes including flood control, water supply, and recreation. The total estimated volume and mass of bottom sediment deposited between 1964 and 2009 in the conservation pool of the reservoir was 1.46 billion cubic feet and 55.8 billion pounds, respectively. The estimated sediment volume occupied about 41 percent of the conservation-pool, water-storage capacity of the reservoir. Water-storage capacity in the conservation pool has been lost to sedimentation at a rate of about 1 percent annually. Mean annual net sediment deposition since 1964 in the conservation pool of the reservoir was estimated to be 1.24 billion pounds per year. Mean annual net sediment yield from the reservoir basin was estimated to be 411,000 pounds per square mile per year Information from sediment cores shows that throughout the history of John Redmond Reservoir, total nitrogen concentrations in the deposited sediment generally were uniform indicating consistent nitrogen inputs to the reservoir. Total phosphorus concentrations in the deposited sediment were more variable than total nitrogen indicating the possibility of changing phosphorus inputs to the reservoir. As the principal limiting factor for primary production in most freshwater environments, phosphorus is of particular importance because increased inputs can contribute to accelerated reservoir eutrophication and the production of algal toxins and taste-and-odor compounds. The mean annual net loads of total nitrogen and total phosphorus deposited in the bottom sediment of the reservoir were estimated to be 2,350,000 pounds per year and 1,030,000 pounds per year, respectively. The estimated mean annual net yields of total nitrogen and total phosphorus from the reservoir basin were 779 pounds per square mile per year and 342 pounds per square mile per year, respectively. Trace element concentrations in the bottom sediment of John Redmond Reservoir generally were uniform over time. As is typical for eastern Kansas reservoirs, arsenic, chromium, and nickel concentrations typically exceeded the threshold-effects guidelines, which represent the concentrations above which toxic biological effects occasionally occur. Trace element concentrations did not exceed the probable-effects guidelines (available for eight trace elements), which represent the concentrations above which toxic biological effects usually or frequently occur. Organochlorine compounds either were not detected or were detected at concentrations that were less than the threshold-effects guidelines. Stream channel banks, compared to channel beds, likely are a more important source of sediment to John Redmond Reservoir from the upstream basin. Other sediment sources include surface-soil erosion in the basin and shoreline erosion in the reservoir.

Scientific Investigations Report↗

USGS2011

Simulation of specific conductance and chloride concentration in Abercorn Creek, Georgia, 2000-2009

The City of Savannah operates an industrial and domestic water-supply intake on Abercorn Creek approximately 2 miles from the confluence with the Savannah River upstream from the Interstate 95 bridge. Chloride concentrations are a major concern for the city because industrial customers require water with low chloride concentrations, and elevated chloride concentrations require additional water treatment in order to meet those needs. The proposed deepening of Savannah Harbor could increase chloride concentrations (the major ion in seawater) in the upper reaches of the lower Savannah River estuary, including Abercorn Creek. To address this concern, mechanistic and empirical modeling approaches were used to simulate chloride concentrations at the city's intake to evaluate potential effects from deepening the Savannah Harbor. The first approach modified the mechanistic Environmental Fluid Dynamics Code (EFDC) model developed by Tetra Tech and used for evaluating proposed harbor deepening effects for the Environmental Impact Statement. Chloride concentrations were modeled directly with the EFDC model as a conservative tracer. This effort was done by Tetra Tech under a separate funding agreement with the U.S. Army Corps of Engineers and documented in a separate report. The second approach, described in this report, was to simulate chloride concentrations by developing empirical models from the available data using artificial neural network (ANN) and linear regression models. The empirical models used daily streamflow, specific conductance (field measurement for salinity), water temperature, and water color time series for inputs. Because there are only a few data points that describe the relation between high specific conductance values at the Savannah River at Interstate 95 and the water plant intake, there was a concern that these few data points would determine the extrapolation of the empirical model and potentially underestimate the effect of deepening the harbor on chloride concentrations at the intake. To accommodate these concerns, two ANN chloride models were developed for the intake. The first model (ANN M1e) used all the data. The second model (ANN M2e) only used data when specific conductance at Interstate 95 was less than 175 microsiemens per centimeter at 25 degrees Celsius. Deleting the conductivity data greater than 175 microsiemens per centimeter removed the "plateau" effect observed in the data. The chloride simulations with the ANN M1 model have a low sensitivity to specific conductance (salinity) at Interstate 95, whereas the chloride simulations with the ANN M2 model have a high sensitivity to salinity at Interstate 95. The two modeling approaches (Tetra Tech's EFDC model and the one described in this report) were integrated into a decision support system (DSS) that combines the historical database, output from EFDC, ANN models, ANN model simulation controls, streaming graphics, and model output. The DSS was developed as a Microsoft ExcelTM/Visual Basic for Applications program, which allowed the DSS to be prototyped, easily modified, and distributed in a familiar spreadsheet format. The EFDC and ANN models were used to simulate various harbor deepening scenarios. To accommodate the geometry changes in the harbor, the ANN models used the EFDC model-simulated salinity changes for a historical condition as input. The DSS uses a graphical user interface and allows the user to interrogate the ANN models and EFDC output. Two scenarios were simulated using the Savannah Chloride Model DSS to demonstrate different input options. One scenario decreased winter streamflows to a constant streamflow for 45 days. Streamflows during the period January 1 to February 15 were set to a constant 3,600 cubic feet per second for the simulation period of October 1, 2006, to October 1, 2009. The decreased winter streamflow resulted in predictions of increased specific conductance by as much as 50 microsiemens per centimeter and chlorid

Georgia↗

USGS2011

Simulation of streamflow, evapotranspiration, and groundwater recharge in the Lower Frio River watershed, south Texas, 1961-2008

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, Fort Worth District; the City of Corpus Christi; the Guadalupe-Blanco River Authority; the San Antonio River Authority; and the San Antonio Water System, configured, calibrated, and tested a watershed model for a study area consisting of about 5,490 mi 2 of the Frio River watershed in south Texas. The purpose of the model is to contribute to the understanding of watershed processes and hydrologic conditions in the lower Frio River watershed. The model simulates streamflow, evapotranspiration (ET), and groundwater recharge by using a numerical representation of physical characteristics of the landscape, and meteorological and streamflow data. Additional time-series inputs to the model include wastewater-treatment-plant discharges, surface-water withdrawals, and estimated groundwater inflow from Leona Springs. Model simulations of streamflow, ET, and groundwater recharge were done for various periods of record depending upon available measured data for input and comparison, starting as early as 1961. Because of the large size of the study area, the lower Frio River watershed was divided into 12 subwatersheds; separate Hydrological Simulation Program-FORTRAN models were developed for each subwatershed. Simulation of the overall study area involved running simulations in downstream order. Output from the model was summarized by subwatershed, point locations, reservoir reaches, and the Carrizo-Wilcox aquifer outcrop. Four long-term U.S. Geological Survey streamflow-gaging stations and two short-term streamflow-gaging stations were used for streamflow model calibration and testing with data from 1991-2008. Calibration was based on data from 2000-08, and testing was based on data from 1991-99. Choke Canyon Reservoir stage data from 1992-2008 and monthly evaporation estimates from 1999-2008 also were used for model calibration. Additionally, 2006-08 ET data from a U.S. Geological Survey meteorological station in Medina County were used for calibration. Streamflow and ET calibration were considered good or very good. For the 2000-08 calibration period, total simulated flow volume and the flow volume of the highest 10 percent of simulated daily flows were calibrated to within about 10 percent of measured volumes at six U.S. Geological Survey streamflow-gaging stations. The flow volume of the lowest 50 percent of daily flows was not simulated as accurately but represented a small percent of the total flow volume. The model-fit efficiency for the weekly mean streamflow during the calibration periods ranged from 0.60 to 0.91, and the root mean square error ranged from 16 to 271 percent of the mean flow rate. The simulated total flow volumes during the testing periods at the long-term gaging stations exceeded the measured total flow volumes by approximately 22 to 50 percent at three stations and were within 7 percent of the measured total flow volumes at one station. For the longer 1961-2008 simulation period at the long-term stations, simulated total flow volumes were within about 3 to 18 percent of measured total flow volumes. The calibrations made by using Choke Canyon reservoir volume for 1992-2008, reservoir evaporation for 1999-2008, and ET in Medina County for 2006-08, are considered very good. Model limitations include possible errors related to model conceptualization and parameter variability, lack of data to better quantify certain model inputs, and measurement errors. Uncertainty regarding the degree to which available rainfall data represent actual rainfall is potentially the most serious source of measurement error. A sensitivity analysis was performed for the Upper San Miguel subwatershed model to show the effect of changes to model parameters on the estimated mean recharge, ET, and surface runoff from that part of the Carrizo-Wilcox aquifer outcrop. Simulated recharge was most sensitive to the changes in the lower-zone ET (LZ

Texas↗

USGS2014

Linkage of the Soil and Water Assessment Tool and the Texas Water Availability Model to simulate the effects of brush management on monthly storage of Canyon Lake, south-central Texas, 1995-2010

The U.S. Geological Survey (USGS), in cooperation with the Texas State Soil and Water Conservation Board, developed and applied an approach to create a linkage between the published upper Guadalupe River Soil Water Assessment Tool (SWAT) brush-management (ashe juniper [ Juniperus ashei ]) model and the full authorization version Guadalupe River Water Availability Model (WAM). The SWAT model was published by the USGS, and the Guadalupe River WAM is available from the Texas Commission on Environmental Quality. The upper Guadalupe River watershed is a substantial component of the Guadalupe River WAM. This report serves in part as documentation of a proof of concept on the feasibility of linking these two water-resources planning models for the purpose of simulating possible increases in water storage in Canyon Lake as a result of different brush-management scenarios. The SWAT-WAM linkage for the upper Guadalupe River is documented with a principal objective to evaluate the distributional characteristics of the monthly water storage of Canyon Lake during selected drought conditions. Focus is on the relative evaluation of select scenarios of large-scale or “extensive” brush management within the upper Guadalupe River watershed. There are six SWAT simulations for the upper Guadalupe River watershed that include a baseline (0-percent management of treatable ashe juniper, the baseline scenario from a previous study in which no percentage of ashe juniper is numerically replaced with grassland) along with five scenarios (extensions of SWAT simulations from a previous study) of 20-, 40-, 60-, 80-, and 100-percent random (numerical) replacement of treatable ashe juniper with grasslands throughout the upper Guadalupe River watershed in south-central Texas. SWAT is a process-based, semidistributed, water-balance model designed to predict the effects of landscape management decisions on water yields. A watershed is subdivided into subbasins, and each subbasin is associated with a single reach on the stream network. In general a WAM, such as the Guadalupe River WAM, provides analysis of generalized water rights in a river and reservoir framework. A WAM accommodates hydrology and water usage through several input files containing water rights, watershed parameters, and naturalized streamflow time series. A WAM is generalized for application to rivers and reservoir systems, and input datasets are uniquely developed for a river basin of concern. The extractions of SWAT output for the five extensive brush-management and baseline scenarios were offset by –21 years and, in general, the results were then mapped to the WAM input-flow file. The offset of –21 years was chosen arbitrarily for technical reasons and means that the period of monthly record 1995–2010 of the upper Guadalupe River SWAT became the synthetic period of monthly record 1974–89, hereinafter 1974–89 (synthetic) period, of the Guadalupe River WAM. The relative (between scenario to baseline) effects of extensive brush-management scenarios by using the SWAT-WAM linkage were evaluated, and two critical intermediate results were total inflow to Canyon Lake from 1995 to 2010 and the monthly storage of Canyon Lake from 1974 to 1989 (synthetic). The first quartile or lower 25th percentile of monthly storage of Canyon Lake for the baseline scenario is 381,000 acre-feet (acre-ft) for the hereinafter 1974–89 (synthetic) period. This lower quartile was chosen for analysis for two critical purposes. First, Canyon Lake is managed with a conservation pool of about 386,200 acre-ft capacity (as recognized by the WAM) and is at or near conservation capacity about 50 percent or more of the time; further, there is intrinsic data censoring that occurs for the monthly storage distribution because Canyon Lake is at or near conservation pool elevation the majority of the time. This intrinsic censoring has the effect of creating a bounded distribution with a left or low-volume tail. Statistical assessment of the brush-management scenarios beginning with the 381,000 acre-ft censoring threshold provides readily interpretable results. Second, the quantification of brush management during periods lacking abundant rainfall, which were defined in this study as months for which Canyon Lake storage was below the 25th percentile for the simulation period, are of substantial interest to water-resource managers and stakeholders in the context of water-supply enhancement. A statistical assessment of the SWAT-WAM linkage for the low-volume tail of the distribution of monthly storage of Canyon Lake is the focus of analysis and interpretation. Drought periods for the analysis are defined as the months (consecutive or not) during which Canyon Lake is below the 25th percentile of storage (381,000 acre-ft) for the baseline scenario. Such months are referred to as being within the “Drought Quartile.” The Drought Quartile is a conceptual and heuristically determined waypoint for the analysis and is not related to any administrative definition of drought by stakeholders or policy makers. The five scenarios and the baseline scenario simulated in the upper Guadalupe River SWAT were all passed through the Guadalupe River WAM by the SWAT-WAM linkage described in this report. A comparison of the mean increase per month in reservoir storage for Canyon Lake conditioned for the Drought Quartile was made. For each of the five brush-management and baseline scenarios, the months with storage below 381,000 acre-ft were extracted. The mean monthly storages during the Drought Quartile were computed for each of the five scenarios and the baseline scenario. The mean of the baseline scenario was 376,458 acre-ft and subsequently was subtracted from the mean monthly storage during the Drought Quartile for each of the five scenarios. The mean monthly offset storages of Canyon Lake during the Drought Quartile were 110 acre-ft (20 percent); 448 acre-ft (40 percent); 754 acre-ft (60 percent); 1,080 acre-ft (80 percent); and 1,090 acre-ft (100 percent). A particular mean was interpreted as follows: the value of 754 acre-ft for the 60-percent brush-management scenario implies that, on average, this scenario indicates an additional 754 acre-ft per month of storage in Canyon Lake relative to the baseline during the Drought Quartile. All of the five scenarios resulted in an increase on average to water supply relative to the baseline scenario during the Drought Quartile through the SWAT-WAM linkage.

Texas↗

USGS2014

Assessing potential effects of highway runoff on receiving-water quality at selected sites in Oregon with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2012, the U.S. Geological Survey and the Oregon Department of Transportation began a cooperative study to demonstrate use of the Stochastic Empirical Loading and Dilution Model (SELDM) for runoff-quality analyses in Oregon. SELDM can be used to estimate stormflows, constituent concentrations, and loads from the area upstream of a stormflow discharge site, from the site of interest and in the receiving waters downstream of the discharge. SELDM also can be used to assess the potential effectiveness of best management practices (BMP) for mitigating potential effects of runoff in receiving waters. Nominally, SELDM is a highway-runoff model, but it is well suited for analysis of runoff from other land uses as well. This report provides case studies and examples to demonstrate stochastic-runoff modeling concepts and to demonstrate application of the model. Basin characteristics from six Oregon highway study sites were used to demonstrate various applications of the model. The highway catchment and upstream basin drainage areas of these study sites ranged from 3.85 to 11.83 acres and from 0.16 to 6.56 square miles, respectively. The upstream basins of two sites are urbanized, and the remaining four sites are less than 5 percent impervious. SELDM facilitates analysis by providing precipitation, pre-storm streamflow, and other variables by region or from hydrologically similar sites. In Oregon, there can be large variations in precipitation and streamflow among nearby sites. Therefore, spatially interpolated geographic information system data layers containing storm-event precipitation and pre-storm streamflow statistics specific to Oregon were created for the study using Kriging techniques. Concentrations and loads of cadmium, chloride, chromium, copper, iron, lead, nickel, phosphorus, and zinc were simulated at the six Oregon highway study sites by using statistics from sites in other areas of the country. Water‑quality datasets measured at hydrologically similar basins in the vicinity of the study sites in Oregon were selected and compiled to estimate stormflow-quality statistics for the upstream basins. The quality of highway runoff and some upstream stormflow constituents were simulated by using statistical moments (average, standard deviation, and skew) of the logarithms of data. Some upstream stormflow constituents were simulated by using transport curves, which are relations between stormflow and constituent concentrations. Stochastic analyses were done by using SELDM to demonstrate use of the model and to illustrate the types of information that stochastic analyses may provide: 1. An analysis was done to demonstrate use of dilution factors as an initial reconnaissance tool for comparing relative risk among sites. 2. An analysis of hardness-dependent, water-quality criteria was done to illustrate the effects of variations in hardness and flow on the application and interpretation of such criteria. This analysis shows that hardness-dependent criteria can vary by an order of magnitude among storm events because hardness is diluted by stormflows. 3. An analysis of uncertainties in input and output values was done to demonstrate that properly selected robust datasets are needed to represent conditions at a site of interest. This analysis shows that the rate of water-quality exceedances that are measured or simulated may depend on sample size and the luck of the draw. 4. An analysis was done to demonstrate that SELDM and other Monte Carlo models may generate extreme values from input statistics, which may or may not be feasible based on physicochemical or hydrological limits. 5. An analysis of BMP modeling methods was done to demonstrate use of the model for estimating treatment requirements for meeting water-quality objectives. 6. An analysis of the use of grab sampling and nonstochastic upstream modeling methods was done to evaluate the potential effects on modeling outcomes. Additional analyses using surrogate water-quality datasets for the upstream basin and highway catchment were provided for six Oregon study sites to illustrate the risk-based information that SELDM will produce. These analyses show that the potential effects of highway runoff on receiving-water quality downstream of the outfall depends on the ratio of drainage areas (dilution), the quality of the receiving water upstream of the highway, and the concentration of the criteria of the constituent of interest. These analyses also show that the probability of exceeding a water-quality criterion may depend on the input statistics used, thus careful selection of representative values is important.

Oregon↗

USGS2014

Occurrence and transport of nitrogen in the Big Sunflower River, northwestern Mississippi, October 2009-June 2011

The Big Sunflower River Basin, located within the Yazoo River Basin, is subject to large annual inputs of nitrogen from agriculture, atmospheric deposition, and point sources. Understanding how nutrients are transported in, and downstream from, the Big Sunflower River is key to quantifying their eutrophying effects on the Gulf. Recent results from two Spatially Referenced Regressions on Watershed attributes (SPARROW models), which include the Big Sunflower River, indicate minimal losses of nitrogen in stream reaches typical of the main channels of major river systems. If SPARROW assumptions of relatively conservative transport of nitrogen are correct and surface-water losses through the bed of the Big Sunflower River are negligible, then options for managing nutrient loads to the Gulf of Mexico may be limited. Simply put, if every pound of nitrogen entering the Delta is eventually delivered to the Gulf, then the only effective nutrient management option in the Delta is to reduce inputs. If, on the other hand, it can be shown that processes within river channels of the Mississippi Delta act to reduce the mass of nitrogen in transport, other hydrologic approaches may be designed to further limit nitrogen transport. Direct validation of existing SPARROW models for the Delta is a first step in assessing the assumptions underlying those models. In order to characterize spatial and temporal variability of nitrogen in the Big Sunflower River Basin, water samples were collected at four U.S. Geological Survey gaging stations located on the Big Sunflower River between October 1, 2009, and June 30, 2011. Nitrogen concentrations were generally highest at each site during the spring of the 2010 water year and the fall and winter of the 2011 water year. Additionally, the dominant form of nitrogen varied between sites. For example, in samples collected from the most upstream site (Clarksdale), the concentration of organic nitrogen was generally higher than the concentrations of ammonia and nitrate plus nitrite; conversely, at sites farther downstream (that is, at Sunflower and Anguilla), nitrate plus nitrite concentrations were generally higher than concentrations of organic nitrogen and ammonia. In addition to the routinely collected samples, water samples from the Big Sunflower River Basin were collected using a Lagrangian sampling scheme, which attempts to follow a single mass of water through time in order to determine how it changes through processing or other pathways as the water moves downstream. Lagrangian sampling was conducted five times during the study period: (1) April 8–21, 2010, (2) May 12–June 3, 2010, (3) June 15–July 1, 2010, (4) August 23–30, 2010, and (5) May 16–20, 2011. Streamflow conditions were variable for each sampling event because of input from local precipitation and irrigation return flow, and streamflow losses through the streambed. Streamflow and total nitrogen flux increased with drainage area, and the dominant form of nitrogen varied with drainage area size and temporally across sampling events. Results from each method indicate relatively conservative transport of nitrogen within the 160 miles between Clarksdale and Anguilla, providing further validation of the SPARROW models. Furthermore, these results suggest relatively conservative transport of nitrogen from the Big Sunflower River to the Gulf of Mexico and, therefore, imply a fairly close association of nutrient application and export from the Big Sunflower River Basin to the Mississippi River. However, within the Big Sunflower River Basin, two potential nitrogen sinks were identified and include the transport and potential transformation of nitrogen through the streambed and the sequestration and potential transformation of nitrogen above the drainage control structures downstream of Anguilla. By coupling these potential loss mechanisms with nitrogen transport dynamics, it may be possible to further reduce the amount of nitrogen leaving the Big Sunflower River Basin and ultimately arriving at the Gulf of Mexico.

Arkansas;Louisiana;Mississippi↗

USGS2014

Hydrogeology, water resources, and water budget of the upper Rio Hondo Basin, Lincoln County, New Mexico, 2010

The upper Rio Hondo Basin occupies a drainage area of 585 square miles in south-central New Mexico and comprises three general hydrogeologic terranes: the higher elevation “Mountain Block,” the “Central Basin” piedmont area, and the lower elevation “Hondo Slope.” As many as 12 hydrostratigraphic units serve as aquifers locally and form a continuous aquifer on the regional scale. Streams and aquifers in the basin are closely interconnected, with numerous gaining and losing stream reaches across the study area. In general, the aquifers are characterized by low storage capacity and respond to short-term and long-term variations in recharge with marked water-level fluctuations on short (days to months) and long (decadal) time scales. Droughts and local groundwater withdrawals have caused marked water-table declines in some areas, whereas periodically heavy monsoons and snowmelt events have rapidly recharged aquifers in some areas. A regional-scale conceptual water budget was developed for the study area in order to gain a basic understanding of the magnitude of the various components of input, output, and change in storage. The primary input is watershed yield from the Mountain Block terrane, supplying about 38,200 to 42,300 acre-feet per year (acre-ft/yr) to the basin, as estimated by comparing the residual of precipitation and evapotranspiration with local streamgage data. Streamflow from the basin averaged about 21,200 acre-ft/yr, and groundwater output left the basin at an estimated 2,300 to 5,700 acre-ft/yr. The other major output (about 13,500 acre-ft/yr) was by public water supply, private water supply, livestock, commercial and industrial uses, and the Bonito Pipeline. The residual in the water budget, the difference between the totals of the input and output terms or the potential change in storage, ranged from -2,200 acre-ft/yr to +5,300 acre-ft/yr. There is a high degree of variability in precipitation and consequently in the water supply; small variations in annual precipitation can result in major changes in overall watershed yield. Changing water-use patterns, concentrated areas of groundwater withdrawal, and variations in precipitation have created localized areas where water-table declines and diminished surface flow are of concern.

Lincoln, New Mexico↗

USGS2015

Chloride concentrations, loads, and yields in four watersheds along Interstate 95, southeastern Connecticut, 2008-11: factors that affect peak chloride concentrations during winter storms

Chloride (Cl - ) concentrations and loads and other water chemistry characteristics were assessed to evaluate potential effects of road-deicer applications on streamwater quality in four watersheds along Interstate 95 (I–95) in southeastern Connecticut from November 1, 2008, through September 30, 2011. Streamflow and water quality were studied in the Four Mile River, Oil Mill Brook, Stony Brook, and Jordan Brook watersheds, where developed land ranged from 9 to 32 percent. Water-quality samples were collected and specific conductance was measured continuously at paired water-quality monitoring sites, upstream and downstream from I–95. Specific conductance values were related to Cl - concentrations to assist in determining the effects of road-deicing operations on the levels of Cl - in the streams. Streamflow and water-quality data were compared with weather data and with the timing, amount, and composition of deicers applied to State highways. Grab samples were collected during winter stormwater-runoff events, such as winter storms or periods of rain or warm temperatures in which melting takes place. Grab samples were also collected periodically during the spring and summer and during base-flow conditions. The estimated Cl - concentrations at the eight water-quality monitoring sites during winter storms peaked as high as 270 milligrams per liter (mg/L) and were well below the U.S. Environmental Protection Agency (EPA) recommended acute chloride toxicity criterion of 860 mg/L and the chronic 4-day average toxicity criterion of 230 mg/L. Cl - concentrations in streams, particularly at sites downstream from I–95, peaked during increased streamflow in the winter and early spring as a result of deicers applied to roads and washed off by stormwater or meltwater. Cl - concentrations during most of the nonwinter seasons decreased during increases in streamflow because storm runoff was more dilute than base flow. However, peaks in specific conductance and estimated chloride concentrations at streams in more urbanized watersheds corresponded to peaks in streamflow well after winter snow or ice events; these delayed peaks in Cl - concentration likely resulted from deicer residue that remained in melting snow piles and on roadsides and (or) that were flushed from soils and shallow groundwater, then discharged downstream. Estimated peak Cl - concentrations varied with the type of winter storm event and were highest during or after winter storms of frozen precipitation and rain, in which the rain or meltwater effectively washed off the deicers. Estimated peak Cl - concentrations correlated positively with the duration of deicer application but generally not with streamflow. Estimated peak Cl - concentrations during the winter season were highest during low streamflow at most sites. Chloride concentrations varied considerably in shallow groundwater as a result of land-use differences. Cl - concentrations were very high (as high as 800 mg/L) in shallow groundwater downstream from I–95 at the Four Mile River site. Chloride/bromide mass concentration ratios and the proximity of a former landfill and sewage lagoon upstream indicate a likely source of Cl - is landfill leachate and possibly sewage leachate. Cl - loads in streams generally were highest in the winter and early spring. The estimated daily Cl - yield for the four monitoring sites downstream from I–95 ranged from 0.0004 ton per day per square mile for one of the least developed watersheds to 0.052 ton per day per square mile for the watershed with the highest percentage of urban development and impervious area. The estimated median contribution of Cl - load from atmospheric deposition was small and ranged from 0.07 percent of Cl - load at the Jordan Brook watershed to 0.57 percent at the Oil Mill Brook watershed. The Cl - loads in streams (outputs) were compared with Cl - load inputs, which include atmospheric deposition and deicer applications; Cl - load inputs were slightly larger than the Cl - load outputs at most of the sites during most years but do not account for the Cl - load in groundwater leaving the watersheds. A multiple linear regression model was developed to describe the variability of the natural log of peak specific conductance, as well as estimated Cl - concentrations. Five significant variables best explained the variability in the natural log of the peak specific conductance after deicing events: (1) snow on ground before deicing event; (2) winter precipitation with rain; (3) specific conductance in base flow; (4) State-operated road lane miles divided by watershed area; and (5) amount of Cl - from deicers applied to State-operated roads per lane mile. In this report, winter precipitation is defined as any type of precipitation, including frozen precipitation and rain, that occurs during the active deicing season, typically November through March. Frozen precipitation is defined here as snow, sleet, freezing rain, or any winter precipitation except rain. The addition of a lane mile in both directions on I–95 would result in an estimate of approximately 2 to 11 percent increase in Cl - input from deicers applied to I–95 and other roads maintained by Connecticut Department of Transportation. The largest estimated increase in Cl - load was in the watersheds with the greatest number miles of I–95 corridor relative to the total lane miles maintained by Connecticut Department of Transportation. On the basis of these estimates and the estimated peak Cl - concentrations during the study period, it is unlikely that the increased use of deicers on the additional lanes would lead to Cl - concentrations that exceed the aquatic habitat criteria.

Connecticut↗

USGS2015

Groundwater/surface-water interactions in the Bad River Watershed, Wisconsin

A groundwater-flow model was developed for the Bad River Watershed and surrounding area by using the U.S. Geological Survey (USGS) finite-difference code MODFLOW-NWT. The model simulates steady-state groundwater-flow and base flow in streams by using the streamflow routing (SFR) package. The objectives of this study were to: (1) develop an improved understanding of the groundwater-flow system in the Bad River Watershed at the regional scale, including the sources of water to the Bad River Band of Lake Superior Chippewa Reservation (Reservation) and groundwater/surface-water interactions; (2) provide a quantitative platform for evaluating future impacts to the watershed, which can be used as a starting point for more detailed investigations at the local scale; and (3) identify areas where more data are needed. This report describes the construction and calibration of the groundwater-flow model that was subsequently used for analyzing potential locations for the collection of additional field data, including new observations of water-table elevation for refining the conceptualization and corresponding numerical model of the hydrogeologic system. The study area can be conceptually divided into three primary hydrogeologic environments. The first encompasses the southern uplands with relatively low topographic relief, where groundwater-flow is unconfined and occurs primarily in sandy till and glacial outwash overlying Archean-aged crystalline bedrock. The second includes a transitional area of higher topographic relief and shallow depth to bedrock, in the vicinity of ridges formed by steeply dipping, early-Proterozoic aged metasedimentary units of the Marquette Range Supergroup (including the Ironwood Formation), and late-Proterozoic igneous units associated with the Midcontinent Rift System (MRS). Groundwater-flow in this area likely occurs primarily through connected networks of bedrock fractures that are not well characterized, and also in isolated pockets of Quaternary deposits. The third and last hydrogeologic environment includes lowlands along Lake Superior where a deep sandstone aquifer is confined by thick deposits of clay-rich till. Model input was compiled by using both published and unpublished data. Constant flux boundary conditions for the model perimeter were developed from a regional analytic element model described in appendix 1 of this report. Pumping from 26 high-capacity wells within the model area was included. The SFR stream network was developed from the National Hydrography Dataset (NHDPlus Version 2) and hydrography from the Wisconsin Department of Natural Resources (WDNR). Hydraulic conductivity values were determined for each model cell by interpolation from a network of pilot points, within zones representing major hydrogeologic units. Recharge to the groundwater system was estimated on a cell-by-cell basis by using the Soil Water Balance code (SWB), with gridded daily temperature and precipitation data for the period 1980–2011, and GIS coverages of soil and land-surface conditions. Estimated recharge varies considerably, following spatial patterns in the precipitation and soil hydrologic group inputs. The lowest recharge values occur in the Superior lowlands, whereas the highest values occur in the upland areas, especially those underlain by sandy soils, and in the vicinity of bedrock hills. The model was calibrated to groundwater-levels and base flows obtained from the USGS National Water Information System (NWIS) database, and groundwater-levels obtained from the WDNR and Band River Band well-construction databases. Calibration was performed via nonlinear regression by using the parameter-estimation software suite PEST. Groundwater levels and base-flow observations in the calibration dataset were well simulated by the calibrated model, with reasonable values of hydraulic conductivity. The pilot-point parameters that were most constrained by observations during model calibration coincided with the locations containing the most wells (head observations)—especially the population centers of Ashland, Mellen, and other communities along the major highway corridors. Results from the calibrated model illustrate differences in the nature of groundwater-ﬂow within the watershed. In the southern part of the watershed, where bedrock is shallow, groundwater ﬂow paths are relatively short, extending from local recharge areas to adjacent ﬁrst and second-order streams. In contrast, laterally continuous deposits of clay-rich till covering the Superior Lowlands isolate most smaller streams from the sandstone aquifer, allowing for longer ﬂow paths toward larger streams such as the Bad, Marengo, and White Rivers. Approximately three-quarters of all ﬁrst-order stream cells were dry in the Superior Lowlands, compared to only half of ﬁrst-order stream cells in the southern bedrock uplands. The model was used to delineate the groundwatershed for the Bad and Kakagon Rivers. “Groundwatershed” is deﬁned as the area contributing groundwater discharge to one of these streams and their tributaries. The groundwatershed was found to align closely with the surface-watershed, with the most notable exception occurring along the southwestern half of Birch Hill, where surface water drains southwest towards the Potato River, and groundwater ﬂows north and east towards Lake Superior. Similarly, the contributing area of groundwater-ﬂow to the Reservation was delineated. Results indicate the off-Reservation groundwater contributing area to be limited in comparison to the extent of the watershed, extending southward into the highlands underlain by MRS igneous rock units, but not further into the area underlain by the Marquette Range Supergroup. Stable isotope samples were collected from 54 wells within the watershed, to investigate sources of groundwater. Oxygen-18 (δ 18 O) values lower than -13.0 per mil were documented in the sampling, and likely indicate the presence of recharge water from the last glacial period (>9,500 years old) beneath the northern portion of the Reservation, in the vicinity of Odanah, Wisconsin. Finally, a new data-worth analysis of potential new monitoring-well locations was performed by using the model. The relative worth of new measurements was evaluated based on their ability to increase conﬁdence in model predictions of groundwater levels and base ﬂows at 35 locations, under the condition of a proposed open-pit iron mine. Results of the new data-worth analysis, and other inputs and outputs from the Bad River model, are available through an online dynamic web mapping service at ( http://wim.usgs.gov/badriver/ ).

Wisconsin↗

USGS2018

Water pressure and ground vibrations induced by water guns near Brandon Road Lock and Dam and Lemont, Illinois

Multiple geophysical sensors were used to characterize the underwater pressure field and ground vibrations of a seismic water gun and its suitability to deter the movement of Asian carps (particularly the silver [ Hypophthalmichthys molitrix ] and bighead [ Hypophthalmichthys nobilis ] carps) while ensuring the integrity of surrounding structures. The sensors used to collect this information were blast-rated hydrophones, surface- and borehole-mounted geophones, and fixed accelerometers. Results from two separate studies are discussed in this report. The Brandon Road study took place in May 2014, in the Des Plaines River, in a concrete-walled channel downstream of the Brandon Road Lock and Dam near Joliet, Illinois. The Lemont study took place in June 2014, in a segment of the dolomite setblock-lined Chicago Sanitary and Ship Canal near Lemont, Illinois. Two criteria were evaluated to assess the potential deterrence to carp migration, and to minimize the expected effect on nearby structures from discharge of the seismic water gun. The first criterion was a 5-pound-per-square-inch (lb/in 2 ) limit for dynamic underwater pressure variations. The second criterion was a maximum velocity and acceleration disturbance of 0.75 inch per second (in/s) for sensitive machinery (such as the lock gates and pumps) and 2.0 in/s adjacent to canal walls, respectively. The criteria were based on previous studies of fish responses to dynamic pressure variations, and effects of vibrations on the structural integrity of concrete walls. The Brandon Road study evaluated the magnitude and extent of the pressure field created by two water gun configurations in the concrete-walled channel downstream of the lock where channel depths ranged from 11 to 14 feet (ft). Data from a single 80-cubic-inch (in³) water gun set at 6 ft below water surface (bws) produced a roughly cylindrical 5-lb/in 2 pressure field 20 ft in radius, oriented vertically, with the radius decreasing to less than 15 ft at the water surface. A combination of two 80-in 3 water guns set at 6 and 8 ft, respectively, produced a similarly shaped 5 lb/in 2 pressure field 30 ft in radius. Neither of the water gun configurations exceeded the given threshold of 5 lb/in 2 above the static pressure along the walls of the canal at the 700 lb/in 2 water gun input pressure. Velocity and acceleration data were collected simultaneously with the underwater pressure data to understand the response of adjacent canal walls to the water gun firings. Maximum velocity and acceleration were 0.239 in/s and 0.0188 feet per second squared (ft/s 2 ), respectively. The Lemont study replicated and expanded upon work done in 2011. The pressure field created by the water gun was evaluated in a deeper environment (about 25 ft of water depth) than that of the Brandon Road study. To replicate the 2011 study, data were collected with the same water gun placements and input pressure, but static underwater pressure monitoring was added. Two 80-in 3 water guns were suspended below a platform at depths of 4 and 14 ft bws. Pressure was lower when the gun suspended at 4 ft bws was fired as compared to firing the single gun suspended at 14 ft bws. Firing both guns simultaneously produced similar pressures to the single gun suspended at 14 ft bws. Data were collected to assess the pressure field produced by two 80-in 3 water guns separated by 80 ft and suspended at a depth of 14 ft bws. The spatial extent of the 5-lb/in 2 threshold varied substantially with gun input air pressure. Firing the water gun with an air pressure of 2,000 lb/in 2 generated a pressure field greater than the threshold at all but one location in the measured region. Additionally, the water gun with an air pressure of 1,000 lb/in 2 did not reach the threshold anywhere in the measured region. Maximum velocity and acceleration were 0.304 in/s and 0.015 ft/s 2 , respectively.

Illinois↗

USGS2018

Simulation of groundwater flow, 1895–2010, and effects of additional groundwater withdrawals on future stream base flow in the Elkhorn and Loup River Basins, central Nebraska—Phase three

The U.S. Geological Survey, in cooperation with the Lewis and Clark, Lower Elkhorn, Lower Loup, Lower Platte North, Lower Niobrara, Middle Niobrara, Upper Elkhorn, and the Upper Loup Natural Resources Districts, designed a study to refine the spatial and temporal discretization of a previously modeled area. This updated study focused on a 30,000-square-mile area of the High Plains aquifer and constructed regional groundwater-flow models to evaluate the effects of groundwater withdrawal on stream base flow in the Elkhorn and Loup River Basins, Nebraska. The model was calibrated to match groundwater-level and base-flow data from the stream-aquifer system from pre-1940 through 2010 (including predevelopment [pre-1895], early development [1895–1940], and historical development [1940 through 2010] conditions) using an automated parameter-estimation method. The calibrated model then was used to simulate hypothetical development conditions (2011 through 2060). Predicted changes to stream base flow based on simulated changes to groundwater withdrawal will aid in developing strategies for management of hydrologically connected water supplies. Additional wells were simulated throughout the model domain and pumped for 50 years to assess the effect of wells on aquifer depletions, including stream base flow. The percentage of withdrawal for each well after 50 years, which was compensated by aquifer reductions to stream base flow, storage, or evapotranspiration, was computed and mapped. These depletions are influenced by aquifer properties, time, and distance from the well. Stream base-flow depletion results showed that the closer the added well was to a stream, the greatest the effect on the stream base flow. Areas of stream base-flow depletion percentages greater than 80 percent were generally within 1 mile (mi) from the stream. The distance increased to 6 mi near the confluence of the Dismal and Middle Loup Rivers, and the North Loup and Calamus Rivers. The percentage of stream base-flow depletion decreased as the distance from the stream increased. Areas more than 10 mi from the stream generally had a stream base-flow depletion of 10 percent or less. Evapotranspiration depletion was largest in areas closest to streams, specifically in the Elkhorn River watershed. It was also larger in areas of interdunal wetlands within the Sand Hills. Evapotranspiration depletion was negligible in areas greater than 5 mi from a stream, with the exception of interdunal areas in Cherry, Grant, and Arthur Counties. The storage depletion percentage increased as the distance from a stream increased. Storage depletion was largest in areas between streams. Areas experiencing the smallest amount of storage depletion were adjacent to streams. Calibrated model outputs and streamflow depletion analysis are publicly available online. Accuracy of the simulations is affected by input data limitations, system simplifications, assumptions, and resources available at the time of the simulation construction and calibration. Most of the important limitations relate either to data used as simulation inputs or to data used to estimate simulation inputs. Development of the regional simulations focused on generalized hydrogeologic characteristics within the study area and did not attempt to describe variations important to local-scale conditions. These simulations are most appropriate for analyzing groundwater-management scenarios for large areas and during long periods and are not suitable for analysis of small areas or short periods.

Nebraska↗

USGS2019

Water-balance modeling of selected lakes for evaluating viability as long-term fisheries in Kidder, Logan, and Stutsman Counties, North Dakota

Water levels in lakes and wetlands in the central North Dakota Missouri Coteau region that were either dry or only sporadically held water since before the 1930s have been rising since the early 1990s in response to an extended wet period. The lakes have remained full since the mid-1990s, which has provided benefits to migratory waterfowl, fisheries, and wildlife. A small shift in climate conditions, either to drier or wetter conditions, can have a large effect on the lake levels of these water bodies. The North Dakota Game and Fish Department identified five lakes as candidates for sustaining long-term fisheries. The lakes are in Kidder, Stutsman, and Logan Counties, and some lakes might receive inflow from mostly freshwater aquifers, such as the Central Dakota and Streeter aquifers, and were mostly dry during the early 1990s. After about 1995, the lakes had filled up and were deep enough to sustain populations of game fish such as walleye, perch, and northern pike. Before investing in development of permanent fisheries and associated infrastructure, such as campgrounds and boat ramps, fisheries biologists needed to know if the lake levels are likely to remain high in coming decades. The U.S. Geological Survey, in cooperation with the North Dakota Game and Fish Department, developed a water-balance model to determine the effects of precipitation, evapotranspiration, and groundwater interaction on lake volumes. The model was developed using climate input data and lake volumes for the calibration period 1992 through 2016, during which historical lake volumes could be estimated using land surface elevation data and Landsat images. Long-term (1940–2018) climate input data were used with the water-balance model to reconstruct historical lake volumes prior to the calibration period, and block-bootstrapping was used to simulate potential future climate input data and lake volumes for 2017 through 2067. The simulated future lake volumes were used to estimate the likelihood of annual lake volumes remaining consistent, increasing, or decreasing through the year 2067. Of the five lakes, Sibley Lake was the most likely to sustain a long-term fishery for a period longer than 50 years. The simulated lake volumes for Alkaline Lake, Big Mallard Marsh, and Remmick Lake indicated the lakes have a 50-percent chance to fall below 75 percent of their 2016 volume by about 2030, 2067, and 2025, respectively. Simulation results for Marvin Miller Lake were substantially different compared to the other four lakes and indicated the lake has a 50-percent chance to fall below 75 percent of its 2016 volume prior to 2025.

North Dakota↗

USGS2019

Pecos River Basin salinity assessment, Santa Rosa Lake, New Mexico, to the confluence of the Pecos River and the Rio Grande, Texas, 2015

The elevated salinity of the Pecos River throughout much of its length is of paramount concern to water users and water managers. Dissolved-solids concentrations in the Pecos River exceed 3,000 milligrams per liter in many of its reaches in the study area, from Santa Rosa Lake, New Mexico, to the confluence of the Pecos River with the Rio Grande, Texas. The salinity of the Pecos River increases downstream and affects the availability of useable water in the Pecos River Basin. In this report, “salinity” and “dissolved-solids concentration” are considered synonymous; both terms are used to refer to the total ionic concentration of dissolved minerals in water. The sources of salinity in the Pecos River Basin are natural (geologic) and anthropogenic, including but not limited to groundwater discharge, springs, and irrigation return flows. Previous studies in the Pecos River Basin were project specific and designed to address salinity issues in specific parts of the basin; therefore, in 2015, the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, New Mexico Interstate Stream Commission, Texas Commission on Environmental Quality, and Texas Water Development Board assessed the major sources of salinity throughout the extent of the basin where elevated salinity in the Pecos River is well documented (that is, in the drainage area of the Pecos River from Santa Rosa Lake to the confluence of the Pecos River and the Rio Grande). The goal was to gain a better understanding of how specific areas might be contributing to the elevated salinity in the Pecos River and how salinity of the Pecos River has changed over time. This assessment includes a literature review and compilation of previously published salinity-related data, which guided the collection of additional water-quality samples and streamflow gain-loss measurements. Differences in water quality of surface-water and groundwater samples, streamflow measurements, and geophysical data were assessed to gain new insights regarding sources of salinity in the Pecos River Basin and a more detailed assessment of potential areas of elevated salinity in the basin. The datasets compiled for this assessment are available in a companion data release. The literature review identified several potential sources of salinity inputs to the Pecos River in New Mexico and Texas. In New Mexico, sources of salinity inputs included sinkhole springs discharging into El Rito Creek, the Bitter Lake National Wildlife Refuge inflow to the Pecos River, inflow from the Rio Hondo, including the main channel and a restored channel at the Bitter Lake National Wildlife Refuge referred to as the “Rio Hondo spring channel,” the outflow from Lea Lake at Bottomless Lakes State Park, and the Malaga Bend region of the Pecos River. In Texas, sources of salinity inputs included Salt Creek downstream from Red Bluff Reservoir and the area near the Horsehead Crossing ford on the Pecos River. The compilation of historical water-quality data revealed a lack of consistent sampling of the same constituents at the same sites along the main stem of the Pecos River, which results in data gaps that hinder the ability to effectively analyze long-term changes in water quality that may help with the understanding of how salinity in the Pecos River has changed over time and identifying the sources of salinity in the Pecos River Basin. To help fill these data gaps, water-quality and streamflow data were collected in the study area in February 2015 by the U.S. Geological Survey. Historical water-quality data and newly collected data from February 2015 were evaluated for selected major-ion concentrations, dissolved-solids concentrations, and deuterium, oxygen, and strontium isotopes. Analysis of the data indicated several areas of increasing salinity in the Pecos River. Most notable increases were in two subreaches of the river, between Acme, N. Mex., and Artesia, N. Mex., and between Orla, Tex., and Grandfalls, Tex. Increasing sodium and chloride concentrations from Acme to Artesia coincided with changes in isotopic ratios within the Pecos River Basin. Changes in isotopic ratios in this reach indicate a likely inflow from an isotopically different source of water compared to the water in the main stem of the Pecos River, such as groundwater inflow, inflow from surface-water features distinct from the main stem of the Pecos River, or both. In the subreach between Orla and Grandfalls, an increase in dissolved-solids concentrations was observed along with a shift in isotope values, indicating that neither evaporative processes in Red Bluff Reservoir nor inflow from Salt Creek likely solely influences the salinity of the Pecos River in this subreach. The highest dissolved-solids concentrations in the Pecos River Basin were measured downstream from Grandfalls, where dissolved-solids concentrations are greater than 16,000 milligrams per liter near Iraan, Tex. Changes in isotopic values (deuterium, oxygen, and strontium) indicate mixing of different waters at several areas along the main stem of the Pecos River. The spatial distribution of the areas of interest from the literature review and the water-quality data are available in the companion data release.

Texas, New Mexico↗

USGS2020

Spatially referenced models of streamflow and nitrogen, phosphorus, and suspended-sediment loads in the southeastern United States

Spatially Referenced Regression On Watershed attributes (SPARROW) models were applied to describe and estimate mean-annual streamflow and transport of total nitrogen (TN), total phosphorus (TP), and suspended sediment (SS) in streams and delivered to coastal waters of the southeastern United States on the basis of inputs and management practices centered near 2012, the base year of the model. Previously published TN and TP models for 2002 served as a starting point and reference for comparison. The datasets developed for the 2012 models not only represent updates of previous conditions but also incorporate new approaches for characterizing sources and transport processes that were not available for previous models. Variability in streamflow across the southeastern United States was explained as a function of precipitation adjusted for evapotranspiration, spring discharge, and municipal and domestic wastewater discharges to streams. Results from the streamflow model were used as input to the water-quality SPARROW models, and areas with large streamflow prediction errors—urban areas and karst areas—were used to provide guidance on where additional data are needed to improve routing of flow. Variability in TN transport in Southeast streams was explained by the following five sources in order of decreasing mass contribution to streams: atmospheric deposition, agricultural fertilizer, municipal wastewater, manure from livestock, and urban land. Variable rates of TN delivery from source to stream were attributed to variation among catchments in climate, soil texture, and vegetative cover, including the extent of cover crops in the watershed. Variability in TP transport in Southeast streams was explained by the following six sources in order of decreasing mass contribution to streams: parent-rock minerals, urban land, manure from livestock, municipal wastewater, agricultural fertilizer, and phosphate mining. Varying rates of TP delivery were attributed to variation in climate, soil erodibility, depth to water table, and the extent of conservation tillage practices in the watershed. Variability in SS transport in Southeast streams was explained by variable sediment export rates for different combinations of land cover and geologic setting (for upland sources of sediment) and by gains in stream power caused by longitudinal changes in channel hydraulics (for channel sources of sediment). Sediment yields for the transitional land cover (shrub, scrub, herbaceous, and barren) varied widely depending on geologic setting and on agricultural land cover. Varying rates of SS delivery, like those for TP, were attributed to variation in climate, soil erodibility, and the extent of conservation tillage practices in the watershed, as well as to areal extent of canopy land cover in the 100-meter buffer along the channel. Relatively large uncertainty, compared to the other three models, for almost all the SS source coefficients indicates the need for caution when interpreting the results from the sediment model. TN, TP, and SS inputs to streams from sources were balanced in the models with losses from physical processes in streams and reservoirs and with water withdrawals. The losses in streams and reservoirs along with withdrawals removed 35, 44, and 65 percent of the TN, TP, and SS load, respectively, that entered streams before reaching coastal waters.

Scientific Investigations Report↗

USGS2021

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗

USGS2023

Simulation of flow and eutrophication in the central Salem River, New Jersey

The central Salem River in New Jersey is subject to periods of water-quality impairment, marked by elevated concentrations of phosphorus and chlorophyll- a , and low concentrations of and large diurnal swings in concentrations of dissolved oxygen. These seasonal eutrophic conditions are controlling factors for water quality in lower reaches, where the river is more lacustrine than in upper reaches, as a result of downstream damming. This biological productivity is supported by nutrient wash-off from agricultural areas in the surrounding watershed. To investigate this impairment, flow measurement and water-quality sampling were conducted during 2007–08 in support of development of a one-dimensional surface-water-quality model that simulates nutrient cycling and transformation processes. The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection, used the U.S. Environmental Protection Agency Water Quality Analysis Simulation Program (WASP) to develop a receiving-water-quality model of the central Salem River between Woodstown and Deepwater, New Jersey, from April 2007 to October 2008. The main-stem river and largest tributary were simulated. In the flow model, kinematic wave flow is used to simulate flow in upper reaches and ponded weir flow is used to simulate flow in lower reaches. The water-quality model makes use of a mass-balance equation to simulate the fate and transport of nutrients, phytoplankton chlorophyll- a , dissolved oxygen, and oxygen demands (an indicator rather than a substance) in the river. Model input included channel characteristics, boundary conditions for flow and water quality, environmental parameters, vertical dispersion coefficients, settling rates, and kinetic constants. Inputs were estimated where field data were lacking, notably for tributary flows and nutrient loads. The model was calibrated to observed flow variables and concentrations of dissolved oxygen, chlorophyll- a , and nutrients at sampling locations, with emphasis on growing-season conditions. Calibration was achieved through graphical and statistical comparison of simulated results to observed data. Sensitivity analyses were performed, and model limitations and applicability were evaluated. Simulated results closely matched observed data in most cases, although some were overpredicted slightly. The most important causes of overprediction were estimated tributary flows for the flow model and estimated tributary watershed loads for the water-quality model. Calibration of dissolved-oxygen concentrations was closer, and predicted diurnal variations were consistent with high algal photosynthesis/respiration, although lack of continuous dissolved-oxygen data precluded verifying these predictions. A similar caveat applies to predicted diurnal variations in chlorophyll- a . Simulated limitations on algal growth were consistent with those based on observed data and indicated phosphorus was the main limiting nutrient, except during certain periods when nitrogen was limiting. Two water-quality management scenarios were simulated with the model to assess the effect of point- and nonpoint-source nutrient reductions on water-quality conditions in the river. Scenarios involved (1) a return of watershed land use to predevelopment natural conditions and (2) an extreme reduction in nutrient input. Although the extreme-nutrient-reduction scenario yielded improvements in water quality, the natural-conditions scenario yielded the largest improvements as indicated by minimal violations of surface-water-quality standards or thresholds. However, years may be needed to attain the full benefit of these management scenarios as a result of accumulation of phosphorus and organic carbon in riverbed sediments in lacustrine reaches. The results of this study indicate that the quality of water in the central Salem River will improve if management policies that mitigate the effects of nutrient-loading practices in the watershed, particularly those related to agriculture, are implemented.

New Jersey↗

USGS2022

Effect of uncertainty of discharge data on uncertainty of discharge simulation for the Lake Michigan Diversion, northeastern Illinois and northwestern Indiana

Simulation models of watershed hydrology (also referred to as “rainfall-runoff models”) are calibrated to the best available streamflow data, which are typically published discharge time series at the outlet of the watershed. Even after calibration, the model generally cannot replicate the published discharges because of simplifications of the physical system embedded in the model structure and uncertainties of the input data and of the estimated model parameters, which, although optimized for the given calibration data, remain uncertain. The input data errors are caused by uncertainties in the forcing data, such as precipitation and other climatological data, and in the published discharges used for calibration. In the numerical algorithms used for calibration, the published discharges are often assumed to be without error, but they are themselves uncertain, typically having been computed using ratings, which are models fitted to uncertain discharge measurements. In this study, uncertainty of published daily discharge data and how the discharge uncertainty is transmitted to the parameter values of the Hydrological Simulation Program–FORTRAN (HSPF) rainfall-runoff model and to the simulated discharge at both calibration and prediction locations were investigated for the Lake Michigan diversion in northeastern Illinois and northwestern Indiana. The HSPF model used in this study is used by the U.S. Army Corps of Engineers as part of quantifying the diversion of water from Lake Michigan by the State of Illinois. In this study, the model is calibrated jointly at two watersheds in the study area; the resulting model is considered the base model in this study. Seven other gaged watersheds in the study area are used for testing predictive simulations. A Bayesian rating curve estimation (BaRatin) approach, the BaRatin stage-period-discharge (SPD) method, was used to estimate the uncertainty of the published discharge from the calibration watersheds. To characterize the effect of the discharge uncertainty on parameter values, the HSPF model parameters were recalibrated to 17 nonrandomly selected pairs of discharge series from the BaRatin SPD analysis. To provide an indicator of the effect of parameter uncertainty to compare to the effect of discharge uncertainty, 1,000 parameter sets also were randomly generated from the estimated parameter covariance matrix of the base model. The recalibrated and random parameter sets were then used in HSPF simulations of discharge at the two calibration watersheds and at the seven prediction watersheds. Selected discharge summary statistics—the period-of-study (POS, water years 1997 to 2015) mean discharge, selected flow-duration curve (FDC) quantiles, and water year mean discharges—are used to characterize the variability between simulated and published discharge. A normalized variability index ( V N ) is used as a measure of the uncertainty of flow statistics arising from the uncertainty of the sources considered in this study. When this index is at least 1, the variability of the simulations is large enough to explain the median error between simulated and published values, although offsetting errors from other sources are also likely. When the index is appreciably less than 1, the variability of the simulations is clearly insufficient to explain the median error between simulated and published values. At the two calibration watersheds and for results of the two simulation sets considered together, the V N values ranged from 0.2 to 0.8 for POS mean discharge, from 0.3 to 0.6 in the median for a set of FDC quantiles, and from 0.1 to 0.2 in the median for water year mean discharges. These values indicate that substantial uncertainty remains unexplained. Even though two watersheds were used in calibration, that calibration was highly constrained because it was applied to the watersheds simultaneously and was subject to parameter regularization that constrained the adjustment of the parameters from their initial values. These constraints were applied to avoid overfitting to the calibration watersheds and thus to increase the likelihood that the resulting parameters would give accurate results at watersheds not used in the calibration, but they created a parameter transfer error in the calibration watershed results shown by the balancing of errors between the two watersheds. Additional remaining error sources include model structural error and meteorological forcing error to the degree that the calibration was unable to adjust the parameters to account for these errors. At the prediction watersheds, the corresponding V N values were almost always substantially lower than those values at the calibration watersheds. This result is expected because the prediction watersheds have additional uncertainty, including parameter transfer error. The work described in this report provides preliminary estimates of a limited range of sources of error in predicted discharge uncertainty. Future work would be beneficial to obtain a better statistical characterization of the effect of the uncertainty of calibration discharge series and to address additional sources of uncertainty, such as from precipitation input data used in calibration and prediction and from structural (model) errors.

Illinois, Indiana↗

USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗