Exercise 3B - Visualising groundwater
In this exercise you will investigate what controls the direction and rate of groundwater flow, and investigate groundwater storage and flow within a layered groundwater system.
Learning Objectives:
At the completion of this exercise students will be able to
calculate groundwater flow rates
describe the range of flow rates within natural groundwater systems and explain the key factors controlling these rates
identify patterns of groundwater flow and hydraulic head contours within layered groundwater systems
Groundwater flow rates
As a first step lets consider the flux law for groundwater flow. Groundwater flow is driven by the distribution of water pressure (hydraulic head, which has units of length) within the groundwater system. The flow of groundwater follows a flux law;
Where \(q\) is the specific discharge (or Darcy Flux) with units of L T-1, \(Q\) is the volumetric flow rate (L3 T-1), \(A\) is the cross-sectional area through which the water flow (L2), and \(\frac{\Delta h}{ \Delta l}\) is the hydraulic gradient. This gradient is simply the change in water pressure (\(h_1-h_2\)) over a given distance (L/L, unitless). The key system property in this flux law is the hydraulic conductivity (\(K\)), which has units of L T-1. In natural groundwater systems the value of hydraulic conductivity can vary over 12 orders of magnitude. Even though this water flux \(q\) has units of L T-1, it’s not actually the velocity of the groundwater, because the water wasn’t flowing through the whole area, just the pore space. To correct for this and calculate the average velocity we need to use the area of the pore space (\(An\), where \(n\) is porosity), instead of the total area. This gives us;
Where \(\bar{v}\) is the average linear velocity of the water flowing through the porous media.
Now that we know these simple equations let’s look at some experimental data from simple “Darcy tubes” that allow us to measure the relationships between rates of groundwater flow, hydraulic gradients and the hydraulic conductivity of the porous media that it flows through. This is how the idea of a “hydraulic conductivity” was first defined, empirically, and it was only decades later that the physical basis for hydraulic conductivity and Darcy’s Law was fully elucidated from first principles (see Hubbert 1940 for some light bedtime reading).
Figure 7: a) Photo of simple Darcy Tube and b) diagram showing terms in Darcy’s Law Equation
Let’s have a go at using this experimental setup to estimate the hydraulic conductivity of different sedimentary materials. To do this all we need to do is measure the flow rate through the tube for our chosen porous media. For each porous material, the slope determines the hydraulic gradient and we just need to measure the volume that flows out of the tube over a set time interval.
During this experiment;
the diameter of the tubes was 40 mm
the inlet at the top of the tube was 25 cm above the outlet at the bottom of the tube, and the lateral distance between the inlet and the outlet was 0.5 m
Over a period 24 hours 168 mL was collected from the tube filled with sand, 3.6 mL was collected from the tube filled with silty sand, and 814 mL of water was collected from the tube filled with gravel.
- Based on this information you can rearrange Eq (10) to calculate the hydraulic conductivity of each porous media. Are you results consistent with the range of values for similar porous media in the literature? Could you measure the hydraulic conductivity of a clay sample using this experimental setup? What could you do to improve this experiment?
- Given these hydraulic conductivities, what would the range of linear velocity we could expect in a groundwater system like the Gnangara Mound? Use the following watertable map (taken from Davidson and Yu 2008, HG20) to estimate a hydraulic conductivity, and then calculate linear velocity using Eq (11). You will need to assume a porosity - look these up, or just use 0.25 (which is a good indicative value for back-of-the-envelope calculations in aquifer materials - why not in clays?). Note that it’s probably ok to assume a porosity value because the range of porosity is within one order of magnitude, whereas hydraulic conductivity will vary by multiple orders of magnitude, even within the same type of material, so this will be the larger source of uncertainty in your calculation.
- How does your estimate of groundwater velocity compare to the streamflow rates you have been working with in previous Exercises? What is the implication for transport and storage of pollutants in the environment?
Figure 8: Section of a map of the watertable across Perth
Physical model of groundwater system
In this part of the exercise we will look at a model of a layered groundwater system. Each layer has a different hydraulic conductivity, \(K\), forming two aquifers and an aquitard. Within each layer there will also be variation in hydraulic conductivity. This spatial variation is called “heterogeneity,” and it is an important feature of real groundwater systems.
Start by annotating the photo of this model (Figure 8) to show the following;
Aquifers and aquitards
Recharge and discharge zones
Water levels in each of the wells
The water table
Hydraulic head boundaries
Figure 9: Physical model of a layered groundwater system.
Now watch this video to see how water moves through the model when dye is injected to see the direction of groundwater flow. Does the flow direction vary spatially? Does the flow direction change over time?
Draw arrows on the following screen-shot to show the direction of groundwater flow within the model domain. We have learnt that the measure of potential for groundwater flow is “hydraulic head” and we also know that groundwater flows from high hydraulic head to low hydraulic head. In fact, groundwater flow will be perpendicular to the equipotentials (lines of equal hydraulic head). Based on your inferred flow paths, try now to add in the hydraulic head contours (or equipotentials).
It might help for you to think about topographic contours on a landscape map - in this case you have elevation contours, and if you placed a ball at the top of a hill it would roll (or flow) from the higher elevation to the lower elevation. As a general rule, we also know that flow through high conductivity layers will be primarily horizontal, and flow through low conductivity layers will be primarily vertical.
Figure 10: Physical model of a layered groundwater system.
What can you learn from this physical model about groundwater systems? By using the physical model we have been able to see directly in the subsurface, but in the real world, the only information we have would be the wells.
If you only had the hydraulic head measurements in the wells what would you know about the groundwater system?
What else would you need to know about this system to calculate groundwater flow velocities (and therefore contaminant transport)?
How important is it that we understand the geology as well as having measurements of hydraulic heads?
What else might you measure in a real groundwater system to get a fuller picture of what is happening?
Computer model of groundwater system
We can also run a “digital twin” of the physical model that we have built in the lab. For this activity we will now use the ParFlow Sandtank model. This is a computer model simulating groundwater flow, and allows users to change aquifer configurations and explore the essential features of groundwater and contaminant flow.
Once you have watched the introductory video lets take a look at running the model. Navigate to this model website. You can learn how to customise the sandtank and run a simulation here. Have a go at manipulating the soil conditions and the hydraulic head sliders, and run simulations pumping or injecting groundwater. Can you reproduce our lab-based physical model?