Geology Reports⌕ Search

Geology topics

Richard G. Niswonger

Publications and source records attributed to Richard G. Niswonger.

63 records · Page 4Linked to original sources
USGS2008

GSFLOW - Coupled Ground-Water and Surface-Water Flow Model Based on the Integration of the Precipitation-Runoff Modeling System (PRMS) and the Modular Ground-Water Flow Model (MODFLOW-2005)

The need to assess the effects of variability in climate, biota, geology, and human activities on water availability and flow requires the development of models that couple two or more components of the hydrologic cycle. An integrated hydrologic model called GSFLOW (Ground-water and Surface-water FLOW) was developed to simulate coupled ground-water and surface-water resources. The new model is based on the integration of the U.S. Geological Survey Precipitation-Runoff Modeling System (PRMS) and the U.S. Geological Survey Modular Ground-Water Flow Model (MODFLOW). Additional model components were developed, and existing components were modified, to facilitate integration of the models. Methods were developed to route flow among the PRMS Hydrologic Response Units (HRUs) and between the HRUs and the MODFLOW finite-difference cells. This report describes the organization, concepts, design, and mathematical formulation of all GSFLOW model components. An important aspect of the integrated model design is its ability to conserve water mass and to provide comprehensive water budgets for a location of interest. This report includes descriptions of how water budgets are calculated for the integrated model and for individual model components. GSFLOW provides a robust modeling system for simulating flow through the hydrologic cycle, while allowing for future enhancements to incorporate other simulation techniques.

Techniques and Methods↗
USGS2007

Streambed infiltration and ground-water flow from the Trout Creek drainage, an intermittent tributary to the Humboldt River, north-central Nevada

Ground water is abundant in many alluvial basins of the Basin and Range Physiographic Province of the western United States. Water enters these basins by infiltration along intermittent and ephemeral channels, which originate in the mountainous regions before crossing alluvial fans and piedmont alluvial plains. Water also enters the basins as subsurface ground-water flow directly from the mountains, where infiltrated precipitation recharges water-bearing rocks and sediments at these higher elevations. Trout Creek, a typical intermittent stream in the Middle Humboldt River Basin in north-central Nevada, was chosen to develop methods of estimating and characterizing streambed infiltration and ground-water recharge in mountainous terrains. Trout Creek has a drainage area of about 4.8 × 10 7 square meters. Stream gradients range from more than 1 × 10 –1 meter per meter in the mountains to 5 × 10 –3 meter per meter at the foot of the piedmont alluvial plain. Trout Creek is perennial in short reaches upstream of a northeast-southwest trending normal fault, where perennial springs discharge to the channel. Downstream from the fault, the water table drops below the base of the channel and the stream becomes intermittent. Snowmelt generates streamflow during March and April, when streamflow extends onto the piedmont alluvial plain for several weeks in most years. Rates of streambed infiltration become highest in the lowest reaches, at the foot of the piedmont alluvial plain. The marked increases in infiltration are attributed to increases in streambed permeability together with decreases in channel-bed armoring, the latter which increases the effective area of the channel. Large quartzite cobbles cover the streambed in the upper reaches of the stream and are absent in the lowest reach. Such changes in channel deposits are common where alluvial fans join piedmont alluvial plains. Poorly sorted coarse and fine sediments are deposited near the head of the fan, while finer-grained but better sorted gravels and sands are deposited near the foot. All flow in Trout Creek is lost to infiltration in the upper and middle reaches of the channel during years of normal to below-normal precipitation. During years of above-normal precipitation, streamflow extends beyond the piedmont alluvial plain to the lower reaches of the channel, where high rates of infiltration result in rapid stream loss. The frequency and duration of streambed infiltration is sufficient to maintain high water contents and low chloride concentrations, compared with interchannel areas, to depths of at least 6 m beneath the channel. Streamflow, streambed infiltration, and unsaturated-zone thickness are all highly variable along intermittent streams, resulting in recharge that is highly variable as well. Average annual ground-water recharge in the mountainous part of the Trout Creek drainage upstream of Marigold Mine was estimated on the basis of chloride balance to be 5.2 × 10 5 cubic meters. Combined with an average annual surface runoff exiting the mountains of 3.4 × 10 5 cubic meters, the total annual volume of inflow to alluvial-basin sediments from the mountainous part of the Trout Creek is 8.6 × 10 5 cubic meters, assuming that all runoff infiltrates the stream channel. This equates to about 7 percent of average annual precipitation, which is about the same percentage estimated for ground-water recharge using the original Maxey-Eakin method.

