Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature, Scientific Data”

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

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

At least 1,009 records · Page 56Linked to original sources

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↗

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in northeastern Missouri, 2021

Bathymetric data were collected at 12 water-supply lakes in northeastern Missouri by the U.S. Geological Survey (USGS) in cooperation with the Missouri Department of Natural Resources (MoDNR) and various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. The lakes were surveyed in March through May 2021. Ten of the lakes had been surveyed previously by the USGS, and the recent surveys were compared to the earlier surveys to document the changes in the bathymetric surface and capacity of the lakes. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data at five of the lakes were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. Data points from the various sources were exported at a gridded data resolution appropriate to each lake, either 0.82 foot, 1.64 feet, or 3.28 feet. Data outside the multibeam survey extent and greater than the surveyed water-surface elevation were obtained from data collected using aerial light detection and ranging (lidar) point cloud data. A linear enforcement technique was used to add points to the dataset in areas of sparse data (the upper ends of coves where the water was shallow or aquatic vegetation precluded data acquisition) based on surrounding multibeam and upland data values. The various point datasets were used to produce a three-dimensional triangulated irregular network surface of the lake-bottom elevations for each lake. A surface area and capacity table was produced from the three-dimensional surface for each lake showing surface area and capacity at specified lake water-surface elevations. Various quality-assurance tests were conducted to ensure quality data were collected with the multibeam, including beam angle checks and patch tests. Additional quality-assurance tests were conducted on the gridded bathymetric data from the survey, the bathymetric surface created from the gridded data, and the contours created from the bathymetric survey. If there were data from a previous bathymetric survey for a given lake, a bathymetric change map was generated from the elevation difference between the previous survey and the 2021 bathymetric survey data points. After reconciling any vertical datum disagreement between the previous survey data and the 2021 survey datum, coincident points between the surveys were identified, and a bathymetric change map was generated using the coincident point data. The mean elevation change between all repeat surveys at most lakes was positive, indicating sedimentation. Relative to previous surveys, the change in capacity at the primary spillway elevation ranged from a 7.7-percent decrease at Memphis Reservoir to a 3.9-percent increase at Old Lake (Bowling Green West). The mean bathymetric change ranged from 0.03 foot at Hazel Creek and 0.07 foot at Shelbina Lake and Bowling Green Reservoir (Jack Floyd Memorial Lake) to 0.63 at Memphis Lake (Lake Showme) and 0.88 at Memphis Reservoir. The time-averaged mean bathymetric change ranged from 0.002 foot per year at Hazel Creek Lake to 0.044 foot per year at Memphis Reservoir. The computed volumetric sedimentation rate generally ranged from 0.14 to 6.80 acre-feet per year at Shelbina Lake and Memphis Lake (Lake Showme), respectively; however, Forest Lake had a substantially larger sedimentation rate of 17.0 acre-feet per year. Some changes observed in some bathymetric change maps are believed to result from the difference in data collection equipment and techniques between the previous and present bathymetric surveys, whereas other erosional features around the perimeter of certain lakes may be the result of wave action during low-water years.

Missouri↗

Geologic map and borehole stratigraphy of Hinkley Valley and vicinity, San Bernardino County, California

Hinkley Valley, in the central to western Mojave Desert of southeastern California, has a long historical record owing to its position as a crossroads for rail and road traffic and its position adjacent to the Mojave River. Subflow in the Mojave River provided groundwater recharge that maintained water consumption and demand by way of shallow wells for local agriculture in the valley. Its crossroads position led to construction of several power-transmission lines, pipeline, and communications cable routes that transect Hinkley Valley. One of these, a natural gas pipeline and its associated compressor station, was the locus of hexavalent chromium, Cr(VI), released into, and consequent contamination of, groundwater. Understanding the movement and fate of the contaminants is a complex hydrologic and geochemical problem. Geologic mapping of the Hinkley Valley area provides framework elements for use in resolving this problem. This report provides new information on surface and subsurface geology to better constrain the origin and geometry of hydrologically important deposits in the Hinkley Valley area and describes youthful faults that may control sediment distribution and groundwater flow. The geologic map (sheet 1) presents substantial new information on surficial geology, including Pliocene deposits, but does not contain significant new work on bedrock. Bedrock investigations were specific to identifying youthful faults and representative outcrops for rocks that were penetrated by boreholes in the valley. Special attention was placed on locating and describing youthful faults. In addition, we analyzed gravity data to (1) map horizontal gradients that we interpret to reflect long-term fault traces and to (2) estimate the depth to bedrock, which is defined as Miocene and older intrusive and metamorphic rocks for the purposes of this report. The subsurface geology of Hinkley Valley was investigated by examining borehole sediment cores and rock encountered at the base of the sediment section. We analyzed the core to determine depositional environments, provenance, and age of the sediment that infilled the valley. Valleys, mountains, and basins in the Hinkley Valley area are topographically complex and incompletely named. The nearly flat floored Hinkley Valley slopes gently northward. It is framed by Mount General and the informally named “Hinkley hills” (southeast of Mount General) on the northeast and by Iron Mountain and Lynx Cat Mountain on the southwest, although breaks in the western mountains allow stream connections between Hinkley Valley and another valley to the west that is herein referred to as Hawes valley. At its south end, Hinkley Valley is traversed by the entrenched Mojave River, which passes east out of the valley past Barstow. North of Hinkley Valley, a few low hills (including Red Hill) separate the valley from a broad west-sloping piedmont that is part of the physiographic Harper Basin (of which the Harper Lake playa is the center). The lower part of this piedmont, however, is referred to as Water Valley, although it is not a distinct valley. The name derives from groundwater sourced from subflow in the Mojave River, which caused shallow water and even artesian flow in Water Valley but not in other parts of the Harper Basin. When water filled the Harper Basin to form Pleistocene Lake Harper it not only submerged Water Valley but also northern Hinkley Valley.

California↗

The Connecticut Streamflow and Sustainable Water Use Estimator—A decision-support tool to estimate water availability at ungaged stream locations in Connecticut

