USGS ScienceSearch

SEARCH · USGS Science

Results for “Water Surface Techniques”

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 127 records · Page 7Linked to original sources

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California

Numerical-simulation and conjunctive-management models of the Hunt-Annaquatucket-Pettaquamscutt stream-aquifer system, Rhode Island

Numerical-simulation and optimization techniques were used to evaluate alternatives for the conjunctive management of ground- and surface-water resources of the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system in central Rhode Island. Ground-water withdrawals from the Hunt-Annaquatucket-Pettaquamscutt aquifer exceeded 8 million gallons per day during months of peak water use during 1993­98, and additional withdrawals have been proposed to meet growing demands from within and outside of the system boundary. The system is defined by the Hunt­-Annaquatucket­-Pettaquamscutt aquifer, which is composed of glacial stratified deposits, and the network of rivers, brooks, and ponds that overlie and are in hydraulic connection with the aquifer. Nearly all of the water withdrawn, however, is derived from depletions of flow in the rivers, brooks, and ponds that overlie the aquifer. Streamflow depletions are of concern to environmental agencies because of the adverse effects that reductions in streamflow can have on aquatic and riparian ecosystems. A conjunctive-management model of the stream-aquifer system was developed to simultaneously address the water-demand and streamflow-depletion issues. The objective of the model was to maximize total ground-water withdrawal from the aquifer during July, August, and September. These three months are generally the time of year when water-supply demands are largest and streamflows are simultaneously lowest. Total withdrawal from the aquifer was limited by a set of constraints specified in the model. These constraints were (1) maximum rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers; (2) minimum monthly water demands of each of three water-supply systems that withdraw water from the aquifer; and (3) minimum and maximum withdrawal rates at each supply well. The conjunctive-management model was formulated mathematically as a linear program. The model was solved by a response-matrix technique that incorporates the results of transient, numerical simulation of the stream-aquifer system into the constraint set of the linear program. The basis of the technique was the assumption that streamflow-depletion rates in each river were a linear function of ground-water-withdrawal rates at each well. This assumption was shown to be valid for the conditions evaluated in this study, primarily because of the very high transmissivity of the aquifer near many of the wells pumped for water supply. A transient, numerical model of the system was developed to simulate an average annual cycle of monthly withdrawal and hydrologic conditions representative of the 56-year period 1941­96. The transient model was used to generate characteristic streamflow-depletion responses in each river to simulated withdrawals at each well; these characteristic responses, or response coefficients, were then incorporated directly into the streamflow-depletion constraints of the linear program. Four sets of applications of the conjunctive-management model were made to determine whether total ground-water withdrawal from the aquifer during July, August, and September could be increased over the current total withdrawal for alternative definitions of the maximum rates of streamflow depletion allowed in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Current conditions were defined as the average monthly withdrawal rates at each supply well, water demands of each of the three water-supply systems, and estimated streamflow-depletion rates during the 6-year period 1993­98. Total withdrawal from all wells in the system from July through September during 1993­98 was 506.5 million gallons. Estimated streamflow-depletion rates for 1993­98 were calculated by use of the transient model, with the 1993­98 average monthly withdrawal rates specified at each supply well. Streamflow-depletion rates calculated for July, August, and September averaged 25 percent of the model-calculated pre-withdrawal streamflow rates for the Hunt River, 19 percent for the Annaquatucket River, and 7 percent for the Pettaquamscutt River. The first set of applications of the model were made with the current estimated rates of streamflow depletion in the Hunt, Annaquatucket, and Pettaquamscutt Rivers. Results of these applications indicated that total withdrawal from the aquifer during July, August, and September could be increased from about 8 to 18 percent (from 546.0 to 596.3 million gallons) over the current total withdrawal. The increased withdrawal would require modifications to the current annual withdrawal schedule of each supply well and, for the 18-percent increase, a modified network of supply wells that would include two new wells in the Annaquatucket River Basin. A second set of model applications then was made to determine if current estimated rates of streamflow depletion in the Hunt River could be reduced without increasing current estimated rates of streamflow depletion in the Annaquatucket or Pettaquamscutt Rivers. Decreases in the current rates of streamflow depletion in the Hunt River would result in increased streamflow in the river during these three months. Results showed that current rates of streamflow depletion in the Hunt River during July, August, and September could be decreased from 5 to 15 percent, depending on whether the existing or modified well network was used. Subsequent model applications indicated that substantial increases in total ground-water withdrawal from the aquifer are possible, but would require increased rates of streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers. Maximum increases in the July through September withdrawal from the aquifer of about 39 to 50 percent (from 705.1 to 760.3 million gallons) over the current total withdrawal were calculated when streamflow-depletion rates in the Annaquatucket and Pettaquamscutt Rivers were allowed to increase from current estimated rates to a maximum of 25 percent of the model-calculated pre-withdrawal streamflow for each river during July, August, and September. Alternatively, it was shown that current estimated rates of streamflow depletion in the Hunt River during July, August, and September could be reduced by as much as 35 percent for the maximum allowed increases in streamflow depletion in the Annaquatucket and Pettaquamscutt Rivers; maximum increased withdrawal from the aquifer, however, would range from 8 to 18 percent over the current total withdrawal for the 35-percent reduction in streamflow-depletion rates in the Hunt River. Results of the different applications of the model demonstrate the usefulness of coupling numerical-simulation and optimization techniques for regional-scale evaluation of water-resource management alternatives. The results of the evaluation must be viewed, however, within the limitations of the quality of data available for the Hunt­-Annaquatucket­-Pettaquamscutt stream-aquifer system and representation of the system by a simulation model. An additional limitation of the analysis was the use of an average annual cycle of monthly withdrawal and hydrologic conditions. Ground-water withdrawal strategies may need to be modified to meet streamflow-depletion constraints during extreme hydrologic events, such as droughts. Contributing areas and sources of water to the supply wells also were delineated by use of a steady-state model of the stream-aquifer system. The model was developed to simulate long-term-average ground-water flow and ground-water/ surface-water interactions in the system during the 56-year period 1941­96. Sources of water to the wells consisted of precipitation and wastewater recharge to the aquifer, streamflow leakage from natural stream-channel losses, streamflow leakage caused by induced infiltration, and lateral ground-water inflow from till and bedrock upland areas.

