Unit 03a : Advanced Hydrogeology Basin Hydrogeologic Cycle.

Unit 03a : Advanced Hydrogeology

Basin Hydrogeologic Cycle

Groundwater Basin

• A groundwater basin is a hydrologic unit of groundwater storage defined as an area more or less separate from neighboring groundwater storage areas.

• Basins are normally delimited by natural physical boundaries such as rivers, flow divides, and flow barriers.

Driving Forces

• The principal driving forces for groundwater flow were recognized very early by Chamberlin (1895) and King (1899):– Topographic (gravitational)– Thermal– Capillary

• Later tectonic strain (compaction) was recognized as an additional driving force.

Topographic Driving Forces

• Thomas S. Chamberlin (1885) identified elevation differences as the topographic drive for groundwater movement

• “the pressure-producing mechanism is the hydrostatic weight of the body of water….”

T.C. Chamberlin

F.H.King

• Franklin H. King (1899) stated that the water table everywhere is a subdued replica of the topography and that water moves from topographically high areas to topographically low areas

• “the dynamic mechanism required to maintain flow is continuous replenishment by precipitation….”

Chamberlin and later King defined and documented the confined and unconfined behaviour of groundwater systems over 100 years ago.

Conceptual Model of Regional Flow• M. King Hubbert (1940) was the first to publish a

conceptual model for a flow field based on head potential

• Hubbert’s model predicts the characteristics of the flow net in both recharge and discharge areas and is consistent with the water table as a subdued replica of the topographic surfaceM. King Hubbert

Hubbert’s model is essentially a 2D solution of Laplace’s equation where all boundaries except the upper surface are specified as no flow. The upper surface can be specified as any function representing the form of the water table. A general solution is:

(x,z) = ao + an cosh(nz/L).cos(nx/L) n =1

where is the hydraulic potential and L is the length of the flow cell

Hubbert’s ModelOne specific solution when the water table (upper surface boundary condition) is a simple cosine function:

(x,zo) = A – B cos(nx/L)has the form:

(x,zo) = A – B cosh(nz/L).cos(nx/L) cosh(nzo/L)

Multi-Cell Model of Regional Flow• J. Tóth (1962, 1963) was the first to significantly

extend the conceptual work of Hubbert• Tóth investigated a more complex flow system

with a sinusoidal water table superimposed on a regional slope

• Tóth identified local, intermediate and regional flow systems based on this simple topographic model.

Tóth’s solution uses a slightly more complex boundary condition than Hubbert:

(x,zo) = B’x/L – b sin(2x/)

The first (B’) term represents the regional slope and the second (b) term the more local sinusoidal relief. The parameter b is the amplitude of the topography and L/ is the number of flow cells

Tóth’s Model

Tóth’s Model• If b=0 (no local-scale topography) only a regional flow

system develops• If B’=0 (no regional-scale topography) only local flow

systems develop• If b=B’=0 (no topography) waterlogged conditions will

develop with the water table near the surface discharging by evapotranspiration

• If B’ and b > 0 and L/ >> 0 then regional, intermediate and local flow systems will develop.

Later Developments• Joe Tóth only considered an isotropic

homogeneous aquifer in his early work and used simple analytical models.

• Freeze and Witherspoon (1966,1967) used Tóth’s basic model but extended it’s application to layered systems using numerical models.

• Freeze and Witherspoon dealt primarily with the role of permeability contrasts in influencing flow lines in layered systems.

• These effects are important in understanding regional flow.

Flow Line Refraction• If k decreases with depth,

equipotentials crowd together and flow becomes more vertical

• If k increases with depth, equipotentials spread apart and flow becomes more horizontal

• If k increases significantly with depth, equipotentials become more widely spaced and flow becomes sub-horizontal

k1

k2

k1

k2

k1

k2

k1< k2

k1> k2

k1<< k2

Aquitard and Aquifer

• For regional flow systems where k2 /k1 or kaquifer/kaquitard tends to be 100 or greater

• Flows in aquitards (k1) are subvertical• Flows in aquifers (k2) are subhorizontal• The spacing of flow lines is a measure of flux so that

the aquifer is acting as a collector, concentrating flow.

k1

k2

k1<< k2

Regional Flow• To predict flow patterns in regional flow systems the

data requirements include:– Permeability distributions– Geometry of basin boundaries (including the surface

topography)

• The most influential factors on flow patterns are:– Ratio of basin depth to lateral extent– Configuration of the water table (topography)– Permeability distribution

• Major valleys collect flow and concentrate discharge• Deep permeable aquifers act as conduits and control

recharge rates and the location of recharge areas• Stratigraphic pinchouts at depth exert an influence

the location of recharge and discharge areas

Effect of Topography on Flow Patterns

Low relief topography

Moderate relief topography

High relief topography