Freshwater streams in Connecticut are subject to many competing demands, including public water supply; agricultural, commercial, and industrial water use; and ecosystem and habitat needs. In recent years, drought has further stressed Connecticut’s water resources. To sustainably allocate and manage water resources among these competing uses, Federal, State, and local water-resource managers require data and modeling tools to estimate the water availability at a variety of temporal and spatial scales for planning purposes. The Connecticut Streamflow and Sustainable Water Use Estimator (CT SSWUE), developed by the U.S. Geological Survey in cooperation with the Connecticut Department of Energy and Environmental Protection, is a decision-support tool for estimating daily unaltered streamflow and sustainable water use at ungaged sites in Connecticut. The CT SSWUE estimates unaltered daily mean streamflow and water-use-adjusted streamflow for the period from October 1, 1960, to September 30, 2015, and the monthly sustainable net withdrawal at ungaged sites in Connecticut. Unaltered streamflow is the estimated daily mean streamflow in a drainage basin in the absence of any water withdrawals or wastewater discharges and with minimal human development. Sustainable net withdrawal is the maximum net withdrawal (withdrawal minus wastewater discharges) that can be drawn from a basin without critically depleting the water available through natural streamflow patterns. Sustainable net withdrawal is defined for this study as the difference between the unaltered daily mean streamflow and a user-defined target minimum streamflow. Weighted least squares and Tobit regression techniques were used to develop equations for estimating streamflow at ungaged sites at 19 streamflow quantiles with exceedance probabilities ranging from 0.005 to 99.995 percent. Regressions were based on streamflow quantiles and basin characteristics from 36 reference streamgages in and around Connecticut. Four basin characteristics—drainage area, mean of the soil permeability, mean of the average annual precipitation, and ratio of the length of streams that overlay sand and gravel deposits to the total length of streams in the basin—are used as explanatory variables in the equations. At an ungaged site, interpolation between the streamflow quantiles estimated from the regression equations produces a continuous flow-duration curve. A time series of daily mean streamflow at an ungaged site is then estimated by assuming that for each day, the streamflow quantile occurs on the same date at both a reference streamgage and the ungaged site. In a remove-one cross validation, estimated unaltered daily mean streamflow agreed well with observed values at reference streamgages, with a few exceptions. Nash Sutcliffe efficiency ranged from −0.43 to 0.97 with a median value of 0.88. The normalized root-mean-square error ranged from 16.6 to 120.4 percent with a median value of 34.5 percent. An empirical method for estimating 95-percent prediction intervals for unaltered daily and monthly mean streamflow was developed and tested by using the cross-validation data. Prediction intervals for unaltered daily mean streamflow at the cross-validation reference streamgages performed well in most cases. Gaged streamflow values from the cross-validation data fell within the prediction intervals a median 96.6 percent of the time for daily mean time series and 93.9 percent of the time for monthly mean time series. The CT SSWUE computes water-use-adjusted streamflow using spatially referenced water-use information provided by the Connecticut Department of Energy and Environmental Protection. Available water-use information included permitted and registered water withdrawals and permitted wastewater discharges during 1998 to 2015 for the Thames River Basin and central coastal drainage basins. Water-use information was incorporated into the U.S. Geological Survey StreamStats web application for Connecticut and can be used for computing water-use-adjusted streamflow and sustainable net withdrawal at selected points of interest. Altered daily streamflow is computed by applying average daily withdrawals and wastewater discharges to the water balance equation. Average daily surface water withdrawals and wastewater discharges are applied directly to the daily water balance equation. Time-lagged alterations on streamflow from groundwater withdrawals or wastewater discharges are estimated by using a response-coefficient method developed from results of previously published, calibrated groundwater models.

Connecticut↗

Arsenic, Boron, and Fluoride Concentrations in Ground Water in and Near Diabase Intrusions, Newark Basin, Southeastern Pennsylvania

During an investigation in 2000 by the U.S. Environmental Protection Agency (USEPA) of possible contaminant releases from an industrial facility on Congo Road near Gilbertsville in Berks and Montgomery Counties, southeastern Pennsylvania, concentrations of arsenic and fluoride above USEPA drinking-water standards of 10 ?g/L and 4 mg/L, respectively, and of boron above the USEPA health advisory level of 600 ?g/L were measured in ground water in an area along the northwestern edge of the Newark Basin. In 2003, the USEPA requested technical assistance from the U.S. Geological Survey (USGS) to help identify sources of arsenic, boron, and fluoride in the ground water in the Congo Road area, which included possible anthropogenic releases and naturally occurring mineralization in the local bedrock aquifer, and to identify other areas in the Newark Basin of southeastern Pennsylvania with similarly elevated concentrations of these constituents. The USGS reviewed available data and collected additional ground-water samples in the Congo Road area and four similar hydrogeologic settings. The Newark Basin is the largest of the 13 major exposed Mesozoic rift basins that stretch from Nova Scotia to South Carolina. Rocks in the Newark Basin include Triassic through Jurassic-age sedimentary sequences of sandstones and shales that were intruded by diabase. Mineral deposits of hydrothermal origin are associated with alteration zones bordering intrusions of diabase and also occur as strata-bound replacement deposits of copper and zinc in sedimentary rocks. The USGS review of data available in 2003 showed that water from about 10 percent of wells throughout the Newark Basin of southeastern Pennsylvania had concentrations of arsenic greater than the USEPA maximum contaminant level (MCL) of 10 ?g/L; the highest reported arsenic concentration was at about 70 ?g/L. Few data on boron were available, and the highest reported boron concentration in well-water samples was 60 ?g/L in contrast to concentrations over 5,000 ?g/L in the Congo Road area. Although concentrations of fluoride up to 4 mg/L were reported for a few well-water samples collected throughout the Newark Basin, about 90 percent of the samples had concentrations of 0.5 mg/L or less. The USGS sampled 58 wells primarily in 5 areas in the Newark Basin, southeastern Pennsylvania, from February 2004 through April 2005 to identify other possible areas of elevated arsenic, boron, and fluoride and to characterize the geochemical environment associated with elevated concentrations of these constituents. Sampled wells included 12 monitor wells at an industrial facility near Congo Road, 45 private-supply wells in Berks, Montgomery, and Bucks Counties, and 1 private-supply well near Dillsburg, York County, an area where elevated fluoride in ground water had been reported in the adjacent Gettysburg Basin. Wells were sampled in transects from the diabase through the adjacent hornfels and into the unaltered shales of the Brunswick Group. Field measurements were made of pH, temperature, dissolved oxygen concentration, and specific conductance. Samples were analyzed in the laboratory for major ions, nutrients, total organic carbon, dissolved and total concentrations of selected trace elements, and boron isotopic composition. Generally, the ground water from the 46 private-supply wells had relatively neutral to alkaline pH (ranging from 6.1 to 9.1) and moderate concentrations of dissolved oxygen. Most water samples were of the calcium-bicarbonate type. Concentrations of arsenic up to 60 ?g/L, boron up to 3,950 ?g/L, and fluoride up to 0.70 mg/L were measured. Drinking-water standards or health advisories (for constituents that do not have standards established) were exceeded most frequently (about 20 percent of samples) for arsenic and boron and less frequently (6 percent or less of samples) for total iron, manganese, sulfate, nitrate, lead, molybdenum, and strontium. In water from 12 monitor

