SEARCH · Search USGS

Results for “Scientific Investigations Report”

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 1,621 records · Page 90Linked to original sources

USGS2006

Application of a stream-aquifer model to Monument Creek for development of a method to estimate transit losses for reusable water, El Paso County, Colorado

The U.S. Geological Survey, in cooperation with Colorado Springs Utilities, the Colorado Water Conservation Board, and the El Paso County Water Authority, began a study in 2004 to (1) apply a stream-aquifer model to Monument Creek, (2) use the results of the modeling to develop a transit-loss accounting program for Monument Creek, (3) revise the existing transit-loss accounting program for Fountain Creek to incorporate new water-management strategies and allow for incorporation of future changes in water-management strategies, and (4) integrate the two accounting programs into a single program with a Web-based user interface. The purpose of this report is to present the results of applying a stream-aquifer model to the Monument Creek study reach. Transit losses were estimated for reusable-water flows in Monument Creek that ranged from 1 to 200 cubic feet per second (ft 3 /s) and for native streamflows that ranged from 0 to 1,000 ft 3 /s. Transit losses were estimated for bank-storage, channel-storage, and evaporative losses. The same stream-aquifer model used in the previously completed (1988) Fountain Creek study was used in the Monument Creek study. Sixteen model nodes were established for the Monument Creek study reach, defining 15 subreaches. Channel length, aquifer length, and aquifer width for the subreaches were estimated from available topographic and geologic maps. Thickness of alluvial deposits and saturated thickness were estimated using lithologic and water-level data from about 100 wells and test holes in or near the Monument Creek study reach. Estimated average transmissivities for the subreaches ranged from 2,000 to 12,000 feet squared per day, and a uniform value of 0.20 was used for storage coefficient. Qualitative comparison of recorded and simulated streamflow at the downstream node for the calibration and verification simulations indicated that the two streamflows compared reasonably well. No adjustments were made to the model parameters. Differences between recorded and simulated streamflow volumes for all calibration and verification simulations ranged from about –8.8 to 7.5 percent; the total error for all simulations was about –0.7 percent. The model was used to estimate bank-storage losses for 10 to 15 native streamflows for each reusable-water flow of 1, 3, 5, 7, 10, 15, 20, 30, 40, 50, 100, and 200 ft 3 /s. Then the 10 to 15 bank-storage loss values were used in least-squares linear regression to estimate a relation between bank-storage loss and native streamflow for each of the 12 reusable-water flow rates. The 12 regression relations then were used to develop “look-up” tables of bank-storage loss for reusable-water flows ranging from 1 to 200 ft 3 /s (in 1-ft 3 /s increments). Additional model simulations indicated that (1) when the ratio of downstream native streamflow to upstream native streamflow was less than 1, bank-storage loss generally increased and (2) when the ratio of downstream native streamflow to upstream native streamflow was larger than 1, bank-storage loss generally decreased. These results were used to develop a bank-storage loss adjustment factor based on the ratio of native streamflow at the downstream node to native streamflow at the upstream node. The model also was used to estimate a recovery period, which is the length of time needed for the bank-storage loss to return to the stream. The recovery period was 1 day for six subreaches; 2 days for four subreaches; between 3 and 12 days for four subreaches; and 28 days for one subreach. Channel-storage losses are about 10 percent of the reusable-water flow for most of the subreaches, except for two subreaches, where the channel-storage losses are about 20 percent, and one subreach, where the losses are about 30 percent, owing to the greater channel lengths. Evaporative losses were estimated by the use of monthly pan-evaporation data and the incremental increase in stream width resulting from any reusable-water flows. Monthly pan-evaporation data were converted to a daily rate. The daily rate, when multiplied by the stream-width increase (in feet) that results from reusable-water flow and by the subreach length (in miles) gives the daily evaporative loss in cubic feet per second.

Colorado↗

USGS2009

Variations in Withdrawal, Return Flow, and Consumptive Use of Water in Ohio and Indiana, with Selected Data from Wisconsin, 1999-2004

This report contains an analysis of water withdrawal and return-flow data for Ohio and withdrawal data for Indiana and Wisconsin to compute consumptive-use coefficients and to describe monthly variability of withdrawals and consumptive use. Concurrent data were available for most water-use categories from 1999 through 2004. Average monthly water withdrawals are discussed for a variety of water-use categories, and average water use per month is depicted graphically for Ohio, Indiana, and Wisconsin (public supply only). For most water-use categories, the summer months were those of highest withdrawal and highest consumptive use. For public supply, average monthly withdrawals ranged from 1,380 million gallons per day (Mgal/d) (November) to 1,620 Mgal/d (July) in Ohio, 621 Mgal/d (December) to 816 Mgal/d (July) in Indiana, and 515 Mgal/d (December) to 694 Mgal/d (July) in Wisconsin. Ohio and Indiana thermoelectric facilities had large increases in average monthly withdrawals in the summer months (5,520 Mgal/d in March to 7,510 Mgal/d in August for Indiana; 7,380 Mgal/d in February to 10,040 Mgal/d in July for Ohio), possibly because of increased electricity production in the summer, a need for additional cooling-water withdrawals when intake-water temperature is high, or use of different types of cooling methods during different times of the year. Average industrial withdrawals ranged from 2,220 Mgal/d (December) to 2,620 Mgal/d (August) in Indiana and from 707 Mgal/d (January) to 787 Mgal/d (August) in Ohio. The Ohio and Indiana irrigation data showed that most withdrawals were in May through October for golf courses, nurseries, and crop irrigation. Commercial water withdrawals ranged from 30.4 Mgal/d (January) to 65.0 Mgal/d (September) in Indiana and from 23.2 Mgal/d (November) to 49.5 Mgal/d (August) in Ohio; commercial facilities that have high water demand in Ohio and Indiana are medical facilities, schools, amusement facilities, wildlife facilities, large stores, colleges, correctional institutions, and national security facilities. Monthly livestock withdrawals were constant for Ohio but were more variable in Indiana and depended on whether the livestock facility operated on a seasonal schedule. Aquaculture withdrawals appeared to correlate with growing seasons and with aeration of ponds during the winter months. Mining withdrawals - specifically, those for nonmetallic mining - tended to be highest in April and may be related to dewatering. Consumptive use and consumptive-use coefficients were computed by two principal methods in this study: the return-flow and withdrawal method (RW; Ohio only) and the winter-base-rate method (WBR; Ohio, Indiana and Wisconsin). The WBR method was not suitable for the thermoelectric, industrial, irrigation, livestock, aquaculture, and mining water-use categories. The RW method was not used for public-supply facilities. A third method, the Standard Industrial Classification code method (SIC), was used only for certain industrial facilities. The public-supply annual average consumptive-use coefficient derived by use of the WBR methods ranged from 6 to 8 percent among Ohio, Indiana, and Wisconsin; the summer average consumptive-use coefficient was considerably higher, ranging from 16 to 20 percent. The commercial annual consumptive-use coefficient for both Ohio and Indiana was 30 percent by the WBR method, which fell within the Ohio annual median (17 percent) and annual average (42 percent) by the RW method. Thermoelectric consumptive use differs greatly by the type of cooling the facility uses; the Ohio annual median consumptive-use coefficient (RW method) was 2 percent for all thermoelectric facilities and facilities with multiple types of cooling, but exclusively once-through-cooling facilities had a median of 0 percent and exclusively closed-loop-cooling facilities had a median of 25 percent. Industrial consumptive-use coefficients varied by type of industry, as reflected by SIC code