Rhode Island

Effects of surface-water and groundwater inflows and outflows on the hydrology of the Tsala Apopka Lake Basin in Citrus County, Florida

The U.S. Geological Survey, in cooperation with the Southwest Florida Water Management District, initiated a study to quantify the inflows and outflows in the Floral City, Inverness, and Hernando pools of the Tsala Apopka Lake Basin in Citrus County, Florida. This study assesses hydrologic changes in pool stages, groundwater levels, spring flows, and streamflows caused by the diversion of streamflow from the Withlacoochee River to the Tsala Apopka Lake Basin through water-control structures. A surface-water/groundwater flow model was developed using hydraulic parameters for lakes, streams, the unsaturated zone, and the underlying surficial and Upper Floridan aquifers estimated using an inverse modeling calibration technique. After calibration, the model was used to assess the relation between inflows and outflows in the Tsala Apopka Lake Basin and changes in pool stages. Simulation results using the calibrated surface-water/groundwater flow model showed that leakage rates from the pools to the Upper Floridan aquifer were largest at the deep lake cells and that these leakage rates to the Upper Floridan aquifer were the highest in the model area. Downward leakage to the Upper Floridan aquifer occurred beneath most of the extent of the Floral City, Inverness, and Hernando pools. These leakage rates depended on the lakebed leakance and the difference between lake stages and heads in the Upper Floridan aquifer. Leakage rates were higher for the Floral City pool than for the Inverness pool, and higher for the Inverness pool than for the Hernando pool. Lakebed leakance was higher for the Floral City pool than for the Hernando pool, and higher for the Hernando pool than for the Inverness pool. Simulation results showed that the average recharge rate to the surficial aquifer was 10.3 inches per year for the 2004 to 2012 simulation period. Areas that recharge the surficial aquifer covered about 86 percent of the model area. Simulations identified areas along segments of the Withlacoochee River and within land-surface depressions that receive water from the surficial aquifer. Recharge rates were largest in physiographic regions having a deep water table. Simulated heads in the Upper Floridan aquifer indicated the general flow directions in the active flow model area were from the northeast toward the southwest and then westward toward the coast, and from the southeast toward the northwest and then westward toward the coast, consistent with flow directions inferred from the estimated potentiometric surface map for May 2010. The largest inflow in the water budget of the Upper Floridan aquifer was downward leakage from the overlying hydrogeologic unit. The largest outflow in the water budget of the Upper Floridan aquifer was spring flow. The calibrated surface-water and groundwater flow model was used to simulate hydrologic scenarios that included changes in rainfall rates, projected increases in groundwater pumping rates for 2025 and 2035, no flow for the 2004–12 period through the eight water-control structures in the Tsala Apopka Lake Basin, and the removal of the Inglis Dam and the Inglis Bypass Spillway on Lake Rousseau. Scenario simulation results were compared to annual average calibrated water levels and flows from 2004 to 2012. Simulated declines in the Tsala Apopka Lake pool stages under the 10-percent lower rainfall scenario were about 0.8, 0.3, and 1.3 feet (ft) for the Floral City, Inverness, and Hernando pools, respectively. Simulated groundwater levels under the same scenario declined up to 5.4 ft in the surficial aquifer and up to 2.9 ft in the Upper Floridan aquifer. Under the projected increases in groundwater pumping rates for 2035 that represented an increase of 36 percent from average 2004 to 2012 pumping rates, the simulated declines in the Floral City, Inverness, and Hernando pool stages were, in downstream order, 0.02, 0.06, and 0.04 ft. The largest drawdown under the projected increases in groundwater pumping rates for 2035 was 2.1 ft in the surficial aquifer and about 1.8 ft in the Upper Floridan aquifer. A scenario of decreased rainfall by 10 percent caused greater declines in water levels and pool stages than projected increases in groundwater pumping rates. The simulation with no flow through the eight Tsala Apopka Lake water-control structures resulted in simulated declines in average pool stage of 1.8, 1.9, and 0.5 ft in the Floral City, Inverness, and Hernando pools, respectively. The simulated removal of the two water-control structures in Lake Rousseau caused flow to increase at Rainbow Springs by 28 cubic feet per second, an increase of 4.7 percent from the average calibrated flow for 2004 to 2012.

Florida

Analysis of temperature profiles for investigating stream losses beneath ephemeral channels

