Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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

Volcanic aquifers of Hawaiʻi—Contributions to assessing groundwater availability on Kauaʻi, Oʻahu, and Maui

The volcanic aquifers of the Hawaiian Islands supply water to 1.46 million residents, diverse industries, and a large component of the U.S. military in the Pacific. Groundwater also supplies fresh water that supports ecosystems in streams and near the coast. Hawaii’s aquifers are remarkably productive given their small size, but the capacity of the islands to store fresh groundwater is limited because each island is surrounded by seawater, and salt water underlies much of the fresh groundwater. The amount of fresh groundwater available for human use from Hawai‘i’s volcanic aquifers is constrained by the consequences of groundwater withdrawal. Restrictions placed on these consequences can translate to limitations on groundwater availability. Changes in recharge resulting from changes in land cover or climate can alter the effect of withdrawals. This study uses numerical models of the volcanic aquifers of the islands of Kaua‘i, O‘ahu, and Maui to quantify the consequences of historical and plausible future withdrawals and changes in recharge. The study compares the results of model simulations of multiple scenarios of historical and projected future withdrawal and recharge. Results of the simulations using the groundwater models of the islands of Kaua‘i, O‘ahu, and Maui have implications for other islands in Hawai‘i. Since the first modern water well was drilled in Hawai‘i in 1879, total groundwater withdrawals on Kaua‘i, O‘ahu, and Maui have risen to nearly 400 million gallons per day. Model simulations indicate that these withdrawals have caused reductions in groundwater discharge to streams and springs, reductions in groundwater discharge to the ocean, changes in subsurface flow between sectors within an island, lowering of groundwater levels, and rise of the interface between fresh water and salt water in the aquifers. Future increases in withdrawals will increase the severity of the consequences. Changes in recharge can alter the effect of withdrawals—increases in recharge can offset the consequences of withdrawals, whereas decreases in recharge can exacerbate the effects of withdrawals. This study quantifies the consequences of withdrawals for past and plausible future circumstances. The models can be used to test other circumstances. Limits placed on the consequences of withdrawals—such as restrictions to protect stream or coastal ecosystems that rely on groundwater discharge and limitations on water-level decline and rise of the freshwater-saltwater interface to protect the productivity of existing wells—can translate to limits on groundwater availability from Hawai‘i’s volcanic aquifers. Setting acceptable limits to the consequences of groundwater withdrawal is also a critical part of assessing groundwater availability. Once these limits are set, numerical models can be used to quantify the amount of water that can be withdrawn within those limits and thereby inform management decisions that seek to balance the need to limit the consequences of groundwater withdrawals with the need to develop water for human use.

Hawaii↗

Effects of aquifer heterogeneity on ground-water flow and chloride concentrations in the Upper Floridan aquifer near and within an active pumping well field, west-central Florida

Chloride concentrations have been increasing over time in water from wells within and near the Eldridge-Wilde well field, near the coast in west-central Florida. Variable increases in chloride concentrations from well to well over time are the combined result of aquifer heterogeneity and ground-water pumping within the Upper Floridan aquifer. Deep mineralized water and saline water associated with the saltwater interface appear to move preferentially along flow zones of high transmissivity in response to ground-water withdrawals. The calcium-bicarbonate-type freshwater of the Upper Floridan aquifer within the study area is variably enriched with ions by mixing with introduced deep and saline ground water. The amount and variability of increases in chloride and sulfate concentrations at each well are related to well location, depth interval, and permeable intervals intercepted by the borehole. Zones of high transmissivity characterize the multilayered carbonate rocks of the Upper Floridan aquifer. Well-developed secondary porosity within the Tampa/Suwannee Limestones and the Avon Park Formation has created producing zones within the Upper Floridan aquifer. The highly transmissive sections of the Avon Park Formation generally are several orders of magnitude more permeable than the Tampa/Suwannee Limestones, but both are associated with increased ground-water flow. The Ocala Limestone is less permeable and is dominated by primary, intergranular porosity. Acoustic televiewer logging, caliper logs, and borehole flow logs (both electromagnetic and heat pulse) indicate that the Tampa/Suwannee Limestone units are dominated by porosity owing to dissolution between 200 and 300 feet below land surface, whereas the porosity of the Avon Park Formation is dominated by fractures that occur primarily from 600 to 750 feet below land surface and range in angle from horizontal to near vertical. Although the Ocala Limestone can act as a semiconfining unit between the Avon Park Formation and the Tampa/Suwannee Limestones, seismic-reflection data and photolinear analyses indicate that fractures and discontinuities in the Ocala Limestone are present within the southwestern part of the well field. It is possible that some fracture zones extend upward from the Avon Park Formation through the Ocala, Suwannee, and Tampa Limestones to land surface. These fractures may provide a more direct hydrologic connection between transmissive zones that are vertically separated by less permeable stratigraphic units. Ground water moves along permeable zones within the Upper Floridan aquifer in response to changes in head gradients as a result of pumping. Borehole geophysical measurements, including flow logs, specific conductance logs, and continuous monitoring of specific conductance at selected fixed depths, indicate that borehole specific conductance varies substantially with time and in response to pumping stresses. Ground-water mixing between hydrogeologic units likely occurs along highly transmissive zones and within boreholes of active production wells. Ground-water movement and water-quality changes were greatest along the most transmissive zones. Variable mixing of three water-type end members (freshwater, deepwater, and saltwater) occurs throughout the study area. Both deepwater and saltwater are likely sources for elevated chloride and sulfate concentrations in ground water. Mass-balance calculations of mixtures of the three end members indicate that deepwater is found throughout the aquifer units. Samples from wells within the southwestern part of the well field indicate that deepwater migrates into the shallow permeable units in the southwestern part of the well field. Deepwater contributes to elevated sulfate and chloride concentrations, which increase with depth and are elevated in wells less than 400 feet deep. The greatest increases in chloride concentrations over time are found in water from wells closest to the saltwater interface. Gro