Scientific Investigations Report↗

Surface-Water Quality Conditions and Long-Term Trends at Selected Sites within the Ambient Water-Quality Monitoring Network in Missouri, Water Years 1993-2008

The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, collects data pertaining to the surface-water resources of Missouri. These data are collected as part of the Missouri Ambient Water-Quality Monitoring Network and constitute a valuable source of reliable, impartial, and timely information for developing an improved understanding of water resources in the State. Six sites from the Ambient Water-Quality Monitoring Network, with data available from the 1993 through 2008 water years, were chosen to compare water-quality conditions and long-term trends of dissolved oxygen, selected physical properties, total suspended solids, dissolved nitrate plus nitrite as nitrogen, total phosphorous, fecal indicator bacteria, and selected trace elements. The six sites used in the study were classified in groups corresponding to the physiography, main land use, and drainage basin size, and represent most stream types in Missouri. Long-term trends in this study were analyzed using flow-adjusted and non-flow adjusted models. Highly censored datasets (greater than 5 percent but less than 50 percent censored values) were not flow-adjusted. Trends that were detected can possibly be related to changes in agriculture or urban development within the drainage basins. Trends in nutrients were the most prevalent. Upward flow-adjusted trends in dissolved nitrate plus nitrite (as nitrogen) concentrations were identified at the Elk River site, and in total phosphorus concentrations at the South Fabius and Grand River sites. A downward flow-adjusted trend was identified in total phosphorus concentrations from Wilson Creek, the only urban site in the study. The downward trend in phosphorus possibly was related to a phosphorus reduction system that began operation in 2001 at a wastewater treatment plant upstream from the sampling site. Total suspended solids concentrations indicated an upward non-flow adjusted trend at the two northern sites (South Fabius and Grand Rivers). The increase in total suspended solids concentrations could be because of soil erosion from land cultivated for row crops. Most trace element data examined in the study were highly censored and could not be used for flow-adjusted trend analyses. Water-quality conditions were assessed to explore relations between data from sites and to the State water-quality standards where applicable for selected constituents. Streamflow varied at each site because of drainage area, land use, and groundwater inputs. Dissolved oxygen and water temperature were similar at all sites except the urban site located on Wilson Creek. Specific conductance was similar between the most northern (South Fabius and Grand River sites) and the most southern sites (Current and Elk River sites). Total suspended solids concentrations were near the method reporting level at all sites, except the northern sites. Streams in northern Missouri are more turbid than streams in southern Missouri and are affected by large volumes of sediment deposition because of soil erosion from land cultivated for row crops. Geometric means of Escherichia coli were calculated from the recreational seasons within the study period. Only the Grand River site exceeded the whole-body-contact standard for frequently used waters. The South Fabius and Grand River sites and the Wilson Creek site had statistically larger densities of both fecal indicator bacteria types than the remaining sites.

Scientific Investigations Report↗

Streamflow Simulations and Percolation Estimates Using the Soil and Water Assessment Tool for Selected Basins in North-Central Nebraska, 1940-2005

The U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, used the Soil and Water Assessment Tool to simulate streamflow and estimate percolation in north-central Nebraska to aid development of long-term strategies for management of hydrologically connected ground and surface water. Although groundwater models adequately simulate subsurface hydrologic processes, they often are not designed to simulate the hydrologically complex processes occurring at or near the land surface. The use of watershed models such as the Soil and Water Assessment Tool, which are designed specifically to simulate surface and near-subsurface processes, can provide helpful insight into the effects of surface-water hydrology on the groundwater system. The Soil and Water Assessment Tool was calibrated for five stream basins in the Elkhorn-Loup Groundwater Model study area in north-central Nebraska to obtain spatially variable estimates of percolation. Six watershed models were calibrated to recorded streamflow in each subbasin by modifying the adjustment parameters. The calibrated parameter sets were then used to simulate a validation period; the validation period was half of the total streamflow period of record with a minimum requirement of 10 years. If the statistical and water-balance results for the validation period were similar to those for the calibration period, a model was considered satisfactory. Statistical measures of each watershed model's performance were variable. These objective measures included the Nash-Sutcliffe measure of efficiency, the ratio of the root-mean-square error to the standard deviation of the measured data, and an estimate of bias. The model met performance criteria for the bias statistic, but failed to meet statistical adequacy criteria for the other two performance measures when evaluated at a monthly time step. A primary cause of the poor model validation results was the inability of the model to reproduce the sustained base flow and streamflow response to precipitation that was observed in the Sand Hills region. The watershed models also were evaluated based on how well they conformed to the annual mass balance (precipitation equals the sum of evapotranspiration, streamflow/runoff, and deep percolation). The model was able to adequately simulate annual values of evapotranspiration, runoff, and precipitation in comparison to reported values, which indicates the model may provide reasonable estimates of annual percolation. Mean annual percolation estimated by the model as basin averages varied within the study area from a maximum of 12.9 inches in the Loup River Basin to a minimum of 1.5 inches in the Shell Creek Basin. Percolation also varied within the studied basins; basin headwaters tended to have greater percolation rates than downstream areas. This variance in percolation rates was mainly was because of the predominance of sandy, highly permeable soils in the upstream areas of the modeled basins.

Scientific Investigations Report↗

Hydrology and water quality of Shell Lake, Washburn County, Wisconsin, with special emphasis on the effects of diversion and changes in water level on the water quality of a shallow terminal lake