Effect of Layers on Flow Patterns

Moderate K ratio

High K ratio

Effect of Dip on Flow Patterns

Dip towards recharge area – high relief

Dip towards discharge area – low relief

Effect of Pinchouts on Flow Patterns

High K below discharge area

High K below recharge area

Mountainous Terrain

• 20% of global flow systems• K distribution dominated by fracturing• Deep circulation –relatively high Kv for

fractured rock• Water circulates to depths where elevated

temperatures exist• Water table free surface relatively unrelated

to topography• Relationship between K and infiltration has

strong influence on flow patterns

Carbonate Terrain

• 10% of global systems supplying >25% global population with water supply

• Self-organized networks of solution enlarged fractures

• Dual porosity systems – matrix dominates storage – fractures dominate flow

• Tracer velocities up to 21 km/d (0.25 m/s) have been measured

• A few springs integrate flow from large areas

Groundwater in Coastal Regions

• Freshwater flow limits saltwater encroachment• Development disturbs natural balance and can

lead to major seawater intrusion• Control of saline intrusion achieved by:

– Artificial recharge of freshwater– Reduction and rearrangement of wells– Development of coastal trough to limit intrusion– Development of coastal pressure ridge– Installation of subsurface flow barriers

Saline Intrusion

• Small seawater wedge balanced by net outflow to ocean

• Seawater wedge advances inland as a result of pumping

• Mixing as a result of fluctuations in aquifer recharge

Ghyben-Hertzberg

• Interface assumed rigid• Fluids assumed immiscible• Segregated flow

• Pressures for fluid columns onshore and offshore are assumed to balance

sz =f(hf + z)

z = hf f / (s - f)

• Assuming f = 1000 kN/m3

and s = 1025 kN/m3 the depth to the saline interface

z 40 hf

• This is called the Ghyben-Hertzberg formula.

zz

hf

Glover Analysis

• Seepage face below sea level

• (x,z) is point on interface

• Q’ = Q/L is discharge per unit length of coast

z2 = 2Q’xf + { Q’f } 2

K(s – f) K(s - f)

• The depth to the interface at the coast (x=0) is:

z = Q’f K(s – f)

• The height of the water table for any x is:

hf = 2Q’x (s – f)

K f

• The greater the flow, the deeper the interface and the greater the gap x.

InterfaceFresh water

Salt water

m

hfx

z

Upconing

• When pump wells lower the head above a saline interface, the interface rises. This is phenomenon is called “upconing”.

• There is an analytical solution that estimates the rise due to upconing for small perturbations of the saline interface.

z = Q’f _ 2dK(s – f)

• The displacement process is unstable and premature breakthrough of salt water at the pump well occurs if z/d is greater than about 0.5.

dz

Surface Features

• Surface features of groundwater flow (mainly discharge) include:– Springs– Seeps– Saline Soils– Permanent and Ephemeral Streams– Marshes, Swamps, Bogs and Fens (Wetlands)– Ponds, Sloughs and Lakes

Prairie Profile

• A central topographic high bounded by areas of natural discharge (Meyboom, 1966)

• Geology is generaly low-K tills over higher-K intertill sands and gravels

• No streams - most discharge occurs by evapotranspiration

Recharge Recharge

SalineSeepage

SalineSeepage

Groundwater and Vegetation• Willow rings occur around

recharge sloughs and local discharge sloughs – phreatophytes with low alkali tolerance

• Saline soils associated with intermediate to regional discharge systems – halophytes with high salt tolerance (foxtails, salicornia)

• Discharge sloughs associated with region flow systems have high TDS and precipitate salts

Bogs and Fens

• Bogs are characterized by spongy peat deposits, acidic waters, and a floor covered by a thick carpet of sphagnum moss.

• Bogs receive all or most of their water from precipitation rather than from runoff, groundwater or streams.

• Fens are peat-forming wetlands that receive nutrients from sources other than precipitation: usually from upslope drainage from soils and from groundwater.

• Fens differ from bogs because they are less acidic and have higher nutrient levels. They support a much more diverse ecology.

Slough and Lake InteractionsRecharge

Discharge

Flow-through

Sloughs can provide recharge or receive discharge or act as both source and sink.

Recharge and Discharge Areas

• Water entering the ground in recharge areas is transmitted to distant points and results in a soil moisture deficit in the soils overlying the recharge area

• Water entering the soil in a discharge areas cannot overcome the upward gradient and is returned to the surface by evapotranspiration locally.

dry soil salinization

Basin Hydrologic Cycle• How much groundwater participates in the

basin hydrologic cycle?• Tóth estimates 90% of recharge never

penetrates deeper than about 80 m.• Tritium studies confirm this theoretical estimate.• Many streams receive much of their baseflow

component from the area within the nearest topographic high

• Regional flow components are small compared with locally derived flows.