Search USGS⌕ Search

Geology topics

Hedeff I. Essaid

Publications and source records attributed to Hedeff I. Essaid.

At least 19 recordsLinked to original sources

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research↗

Simulation of groundwater-flow dynamics in the U.S. Northern High Plains driven by multi-model estimates of surficial aquifer recharge

There is growing interest in incorporating higher-resolution groundwater modeling within the framework of large-scale land surface models (LSMs), including processes such as three-dimensional flow, variable soil saturation, and surface water/groundwater interactions. Conversely, complex groundwater models (e.g., the U.S. Geological Survey Groundwater-Flow Model, MODFLOW) often use simpler representations of land surface dynamics (e.g., surface vegetation, evapotranspiration, recharge) and may benefit from higher process fidelity and temporal resolutions in these inputs. This study investigates the potential of improving groundwater representation in LSMs and land surface dynamics in MODFLOW through forcing MODFLOW with recharge from LSMs. Groundwater simulations build on an existing and well-calibrated MODFLOW model of the U.S. Northern High Plains aquifer, a hydrologically complex basin under the dual impacts of conversion of native vegetation to intense irrigated agricultural fields and climate change. Simulated groundwater recharge from four different land models are used to drive MODFLOW groundwater simulations. Results show relatively large discrepancies between recharge estimates among simulations. Forcing MODFLOW using recharge simulated by some of the LSMs in place of a simple water balance model marginally improves MODFLOW groundwater simulation. Further, our results support the efficacy of coupling LSMs to a sophisticated groundwater model such as MODFLOW. The coupling results in notable improvements in matching the historical groundwater levels through reduction of the skewness coefficient in percent bias histogram (from 1.50 and 1.41 in original LSMs to 0.44 and 0.27, respectively, when MODFLOW is forced by groundwater recharge from LSMs) and reduction of bias. This modeling effort seeks to identify the best compromise between comprehensive land surface processes from global LSMs and advanced representation of groundwater from regional models.

Journal of Hydrology↗

Benchmarking high-resolution hydrologic model performance of long-term retrospective streamflow simulations in the contiguous United States

Because use of high-resolution hydrologic models is becoming more widespread and estimates are made over large domains, there is a pressing need for systematic evaluation of their performance. Most evaluation efforts to date have focused on smaller basins that have been relatively undisturbed by human activity, but there is also a need to benchmark model performance more comprehensively, including basins impacted by human activities. This study benchmarks the long-term performance of two process-oriented, high-resolution, continental-scale hydrologic models that have been developed to assess water availability and risks in the United States (US): the National Water Model v2.1 application of WRF-Hydro (NWMv2.1) and the National Hydrologic Model v1.0 application of the Precipitation–Runoff Modeling System (NHMv1.0). The evaluation is performed on 5390 streamflow gages from 1983 to 2016 ( ∼ 33 years) at a daily time step, including both natural and human-impacted catchments, representing one of the most comprehensive evaluations over the contiguous US. Using the Kling–Gupta efficiency as the main evaluation metric, the models are compared against a climatological benchmark that accounts for seasonality. Overall, the model applications show similar performance, with better performance in minimally disturbed basins than in those impacted by human activities. Relative regional differences are also similar: the best performance is found in the Northeast, followed by the Southeast, and generally worse performance is found in the Central and West areas. For both models, about 80 % of the sites exceed the seasonal climatological benchmark. Basins that do not exceed the climatological benchmark are further scrutinized to provide model diagnostics for each application. Using the underperforming subset, both models tend to overestimate streamflow volumes in the West, which could be attributed to not accounting for human activities, such as active management. Both models underestimate flow variability, especially the highest flows; this was more pronounced for NHMv1.0. Low flows tended to be overestimated by NWMv2.1, whereas there were both over and underestimations for NHMv1.0, but they were less severe. Although this study focused on model diagnostics for underperforming sites based on the seasonal climatological benchmark, metrics for all sites for both model applications are openly available online.

Hydrology and Earth System Sciences↗

Evaluating the dynamics of groundwater, lakebed transport, nutrient inflow and algal blooms in Upper Klamath Lake, Oregon, USA