Scientific Investigations Report↗

Hydrology of Northern Utah Valley, Utah County, Utah, 1975-2005

The ground-water resources of northern Utah Valley, Utah, were assessed during 2003-05 to describe and quantify components of the hydrologic system, determine a hydrologic budget for the basin-fill aquifer, and evaluate changes to the system relative to previous studies. Northern Utah Valley is a horst and graben structure with ground water occurring in both the mountain-block uplands surrounding the valley and in the unconsolidated basin-fill sediments. The principal aquifer in northern Utah Valley occurs in the unconsolidated basin-fill deposits where a deeper unconfined aquifer occurs near the mountain front and laterally grades into multiple confined aquifers near the center of the valley. Sources of water to the basin-fill aquifers occur predominantly as either infiltration of streamflow at or near the interface of the mountain front and valley or as subsurface inflow from the adjacent mountain blocks. Sources of water to the basin-fill aquifers were estimated to average 153,000 (+/- 31,500) acre-feet annually during 1975-2004 with subsurface inflow and infiltration of streamflow being the predominant sources. Discharge from the basin-fill aquifers occurs in the valley lowlands as flow to waterways, drains, ditches, springs, as diffuse seepage, and as discharge from flowing and pumping wells. Ground-water discharge from the basin-fill aquifers during 1975-2004 was estimated to average 166,700 (+/- 25,900) acre-feet/year where discharge to wells for consumptive use and discharge to waterways, drains, ditches, and springs were the principal sources. Measured water levels in wells in northern Utah Valley declined an average of 22 feet from 1981 to 2004. Water-level declines are consistent with a severe regional drought beginning in 1999 and continuing through 2004. Water samples were collected from 36 wells and springs throughout the study area along expected flowpaths. Water samples collected from 34 wells were analyzed for dissolved major ions, nutrients, and stable isotopes of hydrogen and oxygen. Water samples from all 36 wells were analyzed for dissolved-gas concentration including noble gases and tritium/helium-3. Within the basin fill, dissolved-solids concentration generally increases with distance along flowpaths from recharge areas, and shallower flowpaths tend to have higher concentrations than deeper flowpaths. Nitrate concentrations generally are at or below natural background levels. Dissolved-gas recharge temperature data support the conceptual model of the basin-fill aquifers and highlight complexities of recharge patterns in different parts of the valley. Dissolved-gas data indicate that the highest elevation recharge sources for the basin-fill aquifer are subsurface inflow derived from recharge in the adjacent mountain block between the mouths of American Fork and Provo Canyons. Apparent ground-water ages in the basin-fill aquifer, as calculated using tritium/helium-3 data, range from 2 to more than 50 years. The youngest waters in the valley occur near the mountain fronts with apparent ages generally increasing near the valley lowlands and discharge area around Utah Lake. Flowpaths are controlled by aquifer properties and the location of the predominant recharge sources, including subsurface inflow and recharge along the mountain front. Subsurface inflow is distributed over a larger area across the interface of the subsurface mountain block and basin-fill deposits. Subsurface inflow occurs at a depth deeper than that at which mountain-front recharge occurs. Recharge along the mountain front is often localized and focused over areas where streams and creeks enter the valley, and recharge is enhanced by the associated irrigation canals.

Utah↗

StreamStats in Oklahoma– Drainage-basin characteristics and peak-flow frequency statistics for ungaged streams

