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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27Linked to original sources

USGS2011

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↗

USGS2012

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↗

USGS2018

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↗

USGS2019

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↗

USGS2020

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↗

USGS2020

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↗

USGS2024

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↗

USGS2026

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↗

USGS2008

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↗

USGS2006

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↗

USGS2008

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↗

USGS2006

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↗

USGS2013

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↗

USGS2017

Documentation for the “XT3D” option in the Node Property Flow (NPF) Package of MODFLOW 6

This report describes the “XT3D” option in the Node Property Flow (NPF) Package of MODFLOW 6. The XT3D option extends the capabilities of MODFLOW by enabling simulation of fully three-dimensional anisotropy on regular or irregular grids in a way that properly takes into account the full, three-dimensional conductivity tensor. It can also improve the accuracy of groundwater-flow simulations in cases in which the model grid violates certain geometric requirements. Three example problems demonstrate the use of the XT3D option to simulate groundwater flow on irregular grids and through three-dimensional porous media with anisotropic hydraulic conductivity. Conceptually, the XT3D method of estimating flow between two MODFLOW 6 model cells can be viewed in terms of three main mathematical steps: construction of head-gradient estimates by interpolation; construction of fluid-flux estimates by application of the full, three-dimensional form of Darcy’s Law, in which the conductivity tensor can be heterogeneous and anisotropic; and construction of the flow expression by enforcement of continuity of flow across the cell interface. The resulting XT3D flow expression, which relates the flow across the cell interface to the values of heads computed at neighboring nodes, is the sum of terms in which conductance-like coefficients multiply head differences, as in the conductance-based flow expression the NPF Package uses by default. However, the XT3D flow expression contains terms that involve “neighbors of neighbors” of the two cells for which the flow is being calculated. These additional terms have no analog in the conductance-based formulation. When assembled into matrix form, the XT3D formulation results in a larger stencil than the conductance-based formulation; that is, each row of the coefficient matrix generally contains more nonzero elements. The “RHS” suboption can be used to avoid expanding the stencil by placing the additional terms on the right-hand side of the matrix equation and evaluating them at the previous iteration or time step. The XT3D option can be an alternative to the Ghost-Node Correction (GNC) Package. However, the XT3D formulation is typically more computationally intensive than the conductance-based formulation the NPF Package uses by default, either with or without ghost nodes. Before deciding whether to use the GNC Package or XT3D option for production runs, the user should consider whether the conductance-based formulation alone can provide acceptable accuracy for the particular problem being solved.

Techniques and Methods↗

USGS2012

Approaches in highly parameterized inversion - GENIE, a general model-independent TCP/IP run manager

GENIE is a model-independent suite of programs that can be used to generally distribute, manage, and execute multiple model runs via the TCP/IP infrastructure. The suite consists of a file distribution interface, a run manage, a run executer, and a routine that can be compiled as part of a program and used to exchange model runs with the run manager. Because communication is via a standard protocol (TCP/IP), any computer connected to the Internet can serve in any of the capacities offered by this suite. Model independence is consistent with the existing template and instruction file protocols of the widely used PEST parameter estimation program. This report describes (1) the problem addressed; (2) the approach used by GENIE to queue, distribute, and retrieve model runs; and (3) user instructions, classes, and functions developed. It also includes (4) an example to illustrate the linking of GENIE with Parallel PEST using the interface routine.

Techniques and Methods↗

USGS1996

A coupled surface-water and ground-water flow model (MODBRANCH) for simulation of stream-aquifer interaction

Ground-water and surface-water flow models traditionally have been developed separately, with interaction between subsurface flow and streamflow either not simulated at all or accounted for by simple formulations. In areas with dynamic and hydraulically well-connected ground-water and surface-water systems, stream-aquifer interaction should be simulated using deterministic responses of both systems coupled at the stream-aquifer interface. Accordingly, a new coupled ground-water and surface-water model was developed by combining the U.S. Geological Survey models MODFLOW and BRANCH; the interfacing code is referred to as MODBRANCH. MODFLOW is the widely used modular three-dimensional, finite-difference ground-water model, and BRANCH is a one-dimensional numerical model commonly used to simulate unsteady flow in open- channel networks. MODFLOW was originally written with the River package, which calculates leakage between the aquifer and stream, assuming that the stream's stage remains constant during one model stress period. A simple streamflow routing model has been added to MODFLOW, but is limited to steady flow in rectangular, prismatic channels. To overcome these limitations, the BRANCH model, which simulates unsteady, nonuniform flow by solving the St. Venant equations, was restructured and incorporated into MODFLOW. Terms that describe leakage between stream and aquifer as a function of streambed conductance and differences in aquifer and stream stage were added to the continuity equation in BRANCH. Thus, leakage between the aquifer and stream can be calculated separately in each model, or leakages calculated in BRANCH can be used in MODFLOW. Total mass in the coupled models is accounted for and conserved. The BRANCH model calculates new stream stages for each time interval in a transient simulation based on upstream boundary conditions, stream properties, and initial estimates of aquifer heads. Next, aquifer heads are calculated in MODFLOW based on stream stages calculated by BRANCH, aquifer properties, and stresses. This process is repeated until convergence criteria are met for head and stage. Because time steps used in ground-water modeling can be much longer than time intervals used in surface- water simulations, provision has been made for handling multiple BRANCH time intervals within one MODFLOW time step. An option was also added to BRANCH to allow the simulation of channel drying and rewetting. Testing of the coupled model was verified by using data from previous studies; by comparing results with output from a simpler, four-point implicit, open-channel flow model linked with MODFLOW; and by comparison to field studies of L-31N canal in southern Florida.