Scientific Investigations Report↗

USGS2021

Evaluation of actual evapotranspiration rates from the Operational Simplified Surface Energy Balance (SSEBop) model in Florida and parts of Alabama and Georgia, 2000–17

Evapotranspiration (ET) is the water-vapor flux transported from the surface of the Earth into the atmosphere and is the sum of surface water directly evaporated and subsurface water transpired by plants. ET rates are commonly estimated by using potential or reference ET, which might differ from actual ET rates. Actual evapotranspiration (ETa) rates can be estimated by using the Operational Simplified Surface Energy Balance (SSEBop) model. This report evaluates SSEBop ETa rates at the point and basin scales in Florida and parts of Alabama and Georgia for 2000–17. ETa rates computed by using data from 24 micrometeorological stations in Florida are referred to as mETa rates and were used to quantify biases in the SSEBop ETa rates, stratified by generalized land-use type. Bias was computed as mETa minus SSEBop ETa rates for given generalized land-use types, and bias-correction equations were computed by using least-squares regressions. In addition to mETa rates at station locations, annual average ETa rates calculated from the application of a water-balance method to 55 basins in Florida and parts of Alabama and Georgia were used to assess the accuracy of the annual SSEBop ETa rates at the basin scale. Another independent model used to simulate ETa rates was based on monthly reference ET from the statewide daily reference evapotranspiration (ETo) gridded dataset for Florida computed by using Geostationary Operational Environmental Satellite estimates of solar radiation (GOES ETo). ETa at grid points was computed as monthly GOES ETo multiplied by ratios of monthly mETa to GOES ETo, computed at micrometeorological stations and stratified by each generalized land-use type. The coefficient of determination (R 2 ) between monthly mETa and SSEBop ETa rates for all stations combined improved from 0.37 before bias correction of SSEBop ETa rates to 0.79 after the bias correction stratified by land-use type. For individual land-uses types, R 2 varied from 0.59 for the monthly mETa at a station in the land-use type forest to 0.82 for the monthly mETa at stations in the land-use type shallow-water-table pasture. Root-mean-square error (RMSE) was computed as a function of the difference between SSEBop ETa rates and mETa rates. RMSE of monthly SSEBop ETa rates was 1.27 inches per month before the bias corrections improved to 0.73 inch per month after the bias corrections. RMSE for bias-corrected annual SSEBop ETa rates based on micrometeorological stations with complete years of records ranged from 2.01 inches per year (in/yr) for the land-use type of agriculture to 5.73 in/yr for the land-use type of deep water-table pasture, or 4.96 and 21.21 percent errors relative to annual mETa rates, respectively. Bias-corrected annual SSEBop ETa rates were also compared to annual ETa rates computed by using a water-balance method (wbETa) for 55 basins in Florida. Differences in bias-corrected average annual SSEBop ETa rates and average annual wbETa rates for the 55 basins ranged from −3.67 to 5.29 in/yr (−9.24 to 17.36 percent). RMSE when computed as a function of the differences between annual SSEBop ETa rates and wbETa rates decreased, on average, from 4.13 in/yr for the uncorrected bias SSEBop ETa rates to 1.95 in/yr for the bias-corrected SSEBop rates. The average annual bias-corrected SSEBop ETa rates, from all basins, was 36.46 in/yr or 3.41 percent lower than the average annual wbETa rate of 37.79 inches. Bias in SSEBop ETa rates varies based on time step (monthly versus annual), scale (point, basin, statewide), and land-use type. Applications to hydrologic models should consider bias relative to the inherent error in models. Bias-corrected SSEBop ETa rates could be used as calibration targets in models of hydrologic processes, such as groundwater models. Annual bias in SSEBop ETa introduced to the model calibration is typically below the margin of error associated with typical residuals in model simulations, depending on scale. Surface-water and groundwater-flow models with RMSEs on the order of a few feet could benefit from bias-corrected SSEBop values of ETa.

Alabama, Florida, Georgia↗

USGS2024

Peak streamflow trends in Missouri and their relation to changes in climate, water years 1921–2020

This report characterizes changes in peak streamflow in Missouri and the relation of these changes to climatic variability, and provides a foundation for future studies that can address nonstationarity in peak-streamflow frequency analysis in Missouri. Records of annual peak and daily streamflow at streamgages and gridded monthly climatic data (observed and modeled) were examined across four trend periods (100 years, water years 1921–2020; 75 years, 1946–2020; 50 years, 1971–2020; and 30 years, 1991–2020) for trends, change points (abrupt changes in the streamflow time series), and other statistical properties indicative of changing conditions. Peak streamflow magnitudes generally exhibit upward trends across the State for the 100-, 75-, and 50-year trend periods and only in southern Missouri for the 30-year trend period. The medians of the trend magnitudes (normalized by median peak streamflow) range from a 10-percent increase during the 30-year trend period to a 40-percent increase during the 100-year trend period. Changes in the 90-percent quantile of peak streamflow, which correspond to the 10-percent exceedance probability often used for the design of drainage structures, are not as substantial or widespread, showing consistent increases mainly in the southern part of the State in the 50- and 30-year trend periods. Streamgages with trends in peak streamflow often also have change points, or abrupt changes, in streamflow magnitude. Change points in peak streamflows generally follow that of the peak streamflow trends, with upward change points throughout most of the State at the 100- and 75-year trend periods and in southern Missouri at the 30-year trend period. Temporally, clusters upward of change points are observed in the 1970s through 1980s for the 100-, 75-, and 50-year trend periods and around 2006 and 2007 for the 50- and 30-year trend periods. A peaks-over-threshold analysis, which evaluates changes in the frequency of peak streamflows over a certain threshold, indicates that high flows have increased in frequency at 50 to 64 percent of streamgages in the 100- and 75-year trend periods. Most streamgages in the 50- and 30-year trend periods exhibit no change. Although the frequency of high flows has increased at some streamgages and trend periods in Missouri, these increases are not as widespread as the increases in the magnitude of peak streamflow. Upward trends in observed temperature and observed annual precipitation dominate in all trend periods, with no downward trends in precipitation and only two somewhat likely downward trends in temperature for the 100-year trend period. Increases in annual precipitation mostly are limited to southern Missouri for the 30-year trend period. The proportion of precipitation falling as snow has largely decreased in the study basins across the State, which is expected in response to increasing temperature. Upward trends in modeled annual runoff, which in this study incorporates only the effects of climatic variation, are observed in the same geographic areas where there are increases in observed annual precipitation. When peak streamflow and climatic trends are considered together, widespread upward trends in peak streamflows for the 100-, 75-, and 50-year trend periods and for the 30-year trend period mainly in southern Missouri (encompassing both trends and abrupt change) appear to be driven largely by increases in precipitation based on spatial patterns and statistical relations. The prevalence of nonstationarity in peak streamflow in Missouri has important implications for peak-flow frequency analysis. Winter and spring precipitation and the occurrence of extreme precipitation events are expected to increase across the State. If precipitation continues to increase as expected, peak-flow frequency estimates based on older records may no longer represent the hydrologic regime of today, and methods for nonstationary peak-flow frequency analysis may be needed.

Missouri↗

USGS2020

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and northern Chihuahua, Mexico

Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model, and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model version 2 (referred to as MF-OWHM2). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 1,090 acre-feet per year.