Continuous estimates of streamflow are challenging in ephemeral channels. The extremely transient nature of ephemeral streamflows results in shifting channel geometry and degradation in the calibration of streamflow stations. Earlier work suggests that analysis of streambed temperature profiles is a promising technique for estimating streamflow patterns in ephemeral channels. The present work provides a detailed examination of the basis for using heat as a tracer of stream/groundwater exchanges, followed by a description of an appropriate heat and water transport simulation code for ephemeral channels, as well as discussion of several types of temperature analysis techniques to determine streambed percolation rates. Temperature‐based percolation rates for three ephemeral stream sites are compared with available surface water estimates of channel loss for these sites. These results are combined with published results to develop conclusions regarding the accuracy of using vertical temperature profiles in estimating channel losses. Comparisons of temperature‐based streambed percolation rates with surface water‐based channel losses indicate that percolation rates represented 30% to 50% of the total channel loss. The difference is reasonable since channel losses include both vertical and nonvertical component of channel loss as well as potential evapotranspiration losses. The most significant advantage of the use of sediment‐temperature profiles is their robust and continuous nature, leading to a long‐term record of the timing and duration of channel losses and continuous estimates of streambed percolation. The primary disadvantage is that temperature profiles represent the continuous percolation rate at a single point in an ephemeral channel rather than an average seepage loss from the entire channel.

Water Resources Research

Surface-geophysical techniques used to detect existing and infilled scour holes near bridge piers

Surface-geophysical techniques were used with a position-recording system to study riverbed scour near bridge piers. From May 1989 to May 1993, Fathometers, fixed- and swept-frequency continuous seismic- reflection profiling (CSP) systems, and a ground-penetrating radar (GPR) system were used with a laser-positioning system to measure the depth and extent of existing and infilled scour holes near bridge piers. Equipment was purchased commercially and modified when necessary to interface the components and (or) to improve their performance. Three 200-kHz black-and-white chart-recording Fathometers produced profiles of the riverbed that included existing scour holes and exposed pier footings. The Fathometers were used in conjunction with other geophysical techniques to help interpret the geophysical data. A 20-kHz color Fathometer delineated scour-hole geometry and, in some cases, the thickness of fill material in the hole. The signal provided subbottom information as deep as 10 ft in fine-grained materials and resolved layers of fill material as thin as 1 ft thick. Fixed-frequency and swept-frequency CSP systems were evaluated. The fixed-frequency system used a 3.5-, 7.0-, or 14-kHz signal. The 3.5-kHz signal penetrated up to 50 ft of fine-grained material and resolved layers as thin as 2.5-ft thick. The 14-kHz signal penetrated up to 20 ft of fine-grained material and resolved layers as thin as 1-ft thick. The swept-frequency systems used a signal that swept from 2- to 16-kHz. With this system, up to 50 ft of penetration was achieved, and fill material as thin as 1 ft was resolved. Scour-hole geometry, exposed pier footings, and fill thickness in scour holes were detected with both CSP systems. The GPR system used an 80-, 100-, or 300-megahertz signal. The technique produced records in water up to 15 ft deep that had a specific conductance less than 200 µS/cm. The 100-MHz signal penetrated up to 40 ft of resistive granular material and resolved layers as thin as 2-ft thick. Scour-hole geometry, the thickness of fill material in scour holes, and riverbed deposition were detected using this technique. Processing techniques were applied after data collection to assist with the interpretation of the data. Data were transferred from the color Fathometer, CSP, and GPR systems to a personal computer, and a commercially available software package designed to process GPR data was used to process the GPR and CSP data. Digital filtering, predictive- deconvolution, and migration algorithms were applied to some of the data. The processed data were displayed and printed as color amplitude or wiggle-trace plots. These processing methods eased and improved the interpretation of some of the data, but some interference from side echoes from bridge piers and multiple reflections remained in the data. The surface-geophysical techniques were applied at six bridge sites in Connecticut. Each site had different water depths, specific conductance, and riverbed materials. Existing and infilled scour holes, exposed pier footings, and riverbed deposition were detected by the surveys. The interpretations of the geophysical data were confirmed by comparing the data with lithologic and (or) probing data.

Connecticut

Deep aquifer recharge in the Columbia River Basalt Group, upper Umatilla River Basin, northeastern Oregon