The USGS Streamflow Statistics (StreamStats) Program was created to make geographic information systems-based estimation of streamflow statistics easier, faster, and more consistent than previously used manual techniques. The StreamStats user interface is a map-based internet application that allows users to easily obtain streamflow statistics, basin characteristics, and other information for user-selected U.S. Geological Survey data-collection stations and ungaged sites of interest. The application relies on the data collected at U.S. Geological Survey streamflow-gaging stations, computer aided computations of drainage-basin characteristics, and published regression equations for several geographic regions comprising the United States. The StreamStats application interface allows the user to (1) obtain information on features in selected map layers, (2) delineate drainage basins for ungaged sites, (3) download drainage-basin polygons to a shapefile, (4) compute selected basin characteristics for delineated drainage basins, (5) estimate selected streamflow statistics for ungaged points on a stream, (6) print map views, (7) retrieve information for U.S. Geological Survey streamflow-gaging stations, and (8) get help on using StreamStats. StreamStats was designed for national application, with each state, territory, or group of states responsible for creating unique geospatial datasets and regression equations to compute selected streamflow statistics. With the cooperation of the Oklahoma Department of Transportation, StreamStats has been implemented for Oklahoma and is available at http://water.usgs.gov/osw/streamstats/. The Oklahoma StreamStats application covers 69 processed hydrologic units and most of the state of Oklahoma. Basin characteristics available for computation include contributing drainage area, contributing drainage area that is unregulated by Natural Resources Conservation Service floodwater retarding structures, mean-annual precipitation at the drainage-basin outlet for the period 1961-1990, 10-85 channel slope (slope between points located at 10 percent and 85 percent of the longest flow-path length upstream from the outlet), and percent impervious area. The Oklahoma StreamStats application interacts with the National Streamflow Statistics database, which contains the peak-flow regression equations in a previously published report. Fourteen peak-flow (flood) frequency statistics are available for computation in the Oklahoma StreamStats application. These statistics include the peak flow at 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for rural, unregulated streams; and the peak flow at 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for rural streams that are regulated by Natural Resources Conservation Service floodwater retarding structures. Basin characteristics and streamflow statistics cannot be computed for locations in playa basins (mostly in the Oklahoma Panhandle) and along main stems of the largest river systems in the state, namely the Arkansas, Canadian, Cimarron, Neosho, Red, and Verdigris Rivers, because parts of the drainage areas extend outside of the processed hydrologic units.

Oklahoma↗

Observed and forecast flood-inundation mapping application-A pilot study of an eleven-mile reach of the White River, Indianapolis, Indiana

Near-real-time and forecast flood-inundation mapping products resulted from a pilot study for an 11-mile reach of the White River in Indianapolis. The study was done by the U.S. Geological Survey (USGS), Indiana Silver Jackets hazard mitigation taskforce members, the National Weather Service (NWS), the Polis Center, and Indiana University, in cooperation with the City of Indianapolis, the Indianapolis Museum of Art, the Indiana Department of Homeland Security, and the Indiana Department of Natural Resources, Division of Water. The pilot project showed that it is technically feasible to create a flood-inundation map library by means of a two-dimensional hydraulic model, use a map from the library to quickly complete a moderately detailed local flood-loss estimate, and automatically run the hydraulic model during a flood event to provide the maps and flood-damage information through a Web graphical user interface. A library of static digital flood-inundation maps was created by means of a calibrated two-dimensional hydraulic model. Estimated water-surface elevations were developed for a range of river stages referenced to a USGS streamgage and NWS flood forecast point colocated within the study reach. These maps were made available through the Internet in several formats, including geographic information system, Keyhole Markup Language, and Portable Document Format. A flood-loss estimate was completed for part of the study reach by using one of the flood-inundation maps from the static library. The Federal Emergency Management Agency natural disaster-loss estimation program HAZUS-MH, in conjunction with local building information, was used to complete a level 2 analysis of flood-loss estimation. A Service-Oriented Architecture-based dynamic flood-inundation application was developed and was designed to start automatically during a flood, obtain near real-time and forecast data (from the colocated USGS streamgage and NWS flood forecast point within the study reach), run the two-dimensional hydraulic model, and produce flood-inundation maps. The application used local building data and depth-damage curves to estimate flood losses based on the maps, and it served inundation maps and flood-loss estimates through a Web-based graphical user interface.

Indiana↗

Collection, processing, and interpretation of ground-penetrating radar data to determine sediment thickness at selected locations in Deep Creek Lake, Garrett County, Maryland, 2007

The U.S. Geological Survey collected geophysical data in Deep Creek Lake in Garrett County, Maryland, between September 17 through October 4, 2007 to assist the Maryland Department of Natural Resources to better manage resources of the Lake. The objectives of the geophysical surveys were to provide estimates of sediment thickness in shallow areas around the Lake and to test the usefulness of three geophysical methods in this setting. Ground-penetrating radar (GPR), continuous seismic-reflection profiling (CSP), and continuous resistivity profiling (CRP) were attempted. Nearly 90 miles of GPR radar data and over 70 miles of CSP data were collected throughout the study area. During field deployment and testing, CRP was determined not to be practical and was not used on a large scale. Sediment accumulation generally could be observed in the radar profiles in the shallow coves. In some seismic profiles, a thin layer of sediment could be observed at the water bottom. The radar profiles appeared to be better than the seismic profiles for the determination of sediment thickness. Although only selected data profiles were processed, all data were archived for future interpretation. This investigation focused on selected regions of the study area, particularly in the coves where sediment accumulations were presumed to be thickest. GPR was the most useful tool for interpreting sediment thickness, especially in these shallow coves. The radar profiles were interpreted for two surfaces of interest-the water bottom, which was defined as the "2007 horizon," and the interface between Lake sediments and the original Lake bottom, which was defined as the "1925 horizon"-corresponding to the year the Lake was impounded. The ground-penetrating radar data were interpreted on the basis of characteristics of the reflectors. The sediments that had accumulated in the impounded Lake were characterized by laminated, parallel reflections, whereas the subsurface below the original Lake bottom was characterized by more discontinuous and chaotic reflections, often with diffractions indicating cobbles or boulders. The reflectors were picked manually along the water bottom and along the interface between the Lake sediments and the pre-Lake sediments. A simple graphic approach was used to convert traveltimes to depth through water and depth through saturated sediments using velocities of the soundwaves through the water and the saturated sediments. Nineteen cross sections were processed and interpreted in 9 coves around Deep Creek Lake, and the difference between the 2007 horizon and the 1925 horizon was examined. In most areas, GPR data indicate a layer of sediment between 1 and 7 feet thick. When multiple cross sections from a single cove were compared, the cross sections indicated that sediment thickness decreased toward the center of the Lake.