Chihuahua, New Mexico, Texas↗

USGS2021

Analysis of Escherichia coli, total recoverable iron, and dissolved selenium concentrations, loading, and identifying data gaps for selected 303(d) listed streams, Grand Valley, western Colorado, 1980–2018

Tributaries to the Colorado River in the Grand Valley in western Colorado (segment COLCLC13b) have been placed on the State of Colorado 303(d) list as impaired for Escherichia coli (E. coli) , total recoverable iron, and dissolved selenium. The Colorado Department of Public Health and Environment Water Quality Control Division is required to develop total maximum daily loads for these constituents in these tributaries. The U.S. Geological Survey, in cooperation with the Grand Valley Drainage District and Colorado Water Conservation Board, conducted a study to (1) characterize concentrations, loads, and load reductions for E. coli , total recoverable iron, and dissolved selenium using existing data and (2) identify water-quality data gaps to inform future monitoring strategies. This study analyzed water-quality and streamflow data for 3 main-stem sites (2 sites along the Colorado River and 1 site along the Gunnison River) and 29 selected sites on tributaries to the Colorado River. Sample data were available at five sites along Adobe Creek and at six sites along Leach Creek, the two tributaries in the study area that are impaired for E. coli . All geometric mean E. coli concentrations at sites along Adobe Creek and Leach Creek exceeded the State recreational use standard of 126 colony forming units per 100 milliliters (CFU/100 mL). In Adobe Creek, E. coli concentrations in samples ranged from 45.7 to more than 2,420 CFU/100 mL (method upper reporting limit for undiluted samples), and geometric mean concentrations at sites ranged from 301 to 1,180 CFU/100 mL. The E. coli concentrations generally increased in the downstream direction in Adobe Creek; however, increases were not seen between all sites. The largest downstream increase in E. coli concentration was measured between the two most upstream sites. In Leach Creek, concentrations of E. coli in samples ranged from 25.9 to more than 2,420 CFU/100 mL, and geometric mean concentrations at sites ranged from 160 to 259 CFU/100 mL. The E. coli concentrations showed no consistent downgradient increase in Leach Creek. In fact, some of the highest E. coli concentrations were measured at the most upstream site, Leach Creek at Summer Hill Drive. Total recoverable iron concentrations and loads were evaluated at 15 tributary sites for samples collected from August 1993 to February 2018. Median total recoverable iron concentrations ranged from 211 to 4,670 micrograms per liter (µg/L). The chronic aquatic-life water-quality standard (1,000 µg/L) was exceeded in most irrigation season (April through October) samples but was rarely exceeded in nonirrigation season (November through March) samples. Concentrations were often an order of magnitude higher in samples collected during irrigation season than in samples collected during nonirrigation season. None of the sites had enough concurrent total recoverable iron and streamflow data to compute annual loads. As with E. coli , the lack of concurrent total recoverable iron and streamflow information represents a data gap, which needs to be addressed to compute annual loads. Dissolved selenium concentrations and loads were evaluated at 20 tributary sites using discrete water-quality data collected 1991–2018. Dissolved selenium concentrations were higher during nonirrigation season than during irrigation season at tributary sites. However, irrigation season dissolved selenium loads were generally higher than nonirrigation selenium loads, because streamflows were higher during irrigation season. Regression analysis was used to estimate daily dissolved selenium concentrations and loads at three main-stem sites for water years (WYs) 1980–2018 (Gunnison River near Grand Junction and Colorado River near Colorado-Utah State Line) and WYs 2002–18 (Colorado River near Cameo). A trend analysis of dissolved selenium concentrations and loads was completed for these sites from the same respective starting dates but ending in 2017. A continuing downward trend in dissolved selenium concentration was observed at all sites and across all seasonal designations of the analysis. The dissolved selenium concentration decreased by 0.12 µg/L from WY 2002 to 2017 at Colorado River near Cameo, representing an 18-percent decrease during the time period. The dissolved selenium concentration at Gunnison River near Grand Junction decreased by 4.2 µg/L from WY 1980 to 2017, representing a 56-percent decrease overall. During the same time period, dissolved selenium concentration at Colorado River near Colorado-Utah State Line decreased by 3.8 µg/L, representing a 56-percent decrease overall. A downward trend in dissolved selenium load was also observed at all sites and across all seasonal designations of the analysis. The relative contribution of dissolved selenium from the Grand Valley near Grand Junction was estimated by comparing loads at main-stem sites bracketing the study area. The two upstream sites, Colorado River near Cameo and Gunnison River near Grand Junction, contributed 60,300 cumulative pounds and 251,000 cumulative pounds, respectively, during WYs 2002–18. At the furthest downstream site, Colorado River near Colorado-Utah State Line, 490,000 cumulative pounds were estimated during the same time period, indicating that the region between Whitewater and State line contributed approximately 179,000 cumulative pounds or a mean annual load of 10,500 lb/yr. Grand Valley dissolved selenium contributions appear to be stable during WYs 2002–18.

Colorado↗

USGS2006

Rates of evapotranspiration, recharge from precipitation beneath selected areas of native vegetation, and streamflow gain and loss in Carson Valley, Douglas County, Nevada, and Alpine County, California

Rapid growth and development in Carson Valley is causing concern over the continued availability of water resources to sustain such growth into the future. A study to address concerns over water resources and to update estimates of water-budget components in Carson Valley was begun in 2003 by the U.S. Geological Survey, in cooperation with Douglas County, Nevada. This report summarizes micrometeorologic, soil-chloride, and streambed-temperature data collected in Carson Valley from April 2003 through November 2004. Using these data, estimates of rates of discharge by evapotranspiration (ET), rates of recharge from precipitation in areas of native vegetation on the eastern and northern sides of the valley, and rates of recharge and discharge from streamflow infiltration and seepage on the valley floor were calculated. These rates can be used to develop updated water budgets for Carson Valley and to evaluate potential effects of land- and water-use changes on the valley's water budget. Data from eight ET stations provided estimates of annual ET during water year 2004, the sixth consecutive year of a drought with average or below average precipitation since 1999. Estimated annual ET from flood-irrigated alfalfa where the water table was from 3 to 6 feet below land surface was 3.1 feet. A similar amount of ET, 3.0 feet, was estimated from flood-irrigated alfalfa where the water table was about 40 feet below land surface. Estimated annual ET from flood-irrigated pasture ranged from 2.8 to 3.2 feet where the water table ranged from 2 to 5 feet below land surface, and was 4.4 feet where the water table was within 2 feet from land surface. Annual ET estimated from nonirrigated pasture was 1.7 feet. Annual ET estimated from native vegetation was 1.9 feet from stands of rabbitbrush and greasewood near the northern end of the valley, and 1.5 feet from stands of native bitterbrush and sagebrush covering alluvial fans along the western side of the valley. Uncertainty in most ET estimates is about 12 percent, but ranged from +30 and +50 percent to -20 and -40 percent for nonirrigated pasture and native bitterbrush and sagebrush. Estimated rates for water year 2004 likely are less than those during years of average, or above average precipitation when the water table would be closer to land surface. Test holes drilled in areas of native vegetation on the northern and eastern sides of Carson Valley had high concentrations of soil chloride at depths ranging from 4 to 18 feet below land surface at six locations on the eastern side of the valley. The high chloride concentrations indicate that modern-day precipitation at the six locations does not percolate deeper than the root zone of native vegetation. Estimates of the time required to accumulate the measured amount of chloride to depths of about 30 feet below land surface at the six test holes ranged from about 3,000 to 12,000 years. Low concentrations of soil chloride in two test holes on the northern end of Carson Valley and in a test hole on the eastern side of Fish Spring Flat indicate that a small amount of recharge from modern-day precipitation is taking place. Estimated annual recharge from precipitation at the two locations was 0.03 and 0.04 foot on the northern end of the valley and 0.02 foot on the eastern side of Fish Spring Flat. Uncertainty in the estimated recharge rates was about ?0.01 foot. Estimates of the time required to accumulate the measured amount of chloride to depths of about 30 feet below land surface at the three test holes ranged from about 100 to 700 years. The two test holes near the northern end of the valley are in gravel and eolian sand deposits and recharge from precipitation may be taking place at similar rates in other areas with gravel and eolian sand deposits. Based on results from other test holes, recharge at the rate estimated for the test hole on the eastern side of Fish Spring Flat is not likely applicable to a large area. Data from 37 site