Groundwater is an important component of the water resources of the upper Umatilla River Basin of northeastern Oregon. As such, understanding the capacity of the resource is vital. Past studies have estimated recharge in the study area. One recent study of the upper Umatilla River Basin indicated that about 80 percent of recharge entering the groundwater system is discharged to streams in the study area through shallow groundwater-flow paths, leaving about 20 percent of recharge to infiltrate deeper parts of the aquifer system. The purpose of this work is to quantify the spatial distribution and variability of deep aquifer recharge in the study area and to understand the reasons for a relatively low percentage of total recharge reaching the deeper parts of the groundwater-flow system. The study area is divided into two distinct physiographic regions—the highly dissected Blue Mountains and the lowland plains. Underlying both regions of the study area are basalts of the Columbia River Basalt Group (CRBG), which is the principal aquifer in the study area. Deep incision by streams in the Blue Mountains disrupts the lateral continuity of the CRBG aquifer units, and infiltrating water is more readily diverted laterally and discharged to streams and springs. In the lowland plains, incision is less pronounced. The shallow CRBG units might be disrupted, but deeper aquifer units retain their lateral continuity and enable groundwater to infiltrate deeper and flow laterally farther downgradient before discharging. Recharge to the deep basalt aquifers is estimated as the difference between total recharge and base flow. Total recharge is the portion of precipitation and applied irrigation water that infiltrates past the root zone to become groundwater recharge. Of this total recharge, a proportion discharges to springs and streams in the study area, and the remaining water infiltrates below the base level of streams and recharges the deep basalt aquifers and contributes to the regional groundwater flow system. The portion of total recharge that recharges the regional flow system is referred to as deep aquifer recharge. Total recharge is the portion of precipitation and applied irrigation water that infiltrates past the root zone to become groundwater recharge. It is the sum of recharge from precipitation and recharge from infiltration of irrigation water. Recharge from precipitation was calculated using a regression method developed for the Columbia Plateau. Recharge from infiltrating irrigation water was obtained from a water balance model developed for the Columbia Plateau. Base flow, the component of streamflow that represents groundwater discharge as opposed to runoff from the land surface, was estimated using the Base Flow Index Modified (BFI-Modified) method, an empirical hydrograph separation technique. Base flow was estimated in eight subbasins with streamgages within the study area. Five of the eight subbasins in which base flow was estimated had permitted water rights for irrigation that specified surface water as the primary source of water. Maximum surface-water withdrawal for irrigation was estimated for all subbasins in which water rights for irrigation occur. The base-flow estimate from BFI-Modified is assumed to be the minimum amount of base flow. The sum of the BFIModified base-flow estimate and the maximum permitted surface-water withdrawal estimate for each subbasin is assumed to be the maximum amount of base flow at the streamgage. These minimum and maximum estimates of base flow were used to calculate minimum and maximum values of deep aquifer recharge in each subbasin analyzed within the study area. Subbasin estimates were scaled up to the Blue Mountains and lowland plains regions, and to the entire study area. Mean annual total recharge for 1981–2010 in the subbasins, analyzed as part of this work, ranged from 6 inches (in.) in the Patawa and Wildhorse Creek subbasins in the lowland plains to as much as 20 in. in the Umatilla River above Meacham Creek subbasin. Mean annual total recharge totaled 4 in. in the lowland plains region and 14 in. in the Blue Mountains. Mean annual total recharge for the entire study area was 11 in. Mean annual base flow ranged from 1 in. in the Patawa and Wildhorse Creek subbasins in the lowland plains to as much as 14 in. in the Umatilla River above Meacham Creek subbasin in the Blue Mountains. Mean annual deep aquifer recharge ranged from 4 in. in the Patawa and Wildhorse Creek subbasins in the lowland plains to as much as 8 in. in the Isqu’ulktpe Creek subbasin in the Blue Mountains. Deep aquifer recharge was 3–4 in. in the lowland plains region and 6 in. in the Blue Mountains. Over the entire study area, mean annual deep aquifer recharge was 5 in. Most groundwater recharge (both total and deep aquifer) in the study area occurred in the Blue Mountains, which highlights the importance of the Blue Mountains as the principal source of groundwater for the study area and for aquifers farther downgradient. Total recharge in the Blue Mountains region represents 86 percent of the mean annual total recharge in the study area in an area that encompasses 65 percent of the study area. However, only 43–44 percent of the mean annual total recharge remains in the system to recharge the deeper, regional aquifer system because the rest is discharged as base flow within the Blue Mountains region. Within the lowland plains region of the study area, an estimated 67–84 percent of the mean annual total recharge remains in the system to recharge the deep, regional aquifer system. Although total recharge in the study area represents only 14 percent of the total recharge across the study area, it contributes 20–24 percent of the water to the deep aquifer. The difference in the percentage of deep groundwater recharge in the Blue Mountains and the lowland plains is attributed to differences in the degree of stream incision. Stream channels are more incised in the Blue Mountains region than they are in the lowland plains. The dissection of the landscape in the Blue Mountains disrupts the lateral continuity of the CRBG aquifer units and allows groundwater to discharge to springs and streams rather than infiltrate more deeply. In the lowland plains region, incision is much less pronounced and deeper CRBG units likely retain their lateral continuity, enabling groundwater to infiltrate more deeply than in the Blue Mountains.

Oregon

Using natural distributions of short-lived radium isotopes to quantify groundwater discharge and recharge

Radium activity in pore water of wetland sediments often differs from the amount expected from local production, decay, and exchange with solid phases. This disequilibrium results from vertical transport of radium with groundwater that flows between the underlying aquifer and surface water. In situations where groundwater recharge or discharge is significant, the rate of vertical water flow through wetland sediment can be determined from the radium disequilibrium by a combined model of transport, production, decay, and exchange with solid phases. We have developed and tested this technique at three sites in the freshwater portion of the Everglades by quantifying vertical advective velocities in areas with persistent groundwater recharge or discharge and estimating a coefficient of dispersion at a site that is subject to reversals between recharge and discharge. Groundwater velocities (v) were determined to be between 0 and -0.5 cm d-1 for a recharge site and 1.5 ?? 0.4 cm d-1 for a discharge site near Levee 39 in the Everglades. Strong gradients in 223Ra and 224Ra usually occurred at the base of the peat layer, which avoided the problems of other tracers (e.g., chloride) for which greatest sensitivity occurs near the peat surface - a zone readily disturbed by processes unrelated to groundwater flow. This technique should be easily applicable to any wetland system with different production rates of these isotopes in distinct sedimentary layers or surface water. The approach is most straightforward in systems where constant pore-water ionic strength can be assumed, simplifying the modeling of radium exchange.

Limnology and Oceanography

Estimating potential wetland extent along selected river reaches in Indiana using streamflow statistics and flood-inundation mapping techniques

