In the field of hydrogeology, storage properties are physical properties that characterize the capacity of an aquifer to release groundwater. These properties are storativity (S), specific storage (Ss) and specific yield (Sy). According to Groundwater, by Freeze and Cherry (1979), specific storage, S s {\displaystyle S_{s}} [m−1], of a saturated aquifer is defined as the volume of water that a unit volume of the aquifer releases from storage under a unit decline in hydraulic head.

They are often determined using some combination of field tests (e.g., aquifer tests) and laboratory tests on aquifer material samples. Recently, these properties have been also determined using remote sensing data derived from Interferometric synthetic-aperture radar.

Storativity

Definition

Storativity or the storage coefficient is the volume of water released from storage per unit decline in hydraulic head in the aquifer, per unit area of the aquifer. Storativity is a dimensionless quantity, and is always greater than 0.

S = d V w d h 1 A = S s b + S y {\displaystyle S={\frac {dV_{w}}{dh}}{\frac {1}{A}}=S_{s}b+S_{y}\,}

  • V w {\displaystyle V_{w}} is the volume of water released from storage ([L3]);
  • h {\displaystyle h} is the hydraulic head ([L])
  • S s {\displaystyle S_{s}} is the specific storage
  • S y {\displaystyle S_{y}} is the specific yield
  • b {\displaystyle b} is the thickness of aquifer
  • A {\displaystyle A} is the area ([L2])

Confined

For a confined aquifer or aquitard, storativity is the vertically integrated specific storage value. Specific storage is the volume of water released from one unit volume of the aquifer under one unit decline in head. This is related to both the compressibility of the aquifer and the compressibility of the water itself. Assuming the aquifer or aquitard is homogeneous:

S = S s b {\displaystyle S=S_{s}b\,}

Unconfined

For an unconfined aquifer, storativity is approximately equal to the specific yield (S y {\displaystyle S_{y}}) since the release from specific storage (S s {\displaystyle S_{s}}) is typically orders of magnitude less (S s b ≪ S y {\displaystyle S_{s}b\ll \!\ S_{y}}).

S = S y {\displaystyle S=S_{y}\,}

The specific storage is the amount of water that a portion of an aquifer releases from storage, per unit mass or volume of the aquifer, per unit change in hydraulic head, while remaining fully saturated.

Mass specific storage is the mass of water that an aquifer releases from storage, per mass of aquifer, per unit decline in hydraulic head:

( S s ) m = 1 m a d m w d h {\displaystyle (S_{s})_{m}={\frac {1}{m_{a}}}{\frac {dm_{w}}{dh}}}

where

( S s ) m {\displaystyle (S_{s})_{m}} is the mass specific storage ([L−1]);

m a {\displaystyle m_{a}} is the mass of that portion of the aquifer from which the water is released ([M]);

d m w {\displaystyle dm_{w}} is the mass of water released from storage ([M]); and

d h {\displaystyle dh} is the decline in hydraulic head ([L]).

Volumetric specific storage (or volume-specific storage) is the volume of water that an aquifer releases from storage, per volume of the aquifer, per unit decline in hydraulic head (Freeze and Cherry, 1979):

S s = 1 V a d V w d h = 1 V a d V w d p d p d h = 1 V a d V w d p γ w {\displaystyle S_{s}={\frac {1}{V_{a}}}{\frac {dV_{w}}{dh}}={\frac {1}{V_{a}}}{\frac {dV_{w}}{dp}}{\frac {dp}{dh}}={\frac {1}{V_{a}}}{\frac {dV_{w}}{dp}}\gamma _{w}}

where

S s {\displaystyle S_{s}} is the volumetric specific storage ([L−1]);

V a {\displaystyle V_{a}} is the bulk volume of that portion of the aquifer from which the water is released ([L3]);

d V w {\displaystyle dV_{w}} is the volume of water released from storage ([L3]);

d p {\displaystyle dp} is the decline in pressure(N•m−2 or [ML−1T−2]);

d h {\displaystyle dh} is the decline in hydraulic head ([L]) and

γ w {\displaystyle \gamma _{w}} is the specific weight of water (N•m−3 or [ML−2T−2]).

In hydrogeology, volumetric specific storage is much more commonly encountered than mass specific storage. Consequently, the term specific storage generally refers to volumetric specific storage.

In terms of measurable physical properties, specific storage can be expressed as

S s = γ w ( β p + n ⋅ β w ) {\displaystyle S_{s}=\gamma _{w}(\beta _{p}+n\cdot \beta _{w})}

where

γ w {\displaystyle \gamma _{w}} is the specific weight of water (N•m−3 or [ML−2T−2])

n {\displaystyle n} is the porosity of the material (dimensionless ratio between 0 and 1)