Scientific Investigations Report↗

USGS2019

Comparison of groundwater-model construction methods, representations of glacial geology, model designs, and groundwater-model flow simulations within Elkhart County, Indiana

Automated data-processing methods allow hydrologists to efficiently incorporate digital well-record datasets into the construction of hydrostratigraphic frameworks for groundwater-flow models. The method selected to construct the hydrostratigraphic framework can affect the extent of geologic heterogeneity that can be included in the model. The detail generated from a hydrostratigraphic framework can affect groundwater simulation results. The effects of detail on model accuracy, groundwater-flow simulations, and particle-tracking simulations are described in this study. This report compares differences in hydrostratigraphic frameworks and results of groundwater models using (1) a method that incorporates more hydrologic judgment at the expense of using limited lithologic data and (2) a method that is more automated and uses all available lithologic data. The study additionally evaluates the effect of model discretization and inclusion of more (or less) geologic detail on simulation results. Two methods were used to create hydrostratigraphic frameworks of glacial deposits in the St. Joseph River Basin. One method, referred to as the subjective method, manually identifies stratigraphic boundaries using a sample of well logs from State databases and uses two-dimensional kriging to create three model layers of the study area. Indicator kriging is used to define aquifer extent in each layer. The second method, referred to as the objective method, uses three-dimensional kriging to automatically create a detailed heterogeneous model of the study area using all wells logs from the State database. The objective method increases detail in the vertical by greatly increasing the number of computer groundwater model layers from 3 to 30. In Elkhart County, Indiana, a previously published model represents the product of the subjective method, and a newly calibrated model of the same area represents the product of the objective method. An automated calibration procedure was used with the objective model (derived from the objective method) for Elkhart County. The two most-sensitive parameters for the Elkhart County objective model are horizontal hydraulic conductivity of the sand and the combined sand and gravel/gravel deposits. Vertical hydraulic conductivity of the fine-grained and intermediate-sized deposits could not be estimated, possibly indicating major flow paths are along a continuously connected series of sand and gravel deposits and not through a confining layer. The statistics measuring model calibration accuracy for the objective model were slightly better than statistics for the subjective model (model derived from the subjective method) of Elkhart County, but the hydraulic conductivities and flow rates for the two models were different. The mean absolute errors between simulated and measured groundwater levels are 2.04 and 2.16 feet for the objective and subjective models, respectively. Simulated seepage losses from and groundwater discharges to measured stream reaches in the objective model were evenly balanced in terms of over and under simulations of measured values; the subjective model tended to overpredict measured groundwater discharge to streams. The overprediction may be related to the 58 percent greater total inflow and outflow through the subjective model. The greater flow rate through the subjective model results from higher horizontal hydraulic conductivities in the subjective model than in the objective model. Horizontal hydraulic conductivity ranged from 23.9 to 111 feet per day in the objective model and generally ranged from 170 to 370 feet per day in the subjective model. The improvement in calibration statistics for the objective model relative to the subjective model may be from increased detail in how the objective model represents the distribution of fine- and coarse-grained deposits. The improvement also could be associated with the difference in methods used to represent the continuity of the confining unit. The effect of differences in horizontal hydraulic conductivity distributions between the two models for Elkhart County is evident in the groundwater-flow paths simulated by the objective and subjective models. At a withdrawal well location, the flow lines produced by the objective model indicate a wider contributing area than that for the subjective model. The discontinuous confining unit represented in the objective model provided the opportunity for groundwater flow to split into an upper and lower path. The split in flow simulated by the objective model at one location was independently supported by bromide concentrations in groundwater; the subjective model did not duplicate the split in flow.

Indiana↗

USGS2025

Quality of groundwater used for domestic supply in the Gilroy-Hollister basin and surrounding areas, California, 2022

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap. To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs. Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface. Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

California↗

USGS2008

Environmental factors affecting mercury in Camp Far West Reservoir, California, 2001-03

This report documents water quality in Camp Far West Reservoir from October 2001 through August 2003. The reservoir, located at approximately 300 feet above sea level in the foothills of the northwestern Sierra Nevada, California, is a monomictic lake characterized by extreme drawdown in the late summer and fall. Thermal stratification in summer and fall is coupled with anoxic conditions in the hypolimnion. Water-quality sampling was done at approximately 3-month intervals on eight occasions at several stations in the reservoir, including a group of three stations along a flow path in the reservoir: an upstream station in the Bear River arm (principal tributary), a mid-reservoir station in the thalweg (prereservoir river channel), and a station in the deepest part of the reservoir, in the thalweg near Camp Far West Dam. Stations in other tributary arms of the reservoir included those in the Rock Creek arm of the reservoir, a relatively low-flow tributary, and the Dairy Farm arm, a small tributary that receives acidic, metal-rich drainage seasonally from the inactive Dairy Farm Mine, which produced copper, zinc, and gold from underground workings and a surface pit. Several water-quality constituents varied significantly by season at all sampling stations, including major cations and anions, total mercury (filtered and unfiltered samples), nitrogen (ammonia plus organic), and total phosphorus. A strong seasonal signal also was observed for the sulfurisotope composition of aqueous sulfate from filtered water. Although there were some spatial differences in water quality, the seasonal variations were more profound. Concentrations of total mercury (filtered and unfiltered water) were highest during fall and winter; these concentrations decreased at most stations during spring and summer. Anoxic conditions developed in deep parts of the reservoir during summer and fall in association with thermal stratification. The highest concentrations of methylmercury in unfiltered water were observed in samples collected during summer from deepwater stations in the anoxic hypolimnion. In the shallow (less than 14 meters depth) oxic epilimnion, concentrations of methylmercury in unfiltered water were highest during the spring and lowest during the fall. The ratio of methylmercury to total mercury (MeHg/HgT) increased systematically from winter to spring to summer, largely in response to the progressive seasonal decrease in total mercury concentrations, but also to some extent because of increases in MeHg concentrations during summer. Water-quality data for Camp Far West Reservoir are used in conjunction with data from linked studies of sediment and biota to develop and refine a conceptual model for mercury methylation and bioaccumulation in the reservoir and the lower Bear River watershed. It is hypothesized that MeHg is produced by sulfate-reducing bacteria in the anoxic parts of the water column and in shallow bed sediment. Conditions were optimal for this process during late summer and fall. Previous work has indicated that Camp Far West Reservoir is a phosphate-limited system - molar ratios of inorganic nitrogen to inorganic phosphorus in filtered water were consistently greater than 16 (the Redfield ratio), sometimes by orders of magnitude. Therefore, concentrations of orthophosphate were expectedly very low or below detection at all stations during all seasons. It is further hypothesized that iron-reducing bacteria facilitate release of phosphorus from iron-rich sediments during summer and early fall, stimulating phytoplankton growth in the fall and winter, and that the MeHg produced in the hypolimnion and metalimnion is released to the entire water column in the late fall during reservoir destratification (vertical mixing). Mercury bioaccumulation factors (BAF) were computed using data from linked studies of biota spanning a range of trophic position: zooplankton, midge larvae, mayfly nymphs, crayfish, threadfin shad, bluegill,