In this study potential wetland extents were estimated for 12 river reaches covering about 750 river miles in Indiana and parts of Illinois and Ohio. The study was completed by the U.S. Geological Survey in cooperation with the U.S. Department of Agriculture, Natural Resources Conservation Service. This study follows and adds to the work completed in a pilot study and determines that potential wetland extents can be estimated using streamflow statistics, streamgage data, and flood-inundation mapping techniques. The study was designed to assist in the Agricultural Conservation Easement Program. The Agricultural Conservation Easement Program is a voluntary program administered by the Natural Resources Conservation Service that provides technical and financial assistance to private landowners and Tribes to restore, protect, and enhance wetlands in exchange for retiring eligible land from agriculture. For a site to be eligible for wetland restoration, it should be in a zone with sustained or frequent flooding. This study calculated the flows that lasted for a period of 7 consecutive days on average at least once every 2 years (a value termed the “7MQ2”) for all the U.S. Geological Survey streamgages within the selected river reaches. These 7MQ2 flows were related to the stage-discharge tables for each streamgage, and a corresponding water-surface elevation was determined. Maps of estimated wetland extent were prepared using the 7MQ2 inundation elevation data in conjunction with bare-earth land-surface elevation data made publicly available through the online geospatial data clearinghouses of Indiana, Illinois, and Ohio. Flood-inundation mapping techniques were applied with the aid of geographic information system software to generate water-surface planes that represent inundation elevations associated with the 7MQ2 streamflow. Land-surface elevation data from high-resolution digital elevation models were subtracted from the water-surface planes to produce maps of wetland extent. The 12 map products, including datasets and geoprocessing tools, produced from this study will aid the National Resources Conservation Service and its partners with the onsite inundation-zone verification in agricultural land for potential restoration.

Illinois, Indiana, Ohio

Trench infiltration for managed aquifer recharge to permeable bedrock

Managed aquifer recharge to permeable bedrock is increasingly being utilized to enhance resources and maintain sustainable groundwater development practices. One such target is the Navajo Sandstone, an extensive regional aquifer located throughout the Colorado Plateau of the western United States. Spreading-basin and bank-filtration projects along the sandstone outcrop's western edge in southwestern Utah have recently been implemented to meet growth-related water demands. This paper reports on a new cost-effective surface-infiltration technique utilizing trenches for enhancing managed aquifer recharge to permeable bedrock. A 48-day infiltration trench experiment on outcropping Navajo Sandstone was conducted to evaluate this alternative surface-spreading artificial recharge method. Final infiltration rates through the bottom of the trench were about 0.5 m/day. These infiltration rates were an order of magnitude higher than rates from a previous surface-spreading experiment at the same site. The higher rates were likely caused by a combination of factors including the removal of lower permeability soil and surficial caliche deposits, access to open vertical sandstone fractures, a reduction in physical clogging associated with silt and biofilm layers, minimizing viscosity effects by maintaining isothermal conditions, minimizing chemical clogging caused by carbonate mineral precipitation associated with algal photosynthesis, and diminished gas clogging associated with trapped air and biogenic gases. This pilot study illustrates the viability of trench infiltration for enhancing surface spreading of managed aquifer recharge to permeable bedrock. ?? 2010.

Hydrological Processes

Transport and transformations of chlorinated-solvent contamination in a saprolite and fractured rock aquifer near a former wastewater-treatment plant, Greenville, South Carolina

The transport and fate of chlorinated-ethene contamination was investigated in a fractured-rock aquifer downgradient from a wastewater-treatment plant at a gas-turbine manufacturing facility in Greenville, South Carolina. A vapor-diffusion-sampler technique, developed for this investigation, located fracture zones that discharged contaminated ground water to surface water. The distribution of chlorinated compounds and sulfate, comparison of borehole geophysical data, driller's logs, and the aquifer response to pumpage allowed subsurface contaminant-transport pathways to be delineated.The probable contaminant-transport pathway from the former aeration lagoon was southward. The probable pathway of contaminant transport from the former sludge lagoon was southward to and beneath Little Rocky Creek. South of the creek, the major pathway of contaminant transport appeared to be at a depth of approximately 80 to 107 feet below land surface. The contaminant-transport pathway from the former industrial lagoon was not readily discernible from existing data. A laboratory investigation, as well as examination of ground- water-chemistry data collected during this investigation and concentrations of chlorinated compounds collected during previous investigations,indicates that higher chlorinated compounds are being degraded to lower-chlorinated compounds in the contaminated aquifer. The approaches used in this investigation, as well as the findings, have potential application to other fractured-rock aquifers contaminated by chlorinated ethenes.

South Carolina

Simulation of zones of contribution to wells at site GM–38, Naval Weapons Industrial Reserve Plant, Bethpage, New York

A three-dimensional groundwater-flow model is coupled with the particle-tracking program MODPATH to delineate zones of contribution to wells pumping from the Magothy aquifer and supplying water to a chlorinated volatile organic compound removal plant at site GM–38, Naval Weapons Industrial Reserve Plant, Bethpage, New York. By use of driller’s logs, a transitional probability approach generated three alternative realizations of heterogeneity within the Magothy aquifer to assess uncertainty in model representation. Finer-grained sediments with low hydraulic conductivity were realized as laterally discontinuous, thickening towards the south, and comprising about 17 percent of the total aquifer volume. Particle-tracking evaluations of a steady state present conditions model with alternative heterogeneity realizations were used to develop zones of contribution of remedial pumping wells. Because of heterogeneity and high rates of advection within the coarse-grained sediments, transport by dispersion and (or) diffusion was assumed to be negligible. Resulting zones of contribution of existing remedial wells are complex shapes, influenced by heterogeneity of each realization and other nearby hydrologic stresses. The use of two particle tracking techniques helped identify zones of contribution to wells. Backtracking techniques and observations of points of intersection of backward-tracked particles at shells of the GM–38 Hot Spot, as defined by surfaces of equal total volatile organic compound concentration, identified the source of water within the GM–38 Hot Spot to simulated wells. Forward-tracking techniques identified the fate of water within the GM–38 Hot Spot, including well capture and discharge to model constant head and drain boundaries. The percentage of backward-tracked particles, started at GM–38 wells that were sourced from within the Hot Spot, varied from 72.0 to 98.2, depending on the Hot Spot delineation used (present steady state model and Magothy aquifer heterogeneity realization A). The percentage of forward-tracked particles that were captured by GM–38 wells varied from 81.1 to 94.6, depending on the Hot Spot delineation used, with the remainder primarily captured by Bethpage Water District Plant 4 production wells (present steady state model and Magothy aquifer heterogeneity realization A). Less than 1 percent of forward-tracked particles ultimately discharge at model constant head and drain boundaries. The differences between forward- and backward-tracked particle percentage ranges are due to some forward-tracked particles not being captured by GM–38 wells, and some backward-tracked particles not intersecting specific regions of the Hot Spot. During 2013, an aquifer test generated detailed time series of well pumping rates and corresponding water-level responses were recorded at numerous locations. These data were used to verify the present conditions steady state model and demonstrate the sensitivity of model results to transient-state changes.