Maryland↗

Well network installation and hydrogeologic data collection, Assateague Island National Seashore, Worcester County, Maryland, 2010

The U.S. Geological Survey, as part of its Climate and Land Use Change Research and Development Program, is conducting a multi-year investigation to assess potential impacts on the natural resources of Assateague Island National Seashore, Maryland that may result from changes in the hydrologic system in response to projected sea-level rise. As part of this effort, 26 monitoring wells were installed in pairs along five east-west trending transects. Each of the five transects has between two and four pairs of wells, consisting of a shallow well and a deeper well. The shallow well typically was installed several feet below the water table—usually in freshwater about 10 feet below land surface (ft bls)—to measure water-level changes in the shallow groundwater system. The deeper well was installed below the anticipated depth to the freshwater-saltwater interface—usually in saltwater about 45 to 55 ft bls—for the purpose of borehole geophysical logging to characterize local differences in lithology and salinity and to monitor tidal influences on groundwater. Four of the 13 shallow wells and 5 of the 13 deeper wells were instrumented with water-level recorders that collected water-level data at 15-minute intervals from August 12 through September 28, 2010. Data collected from these instrumented wells were compared with tide data collected north of Assateague Island at the Ocean City Inlet tide gage, and precipitation data collected by National Park Service staff on Assateague Island. These data indicate that precipitation events coupled with changes in ambient sea level had the largest effect on groundwater levels in all monitoring wells near the Atlantic Ocean and Chincoteague and Sinepuxent Bays, whereas precipitation events alone had the greatest impact on shallow groundwater levels near the center of the island. Daily and bi-monthly tidal cycles appeared to have minimal influence on groundwater levels throughout the island and the water-level changes that were observed appeared to vary among well sites, indicating that changes in lithology and salinity also may affect the response of water levels in the shallow and deeper groundwater systems throughout the island. Borehole geophysical logs were collected at each of the 13 deeper wells along the 5 transects. Electromagnetic induction logs were collected to identify changes in lithology; determine the approximate location of the freshwater-saltwater interface; and characterize the distribution of fresh and brackish water in the shallow aquifer, and the geometry of the fresh groundwater lens beneath the island. Natural gamma logs were collected to provide information on the geologic framework of the island including the presence and thickness of finer-grained deposits found in the subsurface throughout the island during previous investigations. Results of this investigation show the need for collection of continuous water-level data in both the shallow and deeper parts of the flow system and electromagnetic induction and natural gamma geophysical logging data to better understand the response of this groundwater system to changes in precipitation and tidal forcing. Hydrologic data collected as part of this investigation will serve as the foundation for the development of numerical flow models to assess the potential effects of climate change on the coastal groundwater system of Assateague Island.

Maryl↗

Groundwater contributions to excessive algal growth in the East Fork Carson River, Carson Valley, west-central Nevada, 2010 and 2012