Transport of nutrients to lakes can occur via surface-water inflow, atmospheric deposition, groundwater (GW) inflow and benthic processes. Identifying and quantifying within-lake nutrient sources and recycling processes is challenging. Prior studies in hypereutrophic Upper Klamath Lake, Oregon, USA, indicated that ~60% of the early summer phosphorus (P) load to the lake was internal and hypothesized to be lakebed sediment release. Dynamic nutrient transport processes were examined to better characterize the nutrient sources. One-dimensional heat transport models calibrated to observed lakebed temperatures and a cross-sectional GW flow model provided estimates of GW-inflow rates that were greatest in spring and decreased through summer. One-dimensional solute transport models calibrated to observed lakebed pore-water dissolved silica (Si) and dissolved phosphate-phosphorus (DP) concentrations indicated that nutrients were transported from the lakebed by advection, diffusion, and enhanced mixing by benthic organisms and waves, and that DP removal occurred near the lakebed interface. Estimated water, Si, DP and total-phosphorus (TP) budgets indicated that GW contributed 21% of lake water inflow and at least 26, 20 and 16% of total Si, DP and TP inflow, respectively, when conservatively assuming background GW nutrient concentrations. However, lakebed GW (LGW) is enriched in nutrients during flow through lakebed sediment and the estimated GW contribution increased to 29 (33), 49 (67) and 43% (61%) of total Si, DP and TP inflow, respectively, if 20% (50%) of GW inflow to the lake was assumed to have LGW concentrations. Net nutrient inflow to the lake was greatest in spring and coincident with the annual diatom bloom. Inflowing dissolved nutrients appear to be assimilated by diatoms during the spring and become available for the summer Aphanizomenon flos-aquae bloom when the diatoms senesce. Thus, nutrient-enriched GW inflow and nutrient recycling by successive algal blooms must be considered when evaluating internal nutrient loading to lakes.

Oregon↗

Evaluating the impact of irrigation on surface water – groundwater interaction and stream temperature in an agricultural watershed

Changes in groundwater discharge to streams caused by irrigation practices can influence stream temperature. Observations along two currently flood-irrigated reaches in the 640-square-kilometer upper Smith River watershed, an important agricultural and recreational fishing area in west-central Montana, showed a downstream temperature decrease resulting from groundwater discharge to the stream. A watershed-scale coupled surface water and groundwater flow model was used to examine changes in streamflow, groundwater discharge to the stream and stream temperature resulting from irrigation practices. The upper Smith River watershed was used to develop the model framework including watershed climate, topography, hydrography, vegetation, soil properties and current irrigation practices. Model results were used to compare watershed streamflow, groundwater recharge, and groundwater discharge to the stream for three scenarios: natural, pre-irrigation conditions (PreIrr); current irrigation practices involving mainly stream diversion for flood and sprinkler irrigation (IrrCurrent); and a hypothetical scenario with only groundwater supplying sprinkler irrigation (IrrGW). Irrigation increased groundwater recharge relative to natural PreIrr conditions because not all applied water was removed by crop evapotranspiration. Groundwater storage and groundwater discharge to the stream increased relative to natural PreIrr conditions when the source of irrigation water was mainly stream diversion as in the IrrCurrent scenario. The hypothetical IrrGW scenario, in which groundwater withdrawals were the sole source of irrigation water, resulted in widespread lowering of the water table and associated decreases in groundwater storage and groundwater discharge to the stream. A mixing analysis using model predicted groundwater discharge along the reaches suggests that stream diversion and flood irrigation, represented in the IrrCurrent scenario, has led to cooling of stream temperatures relative to natural PreIrr conditions improving fish thermal habitat. However, the decrease in groundwater discharge in the IrrGW scenario resulting from large-scale groundwater withdrawal for irrigation led to warmer than natural stream temperatures and possible degradation of fish habitat.

Montana↗

Contrasting nitrogen fate in watersheds using agricultural and water quality information