New York

Delineation of areas having elevated electrical conductivity, orientation and characterization of bedrock fractures, and occurrence of groundwater discharge to surface water at the U.S. Environmental Protection Agency Barite Hill/Nevada Goldfields Superfund site near McCormick, South Carolina

During October 2012 through March 2013, the U.S. Geological Survey (USGS), in cooperation with the U.S. Environmental Protection Agency (EPA) Region 4, Superfund Section, conducted borehole geophysical logging, surface geophysical surveys, and water-quality profiling in selected wells and areas to characterize or delineate the extent of elevated subsurface electrical conductivity at the EPA Barite Hill/Nevada Goldfields Superfund site near McCormick, South Carolina. Elevated electrical conductivity measured at the site may be related to native rock materials, waste rock disposal areas used in past operations, and (or) groundwater having elevated dissolved solids (primarily metals and major ions) related to waste migration. Five shallow screened wells and four open-borehole bedrock wells were logged by using a suite of borehole tools, and downhole water-quality profiles were recorded in two additional wells. Well depths ranged from about 26 to 300 feet below land surface. Surface geophysical surveys based on frequency-domain electromagnetic and distributed temperature sensing (DTS) techniques were used to identify areas of elevated electrical conductivity (Earth materials and groundwater) and potential high dissolved solids in groundwater and surface water on land and in areas along the northern unnamed tributary at the site. Results from the electromagnetic-induction logging of four selected wells near the Main Pit and one well located about 800 feet southeast of the Main Pit lake indicate that elevated electrical conductivity extends to a depth of about 110 feet below land surface. Groundwater-quality properties recorded in eight selected wells were highly variable, suggesting a broad spectrum of geochemical conditions and contaminant concentrations within the groundwater system. Ranges of field water-quality properties recorded from water-profiling of groundwater in all wells logged were as follows: pH, 3.1 to 9.2; specific conductance, 48 to 5,300 microsiemens per centimeter; dissolved oxygen, 0.2 to 4.4 milligrams per liter; and water temperature, 17.0 to 18.0 degrees Celsius. The highest specific conductance and lowest pH measurements were made in boreholes located between the Main Pit lake and the northern unnamed tributary. Conceptually, these wells may intercept elevated dissolved solids in groundwater leaking from the Main Pit lake along a flow path that discharges into the unnamed tributary to the north. Results from surface geophysical electromagnetic and fiber-optics surveys confirm areas of focused discharge of groundwater near the Main Pit lake along the northern unnamed tributary. The frequency-domain surface electromagnetic surveys also identified an area with higher levels of elevated electrical conductivity located northwest of the former Rainsford Pit area. Bedrock properties were characterized from borehole geophysical logs collected from three open-borehole bedrock wells. The mean strike azimuth of the borehole foliation data measured in bedrock well IR-1 was 221° (N. 41° E.), and the mean dip angle was 78° to the northwest. Dominant strike azimuth orientations of primary fractures measured in three boreholes were from 210° to 250° (N. 30° E. to N. 70° E.) with a mean dip of 68° northwest. Transmissivity estimates interpreted from the heat-pulse flowmeter data from bedrock well IR-1 were about 69 feet squared per day, and the radius of influence was estimated at about 640 feet.

South Carolina

The Phosphoria Formation at the Hot Springs Mine in Southeast Idaho: A source of selenium and other trace elements to surface water, ground water, vegetation, and biota

Major-element oxides and trace elements in the Phosphoria Formation at the Hot Springs Mine, Idaho were determined by a series of techniques. In this report, we examine the distribution of trace elements between the different solid components aluminosilicates, apatite, organic matter, opal, calcite, and dolomite that largely make up the rocks. High concentrations of several trace elements throughout the deposit, for example, As, Cd, Se, Tl, and U, at this and previously examined sites have raised concern about their introduction into the environment via weathering and the degree to which mining and the disposal of mined waste rock from this deposit might be accelerating that process. The question addressed here is how might the partitioning of trace elements between these solid host components influence the introduction of trace elements into ground water, surface water, and eventually biota, via weathering? In the case of Se, it is partitioned into components that are quite labile under the oxidizing conditions of subaerial weathering. As a result, it is widely distributed throughout the environment. Its concentration exceeds the level of concern for protection of wildlife at virtually every trophic level.

Idaho

Examples of deep-water-bottom multiple dereverberation techniques applied to seismic-reflection data from the Atlantic Outer Continental Shelf

