Groundwater  Processes and Concepts

1. Groundwater  Processes and Concepts S. Hughes, 2003

2. Porosity Porosity: total volume of soil that can be filled with water V = Total volume of element Vi = Volume of Pores Vs = Volume of solids Soil volume V (Saturated) Pore with water solid rm = particles density (grain density) rd = bulk density Void Ratio:

3. Porosity Porosity: total volume of soil that can be filled with water V = Total volume of element Vi = Volume of Pores Vs = Volume of solids Soil volume V (Saturated) Pore with water solid rm = particles density (grain density) rd = bulk density Void Ratio:

4. Typical Values of Porosity

5. Moisture Content Soil volume V (Unsaturated) Moisture Content (volumetric water content): volume of water in total volume V = Total volume of element Vi = Volume of Pores Vs = Volume of solids Saturation (% water content): Pore with air

6. Groundwater The major source of all fresh water drinking supplies in some countries is groundwater. Groundwater is stored underground in aquifers, and is highly vulnerable to pollution. Understanding groundwater processes and aquifers is crucial to the management and protection of this precious resource.

7. Groundwater Groundwater comes from precipitation. Precipitated water must filter down through the vadose zoneto reach the zone of saturation, where groundwater flow occurs. The vadose zonehas an important environmental role in groundwater systems. Surface pollutants must filter through the vadose zone before entering the zone of saturation. Subsurface monitoringof the vadose zone is used to locate plumes of contaminated water, tracking the direction and rate of plume movement.

8. Components of groundwater The rate of infiltration is a function of soil type, rock type, antecedent water, and time. The vadose zoneincludes all the material between the Earth’s surface and the zone of saturation. The upper boundary of the zone of saturation is called the water table. The capillary fringeis a layer of variable thickness that directly overlies the water table. Water is drawn up into this layer by capillary action.

9. Distribution of Water in Subsurface Moisture Profile Soil Profile Description • Different zones • depend on % of pore space filled with water • Unsaturated Zone • Water held by capillary forces, water content near field capacity except during infiltration • Soil zone • Water moves down (up) during infiltration (evaporation) • Capillary fringe • Saturated ar base • Field capacity at top • Saturated Zone • Fully saturated pores Field capacity - Water remaining after gravity drainage Wilting point - Water remaining after gravity drainage & evapotranspiration

10. Field Capacity • After infiltration - saturated • Drainage • Coarse soils: a few hours • Fine soils: a 2-3 days • After drainage • Large pores: air and water • Smaller pores: water only • Soil is at field capacity • Water and air are ideal for plant growth After drying After drainage Saturated

11. Wilting Point • After drainage • Root suction & evaporation • Soil dries out • More difficult for roots • After drying • Root suction not sufficient for plants – plant wilts • Wilting Point = soil water content when plant dies After drying After drainage Saturated

12. air Interface water Net forceinward No net force Surface Tension • Below interface • Forces act equally in all directions • At interface • Some forces are missing • Pulls molecules down and together • Like membrane exerting tension on the surface • Curved interface • Higher pressure on concave side • Pressure increase is balanced by surface tension • s = 0.073 N/m (@ 20oC) • Capillary pressure • Relates pressure on both sides of interface

13. Capillary Pressure Air Negative pressure Solid Solid Water Positive pressure r

14. Capillary Pressure Air Solid Solid Water r Liquid rises due to attractive force of pore, until gravity force stops it.

15. Subsurface Pressure Distribution • Capillary pressure head in zone above water table • Hydrostatic pressure distribution exists below the water table (p = 0). Ground surface Pressure is negative above water table Unsaturated zone Water table Pressure is positive below water table Saturated zone

16. Soil Water Characteristic Curves • Capillary pressure head • Function of: • Pore size distribution • Moisture content Vadose Zone Porosity Capillary Zone

17. Capillary Rise in Soils

18. Groundwater Movement SPECIFIC YIELD (Sy) is the ratio of the volume of water drained from a rock (due to gravity) to the total rock volume. Grain size has a definite effect on specific yield. Smaller grains have larger surface area/volume ratio, which means more surface tension. Fine-grained sediment will have a lower Sy than coarse-grained sediment. SPECIFIC RETENTION (Sr) is the ratio of the volume of water a rock can retain (in spite of gravity) to the total volume of rock. Specific yield plus specific retention equals porosity

19. Porosity, Specific Yield, & Specific Retention

20. Aquifers An aquifer is a formation that allows water to be accessible at a usable rate. Aquifers are permeable layers such as sand, gravel, and fractured rock. Confined aquifershave non-permeable layers, above and below the aquifer zone, referred to as aquitards or aquicludes. These layers restrict water movement. Clay soils, shales, and non-fractured, weakly porous igneous and metamorphic rocks are examples of aquitards.

