Groundwater model

Last updated

Groundwater models are computer models of groundwater flow systems, and are used by hydrologists and hydrogeologists. Groundwater models are used to simulate and predict aquifer conditions.

Contents

Characteristics

Typical aquifer cross-section Aquifer en.svg
Typical aquifer cross-section

An unambiguous definition of "groundwater model" is difficult to give, but there are many common characteristics.

A groundwater model may be a scale model or an electric model of a groundwater situation or aquifer. Groundwater models are used to represent the natural groundwater flow in the environment. Some groundwater models include (chemical) quality aspects of the groundwater. Such groundwater models try to predict the fate and movement of the chemical in natural, urban or hypothetical scenario.

Groundwater models may be used to predict the effects of hydrological changes (like groundwater pumping or irrigation developments) on the behavior of the aquifer and are often named groundwater simulation models. Groundwater models are used in various water management plans for urban areas.

As the computations in mathematical groundwater models are based on groundwater flow equations, which are differential equations that can often be solved only by approximate methods using a numerical analysis, these models are also called mathematical, numerical, or computational groundwater models. [1]

The mathematical or the numerical models are usually based on the real physics the groundwater flow follows. These mathematical equations are solved using numerical codes such as MODFLOW, ParFlow, HydroGeoSphere, OpenGeoSys etc. Various types of numerical solutions like the finite difference method and the finite element method are discussed in the article on "Hydrogeology".

Inputs

For the calculations one needs inputs like:

The model may have chemical components like water salinity, soil salinity and other quality indicators of water and soil, for which inputs may also be needed.

Hydrological inputs

Hydrological factors at the soil surface determining the recharge Surface water cycle.svg
Hydrological factors at the soil surface determining the recharge

The primary coupling between groundwater and hydrological inputs is the unsaturated zone or vadose zone. The soil acts to partition hydrological inputs such as rainfall or snowmelt into surface runoff, soil moisture, evapotranspiration and groundwater recharge. Flows through the unsaturated zone that couple surface water to soil moisture and groundwater can be upward or downward, depending upon the gradient of hydraulic head in the soil, can be modeled using the numerical solution of Richards' equation [2] partial differential equation, or the ordinary differential equation Finite Water-Content method [3] as validated for modeling groundwater and vadose zone interactions. [4]

Operational inputs

The operational inputs concern human interferences with the water management like irrigation, drainage, pumping from wells, watertable control, and the operation of retention or infiltration basins, which are often of an hydrological nature.
These inputs may also vary in time and space.

Many groundwater models are made for the purpose of assessing the effects hydraulic engineering measures.

Boundary and initial conditions

Boundary conditions Boundary value problem-en.svg
Boundary conditions

Boundary conditions can be related to levels of the water table, artesian pressures, and hydraulic head along the boundaries of the model on the one hand (the head conditions), or to groundwater inflows and outflows along the boundaries of the model on the other hand (the flow conditions). This may also include quality aspects of the water like salinity.

The initial conditions refer to initial values of elements that may increase or decrease in the course of the time inside the model domain and they cover largely the same phenomena as the boundary conditions do.

The initial and boundary conditions may vary from place to place. The boundary conditions may be kept either constant or be made variable in time.

Parameters

Example of parameters of an irrigation cum groundwater model File-Saltmod8 (2).JPG
Example of parameters of an irrigation cum groundwater model

The parameters usually concern the geometry of and distances in the domain to be modelled and those physical properties of the aquifer that are more or less constant with time but that may be variable in space.

Important parameters are the topography, thicknesses of soil / rock layers and their horizontal/vertical hydraulic conductivity (permeability for water), aquifer transmissivity and resistance, aquifer porosity and storage coefficient, as well as the capillarity of the unsaturated zone. For more details see the article on hydrogeology.

Some parameters may be influenced by changes in the groundwater situation, like the thickness of a soil layer that may reduce when the water table drops and/the hydraulic pressure is reduced. This phenomenon is called subsidence. The thickness, in this case, is variable in time and not a parameter proper.

Applicability

The applicability of a groundwater model to a real situation depends on the accuracy of the input data and the parameters. Determination of these requires considerable study, like collection of hydrological data (rainfall, evapotranspiration, irrigation, drainage) and determination of the parameters mentioned before including pumping tests. As many parameters are quite variable in space, expert judgment is needed to arrive at representative values.

The models can also be used for the if-then analysis: if the value of a parameter is A, then what is the result, and if the value of the parameter is B instead, what is the influence? This analysis may be sufficient to obtain a rough impression of the groundwater behavior, but it can also serve to do a sensitivity analysis to answer the question: which factors have a great influence and which have less influence. With such information one may direct the efforts of investigation more to the influential factors.