Surplus nitrogen (N) estimates, principal component analysis (PCA), and end-member mixing analysis (EMMA) were used in a multisite comparison contrasting the fate of N in diverse agricultural watersheds. We applied PCA-EMMA in 10 watersheds located in Indiana, Iowa, Maryland, Nebraska, Mississippi, and Washington ranging in size from 5 to 1254 km 2 with four nested watersheds. Watershed Surplus N was determined by subtracting estimates of crop uptake and volatilization from estimates of N input from atmospheric deposition, plant fixation, fertilizer, and manure for the period from 1987 to 2004. Watershed average Surplus N ranged from 11 to 52 kg N ha −1 and from 9 to 32% of N input. Solute concentrations in streams, overland runoff, tile drainage, groundwater (GW), streambeds, and the unsaturated zone were used in the PCA-EMMA procedure to identify independent components contributing to observed stream concentration variability and the end-members contributing to streamflow and NO 3 load. End-members included dilute runoff, agricultural runoff, benthic-processing, tile drainage, and oxic and anoxic GW. Surplus N was larger in watersheds with more permeable soils (Washington, Nebraska, and Maryland) that allowed greater infiltration, and oxic GW was the primary source of NO 3 load. Subsurface transport of NO 3 in these watersheds resulted in some removal of Surplus N by denitrification. In less permeable watersheds (Iowa, Indiana, and Mississippi), NO 3 was rapidly transported to the stream by tile drainage and runoff with little removal. Evidence of streambed removal of NO 3 by benthic diatoms was observed in the larger watersheds.

Indiana, Iowa, Maryland, Nebraska, Mississippi, Wa↗

The role of dynamic surface water-groundwater exchange on streambed denitrification in a first-order, low-relief agricultural watershed

The role of temporally varying surface water-groundwater (SW-GW) exchange on nitrate removal by streambed denitrification was examined along a reach of Leary Weber Ditch (LWD), Indiana, a small, first-order, low-relief agricultural watershed within the Upper Mississippi River basin, using data collected in 2004 and 2005. Stream stage, GW heads (H), and temperatures (T) were continuously monitored in streambed piezometers and stream bank wells for two transects across LWD accompanied by synoptic measurements of stream stage, H, T, and nitrate (NO 3 ) concentrations along the reach. The H and T data were used to develop and calibrate vertical two-dimensional, models of streambed water flow and heat transport across and along the axis of the stream. Model-estimated SW-GW exchange varied seasonally and in response to high-streamflow events due to dynamic interactions between SW stage and GW H. Comparison of 2004 and 2005 conditions showed that small changes in precipitation amount and intensity, evapotranspiration, and/or nearby GW levels within a low-relief watershed can readily impact SW-GW interactions. The calibrated LWD flow models and observed stream and streambed NO 3 concentrations were used to predict temporal variations in streambed NO 3 removal in response to dynamic SW-GW exchange. NO 3 removal rates underwent slow seasonal changes, but also underwent rapid changes in response to high-flow events. These findings suggest that increased temporal variability of SW-GW exchange in low-order, low-relief watersheds may be a factor contributing their more efficient removal of NO 3 .

Indiana↗

Organic contaminant transport and fate in the subsurface: evolution of knowledge and understanding

Toxic organic contaminants may enter the subsurface as slightly soluble and volatile nonaqueous phase liquids (NAPLs) or as dissolved solutes resulting in contaminant plumes emanating from the source zone. A large body of research published in Water Resources Research has been devoted to characterizing and understanding processes controlling the transport and fate of these organic contaminants and the effectiveness of natural attenuation, bioremediation, and other remedial technologies. These contributions include studies of NAPL flow, entrapment, and interphase mass transfer that have advanced from the analysis of simple systems with uniform properties and equilibrium contaminant phase partitioning to complex systems with pore-scale and macroscale heterogeneity and rate-limited interphase mass transfer. Understanding of the fate of dissolved organic plumes has advanced from when biodegradation was thought to require oxygen to recognition of the importance of anaerobic biodegradation, multiple redox zones, microbial enzyme kinetics, and mixing of organic contaminants and electron acceptors at plume fringes. Challenges remain in understanding the impacts of physical, chemical, biological, and hydrogeological heterogeneity, pore-scale interactions, and mixing on the fate of organic contaminants. Further effort is needed to successfully incorporate these processes into field-scale predictions of transport and fate. Regulations have greatly reduced the frequency of new point-source contamination problems; however, remediation at many legacy plumes remains challenging. A number of fields of current relevance are benefiting from research advances from point-source contaminant research. These include geologic carbon sequestration, nonpoint-source contamination, aquifer storage and recovery, the fate of contaminants from oil and gas development, and enhanced bioremediation.

