In the first part of this module, we will focus on the properties of aquifers: What characteristics of a rock or sediment make it a good aquifer? What are the different kinds of aquifers? Fundamentally, the ability to store and transmit water are the two key ingredients that make a subsurface geological formation useful as an aquifer. In Module 6.1, we will explore the detailed physical properties of rocks and sediments that ultimately affect the storage and movement of groundwater. We'll also illustrate with a series of well-known examples of large aquifers tapped for drinking, industrial, and agricultural uses.
In completing this lesson, you will:
Aquifers come in many shapes, sizes, and “flavors”. For example, some aquifer systems span hundreds – or even thousands – of kilometers across several states or nations and may include multiple individual layers of rock or sediment that total thousands of meters in thickness. Other aquifers may be restricted to a few kilometers within a stream valley, and be only a few meters thick.
Nonetheless, there are some important common features of different kinds of aquifers. In this section of the module, we provide a basic overview of aquifers: What is an aquifer? What are the different types of aquifers, and what is their anatomy?
Aquifers are geologic formations in the subsurface that can store & transmit water (Figures 1 and 2). As we will see later in the section titles Basic Aquifer Properties, there are specific rock and soil properties that govern these two functions. Adequate storage requires that there be sufficient void space between particles, in fractures, or generated by compressing the aquifer under pressure, to provide usable quantities of water. Adequate transmission requires that the void spaces where water occurs be well enough connected that it can percolate or flow under either natural or pumping-driven conditions, at a rate that will support sustained use.
These definitions are intentionally vague because they depend on the scale of intended use. For example, an aquifer that provides water for a large city will need to sustain higher pumping rates at wells (on order of tens of thousands of gallons per minute) than one that provides for a single-family (a few to perhaps ten gallons per minute). For reference, Penn State University relies almost exclusively on groundwater pumping for its water supply at the University Park campus, with a total extraction of ~2.5 million gallons per day (about 1750 gallons per minute) distributed among several pumping well fields.
In contrast to an aquifer, an aquitard, often also termed a confining layer or aquiclude, is a geologic formation in the subsurface that does not transmit water effectively – and therefore acts as a barrier to groundwater flow. In general, aquifers are usually composed of sediments or sedimentary rocks with grain sizes larger than fine- to medium sand (>~125 µm diameter), or of fractured rock. Aquitards are typically composed of fine-grained sediments or sedimentary rocks or those in which the pore spaces have been filled by mineral cements (silts, siltstones, shales, clays, cemented sandstones, or unfractured limestones).
In the simplest sense, you might imagine an aquifer formation that may be covered by a veneer of soil, and which extends downward from a few feet below the surface for several tens or even hundreds of feet (Figure 3). At some depth below the land surface, the interstices between soil or sediment particles, or the fractures in the rock, will be water-filled or saturated. Shallower than that depth, these interstices, or pore spaces, will be filled with air, water vapor, and some liquid water bound to the surfaces of the rock (Figures 3-5). This zone is the unsaturated zone, also known as the “zone of soil moisture” or the vadose zone. The water table marks the top of the ground water system and is formally defined as the depth at which the pressure in the subsurface is equal to the atmospheric pressure. Immediately above the water table, there is a narrow zone of saturation termed the capillary fringe. In this zone, water is wicked upward in pore spaces due to capillary forces. This is analogous to capillary tube experiments you may have seen or performed in physics or chemistry classes in high school; it occurs due to interaction between the polar water molecule and the surfaces of the solids, and is directly related to the fact that water has a surface tension (as you may remember from Module 1!). In the capillary fringe, pores are saturated, but pressures are sub-atmospheric, meaning that the water is under suction as it is pulled or wicked upward. The nomenclature “fringe” reflects the fact that slight variations in grain size lead to variations in the height that water is drawn (again, think to the capillary tube experiment, and effects of different tube sizes) (Figure 5).
Any precipitation or surface water that infiltrates to the water table must percolate through the vadose zone in order to recharge the aquifer. As we will see later in this module infiltration and recharge typically constitute only a small fraction – rarely more than 10% - of precipitation, because most water that falls in events is returned to the atmosphere by transpiration or evaporation, becomes runoff (i.e., if the capacity for infiltration is exceeded by the rate of precipitation), or is bound by soils in the vadose zone.
In more detail, there are three main classifications of aquifers, defined by their geometry and relationship to topography and the subsurface geology (Figures 6-9). The simple aquifer shown in Figure 6 is termed an unconfined aquifer because the aquifer formation extends essentially to the land surface. As a result, the aquifer is in pressure communication with the atmosphere. Unconfined aquifers are also known as water table aquifers because the water table marks the top of the groundwater system.
A second common type of aquifer is a confined aquifer, which is isolated from pressure communication with overlying or underlying geologic formations – and with the land surface and atmosphere – by one or more confining layers or confining units. Confined aquifers differ from unconfined aquifers in two fundamental and important ways. First, confined aquifers are typically under considerable pressure, which may be derived from recharge at a higher elevation or from the weight of the overlying rock and soil (known as the overburden). In some cases, the pressure is high enough that wells drilled into the aquifer are free-flowing. This condition requires that the water pressure in the aquifer is sufficient to drive water up the wellbore and above the land surface, and such wells are called artesian wells (Figure 7). Second, confined aquifers typically remain saturated over their entire thickness, even as water is removed by pumping wells. Water extracted from the aquifer comes only from the depressurization of the aquifer – a combination of depressurization and expansion of the water itself, and relaxation of the aquifer formation upon reduction in pressure (Figure 8).
The third main type of aquifer is a perched aquifer (Figure 6). Perched aquifers occur above discontinuous aquitards, which allow groundwater to “mound” above them. Thee aquifers are perched, in that they sit above the regional water table, and within the regional vadose zone (i.e. there is an unsaturated zone below the perched aquifer). The dimensions of perched aquifers are typically small (dictated by climate conditions and the size of aquitard layers), and the volume of water they contain is sensitive to climate conditions and therefore highly variable in time.
In the previous section, you’ve learned about the different types of aquifers, and the basic characteristics that define an aquifer – namely the ability to store and transmit water. But what, exactly, about a rock or sediment beneath the ground determines whether the rock can hold water, or whether water can percolate through it? In the following section, we will explore this question in more detail, to define the important individual properties of the rock.
Porosity (usually denoted by the symbol η, which is Greek letter 'eta') is the primary aquifer property that controls water storage, and is defined as the volume of void space (i.e., that can hold water in the zone of saturation) as a proportion of the total volume (Figure 10).
Porosity is expressed as either a fraction, or a percentage:
, or if reported as a percentage,
For aquifers composed of sedimentary rocks or sediments, porosity is usually in the range of ~10-35%. For unfractured crystalline rock, porosity is quite a bit lower - on order of a few percent - because there is little porosity between individual grains other than the vary narrow interfaces along their boundaries (Figure 10).
Several factors can affect porosity. In sedimentary rock and sediments, controls on porosity include sorting, cementation, overburden stress (related to burial depth), and grain shape. Poorly sorted sedimentary deposits, in which there is a wide distribution of grain sizes, typically have lower porosity than well-sorted ones (Figure 11). This is because the finer particles are able to fill in spaces between the larger grains. Cementation caused by precipitation of minerals (typically calcium carbonate or silica) at grain boundaries also reduces porosity (Figure 11). Angular grains generally allow more efficient packing of particles than rounded or spherical ones, also leading to slightly lower porosity. Finally, the more deeply sediments or sedimentary rock are buried, the larger the weight of the overburden; the higher stress leads to compaction, tighter packing of the grains, and lower overall porosity.
Secondary processes that act on the rock or sediment after its formation, primarily weathering by physical or chemical mechanisms, can also affect porosity. Physical weathering by wind or water movement can remove fine clay-sized particles from the sediment (a process termed winnowing), leading to increased porosity near the Earth’s surface. Chemical weathering of certain rock types can lead to clay and oxide formation; depending on the environment and initial composition of the aquifer grains, the clays and oxides may subsequently be removed (porosity increase), or they may grow at the boundaries of other particles and reduce the porosity.
In fractured rock (whether fractured crystalline or cemented sedimentary rock), porosity is typically ~2-5%. The pore space is almost entirely composed of the fractures or cracks themselves, which are typically a millimeter or less in aperture (Figure 12). Two primary factors control porosity – and the connectedness of porosity – in fractured rock. First, increased stress, related to the depth of burial and the weight of the overburden, exerts a clamping force that causes the closure of cracks or fractures (decreases fracture aperture). In some types of rock – most notably limestones – chemical weathering occurs via dissolution as water flows along and through fractures. This leads to increased fracture aperture. Significant enlargement of fractures can lead to the development of karst, typified by large open fractures, caves, and caverns, as well as sinkholes and hummocky topography that ensue as the underlying rock is gradually dissolved (Figure 13).
Interestingly, grain size does not affect porosity. For example, consider a box filled with spherical particles packed as tightly as possible. The proportion of empty space (porosity) would be the same whether the particles are marble-sized, pea-sized, or golf ball-sized; the porosity is controlled entirely by the geometry of the particles – not their dimensions. As we’ll see in the next section, however, grain size does strongly affect the ability of aquifers to transmit water because it directly controls the size of the pore spaces where the water percolates. For example, unconsolidated clays (grain sizes of a few to tens of microns) commonly have porosities of over 50-60%, but they transmit water only one-thousandth to one millionth as well as sands with porosities of 20-30% (grain sizes of a few hundred microns).
Specific Yield (denoted as Sy) is another important quantity for water storage in unconfined aquifers. Sy is defined as the proportion of water occupying void spaces that drains under gravity. Because some water is bound, or adsorbed, to the aquifer particles or fractures, the specific yield is always lower than the porosity (take a look at Figure 8 [7], inset at top left). The attraction between water molecules and the aquifer is due – you guessed it! – to the polar nature of water and surface tension. In sands and fractured rock, Sy is typically a large fraction (>90%) of the porosity, whereas in fine-grained sedimentary deposits Sy may be as low as a few percent because the surface area interacting with water molecules is higher, and pores are smaller, allowing the aquifer to retain more water. In unconfined aquifers, Sy controls the amount of water that can be extracted by pumping.
In confined aquifers, the compressibility of the aquifer is the dominant control on water storage and release. As described above (see Figure 8, right panel), when water is extracted from confined aquifers by pumping or flow to natural springs, the aquifer remains saturated, but the water pressure decreases. Upon depressurization, the aquifer itself can compress slightly. If water is recharged or injected, the opposite occurs: pressure increases and the aquifer expands very slightly. Essentially, by increasing water pressure, more water mass is being “crammed” into the pore space in the aquifer, and vice versa. Although exaggerated, one way to visualize this is to think of pores in the rock or aquifer as a juice box. By changing the pressure inside the box, it will expand or contract. In the same way that a soft juice box will deform more than a stiff one for a given change in pressure, a more compressible aquifer will yield more water than a stiffer aquifer, for the same depressurization. The storage of water in confined aquifers is termed the specific storage, and reported in units of Volume of water/Volume of aquifer per change in water level (so the units are 1/length; e.g., 1/m or 1/ft).
The ability of an aquifer to transmit water – or of an aquitard to slow the flow of water – is the second essential ingredient controlling groundwater movement. It is also the most variable in natural materials; distances in astronomy are the only other quantity in nature that varies over a similar range! For example, the difference in groundwater flow rate for shale vs. gravel is a factor of 1,000,000,000,000 (yup…one trillion). That’s the difference between the size of an iPhone and the distance from the Earth to the Sun.
Groundwater transport properties are described by two related quantities. Hydraulic Conductivity, denoted by K, is a measure of the ability of a particular fluid (usually water) to flow through the rock or sediment. Permeability, denoted by a lower-case k, is often also termed “intrinsic permeability” and describes the ability of the geologic formation alone to transmit fluid. Although related, the key difference is that hydraulic conductivity combines properties of the geologic formation and the fluid, whereas permeability describes only the rock properties. As described in Sidebar 1, the basic concept of hydraulic conductivity emerged from a series of ingenious experiments conducted in the mid-1880s by Henri Darcy, a French Engineer. These experiments led to Darcy’s Law, which forms the foundation for much of modern hydrogeology and petroleum engineering.
To illustrate the difference between K and k, consider the sandstone in Figure 14 below. The sandstone itself has a permeability, which is controlled by the size of the grains and pore spaces through which water can percolate, and the connectedness and geometry of the pores (more on that in a moment!). That permeability is a characteristic of the sandstone, regardless of whatever fluid might be moving through it, the temperature, or anything else. But the flow rate of water through this sandstone will be different than for oil, or for air, or any other fluid. So the same sandstone also has a hydraulic conductivity specific to a given fluid of interest.
More specifically, it is the viscosity and density of the fluid that matter. More viscous fluids will flow more slowly through the same rock than less viscous ones. This is important for comparing different fluids (say, oil vs. water – whether you are thinking about an oil reservoir or contamination of groundwater by a gasoline spill). It is also important in considering the effects of temperature, because water viscosity decreases with increasing temperature: it’s less than half as viscous at 90° than at 32° F. So even for the same aquifer, the hydraulic conductivity goes up if it is warmer! This makes some sense – if the water is less viscous (i.e. “thinner”), it will flow more easily through the aquifer.
So…that’s how we define permeability and hydraulic conductivity. But what controls their magnitude? The main factors are grain size and shape, sorting, porosity (degree of compaction or fracture aperture), particle orientation or alignment that affects the tortuosity of the flow path, and cementation. Tortuosity is a measure of how far fluid must go to “circumnavigate” its way around particles: higher tortuosity indicates that water must go farther to get to its destination (a more tortuous path). For all of these mechanisms, the key underlying control on groundwater movement is the viscous resistance resulting from the interaction of the fluid with solid surfaces in the aquifer (grain edges or fracture walls).
Fractured aquifers are one important and widely used class of aquifer because they are commonly both highly permeable and rapidly recharged. For example, groundwater recharge to the limestone aquifer beneath Nittany Valley in the Spring Creek watershed is around 30-45% of the annual precipitation (in comparison to typical recharge of <10% of precipitation). Fractured aquifers are permeable despite their overall low porosity (usually <5%) because natural fractures usually form in consistent orientations and are well connected in networks over hundreds of meters to tens of kilometers or more (Figures 15-16). The preferred orientation of major fractures leads to anisotropy in permeability, in which the aquifer may be more permeable parallel to the dominant fracture directions than in other orientations.
The rapid flow rates and direct pathways for recharge from the land surface also lead to concerns specific to fractured aquifers. In the absence of confining layers or thick soils, rapid recharge along fractures that extend to or nearly to the land surface increase vulnerability of contamination by surface activity, including fertilization of fields, pesticide application, or spills. Direct connections between surface water bodies and groundwater through major fracture systems also increase the potential for water-borne pathogens to enter the groundwater system, especially during periods of high flow or if confining layers along stream beds are breached. Compounding this risk, if contamination does occur, flow along fracture networks can be very rapid and the direction and rates of contaminant transport difficult to predict - unless the fracture network in the subsurface is extraordinarily well known, which is rarely the case. Because of their potential for contamination, fractured aquifers are a subject of highly active research, including dedicated large-scale field programs (e.g., check out the U.S. Geological Survey’s Mirror Lake project [9]).
Ground water flow systems extend over a wide range of scales, from small perched aquifers that may supply water for a single-family, to regional rock formations that span thousands of km and cross several states (Figure 18). These regional systems supply water for irrigation and domestic uses in many areas, especially in semi-arid and arid parts of the American West and coastal population centers along the East coast (remember Module 1, figures 10-12?). These regional systems commonly consist of several layered sedimentary formations and may extend to several kilometers in depth. The U.S. Geological Survey has compiled detailed studies of regional aquifer systems across the U.S., with useful information about climate, recharge, subsurface geology, use, and problems related to water quality or quantity (a list and links for each of the principal regional aquifers in the U.S can be found at USGS Groundwater Information [11]. A detailed atlas with information about the major aquifer systems in particular regions of the U.S. can be viewed at USGS Ground Water Atlas of the United States [12]. In this module, we will focus on a few example regional aquifer systems of particular relevance to the Northeastern and mid-Atlantic U.S. and the Central Valley of CA.
The Valley and Ridge aquifer system extends SW-NE across Central PA, West VA, and VA (Figure 19, purple area), and is the main water supply for much of this region. It is composed of layered Paleozoic sedimentary strata (shales, sandstones, and limestones) that were folded and deformed by a series tectonic collisions over 200 million years ago. The modern valleys in the Valley and Ridge province have formed where limestone, which is most susceptible to erosion, was exposed in the core of anticlines, or upfolds (Figure 20). The more resistant sandstones and shales form the regional ridges, like Mt. Nittany and Bald Eagle Ridge.
The principal aquifer unit in this system is the fractured limestone that underlies the valleys. As noted above, because it is fractured, it recharges rapidly, has a high fracture permeability, and wells drilled along the fractures are highly productive (c.f. Figure 17). Recharge is focused on the flanks of the ridges, where runoff flows over the less permeable shale and sandstone units and enters the groundwater through fractures or sinkholes above the limestone at the valley edges. Groundwater flow is generally toward the center of the valleys, and springs commonly feed the surface water systems. The water is characterized by a high hardness (Mg and Ca content; we’ll cover this in more detail in Module 7), derived from limestone dissolution. Dissolution of the limestone has formed extensive karst features (caves, caverns, sinkholes) throughout the region.
The Atlantic Coastal Plain aquifer system extends North-South along much of the Eastern portions of New Jersey, Delaware, Pennsylvania, Virginia, and North Carolina (Figure 19). It consists of a sequence of layered sedimentary aquifers (sands and gravels) separated by series of aquitards, all deposited starting around 100 million years ago and continuing today. The layers slope, or dip, to the East and extend offshore for tens of km beneath the continental shelf (Figure 21).
Recharge occurs by both natural and managed infiltration on land across much of the coastal plain; groundwater flow in the subsurface is mainly to the East along the sediment layers. One interesting consequence of this flow pattern is that there may be a sizable freshwater resource offshore that could be accessed by drilling in relatively shallow water on the continental shelf. During the last ice age, when conditions were substantially wetter than today and a nearly mile-thick ice sheet covered the northern extent of the aquifer system, recharge was probably even larger - and thus may have forced fresh water several tens of km offshore, where that “fossil” water may remain today!
The Atlantic Coastal Plain system is an important water source for domestic/municipal supply and industry in population centers throughout coastal North Carolina, Maryland, Virginia, Delaware, and New Jersey. However, concentrated, localized pumping has led to a reversal of flow direction (toward the wells instead of Eastward) in some of the aquifer units throughout the region. In addition to overarching concerns about the sustainability of withdrawals that exceed recharge rates, the flow reversal has led to local salt-water intrusion, whereby saline ocean water infiltrates the aquifer and in some cases renders it non-potable.
The Central Valley Aquifer system of Central CA lies in a large structural basin running approximately North-South, between the Coast Ranges to the West and the Sierra Nevada mountains to the East (Figure 22). The deep elongate basin is infilled with marine and continental sediments, primarily composed of interlayered sands and clays. The basin itself is formed by tectonic processes caused by East-West extension (these are the same forces that are causing continued uplift of the major mountain chains throughout the Basin and Range province of the Southwestern US, and which are one major control on orographic precipitation patterns in that region).
The continental deposits (Figure 22, orange) comprise the main aquifer units and range from one-half to over two miles in thickness. As is the case in the Valley and Ridge, the recharge is primarily focused around the valley perimeter as runoff over the flanks of surrounding high topography infiltrates and enters the groundwater system. Groundwater flow is primarily inward, toward valley center, with a component of flow down-valley to the North, parallel to surface water flow in the San Joaquin River.
The thick sedimentary sequence has formed a vast expanse of flat topography on the natural floodplain of the San Joaquin River. This, in combination with a mild climate that allows a year-round growing season, has made the Central Valley one of the most productive and largest agricultural centers in the world. The Central Valley aquifer system is highly utilized, primarily to augment limited allocations of surface water for irrigation. Since the mid-1920s, groundwater withdrawals have generally outpaced natural recharge to the aquifer, leading to dropping water levels, irreversible aquifer compaction, and land subsidence (as will be discussed in more detail next week, in Module 6.2). Until recently, groundwater withdrawals were neither heavily monitored nor regulated. However, in the face of an ongoing multi-year drought, a 2014 bill was signed into law that restricts pumping and implements groundwater sustainability plans (See What to Know about California's New Groundwater Law [14]; see also New California Groundwater Pumping Rules Signed Into Law [15]). Shallow aquifer units in the valley are also plagued by a wide range of water quality concerns associated with irrigation and return flow of irrigation water to the aquifer via infiltration; these include leaching of selenium, boron, and other constituents from soils; salinization; and high concentrations of pesticides and fertilizers. We’ll discuss all of these issues in more detail in upcoming modules about water quality and the effects of climate change.
In 1855, Henri Darcy, a French hydraulic engineer (Figure 24), oversaw a series of experiments aimed to understand the rates of water flow through sand layers, and their relationship to pressure loss along the flow paths. Darcy’s experiments consisted of a vertical steel column, with a water inlet at one end and an outlet at the other. The water pressure was controlled at the inlet and outlet ends of the column using reservoirs with constant water levels (Figure 25) (denoted h1 and h2). The experiments included a series of tests with different packings of river sand, and a suite of tests using the same sand pack and column, but for which the inlet and outlet pressures were varied. For part of our in-class activity this week, we will perform our own set of “Darcy Tube” experiments, and also work with the original dataset generated by Darcy in his experiments.
Darcy’s findings laid the foundation for the modern science of hydrogeology by quantifying the relationships between volumetric groundwater flow rate, driving forces, and aquifer properties. Specifically, Darcy’s experiments revealed proportionalities between the flux of water, Q, through the laboratory “aquifer” and different characteristics of the experimental system (refer to Figure 25 above):
Q was directly proportional to the difference in water levels from inlet to outlet, :
- Q was directly proportional to the cross sectional area of the tube:
- Q was inversely proportional to the length of the column:
Combining these proportionalities leads to Darcy’s Law, the empirical law that describes groundwater flow:where K is a constant of proportionality that defines the water flux for a given hydraulic gradient . The above equation can also be recast in terms of the water volume flux per unit area, Q/A (also called "Darcy flux" or "Darcy velocity" with units of length per time):
Links
[1] http://pubs.usgs.gov/circ/circ1186/index.html
[2] http://water.usgs.gov/edu/capillaryaction.html
[3] http://water.usgs.gov/edu/earthgwaquifer.html
[4] http://water.usgs.gov/edu/gallery/artesian-well-georgia.html
[5] http://pubs.usgs.gov/sir/2008/5059/images/or19-0129_fig02.png
[6] http://earthsci.org/education/teacher/basicgeol/groundwa/groundwa.html#Porosity
[7] https://www.e-education.psu.edu/earth111/node/911
[8] http://dpipwe.tas.gov.au/water/groundwater/aquifers/dual-porosity-aquifers
[9] https://www.usgs.gov/publications/mirror-lake-past-present-and-future
[10] http://toxics.usgs.gov
[11] http://water.usgs.gov/ogw/aquiferbasics/alphabetical.html
[12] http://pubs.usgs.gov/ha/ha730/gwa.html
[13] http://water.usgs.gov/GIS/metadata/usgswrd/XML/aquifers_us.xml
[14] https://www.kqed.org/science/21706/what-to-know-about-californias-new-groundwater-law
[15] https://www.waterboards.ca.gov/water_issues/programs/sgma/about_sgma.html#:~:text=The%20Sustainable%20Groundwater%20Management%20Act%20(SGMA)%20was%20enacted%20in%20order,and%20medium-priority%20groundwater%20basins