21. Aquifers Sometimes a lens of non-permeable materialwill be found within more permeable material. Water percolating through the unsaturated zone will be intercepted by this layer and will accumulate on top of the lens. This water is a perched aquifer. An unconfined aquiferhas no confining layers that retard vertical water movement. Artesian aquifersare confined under hydraulic pressure, resulting in free-flowing water, either from a spring or from a well.

22. Aquifer Types

23. Aquifer Types

24. Aquifer Storage Coefficient • Fluid Compressibility (b) • Porous Medium Compressibility (a) • Confined Aquifer • Water produced by 2 mechanisms • Aquifer compaction due to increasing effective stress, DV = arg • Water expansion due to decreasing pressure, DV = bfrg S = rg(a+fb) • Unconfined aquifer • Water produced by draining pores, S = Sy

25. Recharge and Discharge Water is continually recycledthrough aquifer systems.Groundwater recharge is any water added to the aquifer zone. Processes that contribute to groundwater recharge include precipitation, streamflow, leakage (reservoirs, lakes, aqueducts), and artificial means (injection wells). Groundwater discharge is any process that removes water from an aquifer system. Natural springs and artificial wells are examples of discharge processes. S. Hughes, 2003

26. Recharge and Discharge Groundwater supplies 30% of the water present in our streams. Effluent streamsact as discharge zones for groundwater during dry seasons. This phenomenon is known as base flow. Groundwater overdraftreduces the base flow, which results in the reduction of water supplied to our streams. S. Hughes, 2003

27. Perennial Stream (effluent) • Humid climate • Flows all year -- fed by groundwater base flow (1) • Discharges groundwater S. Hughes, 2003

28. Ephemeral Stream (influent) • Semiarid or arid climate • Flows only during wet periods (flashy runoff) • Recharges groundwater S. Hughes, 2003

29. Springs • Discharge of groundwater • from a spring. • Springs generally emerge at • the base of a hillslope. • Some springs produce water • that has traveled for many • kilometers; while others emit • water that has traveled only • a few meters. • Springs represent places • where the saturated zone • (below the water table) • comes in contactwith the • land surface. S. Hughes, 2003

30. (from Keller, 2000, Figure 10.9) Summary of Groundwater Systems S. Hughes, 2003

31. Groundwater Movement -- General Concepts The water tableis actually a sloping surface. Slope (gradient)is determined by the difference in water table elevation (h) over a specified distance (L). Direction of flow is downslope. Flow ratedepends on the gradient and the properties of the aquifer.

32. Darcy’s Law Q = KIA -- Henry Darcy, 1856, studied water flowing through porous material. His equation describesgroundwater flow. • Darcy’s experiment: • Water is applied under pressure through end A, flows through the pipe, and discharges at end B. • Water pressure is measured using piezometer tubes Hydraulic head = dh(change in height between A and B) Flow length = dL (distance between the two tubes) Hydraulic gradient (I) = dh / dL

33. Darcy’s Law The velocity of groundwater is based on hydraulic conductivity (K), as well as the hydraulic head (I). The equation to describe the relations between subsurface materials and the movement of water through them is Q = KIA Q = Discharge = volumetric flow rate, volume of water flowing through an aquifer per unit time (m3/day) A = Areathrough which the groundwater is flowing, cross-sectional area of flow (aquifer width x thickness, in m2) Rearrange the equation to Q/A = KI, known as the flux (v), which is an apparent velocity Actual groundwater velocityis higher than that determined by Darcy’s Law.

34. Darcy’s Law FLUX given by v = Q/A = KIis the IDEAL velocityof groundwater; it assumes that water molecules can flow in a straight line through the subsurface. NOTE:Flux doesn't account for the water molecules actually following a tortuous path in and out of the pore spaces. They travel quite a bit farther and faster in reality than the flux would indicate. DARCY FLUXgiven by vx = Q/An = KI/n (m/sec)is the ACTUAL velocityof groundwater, which DOES account for tortuosity of flow paths by including porosity (n) in the calculation. Darcy velocity is higher than ideal velocity. Darcy’s Lawis used extensively in groundwater studies. It can help answer important questions such as the direction a pollution plume is moving in an aquifer, and how fast it is traveling. S. Hughes, 2003

35. Groundwater Movement The tortuous path of groundwater moleculesthrough an aquifer affects the hydraulic conductivity. How do the following properties contribute to the rate of water movement? • Clay content and • adsorptive properties • Packing density • Friction • Surface tension • Preferred orientation • of grains • Shape (angularity or • roundness) of grains • Grain size • Hydraulic gradient