Water Resources Research↗

Watershed-scale modeling of streamflow change in incised montane meadows

Land use practices have caused stream channel incision and water table decline in many montane meadows of the Western United States. Incision changes the magnitude and timing of streamflow in water supply source watersheds, a concern to resource managers and downstream water users. The hydrology of montane meadows under natural and incised conditions was investigated using watershed simulation for a range of hydrologic conditions. The results illustrate the interdependence between: watershed and meadow hydrology; bedrock and meadow aquifers; and surface and groundwater flow through the meadow for the modeled scenarios. During the wet season, stream incision resulted in less overland flow and interflow and more meadow recharge causing a net decrease in streamflow and increase in groundwater storage relative to natural meadow conditions. During the dry season, incision resulted in less meadow evapotranspiration and more groundwater discharge to the stream causing a net increase in streamflow and a decrease in groundwater storage relative to natural meadow conditions. In general, for a given meadow setting, the magnitude of change in summer streamflow and long-term change in watershed groundwater storage due to incision will depend on the combined effect of: reduced evapotranspiration in the eroded meadow; induced groundwater recharge; replenishment of dry season groundwater storage depletion in meadow and bedrock aquifers by precipitation during wet years; and groundwater storage depletion that is not replenished by precipitation during wet years.

Water Resources Research↗

Design of Cycle 3 of the National Water-Quality Assessment Program, 2013-23: Part 2: Science plan for improved water-quality information and management

This report presents a science strategy for the third decade of the National Water-Quality Assessment (NAWQA) Program, which since 1991, has been responsible for providing nationally consistent information on the quality of the Nation's streams and groundwater; how water quality is changing over time; and the major natural and human factors that affect current water quality conditions and trends. The strategy is based on an extensive evaluation of the accomplishments of NAWQA over its first two decades, the current status of water-quality monitoring activities by USGS and its partners, and an updated analysis of stakeholder priorities. The plan is designed to address priority issues and national needs identified by NAWQA stakeholders and the National Research Council (2012) irrespective of budget constraints. This plan describes four major goals for the third decade (Cycle 3), the approaches for monitoring, modeling, and scientific studies, key partnerships required to achieve these goals, and products and outcomes that will result from planned assessment activities. The science plan for 2013–2023 is a comprehensive approach to meet stakeholder priorities for: (1) rebuilding NAWQA monitoring networks for streams, rivers, and groundwater, and (2) upgrading models used to extrapolate and forecast changes in water-quality and stream ecosystem condition in response to changing climate and land use. The Cycle 3 plan continues approaches that have been central to the Program’s long-term success, but adjusts monitoring intensities and study designs to address critical information needs and identified data gaps. Restoration of diminished monitoring networks and new directions in modeling and interpretative studies address growing and evolving public and stakeholder needs for water-quality information and improved management, particularly in the face of increasing challenges related to population growth, increasing demands for water, and changing land use and climate. However, a combination of funding growth and extensive collaboration with other USGS programs and other Federal, State, and local agencies, public interest groups, professional and trade associations, academia, and private industry will be needed to fully realize the monitoring and modeling goals laid out in this plan (USGS Fact Sheet 2013-3008).

Open-File Report↗

VS2DI: Model use, calibration, and validation

VS2DI is a software package for simulating water, solute, and heat transport through soils or other porous media under conditions of variable saturation. The package contains a graphical preprocessor for constructing simulations, a postprocessor for displaying simulation results, and numerical models that solve for flow and solute transport (VS2DT) and flow and heat transport (VS2DH). Flow is described by the Richards equation, and solute and heat transport are described by advection-dispersion equations; the finite-difference method is used to solve these equations. Problems can be simulated in one, two, or three (assuming radial symmetry) dimensions. This article provides an overview of calibration techniques that have been used with VS2DI; included is a detailed description of calibration procedures used in simulating the interaction between groundwater and a stream fed by drainage from agricultural fields in central Indiana. Brief descriptions of VS2DI and the various types of problems that have been addressed with the software package are also presented.