Seismic-reflection data recorded in deep water over the Atlantic Outer Continental Shelf are often dominated by multiply-reflected seismic energy. This energy reverberates between the surface of the water and the seafloor (or other strong reflectors), and makes portions of the seismic data completely useless. Several different data-processing techniques can be applied to partially suppress these multiples and enhance the interpretability of the data. These techniques include (1) the three-point operator, (2) predictive deconvolution, (3) near-trace muting, (4) spatially variant bandpass filters, (5) Nth root stack, and (6) trace distance weighting. Application of these methods to several seismic lines indicates that trace distance weighting is the most useful method studied for suppressing deep-water-bottom multiples for data from the Atlantic Outer Continental Shelf.

Open-File Report

Summary of exploration geochemical and mineralogical studies at the giant Pebble porphyry Cu-Au-Mo deposit, Alaska: Implications for exploration under cover

Exploration geochemical and mineralogical studies by the U.S. Geological Survey at the Pebble porphyry Cu-Au-Mo deposit were designed to (1) determine whether the concealed deposit can be detected with surface samples, (2) better understand the processes of metal migration from the deposit to the surface, and (3) test existing methods for assessing concealed mineral resources and/or develop new ones. Surface water (ponds, streams, and springs), pond and stream sediment, soils subjected to various leaching techniques, and glacial till samples were collected. The tilted nature of the undisturbed orebody, varying depth of cover, and later glacial processes, strongly influence the geochemical responses and processes active on the various sample media. The multimedia approach aids in identifying possible processes that caused the significant geochemical variations within and among the various media. These processes include the following: In the Pebble West zone, thin cover and local exposure of the orebody have facilitated the oxidation of pyrite and other sulfides, and associated ferrous-ferric iron reactions, resulting in the local natural acidification of ponds observed in the West zone, and in associated metal anomalies in waters, sediments, and soils. In contrast, the East zone is concealed by both glacial deposits and underlying thick cover rocks, which precludes the oxidation of sulfides in the underlying orebody. Low-level geochemical anomalies in circumneutral spring and pond waters from the East zone are discernible only by using high resolution-inductively coupled plasma-mass spectrometry with lower limits of determination two and perhaps three orders of magnitude lower than traditional methods. A variety of partial leaches of soils over the East zone reveal geochemical anomalies in a similar suite of elements that may be related to upwelling waters from depth along graben-bounding faults. The indicator minerals gold, jarosite, and andradite in till reveal a displaced mineralogical anomaly to the west and south of the Pebble orebody, as ore-related minerals were scraped from the orebody and deposited in till downice of the deposit. Geochemical anomalies in pond water and sediment over the displaced till are attributed to the ore-related minerals in till. This orientation study demonstrates the strong control of local geologic and geochemical settings on the effectiveness of different traditional and newer reconnaissance geochemical exploration techniques and thus has important implications for exploration.

Alaska

Low-flow characteristics of streams under natural and diversion conditions, Waipiʻo Valley, Island of Hawaiʻi, Hawaiʻi

Over the past 100 years, natural streamflow in Waipiʻo Valley has been reduced by the transfer of water out of the valley by Upper and Lower Hāmākua Ditches. The physical condition and diversion practices along the two ditch systems have varied widely over the years, and as a result, so have their effects on natural streamflow in Waipiʻo Valley. Recent renovation and improvements to Lower Hāmākua Ditch system, along with proposals for its future operation and water-diversion strategies, have unknown implications. The purpose of this report is to quantify the availability of streamflow and to determine the effects of current and proposed diversion strategies on the low-flow hydrology in Waipiʻo Valley. In this report, the low-flow hydrology of Waipiʻo Valley is described in terms of flow-duration statistics. Flow-duration statistics were computed for three locations in the Waipiʻo Valley study area where long-term surface-water gaging stations have been operated. Using a variety of streamflow record-extension techniques, flow-duration statistics were estimated at an additional 13 locations where only few historical data are available or where discharge measurements were made as part of this study. Flow-duration statistics were computed to reflect natural conditions, current (2000-2005) diversion conditions, and proposed future diversion conditions at the 16 locations. At the downstream limit of the study area, on Wailoa Stream at an altitude of 190 feet, a baseline for evaluating the availability of streamflow is provided by computed flow-duration statistics that are representative of natural, no-diversion conditions. At the Wailoa gaging station, 95- and 50-percentile discharges under natural conditions were determined to be 86 and 112 cubic feet per second, respectively. Under 1965-1969 diversion conditions, natural 95- and 50-percentile discharges were reduced by 52 and 53 percent, to 41 and 53 cubic feet per second, respectively. Under current (2000-2005) diversion conditions, natural 95- and 50-percentile discharges were reduced by 21 and 24 percent, to 68 and 85 cubic feet per second, respectively. Under proposed future diversion conditions, natural 95- and 50-percentile discharges would be reduced by 33 and 24 percent, to 58 and 85 cubic feet per second, respectively. Compared to discharges that reflect current (2000-2005) diversion conditions, proposed future diversion conditions would reduce 95-percentile discharges, which are representative of moderate drought levels in the stream, by 15 percent. No change would be expected in 50-percentile discharges, which are representative of normal conditions. The effects of current (2000-2005) and proposed future diversion conditions on the natural flow of streams in the Waipiʻo Valley study area differ, depending on the location. Under current (2000-2005) diversion conditions, reductions in natural 95- or 50-percentile discharges of greater than 30 percent were found in Kawainui Stream downstream from Upper Hamakua Ditch to an altitude of about 1,435 feet and in the reach of Waimā Stream between Upper and Lower Hāmākua Ditches. Under proposed future diversion conditions, reductions in natural 95- or 50-percentile discharges of greater than 30 percent were found in Kawainui Stream downstream from Upper Hamakua Ditch to an altitude of about 1,435 feet, in the reach of Waimā Stream between Upper and Lower Hāmākua Ditches, and along most stream reaches downstream from Lower Hāmākua Ditch, except for Waimā Stream.