36. Hydraulic Conductivity • A combined property of the medium and the fluid • Ease with which fluid moves through the medium k = intrinsic permeability ρ = density µ = dynamic viscosity g = gravitational constant Porous medium property Fluid properties

37. Hydraulic Conductivity • Specific discharge (q) per unit hydraulic gradient • Ease with which fluid it transorted through porous medium • Depends on both matrix and fluid properties • Fluid properties: • Density r, and • Viscosity m • Matrix properties • Pore size distribution • Pore shape • Tortuosity • Specific surface area • Porosity k= intrinsic permeability [L2]

38. Estimating ConductivityKozeny – Carman Equation • A combined property of the medium and the fluid • Ease with which fluid moves through the medium k = intrinsic permeability ρ = density g = gravitational constant µ = dynamic viscosity d = mean particle size f = porosity Kozeny – Carman eq.

39. Lab Measurement of ConductivityPermeameters • Darcy’s Law is useless unless we can measure the parameters • Set up a flow pattern such that • We can derive a solution • We can produce the flow pattern experimentally • Hydraulic Conductivity is measured in the lab with a permeameter • Steady or unsteady 1-D flow • Small cylindrical sample of medium

40. Lab Measurement of ConductivityConstant Head Permeameter • Flow is steady • Sample: Right circular cylinder • Length, L • Area, A • Constant head difference (b) is applied across the sample producing a flow rate Q • Darcy’s Law Continuous Flow Overflow Outflow Q A

41. Lab Measurement of ConductivityFalling Head Permeameter • Flow rate in the tube must equal that in the column Outflow Q

42. WELL SORTED Coarse (sand-gravel) WELL SORTED Fine (silt-clay) POORLY SORTED Coarse - Fine Permeability and Hydraulic Conductivity High Low Sorting of material affects groundwater movement. Poorly sorted (well graded) material is less porous than well-sorted material.

43. Groundwater Movement Porosity and hydraulic conductivity of selected earth materials • Hydraulic • Porosity Conductivity • Material (%) (m/day) • Unconsolidated • Clay 450.041 • Sand 35 32.8 • Gravel 25 205.0 • Gravel and sand 20 82.0 • Rock • Sandstone 15 28.7 • Dense limestone or shale 5 0.041 • Granite 1 0.0041

44. Aquifer Transmissivity • Transmissivity (T) • Discharge through thickness of aquifer per unit width per unit head gradient • Product of conductivity and thickness Hydraulic gradient = 1 m/m Potentiometric Surface Confining Bed Transmissivity, T, volume of water flowing an area 1 m x b under hydraulic gradient of 1 m/m 1 m b 1 m Confined Aquifer 1 m Conductivity, K, volume of water flowing an area 1 m x 1 m under hydraulic gradient of 1 m/m

45. Heterogeneity and Anisotropy • Homogeneous aquifer • Properties are the same at every point • Heterogeneous aquifer • Properties are different at every point • Isotropic aquifer • Properties are same in every direction • Anisotropic aquifer • Properties are different in different directions • Often results from stratification during sedimentation www.usgs.gov

46. Layered Porous Media(Flow Parallel to Layers) W Q b

47. Layered Porous Media(Flow Perpendicular to Layers) Q W b Q

48. Groundwater Flow Nets Water table contour lines are similar to topographic lines on a map. They essentially represent "elevations" in the subsurface. These elevations are the hydraulic head mentioned above. Water table contour lines can be used to determine the direction groundwater will flowin a given region. Many wells are drilled and hydraulic head is measured in each one. Water table contours (called equipotential lines) are constructed to join areas of equal head. Groundwater flow lines, which represent the paths of groundwater downslope, are drawn perpendicular to the contour lines. A map of groundwater contour lines with groundwater flow lines is called a flow net.

49. Well #1 4252’ elevdepth to WT = 120’ (WT elev = 4132’) Well #2 4315’ elevdepth to WT = 78’ (WT elev = 4237’) 4220 4200 4240 4180 • Calculate WT elevations. • Interpolate contour intervals. 4260 4140 4160 4180 4280 4200 N 4220 • Connect contours of equal elevation. • Draw flow lines perpendicular to contours. 4240 4260 4300 4280 Well #3 4397’ elevdepth to WT = 95’ (WT elev = 4302’) Groundwater Movement Determine flow direction from well data:

50. A simple flow net Cross-profile view Qal 100 50 WT Qal Aquitard (granite) • Effect of a producing well • Notice the approximate diameter of the cone of depression Aquitard Qal well Groundwater Flow Nets