Shell Lake is a relatively shallow terminal lake (tributaries but no outlets) in northwestern Wisconsin that has experienced approximately 10 feet (ft) of water-level fluctuation over more than 70 years of record and extensive flooding of nearshore areas starting in the early 2000s. The City of Shell Lake (City) received a permit from the Wisconsin Department of Natural Resources in 2002 to divert water from the lake to a nearby river in order to lower water levels and reduce flooding. Previous studies suggested that water-level fluctuations were driven by long-term cycles in precipitation, evaporation, and runoff, although questions about the lake’s connection with the groundwater system remained. The permit required that the City evaluate assumptions about lake/groundwater interactions made in previous studies and evaluate the effects of the water diversion on water levels in Shell Lake and other nearby lakes. Therefore, a cooperative study between the City and U.S. Geological Survey (USGS) was initiated to improve the understanding of the hydrogeology of the area and evaluate potential effects of the diversion on water levels in Shell Lake, the surrounding groundwater system, and nearby lakes. Concerns over deteriorating water quality in the lake, possibly associated with changes in water level, prompted an additional cooperative project between the City and the USGS to evaluate efeffects of changes in nutrient loading associated with changes in water levels on the water quality of Shell Lake. Numerical models were used to evaluate how the hydrology and water quality responded to diversion of water from the lake and historical changes in the watershed. The groundwater-flow model MODFLOW was used to simulate groundwater movement in the area around Shell Lake, including groundwater/surface-water interactions. Simulated results from the MODFLOW model indicate that groundwater flows generally northward in the area around Shell Lake, with flow locally converging toward the lake. Total groundwater inflow to Shell Lake is small (approximately 5 percent of the water budget) compared with water entering the lake from precipitation (83 percent) and surface-water runoff (13 percent). The MODFLOW model also was used to simulate average annual hydrologic conditions from 1949 to 2009, including effects of the removal of 3 billion gallons of water during 2003–5. The maximum decline in simulated average annual water levels for Shell Lake due to the diversion alone was 3.3 ft at the end of the diversion process in 2005. Model simulations also indicate that although water level continued to decline through 2009 in response to local weather patterns (local drought), the effects of the diversion decreased after the diversion ceased; that is, after 4 years of recovery (2006–9), drawdown attributable to the diversion alone decreased by about 0.6 ft because of increased groundwater inflow and decreased lake-water outflow to groundwater caused by the artificially lower lake level. A delayed response in drawdown of less than 0.5 ft was transmitted through the groundwater-flow system to upgradient lakes. This relatively small effect on upgradient lakes is attributed in part to extensive layers of shallow clay that limit lake/groundwater interaction in the area. Data collected in the lake indicated that Shell Lake is polymictic (characterized by frequent deep mixing) and that its productivity is limited by the amount of phosphorus in the lake. The lake was typically classified as oligotrophic-mesotrophic in June, mesotrophic in July, and mesotrophic-eutrophic in August. In polymictic lakes like Shell Lake, phosphorus released from the sediments is not trapped near the bottom of the lake but is intermittently released to the shallow water, resulting in deteriorating water quality as summer progresses. Because the productivity of Shell Lake is limited by phosphorus, the sources of phosphorus to the lake were quantified, and the response in water quality to changes in phosphorus inputs were evaluated by means of eutrophication models. During 2009, the total input of phosphorus to Shell Lake was 1,730 pounds (lb), of which 1,320 lb came from external sources (76 percent) and 414 lb came from internal loading from sediments in the lake (24 percent). The largest external source was from surface-water runoff, which delivered about 52 percent of the total phosphorus load compared with about 13 percent of the water input. The second largest source was from precipitation (wetfall and dryfall), which delivered 19 percent of the load compared to about 83 percent of the water input. Contributions from septic systems and groundwater accounted for about 3 and 2 percent, respectively. Increased runoff raises water levels in the lake but does not necessarily increase phosphorus loading because phosphorus concentrations in the tributaries decline during increased flow, possibly because of shorter retention times in upstream wetlands. Phosphorus loading to the lake in 2009 represented what occurred after a series of dry years; therefore, this information was combined with data from 2011, a wet year, to estimate phosphorus loading during a range of hydrologic conditions by estimating loading from each component of the phosphorus budget for each year from 1949 to 2011. Comparisons of historical water-quality records with historical water levels and applications of a hydrodynamic model (Dynamic Lake Model, DLM) and empirical eutrophication models were used to understand how changes in water level and the coinciding changes in phosphorus loading affect the water quality of Shell Lake. DLM simulations indicate that large changes in water level (approximately 10 ft) affect the persistence of stratification in the lake. During periods with low water levels, the lake is a well-mixed, polymictic system, with water quality degrading slightly as summer progresses. During periods with high water levels, the lake is more stratified, and phosphorus from internal loading is trapped in the hypolimnion and released later in summer, which results in more extreme seasonality in water quality and better clarity in early summer. Results of eutrophication model simulations using a range in external phosphorus inputs illustrate how water quality in Shell Lake (phosphorus and chlorophyll a concentrations and Secchi depths) responds to changes in external phosphorus loading. Results indicate that a 50-percent reduction in external loading from that measured in 2009 would be required to change phosphorus concentrations from 0.018 milligram per liter (mg/L) (measured in 2009) to 0.012 mg/L (estimated for the mid-1800s from analysis of diatoms in sediment cores). Such reductions in phosphorus loading cannot be accomplished by targeting septic systems or internal loading alone because septic systems contribute only about 3 percent of the phosphorus input to the lake, and internal loading from the sediments of Shell Lake contributes only about 25 percent of phosphorus input. Complete elimination of phosphorus from septic systems and internal loading would decrease the phosphorus concentrations in the lake by 0.003–0.004 mg/L. Therefore, reducing phosphorus concentration in the lake more than by 0.004 mg/L requires decreasing phosphorus loading from surface-water contributions, primarily runoff to the lake. Reconstructed changes in water quality from 1860 to 2010, based on changes in the diatom communities archived in the sediments and eutrophication model simulations, suggest that anthropogenic changes in the watershed (sawmill construction in 1881; the establishment of the village of Shell Lake; and land-use changes in the 1920s, including increased agriculture) had a much larger effect on water quality than the natural changes associated with fluctuations in water level. Although the effects of natural changes in water level on water quality appear to be small, changes in water level do have a modest effect on water quality, primarily manifested as small improvements during higher water levels. Fluctuations in water level, however, have a larger effect on the seasonality of water-quality patterns, with better water quality, especially increased Secchi depths, in early summer during years with high water levels.

Wisconsin↗

Use of the Soil and Water Assessment Tool (SWAT) for simulating hydrology and water quality in the Cedar River Basin, Iowa, 2000--10