Nevada↗
USGS2006

Documentation of the Unsaturated-Zone Flow (UZF1) Package for modeling Unsaturated Flow Between the Land Surface and the Water Table with MODFLOW-2005

Percolation of precipitation through unsaturated zones is important for recharge of ground water. Rain and snowmelt at land surface are partitioned into different pathways including runoff, infiltration, evapotranspiration, unsaturated-zone storage, and recharge. A new package for MODFLOW-2005 called the Unsaturated-Zone Flow (UZF1) Package was developed to simulate water flow and storage in the unsaturated zone and to partition flow into evapotranspiration and recharge. The package also accounts for land surface runoff to streams and lakes. A kinematic wave approximation to Richards? equation is solved by the method of characteristics to simulate vertical unsaturated flow. The approach assumes that unsaturated flow occurs in response to gravity potential gradients only and ignores negative potential gradients; the approach further assumes uniform hydraulic properties in the unsaturated zone for each vertical column of model cells. The Brooks-Corey function is used to define the relation between unsaturated hydraulic conductivity and water content. Variables used by the UZF1 Package include initial and saturated water contents, saturated vertical hydraulic conductivity, and an exponent in the Brooks-Corey function. Residual water content is calculated internally by the UZF1 Package on the basis of the difference between saturated water content and specific yield. The UZF1 Package is a substitution for the Recharge and Evapotranspiration Packages of MODFLOW-2005. The UZF1 Package differs from the Recharge Package in that an infiltration rate is applied at land surface instead of a specified recharge rate directly to ground water. The applied infiltration rate is further limited by the saturated vertical hydraulic conductivity. The UZF1 Package differs from the Evapotranspiration Package in that evapotranspiration losses are first removed from the unsaturated zone above the evapotranspiration extinction depth, and if the demand is not met, water can be removed directly from ground water whenever the depth to ground water is less than the extinction depth. The UZF1 Package also differs from the Evapotranspiration Package in that water is discharged directly to land surface whenever the altitude of the water table exceeds land surface. Water that is discharged to land surface, as well as applied infiltration in excess of the saturated vertical hydraulic conductivity, may be routed directly as inflow to specified streams or lakes if these packages are active; otherwise, this water is removed from the model. The UZF1 Package was tested against the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model for a vertical unsaturated flow problem that includes evapotranspiration losses. This report also includes an example in which MODFLOW-2005 with the UZF1 Package was used to simulate a realistic surface-water/ground-water flow problem that includes time and space variable infiltration, evapotranspiration, runoff, and ground-water discharge to land surface and to streams. Another simpler problem is presented so that the user may use the input files as templates for new problems and to verify proper code installation.

Techniques and Methods↗
USGS2005

Documentation of the Streamflow-Routing (SFR2) Package to Include Unsaturated Flow Beneath Streams - A Modification to SFR1

Many streams in the United States, especially those in semiarid regions, have reaches that are hydraulically disconnected from underlying aquifers. Ground-water withdrawals have decreased water levels in valley aquifers beneath streams, increasing the occurrence of disconnected streams and aquifers. The U.S. Geological Survey modular ground-water model (MODFLOW-2000) can be used to model these interactions using the Streamflow-Routing (SFR1) Package. However, the approach does not consider unsaturated flow between streams and aquifers and may not give realistic results in areas with significantly deep unsaturated zones. This documentation describes a method for extending the capabilities of MODFLOW-2000 by incorporating the ability to simulate unsaturated flow beneath streams. A kinematic-wave approximation to Richards' equation was solved by the method of characteristics to simulate unsaturated flow beneath streams in SFR1. This new package, called SFR2, includes all the capabilities of SFR1 and is designed to be used with MODFLOW-2000. Unlike SFR1, seepage loss from the stream may be restricted by the hydraulic conductivity of the unsaturated zone. Unsaturated flow is simulated independently of saturated flow within each model cell corresponding to a stream reach whenever the water table (head in MODFLOW) is below the elevation of the streambed. The relation between unsaturated hydraulic conductivity and water content is defined by the Brooks-Corey function. Unsaturated flow variables specified in SFR2 include saturated and initial water contents; saturated vertical hydraulic conductivity; and the Brooks-Corey exponent. These variables are defined independently for each stream reach. Unsaturated flow in SFR2 was compared to the U.S. Geological Survey's Variably Saturated Two-Dimensional Flow and Transport (VS2DT) Model for two test simulations. For both test simulations, results of the two models were in good agreement with respect to the magnitude and downward progression of a wetting front through an unsaturated column. A third hypothetical simulation is presented that includes interaction between a stream and aquifer separated by an unsaturated zone. This simulation is included to demonstrate the utility of unsaturated flow in SFR2 with MODFLOW-2000. This report includes a description of the data input requirements for simulating unsaturated flow in SFR2.