California↗

USGS2016

Estimating flood magnitude and frequency at gaged and ungaged sites on streams in Alaska and conterminous basins in Canada, based on data through water year 2012

Estimates of the magnitude and frequency of floods are needed across Alaska for engineering design of transportation and water-conveyance structures, flood-insurance studies, flood-plain management, and other water-resource purposes. This report updates methods for estimating flood magnitude and frequency in Alaska and conterminous basins in Canada. Annual peak-flow data through water year 2012 were compiled from 387 streamgages on unregulated streams with at least 10 years of record. Flood-frequency estimates were computed for each streamgage using the Expected Moments Algorithm to fit a Pearson Type III distribution to the logarithms of annual peak flows. A multiple Grubbs-Beck test was used to identify potentially influential low floods in the time series of peak flows for censoring in the flood frequency analysis. For two new regional skew areas, flood-frequency estimates using station skew were computed for stations with at least 25 years of record for use in a Bayesian least-squares regression analysis to determine a regional skew value. The consideration of basin characteristics as explanatory variables for regional skew resulted in improvements in precision too small to warrant the additional model complexity, and a constant model was adopted. Regional Skew Area 1 in eastern-central Alaska had a regional skew of 0.54 and an average variance of prediction of 0.45, corresponding to an effective record length of 22 years. Regional Skew Area 2, encompassing coastal areas bordering the Gulf of Alaska, had a regional skew of 0.18 and an average variance of prediction of 0.12, corresponding to an effective record length of 59 years. Station flood-frequency estimates for study sites in regional skew areas were then recomputed using a weighted skew incorporating the station skew and regional skew. In a new regional skew exclusion area outside the regional skew areas, the density of long-record streamgages was too sparse for regional analysis and station skew was used for all estimates. Final station flood frequency estimates for all study streamgages are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Regional multiple-regression analysis was used to produce equations for estimating flood frequency statistics from explanatory basin characteristics. Basin characteristics, including physical and climatic variables, were updated for all study streamgages using a geographical information system and geospatial source data. Screening for similar-sized nested basins eliminated hydrologically redundant sites, and screening for eligibility for analysis of explanatory variables eliminated regulated peaks, outburst peaks, and sites with indeterminate basin characteristics. An ordinary least‑squares regression used flood-frequency statistics and basin characteristics for 341 streamgages (284 in Alaska and 57 in Canada) to determine the most suitable combination of basin characteristics for a flood-frequency regression model and to explore regional grouping of streamgages for explaining variability in flood-frequency statistics across the study area. The most suitable model for explaining flood frequency used drainage area and mean annual precipitation as explanatory variables for the entire study area as a region. Final regression equations for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probability discharge in Alaska and conterminous basins in Canada were developed using a generalized least-squares regression. The average standard error of prediction for the regression equations for the various annual exceedance probabilities ranged from 69 to 82 percent, and the pseudo-coefficient of determination (pseudo-R 2 ) ranged from 85 to 91 percent. The regional regression equations from this study were incorporated into the U.S. Geological Survey StreamStats program for a limited area of the State—the Cook Inlet Basin. StreamStats is a national web-based geographic information system application that facilitates retrieval of streamflow statistics and associated information. StreamStats retrieves published data for gaged sites and, for user-selected ungaged sites, delineates drainage areas from topographic and hydrographic data, computes basin characteristics, and computes flood frequency estimates using the regional regression equations.

Alaska↗

USGS2007

Evaluation of Approaches for Managing Nitrate Loading from On-Site Wastewater Systems near La Pine, Oregon

This report presents the results of a study by the U.S. Geological Survey, done in cooperation with the Oregon Department of Environmental Quality and Deschutes County, to develop a better understanding of the effects of nitrogen from on-site wastewater disposal systems on the quality of ground water near La Pine in southern Deschutes County and northern Klamath County, Oregon. Simulation models were used to test the conceptual understanding of the system and were coupled with optimization methods to develop the Nitrate Loading Management Model, a decision-support tool that can be used to efficiently evaluate alternative approaches for managing nitrate loading from on-site wastewater systems. The conceptual model of the system is based on geologic, hydrologic, and geochemical data collected for this study, as well as previous hydrogeologic and water quality studies and field testing of on-site wastewater systems in the area by other agencies. On-site wastewater systems are the only significant source of anthropogenic nitrogen to shallow ground water in the study area. Between 1960 and 2005 estimated nitrate loading from on-site wastewater systems increased from 3,900 to 91,000 pounds of nitrogen per year. When all remaining lots are developed (in 2019 at current building rates), nitrate loading is projected to reach nearly 150,000 pounds of nitrogen per year. Low recharge rates (2-3 inches per year) and ground-water flow velocities generally have limited the extent of nitrate occurrence to discrete plumes within 20-30 feet of the water table; however, hydraulic-gradient and age data indicate that, given sufficient time and additional loading, nitrate will migrate to depths where many domestic wells currently obtain water. In 2000, nitrate concentrations greater than 4 milligrams nitrogen per liter (mg N/L) were detected in 10 percent of domestic wells sampled by Oregon Department of Environmental Quality. Numerical simulation models were constructed at transect (2.4 square miles) and study-area (247 square miles) scales to test the conceptual model and evaluate processes controlling nitrate concentrations in ground water and potential ground-water discharge of nitrate to streams. Simulation of water-quality conditions for a projected future build-out (base) scenario in which all existing lots are developed using conventional on-site wastewater systems indicates that, at equilibrium, average nitrate concentrations near the water table will exceed 10 mg N/L over areas totaling 9,400 acres. Other scenarios were simulated where future nitrate loading was reduced using advanced treatment on-site systems and a development transfer program. Seven other scenarios were simulated with total nitrate loading reductions ranging from 15 to 94 percent; simulated reductions in the area where average nitrate concentrations near the water table exceed 10 mg N/L range from 22 to 99 percent at equilibrium. Simulations also show that the ground-water system responds slowly to changes in nitrate loading due to low recharge rates and ground-water flow velocity. Consequently, reductions in nitrate loading will not immediately reduce average nitrate concentrations and the average concentration in the aquifer will continue to increase for 25-50 years depending on the level and timing of loading reduction. The capacity of the ground-water system to receive on-site wastewater system effluent, which is related to the density of homes, presence of upgradient residential development, ground-water recharge rate, ground-water flow velocity, and thickness of the oxic part of the aquifer, varies within the study area. Optimization capability was added to the study-area simulation model and the combined simulation-optimization model was used to evaluate alternative approaches to management of nitrate loading from on-site wastewater systems to the shallow alluvial aquifer. The Nitrate Loading Management Model (NLMM) was formulated to find the minimum red