Excessive algal growth and low dissolved oxygen concentrations were observed during low streamflow conditions during summer months along a 5,800-foot reach of the East Fork Carson River in Carson Valley, west-central Nevada. Algal growth from nutrient enrichment of a stream reduces aquatic diversity, threatens fish ecology and stream health, and can be a recreational nuisance. In response to concerns that groundwater discharging to the 5,800-foot reach of the East Fork Carson River may be a source of nutrients to the stream, the U.S. Geological Survey, in cooperation with the Carson Water Subconservancy District and the Nevada Division of Environmental Protection, conducted studies during the summers of 2010 and 2012 to gain an improved understanding of the contributions of nutrients to the stream from groundwater, characterize algal conditions and algal effects on water quality, assess potential sources of nitrate in groundwater discharging to the stream, and evaluate nitrate reduction in groundwater from denitrification. A reconnaissance study in the summer of 2010 along the 5,800-foot study reach located a subreach with clear evidence of nutrient-rich groundwater discharging to the stream. At the subreach, nitrate plus nitrite (referred to hereafter as nitrate) concentrations in groundwater discharging to the stream were high (average 2.75 milligrams per liter as nitrogen) along the right bank. The stream at this location had the highest stream nitrate concentrations (average 0.056 milligrams per liter as nitrogen) compared to other locations upstream and downstream of the subreach. As a result, the 2012 study focused on a 405-foot subreach of the East Fork Carson River centered where results from the 2010 study found the highest stream and groundwater concentrations of nitrate, as well as the greatest observed contributions of groundwater discharge to the stream. Groundwater nutrient concentrations were much higher than stream nutrient concentrations during the summer of 2012 during low streamflow conditions at the 405-foot subreach of the East Fork Carson River. Average groundwater nitrate and orthophosphate concentrations along the right bank of the 405‑foot subreach were 9 and 12 times higher, respectively, than in the stream at this subreach. Groundwater discharge rates to the study reach based on different methods varied from 0.09 to 1.2 cubic feet per second per mile. Estimated groundwater discharge rates to the right bank of the study subreach were used to calculate groundwater nutrient load estimates to the subreach right bank, which were found to be low when compared to stream nutrient loads. Elevated algal biomass levels above nuisance thresholds were observed during the summers of 2010 and 2012. The study reach was characterized as mesotrophic-eutrophic during the 2010 study and eutrophic during the 2012 study. The presence of algae caused daily dissolved oxygen and pH fluctuations in the stream, resulting in exceedances of the State of Nevada water-quality standards owing to low dissolved oxygen concentrations and high pH concentrations, although the standards might not have been applicable during 2012 because of extremely low streamflow. The addition of nutrients to the stream from the constant supply in groundwater discharge sustains the growth of algae during low streamflow conditions. In the summer when streamflow is low or very low, nutrient-rich groundwater discharge enters the stream through the sediment-water interface at the streambed. Because the attached algae is thick and stream velocity is low, the nutrient-rich water pools at the sediment-water interface. Higher nutrient concentrations at the streambed create a favorable microenvironment for algae attached to the substrate to consume available nutrients from the groundwater before the groundwater mixes with overlying stream water. The source of nitrate in groundwater in this subreach is anthropogenic because nitrate concentrations are greater than background groundwater nitrate concentrations in Douglas County, high groundwater nitrate concentrations are only found at the right bank of the stream near a housing development, and organic wastewater compounds indicative of human-derived sources were also detected in groundwater wells on the right bank of the stream. Nitrogen and oxygen isotope concentrations of nitrate in shallow groundwater were used to determine the specific source of the nitrate, but the isotopic values indicated denitrification was occurring. Further investigation is needed to determine the specific anthropogenic source of the nitrate in the groundwater because the denitrification present in all samples obscures the original source of nitrogen.

Nevada↗

Calibration of Precipitation-Runoff Modeling System (PRMS) to simulate prefire and postfire hydrologic response in the upper Rio Hondo Basin, New Mexico

The Precipitation-Runoff Modeling System (PRMS) is widely used to simulate the effects of climate, topography, land cover, and soils on landscape-level hydrologic responses and streamflow. The U.S. Geological Survey (USGS), in cooperation with the New Mexico Department of Homeland Security and Emergency Management, developed procedures to apply the PRMS model to simulate the effects of fire on hydrologic responses. A PRMS model was built of the upper Rio Hondo Basin from the headwaters to approximately 19 miles downstream from the USGS streamgage Rio Hondo above Chavez Canyon near Hondo, New Mexico, by using 24 hydrologic response units (HRUs), or hydrologically similar subareas, from the National Hydrologic Model. A quasi-graphical user interface was created to easily query and analyze published PRMS sensitivity-analysis data. Simulation of mean daily streamflow was most sensitive to parameters related to snowmelt or infiltration throughout the upper Rio Hondo Basin. In the basin’s eastern and northern HRUs, flashiness and timing of streamflow were most sensitive to interflow; in many western-basin HRUs (higher elevations), flashiness of streamflow was most sensitive to soil moisture parameters, and timing of streamflow was most sensitive to infiltration and evapotranspiration parameters. The PRMS model was calibrated for the fire-affected North Fork Eagle Creek subwatershed by comparing modeled to observed daily streamflow for the nonfrozen (May through October) period for a prefire and postfire time period. The prefire model was calibrated for the period 2007–12 before the 2012 fire, and the postfire model was calibrated for a 2-year (2014–15) period after the fire. Model parameterization combined manual adjustment of 8 parameters on the basis of prior knowledge and automated adjustment of the most sensitive parameters by using the Let Us Calibrate interface. A gridded, daily precipitation dataset that captured the spatial heterogeneity across the study watershed was used as the precipitation input for calibration. Model performance was assessed as satisfactory by using standard statistical measures for prefire and postfire periods. The calibrated model was run by using data from a single precipitation gage to better represent the effect of localized, extreme storms on postfire hydrologic response. The calibrated models for prefire and postfire conditions simulated streamflows with greater consistency than the uncalibrated model for the corresponding (prefire or postfire) period of hydrographic record. The effect of fire on streamflow was found to be primarily a shift from streamflow dominated by base flow prior to fire to streamflow dominated by surface runoff after fire.

New Mexico↗