Hawai?i

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in northwestern Missouri, 2019 and 2020

Bathymetric data were collected at 12 water-supply lakes in northwestern Missouri by the U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources and in collaboration with various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. Ten of the lakes were surveyed from July to September 2019, one of the original 10 was resurveyed in March 2020, and two lakes of high interest near Maryville were surveyed in June 2020. Six of the lakes had been surveyed by the U.S. Geological Survey before, and the recent surveys were compared to the earlier surveys to document the changes in the bathymetric surface and capacity of the lake and to produce a bathymetric change map. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. At Hamilton Reservoir, a Global Navigation Satellite System survey receiver was used to collect additional bathymetric data at several points across four transects and around the perimeter of a substantial shallow area filled with aquatic vegetation upstream from a low-clearance bridge on the northern arm. Data points from the various sources were exported at a gridded data resolution appropriate to each lake. Data outside the multibeam echosounder survey extent and greater than the surveyed water-surface elevation generally were obtained from data collected using aerial light detection and ranging point cloud data, 1/9 arc-second National Elevation Dataset data based on aerial light detection and ranging data, or both. 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 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 data from a previous bathymetric survey existed at a given lake, a bathymetric change map was generated from the elevation difference between the previous survey and the 2019 bathymetric survey data points. After applying any vertical elevation changes to the previous survey data to ensure a match to the 2019 survey datum, coincident points between the surveys were found, and a bathymetric change map was generated using the coincident point data. A decrease in capacity was observed at all the lakes for which a previous survey existed. The decrease in capacity at the primary spillway or intake elevation ranged from 0.8 percent at Lake Viking to 21.4 percent at Middle Fork Grand River Reservoir. The mean bathymetric change ranged from 0.33 foot at Willow Brook Lake to 1.18 feet at Middle Fork Grand River Reservoir. The computed sedimentation rate generally ranged from 0.54 to 4.19 acre-feet per year at Maysville Lake and Middle Fork Grand River Reservoir, respectively; however, Lake Viking had the largest sedimentation rate of 14.9 acre-feet per year, despite having the smallest decrease in capacity at the spillway elevation of only 0.8 percent and a mean bathymetric change of only 0.4 foot. Evidence of dredging was observed in the bathymetric surface for Lake Viking. Some changes observed in some bathymetric change maps are hypothesized to result from the difference in data collection equipment and techniques between the previous and present bathymetric surveys. Certain erosional features around the perimeter of certain lakes may be the result of wave action during low-water years.

Missouri

Hydrologic description of the Braden River watershed, west-central Florida

The Braden River watershed drains an 83-square mile area in west-central Florida and is the largest tributary to the Manatee River. The hydrology of the Braden River was altered in 1936 when the city of Bradenton created Ward Lake, a reservoir with an 838-foot broad-crested weir 6 miles upstream from the mouth. In 1985 the reservoir, which is the sole source of drinking water for the city of Bradenton, was expanded and supplies an annual average of 5.7 million gallons of water per day. The Braden River can be hydrologically divided into three distinct sections that include an 8.6-mile reach of naturally incised, free-flowing channel; a 6.4-mile reach of impounded river created by the Ward Lake reservoir and weir; and a 6-mile reach of tidal estuary. Ten first-order and two second-order tributaries that flow into the Braden River were examined in this report. The Braden River watershed is dominated by low topographic relief. The two physiographic zones that contain the Braden River watershed, the Gulf Coast Lowlands and De Soto Plain, are both poorly drained and have numerous depressional features. The climate is subtropical with an annual average rainfall of 56 inches, annual average temperatures of 72 degrees Fahrenheit, and estimated annual lake evaporation of 52 inches. The soil series in the watershed are predominantly Myakka-Cassia and the EauGallie-Floridana; these series are characterized as nearly level and poorly drained soils. Land use within the watershed is the fastest changing characteristic that affects the hydrology of the system. The western half of the watershed is typically urban and includes parts of the city of Bradenton. Land use in the eastern half of the watershed is predominantly agricultural, but the explosive population growth of the area is driving the development of medium to high-density residential communities. The three major aquifers underlying the Braden River watershed are the surficial, intermediate, and Floridan aquifer systems. The surficial aquifer generally is underlain in places by a clay layer that enhances the ground-water flow of the surficial aquifer to surface-water bodies. The intermediate aquifer system has discontinuous water-bearing units, but retards ground-water movement between the surficial and Floridan aquifer system. The Floridan aquifer system consists of the Upper and Lower Floridan aquifers separated by a middle confining unit. The Upper Floridan aquifer is the primary source for ground-water withdrawals in the watershed and has, at times, heads 20 feet higher than land-surface elevation. Discharge over the Ward Lake weir into the tidal estuary was measured using volumetric and standard discharge measurement techniques. Annual mean flow for water years 1993 and 1994 were 59.7 and 57.3 cubic feet per second, respectively. Weir coefficients, calculated from discharge measurements, ranged from 0.023 to 2.99, depending on the head of water over the weir, and the method of determining length of flow on the weir. Weir coefficients calculated from the theoretical rating ranged from 0.032 to 3.11. No significant seepage was found around the ends of the weir, and no leakage was detected through the weir.

Florida