Scientific Investigations Report↗

USGS2021

Evaluation of hydrologic simulation models for fields with subsurface drainage to mitigated wetlands in Barnes, Dickey, and Sargent Counties, North Dakota

Proper identification of wetlands, along with a better understanding of the hydrology of mitigated wetlands, is needed to assist with conservation efforts aimed at maintaining the productivity and ecological function (wetland mitigation) of agricultural lands. The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service, completed a study to evaluate two models for simulating hydrologic conditions in fields with subsurface drainage to mitigated wetlands at several sites in North Dakota. These two models were evaluated as possible tools for water resource managers to use for designing wetland mitigation projects in the area in the future. The Soil-Plant-Atmosphere-Water (SPAW) model simulates the daily hydrologic water budgets of agricultural landscapes by two linked routines, one for farm fields (field hydrology) and one for impoundments such as wetlands and ponds (pond model). The DRAINMOD model was used in conjunction with the SPAW model because although the SPAW model can be used to simulate the hydrology of small drainage basins containing wetlands, the SPAW model does not contain routines to simulate drainage, either subsurface drainage or surface (drainage ditches), that can directly affect the wetland hydrology. The wetlands in the study areas in this report are all downstream from and adjacent to drained agricultural fields. SPAW and DRAINMOD models were developed and calibrated at three study areas (study areas B, D, and S) to evaluate how the models simulated field-scale hydrologic characteristics and the water balance in wetlands from January 1, 2003, through December 31, 2018. The SPAW model developed for study area B included five modeled fields in the field hydrology portion of SPAW that contributed inflow to one wetland simulated in the pond model portion of SPAW. Simulated wetland water depths were most similar to water depths measured at site BWET1, with an absolute mean error of 0.10 foot and a root mean square error of 0.14 foot. Site BWET2 had slightly larger errors, with an absolute mean error of 0.22 foot and a root mean square error of 0.28 foot. Simulated water depths were similar to the pattern of measured water depths at BWET1 and BWET2 from about mid-April 2018 through about mid-September 2018, but overpredicted water depths in the fall from about mid-September 2018 through about mid-October 2018. The SPAW model developed for study area D included six modeled fields in the field hydrology portion of SPAW that contributed inflow to five wetlands connected in series in the pond model portion of SPAW. Simulated water depths compared relatively well to water depths in the five wetlands, with the absolute mean error ranging from 0.17 foot (DWET1) to 0.39 foot (DWET2), and the root mean square error ranging from 0.28 foot (DWET1) to 0.56 foot (DWET5). The SPAW model developed for study area S included one modeled field in the field hydrology portion of SPAW that contributed inflow to one wetland in the pond model portion of SPAW. Among the SPAW models developed for the three study areas, the model for study area S had the best comparison between simulated and measured water depths, with an absolute mean error of 0.06 foot and a root mean square error of 0.10 foot. DRAINMOD models were developed and calibrated at the three study areas and provided inflow from subsurface drainage discharge to the SPAW models for simulating water levels in wetlands in the study areas. The calibrated DRAINMOD model for study area B showed the variability of hydrologic processes in the modeled field throughout the wide range of hydrologic conditions from January 1, 2003, through December 31, 2018. In general, the discharge through the modeled subsurface drainage system was in the spring and early summer (April through June) most years, with little to no discharge later in the year. Although the subsurface drainage system in study area D was the most complex among the three study areas and was simplified into a uniform system within DRAINMOD, simulated water table depths at study area D correlated better to measured water table depths compared to results from the model applications at the other two study areas. Simulated water table depths had an absolute mean error of 0.30 foot and root mean square error of 0.37 foot at site DGW1 and an absolute mean error of 0.29 foot and a root mean square error of 0.34 foot at site DGW2. Although the subsurface drainage system in study area S was the simplest and the modeled field was the smallest among the three study areas, simulated water table depths at study area S did not correlate as well to measured water table depths compared to results from the model applications at the other two study areas. The SPAW and DRAINMOD model applications at the three study areas in southeast North Dakota adequately simulated the hydrologic processes for fields with subsurface drainage that are connected to adjacent wetlands. However, more measured data would be needed to fully evaluate the models throughout the range of possible climatic conditions.

North Dakota↗

USGS2008

Estimated Loads of Suspended Sediment and Selected Trace Elements Transported through Milltown Reservoir in the Upper Clark Fork Basin, Montana, Water Years 2004-07

The purpose of this report is to present estimated daily and annual loads of suspended sediment and selected trace elements for water years 2004-07 at two sites upstream and one site downstream from Milltown Reservoir. Milltown Reservoir is a National Priorities List Superfund site in the upper Clark Fork basin of western Montana where sediments enriched in trace elements from historical mining and ore processing have been deposited since the construction of Milltown Dam in 1907. The estimated loads were used to quantify annual net gains and losses (mass balance) of suspended sediment and trace elements within Milltown Reservoir before and after June 1, 2006, which was the start of Stage 1 of a permanent drawdown of the reservoir in preparation for removal of Milltown Dam. This study was done in cooperation with the U.S. Environmental Protection Agency. Daily loads of suspended sediment were estimated for water years 2004-07 by using either high-frequency sampling as part of daily sediment monitoring or regression equations relating suspended-sediment discharge to streamflow. Daily loads of unfiltered-recoverable arsenic, cadmium, copper, iron, lead, manganese, and zinc were estimated by using regression equations relating trace-element discharge to suspended-sediment discharge. Regression equations were developed from data for eriodic water-quality samples collected during water years 2004-07. The equations were applied to daily records of either streamflow or suspended-sediment discharge to produce estimated daily loads. Variations in daily suspended-sediment and trace-element loads generally coincided with variations in streamflow. For most of the period before June 1, 2006, differences in daily loads transported to and from Milltown Reservoir were minor or indicated small amounts of deposition; however, losses of suspended sediment and trace elements from the reservoir occurred during temporary drawdowns in July-August 2004 and October-December 2005. After the start of Stage 1 of the permanent drawdown on June 1, 2006, losses of suspended sediment and trace elements from the reservoir persisted for all streamflow conditions during the entire interval of the Stage 1 drawdown (June 1, 2006-September 30, 2007) within the study period. Estimated daily loads of suspended sediment and trace elements were summed for each year to produce estimated annual loads used to determine the annual net gains (deposition) or losses (erosion) of each constituent within Milltown Reservoir during water years 2004-07. During water year 2004, there was an annual net gain of suspended sediment in the reservoir. The annual net gains and losses of trace elements were inconsistent in water year 2004, with gains occurring for arsenic ad iron, but losses occurring for cadmium, copper, lead, manganese, and zinc. In water year 2005, there were annual net gains of suspended sediment and all the trace elements within the reservoir. In water year 2006, there were annual net losses of all constituents from the reservoir, likely as the result of sediment erosion from the reservoir during both a temporary drawdown in October-December 2005 and Stage 1 of the permanent drawdown that continued after June 1, 2006. In water year 2007, when the Stage 1 drawdown was in effect for the entire year, there were large annual net losses of suspended sediment and trace elements from the reservoir. The annual net losses of constituents from Milltown Reservoir in water year 2007 were the largest of any year during the 2004-07 study period. In water year 2007, the annual net loss of suspended sediment from the reservoir was 130,000 tons, which was more than double (about 222 percent) the combined inflow to the reservoir. The largest annual net losses of trace elements in water year 2007, in percent of the combined inflow to the reservoir, occurred for cadmium, copper, lead, and zinc-about 190 percent for cadmium, 170 percent for copper, 150 percent for lead, and 238 p