A distributed temperature sensing investigation of groundwater discharge to Haskell Lake, Lac du Flambeau Reservation, Wisconsin, July 27–August 1, 2016

Haskell Lake is a shallow, 89-acre drainage lake in the headwaters of the Squirrel River, on the Lac du Flambeau Reservation in northern Wisconsin. Historically, this lake was an important producer of wild rice for the Lac du Flambeau Band of Lake Superior Chippewa Indians (LDF Tribe); but, beginning in the late 1970s, the rice began to diminish and by the late 1990s, the lake no longer had harvestable stands. Restoring wild rice to Haskell Lake is a long-term priority for the LDF Tribe. A first step towards that effort is the cleanup of a petroleum-contamination plume in the shallow aquifer upgradient of the northern end of the lake. Knowledge of the downgradient extent of the plume and the locations where contaminated water is discharging to the lake is needed to inform cleanup efforts. A cooperative study between the U.S. Geological Survey and the LDF Tribe was initiated to characterize the distribution of groundwater discharge to Haskell Lake in the areas downgradient of the contamination plume. A fiber optic distributed temperature sensing system was used to monitor temperatures at the sediment-water interface for a 7-day period in July and August 2016. Challenges during the investigation included data storage and power supply limitations, maintenance of calibration baths, accurate location of the cable in space, cable placement in weeds and soft sediment, the confounding effects of solar radiation, and contamination of the data by multiple sources of instrument noise. The problem of instrument noise was overcome by solving the fiber optic distributed temperature sensing calibration equation for two parameters that describe temporal variation in the source laser and the photon detectors that observe the backscatter. Early morning temperatures, when the influence of solar radiation via direct warming of the sediment-water interface is minimized, were used to evaluate groundwater discharge, similar to other studies. The results indicate a persistent, horizontal variation in temperature of as much as 5.5 degrees Celsius across the study area, with cooler temperatures interpreted to indicate spatially discrete preferential groundwater discharge. Results of the study can be used to determine locations for collecting lakebed pore water samples to better define the extent of contamination discharging to the lake.

Wisconsin↗

Benthic vertical hydraulic gradients in Upper Klamath Lake, Oregon, 2017

Groundwater piezometers and lake stilling wells were deployed as paired sets at 10 locations in Upper Klamath Lake in south-central Oregon from May to October 2017 to measure hydraulic heads in and beneath the lake. Continuous water-level data from piezometers and stilling wells were then used to calculate the vertical hydraulic gradient (VHG) across the sediment-water interface to determine the direction and relative magnitude of the movement of water between the lake and underlying sediments. Over the study period, heads in lake-bed sediments closely tracked lake levels, both decreasing from spring into autumn. Instantaneous VHG was highly dynamic at all sites and exhibited high-frequency (less than 1 day to less than 1 hour) variations in magnitude and direction. Instantaneous and weekly mean VHG values often exceeded, but were commonly within, the range of measurement uncertainty (VHG less than +0.009 foot per foot [ft/ft] and greater than -0.009 ft/ft). 63 percent of instantaneous VHG values and 66 percent of weekly mean VHG values were within this range. Study period mean VHG was within measurement uncertainty at seven of the nine sites that had continuous water-level data, but two littoral sites (LC03 and LS01) had positive (upward) values greater than measurement uncertainty and are likely locations of vertical groundwater seepage. Data collected in this study provide new information about the hydraulic conditions at the sediment-water interface in UKL and demonstrate that sediment-groundwater exchange in UKL is spatially and temporally heterogeneous.

Oregon↗

Assessment of fresh groundwater discharge and saline surface-water intrusion at Operable Unit 2, North Chevalier Field Disposal Area (Site 11), Naval Air Station Pensacola, Florida, 2018–22

Site 11 is a former landfill at North Chevalier Field Disposal Area in Operable Unit 2 at Naval Air Station Pensacola, in northwest Florida. Site 11 is adjacent to Bayou Grande, a shallow, tidally influenced, saline estuary of the Pensacola Bay watershed. Federal and Florida regulators have expressed concern that contaminants detected in groundwater beneath the inland parts of Site 11 may discharge to Bayou Grande. In 2017, the Department of Defense, U.S. Navy, Naval Facilities Engineering Systems Command Southeast asked the U.S. Geological Survey to assess the occurrence of fresh groundwater discharge to Bayou Grande at Site 11 and to delineate to the extent practicable the location of groundwater discharge. Between 2018 and 2022, the U.S. Geological Survey used a multiple-lines-of-evidence approach that included a visual method and three physical methods based on the temperature difference between groundwater and surface water to assess groundwater discharge. One of the physically based methods also used the difference in specific conductance between fresh groundwater and brackish to saline surface water. Combined, the data indicate that fresh groundwater from across Site 11 discharges primarily along the shoreline of the northern and northeastern part of Site 11. The data also indicate that saline surface water from Bayou Grande intrudes tens of feet into the shallow aquifer beneath Site 11. The combined data indicate that the interface between fresh groundwater and saline surface water changes over space and time. Any new monitoring wells proposed for installation near the shoreline of Site 11 should include approaches to monitor the changes in the location of the freshwater/saltwater interface. Care would need to be taken to collect any groundwater samples at the correct season and tidal period to provide the highest probability of collecting a representative sample of Site 11 groundwater unaffected by saltwater.