Techniques of Water-Resources Investigations↗

USGS2002

Near-real-time simulation and internet-based delivery of forecast-flood inundation maps using two-dimensional hydraulic modeling--A pilot study for the Snoqualmie River, Washington

A system of numerical hydraulic modeling, geographic information system processing, and Internet map serving, supported by new data sources and application automation, was developed that generates inundation maps for forecast floods in near real time and makes them available through the Internet. Forecasts for flooding are generated by the National Weather Service (NWS) River Forecast Center (RFC); these forecasts are retrieved automatically by the system and prepared for input to a hydraulic model. The model, TrimR2D, is a new, robust, two-dimensional model capable of simulating wide varieties of discharge hydrographs and relatively long stream reaches. TrimR2D was calibrated for a 28-kilometer reach of the Snoqualmie River in Washington State, and is used to estimate flood extent, depth, arrival time, and peak time for the RFC forecast. The results of the model are processed automatically by a Geographic Information System (GIS) into maps of flood extent, depth, and arrival and peak times. These maps subsequently are processed into formats acceptable by an Internet map server (IMS). The IMS application is a user-friendly interface to access the maps over the Internet; it allows users to select what information they wish to see presented and allows the authors to define scale-dependent availability of map layers and their symbology (appearance of map features). For example, the IMS presents a background of a digital USGS 1:100,000-scale quadrangle at smaller scales, and automatically switches to an ortho-rectified aerial photograph (a digital photograph that has camera angle and tilt distortions removed) at larger scales so viewers can see ground features that help them identify their area of interest more effectively. For the user, the option exists to select either background at any scale. Similar options are provided for both the map creator and the viewer for the various flood maps. This combination of a robust model, emerging IMS software, and application interface programming should allow the technology developed in the pilot study to be applied to other river systems where NWS forecasts are provided routinely.

Water-Resources Investigations Report↗

USGS2003

Water movement through thick unsaturated zones overlying the central High Plains aquifer, southwestern Kansas, 2000-2001

The role of irrigation as a driving force for water and chemical movement to the central High Plains aquifer is uncertain because of the thick unsaturated zone overlying the aquifer. Water potentials and profiles of tritium, chloride, nitrate, and pesticide concentrations were used to evaluate water movement through thick unsaturated zones overlying the central High Plains aquifer at three sites in southwestern Kansas. One site was located in rangeland and two sites were located in areas dominated by irrigated agriculture. In 2000?2001, the depth to water at the rangeland site was 50 meters and the depth to water at the irrigated sites was about 45.4 meters. Irrigation at the study sites began in 1955?56. Measurements of matric potential and volumetric water content indicate wetter conditions existed in the deep unsaturated zone at the irrigated sites than at the rangeland site. Total water potentials in the unsaturated zone at the irrigated sites systematically decreased with depth to the water table, indicating a potential existed for downward water movement from the unsaturated zone to the water table at those sites. At the rangeland site, total water potentials in the deep unsaturated zone indicate small or no potential existed for downward water movement to the water table. Postbomb tritium was not detected below a depth of 1.9 meters in the unsaturated zone or in ground water at the rangeland site. In contrast, postbomb tritium was detected throughout most of the unsaturated zone and in ground water at both irrigated sites. These results indicate post-1953 water moved deeper in the unsaturated zone at the irrigated sites than at the rangeland site. The depth of the interface between prebomb and postbomb tritium and a tritium mass-balance method were used to estimate water fluxes in the unsaturated zone at each site. The average water fluxes at the rangeland site were 5.4 and 4.4 millimeters per year for the two methods, which are similar to the average water flux (5.1 millimeters per year) estimated using a chloride mass-balance method. Tritium profiles in the unsaturated zone at the irrigated sites were complicated by the presence of tritium-depleted intervals separating upper and lower zones containing postbomb tritium. If the interface between prebomb and postbomb tritium was at the top of the tritium-depleted interval and postbomb tritium detected beneath that interval was from the declining water table in the area, then the average water flux at the irrigated sites was estimated to be 21 to 54 millimeters per year. If postbomb tritium detected beneath the tritium-depleted interval was from bypass or preferential water movement through the local unsaturated zone instead of the declining water table, then the minimum water flux at the irrigated sites was estimated to be 106 to 116 millimeters per year. In either case, water fluxes at the irrigated sites were at least 4 to 12 times larger than the flux at the rangeland site, indicating irrigation was an important driving force for water movement through the unsaturated zone. The presence of postbomb tritium and large nitrate and total pesticide concentrations (24 milligrams per liter as nitrogen and 0.923 microgram per liter, respectively) in ground water at the irrigated sites indicates irrigation water also was an important driving force for chemical movement to the water table. The persistence of a downward hydraulic gradient from the deep unsaturated zone to the water table at the irrigated sites, in addition to large nitrate and atrazine concentrations in deep soil water (34 milligrams per liter as nitrogen and 0.79 microgram per liter, respectively), indicate that the deep unsaturated zone will be a source of nitrate and atrazine to the aquifer in the future.

Water-Resources Investigations Report↗