Scientific Investigations Report↗

USGS2009

Estimated Loads of Suspended Sediment and Selected Trace Elements Transported through the Milltown Reservoir Project Area Before and After the Breaching of Milltown Dam in the Upper Clark Fork Basin, Montana, Water Year 2008

This report presents estimated daily and cumulative loads of suspended sediment and selected trace elements transported during water year 2008 at three streamflow-gaging stations that bracket the Milltown Reservoir project area in the upper Clark Fork basin of western Montana. Milltown Reservoir is a National Priorities List Superfund site where sediments enriched in trace elements from historical mining and ore processing have been deposited since the construction of Milltown Dam in 1907. Milltown Dam was breached on March 28, 2008, as part of Superfund remedial activities to remove the dam and contaminated sediment that had accumulated in Milltown Reservoir. The estimated loads transported through the project area during the periods before and after the breaching of Milltown Dam, and for the entire water year 2008, were used to quantify the net gain or loss (mass balance) of suspended sediment and trace elements within the project area during the transition from a reservoir environment to a free-flowing river. This study was done in cooperation with the U.S. Environmental Protection Agency. Streamflow during water year 2008 compared to long-term streamflow, as represented by the record for Clark Fork above Missoula (water years 1930-2008), generally was below normal (long-term median) from about October 2007 through April 2008. Sustained runoff started in mid-April, which increased flows to near normal by mid-May. After mid-May, flows sharply increased to above normal, reaching a maximum daily mean streamflow of 16,800 cubic feet per second (ft3/s) on May 21, which essentially equaled the long-term 10th-exceedance percentile for that date. Flows substantially above normal were sustained through June, then decreased through the summer and reached near-normal by August. Annual mean streamflow during water year 2008 (3,040 ft3/s) was 105 percent of the long-term mean annual streamflow (2,900 ft3/s). The annual peak flow (17,500 ft3/s) occurred on May 21 and was 112 percent of the long-term mean annual peak flow (15,600 ft3/s). About 81 percent of the annual flow volume was discharged during the post-breach period. Daily loads of suspended sediment were estimated directly by using high-frequency sampling of the daily sediment monitoring. Daily loads of unfiltered-recoverable arsenic, cadmium, copper, iron, lead, manganese, and zinc were estimated by using regression equations relating trace-element discharge to either streamflow or suspended-sediment discharge. Regression equations for estimating trace-element discharge in water year 2008 were developed from instantaneous streamflow and concentration data for periodic water-quality samples collected during all or part of water years 2004-08. The equations were applied to records of daily mean streamflow or daily suspended-sediment loads to produce estimated daily trace-element loads. Variations in daily suspended-sediment and trace-element loads generally coincided with variations in streamflow. Relatively small to moderately large daily net losses from the project area were common during the pre-breach period when low-flow conditions were prevalent. Outflow loads from the project area sharply increased immediately after the breaching of Milltown Dam and during the rising limb and peak flow of the annual hydrograph. Net losses of suspended sediment and trace elements from the project area decreased as streamflow decreased during the summer, eventually becoming small or reaching an approximate net balance between inflow and outflow. Estimated daily loads of suspended sediment and trace elements for all three stations were summed to determine cumulative inflow and outflow loads for the pre-breach and post-breach periods, as well as for the entire water year 2008. Overall, the mass balance between the combined inflow loads from two upstream source areas (upper Clark Fork and Blackfoot River basins) and the outflow loads at Clark Fork above Missoula indicates net losses

Scientific Investigations Report↗

USGS2013

Groundwater quality and the relation between pH values and occurrence of trace elements and radionuclides in water samples collected from private wells in part of the Kickapoo Tribe of Oklahoma Jurisdictional Area, central Oklahoma, 2011

From 1999 to 2007, the Indian Health Service reported that gross alpha-particle activities and concentrations of uranium exceeded the Maximum Contaminant Levels for public drinking-water supplies in water samples from six private wells and two test wells in a rural residential neighborhood in the Kickapoo Tribe of Oklahoma Jurisdictional Area, in central Oklahoma. Residents in this rural area use groundwater from Quaternary-aged terrace deposits and the Permian-aged Garber-Wellington aquifer for domestic purposes. Uranium and other trace elements, specifically arsenic, chromium, and selenium, occur naturally in rocks composing the Garber-Wellington aquifer and in low concentrations in groundwater throughout its extent. Previous studies have shown that pH values above 8.0 from cation-exchange processes in the aquifer cause selected metals such as arsenic, chromium, selenium, and uranium to desorb (if present) from mineral surfaces and become mobile in water. On the basis of this information, the U.S. Geological Survey, in cooperation with the Kickapoo Tribe of Oklahoma, conducted a study in 2011 to describe the occurrence of selected trace elements and radionuclides in groundwater and to determine if pH could be used as a surrogate for laboratory analysis to quickly and inexpensively identify wells that might contain high concentrations of uranium and other trace elements. The pH and specific conductance of groundwater from 59 private wells were measured in the field in an area of about 18 square miles in Lincoln and Pottawatomie Counties. Twenty of the 59 wells also were sampled for dissolved concentrations of major ions, trace elements, gross alpha-particle and gross beta-particle activities, uranium, radium-226, radium-228, and radon-222 gas. Arsenic concentrations exceeded the Maximum Contaminant Level of 10 micrograms per liter in one sample having a concentration of 24.7 micrograms per liter. Selenium concentrations exceeded the Maximum Contaminant Level of 50 micrograms per liter in one sample having a concentration of 147 micrograms per liter. Both samples had alkaline pH values, 8.0 and 8.4, respectively. Uranium concentrations ranged from 0.02 to 383 micrograms per liter with 5 of 20 samples exceeding the Maximum Contaminant Level of 30 micrograms per liter; the five wells with uranium concentrations exceeding 30 micrograms per liter had pH values ranging from 8.0 to 8.5. Concentrations of uranium and radon-222 and gross alpha-particle activity showed a positive relation to pH, with the highest concentrations and activity in samples having pH values of 8.0 or above. The groundwater samples contained dissolved oxygen and high concentrations of bicarbonate; these characteristics are also factors in increasing uranium solubility. Concentrations of radium-226 and radium-228 (combined) ranged from 0.03 to 1.7 picocuries per liter, with a median concentration of 0.45 picocuries per liter for all samples. Radon-222 concentrations ranged from 95 to 3,600 picocuries per liter with a median concentration of 261 picocuries per liter. Eight samples having pH values ranging from 8.0 to 8.7 exceeded the proposed Maximum Contaminant Level of 300 picocuries per liter for radon-222. Eight samples exceeded the 15 picocuries per liter Maximum Contaminant Level for gross alpha-particle activity at 72 hours (after sample collection) and at 30 days (after the initial count); those samples had pH values ranging from 8.0 to 8.5. Gross beta-particle activity increased in 15 of 21 samples during the interval from 72 hours to 30 days. The increase in gross beta-particle activity over time probably was caused by the ingrowth and decay of uranium daughter products that emit beta particles. Water-quality data collected for this study indicate that pH values above 8.0 are associated with potentially high concentrations of uranium and radon-222 and high gross alpha-particle activity in the study area. High pH values also are associated with potentially high concentrations of arsenic, chromium, and selenium in groundwater when these elements occur in the aquifer matrix along groundwater-flow paths.