When sufficient data have been assembled, it is possible to determine some of missing information by calibration. This implies that one assumes a range of values for the unknown or doubtful value of a certain parameter and one runs the model repeatedly while comparing results with known corresponding data. For example, if salinity figures of the groundwater are available and the value of hydraulic conductivity is uncertain, one assumes a range of conductivities and the selects that value of conductivity as "true" that yields salinity results close to the observed values, meaning that the groundwater flow as governed by the hydraulic conductivity is in agreement with the salinity conditions. This procedure is similar to the measurement of the flow in a river or canal by letting very saline water of a known salt concentration drip into the channel and measuring the resulting salt concentration downstream.

Dimensions

Groundwater models can be one-dimensional, two-dimensional, three-dimensional and semi-three-dimensional. Two and three-dimensional models can take into account the anisotropy of the aquifer with respect to the hydraulic conductivity, i.e. this property may vary in different directions.

One-, two- and three-dimensional

Two-dimensional model of subsurface drainage in a vertical plane DrainSection2.png
Two-dimensional model of subsurface drainage in a vertical plane
Three-dimensional grid, Modflow MODFLOW 3D grid.png
Three-dimensional grid, Modflow
  1. One-dimensional models can be used for the vertical flow in a system of parallel horizontal layers.
  2. Two-dimensional models apply to a vertical plane while it is assumed that the groundwater conditions repeat themselves in other parallel vertical planes (Fig. 4). Spacing equations of subsurface drains and the groundwater energy balance applied to drainage equations [5] are examples of two-dimensional groundwater models.
  3. Three-dimensional models like Modflow [6] require discretization of the entire flow domain. To that end the flow region must be subdivided into smaller elements (or cells), in both horizontal and vertical sense. Within each cell the parameters are maintained constant, but they may vary between the cells (Fig. 5). Using numerical solutions of groundwater flow equations, the flow of groundwater may be found as horizontal, vertical and, more often, as intermediate.

Semi three-dimensional

Map of a radial semi 3-dimensional model consisting of vertical concentrical cylinders through which the flow passes radially to the well Well arrangement.svg
Map of a radial semi 3-dimensional model consisting of vertical concentrical cylinders through which the flow passes radially to the well

In semi 3-dimensional models the horizontal flow is described by 2-dimensional flow equations (i. e. in horizontal x and y direction). Vertical flows (in z-direction) are described (a) with a 1-dimensional flow equation, or (b) derived from a water balance of horizontal flows converting the excess of horizontally incoming over the horizontally outgoing groundwater into vertical flow under the assumption that water is incompressible.

There are two classes of semi 3-dimensional models:

Continuous radial model

An example of a non-discretized radial model is the description of groundwater flow moving radially towards a deep well in a network of wells from which water is abstracted. [7] The radial flow passes through a vertical, cylindrical, cross-section representing the hydraulic equipotential of which the surface diminishes in the direction of the axis of intersection of the radial planes where the well is located.

Map of a two-dimensional grid over an alluvial fan for a prismatic semi 3-dimensional model, SahysMod Garmsar.gif
Map of a two-dimensional grid over an alluvial fan for a prismatic semi 3-dimensional model, SahysMod

Prismatically discretized model

Prismatically discretized models like SahysMod [8] have a grid over the land surface only. The 2-dimensional grid network consists of triangles, squares, rectangles or polygons. Hence, the flow domain is subdivided into vertical blocks or prisms. The prisms can be discretized into horizontal layers with different characteristics that may also vary between the prisms. The groundwater flow between neighboring prisms is calculated using 2-dimensional horizontal groundwater flow equations. Vertical flows are found by applying one-dimensional flow equations in a vertical sense, or they can be derived from the water balance: excess of horizontal inflow over horizontal outflow (or vice versa) is translated into vertical flow, as demonstrated in the article Hydrology (agriculture).

In semi 3-dimensional models, intermediate flow between horizontal and vertical is not modelled like in truly 3-dimensional models. Yet, like the truly 3-dimensional models, such models do permit the introduction of horizontal and vertical subsurface drainage systems.

Semiconfined aquifers with a slowly permeable layer overlying the aquifer (the aquitard) can be included in the model by simulating vertical flow through it under influence of an overpressure in the aquifer proper relative to the level of the watertable inside or above the aquitard.

Groundwater modeling software and references

See also

