Introduction and Characteristics of Flow

By James W. LaBaugh and Donald O. Rosenberry Chapter 1 of

Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water
Edited by Donald O. Rosenberry and James W. LaBaugh

Techniques and Methods Chapter 4D2

U.S. Department of the Interior U.S. Geological Survey

Contents
Introduction.....................................................................................................................................................5 Purpose and Scope........................................................................................................................................6 Characteristics of Water Exchange Between Surface Water and Ground Water..............................7 Characteristics of Near-Shore Sediments........................................................................................8 Temporal and Spatial Variability of Flow..........................................................................................10 Defining the Purpose for Measuring the Exchange of Water Between Surface Water and Ground Water...........................................................................................................................12 Determining Locations of Water Exchange. ....................................................................................12 Measuring Direction of Flow.............................................................................................................15 Measuring the Quantity of Flow........................................................................................................15 Measuring Temporal Variations in Flow..........................................................................................15 Methods of Investigation.............................................................................................................................16 Watershed-Scale Rainfall-Runoff Models. ......................................................................................16 Stream Discharge Measurements....................................................................................................17 Ground-Water Flow Modeling...........................................................................................................17 Direct Measurement of Hydraulic Properties.................................................................................20 Examination and Analysis of Aerial Infrared Photography and Imagery...................................20 Dye and Tracer Tests...........................................................................................................................21 Thermal Profiling..................................................................................................................................26 Use of Specific-Conductance Probes..............................................................................................26 Electrical Resistivity Profiling............................................................................................................27 Hydraulic Potentiomanometer (Portable Wells) Measurements. ................................................27 Seepage Meter Measurements........................................................................................................28 Biological Indicators...........................................................................................................................28 References.....................................................................................................................................................30

Figures
1. Summary of techniques that have been used for the measurement or estimation of water fluxes between surface water and ground water. ...................................................6 2. Typical hydraulic conditions in the vicinity of the shoreline of a surface-water body...............8 3. Decrease in seepage discharge with distance from shore...................................................8 4. Generalized hydrologic landscapes...........................................................................................9 5. Flow of salt tracer into a sandy lakebed, Perch Lake, Ontario............................................10 6. The length of the flow path and the direction of flow can vary seasonally and with distance from shore. ...........................................................................................................10 7. Example of ground-water exchange with Mirror Lake wherein lake water flows into ground water near shore, and ground water flows into the lake further from the shoreline.................................................................................................................................11 8. Example of the effect of transpiration on the water table and the direction of water flux between surface water and ground water..........................................................12

9. Example of rise and fall of water-table mounds at the edge of surface-water bodies and changes in flow direction......................................................................................13 10. Example of changes in flow direction related to onset and dissipation of a water-table mound adjacent to a river....................................................................................14 11. Example of flow interaction between surface water and ground water...........................15 12. Fluvial-plain ground-water and stream-channel interactions showing channel cross sections..............................................................................................................................16 13. Example of the use of hydrographs to determine the amount of ground-water discharge to a stream.................................................................................................................17 14. Example of the use of discharge measurements to determine surface-water/ ground-water interaction...........................................................................................................18 15. Example of model grid used for simulation of surface-water/ground-water interaction.....................................................................................................................................19 16. Example of shoreline segment definition for the calculation of water fluxes between surface water and ground water at a lake.............................................................21 17. Example of use of thermal infrared imagery to delineate areas of discrete and diffuse ground-water discharge to surface water. ................................................................22 18. Example of the use of dye to examine water fluxes between surface water and ground water................................................................................................................................23 19. Example of the use of a towed specific-conductance probe to identify groundwater discharge to surface waters..........................................................................................26 20. A, photograph of minipiezometer/hydraulic potentiomanometer in use; B, diagram of hydraulic-potentiomanometer system....................................................................................27 21. Full-section view of seepage meter showing details of placement in the sediment...............28 22. Photographs of marsh marigold in Shingobee Lake, Minnesota. ........................................29

Chapter 1 Introduction and Characteristics of Flow

By James W. LaBaugh and Donald O. Rosenberry

Introduction
Interest in the use and development of our Nations surface- and ground-water resources has increased significantly during the past 50 years (Alley and others, 1999; Hutson and others, 2004). At the same time, a variety of techniques and methods have been developed to examine and monitor these water resources. Quantifying the connection between surface water and ground water also has become more important because the use of one of these resources can have unintended consequences on the other (Committee on Hydrologic Science, National Research Council, 2004). In an attempt to convey the importance of the linkages and interfaces between surface water and ground water, the two have been described as a single resource (Winter and others, 1998). An improved understanding of the connection between surface and ground waters increasingly is viewed as a prerequisite to effectively managing these resources (Sophocleous, 2002). Thus, water-resource managers have begun to incorporate management strategies that require quantifying flow between surface water and ground water (Danskin, 1998; Bouwer and Maddock, 1997; Dokulil and others, 2000; Owen-Joyce and others, 2000; Barlow and Dickerman, 2001; Jacobs and Holway, 2004). The use of surface water (or ground water) can change the location, rate, and direction of flow between surface water and ground water (Stromberg and others, 1996; Glennon, 2002; Galloway and others, 2003). Pumping wells in the vicinity of rivers commonly cause river water to flow into the underlying ground-water body, which can affect the quality of the ground water (Childress and others, 1991; McCarthy and others, 1992; Lindgren and Landon, 1999; Steele and Verstraeten, 1999; Zarriello and Reis, 2000; Sheets and others, 2002). In some cases, ground water is pumped to provide water for cooling industrial equipment and then discharged into lakes, ponds, or rivers (Andrews and Anderson, 1978; Hutson and others, 2004). Ground water also may be pumped specifically to maintain lake levels for recreation purposes, especially during droughts (Stewart and Hughes, 1974; Mcleod, 1980; Belanger and Kirkner, 1994; Metz and Sacks, 2002). Surface water can be directed into surface basins where water percolates to the underlying aquifera process known as artificial recharge (Galloway and others, 2003). Ground-water discharge areas, where ground water flows into

surface water, can be important habitats for fish (Garrett and others, 1998; Power and others, 1999; Malcolm and others, 2003a, 2003b). Water in irrigation canals can flow or seep to an underlying aquifer, which eventually discharges water to rivers, thereby sustaining streamflow essential for the maintenance of fish populations (Konrad and others, 2003). Interest in the interaction of surface water and ground water is not confined to inland waters. This interaction has been studied in coastal areas because fresh ground-water supplies can be affected by intrusion of saltwater (Barlow and Wild, 2002). Beyond the issue of water supply for human consumption, increased attention has been given to the ground water that discharges to oceans and estuaries, both in terms of water quantity and quality (Bokuniewicz, 1980; Moore, 1996, 1999; Linderfelt and Turner, 2001). Discharge of fresh ground water to oceans and estuaries, also referred to as submarine ground-water discharge, is important in maintaining the flora and fauna that have evolved to exploit this source of fresh water in a saline environment (Johannes, 1980; Simmons, 1992; Corbett and others, 1999). Nitrate in submarine groundwater discharge to estuaries and coastal waters can result in eutrophication of those waters (Johannes, 1980; Johannes and Hearn, 1985; Valiela and others, 1990; Taniguchi and others, 2002). Withdrawals or pumping of ground water at near-shore, inland locations can reduce the submarine discharge of ground water offshore and change the environmental conditions of these settings (Simmons, 1992). Some coral reefs may be endangered by diminished submarine ground-water discharge (Bacchus, 2001, 2002). The variety of settings of interest for the examination of the interaction between surface water and ground water makes evident the need for methods to describe and quantify that flow. The exact method chosen for each setting will vary depending on the physical and hydrological conditions present in those settings, as well as the scale of the interaction. Some degree of measurement uncertainty accompanies each method or technique. Thus, it is prudent to consider using more than one method to examine the interaction between surface water and ground water. Because numerous techniques and methods are available to describe and quantify the flow between surface water and ground water, it is useful to provide water-resource investigators an overview of available techniques and methods, as well as their application.

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

Purpose and Scope

Several methods have been developed and applied to the study of the exchange between surface water and ground water (fig. 1). Different methods are better suited for characterizing or measuring flow over large or small areas. If an initial view of a considerable area or distance is needed to determine where measurable ground-water discharge is occurring, aerial infrared photography or imagery can be effective reconnaissance tools. On a smaller scale, some methods may involve direct measurement of sediment temperature or specific conductance

along transects within a surface-water body, or use of dyes or other tracers to indicate the direction and rate of water movement. The measurement of water levels in well networks in the watershed can be used to determine ground-water gradients relative to adjacent surface water, which in turn can indicate the direction and rate of flow between the surface-water body and the underlying aquifer. In streams and rivers, measurement of flow at the endpoints of a channel reach can reveal if the reach is gaining flow from ground water or losing flow to ground water. Addition of tracers to streams also can be used to determine surface-water interaction with ground water over a range of scales. Local interaction of surface water with ground

A E

D C

F E F

Figure 1. Summary of techniques that have been used for the measurement or estimation of water fluxes between surface water and ground water. Techniques illustrated include: (A) aerial infrared photography and imagery, (B) thermal profiling, (C) the use of temperature and specific-conductance probes, (D) dyes and tracers, (E) hydraulic potentiomanometers, (F) seepage meters, (G)well networks, and (H) streamflow measurements. (Artwork by John M. Evans, U.S. Geological Survey, retired.)

Introduction and Characteristics of Flow 7 water is measured by placing devices such as thermistors, minipiezometers, and seepage meters in the sediment, to monitor temperature gradients, hydraulic gradients, or quantity of flow. Determination of how interaction of surface water with ground water changes over time is made possible by using datarecording devices (data loggers) in conjunction with pressure transducers, thermistors, and water-quality probes. This report is designed to make the reader aware of the breadth of approaches (fig. 1) 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 nearsurface environment. 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. To that end, at this point it is useful to examine the content of each chapter in more detail. Chapter 1 provides an overview of typical settings in the landscape where interactions between surface water and ground water occur. The chapter reviews the literature, particularly recent publications, and describes many welldocumented methods for defining the flow between surface and ground waters. A brief overview of the theory behind each method is provided. Information is presented about the field settings where the method has been applied successfully, and, where possible, generalizes the requirements of the physical setting necessary to the success of the method. Strengths and weaknesses of each method are noted, as appropriate. This will aid the investigator in choosing methods to apply to their setting. For those already familiar with some of these methods, the review of recent literature provides information about improvements in these methods. Chapter 2 describes 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. The first method involves measurement of water levels in a network of wells in combination with measurement of the stage of the surface-water body to calculate gradients and then water flow. The second method involves the use of portable piezometers (wells) or hydraulic potentiomanometers to measure gradients. In the third method, seepage meters are used to directly measure flow across the sediment-water interface at the bottom of the surface-water body. Factors that affect measurement scale, accuracy, sources of error in using each of the methods, common problems and mistakes in applying the methods, and conditions under which each method is well- or ill-suited also are described. Chapter 3 presents an overview of methods that are commonly used in the hydrogeologic investigation and characterization of karst aquifers and in the study of water fluxes in karst terranes. Special emphasis is given to describing the techniques involved in conducting water-tracer tests using fluorescent dyes. Dye-tracer test procedures described herein represent commonly accepted practices derived from a variety of published and previously unpublished sources. Methods that are commonly applied to the analysis of karst spring discharge (both flow and water chemistry) also are reviewed and summarized. Chapter 4 reviews early work addressing heat as a tracer in hydrological investigations of the near-surface environment, describes recent advances in the field, and presents selected new results designed to identify the broad application of heat as a tracer to investigate surface-water/ground-water exchanges. An overview of field techniques for estimating water fluxes between surface water and ground water with heat is provided. To familiarize readers with flow conditions that may occur during their studies, the next section of Chapter 1 describes commonly observed interactions between surface water and ground water.

Characteristics of Water Exchange Between Surface Water and Ground Water

Most measurements made for the purpose of quantifying exchange between surface water and ground water are obtained at points within a short distance of the shoreline of the surface-water body. Shorelines represent the horizontal interface between ground water and surface water, an interface that is highly dynamic spatially and temporally. Because of the complex physical processes that occur in precisely the area where measurements are needed, it is important to understand those processes at shorelines and the range of potential changes in conditions at shorelines that occur over time. The following section elaborates these points.

8 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water Typically, a significant break in slope in the water table occurs where the horizontal surface of a lake, stream, or wetland intersects the sloping surface of the ground-water table (fig. 2). Because of this break in slope, ground-water flow lines diverge where they extend beneath and end at the sediment-water interface. Diverging flow lines indicate that the rate of flow per unit area is decreasing. Given homogeneous and isotropic conditions in the porous media adjacent to and beneath the sediment-water interface, seepage across the interface will decrease exponentially with distance from shore (fig. 3) (McBride and Pfannkuch, 1975; Pfannkuch and Winter, 1984). The movement of water between surface water and ground water can occur in a variety of settings or landscapes (fig. 4), each of which can be related to the break in slope of the water table defined by an upland adjacent to a lowland separated by an intervening steeper slope (Winter, 2001). Ground-water flow lines bend substantially beneath the sediment-water interface just before they intersect the surface-water body. Measurements of hydraulic-head gradients typically assume that the flow lines either are horizontal (in the case of comparing heads in near-shore wells with surfacewater stage) or vertical (in the case of inserting the screened intervals of wells to some depth beneath the sediment-water interface). In reality, the orientation of the flow lines are somewhere between horizontal and vertical as shown in figure 5.
Land surface
Water table

Surface water

Ground-water flow path

Figure 3. Decrease in seepage discharge with distance from shore (from Winter and others, 1998).

Characteristics of Near-Shore Sediments