Techniques and Methods↗
USGS2002

Analysis of temperature profiles for investigating stream losses beneath ephemeral channels

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

Water Resources Research↗
USGS2001

Analysis of streambed temperatures in ephemeral channels to determine streamflow frequency and duration

Spatial and temporal patterns in streamflow are rarely monitored for ephemeral streams. Flashy, erosive streamflows common in ephemeral channels create a series of operational and maintenance problems, which makes it impractical to deploy a series of gaging stations along ephemeral channels. Streambed temperature is a robust and inexpensive parameter to monitor remotely, leading to the possibility of analyzing temperature patterns to estimate streamflow frequency and duration along ephemeral channels. A simulation model was utilized to examine various atmospheric and hydrological upper boundary conditions compared with a series of hypothetical temperature‐monitoring depths within the streambed. Simulation results indicate that streamflow events were distinguished from changing atmospheric conditions with greater certainty using temperatures at shallow depths (e.g., 10–20 cm) as opposed to the streambed surface. Three ephemeral streams in the American Southwest were instrumented to monitor streambed temperature for determining the accuracy of using this approach to ascertain the long‐term temporal and spatial extent of streamflow along each stream channel. Streambed temperature data were collected at the surface or at shallow depth along each stream channel, using thermistors encased in waterproof, single‐channel data loggers tethered to anchors in the channel. On the basis of comparisons with site information, such as direct field observations and upstream flow records, diurnal temperature variations successfully detected the presence and duration of streamflow for all sites.

Water Resources Research↗
USGS1999

Numerical simulation of vertical ground-water flux of the Rio Grande from ground-water temperature profiles, central New Mexico

An important gap in the understanding of the hydrology of the Middle Rio Grande Basin, central New Mexico, is the rate at which water from the Rio Grande recharges the Santa Fe Group aquifer system. Several methodologies-including use of the Glover-Balmer equation, flood pulses, and channel permeameters- have been applied to this problem in the Middle Rio Grande Basin. In the work presented here, ground-water temperature profiles and ground-water levels beneath the Rio Grande were measured and numerically simulated at four sites. The direction and rate of vertical ground-water flux between the river and underlying aquifer was simulated and the effective vertical hydraulic conductivity of the sediments underlying the river was estimated through model calibration. Seven sets of nested piezometers were installed during July and August 1996 at four sites along the Rio Grande in the Albuquerque area, though only four of the piezometer nests were simulated. In downstream order, these four sites are (1) the Bernalillo site, upstream from the New Mexico State Highway 44 bridge in Bernalillo (piezometer nest BRN02); (2) the Corrales site, upstream from the Rio Rancho sewage treatment plant in Rio Rancho (COR01); (3) the Paseo del Norte site, upstream from the Paseo del Norte bridge in Albuquerque (PDN01); and (4) the Rio Bravo site, upstream from the Rio Bravo bridge in Albuquerque (RBR01). All piezometers were completed in the inner-valley alluvium of the Santa Fe Group aquifer system. Ground-water levels and temperatures were measured in the four piezometer nests a total of seven times in the 24-month period from September 1996 through August 1998. The flux between the surface- and ground-water systems at each of the field sites was quantified by one-dimensional numerical simulation of the water and heat exchange in the subsurface using the heat and water transport model VS2DH. Model calibration was aided by the use of PEST, a model-independent computer program that uses nonlinear parameter estimation. Mean vertical hydraulic conductivities were estimated by model calibration and range from 1.5x10 -5 to 5.8x10 -6 meters per second (m/s). Mean simulated vertical ground-water flux for the BRN02 piezometer nest is 3.30x10 -7 m/s; for the COR01 piezometer nest is 3.58x10 -7 m/s; for the PDN01 piezometer nest is 4.22x10 -7 m/s; and for the RBR01 piezometer nest is 2.05x10 -7 m/s. Comparison of the simulated vertical fluxes and vertical hydraulic conductivities derived from this study with values from other studies in the Middle Rio Grande Basin indicate agreement between 1 and 3.5 orders of magnitude for hydraulic conductivity and within 1 order of magnitude for vertical flux.

New Mexico↗