Footnotes

  1. Rushton, K.R., 2003, Groundwater Hydrology: Conceptual and Computational Models. John Wiley and Sons Ltd. ISBN   0-470-85004-3
  2. Richards, L. A. (1931), Capillary conduction of liquids through porous mediums, J. Appl. Phys., 1(5), 318–333.
  3. Ogden, F. L., W. Lai, R. C. Steinke, J. Zhu, C. A. Talbot, and J. L. Wilson (2015), A new general 1-D vadose zone solution method, Water Resour. Res., 51, doi:10.1002/2015WR017126.
  4. Ogden, F. L., W. Lai, R. C. Steinke, and J. Zhu (2015), Validation of finite water-content vadose zone dynamics method using column experiments with a moving water table and applied surface flux, Water Resour. Res., 51, doi:10.1002/2014WR016454.
  5. The energy balance of groundwater flow applied to subsurface drainage in anisotropic soils by pipes or ditches with entrance resistance, International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands. Download from : , or directly as PDF : Archived 2009-02-19 at the Wayback Machine Paper based on: R.J. Oosterbaan, J. Boonstra and K.V.G.K. Rao, 1996, The energy balance of groundwater flow. Published in V.P.Singh and B.Kumar (eds.), Subsurface-Water Hydrology, p. 153-160, Vol.2 of Proceedings of the International Conference on Hydrology and Water Resources, New Delhi, India, 1993. Kluwer Academic Publishers, Dordrecht, The Netherlands. ISBN   978-0-7923-3651-8. Download as PDF : . The corresponding EnDrain model can be downloaded from:
  6. Online guide to MODFLOW-2000 and MODFLOW-2005
  7. Subsurface drainage by (tube)wells: well spacing equations for fully and partially penetrating wells in uniform or layered aquifers with or without anisotropy and entrance resistance. Paper explaining the basics of the WellDrain model, International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands. Updated version of: R.J.Oosterbaan, Tube well spacing formulas for subsurface drainage. Published in: Smith, K.V.H. and D.W. Rycroft (eds.) Hydraulic Design in Water Resources Engineering: Land Drainage. Proceedings of the 2nd International Conference, Southampton University, p. 75‑84. Springer‑Verlag, Berlin, 1986. Download as PDF : . Download the WellDrain model from:
  8. ILRI, 1995. SahysMod: Spatial Agro-Hydro-Salinity Model. Description of Principles, User Manual, and Case Studies. International Institute for Land Reclamation and Improvement (ILRI), Wageningen, Netherlands. On line:

Related Research Articles

<span class="mw-page-title-main">Aquifer</span> Underground layer of water-bearing permeable rock

An aquifer is an underground layer of water-bearing material, consisting of permeable or fractured rock, or of unconsolidated materials. Aquifers vary greatly in their characteristics. The study of water flow in aquifers and the characterization of aquifers is called hydrogeology. Related terms include aquitard, which is a bed of low permeability along an aquifer, and aquiclude, which is a solid, impermeable area underlying or overlying an aquifer, the pressure of which could lead to the formation of a confined aquifer. The classification of aquifers is as follows: Saturated versus unsaturated; aquifers versus aquitards; confined versus unconfined; isotropic versus anisotropic; porous, karst, or fractured; transboundary aquifer.

<span class="mw-page-title-main">Hydrogeology</span> Study of the distribution and movement of groundwater

Hydrogeology is the area of geology that deals with the distribution and movement of groundwater in the soil and rocks of the Earth's crust. The terms groundwater hydrology, geohydrology, and hydrogeology are often used interchangeably.

In science and engineering, hydraulic conductivity, is a property of porous materials, soils and rocks, that describes the ease with which a fluid can move through the pore space, or fracture network. It depends on the intrinsic permeability of the material, the degree of saturation, and on the density and viscosity of the fluid. Saturated hydraulic conductivity, Ksat, describes water movement through saturated media. By definition, hydraulic conductivity is the ratio of volume flux to hydraulic gradient yielding a quantitative measure of a saturated soil's ability to transmit water when subjected to a hydraulic gradient.

Used in hydrogeology, the groundwater flow equation is the mathematical relationship which is used to describe the flow of groundwater through an aquifer. The transient flow of groundwater is described by a form of the diffusion equation, similar to that used in heat transfer to describe the flow of heat in a solid. The steady-state flow of groundwater is described by a form of the Laplace equation, which is a form of potential flow and has analogs in numerous fields.

The United States Environmental Protection Agency (EPA) Storm Water Management Model (SWMM) is a dynamic rainfall–runoff–subsurface runoff simulation model used for single-event to long-term (continuous) simulation of the surface/subsurface hydrology quantity and quality from primarily urban/suburban areas.

<span class="mw-page-title-main">Water balance</span> Looks at how water moves in a closed system

The law of water balance states that the inflows to any water system or area is equal to its outflows plus change in storage during a time interval. In hydrology, a water balance equation can be used to describe the flow of water in and out of a system. A system can be one of several hydrological or water domains, such as a column of soil, a drainage basin, an irrigation area or a city.