The U.S. Geological Survey, in cooperation with the Iowa Department of Natural Resources, used the Soil and Water Assessment Tool to simulate streamflow and nitrate loads within the Cedar River Basin, Iowa. The goal was to assess the ability of the Soil and Water Assessment Tool to estimate streamflow and nitrate loads in gaged and ungaged basins in Iowa. The Cedar River Basin model uses measured streamflow data from 12 U.S. Geological Survey streamflow-gaging stations for hydrology calibration. The U.S. Geological Survey software program, Load Estimator, was used to estimate annual and monthly nitrate loads based on measured nitrate concentrations and streamflow data from three Iowa Department of Natural Resources Storage and Retrieval/Water Quality Exchange stations, located throughout the basin, for nitrate load calibration. The hydrology of the model was calibrated for the period of January 1, 2000, to December 31, 2004, and validated for the period of January 1, 2005, to December 31, 2010. Simulated daily, monthly, and annual streamflow resulted in Nash-Sutcliffe coefficient of model efficiency (E NS ) values ranging from 0.44 to 0.83, 0.72 to 0.93, and 0.56 to 0.97, respectively, and coefficient of determination (R 2 ) values ranging from 0.55 to 0.87, 0.74 to 0.94, and 0.65 to 0.99, respectively, for the calibration period. The percent bias ranged from -19 to 10, -16 to 10, and -19 to 10 for daily, monthly, and annual simulation, respectively. The validation period resulted in daily, monthly, and annual E NS values ranging from 0.49 to 0.77, 0.69 to 0.91, and -0.22 to 0.95, respectively; R 2 values ranging from 0.59 to 0.84, 0.74 to 0.92, and 0.36 to 0.92, respectively; and percent bias ranging from -16 for all time steps to percent bias of 14, 15, and 15, respectively. The nitrate calibration was based on a small subset of the locations used in the hydrology calibration with limited measured data. Model performance ranges from unsatisfactory to very good for the calibration period (January 1, 2000, to December 31, 2004). Results for the validation period (January 1, 2005, to December 31, 2010) indicate a need for an increase of measured data as well as more refined documented management practices at a higher resolution. Simulated nitrate loads resulted in monthly and annual E NS values ranging from 0.28 to 0.82 and 0.61 to 0.86, respectively, and monthly and annual R 2 values ranging from 0.65 to 0.81 and 0.65 to 0.88, respectively, for the calibration period. The monthly and annual calibration percent bias ranged from 4 to 7 and 5 to 7, respectively. The validation period resulted in all but two E NS values less than zero. Monthly and annual validation R 2 values ranged from 0.5 to 0.67 and 0.25 to 0.48, respectively. Monthly and annual validation percent bias ranged from 46 to 68 for both time steps. A daily calibration and validation for nitrate loads was not performed because of the poor monthly and annual results; measured daily nitrate data are available for intervals of time in 2009 and 2010 during which a successful monthly and annual calibration could not be achieved. The Cedar River Basin is densely gaged relative to other basins in Iowa; therefore, an alternative hydrology scenario was created to assess the predictive capabilities of the Soil and Water Assessment Tool using fewer locations of measured data for model hydrology calibration. Although the ability of the model to reproduce measured values improves with the number of calibration locations, results indicate that the Soil and Water Assessment Tool can be used to adequately estimate streamflow in less densely gaged basins throughout the State, especially at the monthly time step. However, results also indicate that caution should be used when calibrating a subbasin that consists of physically distinct regions based on only one streamflow-gaging station.

Iowa↗

Determination of dilution factors for discharge of aluminum-containing wastes by public water-supply treatment facilities into lakes and reservoirs in Massachusetts

Dilution of aluminum discharged to reservoirs in filter-backwash effluents at water-treatment facilities in Massachusetts was investigated by a field study and computer simulation. Determination of dilution is needed so that permits for discharge ensure compliance with water-quality standards for aquatic life. The U.S. Environmental Protection Agency chronic standard for aluminum, 87 micrograms per liter (μg/L), rather than the acute standard, 750 μg/L, was used in this investigation because the time scales of chronic exposure (days) more nearly match rates of change in reservoir concentrations than do the time scales of acute exposure (hours). Whereas dilution factors are routinely computed for effluents discharged to streams solely on the basis of flow of the effluent and flow of the receiving stream, dilution determination for effluents discharged to reservoirs is more complex because (1), compared to streams, additional water is available for dilution in reservoirs during low flows as a result of reservoir flushing and storage during higher flows, and (2) aluminum removal in reservoirs occurs by aluminum sedimentation during the residence time of water in the reservoir. Possible resuspension of settled aluminum was not considered in this investigation. An additional concern for setting discharge standards is the substantial concentration of aluminum that can be naturally present in ambient surface waters, usually in association with dissolved organic carbon (DOC), which can bind aluminum and keep it in solution. A method for dilution determination was developed using a mass-balance equation for aluminum and considering sources of aluminum from groundwater, surface water, and filter-backwash effluents and losses caused by sedimentation, water withdrawal, and spill discharge from the reservoir. The method was applied to 13 reservoirs. Data on aluminum and DOC concentrations in reservoirs and influent water were collected during the fall of 2009. Complete reservoir volume was determined to be available for mixing on the basis of vertical and horizontal aluminum-concentration profiling. Losses caused by settling of aluminum were assumed to be proportional to aluminum concentration and reservoir area. The constant of proportionality, as a function of DOC concentration, was established by simulations in each of five reservoirs that differed in DOC concentration. In addition to computing dilution factors, the project determined dilution factors that would be protective with the same statistical basis (frequency of exceedance of the chronic standard) as dilutions computed for streams at the 7-day-average 10-year-recurrence annual low flow (the 7Q10). Low-flow dilutions are used for permitting so that receiving waters are protected even at the worst-case flow levels. The low-flow dilution factors that give the same statistical protection are the lowest annual 7-day-average dilution factors with a recurrence of 10 years, termed 7DF10s. Determination of 7DF10 values for reservoirs required that long periods of record be simulated so that dilution statistics could be determined. Dilution statistics were simulated for 13 reservoirs from 1960 to 2004 using U.S. Geological Survey Firm-Yield Estimator software to model reservoir inputs and outputs and present-day values of filter-effluent discharge and aluminum concentration. Computed settling velocities ranged from 0 centimeters per day (cm/d) at DOC concentrations of 15.5 milligrams per liter (mg/L) to 21.5 cm/d at DOC concentrations of 2.7 mg/L. The 7DF10 values were a function of aluminum effluent discharged. At current (2009) effluent discharge rates, the 7DF10 values varied from 1.8 to 115 among the 13 reservoirs. In most cases, the present-day (2009) discharge resulted in receiving water concentrations that did not exceed the standard at the 7DF10. Exceptions were one reservoir with a very small area and three reservoirs with high concentrations of DOC. Maximum permissible discharges were determined for water-treatment plants by adjusting discharges upward in simulations until the 7DF10 resulted in reservoir concentrations that just met the standard. In terms of aluminum flux, these discharges ranged from 0 to 28 kilograms of aluminum per day.

Massachusetts↗

Microgravity methods for characterization of groundwater-storage changes and aquifer properties in the karstic Madison aquifer in the Black Hills of South Dakota, 2009-12