Transactions of the ASABE↗

Simulation of aerobic and anaerobic biodegradation processes at a crude oil spill site

A two-dimensional, multispecies reactive solute transport model with sequential aerobic and anaerobic degradation processes was developed and tested. The model was used to study the field-scale solute transport and degradation processes at the Bemidji, Minnesota, crude oil spill site. The simulations included the biodegradation of volatile and nonvolatile fractions of dissolved organic carbon by aerobic processes, manganese and iron reduction, and methanogenesis. Model parameter estimates were constrained by published Monod kinetic parameters, theoretical yield estimates, and field biomass measurements. Despite the considerable uncertainty in the model parameter estimates, results of simulations reproduced the general features of the observed groundwater plume and the measured bacterial concentrations. In the simulation, 46% of the total dissolved organic carbon (TDOC) introduced into the aquifer was degraded. Aerobic degradation accounted for 40% of the TDOC degraded. Anaerobic processes accounted for the remaining 60% of degradation of TDOC: 5% by Mn reduction, 19% by Fe reduction, and 36% by methanogenesis. Thus anaerobic processes account for more than half of the removal of DOC at this site.

Water Resources Research↗

A multilayered sharp interface model of coupled freshwater and saltwater flow in coastal systems: Model development and application

A quasi three-dimensional, finite difference model, that simulates freshwater and saltwater flow separated by a sharp interface, has been developed to study layered coastal aquifer systems. The model allows for regional simulation of coastal groundwater conditions, including the effects of saltwater dynamics on the freshwater system. Vertically integrated freshwater and saltwater flow equations incorporating the interface boundary condition are solved within each aquifer. Leakage through confining layers is calculated by Darcy's law, accounting for density differences across the layer. The locations of the interface tip and toe, within grid blocks, are tracked by linearly extrapolating the position of the interface. The model has been verified using available analytical solutions and experimental results. Application of the model to the Soquel-Aptos basin, Santa Cruz County, California, illustrates the use of the quasi three-dimensional, sharp interface approach for the examination of freshwater-saltwater dynamics in regional systems. Simulation suggests that the interface, today, is still responding to long-term Pleistocene sea level fluctuations and has not achieved equilibrium with present day sea level conditions.

Water Resources Research↗

Effects of upstream dams versus groundwater pumping on stream temperature under varying climate conditions

The relative impact of a large upstream dam versus in‐reach groundwater pumping on stream temperatures was analyzed for humid, semiarid, and arid conditions with long dry seasons to represent typical climate regions where large dams are present, such as the western United States or eastern Australia. Stream temperatures were simulated using the CE‐QUAL‐W2 water quality model over a 110 km model grid, with the presence or absence of a dam at the top of the reach and pumping in the lower 60 km of the reach. Measured meteorological data from three representative locations were used as model input to simulate the impact of varying climate conditions on streamflow and stream temperature. For each climate condition four hypothetical streamflow scenarios were modeled: (1) natural (no dam or pumping), (2) large upstream dam present, (3) dam with in‐reach pumping, and (4) no dam with pumping, resulting in 12 cases. Dam removal, in the presence or absence of pumping, resulted in significant changes in stream temperature throughout the year for all three climate conditions. From March to August, the presence of a dam caused monthly mean stream temperatures to decrease on average by approximately 3.0°C, 2.5°C, and 2.0°C for the humid, semiarid, and arid conditions, respectively; however, stream temperatures generally increased from September to February. Pumping caused stream temperatures to warm in summer and cool in winter by generally less than 0.5°C because of a smaller pumping‐induced alteration in streamflow relative to the dam. Though the presence or absence of a large dam led to greater changes in stream temperature than the presence or absence of pumping, ephemeral conditions were increased both temporally and spatially because of pumping.

Water Resources Research↗