Florida↗

Estimation of impounded sediment volume in the Similkameen River upstream of Enloe Dam, Okanogan County, Washington

The Enloe Dam was built in 1920 for hydropower generation and impounds a steep-banked, narrow reach of the Similkameen River in north-central Washington. During the subsequent century, sediment from the Similkameen River watershed, which includes historical mining operations, accumulated within Enloe Dam’s impoundment. Enloe Dam ceased hydropower production in 1958 and is currently under consideration for removal, but there are concerns that the remobilization of sediments may have harmful impacts on downstream water quality. To complement previously published analyses of heavy metal concentrations within sediments and assess the total volume that may be transported following dam removal, this report presents estimates of the volume of sediment trapped behind Enloe Dam which was measured in 2020. The volume of sediment was estimated by comparing a bathymetric survey, collected using an acoustic Doppler continuous profiler, and a survey of the bedrock-sediment interface, collected using a continuous resistivity profiler. The study area spanned a 2.6-kilometer reach of the Similkameen River upstream from Enloe Dam. The volume of impounded sediment was calculated in 2020 by subtracting the elevation of the bedrock-sediment interface measured using the continuous resistivity profiler from a bathymetric surface measured by the acoustic Doppler current profiler. In 2020, the estimated volume of sediment impounded by Enloe Dam was 2.17±0.04 million cubic meters (Mm 3 ) compared to 1.37 Mm 3 measured in 1971. This equates to a deposition rate of approximately 16,300 cubic meters per year from 1971 to 2020. Continuous resistivity profiles revealed that bedrock within Enloe Dam’s impoundment was deepest (about 20 m) in the southern, downstream end of the profiles and shallowest (less than 5 m) in the northern upstream end of the profiles.

Washington↗

Determination of the delta(15N/14N)of Ammonium (NH4+) in Water: RSIL Lab Code 2898

The purpose of the technique described by Reston Stable Isotope Laboratory (RSIL) lab code 2898 is to determine the N isotopic composition, delta(15N/14N), abbreviated as d15N, of ammonium (NH4+) in water (freshwater and saline water). The procedure involves converting dissolved NH4+ into NH3 gas by raising the pH of the sample to above 9 with MgO and subsequently trapping the gas quantitatively as (NH4)2SO4 on a glass fiber (GF) filter. The GF filter is saturated with NaHSO4 and pressure sealed between two gas-permeable polypropylene filters. The GF filter 'sandwich' floats on the surface of the water sample in a closed bottle. NH3 diffuses from the water through the polypropylene filter and reacts with NaHSO4, forming (NH4)2SO4 on the GF filter. The GF filter containing (NH4)2SO4 is dried and then combusted with a Carlo Erba NC 2500 elemental analyzer (EA), which is used to convert total nitrogen in a solid sample into N2 gas. The EA is connected to a continuous-flow isotope-ratio mass spectrometer (CF-IRMS), which determines the relative difference in ratios of the amounts of the stable isotopes of nitrogen (15N and 14N) of the product N2 gas and a reference N2 gas. The filters containing the samples are compressed in tin capsules and loaded into a Costech Zero-Blank Autosampler on the EA. Under computer control, samples then are dropped into a heated reaction tube that contains an oxidant, where combustion takes place in a He atmosphere containing an excess of O2 gas. To remove S-O gases produced from the NaHSO4, a plug of Ag-coated Cu wool is inserted at the bottom of the reaction tube. Combustion products are transported by a He carrier through a reduction furnace to remove excess O2, toconvert all nitrogen oxides to N2, and to remove any remaining S-O gases. The gases then pass through a drying tube to remove water. The gas-phase products, mainly N2 and a small amount of background CO2, are separated by a gas chromatograph (GC). The gas is then introduced into the IRMS through a Finnigan ConFlo II interface. The ConFlo II interface is used to introduce not only sample into the IRMS but also N2 reference gas and He for sample dilution. The flash combustion is quantitative, so no isotopic fractionation is involved. The IRMS is a Finnigan Delta V CF-IRMS with 10 cups and is capable of detecting ion beams with mass/charge (m/z) 28, 29, 30. The ion beams from N2 are as follows: m/z 28 = 14N14N, m/z 29 = 14N15N, and m/z 30 = 15N15N. The ion beam with m/z 30 also represents 14N16O, which may indicate contamination or incomplete reduction.

Techniques and Methods↗

Determination of the δ 15 N and δ 13 C of total nitrogen and carbon in solids; RSIL lab code 1832