A study of groundwater storage in the karstic Madison aquifer in the Black Hills of South Dakota using microgravity methods was conducted by the U.S. Geological Survey in cooperation with West Dakota Water Development District, South Dakota Department of Environment and Natural Resources, and Lawrence County. Microgravity measurements from 2009 to 2012 were used to investigate groundwater-storage changes and effective porosity in unconfined areas of the Madison aquifer. Time-lapse microgravity surveys that use portable high-sensitivity absolute and relative gravimeters indicated temporal-gravity changes as a result of changing groundwater mass. These extremely precise measurements of gravity required characterization and removal of internal instrumental and external environmental effects on gravity from the raw data. The corrected data allowed groundwater-storage volume to be quantified with an accuracy of about plus or minus 0.5 foot of water per unit area of aquifer. Quantification of groundwater-storage change, coupled with water-level data from observation wells located near the focus areas, also was used to calculate the effective porosity at specific altitudes directly beneath gravity stations. Gravity stations were established on bedrock outcrops in three separate focus areas for this study. The first area, the Spring Canyon focus area, is located to the south of Rapid City with one gravity station on the rim of Spring Canyon near the area where Spring Creek sinks into the Madison aquifer. The second area, the Doty focus area, is located on outcrops of the Madison Limestone and Minnelusa Formation to the northwest of Rapid City, and consists of nine gravity stations. The third area, the Limestone Plateau focus area, consists of a single gravity station in the northwestern Black Hills located on an outcrop of the Madison Limestone. An absolute-gravity station, used to tie relative-gravity survey data together, was established on a relatively impermeable bedrock outcrop to minimize groundwater-storage change at the reference location. Data from the three focus areas allow for interpretation of groundwater-storage characteristics using microgravity measurements. Gravity measurements, together with water-level data from an observation well located 2 miles from the Spring Canyon focus area and measured streamflow in Spring Creek, provided evidence that rapid groundwater-storage change, responding to changes in sinking streamflow over the recharge area of the aquifer, occurred in the Madison aquifer directly beneath the gravity station at Spring Canyon. This phenomenon likely was a result of groundwater movement through caverns, conduits, and fractures, which are common in karst aquifers. Spatially and temporally separated microgravity data for the Doty focus area indicated horizontal and vertical heterogeneity of effective porosity for the Madison aquifer. One such example of this was indicated by water-level measurements at an observation well and gravity measurements at four gravity stations in the southeastern part of the Doty area, which were used to estimate effective porosity values ranging from greater than 0 to 0.18. A decrease in groundwater storage determined by microgravity measurements during the spring recharge period for five upgradient stations in the Doty focus area indicated the possibility of rapid release and downgradient cascading of perched groundwater. Evidence for similar phenomena was documented for Wind Cave and Brooks Cave in the Black Hills. Absolute-gravity measurements at the Limestone Plateau focus area confirmed the relation between water levels in an observation well and changes in groundwater storage. Comparison of these gravity measurements with water levels in a nearby observation well resulted in an effective porosity estimate of 0.02 for the Madison aquifer beneath the gravity station.

South Dakota↗

Environmental setting, water budget, and stream assessment for the Broad Run watershed, Chester County, Pennsylvania

The Broad Run watershed lies almost entirely in West Bradford Township, Chester County, Pa., and drains 7.08 square miles to the West Branch Brandywine Creek. Because of the potential effect of encroaching development and other stresses on the Broad Run watershed, West Bradford Township, the Chester County Water Resources Authority, and the Chester County Health Department entered into a cooperative study with the U.S. Geological Survey to complete an annual water budget and stream assessment of overall conditions. The annual water budget quantified the basic parameters of the hydrologic cycle for the climatic conditions present from April 1, 2003, to March 31, 2004. These water-budget data identified immediate needs and (or) deficits that were present within the hydrologic cycle during that period, if present; however, an annual water budget encompassing a single year does not identify long-term trends. The stream assessment was conducted in two parts and assessed the overall condition of the watershed, an overall assessment of the fluvial-geomorphic conditions within the watershed and an overall assessment of the stream-quality conditions. The data collected will document present (2004) conditions and identify potential vulnerabilities to future disturbances. For the annual period from April 1, 2003, to March 31, 2004, determination of an annual water budget indicated that of the 67.8 inches of precipitation that fell on the Broad Run watershed, 38.8 inches drained by way of streamflow to the West Branch Brandywine Creek. Of this 38.8 inches of streamflow, local-minimum hydrograph separation techniques determined that 7.30 inches originated from direct runoff and 31.5 inches originated from base flow. The remaining precipitation went into ground-water storage (1.71 inches) and was lost to evapotranspiration (27.3 inches). Ground-water recharge for this period-35.2 inches-was based on these values and an estimated ground-water evapotranspiration rate of 2 inches. Assessment of fluvial-geomorphic conditions included large-scale mapping of stream classes within the Broad Run watershed and in-depth study of three representative stream reaches also within the Broad Run watershed. Based on the total distance of all stream reaches classified within the Broad Run watershed, 61 percent were classified as C-class, 14 percent as E-class, 13 percent as B-class, 5 percent as F-class, 4 percent as undifferentiated B- and F-class, 2 percent as G-class, and less than 1 percent as A-class. The map of stream classes indicates that the Broad Run watershed currently has no large-scale areas of stream impairment and that, generally, the stream is not entrenched and the main branch of the Broad Run has an available, functioning flood plain. Smaller tributary streams, however, showed signs of localized entrenchment due to site-specific influences such as natural stream-channel evolution, localized channelization, localized contraction due to road and driveway crossings, and (or) increased localized runoff. For example, one small reach along a tributary channel was observed to become entrenched due to runoff originating from a new housing development. Entrenched stream reaches are merely located by large-scale mapping and require individual assessment to determine potential causes of entrenchment and (or) future restorative actions. Three in-depth geomorphic study sites showed that the Broad Run can currently be considered graded or in a state of dynamic equilibrium. The sites did, however, identify certain vulnerabilities to future changes within the watershed. These vulnerabilities included disruption of the present sediment supply, including both increases and (or) reductions in the current sediment loads within the Broad Run; increases in both magnitude and duration of storm-water runoff; encroachment of development onto present flood-plain areas, and (or) alterations to riparian zones. Assessment of stream-quality conditions includ

Scientific Investigations Report↗

Landscape evolution in eastern Chuckwalla Valley, Riverside County, California