Design of Cycle 3 of the National Water-Quality Assessment Program, 2013-2022: Part 1: Framework of Water-Quality Issues and Potential Approaches

In 1991, the U.S. Congress established the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) Program to develop long-term, nationally consistent information on the quality of the Nation's streams and groundwater. Congress recognized the critical need for this information to support scientifically sound management, regulatory, and policy decisions concerning the increasingly stressed water resources of the Nation. The long-term goals of NAWQA are to: (1) assess the status of water-quality conditions in the United States, (2) evaluate long-term trends in water-quality conditions, and (3) link status and trends with an understanding of the natural and human factors that affect water quality. These goals are national in scale, include both surface water and groundwater, and include consideration of water quality in relation to both human uses and aquatic ecosystems. Since 1991, NAWQA assessments and findings have fostered and supported major improvements in the availability and use of unbiased scientific information for decisionmaking, resource management, and planning at all levels of government. These improvements have enabled agencies and stakeholders to cost-effectively address a wide range of water-quality issues related to natural and human influences on the quality of water and potential effects on aquatic ecosystems and human health (http://water.usgs.gov/nawqa/xrel.pdf). NAWQA, like all USGS programs, provides policy relevant information that serves as a scientific basis for decisionmaking related to resource management, protection, and restoration. The information is freely available to all levels of government, nongovernmental organizations, industry, academia, and the public, and is readily accessible on the NAWQA Web site and other diverse formats to serve the needs of the water-resource community at different technical levels. Water-quality conditions in streams and groundwater are described in more than 1,700 publications (available online at http://water.usgs.gov/nawqa/bib/), and are documented by more than 14 million data records representing about 7,600 stream sites, 8,100 wells, and 2,000 water-quality and ecological constituents that are available from the NAWQA data warehouse (http://infotrek.er.usgs.gov/traverse/f?p=NAWQA:HOME:0). The Program promotes collaboration and liaison with government officials, resource managers, industry representatives, and other stakeholders to increase the utility and relevance of NAWQA science to decisionmakers. As part of this effort, NAWQA supports integration of data from other organizations into NAWQA assessments, where appropriate and cost-effective, so that more comprehensive findings are available across geographic and temporal scales.

Open-File Report↗

Crude oil at the Bemidji Site: 25 years of monitoring, modeling, and understanding

The fate of hydrocarbons in the subsurface near Bemidji, Minnesota, has been investigated by a multidisciplinary group of scientists for over a quarter century. Research at Bemidji has involved extensive investigations of multiphase flow and transport, volatilization, dissolution, geochemical interactions, microbial populations, and biodegradation with the goal of providing an improved understanding of the natural processes limiting the extent of hydrocarbon contamination. A considerable volume of oil remains in the subsurface today despite 30 years of natural attenuation and 5 years of pump‐and‐skim remediation. Studies at Bemidji were among the first to document the importance of anaerobic biodegradation processes for hydrocarbon removal and remediation by natural attenuation. Spatial variability of hydraulic properties was observed to influence subsurface oil and water flow, vapor diffusion, and the progression of biodegradation. Pore‐scale capillary pressure‐saturation hysteresis and the presence of fine‐grained sediments impeded oil flow, causing entrapment and relatively large residual oil saturations. Hydrocarbon attenuation and plume extent was a function of groundwater flow, compound‐specific volatilization, dissolution and biodegradation rates, and availability of electron acceptors. Simulation of hydrocarbon fate and transport affirmed concepts developed from field observations, and provided estimates of field‐scale reaction rates and hydrocarbon mass balance. Long‐term field studies at Bemidji have illustrated that the fate of hydrocarbons evolves with time, and a snap‐shot study of a hydrocarbon plume may not provide information that is of relevance to the long‐term behavior of the plume during natural attenuation.

Minnesota↗

Characterization of Interactions between Surface Water and Near-Stream Groundwater along Fish Creek, Teton County, Wyoming, by Using Heat as a Tracer