Although some investigators have found that seepage decreases exponentially with distance from shore (Lee, 1977; Fellows and Brezonik, 1980; Erickson, 1981; Attanayake and Waller, 1988; Rosenberry, 1990), other studies report that the decrease in flow across the sediment-water interface is not

exponential because of heterogeneity of the sediment. One of the early fndings of a departure from what would be expected in a homogeneous, isotropic setting was reported by Woessner and Sullivan (1984) in their study of Lake Mead, Nevada. At many of the transects across which they collected data in Lake Mead, they found seepage did not decrease exponentially, and furthermore, that seepage sometimes decreased and then increased with distance from shore. They reported a large variability in seepage with distance from shore. This variability was attributed to heterogeneity in the sediments in the vicinity of the sediment-water interface. Krabbenhoft and Anderson (1986) also reported that seepage was focused in a gravel lens that intersected the lakebed some distance from shore at Trout Lake, Wisconsin. It now generally is recognized that aquifers adjacent to and beneath surface-water bodies rarely can be considered homogeneous, and usually are not isotropic. Many processes act to create heterogeneity at the sedimentwater interface. A few are listed below. 1. Fluvial processesDepositional and erosional processes occur nearly constantly in streambeds and riverbeds, making heterogeneity a significant feature in these sediment-water interfaces. Organic deposits commonly are buried by deposition of inorganic material, resulting in interlayering of these different sediment types. Channel aggradation and flood scour can cause a shoreline to shift laterally many meters. Seasonal erosion and deposition related to spring floods also create a temporal component to the heterogeneity. 2. Edge effectsShoreline erosion and deposition related to wave action in lakes, large wetlands, and rivers create heterogeneity at the sediment-water interface. Waves erode banks, which subsequently fail as new material slumps into the surface-water body. Fine-grained sediments are moved away from shore, often leaving a cobble- to boulder-sized pavement at the shoreline. Sediment deposition by overland flow commonly results in near-shore, fan-shaped deposits following heavy rainfall. Waves also rework sediments following slump events or sediment transport associated with overland flow, causing movement of fine-grained materials into voids created by movement of cobble- to boulder-sized sediments. In addition, changing surface-water stage causes the position

eg me

nt

Water-table well
es

Break in slope

m A

Groun

d-wate r flow

lines

b L
Flow pa ss

omain es beneath la ke and is outside of local-flow d

Figure 2. Typical hydraulic conditions in the vicinity of the shoreline of a surface-water body. (Artwork by Donald O. Rosenberry, U.S. Geological Survey.)

Sh

ore

lin

Surface-water body

Introduction and Characteristics of Flow 9

MOUNTAIN VALLEY
Water table Seepage face Break in slope
Direction of ground-water flow Direction of regional ground-water flow

RIVERINE VALLEY
Water table Regional upland

Land surface

Water table

Flood levels

Direction of local ground-water flow

A
PLAYA

D
Regional upland

COASTAL TERRAIN

Water table Terrace Land surface Ocean

Water table Direction of regional ground-water flow Direction of local ground-water flow

Direction of local ground-water flow

B
PLATEAU AND HIGH PLAINS
Land surface
Water table Direction of local ground-water flow Direction of regional ground-water flow

E
HUMMOCKY TERRAIN
Water table

Direction of local ground-water flow Direction of regional ground-water flow

Figure 4. Generalized hydrologic landscapes: A, narrow uplands and lowlands separated by a large steep valley side (mountainous terrain); B, large broad lowland separated from narrow uplands by steeper valley sides (playas and basins of interior drainage); C, small narrow lowlands separated from large broad uplands by steeper valley side (plateaus and high plains); D, small fundamental hydrologic landscape units nested within a large fundamental hydrologic landscape unit (large riverine valley with terraces); E, small fundamental hydrologic landscape units superimposed on a larger fundamental hydrologic landscape unit (coastal plain with terraces and scarps); F, small fundamental hydrologic landscape units superimposed at random on large fundamental hydrologic landscape units (hummocky glacial and dune terrain) (from Winter, 2001, copyright the American Water Resources Association, used with permission).

of the shoreline to change over time, resulting in lateral movement of all of the previously mentioned depositional and erosional processes that occur at the shoreline. Accumulation of organic debris, including buried logs and decayed plant matter, also contributes to heterogeneity as it is incorporated with the inorganic sediments, particularly on the downwind shores of surface-water bodies. In surface-water bodies that are ice covered during winter, ice rafting during fall and spring, when ice is forming or when the ice cover is melting, can substantially rework sediments at the downwind shoreline.

3. Biological processesBenthic invertebrates constantly rework sediments, particularly organic sediments, as they carry out their life cycles. Bioturbation and bioirrigation are important processes for organic sediments in deeper water environments, but it can be significant in some nearshore settings also. Aquatic birds disturb the sediment as they search for benthic invertebrates, and fish rework sediments as they create spawning redds. Beavers and muskrats can make large-scale disturbances by removing considerable amounts of sediments for lodges and passageways, and the construction of dams.

10 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water
C'
SM19 SM6 SM3 SM20 SM21

73
INJECTION WELLS
ML10

43

23 10 3

ML39 0 0 2

ML36 1 METER 4 FEET ML29 ML26

ML37

ML34 ML7

ML19

ML15

SM6 - seepage meter ML29 - multi-level well

ML8

ML9

Figure 5. Flow of salt tracer into a sandy lakebed, Perch Lake, Ontario. Numbers in shaded areas (representing center of mass of salt plume) are the days following salt injection (modified from Lee and others, 1980, copyright 1980 by the American Society of Limnology and Oceanography, Inc., used with permission).

Temporal and Spatial Variability of Flow

Flow across the sediment-water interface commonly changes in direction and velocity temporally and spatially. Many occurrences of spatial and temporal variability in the exchange between surface water and ground water are described in the literature; a few examples are provided herein. Some of this variability is summarized in figure 6. In this illustration, water flows from the surface-water body to ground water through the bottom sediments located beyond a low-permeability layer some distance from shore. Yet closer to shore, ground water flows into the surface-water body. Finally, near the shore a depression in the water table created by evapotranspiration causes flow out of the lake. At the south shoreline of Mirror Lake, New Hampshire, water flows from the lake to ground water between the shoreline and approximately 8 meters from shore, and beyond that point, flow from ground water to the lake occurs (fig. 7) (Asbury, 1990; Rosenberry, 2005). In many settings, evapotranspiration during the summer months can depress the water table adjacent to the shoreline of wetlands, streams, and lakes below the level of the surfacewater body (fig. 8) (Meyboom, 1966, 1967; Doss, 1993; Winter and Rosenberry, 1995; Rosenberry and others, 1999; Fraser and others, 2001). As a result, seasonal, and sometimes diurnal, reversals in flow between surface water and ground water may occur at the shoreline. The changes in direction of flow between surface water and ground water result from fluctuations in the amount of water removed from the water table because of evapotranspiration by plants along the margins of the surface-water body. On a seasonal basis, once evapotranspiration ceases to remove water from the near-shore regions, the near-shore depression in the water table dissipates, which then allows ground water to flow into the surface-water body.

On a diurnal basis, more evaportranspiration in the day and less at night can cause the water table to fluctuate between levels below and above the adjacent surface-water level. In many locations, water-table mounds can develop at the edge of surface-water bodies. Many studies have shown transient water-table mounds that form in response to precipitation or snowmelt (fig. 9) (see, for example, Winter, 1986; Rosenberry and Winter, 1997; Lee and Swancar, 1997). Most of these water-table mounds were of short duration and formed in response to large rainfall events. Reversals of flow of longer duration also occur at some settings. Jaquet (1976) reported a reversal of flow along part of the shoreline at Snake Lake, Wisconsin, following spring thaw and considerable rainfall (19 centimeters) over a 5-week period that persisted for several months.

Lan

d su

rfac

Surface-water body

Low-permeability layer

Figure 6. The length of the flow path and the direction of flow can vary seasonally and with distance from shore. (Artwork by Donald O. Rosenberry, U.S. Geological Survey.)

Introduction and Characteristics of Flow 11

A
Transect line

Mirror Lake

Study site Study site Outlet 72

44

Hubbard

Brook

NEW HAMPSHIRE

93

0 0 Line of equal hydraulic head (dashed where interval is variable) 600

300 1200

600 METERS 1800 2400 FEET

260 250

ALTITUDE, IN METERS ABOVE SEA LEVEL

240 230 220 210 200 190 180 170 160 150

Water table Piezometer point Study site Water-table well Mirror Lake 214 213.36 meters Hubbard Brook

5 23
23 0
225
220

Lake sediments 215 215 Glacial drift

2.1 21

3.2 21 .3 213 214

21 0.9

207.9

Direction of ground-water flow Crystalline bedrock

+4 85 120 54 20 153 29 23 0.6 0.4 2 3

+3

+0.3 Seepage rates, in centimeters/day 0 2 METERS 2 4 6 8 FEET

0.6 11 12

2 4 6

Well

Shoreline

Study site

Figure 7. Example of ground-water exchange with Mirror Lake wherein lake water flows into ground water near shore (shown by negative numbers), and ground water flows into the lake farther from the shoreline (shown by positive numbers). (Modified from Rosenberry, 2005, copyright 2005 by the American Society of Limnology and Oceanography, Inc., used with permission.)

12 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water
Transpiration
surfa ce

Land

Figure 8. Example of the effect of transpiration on the water table and the direction of water flux between surface water and ground water (from Winter and others, 1998).

Surface water

Water table during dormant season Water table during growing season

Ground-water levels adjacent to streams also fluctuate in response to the rise and fall of water in the stream (Winter, 1999). An example of the resulting changes in flow direction between a stream and ground water is illustrated by data from the Cedar River in Iowa (fig. 10). In flowing waters, movement of surface water into the subsurface and out again occurs both at the bottom of the stream channel and beneath upland areas between bends in the open channel (fig. 11). This transient flow of surface water into and out of the subsurface is also known as hyporheic flow (Orghidan, 1959). Ground water flowing toward a surface-water body may discharge directly into that body or mix with hyporheic flow prior to emerging into open-water flow. The various interactions with ground water include situations in which flow is parallel to the stream (fig. 12) and does not intersect the surface water (Woessner, 1998, 2000).

Defining the Purpose for Measuring the Exchange of Water Between Surface Water and Ground Water
Water-resource investigators and water-resource managers have many reasons to quantify the flow between surface water and ground water. Perhaps the most common reasons include: calculating hydrological and chemical budgets of surface-water bodies, collecting calibration data for watershed or ground-water models, locating contaminant plumes, locating areas of surface-water discharge to ground water, improving their understanding of processes at the interface between surface water and ground water, and determining the relation of water exchange between surface water and ground water to aquatic habitat. For many investigations, it is sufficient to make a qualitative determination regarding the direction and relative magnitude of flow, either into or out of the surfacewater body. Methods for quantifying flows should be selected to be appropriate for the scale of the study. For a watershed-scale study in which multiple basins may be involved, small-scale flow phenomena, such as near-shore depressions in the water table or spatial variability of flux related to geologic