This study investigates sedimentary and geomorphic processes in eastern Chuckwalla Valley, Riverside County, California, a region of arid, basin-and-range terrain where extensive solar-energy development is planned. The objectives of this study were to (1) measure local weather parameters and use them to model aeolian sediment-transport potential; (2) identify surface sedimentary characteristics in representative localities; and (3) evaluate long-term landscape evolution rates and processes by analyzing stratigraphy in combination with luminescence geochronology. The new stratigraphic and geochronologic data presented in this report demonstrate the varying local significance of aeolian, alluvial fan, lacustrine (playa), and possibly Colorado River influence over a range of time scales. The dominant sand-transport direction in eastern Chuckwalla Valley is toward the northeast, consistent with the recognized regional west-to-east wind direction. However, occasional strong wind events from the north can transport large quantities of sand southward and temporarily reshape local geomorphic features. Influence of a northwest wind direction is also locally dominant around mountain ranges and controls the modern morphology of the Palen dune field. Modeled sand fluxes are on the order of 10 5 kilograms per meter width per year at the site of weather monitoring, 5 kilometers northwest of the Mule Mountains. Aeolian dunes are locally well developed and actively migrating. Their location and activity are determined largely by sediment supply from playa surfaces and ephemeral stream channels, which also control the dunes’ spatial extent and migration potential; stream channels act as both source and sink for aeolian sediment in this environment. Excavations at five sites along a northwest-to-southeast transect reveal that playa deposits formed around 266–226 thousand years ago south of the McCoy Mountains and immediately north of the present location of Interstate 10. The playa material is overlain by late Pleistocene to Holocene alluvial fan deposits. To the southeast (south of Interstate 10, but north of the Mule Mountains), we identified rapid accumulation of alluvial sediment around the time of the Last Glacial Maximum (23–20 thousand years ago), unconformably overlain by a locally varying assemblage of recent aeolian material or Holocene alluvial fan sediment. We have used stratigraphic characteristics and luminescence ages to calculate accumulation rates for sites in eastern Chuckwalla Valley, and thereby to identify spatial variation in landscape stability over decadal and longer time scales. If future solar-energy development plans are to include natural sand-transport corridors, plans would entail retaining the ability for sand to be transported eastward from the ephemeral stream channels and playas that supply sediment to the dunes, sand sheets, and sand ramps of Chuckwalla Valley, and also to allow for southward transport during episodic strong weather events several times per year. The aeolian sediment-transport corridors are dynamic spatially and temporally, reorganizing on the basis of seasonal changes to wind drift potential. Future landscape stability also will be determined by climate-driven changes to vegetation and thereby to aeolian sediment availability. In a warmer, drier climate, aeolian sediment activity is expected to increase, owing to a decrease in stabilizing vegetation cover and more extreme rain that supplies sediment to ephemeral stream channels and playas from which it is remobilized by wind.

California↗

Numerical simulation of groundwater flow in the Columbia Plateau Regional Aquifer System, Idaho, Oregon, and Washington

A three-dimensional numerical model of groundwater flow was constructed for the Columbia Plateau Regional Aquifer System (CPRAS), Idaho, Oregon, and Washington, to evaluate and test the conceptual model of the system and to evaluate groundwater availability. The model described in this report can be used as a tool by water-resource managers and other stakeholders to quantitatively evaluate proposed alternative management strategies and assess the long‑term availability of groundwater. The numerical simulation of groundwater flow in the CPRAS was completed with support from the Groundwater Resources Program of the U.S. Geological Survey Office of Groundwater. The model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model, MODFLOW-NWT. The model uses 3-kilometer (9,842.5 feet) grid cells that subdivide the model domain by 126 rows and 131 columns. Vertically, the model domain was subdivided into six geologic model units. From youngest to oldest, the units are the Overburden, the Saddle Mountains Basalt, the Mabton Interbed, the Wanapum Basalt, the Vantage Interbed, and the Grande Ronde Basalt. Natural recharge was estimated using gridded historical estimates of annual precipitation for the period 1895–2007. Pre-development recharge was estimated to be the average natural recharge for this period. Irrigation recharge and irrigation pumping were estimated using a remote-sensing based soil-water balance model for the period 1985–2007. Pre-1985 irrigation recharge and pumping were estimated using previously published compilation maps and the history of large-scale irrigation projects. Pumping estimates for municipal, industrial, rural, residential, and all other uses were estimated using reported values and census data. Pumping was assumed to be negligible prior to 1920. Two models were constructed to simulate groundwater flow in the CPRAS: a steady-state predevelopment model representing conditions before large-scale pumping and irrigation altered the system, and a transient model representing the period 1900–2007. Automated parameter-estimation techniques (steady-state predevelopment model) and traditional trial-and-error (transient model) methods were used for calibration. To calibrate the steady-state and transient models, 10,525 and 46,460 water level measurements, respectively, and 50 base-flow estimates were used. The steady-state model simulated the shape, slope, and trends of a potentiometric surface that was generally consistent with mapped water levels. For the transient model, the mean and median difference between simulated and measured hydraulic heads is -10 and 4 ft, respectively, with a standard deviation of 164 ft over a 5,648 ft range of measured heads. The residuals for the simulation period show that 52 percent of the simulated heads exceeded measured heads with a median residual value of 43 ft, and 48 percent were less than measured heads with a median residual value of -76 ft. The CPRAS model was constructed to derive components of the groundwater budget and help understand the interactions of stresses, such as recharge, groundwater pumping, and commingling wells on the groundwater and surface-water system. Through these applications, the model can be used to identify trends in groundwater storage and use, and quantify groundwater availability. The annual groundwater budgets showed several patterns of change over the simulation period. Groundwater pumping was negligible until the 1950s and began to increase significantly during the 1970s and 1980s. Recharge was highly variable due to the interannual variability of precipitation, but began to increase in the late 1940s due to the increase in surface-water irrigation projects. Groundwater contributions to streamflow (base flow) followed recharge closely. However, in areas of significant groundwater-level decline, base flow is reduced. Groundwater pumping had the greatest effect on water levels, followed by irrigation enhanced recharge. Commingling was a larger factor in structurally complex upland areas where hydraulic-head gradients are naturally high. Groundwater pumping has increased substantially over the past 40–50 years; this increase resulted in declining water levels at depth and decreased base flows over much of the study area. The effects of pumping are mitigated somewhat by the increase of surface-water irrigation, especially in the shallow Overburden unit, and commingling wells in some areas. During dry to average years, groundwater pumping causes a net loss of groundwater in storage and current condition (2000–2007) groundwater pumping exceeds recharge in all but the wettest of years.

Idaho, Oregon, Washington↗

Borehole analysis, single-well aquifer testing, and water quality for the Burnpit well, Mount Rushmore National Memorial, South Dakota