In geotechnical engineering, watertable control is the practice of controlling the height of the water table by drainage. Its main applications are in agricultural land and in cities to manage the extensive underground infrastructure that includes the foundations of large buildings, underground transit systems, and extensive utilities.

<span class="mw-page-title-main">Soil salinity control</span> Controlling the problem of soil salinity

Soil salinity control refers to controlling the process and progress of soil salinity to prevent soil degradation by salination and reclamation of already salty (saline) soils. Soil reclamation is also known as soil improvement, rehabilitation, remediation, recuperation, or amelioration.

Well drainage means drainage of agricultural lands by wells. Agricultural land is drained by pumped wells to improve the soils by controlling water table levels and soil salinity.

<span class="mw-page-title-main">SahysMod</span>

SahysMod is a computer program for the prediction of the salinity of soil moisture, groundwater and drainage water, the depth of the watertable, and the drain discharge in irrigated agricultural lands, using different hydrogeologic and aquifer conditions, varying water management options, including the use of ground water for irrigation, and several crop rotation schedules, whereby the spatial variations are accounted for through a network of polygons.

<span class="mw-page-title-main">SaltMod</span> Salinity modelling software

SaltMod is computer program for the prediction of the salinity of soil moisture, groundwater and drainage water, the depth of the watertable, and the drain discharge (hydrology) in irrigated agricultural lands, using different (geo)hydrologic conditions, varying water management options, including the use of ground water for irrigation, and several cropping rotation schedules. The water management options include irrigation, drainage, and the use of subsurface drainage water from pipe drains, ditches or wells for irrigation.

<span class="mw-page-title-main">Hydrus (software)</span> Hydrologic simulation software suite

Hydrus is a suite of Windows-based modeling software that can be used for analysis of water flow, heat and solute transport in variably saturated porous media. HYDRUS suite of software is supported by an interactive graphics-based interface for data-preprocessing, discretization of the soil profile, and graphic presentation of the results. While HYDRUS-1D simulates water flow, solute and heat transport in one-dimension, and is a public domain software, HYDRUS 2D/3D extends the simulation capabilities to the second and third dimensions, and is distributed commercially.

The groundwater energy balance is the energy balance of a groundwater body in terms of incoming hydraulic energy associated with groundwater inflow into the body, energy associated with the outflow, energy conversion into heat due to friction of flow, and the resulting change of energy status and groundwater level.

<span class="mw-page-title-main">Hydrological model</span> Predicting and managing water resources

A hydrologic model is a simplification of a real-world system that aids in understanding, predicting, and managing water resources. Both the flow and quality of water are commonly studied using hydrologic models.

<span class="mw-page-title-main">Agricultural hydrology</span>

Agricultural hydrology is the study of water balance components intervening in agricultural water management, especially in irrigation and drainage.

GSSHA is a two-dimensional, physically based watershed model developed by the Engineer Research and Development Center of the United States Army Corps of Engineers. It simulates surface water and groundwater hydrology, erosion and sediment transport. The GSSHA model is used for hydraulic engineering and research, and is on the Federal Emergency Management Agency (FEMA) list of hydrologic models accepted for use in the national flood insurance program for flood hydrograph estimation. Input is best prepared by the Watershed Modeling System interface, which effectively links the model with geographic information systems (GIS).

The following outline is provided as an overview of and topical guide to hydrology:

<span class="mw-page-title-main">Drainage equation</span> Equation used in drainage design

A drainage equation is an equation describing the relation between depth and spacing of parallel subsurface drains, depth of the watertable, depth and hydraulic conductivity of the soils. It is used in drainage design.

Hydrogeophysics is a cross-disciplinary area of research that uses geophysics to determine parameters and monitor processes for hydrological studies of matters such as water resources, contamination, and ecological studies. The field uses knowledge and researchers from geology, hydrology, physics, geophysics, engineering, statistics, and rock physics. It uses geophysics to provide quantitative information about hydrogeological parameters, using minimally invasive methods. Hydrogeophysics differs from geophysics in its specific uses and methods. Although geophysical knowledge and methods have existed and grown over the last half century for applications in mining and petroleum industries, hydrogeological study sites have different subsurface conditions than those industries. Thus, the geophysical methods for mapping subsurface properties combine with hydrogeology to use proper, accurate methods to map shallow hydrological study sites.

<span class="mw-page-title-main">EnDrain</span> Software

EnDrain is software for the calculation of a subsurface drainage system in agricultural land. The EnDrain program computes the water flow discharged by drains, the hydraulic head losses and the distance between drains, also obtaining the curve described by water-table level. Such calculations are necessary to design a drainage system in the framework of an irrigation system for water table and soil salinity control.