Fish Creek, a tributary of the Snake River, is about 25 river kilometers long and is located in Teton County in western Wyoming near the town of Wilson. Local residents began observing an increase in the growth of algae and aquatic plants in the stream during the last decade. Due to the known importance of groundwater to surface water in the area, the U.S. Geological Survey (USGS), in cooperation with the Teton Conservation District, conducted a study to characterize the interactions between surface water and near-stream groundwater along Fish Creek. The study has two main objectives: (1) develop an improved spatial and temporal understanding of water flow (fluxes) between surface water and groundwater, and (2) use a two-dimensional groundwater-flow and heat-transport model to interpret observed temperature and hydraulic-head distributions and to describe groundwater flow near Fish Creek. The study is intended to augment hydrologic information derived from previously published results of a seepage investigation on Fish Creek. Seepage measurements provide spatially averaged gains and losses over an entire reach for one point in time, whereas continuous temperature and water-level measurements provide continuous estimates of gain and loss at a specific location. Stage, water-level, and temperature data were collected from surface water and from piezometers completed in an alluvial aquifer at three cross sections on Fish Creek at Teton Village, Resor's Bridge, and Wilson from October 2004 to October 2006. The flow and energy (heat) transport model VS2DH was used to simulate flow through the streambed of Fish Creek at the Teton Village cross section from April 15 to October 14, 2006, (183 recharge periods) and at the Resor's Bridge and Wilson cross sections from June 6, 2005, to October 14, 2006 (496 recharge periods). A trial-and-error technique was used to determine the best match between simulated and measured data. These results were then used to calibrate the cross-sectional models and determine horizontal and vertical hydraulic conductivities. The fluxes of groundwater into the stream or fluxes of stream water into the alluvial aquifer were estimated by using the calibrated VS2DH model for each cross section. Results of the simulations indicated that surface water/groundwater interaction and hydraulic properties were different at the three cross sections. At the most upstream cross section, Teton Village, Fish Creek flowed intermittently and continually gained relatively large quantities of water from April through September. During other times of the year, the stream was dry near the cross section. Saturated hydraulic conductivity set at 1x10-4 m/s in both the horizontal and vertical directions resulted in the best match between simulated and measured temperatures. The Resor's Bridge cross section, about midway between the other two cross sections, was near the point where perennial flow begins. At this cross section, the stream gained water from groundwater during high flow in late spring and summer, was near equilibrium with groundwater during August and September, and lost water to groundwater during the remainder of the year. Horizontal hydraulic conductivity set at 5x10-5 m/s and vertical hydraulic conductivity set at 1x10-5 m/s resulted in the best match between simulated and measured temperatures. The Wilson cross section, the most downstream site, was at USGS streamflow-gaging station 13016450. This part of the stream is perennial and was almost always gaining a small volume of water from groundwater. Saturated hydraulic conductivity set at 1x10-4 m/s in the horizontal direction and at 5x10-6 m/s in the vertical direction resulted in the best match between simulated and measured temperatures. Quantitative values of the flux from groundwater into surface water were estimated by using VS2DH and ranged from 1.1 to 6.6 cubic meters per day (m3/d) at the Teton Village cross section, from -3.8 to 7.4 m3/d at t

Scientific Investigations Report↗

Influence of groundwater pumping on streamflow restoration following upstream dam removal

We compared streamflow in basins under the combined impacts of an upland dam and groundwater pumping withdrawals, by examining streamflow in the presence and absence of each impact. As a qualitative analysis, inter-watershed streamflow comparisons were performed for several rivers flowing into the east side of the Central Valley, CA. Results suggest that, in the absence of upland dams supporting large reservoirs, some reaches of these rivers might develop ephemeral streamflow in late summer. As a quantitative analysis, we conducted a series of streamflow/groundwater simulations (using MODFLOW-2000 plus the streamflow routing package, SFR1) for a representative hypothetical watershed, with an upland dam and groundwater pumping in the downstream basin, under humid, semi-arid, and arid conditions. As a result of including the impact of groundwater pumping, post-dam removal simulated streamflow was significantly less than natural streamflow. The model predicts extensive ephemeral conditions in the basin during September for both the arid and semi-arid cases. The model predicts continued perennial conditions in the humid case, but spatially weighted, average streamflow of only 71% of natural September streamflow, as a result of continued pumping after dam removal.

California↗