The purpose of the Reston Stable Isotope Laboratory (RSIL) lab code 1832 is to determine the δ( 15 N/ 14 N), abbreviated as δ 15 N, and the δ( 13 C/ 12 C), abbreviated as δ 13 C, of total nitrogen and carbon in a solid sample. A Carlo Erba NC 2500 elemental analyzer (EA) is used to convert total nitrogen and carbon in a solid sample into N 2 and CO 2 gas. The EA is connected to a continuous flow isotope-ratio mass spectrometer (CF-IRMS), which determines the relative difference in stable nitrogen isotope-amount ratio ( 15 N/ 14 N) of the product N 2 gas and the relative difference in stable carbon isotope-amount ratio ( 13 C/ 12 C) of the product CO 2 gas. The combustion is quantitative; no isotopic fractionation is involved. Samples are placed in tin capsules and loaded into a Costech Zero Blank Autosampler on the EA. Under computer control, samples then are dropped into a heated reaction tube that contains an oxidant, where combustion takes place in a helium atmosphere containing an excess of oxygen gas. Combustion products are transported by a helium carrier through a reduction furnace to remove excess oxygen and to convert all nitrous oxides into N 2 and through a drying tube to remove water. The gas-phase products, mainly CO 2 and N 2 , are separated by a gas chromatograph. The gas is then introduced into the IRMS through a Finnigan MAT (now Thermo Scientific) ConFlo II interface. The Finnigan MAT ConFlo II interface is used for introducing not only sample into the IRMS but also N 2 and CO 2 reference gases and helium for sample dilution. The flash combustion is quantitative; no isotopic fractionation is involved. The IRMS is a Thermo Scientific Delta V CF-IRMS. It has a universal triple collector, two wide cups with a narrow cup in the middle; it is capable of measuring mass/charge ( m/z ) 28, 29, 30 or with a magnet current change 44, 45, 46, simultaneously. The ion beams from these m/z values are as follows: m/z 28 = N 2 = 14 N/ 14 N; m/z 29 = N 2 = 14 N/ 15 N primarily; m/z 30 = NO = 14 N/ 16 O primarily, which is a sign of contamination or incomplete reduction; m/z 44 = CO 2 = 12 C 16 O 16 O; m/z 45 = CO 2 = 13 C 16 O 16 O primarily; and m/z 46 = CO 2 = 12 C 16 O 18 O primarily.

Techniques and Methods↗

Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water

This report focuses on measuring the flow of water across the interface between surface water and ground water, rather than the hydrogeological or geochemical processes that occur at or near this interface. The methods, however, that use hydrogeological and geochemical evidence to quantify water fluxes are described herein. This material is presented as a guide for those who have to examine the interaction of surface water and ground water. The intent here is that both the overview of the many available methods and the in-depth presentation of specific methods will enable the reader to choose those study approaches that will best meet the requirements of the environments and processes they are investigating, as well as to recognize the merits of using more than one approach. This report is designed to make the reader aware of the breadth of approaches available for the study of the exchange between surface and ground water. To accomplish this, the report is divided into four chapters. Chapter 1 describes many well-documented approaches for defining the flow between surface and ground waters. Subsequent chapters provide an in-depth presentation of particular methods. Chapter 2 focuses on three of the most commonly used methods to either calculate or directly measure flow of water between surface-water bodies and the ground-water domain: (1) measurement of water levels in well networks in combination with measurement of water level in nearby surface water to determine water-level gradients and flow; (2) use of portable piezometers (wells) or hydraulic potentiomanometers to measure hydraulic gradients; and (3) use of seepage meters to measure flow directly. Chapter 3 focuses on describing the techniques involved in conducting water-tracer tests using fluorescent dyes, a method commonly used in the hydrogeologic investigation and characterization of karst aquifers, and in the study of water fluxes in karst terranes. Chapter 4 focuses on heat as a tracer in hydrological investigations of the near-surface environment.

Techniques and Methods↗

MODFLOW-2005, the U.S. Geological Survey modular ground-water model - documentation of shared node local grid refinement (LGR) and the boundary flow and head (BFH) package

This report documents the addition of shared node Local Grid Refinement (LGR) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference ground-water flow model. LGR provides the capability to simulate ground-water flow using one block-shaped higher-resolution local grid (a child model) within a coarser-grid parent model. LGR accomplishes this by iteratively coupling two separate MODFLOW-2005 models such that heads and fluxes are balanced across the shared interfacing boundary. LGR can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined ground-water systems. Traditional one-way coupled telescopic mesh refinement (TMR) methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled shared-node method of LGR provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR, evaluates LGR accuracy and performance for two- and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR.

Techniques and Methods↗

MODFLOW–LGR—Documentation of ghost node local grid refinement (LGR2) for multiple areas and the boundary flow and head (BFH2) package

This report documents the addition of ghost node Local Grid Refinement (LGR2) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference groundwater flow model. LGR2 provides the capability to simulate groundwater flow using multiple block-shaped higher-resolution local grids (a child model) within a coarser-grid parent model. LGR2 accomplishes this by iteratively coupling separate MODFLOW-2005 models such that heads and fluxes are balanced across the grid-refinement interface boundary. LGR2 can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined groundwater systems. Traditional one-way coupled telescopic mesh refinement methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled ghost-node method of LGR2 provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR2, evaluates accuracy and performance for two-and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH2) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR2.

Techniques and Methods↗