β p {\displaystyle \beta _{p}} is the compressibility of the bulk aquifer material (m2N−1 or [LM−1T2]), and

β w {\displaystyle \beta _{w}} is the compressibility of water (m2N−1 or [LM−1T2])

The compressibility terms relate a given change in stress to a change in volume (a strain). These two terms can be defined as:

β p = − d V t d σ e 1 V t {\displaystyle \beta _{p}=-{\frac {dV_{t}}{d\sigma _{e}}}{\frac {1}{V_{t}}}}

β w = − d V w d p 1 V w {\displaystyle \beta _{w}=-{\frac {dV_{w}}{dp}}{\frac {1}{V_{w}}}}

where

σ e {\displaystyle \sigma _{e}} is the effective stress (N/m2 or [MLT−2/L2])

These equations relate a change in total or water volume (V t {\displaystyle V_{t}} or V w {\displaystyle V_{w}}) per change in applied stress (effective stress — σ e {\displaystyle \sigma _{e}} or pore pressure — p {\displaystyle p}) per unit volume. The compressibilities (and therefore also Ss) can be estimated from laboratory consolidation tests (in an apparatus called a consolidometer), using the consolidation theory of soil mechanics (developed by Karl Terzaghi).

Determination of the storage coefficient of aquifer systems

Aquifer-test analysis

Aquifer-test analyses provide estimates of aquifer-system storage coefficients by examining the drawdown and recovery responses of water levels in wells to applied stresses, typically induced by pumping from nearby wells.

Stress-strain analysis

Elastic and inelastic skeletal storage coefficients can be estimated through a graphical method developed by Riley. This method involves plotting the applied stress (hydraulic head) on the y-axis against vertical strain or displacement (compaction) on the x-axis. The inverse slopes of the dominant linear trends in these compaction-head trajectories indicate the skeletal storage coefficients. The displacements used to build the stress-strain curve can be determined by extensometers, InSAR or levelling.

Laboratory consolidation tests

Laboratory consolidation tests yield measurements of the coefficient of consolidation within the inelastic range and provide estimates of vertical hydraulic conductivity. The inelastic skeletal specific storage of the sample can be determined by calculating the ratio of vertical hydraulic conductivity to the coefficient of consolidation.

Model simulations and calibration

Simulations of land subsidence incorporate data on aquifer-system storage and hydraulic conductivity. Calibrating these models can lead to optimized estimates of storage coefficients and vertical hydraulic conductivity.

Specific yield

Values of specific yield
MaterialSpecific Yield (%)
minavgmax
Unconsolidated deposits
Clay025
Sandy clay (mud)3712
Silt3819
Fine sand102128
Medium sand152632
Coarse sand202735
Gravelly sand202535
Fine gravel212535
Medium gravel132326
Coarse gravel122226
Consolidated deposits
Fine-grained sandstone21
Medium-grained sandstone27
Limestone14
Schist26
Siltstone12
Tuff21
Other deposits
Dune sand38
Loess18
Peat44
Till, predominantly silt6
Till, predominantly sand16
Till, predominantly gravel16

Specific yield, also known as the drainable porosity, is a ratio and is the volumetric fraction of the bulk aquifer volume that a given aquifer will yield when all the water is allowed to drain out of it under the forces of gravity:

S y = V w d V T {\displaystyle S_{y}={\frac {V_{wd}}{V_{T}}}}

where

V w d {\displaystyle V_{wd}} is the volume of water drained, and

V T {\displaystyle V_{T}} is the total rock or material volume

It is primarily used for unconfined aquifers since the elastic storage component, S s {\displaystyle S_{s}}, is relatively small and usually has an insignificant contribution. Specific yield can be close to effective porosity, but several subtleties make this value more complicated than it seems. Some water always remains in the formation, even after drainage; it clings to the grains of sand and clay. Also, the value of a specific yield may not be fully realized for a very long time due to complications caused by unsaturated flow. Problems related to unsaturated flow are simulated using the numerical solution of Richards Equation, which requires estimation of the specific yield, or the numerical solution of the Soil Moisture Velocity Equation, which does not require estimation of the specific yield.

See also

  • Freeze, R.A. and J.A. Cherry. 1979. Groundwater. Prentice-Hall, Inc. Englewood Cliffs, NJ. 604 p.
  • Morris, D.A. and A.I. Johnson. 1967. Summary of hydrologic and physical properties of rock and soil materials as analyzed by the Hydrologic Laboratory of the U.S. Geological Survey 1948-1960. U.S. Geological Survey Water Supply Paper 1839-D. 42 p.
  • De Wiest, R. J. (1966). On the storage coefficient and the equations of groundwater flow. Journal of Geophysical Research, 71(4), 1117–1122.

Specific