Oklahoma↗

USGS2021

Groundwater flow conceptualization of the Pahute Mesa–Oasis Valley Groundwater Basin, Nevada—A synthesis of geologic, hydrologic, hydraulic-property, and tritium data

This report provides a groundwater-flow conceptualization that integrates geologic, hydrologic, hydraulic-property, and radionuclide data in the Pahute Mesa–Oasis Valley (PMOV) groundwater basin, southern Nevada. Groundwater flow in the PMOV basin is of interest because 82 underground nuclear tests were detonated, most near or below the water table. A potentiometric map and nine sets of hydrostratigraphic and hydrologic cross sections supplement the conceptualization. Potentiometric contours indicate that groundwater in the PMOV basin generally flows south-southwest and discharges at Oasis Valley. Groundwater encounters an alternating sequence of low- and high-transmissivity rocks, referred to as dams and pools, respectively, as it moves from east to west across eastern Pahute Mesa. Flow from all Pahute Mesa nuclear tests is to Oasis Valley and is well-constrained by water-level data. Flow converges along a corridor of high transmissivity between Pahute Mesa and Oasis Valley. The location of the lateral PMOV basin boundary is well defined, and this boundary, with a few minor exceptions, represents a no-flow boundary. Some boundary uncertainty exists in the northeastern part of the basin, but potential flow-rate estimates across the northeastern boundary resulting from this uncertainty are small relative to the basin groundwater budget. Recharge in the PMOV basin is derived from episodic pulses of modern water and the diffuse percolation of old water (greater than 1,000 years). Episodic recharge is a minor recharge component observed as a rise in groundwater levels that occurs 3 months to 1 year following a wet winter. Minor amounts of episodic recharge through an unsaturated zone in excess of 1,000 feet (ft) requires preferential flow through faults and fractures. The dominant recharge component is slow, steady, diffuse percolation of old water through the unsaturated zone. A large component of old water recharging the groundwater system is consistent with observations of isotopically light deuterium and oxygen 18 compositions in water from wells on Pahute Mesa and central Oasis Valley. About half the recharge in the PMOV basin is derived from the eastern Pahute Mesa area. The remaining recharge is derived primarily from other highland areas including Timber Mountain, Belted and Kawich Ranges, and Black Mountain. The PMOV groundwater system is nearly steady state, where recharge is balanced by the 5,900 acre-feet per year of natural discharge at Oasis Valley. This assumption is reasonable because the basin is dominated by steady-state conditions, where long-term changes in groundwater storage are minimal. Total groundwater withdrawals from 1963 to 2018 have amounted to less than 10 percent of annual groundwater discharge and less than 0.2 percent of the basin’s groundwater storage. Therefore, present-day (2020) conditions are considered representative of predevelopment (pre-1950) conditions in nearly all areas of the basin. The lower PMOV basin boundary is defined at 4,000 ft below the water table to encompass all underground nuclear tests and tritium plumes. This boundary defines the lower boundary of radionuclide migration. However, nearly all flow and tritium transport occur in the upper 1,600 ft of the saturated zone because a transmissivity-with-depth relation indicates that greater than 90 percent of the transmissivity contributing to groundwater flow occurs within 1,600 ft of the water table. Rocks at deeper depths have low transmissivity because argillic and mineralized alterations plug the fractures. Volcanic rocks form the primary aquifers and confining units in the PMOV basin. Volcanic hydrogeologic units (HGUs) and hydrostratigraphic units (HSUs) have transmissivity distributions that span up to eight orders of magnitude with considerable overlap between distributions. Despite the large overlap between units, mean transmissivities of aquifers are one-to-two orders of magnitude greater than the confining units. However, all volcanic-rock HGUs and HSUs are composite units, meaning that they can function spatially as either an aquifer or confining unit .

Nevada↗

USGS2005

Hydrogeology and simulation of source areas of water to production wells in a colluvium-mantled carbonate-bedrock aquifer near Shippensburg, Cumberland and Franklin Counties, Pennsylvania

This report presents the results of a study by the U.S. Geological Survey in cooperation with the Shippensburg Borough Authority to evaluate the source areas of water to production wells in a colluvium-mantled carbonate-bedrock aquifer in Cumberland and Franklin Counties, Pa. The areal extent of the zone of contribution was simulated for three production wells near Shippensburg, Pa. by use of a ground-water-flow model. A 111-square-mile area was selected as the model area and includes areas of the South Mountain Section and the Great Valley Section of the Valley and Ridge Physiographic Province. Within the model area, the geologic units in the South Mountain area are predominantly metamorphic rocks and the geologic units in the Great Valley are predominantly carbonate rocks. Hydrologic and geologic information were compiled to establish a conceptual model of ground-water flow. Characteristics of aquifer materials were determined, and streamflow and water levels were measured. Streamflow measurements in November 2003 showed all streams lost water as they flowed from South Mountain over the colluvium-mantled carbonate aquifer into the Great Valley. Some streams lost more than 1 cubic foot per second to the aquifer in this area. The Shippensburg Borough Authority owns three production wells in the model area. Two wells, Cu 969 and Fr 823, are currently (2004) used as production wells and produce 500,000 and 800,000 gallons per day, respectively. Well Cu 970 is intended to be brought on line as a production well in the future. Water levels were measured in 43 wells to use for model calibration. Water-level fluctuations and geophysical logs indicated confined conditions in well Cu 970. Ground-water flow was simulated with a model that consisted of two vertical layers, with five zones in each layer. The units were hydrostratigraphic units that initially were based on geologic formations, but boundaries were adjusted during model calibration. Model calibration resulted in a root mean square error of 9.8 feet. A parameter-estimation package was used during model calibration to estimate three parameters. The parameter estimation resulted in a value of 233 feet per day for horizontal hydraulic conductivity of the highly fractured carbonate rocks and sandy colluvium in layer 1; 3.97 feet per day for horizontal hydraulic conductivity of the ridge-forming unit in layer 1; and a value of 1.73 for horizontal anisotropy in both layers. The calibrated model was used to delineate the areal extent of the zone of contribution for wells Cu 969 and Fr 823. Although well Cu 970 is not currently (2004) being used, the areal extent of its zone of contribution also was simulated without additional model calibration. The shape of the areal extent of the zone of contribution was similar for each well and included an area that extended from the well southwest along the Tomstown Formation, and then extended southeast into the metamorphic rocks of South Mountain. The contributing areas from the watersheds of losing streams were also delineated because losing stream reaches bisect the areal extent of the zones of contribution. Spatial uncertainty of the areal extent of the zone of contribution was illustrated using a Monte-Carlo analysis. The model was run 1,000 times using randomly generated parameter sets that were normally distributed within the confidence interval around the optimal values for the three estimated parameters. The model converged and had a reasonable water budget for 980 of the model runs. For each of those 980 model runs, the recharge area was determined, and the results for all runs were compiled and contoured. The results of the Monte-Carlo analysis were compared to the results of the deterministic model, illustrating that the deterministic model has the greatest certainty in the area closest to each well in the Tomstown Formation. The areas farther from the well, upgradient, and in the metamorphic rocks have a higher degree

Pennsylvania↗