variability, likely are of little importance to the overall study goal. In such watershed-scale studies, the net flux integrated over an entire stream reach, or lake, or wetland often is the desired result. Watershed-scale flow modeling, ground-water flow modeling, flow-net analysis, or dye- and geochemicaltracer tests, often are used in such large-scale studies, studies on the order of hundreds of meters or a kilometer or more in length or breadth. If the goal of a study is to identify and (or) delineate zones or areas of flow of surface water to ground water, or flow of ground water to surface water, smaller scale spatial and temporal variations in flow become important, and measurement tools that provide results over an intermediate scale, many tens to hundreds of meters should be selected. In many instances, measurement of surface-water flow at two places some distance apart in a segment of stream, which enables calculation of gains or losses in flow in the segment, is appropriate for these types of studies. For local, small-scale studies in which flow to or from surface water may be focused, small-scale tools such as seepage meters, small portable wells (minipiezometers or hydraulic potentiomanometers), and buried temperature probes may be most appropriate. Devices designed to measure flow in a small area are known as seepage meters because the term seep refers to a small area where water moves slowly to the land surface (USGS Water Basics Glossary http://capp.water.usgs.gov/GIP/h2o_gloss/). Seepage is defined as the slow movement of water through small cracks, pores, interstices, and so forth, of a material into or out of a body of surface or subsurface water (USGS Water Science Glossary of Terms http://ga.water.usgs.gov/edu/dictionary.html#S). Once the water-resource investigator has decided on the purpose of the study and the scale of the investigation, methods of investigation can be chosen to most effectively determine where an exchange between surface water and ground water is taking place, the direction of flow, the rate or quantity of that flow, and whether the rate and direction of flow changes over time.

Determining Locations of Water Exchange

The investigator who wishes to determine where water exchange is taking place between surface water and ground water has many options, particularly in the case of groundwater discharge to surface water. Reconnaissance tools useful over larger areas, such as dye-tracer tests, aerial photography and imagery, temperature and specific-conductance probes, and surface-water discharge measurements, can be supplemented by reconnaissance tools useful in smaller areas of interest, such as seepage meters, minipiezometers, and biological indicators.

Introduction and Characteristics of Flow 13 A

180 160 140 120 100 80 60 40
125

B WT - 15

B'

Pivot point Midlake

ALTITUDE, IN FEET, ABOVE OR BELOW SEA LEVEL

127.24 127.66 2PN-9 No data 2PN-20 126

12 7

128.09

Lake Lucerne
1

125.11

124.90 125.96 2PN-35 125.23 2PN-50

125 124 123 122 121 120

125.77 1PN-7 125.73 1PN-20 125.67 1PN-35 125.60 1PN-50 119.48 1PN-65 119.68 1PN-22

12

126

EXPLANATION
Equipotential lineInterval 1 foot. Datum is sea level Direction of ground-water flow Well and number HeadIn feet above sea level Water table
1

20 Sea level 20 40
1PN-7 119.48 2PN-130 125.77

Value for 8/8/85

1,500 1,750 2,000 2,250 2,500

119.54 1PN-155

250

500

750

1,000

1,250

2,750

3,000 3,250

B
180 160 140 B WT-15

DISTANCE, IN FEET, FROM WELL WT-15 ALONG PART OF TRANSECT BB

VERTICAL SCALE GREATLY EXAGGERATED
B'

Pivot point Midlake

ALTITUDE, IN FEET, ABOVE OR BELOW SEA LEVEL

129.88

127.37 2PN-9 127.50 2PN-20

130.02

120 100 80 60 40 20 Sea Level 20 40

129

Lake Lucerne
125.96 125.70 Midlake
126

127.25 1PN-20

128.37 1PN-7

WT-17
130

2PN

1PN

128

127.45 2PN-50

127

128.14 2PN-35

121

122

124 123

125

127.19 127.12 121.46 120.34

1PN-35 1PN-50

1PN-65 1PN-72

EXPLANATION
121

2PN-35 120.61 2PN-130 128.14

Equipotential lineInterval 1 foot. Datum is sea level Direction of ground-water flow Well and number HeadIn feet above sea level Water table

120.23 1PN-155

250

500

750

1,000

1,250

1,500

1,750

2,000

2,250

2,500

2,750

3,000

WT - 17
7

2PN

1PN

129

128

127

DISTANCE, IN FEET, FROM WELL WT-15 ALONG PART OF TRANSECT BB

VERTICAL SCALE GREATLY EXAGGERATED

Figure 9. Example of rise and fall of water-table mounds at the edge of surfacewater bodies and changes in flow direction (from Lee and Swancar, 1997). A, Vertical distribution of head showing downward head gradient conditions, August 6, 1985. B, Vertical distribution of head during high water-level conditions, October 17, 1985.

126

3,250

14 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water A
210 208 206 204 202 200 198 196 194 192 Cedar River (0.12)*
0.3
0.2 0.3 0.2 0.1 0.1 0.3

Sampling period February 2022

Water table

Geologic boundary

B
210 208 206 204 202 200 198 196 194 192 ALTITUDE, IN METERS Cedar River (0.51)*
0.5

Sampling period March 2022

Water table

0.2

.5 0.6 0.6 0 0.4 0.3

0.5

0.5

0.2 0.3 0.1

0.1

Geologic boundary

C
210 208 206 204 202 200 198 196 194 192 Cedar River (0.21)*
0.5

Sampling period April 35

0.2

0.3

0.4

0.2

Water table

0.1

Geologic boundary
Vertical Exaggeration 4 Datum is sea level

0 0

20 100

40

60 METERS 200 FEET

EXPLANATION
Fine- to coarse-grained silty/clayey sand Fine- to coarse-grained sand Till
0.1

(X.XX)*

Line of equal atrazine concentration, in micrograms per liter Concentration of atrazine in river at time of sampling, in micrograms per liter Well location Direction of ground-water flow

Figure 10. Example of changes in flow direction related to onset and dissipation of a water-table mound adjacent to a river (from Squillace and others, 1993, used in accordance with usage permissions of the American Geophysical Union wherein all authors are U.S. Government employees). Hydrogeologic sections for part of the Cedar River, Iowa, for three periods in 1990. A, Movement of ground water into the river prior to a period of high river stage. B, Movement of river water into the contiguous aquifer during high river stage. C, Return of some of the water from the aquifer during declining river stage.

Introduction and Characteristics of Flow 15

A. View from Above

Bank

Gravel bar

Figure 11. Example of flow interaction between surface water and ground water (from Dumouchelle, 2001). Schematic of flow in Chapman Creek, west-central Ohio. (Arrows indicate direction of flow. Diagrams not to scale.)

Stream

Gravel bar

Bank

B. Sectional View

Water surface

commonly known as a seepage run (Harvey and Wagner, 2000). If the amount of flow in the stream has increased over the reach, the increase may be attributed to ground-water discharge to the stream. If flow in the selected reach of stream has decreased, the decrease may be attributed to surface water flowing into ground water. It is important to recognize, however, that the direction of flow indicated by any change in streamflow is a net direction over the selected reach, and that within the reach, water may be moving into and out of the stream (and conversely, into and out of the underlying aquifer). It is important to account for any inflows or outflows within the stream reach, such as diversions for irrigation or channelized return flows from fields.

Sand and gravel

Measuring the Quantity of Flow

The volume of water flowing between surface water and ground water, either as surface water into ground water or ground water into surface water, can be measured directly with seepage meters. Measurement of changes in water temperatures over time at a specific site above the sediment, at the sediment-water interface, and within the sediment makes possible the determination of the amount of water exchange occurring between surface water and ground water. The exchange of water between surface water and ground water also can be examined and estimated by using dye tracer tests or by using other tracers. Such dyes or tracers are added directly to a stream and then their concentrations are measured at some point or points downstream. Changes in the concentration of the dye or other tracer over time downstream from where they are injected enables calculation of ground-water inputs.

Till

Measuring Direction of Flow

Comparison of surface-water levels and adjacent groundwater levels indicates direction of flow. If the surface-water level is higher than adjacent ground-water levels, the direction of flow is from the surface water to ground water. If the opposite is the caseground-water levels are higher than nearby surface-water levelsthen the direction of flow is from ground water to surface water. In addition to indicating the direction of flow, water-level measurements provide information about the magnitude of the hydraulic gradients between surface water and ground water. In some instances, however, these gradients can be altered locally. For example, vegetation between the wells and the edge of the surface-water body can transpire sufficient water to cause a local depression in the water table close to the edge of the surface water (Meyboom, 1966, 1967; Doss, 1993; Rosenberry and Winter, 1997; Fraser and others, 2001). Thus, it can be important to measure the direction of flow at a local scale using portable wells, minipiezometers, or hydraulic potentiomanometers. Another way to determine if a section of stream or river is receiving ground-water discharge or is losing water to the underlying aquifer is by measurement of surface-water flow at two places some distance apart in a reach of stream, a practice

Measuring Temporal Variations in Flow

In many instances, the rate of exchange between surface water and ground water varies over time scales of hours, days, or months. The direction of flow also may reverse on a seasonal basis or temporarily during a flood, for example. Measuring temporal variation in the rate of water exchange requires multiple measurements over these time periods. Measuring devices equipped with data recorders (data loggers) enable the investigator to record repeated measurements at specified time intervals to document temporal changes.

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

Figure 12. Fluvial-plain ground-water and stream-channel interactions showing channel cross sections classified as: A, gaining; B and C, losing; D, zero exchange; and E, flow-through. The stream is dark blue. The water table and stream stage (thicker lines), ground-water flow (arrows), and equipotential lines (dashed) are shown (from Woessner, 1998, copyright American Institute of Hydrology, used with permission).

Methods of Investigation
Common methods to examine exchange of water between surface-water bodies and ground-water bodies are described below. Some of these methods make use of already installed hydrological instruments and existing data, rather than requiring the investigator to make measurements of hydrologic characteristics. When using such methods, however, the investigator may install wells, stream-gaging equipment, or rain gages, as needed, to obtain sufficient data to make the application of methods possible and the results less uncertain. Other methods require that the investigator make additional, specific measurements or observations of hydrological, physical, chemical, or biological characteristics.

Watershed-Scale Rainfall-Runoff Models

Many analytical and numerical models that relate precipitation, ground-water recharge, and ground-water discharge to temporal variability of flow in a stream have been developed. A fundamental assumption in these models is that streamflow is an integrated response to these processes over the streams watershed, and that ground-water discharge to the

stream provides the steady flow in the stream between rainfall events, commonly referred to as baseflow. Analytical models generally determine baseflow through hydrograph separation techniques. Several automated routines have been developed to assist in this determination (Rutledge, 1992; Rutledge, 1998) (fig. 13). Other analytical methods also have been used to quantify the interaction between ground water and surface water, including an analytic-element method (Mitchell-Bruker and Haitjema, 1996) and a nonparametric regression model (Adamowski and Feluch, 1991). Several numerical models commonly referred to as rainfall-runoff models have been developed; these models areally divide watersheds and subwatersheds and calculate hydrologic parameters for each smaller area (for example, Federer and Lash, 1978; Leavesley and others, 1983, 1996, 2002; Beven and others, 1984; Beven, 1997; Buchtele and others, 1998). Rainfall-runoff models generally are calibrated to match river flow at the outlet of a watershed or subwatershed. Some models include the ground-water component of flow in each area. The current trend is to couple distributed-area watershed-scale models with ground-water flow models in order to better determine the temporal and spatial variability of the interaction between ground water and surface water (for example, Leavesley and Hay, 1998; Beven and Feyen, 2002).

Introduction and Characteristics of Flow 17

A
10,000

FLOW, IN CUBIC FEET PER SECOND

1,000

Streamflow

100
Pulse

10

Figure 13. Example of the use of hydrographs to determine the amount of ground-water discharge to a stream (from Rutledge, 2000). Hydrographs of streamflow for Big Hill Creek near Cherryvale, Kansas, for March 1974 (blue circles and dashed line), and hydrograph of estimated ground-water discharge using the PULSE model (red line). (Note: In each example, the total recharge modeled is 0.73 inch, which is the same as the total recharge estimated from RORA [a recession-curve-displacement method for estimating recharge] for this period. In example A, recharge is modeled as 0.65inch on day 69 and 0.08 inch on day 74. In example B, recharge is modeled as a gradual process that is constant from day 68 to day72).

60

70

80

90

TIME, IN DAYS

B
10,000 Streamflow 1,000

100 Pulse 10

60

70

80

90

TIME, IN DAYS

Stream Discharge Measurements

Measurements of stream discharge (Rantz and others, 1982a, b; Oberg and others, 2005) made as part of seepage runs (described earlier) can be used to determine the occurrence and rate of exchange of water between surface water and ground water in streams and rivers (fig. 14). The results of seepage runs have been used to provide an integrated value for flow between a stream and ground water along a specific stream reach. This method works well in small streams, but for larger streams and rivers, the errors associated with the measurement of flow in the channel often are greater than the net exchange of water to or from the stream or river. This method also requires that any tributaries that discharge to a stream along the reach of interest be measured and subtracted from the downstream discharge measurement. Likewise, withdrawals from the stream, such as that for irrigation, must be measured and added to the downstream discharge measurement.

The application of seepage-run data, however, is limited by the ratio of the net flow of water to or from the stream along a stream reach to the flow of water in the stream. The net exchange of water across the streambed must be greater than the cumulative errors in streamflow measurements. For example, if the errors in the stream discharge measurements are 5 percent of the true, actual flow, then according to the rules of error propagation, in order to be able to detect the net flow of water to or from the stream along the reach of interest, the value of net flow must be greater than 7 percent of the streamflow. Despite these limitations, many hydrologic studies have made use of this method with good results [for example, Ramapo River, New JerseyHill and others (1992); Bear River, Idaho and UtahHerbert and Thomas (1992); Souhegan River, New HampshireHarte and others (1997); Lemhi River, IdahoDonato (1998); constructed stream channel Baden-Wrttemberg, GermanyKaleris (1998)]. The information gained from seepage runs can be enhanced with data obtained by using other techniques such as minipiezometers, seepage meters, temperature and specific-conductance measurements to better define surface-water/ground-water fluxes [for example, creeks and rivers in the Puget Sound area of WashingtonSimonds and others (2004); and Chapman Creek, OhioDumouchelle (2001)]. Seepage-run results also can provide estimates of hydraulic conductivity of the streambed on a scale appropriate for ground-water flow modeling (Hill and others, 1992).

FLOW, IN CUBIC FEET PER SECOND

Ground-Water Flow Modeling

Since 1983, most investigators who have used the numerical modeling approach in the quantification of flows between surface water and ground water have used the U.S. Geological Survey MODFLOW modular modeling code (Harbaugh and others, 2000). This finite-difference model contains an original river package that can simulate flows to or from a river, assuming the river stage does not change during a specified time period (referred to as a stress period in MODFLOW), but can change from one time period to the next. Several other MODFLOW modules or packages also have been developed to simulate fluxes between surface water and ground water. These include streamflow routing packages (Prudic, 1989;

18 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water
11340' 11350' Salmon 4510'

BOR gaging station at L-1 diversion BOR gaging station at L-3a diversion 4510'

EXPLANATION
14 Seepage run reach and identification number Gaging stationRecords kept by U.S. Geological Survey (USGS) or Bureau of Reclamation (BOR) Other discharge measurement site Arbitrary boundary between upper and lower Lemhi River Basin

14 13 12

USGS gaging station below L-5 diversion near Salmon BOR gaging station at Barracks Lane 11330'

11
Le

10
i mh

LOWER BASIN

4500'

4500'

Tendoy

USGS gaging station near Lemhi

Lemhi
4550'

6
4450' 11330'

5 4
Ri ve r

11310'

BOR gaging station at Bureau of Land Management McFarland Campground

UPPER BASIN
2 1
Leadore BOR gaging station at Leadore 4440'

4440' 11350' 11340'

4430'

4430'

11330' 11320' 11310' 0 2 0 2 4 6 4

11300' 8 KILOMETERS 6 8 MILES

Figure 14. Example of the use of discharge measurements to determine surface-water/ground-water interaction (from Donato, 1998). Seepage run reaches, gaging stations, and discharge measurement sites in the Lemhi River Basin, east-central Idaho, August and October 1997.

Introduction and Characteristics of Flow 19 A

13,000 FEET

13,000 FEET PLAN VIEW

B
LAYER 1 2 3 4 5

50 FEET

CROSS-SECTIONAL VIEW

Figure 15. Example of model grid used for simulation of surface-water/ground-water interaction (from Merritt and Konikow, 2000). The lateral and vertical grid discretization for test simulation 1: A, Plan viewShaded area is the surface extent of lake cells in layer 1. Interior grid dimensions are 500 and 1,000 meters. Border row/column cells are 250 feet thick. The locations of hypothetical observations wells are denoted by a and b. B, Cross-sectional viewShaded area is the cross section of the lake. Although a nominal 10-foot thickness is shown for layer 1, the upper surface of layer 1 is not actually specified, and lake stages and aquifer water-table altitudes may rise higher than the nominal surface shown above.

20 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water Prudic and others, 2004), a reservoir package (Fenske and others, 1996), a lake package (Cheng and Anderson, 1993), and more recently, a more elaborate lake package (Merritt and Konikow, 2000) (fig. 15). More advanced MODFLOW-based programs have been developed to couple one-dimensional, unsteady streamflow routing with MODFLOW (Jobson and Harbaugh, 1999; Swain and Wexler, 1996). Advances also are being made in coupling MODFLOW with watershed models that simulate many of the surface-water processes within a basin (Sophocleous and others, 1999; Sophocleous and Perkins, 2000; Niswonger and others, 2006). One of the most challenging aspects of coupling ground-water and surface-water models has been representation of flow through the unsaturated zone beneath a stream. Two new programs have recently been developed for MODFLOW to simulate one-dimensional (Niswonger and Prudic, 2005) and threedimensional (Thoms and others, 2006) flow in the unsaturated zone. Many of these packages require a determination of the transmissivity of the sediments at the interface between the aquifer and the surface-water body. Transmissivity is determined by multiplying hydraulic conductivity by the thickness of the lakebed or riverbed sediments. Shoreline segments are delineated/selected on the assumption that the gradient between a nearby well and the surface-water body, the hydraulic conductivity of the sediments, and the cross-sectional area through which water flows to enter or leave the lake, are uniform along the entire segment (fig. 16). Flows through each segment are summed for the entire surface-water body to compute net flow. The scale of the shoreline segments, and the scale of the study, depend on the scale of the physical setting of interest and the density of monitoring wells. Further detail regarding this method is provided in Chapter 2, in the section Wells and Flow-Net Analysis.

Examination and Analysis of Aerial Infrared Photography and Imagery

Aerial infrared photography and imagery have been used to locate areas of ground-water discharge to surface waters (Robinove, 1965; Fischer and others, 1966; Robinove and Anderson, 1969; Taylor and Stingelin, 1969). This technique is effective only if the temperatures of surface water and ground water are appreciably different. Information obtained from infrared scanners can be captured electronically or transferred to film, on which tonal differences correspond to differences in temperature (Robinove and Anderson, 1969; Banks and others, 1996) (fig. 17). Published studies indicate tonal differences corresponding to a difference in temperature of approximately 2 degrees Celsius are distinguishable (Pluhowski, 1972; Rundquist and others, 1985; Banks and others, 1996). Within the limits of the ability of infrared imagery to distinguish temperature differences between surface water and ground water, the inspection of such imagery enables more rapid identification of gaining reaches in streams over large areas than can be accomplished by stream surveys that measure temperature directly (Pluhowski, 1972). Another advantage of this method in identifying areas of ground-water discharge to surface-water bodies is its application where using other techniques such as dye tracing or direct temperature measurements are impractical, or access on the ground is difficult (Campbell and Keith, 2001) or dangerous (Banks and others, 1996). Using thermal-infrared imagery to distinguish zones of ground-water discharge is practical for locating diffuse and focused ground-water discharge (Banks and others, 1996). This capability has been demonstrated in a variety of environments. Examples for lakes are Crescent Lake, Nebraska (Rundquist and others, 1985), where the flow is diffuse and occurs over a large area, and Great Salt Lake and Utah Lake in Utah (Baskin, 1998), where the ground-water flow into the lake is focused at springs. Campbell and Keith (2001) found the technique useful in locating many springs flowing into streams and reservoirs in northern Alabama. Examples for estuaries are creeks flowing into Chesapeake Bay, Maryland, and the shorelines of the Gunpowder River and the Chesapeake Bay into which the river flows (Banks and others,

Direct Measurement of Hydraulic Properties

The relation between the stage of a surface-water body and the hydraulic head measured in one or more nearby water-table wells can be used to calculate flows of water between surface water and ground water [Williams Lake, MinnesotaLaBaugh and others (1995); Vandercook Lake, WisconsinWentz and others (1995); large saline lakes in central AsiaZekster (1996); Lake Lucerne, FloridaLee and Swancar (1997); Waquoit Bay, Cape Cod, MassachusettsCambareri and Eichner (1998); Otter Tail River, MinnesotaPuckett and others (2002)]. The Darcy equation (eq. 1) is used to calculate flow between ground water and surface water along specific segments of shoreline. where Q Q  KA h1 h2 , L (1)

and

K h1 h2

is flow through a vertical plane that extends beneath the shoreline of a surface-water body (L3/T), is the area of the plane through which all water must pass to either originate from the surface-water body or end up in the surface-water body, depending on the direction of flow (L2), is horizontal hydraulic conductivity (L/T), is hydraulic head at the upgradient well (L), is hydraulic head at the shoreline of the surface-water body (L), is distance from the well to the shoreline (L).

Introduction and Characteristics of Flow 21

EXPLANATION
8 Well location and I3 04 24

number Area of inseepage Area of outseepage

05
Mary Lake
I3
04
03

Williams Lake

31

02
4 I2 8 7

01

9500'

I1

4700'

Williams Lake

25

MINNESOTA

0 0

100 300

200 300 METERS 600 900 1200 FEET

Figure 16. Example of shoreline segment definition for the calculation of water fluxes between surface water and ground water at a lake (modified from LaBaugh and others, 1995, used in accordance with author rights of the National Research Council of Canada Press). Location of wells and shoreline segments used to calculate flow between surface water and ground water at Williams Lake, Minnesota.

1996), as well as creeks and rivers flowing into Long Island Sound, New York (Pluhowski, 1972). Examples of the use of infrared imagery to detect areas of ground-water discharge to marine waters include the delineation of areas of diffuse ground-water flow into Long Island Sound (Pluhowski, 1972) and focused ground-water flow as springs to the ocean, such as around the perimeter of the island of Hawaii (Fischer and others, 1966).

Dye and Tracer Tests

Dyes and other soluble tracers can be added to water and then tracked to provide direct, qualitative information about ground-water movement to streams. Fluorescent dyes that are readily detected at small concentrations and pose little environmental risk make a useful tool for tracing ground-water flow paths, particularly in karst terrane (Aley and Fletcher, 1976; Smart and Laidlaw, 1977; Jones, 1984; Mull and others,

1988). Thus, dye-tracer studies can be used to determine the time-of-travel for ground water to move to and into surface water, as well as hydraulic properties of aquifer systems (Mull and others, 1988). The use of dyes as tracers is described in more detail in Chapter 3. Commonly, a reconnaissance of the ground-water basin is made to identify likely areas of potential surface-water flow into ground water or ground-water flow to the surface. An inventory is made of springs, sinkholes, boreholes or screened wells, and sinking streams. Appropriate sites then are picked for dye injection, and the potential discharge areas, springs, and stream reaches are monitored over an appropriate period of time, hours or days, for appearance of the dye (fig. 18). Solute tracers have been used to aid in the determination of water gains or losses within the channel of a stream or river (Kilpatrick and Cobb, 1985). This technique is known as dilution-gaging. A variety of tracers have been used in such studies, either alone or in combination, usually including a

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

Figure 17. Example of use of thermal infrared imagery to delineate areas of discrete and diffuse ground-water discharge to surface water (from Banks and others, 1996, reprinted from Ground Water with permission from the National Ground Water Association, copyright 1996, thermal imagery of O-Field study area, Aberdeen Proving Ground, Maryland).

Introduction and Characteristics of Flow 23

Figure 18. Example of the use of dye to examine water fluxes between surface water and ground water (from Carter and others, 2002). Dye testing has been done in Boxelder Creek, South Dakota, which can lose as much as 50 cubic feet per second of flow to the bedrock aquifers. In the upper left photograph, nontoxic, red dye is poured into Boxelder Creek upstream from a major loss zone. In the upper right photograph, dye in the stream can be seen disappearing into a sinkhole in the Madison Limestone. In the bottom photograph, dye in the stream emerges downstream at Gravel Spring, which is about 671 meters (2,200 feet) (linear distance) from the major loss zone. The length of time for the first arrival of dye to travel this distance is variable depending on flow conditions but generally is about 1 to 2 hours (Strobel and others, 2000). Thus, the ground-water velocity is about 0.3 to 0.6 kilometer per hour (0.2 to 0.4 mile per hour), which is a very fast rate for ground water. Dye also has been recovered at City Springs, which is in the Rapid Creek Basin, about 30 days after injection. This demonstrates that ground-water flow paths are not necessarily restricted by surface-water drainage basins. (Photographs by Derric L. Isles, South Dakota Department of Environment and Natural Resources.)

24 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water tracer expected to be nonreactive in the waters of the stream to which it is added, such as lithium (for example, in the Snake River, ColoradoBencala and others, 1990) or chloride (for example, Chalk Creek, ColoradoKimball, 1997). In this technique, a known quantity of solute is added at a specified rate for a short interval of time at the upstream cross section of the stream segment of interest, and concentrations of the solute are measured at one or more points downstream over time. Discharge is calculated from the amount of dilution that occurs at the downstream point or points. Kimball (1997) indicated that Qs (the discharge in the stream) is calculated as follows: where and CA is the tracer concentration upstream from the injection point. Qs = (CiQi) / (CB CA) Ci Qi CB is the tracer concentration in the injection solution, is the rate of injection into the stream, is the tracer concentration downstream, (2) where and R is residual. Chemical budgets have been calculated in lake and wetland studies where water exchange between surface-water bodies and underlying ground water was of interest [see, for example, Rawson Lake, OntarioSchindler and others (1976); Thoreaus Bog, MassachusettsHemond (1983); Williams Lake, MinnesotaLaBaugh and others (1995); LaBaugh and others (1997); multiple lakes in Polk and Highlands Counties, FloridaSacks and others (1998); Lake Kinneret, IsraelRimmer and Gideon (2003)]. The equation can be modified to solve for any unknown flow term, provided that the remainder of the flow terms are known. The accuracy of the method depends greatly on the accuracy of the other flow and chemical-concentration measurements. The size of the residual term often is considered a general indicator of the accuracy of the method, but a small residual does not always indicate an accurate chemical balance. LaBaugh (1985) and Choi and Harvey (2000) provide examples of the use of error analysis to quantify the uncertainty associated with water-flux results obtained using this method. The ratios of the isotopes of oxygen and hydrogen present in water have been used for decades to distinguish sources of water, including ground-water discharge to surfacewater bodies (for example, Dincer, 1968). These isotopes are useful because they are part of the water and not solutes dissolved in the water. The method works well when the degree of isotopic fractionation of the water is different for different sources of water. The process of evaporation tends to remove lighter isotopes, leaving the heavier isotopes behind. Thus, the ratio of lighter to heavier isotopes will change over time in the water and the water vapor. More detailed explanation of the isotopic fractionation in catchment water is given in Kendall and others (1995). If the isotopic compositions of different sources of water are distinct, then simple mixing models can be used to identify sources of water. A brief example is presented here, but more detailed explanations and examples of the use of this method can be found in Krabbenhoft and others (1994), LaBaugh and others (1997), Sacks (2002) (applied

P GWI SI E GWO SO VL Cx

is precipitation, is ground-water flux into lake or wetland, is streamflow into lake or wetland, is evaporation, is flux of lake or wetland water to ground water, is streamflow out of lake or wetland, is change in lake or wetland volume, is chemical concentration of hydrologic component,

When coupled with stream-segment discharge measurements, use of solute tracers also enables calculation of the rates of ground-water inflow and outflow within a stream segment (Harvey and Wagner, 2000). At the same time the solute is injected and monitored within the stream segment, physical velocity measurements (streamflow) are made at the upstream and downstream sections of the stream reach. The streamflow measurements provide information on whether or not there was a net loss or gain of flow within the reach due to interaction with ground water. Harvey and Wagner (2000) indicate the solute tracer, or dilution-gaging, values determine groundwater inflow. Thus, ground-water outflow can be calculated by subtracting the net loss or gain from the solute tracer-derived ground-water inflow value. Calculation of chemical budgets for a stream, lake, or wetland is another way in which solutes can be used to make quantitative estimates of surface-water exchange with ground water. Conservative chemicals in a watershed are those that are not altered by the porous media through which they flow, and occur at concentrations for which changes in concentration because of chemical precipitation are not likely to occur. Conservative chemicals can be used to determine the volume of ground water that flows into or out of a surface-water body, provided that all other fluxes are known. A common form of the chemical-budgeting equation for a lake or wetland is P(CP) + GWI(CGWI) + SI(CSI) GWO(CGWO) SO(CSO) = VL(CL) R, (3)

Introduction and Characteristics of Flow 25 to lakes), Kendall and others (1995) (a brief description of the methods), and Kendall and McDonnell (1998) (detailed descriptions of numerous isotopic methods). For determination of ground-water discharge to streams and rivers, a simple two-component mixing model often is used: where Q QSS = QGWGW + QPP is discharge, is the stable-isotopic composition in parts per thousand enrichment or depletion (per mil) relative to a standard, is stream water, is ground water, (4) atmosphere will be lost from that water because of volatilization. Thus, ground water commonly contains higher activities of 222Rn than does surface water, from 3 to 4 orders of magnitude greater (Burnett and others, 2001). 222Rn is a radioactive daughter isotope of radium 226 (226Ra) and has a half life of 3.82 days. Determination of the activity of 222Rn and 226Ra in surface water enables the calculation of the 222Rn excesshow much more 222Rn is present in surface water than would be expected based on the 226Ra content of the water. Determination of the activity of 222Rn and 226Ra in sediment water or ground water is used to determine the 222Rn flux into surface water (Cable and others, 1996; Corbett and others, 1998), which can account for the excess 222Rn in the surface water. The mass balance or flux of 222Rn has been used to determine ground-water discharge to several types of surface-water bodies (Kraemer and Genereaux, 1998): in streams, such as the Bickford watershed, Massachusetts (Genereaux and Hemond, 1990); rivers, such as the Rio Grande de Manati, Puerto Rico (Ellins and others, 1990); lakes, such as Lake Kinneret, Israel (Kolodny and others, 1999); estuaries, such as Chesapeake Bay (Hussain and others, 1999), Charlotte Harbor, Florida (Miller and others, 1990), and Florida Bay (Corbett and others, 1999), as well as the coastal ocean, such as in Kanaha Bay, Oahu, Hawaii (Garrison and others, 2003); the Gulf of Mexico off of Florida (Cable and others, 1996; Burnett and others, 2001); and the Atlantic Ocean off the coast of South Carolina (Corbett and others, 1998). In the radium isotope method, the surface-water activities of the four naturally occurring radium isotopes226Ra, 228 Ra, 223Ra, and 224Raare compared to activities in sediment water or ground water to determine fluxes (Moore, 2000). The source of the radium isotopes is the decay of uranium and thorium in sediments or rocks. Water in contact with solid materials containing the source of the isotopes will accumulate the isotopes. Ground water will accumulate more of the isotopes because of waters presence within the matrix of the sediments or rocks. Surface waters will accumulate less of the isotopes because the sediments or rocks are less abundant relative to the water (Kraemer, 2005). Kraemer indicates the ratio of the longer lived isotopes (226Ra half-life of 1,601 years, 228 Ra half-life of 5.8 years) can be used as an indicator of the types of sediments or rocks through which ground water has traveled, because of differences in uranium and thorium content between rock types. Kraemer (2005) also notes that the short-lived isotopes (223Ra half life of 11.4 days, and 224Ra half-life of 3.7 days) provide some indication of the timing of ground-water discharge. Naturally occurring radium isotopes have been useful in the identification of ground-water inflow to lakes, such as Cayuga Lake, New York (Kraemer, 2005), freshwater wetlands in the Florida Everglades (Krest and Harvey, 2003), estuarine wetlands, North Inlet salt marsh, South Carolina (Krest and others, 2000), and coastal waters in the central South Atlantic Bight (Moore, 2000).

S GW and P is precipitation. For lakes and wetlands, where sources of water are more numerous, slightly more complex mixing models can be used, such as those provided by Krabbenhoft and others (1990): GWI  P (D L D P )  E (D E D L ) D GWI D L (5)

or that provided by Krabbenhoft and others (1994): GWO = P P VL L GWI GWI , L (6)

where GWI is ground-water flux into lake, GWO is flux of water from lake to ground water, P is precipitation, E is evaporation, VL is change in the volume of the lake, and X is per mil value for hydrologic component. Where equations 5 and 6 are derived from equation 3 applied to stable isotopes at steady state: d(VL)/dt = GWI GWI + P P + Si Si GWO L E E So L = 0. (7)

Investigation of ground-water discharge into inland and marine surface water also is feasible through measurement of radon and radium isotopes (Corbett and others, 1998; Moore, 2000). In the radon isotope method, a mass balance is constructed for radon-222 (222Rn), which is a chemically and biologically inert radioactive gas formed by the disintegration of the parent nuclide radium (Corbett and others, 1998, 1999). Because radon is a gas, radon in water in contact with the

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

Thermal Profiling
Measurements of temperature have been used to determine qualitatively the locations of rapid discharge of ground water to surface water; common measurement methods include towing a tethered temperature probe, in-situ measurements of temperature, and thermal imagery (Lee, 1985; Baskin, 1998; Rosenberry and others, 2000). Temperature also can be measured at different depths beneath a stream or lake to determine the rate of vertical flow through a surface-water bed either into or out of the surface-water body (Lapham, 1989). This method is effective when flow through the lakebed is sufficient to allow advective processes to be significant relative to conductive temperature signals. The method requires multiple measurements of temperature over weeks or months and the simultaneous solution of the flow of fluid and heat in one dimension. Depths of temperature measurement typically extend up to 3 to 6 meters beneath the sediment bed. The solution requires the assumption that flow is vertical through the surface-water bed and that the media are homogeneous and isotropic. Taniguchi (1993) used seasonal changes in sediment temperature beneath a surface-water body to develop type curves that can be used to estimate vertical fluxes through the surface-water bed. Subsediment temperature has been used over short distances beneath the surface-water bed and makes use of the temperature response in the sediments to diurnal changes in surface-water temperature (Constantz and others, 1994; Stonestrom and Constantz, 2003). This method adjusts parameters in a one-dimensional, variably saturated heattransport model (VS2DH) [developed in the USGS (Healy and Ronan, 1996)] until the simulation results match the temperature data collected beneath the surface-water body. Sediment temperature measurements and data also have been used to determine streamflow frequency in ephemeral stream channels (Constantz and others, 2001; Stonestrom and Constantz, 2003). This method is described in greater detail in Chapter 4. Conant (2004) used measurements of hydraulic-head gradient and hydraulic conductivity in wells installed in a streambed to determine rates of exchange between ground water and surface water. Conant then developed an empirical relation between streambed flux and streambed temperature relative to stream-water temperature
Figure 19. Example of the use of a towed specific-conductance probe to identify ground-water discharge to surface waters (from Harvey and others, 1997, used in accordance with usage permissions of the American Geophysical Union). The sediment probe is being towed behind a small boat. The probe is used to detect areas of more electrically conductive ground-water inflow. In saline waters, the probe is used to detect less electrically conductive fresh ground-water inflow.

and used the empirical relation to indicate rate of discharge of ground water to the stream at locations where only streambed temperature was measured.

Use of Specific-Conductance Probes

Specific-conductance probes are another tool that can be useful for locating areas of ground-water discharge to surface waters (Lee, 1985; Vanek and Lee, 1991; Harvey and others, 1997) (fig. 19). Such probes are suspended from a boat with a cable connecting the probes to a specificconductance meter on board the boat. The housing for the electrically conductive probes is designed to maintain contact with the sediments (Lee, 1985) so that the probes are dragged through bottom sediments at a depth of 1 to 3 centimeters (Vanek and Lee, 1991). This method depends on having the existence of a difference in the electrical conductance of surface and ground water great enough to be detected by the sensors (probes) and thus to identify points or areas of ground-water discharge to the surface water. Thus, saline waters receiving ground-water discharge of less salinity (Vanek and Lee, 1991) or fresh surface waters receiving more mineralized ground-water discharge (Harvey and others, 1997) are environments where this technique is effective. Changes in electrical conductance, however, also may reflect changes in sediment type so that the technique should be considered only as a reconnaissance tool. The identification of places where ground-water discharge may be occurring should be verified by other, complementary techniques (Vanek and Lee, 1991). Specific-conductance and temperature probes have been used in this way along many kilometers of river reaches in combination with Global Positioning System information to determine the precise locations of groundwater discharge and their variation with season (Vaccaro and Maloy, 2006).

Conductance

Distance

Introduction and Characteristics of Flow 27

Electrical Resistivity Profiling

Electrical geophysical survey methods are applicable to delineating submarine ground-water discharge (Loke and Lane, 2004). Salinity differences in water can be detected by the resistance to the passage of current between electrodes. Fresh or weakly saline waters will be more resistant to movement of an electric current than more saline waters. Two resistivity data acquisition geometries commonly are used: (1) continuous resistivity profiling (Belaval and others, 2003; Manheim and others, 2004; Day-Lewis and others, 2006), and (2) marine resistivity (Taniguchi and others, 2006). In continuous resistivity profiling, a streamer comprising a set of floating electrodes spaced at a regular interval is towed on the water surface behind a boat. As the boat travels along a transect, electrical current is applied at a fixed time interval at one or more electrode pairs, and electrical potentials are measured simultaneously between other electrode pairs. At the same time, water depths along the transect are measured with echo sounding. In marine resistivity profiling, measurements are made using electrodes placed on the water bottom. Regardless A

of the acquisition geometry, data are inverted to produce two-dimensional cross sections, or tomograms, of subsurface resistivity.Such data can indicate locations within a transect where submarine ground-water discharge is occurring, as well as delineating the subsurface saltwater/freshwater interface or geologic structure. Tomograms are commonly interpreted in the context of direct, discrete measurements of conductivity, temperature, and depth.

Hydraulic Potentiomanometer (Portable Wells) Measurements

Many devices have been designed to be installed through the bed of a surface-water body to measure the vertical hydraulic-head gradient beneath the surface-water body. These devices provide a direct measurement of hydraulic head relative to surface-water stage at the depth to which the probe is inserted beneath the surface-water bed. One of the most commonly used devices of this type is the hydraulic potentiomanometer (fig. 20), also sometimes called a minipiezometer B
Valve Tube C Pump

To optional sample bottle Tube A Manometer board Probe Water surface

Tube B

Valve or clamp

Filter

Lake bed

Figure 20. A, Hydraulic potentiomanometer with users demonstrating hydraulic-head difference between surface water and ground water. The hydraulic head of the river (right side of the board) is higher than the hydraulic head of the ground water (left side of the board), indicating a downward gradient (from Simonds and Sinclair, 2002). Photograph by Kirk A. Sinclair, Washington State Department of Ecology, 2002. B, Components of a hydraulic-potentiomanometer system (from Winter and others, 1988, copyright 1988 by the American Society of Limnology and Oceanography, Inc., used with permission).

28 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water (Winter and others, 1988). Although this device does not provide direct measurements of flow across the sediment-water interface, it is useful for making qualitative determinations of direction of flow across the sediment-water interface. Estimates of hydraulic conductivity also can be made on the basis of the amount of suction required to pull water through the screen at the end of the probe. These devices often are used in conjunction with a network of wells to obtain additional hydraulic-head data along shoreline segments where wells are not or cannot be situated. Details about use of these devices is provided in Chapter 2.
Water surface

Bag

Seepage cylinder Water Sediment

Seepage Meter Measurements

Seepage meters are devices that isolate a small area of the bed of a surface-water body and measure the flow of water across that area. References as early as 1944 indicate that early seepage meters were designed to measure water loss through unlined canals (Carr and Winter, 1980). Beginning in the early 1970s, seepage meters have been used in lakes, rivers, wetlands, and estuaries to measure flows between surface water and ground water in natural settings. Use in coastal environments has increased during recent years (Cable and others, 1997). One of the most common devices, called a half-barrel seepage meter, uses a cut-off end of a steel (or plastic) storage drum to isolate a small circular area of the surface-water bed, and a plastic bag is attached to the barrel to register the change in water volume over the time of bag attachment (Lee, 1977; Lee and Cherry, 1978) (fig. 21). More details about the use and interpretation of data from seepage meters are included in Chapter 2. Many device modifications have been made to the basic design of seepage meters for use in deep water, soft sediments, shallow water, and areas exposed to large waves. Many investigations of surface-water and ground-water interaction require integrating point measurements of water flux in order to interpret the total flow between surface water and ground water. The most common method is to average point measurements and apply that average value to all or part of the surface-water body of interest. This method is not appropriate, however, where ground-water flux into a surface water declines with greater distance from shore. If data are collected along transects perpendicular to the shore, a curve can be fit through the data. Because ground-water flow to surface water commonly is distributed exponentially with distance from shore, exponential curves have been fit through transect data, and the equation for the curve has been used to calculate water flux for an area deemed to be representative of the transect point measurements (for example, Fellows and Brezonik, 1980). Where ground-water flux is not exponentially distributed with distance from shore, a plot of point measurements of seepage flux with distance from shore is made, and the area under the curve can be determined to represent total seepage for a unit width of shoreline. That value can

Water surface

Bag

Water Sediment

Figure 21. Full-section view of seepage meter showing details of placement in the sediment (modified from Lee and Cherry, 1978, used with permission of the Journal of Geoscience Education).

then be multiplied by an appropriate shoreline length to determine a flux volume for a specific area of the surface-water body. If point measurements are made at a sufficient number of transects, ground-water flux for an area of a surface-water body can be determined.

Biological Indicators
The biological response to conditions of flow at the sediment-water interface can be an indicator of the direction and relative magnitude of flow. The growing field of ground-water ecology has made frequent use of distributions of specific types of plants and animals as an indicator of ground-water/surface-water interaction (Danielopol, 1984; Danielopol and others, 1997; Lodge and others, 1989; Malard and others, 1996; Goslee and others, 1997; Wetzel, 1999) (fig.22). These methods are useful reconnaissance tools to aid in locating areas in need of more detailed investigations. Typically, these methods involve identifying species or groups of species of plants or animals that are known to thrive in places where ground water discharges to surface water, but some of the indicators also indicate areas where surface water flows into ground water. Identification of specific plant and animal species is necessary for use of these methods, but some of the species are readily identifiable without requiring intensive biological or ecological training (Rosenberry and others, 2000).

Introduction and Characteristics of Flow 29

Figure 22. Marsh marigold in Shingobee Lake, Minnesota. Presence indicates location of ground-water discharge to the lake. (Photographs by Donald O. Rosenberry, U.S. Geological Survey.)

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

References
Adamowski, K., and Feluch, W., 1991, Application of nonparametric regression to groundwater level prediction: Canadian Journal of Civil Engineering, v. 18, no. 4, p.600606. Aley, T.J., and Fletcher, M.W., 1976, The water tracers cookbookA guide to the art and science of water tracing materials with particular emphasis on the use of fluorescent dyes, Lycopodium spores, and bacteriophage in ground water investigations: Missouri Speleology, v. 16, no. 3, 32 p. Alley, W.M., Reilly, T.E., and Franke, O.L., 1999, Sustainability of ground-water resources: U.S. Geological Survey Circular 1186, 79 p. Andrews, C.B., and Anderson, M.P., 1978, Impact of a power plant on the ground-water system of a wetland: Ground Water, v. 16, no. 2, p. 105111. Asbury, C.E., 1990, The role of groundwater seepage in sediment chemistry and nutrient budgets in Mirror Lake, New Hampshire: Ithaca, Cornell University, Ph.D. dissertation, 275 p. Attanayake, M.P., and Waller, D.H., 1988, Use of seepage meters in a groundwater-lake interaction study in a fractured rock basinA case study: Canadian Journal of Civil Engineering, v. 15, p. 984989. Bacchus, S.T., 2001, Knowledge of groundwater responsesA critical factor in saving Floridas threatened and endangered species, Part 1Marine ecological disturbances: Endangered Species UPDATE, v. 18, no. 3, p. 7990. Bacchus, S.T., 2002, The ostrich component of the multiple stressor modelUndermining South Florida, in Porter, J.W., and Porter, K.G., eds., The Everglades, Florida Bay, and Coral Reefs of the Florida KeysAn Ecosystem Sourcebook: Boca Raton, CRC Press, p. 677748. Banks, W.S.L., Paylor, R.L., and Hughes, W.B., 1996, Using thermal-infrared imagery to delineate ground-water discharge: Ground Water, v. 34, no. 3, p. 434443. Barlow, P.M., and Dickerman, D.C., 2001, Numerical simulation and conjunctive-management models of the Hunt-Annaquatucket-Pettaquamscutt stream-aquifer system, Rhode Island: U.S. Geological Survey Professional Paper1636, 88p., one map sheet. Barlow, P.M., and Wild, E.C., 2002, Bibliography on the occurrence and intrusion of saltwater in aquifers along the Atlantic Coast of the United States: U.S. Geological Survey Open-File Report 2002235, 30 p.

Baskin, R.L., 1998, Locating shoreline submarine springs in two Utah lakes using thermal imagery, in Pitman, J.K., and Carroll, A.R., eds., Modern & ancient lake systemsNew problems and perspectives: Salt Lake City, Utah Geological Association, v. 26, p. 5157. Belanger, T.V., and Kirkner, R.A., 1994, Groundwater/surface water interaction in a Florida augmentation lake: Lake and Reservoir Management, v. 8, no. 2, p. 165174. Belaval, Marcel, Lane, J.W., Jr., Lesmes, D.P., and Kineke, G.C., 2003, Continuous-resistivity profiling for coastal groundwater investigationsThree case studies, in Symposium on the Application of Geophysics to Engineering and Environmental Problems (SAGEEP), April 610, 2003, SanAntonio, Texas, Proceedings: Denver, Colorado, Environmental and Engineering Society, CDROM, 14 p. Bencala, K.E., McKnight, D.M., and Zellweger, G.W., 1990, Characterization of transport in an acidic and metal-rich mountain stream based on a lithium tracer injection and simulations of transient storage: Water Resources Research, v. 26, no. 5, p. 9891000. Beven, K.J., ed., 1997, Distributed hydrological modelling Application of the TOPMODEL concept: Chichester, Wiley, 348 p. Beven, K.J., and Feyen, J., 2002, The future of distributed modeling: Hydrological Processes, v. 1, p. 169172. Beven, K.J., Kirkby, M.J., Schofield, N., and Tagg, A.F., 1984, Testing a physically based flood-forecasting model (TOPMODEL) for three U.K. catchments: Journal of Hydrology, v. 69, p. 119143. Bokuniewicz, Henry, 1980, Groundwater seepage into Great South Bay, New York: Estuarine and Coastal Marine Science, v. 10, p. 437444. Bouwer, H., and Maddock, T., III, 1997, Making sense of the interactions between groundwater and streamflowLessons for water masters and adjudicators: Rivers, v. 6, no. 1, p.1931. Buchtele, J., Buchtelov, M., Gallart, F., Latron, J., Llorens, P., Salvany, C., and Herrmann, A., 1998, Rainfall-runoff processes modelling using SACRAMENTO and BROOK models in the Cal Rod catchment (Pyrenees, Spain), Conference on Catchment Hydrological and Biochemical Processes in Changing Environment: Liblice, Czechoslovakia, p. 1316. Burnett, W.C., Kim, G., and Lane-Smith, D., 2001, A continuous monitor for assessment of 222Rn in the coastal ocean: Journal of Radioanalytical and Nuclear Chemistry, v. 249, no. 1, p. 167172.

Introduction and Characteristics of Flow 31 Cable, J.E., Bugna, G.C., Burnett, W.C., and Chanton, J.P., 1996, Application of 222Rn and CH4 for assessment of groundwater discharge to the coastal ocean: Limnology and Oceanography, v. 41, no. 6, p. 13471353. Cable, J.E., Burnett, W.C., Chanton, J.P., Corbett, D.R., and Cable, P.H., 1997, Field evaluation of seepage meters in the coastal marine environment: Estuarine, Coastal and Shelf Science, v. 45, no. 3, p. 367375. Cambareri, T.C., and Eichner, E.M., 1998, Watershed delineation and ground-water discharge to a coastal embayment: Ground Water, v. 36, no. 4, p. 626634. Campbell, C.W., and Keith, A.G., 2001, Karst ground-water hydrologic analyses based on aerial thermography: Hydrological Science and Technology, Special Issue, v. 17, no.14, p. 5968. Carr, M.R., and Winter, T.C., 1980, An annotated bibliography of devices developed for direct measurement of seepage: U.S. Geological Survey Open-File Report 80344, 38 p. Carter, J.M., Driscoll, D.G., Williamson, J.E., and Lindquist, V.A., 2002, Atlas of water resources in the Black Hills area, South Dakota: U.S. Geological Survey Hydrologic Investigations Atlas HA747, 120 p. Cheng, Xiangxue, and Anderson, M.P., 1993, Numerical simulation of ground-water interaction with lakes allowing for fluctuating lake levels: Ground Water, v. 31, no. 6, p.929933. Childress, C.J.O., Sheets, R.A., and Bair, E.S., 1991, Hydrology and water quality near the South Well Field, southern Franklin County, Ohio, with emphasis on the simulation of ground-water flow and transport of Scioto River: U.S. Geological Survey Water-Resources Investigations Report914080, 78 p. Choi, Jungyill, and Harvey, J.W., 2000, Quantifying timevarying ground-water discharge and recharge in wetlands of the northern Florida Everglades: Wetlands, v. 20, no. 3, p.500511. Committee on Hydrologic Science, National Reseach Council, 2004, Groundwater fluxes across interfaces: Washington, D.C., The National Academies Press, 85 p. Conant, Brewster, Jr., 2004, Delineating and quantifying ground-water discharge zones using streambed temperatures: Ground Water, v. 42, no. 2, p. 243257. Constantz, James, Stonestrom, David, Stewart, A.E., Niswonger, Rich, and Smith, T.R., 2001, Analysis of streambed temperatures in ephemeral channels to determine streamflow frequency and duration: Water Resources Research, v. 37, no. 2, p. 317328. Constantz, James, Thomas, C.L., and Zellweger, G., 1994, Influence of diurnal variations in stream temperature on streamflow loss and groundwater recharge: Water Resources Research, v. 30, no. 12, p. 32533264. Corbett, D.R., Burnett, W.C., Cable, P.H., and Clark, S.B., 1998, A multiple approach to the determination of radon fluxes from sediments: Journal of Radioanalytical and Nuclear Chemistry, v. 236, p. 247252. Corbett, D.R., Chanton, Jeffrey, Burnett, William, Dillon, Kevin, Rutkowski, Christine, and Fourqurean, J.W., 1999, Patterns of groundwater discharge into Florida Bay: Limnology and Oceanography, v. 44, no. 4, p. 10451055. Danielopol, D.L., 1984, Ecological investigations on the alluvial sediments of the Danube in the Vienna areaA phreatobiological project: Verhandlungen Internationale Vereiningen Limnologie, v. 22, p. 17551761. Danielopol, D.L., Rouch, R., Pospisil, P., Torreiter, P., and Moeszlacher, F., 1997, Ecotonal animal assemblages; their interest for groundwater studies, in Gibert, J., Mathieu, J., and Fournier, F., eds., Groundwater/surface water ecotonesBiological and hydrological interactions and management options: Cambridge, Cambridge University Press, p. 1120. Danskin, W.R., 1998, Evaluation of the hydrologic system and selected water-management alternatives in the Owens Valley, California: U.S. Geological Survey Water-Supply Paper2370H, 175 p., 6 pls. in pocket. Day-Lewis, F.D., White, E.A., Belaval, M., Johnson, C.D., and Lane, J.W., Jr., 2006, Continuous resistivity profling to delineate submarine ground-water dischargeExamples and limitations: The Leading Edge, v. 25, no. 6, p.724728. Dincer, T., 1968, The use of oxygen 18 and deuterium concentrations in the water balance of lakes: Water Resources Research, v. 4, no. 6, p. 12891306. Dokulil, M.T., Teubner, K., and Donabaum, K., 2000, Restoration of a shallow, ground-water fed urban lake using a combination of internal management strategiesA case study: Advances in Limnology, v. 55, p. 271282. Donato, M.M., 1998, Surface-water/ground-water relations in the Lemhi River Basin, East-Central Idaho: U.S. Geological Survey Water-Resources Investigations Report 984185, 29p. Doss, P.K., 1993, Nature of a dynamic water table in a system of non-tidal freshwater coastal wetlands: Journal of Hydrology, v. 141, p. 107126. Dumouchelle, D.H., 2001, Evaluation of ground-water/ surface-water relations, Chapman Creek, west-central, Ohio, by means of multiple methods: U.S. Geological Survey Water-Resources Investigations Report 20014202, 13 p.

32 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water Ellins, K.K., Roman-Mas, A., and Lee, R., 1990, Using 222Rn to examine groundwater/surface discharge interaction in the Rio Grande de Manati, Puerto Rico: Journal of Hydrology, v. 115, p. 319341. Erickson, D.R., 1981, A study of littoral groundwater seepage at Williams Lake, Minnesota, using seepage meters and wells: Minneapolis, University of Minnesota, M.S. thesis, 135 p. Federer, C.A., and Lash, D., 1978, BROOKA hydrologic simulation model for eastern forests: Water Resources Research Center, University of New Hampshire Research Report no. 19, 84 p. Fellows, C.R., and Brezonik, P.L., 1980, Seepage flow into Florida lakes: Water Resources Bulletin, v. 16, no. 4, p.635641. Fenske, J.P., Leake, S.A., and Prudic, D.E., 1996, Documentation of a computer program (RES1) to simulate leakage from reservoirs using the modular finite-difference groundwater flow model (MODFLOW): U.S. Geological Survey Open-File Report 96364, 51 p. Fischer, W.A., Davis, D.A., and Sousa, T.M., 1966, Fresh-water springs of Hawaii from infrared images: U.S. Geological Survey Hydrologic Atlas 218, 1 map. Fraser, C.J.D., Roulet, N.T., and Lafleur, M., 2001, Ground water flow pattern in a large peatland: Journal of Hydrology, v. 246, p. 142154. Galloway, D.L., Alley, W.M., Barlow, P.M., Reilly, T.E., and Tucci, P., 2003, Evolving issues and practices in managing ground-water resourcesCase studies on the role of science: U.S. Geological Survey Circular 1247, 73 p. Garrett, J.W., Bennett, D.H., Frost, F.O., and Thurow, R.F., 1998, Enhanced incubation success for Kokanee spawning in groundwater upwelling sites in a small Idaho stream: North American Journal of Fisheries Management, v. 18, no. 4, p. 925930. Garrison, G.H., Glenn, C.R., and McMurtry, G.M., 2003, Measurement of submarine groundwater discharge in Kahana Bay, Oahu, Hawaii: Limnology and Oceanography, v. 48, no. 2, p. 920928. Genereux, D.P., and Hemond, H.F., 1990, Naturally occurring radon 222 as a tracer for streamflow generationSteadystate methodology and field application: Water Resources Research, v. 26, no. 12, p. 30653075. Glennon, Robert, 2002, Water folliesGroundwater pumping and the fate of Americas fresh waters: Washington, D.C., Island Press, 314 p. Goslee, S.C., Brooks, R.P., and Cole, C.A., 1997, Plants as indicators of wetland water source: Plant Ecology, v. 131, p.199206. Harbaugh, A.W., Banta, E.R., Hill, M.C., and McDonald, M.G., 2000, MODFLOW-2000, the U.S. Geological Survey modular ground-water modelUser guide to modularization concepts and the ground-water flow process: U.S.Geological Survey Open-File Report 200092, 121 p. Harte, P.T., Flynn, R.J., Kiah, R.G., Severance, T., and Coakley, M.F., 1997, Information on hydrologic and physical properties of water to assess transient hydrology of the MilfordSouhegan glacial-drift aquifer, Milford, New Hampshire: U.S. Geological Survey Open-File Report97414, 45 p. Harvey, F.E., Lee, D.R., Rudolph, D.L., and Frape, S.K., 1997, Locating groundwater discharge in large lakes using bottom sediment electrical conductivity mapping: Water Resources Research, v. 33, no. 11, p. 26092615. Harvey, J.W., and Wagner, B.J., 2000, Quantifying hydrologic interactions between streams and their subsurface hyporheic zones, in Jones, J.B., and Mulholland, P.J., eds., Streams and ground waters: San Diego, Academic Press, p. 344. Healy, R.W., and Ronan, A.D., 1996, Documentation of computer program VS2DH for simulation of energy transport in variably saturated porous mediaModification of the U.S. Geological Survey's computer program VS2DT: U.S. Geological Survey Water-Resources Investigations Report964230, 36 p. Hemond, H.F., 1983, The nitrogen budget of Thoreaus Bog: Ecology, v. 64, p. 94109. Herbert, L.R., and Thomas, B.K., 1992, Seepage study of the Bear River including Cutler Reservoir in Cache Valley, Utah and Idaho: Salt Lake City, Utah Department of Natural Resources, Division of Water Rights, State of Utah Department of Natural Resources Technical Publication no. 105, 18 p. Hill, M.C., Lennon, G.P., Brown, G.A., Hebson, C.S., and Rheaume, S.J., 1992, Geohydrology of simulation of ground-water flow in the valley-fill deposits in the Ramapo River Valley, New Jersey: U.S. Geological Survey WaterResources Investigations Report 904151, 92 p. Hussain, N., Church, T.M., and Kim, G., 1999, Use of 222Rn and 226Ra to trace groundwater discharge into the Chesapeake Bay: Marine Chemistry, v. 65, p. 127134. Hutson, S.S., Barber, N.L., Kenny, J.F., Linsey, K.S., Lumia, D.S., and Maupin, M.A., 2004, Estimated use of water in the United States in 2000: U.S. Geological Survey Circular1268, 52 p.

Introduction and Characteristics of Flow 33 Jacobs, K.L., and Holway, J.M., 2004, Managing for sustainability in an arid climateLessons learned from 20 years of groundwater management in Arizona, USA: Hydrogeology Journal, v. 12, no. 1, p. 5265. Jaquet, N.G., 1976, Ground-water and surface-water relationships in the glacial province of northern WisconsinSnake Lake: Ground Water, v. 14, no. 2, p. 194199. Jobson, H.E., and Harbaugh, A.W., 1999, Modifications to the diffusion analogy surface-water flow model (DAFLOW) for coupling to the modular finite-difference ground-water flow model (MODFLOW): U.S. Geological Survey Open-File Report 99217, 107 p. Jones, W.K., ed., 1984, Water tracing special issue: National Speleological Society Bulletin, v. 46, no. 2, 48 p. Johannes, R.E., 1980, The ecological significance of the submarine discharge of groundwater: Marine Ecology Progress Series, v. 3, p. 365373. Johannes, R.E., and Hearn, C.J., 1985, The effect of submarine groundwater discharge on nutrient and salinity regimes in a coastal lagoon off Perth, Western Australia: Estuarine, Coastal and Shelf Science, v. 21, p. 789800. Kaleris, V., 1998, Quantifying the exchange rate between groundwater and small streams: Journal of Hydraulic Research, v. 36, no. 6, p. 913932. Kilpatrick, F.A., and Cobb, E.D., 1985, Measurement of discharge using tracers: U.S. Geological Survey, Techniques of Water-Resources Investigations, book 3, chap. A16, 73 p. Kendall, Carol, and McDonnell, J.J., eds., 1998, Isotope tracers in catchment hydrology: New York, Elsevier Science Publishers, 839 p. Kendall, Carol, Sklash, M.G., and Bullen, T.K., 1995, Isotope tracers of water and solute sources in catchments, in Trudgill, S.T., ed., Solute modelling in catchment systems: New York, John Wiley & Sons, Ltd., p. 261303. Kimball, Briant, 1997, Use of tracer injections and synoptic sampling to measure metal loading from acid mine drainage: U.S. Geological Survey Fact Sheet 24596, 4 p. Konrad, C.P., Drost, B.W., and Wagner, B.J., 2003, Hydrogeology of the unconsolidated sediments, water quality, and ground-water/surface-water exchanges in the Methow River Basin, Okanogan County, Washington: U.S. Geological Survey Water-Resources Investigations Report 20034244, 137 p. Kolodny, Yehoshua, Katz, Amitai, Starinsky, Abraham, Moise, Tamar, and Simon, Ehud, 1999, Chemical tracing of salinity sources in Lake Kinneret (Sea of Galilee), Israel: Limnology and Oceanography, v. 44, no. 4, p. 10351044. Krabbenhoft, D.P., and Anderson, M.P., 1986, Use of a numerical ground-water flow model for hypothesis testing: Ground Water, v. 24, p. 4955. Krabbenhoft, D.P., Anderson, M.P., and Bowser, C.J., 1990, Estimating groundwater exchange with lakes, 2Calibration of a three-dimensional, solute transport model to a stable isotope plume: Water Resources Research, v. 26, no.10, p.24552462. Krabbenhoft, D.P., Bowser, C.J., Kendall, C., and Gat, J.R., 1994, Use of oxygen-18 and deuterium to assess the hydrology of groundwater-lake systems, in Baker, L.A., ed., Environmental Chemistry of Lakes and Reservoirs: American Chemical Society, Advances in Chemistry Series no. 237, p.6790. Kraemer, T.F., 2005, Radium isotopes in Cayuga Lake, New YorkIndicators of inflow and mixing processes: Limnology and Oceanography, v. 50, no. 1, p. 158168. Kraemer, T.F., and Genereux, D.P., 1998, Applications of uranium and thorium series radionuclides in catchment hydrology studies, in Kendall, Carol, and McDonnell, J.J., eds., Isotope tracers in catchment hydrology: Amsterdam, Elsevier, p. 679722. Krest, J.M., and Harvey, J.W., 2003, Using natural distributions of short-lived radium isotopes to quantify ground water discharge and recharge: Limnology and Oceanography, v. 48, no. 1, p. 290298. Krest, J.M., Moore, W.S., Gardner, L.R., and Morris, J.T., 2000, Marsh nutrient export supplied by groundwater dischargeEvidence from radium measurements: Global Biogeochemical Cycles, v. 14, no. 1, p. 167176. LaBaugh, J.W., 1985, Uncertainty in phosphorus retention, Williams Fork Reservoir, Colorado: Water Resources Research, v. 21, no. 11, p. 16841692. LaBaugh, J.W., Rosenberry, D.O., and Winter, T.C., 1995, Groundwater contribution to the water and chemical budgets of Williams Lake, Minnesota, 19801991: Canadian Journal of Fisheries and Aquatic Science, v. 52, p. 754767. LaBaugh, J.W., Winter, T.C., Rosenberry, D.O., Schuster, P.F., Reddy, M.M., and Aiken, G.A., 1997, Hydrological and chemical estimates of the water balance of a closed-basin lake in north-central Minnesota: Water Resources Research, v. 33, no. 12, p. 27992812. Lapham, W.W., 1989, Use of temperature profiles beneath streams to determine rates of vertical ground-water flow and vertical hydraulic conductivity: U.S. Geological Survey Water-Supply Paper 2337, 35 p.

34 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water Leavesley, G.H., and Hay, L., 1998, The use of coupled atmospheric and hydrological models for water-resources management in headwater basins, in Kovar, K., Tappeiner, U., Peters, N.E., and Craig, R.G., eds., Hydrology, water resources, and ecology in headwaters, Meran, Italy: Meran, Italy IAHS-AISH Publication, p. 259265. Leavesley, G.H., Lichty, R.W., Troutman, B.M., and Saindon, L.G., 1983, Precipitation-runoff modeling systemUser's manual: U.S. Geological Survey Water-Resources Investigations Report 834238, 207 p. Leavesley, G.H., Markstrom, S.L., Restrepo, P.J., and Viger, R.J., 2002, A modular approach to addressing model design, scale, and parameter estimation issues in distributed hydrological modeling: Hydrological Processes, v. 16, p.173187. Leavesley, G.H., Restrepo, P.J., Markstrom, S.L., Dixon, M., and Stannard, L.G., 1996, The modular modeling system MMS, User's manual: U.S. Geological Survey Open-File Report 96151, 142 p. Lee, D.R., 1977, A device for measuring seepage flux in lakes and estuaries: Limnology and Oceanography, v. 22, no. 1, p.140147. Lee, D.R., 1985, Method for locating sediment anomalies in lakebeds that can be caused by groundwater flow: Journal of Hydrology, v. 79, p. 187193. Lee, D.R., and Cherry, J.A., 1978, A field exercise on groundwater flow using seepage meters and minipiezometers: Journal of Geological Education, v. 27, p. 620. Lee, D.R., Cherry, J.A., and Pickens, J.F., 1980, Groundwater transport of a salt tracer through a sandy lakebed: Limnology and Oceanography, v. 25, no. 1, p. 4561. Lee, T.M., and Swancar, Amy, 1997, Influence of evaporation, ground water, and uncertainty in the hydrologic budget of Lake Lucerne, a seepage lake in Polk County, Florida: U.S. Geological Survey Water-Supply Paper 2439, 61 p. Linderfelt, W.R., and Turner, J.V., 2001, Interaction between shallow groundwater, saline surface water and nutrient discharge in a seasonal estuaryThe SwanCanning system: Hydrological Processes, v. 15, p. 26312653. Lindgren, R.J., and Landon, M.K., 1999, Effects of groundwater withdrawals on the Rock River and associated valley aquifer, eastern Rock County, Minnesota: U.S. Geological Survey Water-Resources Investigations Report 994157, 103p. Lodge, D.M., Krabbenhoft, D.P., and Striegl, R.G., 1989, A positive relationship between groundwater velocity and submersed macrophyte biomass in Sparkling Lake, Wisconsin: Limnology and Oceanography, v. 34, no. 1, p. 235239. Loke, M.H., and Lane, J.W., Jr., 2004, Inversion of data from electrical resistivity imaging surveys in water-covered areas: Exploration Geophysics, v. 35, p. 266271. Malard, F., Plenet, S., and Gibert, J., 1996, The use of invertebrates in ground water monitoringA rising research field: Ground Water Monitoring and Remediation, v. 16, no. 2, p.103113. Malcolm, I.A., Soulsby, C., Youngson, A.F., and Petry, J., 2003a, Heterogeneity in ground-water/surface-water interactions in the hyporheic zone of a salmonid spawning stream: Hydrological Processes, v. 17, p. 601617. Malcolm, I.A., Youngson, A.F., and Soulsby, Chris, 2003b, Survival of salmonid eggs in a degraded gravel-bed streamEffects of ground-water/surface-water interactions: Hydrobiologia, v. 444, p. 303316. Manheim, F.T., Krantz, D.E., and Bratton, J.F., 2004, Studying ground water under Delmarva coastal bays using electrical resistivity: Ground Water, v. 42, no. 7, p. 10521068. McBride, M.S., and Pfannkuch, H.O., 1975, The distribution of seepage within lakebeds: U.S. Geological Survey Journal of Research, v. 3, no. 5, p. 505512. McCarthy, K.A., McFarland, W.D., Wilkinson, W.D., and White, L.D., 1992, The dynamic relationship between ground water and the Columbia RiverUsing deuterium and oxygen-18 as tracers: Journal of Hydrology, v. 135, p.112. McLeod, R.S., 1980, The effects of using ground water to maintain water levels of Cedar Lake, Wisconsin: U.S.Geological Survey Water-Resources Investigations Report8023, 35 p. Merritt, M.L., and Konikow, L.F., 2000, Documentation of a computer program to simulate lake-aquifer interaction using the MODFLOW Ground-Water Flow Model and the MOC3D Solute-Transport Model: U.S. Geological Survey Water-Resources Investigations Report 004167, 146 p. Metz, P.A., and Sacks, L.A., 2002, Comparison of the hydrogeology and water quality of a ground-water augmented lake with two non-augmented lakes in northwest Hillsborough County, Florida: U.S. Geological Survey WaterResources Investigations Report 024032, 74 p. Meyboom, P., 1966, Unsteady groundwater flow near a willow ring in hummocky moraine: Journal of Hydrology, v. 4, p.3862. Meyboom, P., 1967, Mass-transfer studies to determine the groundwater regime of permanent lakes in hummocky moraine of western Canada: Journal of Hydrology, v. 5, p.117142.

Introduction and Characteristics of Flow 35 Miller, R.L., Kraemer, T.F., and McPherson, B.F., 1990, Radium and radon in Charlotte Harbor Estuary, Florida: Estuarine, Coastal and Shelf Science, v. 31, p. 439457. Mitchell-Bruker, S., and Haitjema, H.M., 1996, Modeling steady state conjunctive groundwater and surface water flow with analytic elements: Water Resources Research, v. 32, no. 9, p. 27252732. Moore, W.S., 1996, Large groundwater inputs to coastal waters revealed by 226Ra enrichments: Nature, v. 380, p.612614. Moore, W.S., 1999, The subterranean estuaryA reaction zone of ground water and sea water: Marine Geochemistry, v. 65, p. 111125. Moore, W.S., 2000, Determining coastal mixing rates using radium isotopes: Continental Shelf Research, v. 20, p. 19932007. Mull, D.S., Liebermann, T.D., Smoot, J.L., and Woosley, L.H., Jr., 1988, Application of dye-tracing techniques for determining solute transport characteristics of ground water in karst terranes: U.S. Environmental Protection Agency Report 904/688001. Niswonger, R.G., and Prudic, D.E., 2005, Documentation of the streamflow-routing (SFR2) package to include unsaturated flow beneath streamsA modification to SFR1: U.S.Geological Survey Techniques and Methods 6A13, 48 p. Niswonger, R.G., Markstrom, R.L., Regan, R.S., Prudic, D.E., Pohll, Greg, and Viger, R.J., 2006, Modeling ground-water/ surface-water interactions with GSFLOW, a new USGS model: Golden, Colo., Proceedings, MODFLOW and More 2006, Managing Ground Water Systems, May 2124, 2006, p. 99103. Oberg, K.A., Morlock, S.E., and Caldwell, W.S., 2005, Quality-assurance plan for discharge measurements using acoustic Doppler current profilers: U.S. Geological Survey Scientific Investigations Report 20055183, 44 p. Orghidan, T., 1959, Ein neuer Lebensraum des unteriridischen Wassers: Der hyporheische Biotop: Archiv fr Hydrobiologie, v. 55, p. 392414. Owen-Joyce, S.J., Wilson, R.P., Carpenter, M.C., and Fink, J.B., 2000, Method to identify wells that yield water that will be replaced by water from the Colorado River downstream from Laguna Dam in Arizona and California: U.S.Geological Survey Water-Resources Investigations Report 004085, 31 p. 19 maps. Pfannkuch, H.O., and Winter, T.C., 1984, Effect of anisotropy and groundwater system geometry on seepage through lakebeds, 1, Analog and dimensional analysis: Journal of Hydrology, v. 75, p. 213237. Pluhowski, E.J., 1972, Hydrologic interpretations based on infrared imagery of Long Island, New York, Contributions to the hydrology of the United States: U.S. Geological Survey Water-Supply Paper 2009B, 20 p. Power, G., Brown, R.S., and Imhof, J.G., 1999, Groundwater and fishInsights from northern North America: Hydrological Processes, v. 13, p. 401422. Prudic, D.E., 1989, Documentation of a computer program to simulate stream-aquifer relations using the modular finite-difference ground-water flow model: U.S. Geological Survey Open-File Report 88729, 113 p. Prudic, D.E., Konikow, L.F., and Banta, E.R., 2004, A new streamflow-routing (SFR1) package to simulate streamaquifer interaction with MODFLOW-2000: U.S. Geological Survey Open-File Report 20041042, 95 p. Puckett, L.J., Cowdery, T.K., McMahon, P.B., Tornes, L.H., and Stoner, J.D., 2002, Using chemical, hydrologic, and age dating analysis to delineate redox processes and flow paths in the riparian zone of a glacial outwash aquifer-stream system: Water Resources Research, v. 38, no. 8, p. 91 to 920. Rantz, S.E., and others, 1982a, Measurement and computation of streamflowVolume 1, Measurement of stage and discharge: U.S. Geological Survey Water-Supply Paper2175, 284 p. Rantz, S.E., and others, 1982b, Measurement and computation of streamflowVolume 2, Computation of discharge: U.S.Geological Survey Water-Supply Paper 2175, 631 p. Rimmer, Alon, and Gideon, Gal, 2003, Estimating the saline springs component in the solute and water balance of Lake Kinneret, Israel: Journal of Hydrology, v. 284, p. 228243. Robinove, C.J., 1965, Infrared photography and imagery in water resources research: Journal of the American Water Works Association, v. 57, pt. 2, p. 834840. Robinove, C.J., and Anderson, D.G., 1969, Some guidelines for remote sensing in hydrology: Water Resources Bulletin, v. 5, no. 2, p. 1019. Rosenberry, D.O., 1990, Inexpensive groundwater monitoring methods for determining hydrologic budgets of lakes and wetlands, in National Conference on Enhancing the States Lake and Wetland Management Programs: U.S. Environmental Protection Agency, North American Lake Management Society, p. 123131. Rosenberry, D.O., 2005, Integrating seepage heterogeneity with the use of ganged seepage meters: Limnology and Oceanography/Methods, v. 3, p. 131142.

36 Field Techniques for Estimating Water Fluxes Between Surface Water and Ground Water Rosenberry, D.O., Bukaveckas, P.A., Buso, D.C., Likens, G.E., Shapiro, A.M., and Winter, T.C., 1999, Migration of road salt to a small New Hampshire lake: Water Air and Soil Pollution, v. 109, p. 179206. Rosenberry, D.O., Striegl, R.G., and Hudson, D.C., 2000, Plants as indicators of focused ground water discharge to a northern Minnesota lake: Ground Water, v. 38, no. 2, p.296303. Rosenberry, D.O., and Winter, T.C., 1997, Dynamics of watertable fluctuations in an upland between two prairie-pothole wetlands in North Dakota: Journal of Hydrology, v. 191, p.266289. Rundquist, D., Murray, G., and Queen, L., 1985, Airborne thermal mapping of a flow-through lake in the Nebraska sandhills: Water Resources Research, v. 21, no. 6, p. 989994. Rutledge, A.T., 1992, Methods of using streamflow records for estimating total and effective recharge in the Appalachian Valley and Ridge, Piedmont, and Blue Ridge physiographic provinces: American Water Resource Association Monograph Series no. 17, p. 5973. Rutledge, A.T., 1998, Computer programs for describing the recession of ground-water discharge and for estimating mean ground-water recharge and discharge from streamflow dataUpdate: U.S. Geological Survey Water-Resources Investigations Report 984148, 44 p. Rutledge, A.T., 2000, Considerations for use of the RORA program to estimate ground-water recharge from streamflow records: U.S. Geological Survey Open-File Report 2000156, 52 p. Sacks, L.A., 2002, Estimating ground-water inflow to lakes in central Florida using the isotope mass-balance approach: U.S. Geological Survey Water-Resources Investigations Report 024192, 59 p. Sacks, L.A., Swancar, Amy, and Lee, T.M., 1998, Estimating ground-water exchange with lakes using water-budget and chemical mass-balance approaches for ten lakes in ridge areas of Polk and Highlands Counties, Florida: U.S.Geological Survey Water-Resources Investigations Report984133, 52 p. Schindler, D.W., Newbury, R.W., Beaty, K.G., and Campbell, P., 1976, Natural water and chemical budgets for a small Precambrian lake basin in central Canada: Journal of the Fisheries Research Board of Canada, v. 33, no. 11, p.25262543. Sheets, R.A., Darner, R.A., and Whitteberry, B.L., 2002, Lag times of bank filtration at a well field, Cincinnati, Ohio, USA: Journal of Hydrology, v. 266, p. 162174. Simmons, G.M., Jr., 1992, Importance of submarine groundwater discharge (SGWD) and seawater cycling to material flux across sediment/water interfaces in marine environments: Marine Ecology Progress Series, v. 81, p. 173184. Simonds, F.W., and Sinclair, K.A., 2002, Surface-water/ ground-water interactions along the lower Dungeness River and vertical hydraulic conductivity of streambed sediments, Clallam County, Washington, September 1999July 2001: U.S. Geological Survey Water-Resources Investigations Report 024161, 60 p. Simonds, F.W., Longpre, C.I., and Justin, G.B., 2004, Groundwater system in the Chimacum Creek Basin and surfacewater/ground-water interaction in Chimacum and Tarboo Creeks and the Big and Little Quilcene Rivers, Eastern Jefferson County, Washington: U.S. Geological Survey Scientific Investigations Report 20045058, 60 p., one plate. Smart, P.L., and Laidlaw, I.M.S., 1977, An evaluation of some fluorescent dyes for water tracing: Water Resources Research, v. 13, p. 1533. Sophocleous, Marios, 2002, Interactions between groundwater and surface waterThe state of the science: Hydrogeology Journal, v. 10, p. 5267. Sophocleous, M.A., Koelliker, J.K., Govindaraju, R.S., Birdie, T., Ramireddygari, S.R., and Perkins, S.P., 1999, Integrated numerical modeling for basin-wide water management The case of the Rattlesnake Creek basin in south-central Kansas: Journal of Hydrology, v. 214, p. 179196. Sophocleous, Marios, and Perkins, S.P., 2000, Methodology and application of combined watershed and ground-water models in Kansas: Journal of Hydrology, v. 236, p. 185201. Steele, G.V., and Verstraeten, I.M., 1999, Effects of pumping collector wells on riverAquifer interactions at Platte River Island near Ashland, Nebraska, 1998: U.S. Geological Survey Water-Resources Investigations Report 994161, 6 p. Stewart, J.W., and Hughes, G.H., 1974, Hydrologic consequences of using groundwater to maintain lake levels affected by water wells near Tampa, Florida: Tallahassee, Florida, Bureau of Geology Report of Investigation, no. 74, 41 p. Stonestrom, D.A., and Constanz, Jim, eds., 2003, Heat as a tool for studying the movement of ground water near streams: U.S. Geological Survey Circular 1260, 105p. (http://water.usgs.gov/pubs/circ/2003/circ1260/pdf/ Circ1260.pdf) Strobel, M.L., Sawyer, J.F., and Rahn, P.H., 2000, Field trip road logHydrogeology of the central Black Hills of SouthDakota, in Strobel, M.L., and others, eds., Hydrology of theBlack Hills, Proceedings of the 1999 Conference on the Hydrology of the Black Hills: Rapid City, South Dakota School of Mines and Technology Bulletin no. 20, p.239245. Stromberg, J.C., Tiller, R., and Richter, B., 1996, Effects of groundwater decline on riparian vegetation of semiarid regionsThe San Pedro, Arizona: Ecological Applications, v. 6, no. 1, p. 113131.

Introduction and Characteristics of Flow 37 Swain, E.D., and Wexler, E.J., 1996, A coupled surface-water and ground-water flow model (MODBRANCH) for simulation of stream-aquifer interaction: U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, chap. A6, 125 p. Squillace, P.J., Thurman, E.M., and Furlong, E.T., 1993, Groundwater as a nonpoint source of atrazine and deethylatrazine in a river during base flow conditions: Water Resources Research, v. 29, no. 6, p. 17191729. Taniguchi, M., 1993, Evaluation of vertical groundwater fluxes and thermal properties of aquifers based on transient temperature-depth profiles: Water Resources Research, v.29, no. 7, p. 20212026. Taniguchi, Makato, Burnett, W.C., Cable, J.E., and Turner, J.V., 2002, Investigation of submarine groundwater discharge: Hydrological Processes, v. 16, p. 21152129. Taniguchi, M., Ishitobi, T., and Shimada, J., 2006, Dynamics of submarine groundwater discharge and freshwaterseawater interface: Journal of Geophysical Research, v. 111, C01008, doi:10.1029/2005JC002924. Taylor, J.I., and Stingelin, R.W., 1969, Infrared imaging for water resources studies: Journal of the Hydraulics Division, Proceedings of the American Society of Civil Engineers, v.95, no. 1, p. 175189. Thoms, R.B., Johnson, R.L., and Healy, R.W., 2006, Users guide to the Variably Saturated Flow (VSF) Process for MODFLOW: U.S. Geological Survey Techniques and Methods 6A18, 58 p. Vaccaro, J.J., and Maloy, K.J., 2006, A thermal profile method to identify potential ground-water discharge areas and preferred salmonid habitats for long river reaches: U.S.Geological Survey Scientific Investigations Report 20065136, 16 p. (http://pubs.water.usgs.gov/sir2006-5136/) Valiela, Ivan, Costa, Joseph, Foreman, Kenneth, Teal, J.M., Howes, Brian, and Aubrey, David, 1990, Transport of groundwater-borne nutrients from watersheds and their effects on coastal waters: Biogeochemistry, v. 10, p.177197. Vanek, Vladimir, and Lee, D.R., 1991, Mapping submarine groundwater discharge areasAn example from Laholm Bay, southwest Sweden: Limnology and Oceanography, v. 36, no. 6, p. 12501262. Wentz, D.A., Rose, W.J., and Webster, K.E., 1995, Long-term hydrologic and biogeochemical responses of a soft water seepage lake in north central Wisconsin: Water Resources Research, v. 31, no. 1, p. 199212. Wetzel, R.G., 1999, Plants and water in and adjacent to lakes, in Baird, A.J., and Wilby, R.L., eds., Eco-Hydrology: New York, Routledge, p. 269299. Winter, T.C., 1986, Effect of ground-water recharge on configuration of the water table beneath sand dunes and on seepage in lakes in the Sandhills of Nebraska, U.S.A.: Journal of Hydrology, v. 86, p. 221237. Winter, T.C., 1999, Relation of streams, lakes, and wetlands to groundwater flow systems: Hydrogeology Journal, v. 7, p.2845. Winter, T.C., 2001, The concept of hydrologic landscapes: Journal of the American Water Resources Association, v.37, no. 2, p. 335349. Winter, T.C., Harvey, J.W., Franke, O.L., and Alley, W.M., 1998, Ground water and surface water a single resource: U.S. Geological Survey Circular 1139, 79 p. Winter, T.C., LaBaugh, J.W., and Rosenberry, D.O., 1988, The design and use of a hydraulic potentiomanometer for direct measurement of differences in hydraulic head between groundwater and surface water: Limnology and Oceanography, v. 33, no. 5, p. 12091214. Winter, T.C., and Rosenberry, D.O., 1995, The interaction of ground water with prairie pothole wetlands in the Cottonwood Lake area, east-central North Dakota, 19791990: Wetlands, v. 15, no. 3, p. 193211. Woessner, W.W., 1998, Changing views of stream-groundwater interaction, in Van Brahana, J., Eckstein, Y., Ongley, L.K., Schneider, R., and Moore, J.E., eds., Gambling with groundwaterPhysical, chemical, and biological aspects of aquifer-stream relations: Las Vegas, Nev., Proceedings of the Joint Meeting of the XXVIII Congress of the International Association of Hydrogeologists and the annual meeting of the American Institute of Hydrologists, Las Vegas, Nevada, September 28 to October 2, 1998, p. 16. Woessner, W.W., 2000, Stream and fluvial plain ground water interactionsRescaling hydrogeologic thought: Ground Water, v. 38, no. 3, p. 423429. Woessner, W.W., and Sullivan, K.E., 1984, Results of seepage meter and mini-piezometer study, Lake Mead, Nevada: Ground Water, v. 22, no. 5, p. 561568. Zarriello, P.J., and Reis, K.G., III, 2000, A precipitation-runoff model for the analysis of the effects of water withdrawals on streamflow, Ipswich River Basin, Massachusetts: U.S. Geological Survey Water-Resources Investigations Report004029, 99 p. Zekster, I.S., 1996, Groundwater discharge into lakesA review of recent studies with particular regard to large saline lakes in central Asia: International Journal of Salt Lake Research, v. 4, p. 233249.