Mount Rushmore National Memorial (hereafter referred to as “the memorial”), in western South Dakota, is maintained by the National Park Service (NPS) and includes 1,278 acres of land in the east-central part of the Black Hills. An ongoing challenge for NPS managers at the memorial is providing water from sustainable and reliable sources for operations, staff, and the increasing number of visitors. In 2020, the U.S. Geological Survey (USGS) and NPS completed a hydrological study of the Burnpit well (well 5), a 580-foot-deep open hole groundwater well completed in metamorphic (crystalline) rock at the memorial. The purpose of this study was to estimate the geological and hydraulic properties of the aquifer supplying the well and to determine the water quality of the groundwater from the well. The study provides NPS staff and managers background information for assessing future uses for the well. Methods for data collection and analysis for the study included borehole and video camera analysis in 2020, aquifer testing by the NPS in 2009 and the USGS in 2020, and water-quality sampling in 2020. Borehole camera video generally matched the lithology recorded in the well log. Fractures recorded in the well log and observed with the borehole camera, including more than 20 less prominent fractures and rough sidewall areas, indicated a fractured aquifer. The fractures are the primary conduits for groundwater flow through the rock and into the well. Transmissivity was estimated for the upper and lower water-level drawdown zones at the Burnpit well with data from the NPS and USGS using the Theis and Cooper-Jacob methods. Transmissivity for the NPS test using the Theis method was 9.0 and 11 feet squared per day (ft 2 /d) for the upper and lower drawdown zones, respectively. Using the Cooper-Jacob method, the transmissivity was 22 and 14 ft 2 /d for the upper and lower drawdown zones of the aquifer, respectively. Transmissivity estimates from data from the USGS test were similar. The Theis method, applied to the upper and lower drawdown zones of the aquifer, produced transmissivity estimates of 7.7 and 10 ft 2 /d, and the Cooper-Jacob method produced estimates of 9.7 and 12 ft 2 /d, respectively. Storativity (specific yield) estimated using the Theis method for the NPS aquifer-test data was 0.85 and 0.92 for the upper and lower drawdown zones of the aquifer, respectively. The Cooper-Jacob method applied to the NPS aquifer-test data produced storativity estimates of 0.11 and 0.50 for the upper and lower drawdown zones, respectively. The Theis method applied to the USGS aquifer-test data estimated storativity values of 0.77 and 1.0 for the upper and lower drawdown zones, respectively. The Cooper-Jacob method estimated storativity of 0.50 and 0.60 for the upper and lower drawdown zones of the USGS aquifer test, respectively. The estimated storativity values from the NPS and USGS aquifer tests for the upper and lower drawdown zones were higher than expected for limestones and schists. The hypothetical equilibrium drawdown for the Burnpit well was estimated after the NPS test in 2009 at no more, and possibly less, than 35 gallons per minute. The NPS noted that the sustainable yield likely was overestimated because the water level did not stabilize during the NPS aquifer test. The specific capacity for the NPS aquifer test in 2009 was 0.16 gallon per minute per foot ([gal/min]/ft) of drawdown at 3 hours, and the specific capacity for the USGS aquifer test in 2020 was 0.13 (gal/min)/ft of drawdown at 3 hours. The rate of water-level recovery after pumping ceased was 0.017 and 0.013 (gal/min)/ft for the NPS and USGS aquifer tests, respectively. The water-level recovery rate was nearly an order of magnitude less than the specific capacity estimated during pumping, indicating that water levels in the Burnpit well may not recover quickly enough during pumping to provide for a continuous source of water. Water-quality samples were collected at the Burnpit well on June 24 and July 23, 2020, and analyzed for field-measured properties, major ions, metals, nutrients, and perchlorate. Iron, zinc, and lithium concentrations for unfiltered samples in the well were at least three times greater than the mean filtered sample concentrations reported for crystalline aquifers in the Black Hills. Manganese concentrations were less than the mean concentration for crystalline aquifers but exceeded the U.S. Environmental Protection Agency (EPA) secondary drinking-water standards. The iron concentration from the June 24 sample was about 11 times greater than the EPA secondary drinking-water standards and mean concentrations from crystalline aquifers in the Black Hills. Arsenic concentrations in Burnpit well samples collected in 2020 were greater than the EPA primary drinking-water standard and the mean concentration for crystalline aquifers in the Black Hills. Arsenic occurs naturally in the rock of crystalline aquifers, and concentrations from samples in the Black Hills commonly exceed the EPA primary drinking-water standard of 10 micrograms per liter. High concentrations of arsenic, iron, and manganese metals in the Burnpit well make groundwater from the well in its natural state unusable as a drinking-water source, and water treatment would be necessary to reduce the trace element concentrations to less than the EPA primary and secondary drinking-water standards. However, if the memorial has immediate nonpotable water requirements, such as for construction and fire suppression, groundwater from the Burnpit well could provide water without causing additional stress to current (2021) drinking-water sources.

South Dakota↗

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia↗

Fluctuations in groundwater levels related to regional and local withdrawals in the fractured-bedrock groundwater system in northern Wake County, North Carolina, March 2008-February 2009

A study of dewatering of the fractured-bedrock aquifer in a localized area of east-central North Carolina was conducted from March 2008 through February 2009 to gain an understanding of why some privately owned wells and monitoring wells were intermittently dry. Although the study itself was localized in nature, the resulting water-resources data and information produced from the study will help enable resource managers to make sound water-supply and water-use decisions in similar crystalline-rock aquifer setting in parts of the Piedmont and Blue Ridge Physiographic Provinces. In June 2005, homeowners in a subdivision of approximately 11 homes on lots approximately 1 to 2 acres in size in an unincorporated area of Wake County, North Carolina, reported extremely low water pressure and temporarily dry wells during a brief period. This area of the State, which is in the Piedmont Physiographic Province, is undergoing rapid growth and development. Similar well conditions were reported again in July 2007. In an effort to evaluate aquifer conditions in the area of intermittent water loss, a study was begun in March 2008 to measure and monitor water levels and groundwater use. During the study period from March 2008 through February 2009, regular dewatering of the fractured-bedrock aquifer was documented with water levels in many wells ranging between 100 and 200 feet below land surface. Prior to this period, water levels from the 1980s through the late 1990s were reported to range from 15 to 50 feet below land surface. The study area includes three community wells and more than 30 private wells within a 2,000-foot radius of the dewatered private wells. Although groundwater levels were low, recovery was observed during periods of heavy rainfall, most likely a result of decreased withdrawals owing to less demand for irrigation purposes. Similar areal patterns of low groundwater levels were delineated during nine water-level measurement periods from March 2008 through February 2009. Correlation of groundwater-level distribution patterns with orientations of geologic structures obtained from surficial mapping, borehole geophysical measurements, and interpretation of fracture traces suggests two dominant trends striking north-south and N. 65 degrees W. A variation in overall response to groundwater withdrawals was noted in the continuous groundwater-level records for the monitored observation wells and dewatered private wells. The largest overall declines during the study period were observed in an observation well in which the water-level declined as much as 247 feet from mid-July through early August 2008, during a period of heavy usage. A private well had a water-level decline of about 94 feet during the same monitoring period. The large declines recorded in the observation well and the private well indicated a substantial temporary loss of storage in the fractured-bedrock aquifer near the wells, thus reducing the amount of water available to shallow wells in the area (those wells with total depths of about 300 feet), and resulting in temporary well failures until such time as the aquifer recovered.

North Carolina↗

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