{"id":1136,"date":"2026-05-12T21:46:21","date_gmt":"2026-05-12T21:46:21","guid":{"rendered":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/chapter\/groundwater-flow-equations\/"},"modified":"2026-05-13T16:13:37","modified_gmt":"2026-05-13T16:13:37","slug":"groundwater-flow-equations","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/chapter\/groundwater-flow-equations\/","title":{"raw":"4.4 Groundwater Flow Equations","rendered":"4.4 Groundwater Flow Equations"},"content":{"raw":"<div id=\"groundwater-flow-equations\" class=\"level1\">\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nDarcy\u2019s law is the basic equation that governs the one\u2011dimensional flow of water in a porous medium. So far, we have treated the gradient term as a simple quantity that can be calculated by dividing the difference in head from one point to the next by the distance between the two points (<a href=\"#X8ccfdf9def3efe9523aca1224e0b7c264c65f53\"><span data-custom-style=\"Hyperlink\">Figure\u00a017<\/span><\/a>). Also, we have assumed that wells A and B are along a line oriented with steepest head decline for the point of interest.\r\n\r\n<\/div>\r\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\r\n\r\n<img style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.5in; height: 1.86in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file19-2.png\" alt=\"\" \/>\r\n<div class=\"bc-figcaption figcaption\">\r\n<div data-custom-style=\"Image Caption\">\r\n<div class=\"figcaption-text\">\r\n\r\n<span id=\"X8ccfdf9def3efe9523aca1224e0b7c264c65f53\" class=\"anchor\"><strong>Figure\u00a017\u00a0\u2011\u00a0<\/strong>Cross section showing the overall gradient between two wells.<\/span>\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nIt is not inappropriate to do this in practice; however, it may not be entirely accurate. The water table could look something like <a href=\"#Xb8252a415e1d32378897944ada1e65db1372435\"><span data-custom-style=\"Hyperlink\">Figure\u00a018<\/span><\/a>. In this case, the gradient is smaller near well A and larger near well B.\r\n\r\n<\/div>\r\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\r\n\r\n<img style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.5in; height: 1.86in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file20-2.png\" alt=\"\" \/>\r\n<div class=\"bc-figcaption figcaption\">\r\n<div data-custom-style=\"Image Caption\">\r\n<div class=\"figcaption-text\">\r\n\r\n<span id=\"Xb8252a415e1d32378897944ada1e65db1372435\" class=\"anchor\"><strong>Figure\u00a018\u00a0\u2011\u00a0<\/strong>Cross section showing the water table between two wells.<\/span>\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nThe hydraulic gradient represents the slope of the water table or potentiometric surface at a point, and that slope varies from one point to the next. In practice, we can occasionally ignore the type of variations shown in <a href=\"#Xb8252a415e1d32378897944ada1e65db1372435\"><span data-custom-style=\"Hyperlink\">Figure\u00a018<\/span><\/a> and assume that they are not significant. However, it is important to have a basic understanding of the theory and mathematics underlying the hydraulic gradient so that you can decide when it is appropriate to ignore those variations and when it is important to acknowledge them.\r\n\r\n<\/div>\r\n<div id=\"differential-equations\" class=\"level3\">\r\n<h3>Differential Equations<\/h3>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nNatural, dynamic, three\u2011dimensional systems (like aquifers) contain an infinite number of points, and the conditions in those systems vary from point to point. Groundwater flow is driven by a gradient. at each point. A gradient is a slope at a point (i.e., a rise of a line on a graph divided by the run), which is also known as a derivative. Derivatives are the underlying basis of differential equations. So, we use differential equations to describe and understand groundwater systems.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nA derivative is the instantaneous change in a parameter with respect to some independent variable. In Darcy\u2019s law, the hydraulic gradient is a derivative. <a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a> graphically shows the slope of the water table over progressively smaller intervals (<a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a>a through c) with <a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a>d showing the slope at an infinitesimally small distance (i.e., at a point). We use the symbol [latex]\\mathbf{\\mathrm{\\Delta}}[\/latex] to indicate a finite amount of change in a parameter; for example, [latex]\\mathrm{\\Delta}h[\/latex] refers to a specific change in head over a given lateral distance ([latex]\\mathrm{\\Delta}l)[\/latex].\r\n\r\n<\/div>\r\n<div data-custom-style=\"Captioned Figure\">\r\n\r\n<img style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.38961in; height: 8.25247in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file21-2.png\" alt=\"A diagram of a straight line Description automatically generated with medium confidence\" \/>\r\n\r\n<\/div>\r\n<div data-custom-style=\"Block Text\">\r\n<blockquote>\r\n<div class=\"figcaption-text\">\r\n\r\n<span id=\"X4068acdddafd92e7d11d8318596b42f20b720f6\" class=\"anchor\"><strong>Figure\u00a019\u00a0-<\/strong>\u00a0a), b), c) Illustration of hydraulic gradient as a slope along progressively smaller finite distances ([latex]\\mathrm{\\Delta}h\/\\mathrm{\\Delta}l[\/latex]) and d) as an infinitesimal derivative at a single point ([latex]dh\/dl[\/latex]).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nIn a differential equation, we use the symbol d to signify the infinitesimal change in magnitude of a parameter (e.g., head) with respect to some independent variable (e.g., length along the flow path). So, for Darcy\u2019s law, we are looking at the hydraulic gradient at a single point in a single direction (along the length of the flow path), and the derivative dh\/dl represents the gradient in that direction at that single point.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nSince groundwater flow (i.e., [latex]q[\/latex]) is driven by the gradient, q also varies from one point to the next. So, we can think about the change in flow with respect to length as [latex]\\frac{dq}{dl}[\/latex].\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nAn aquifer, however, exists in three dimensions and contains an infinite number of points, and at each point groundwater is flowing in some direction. Let\u2019s choose a single point in our aquifer and consider how groundwater is flowing at that point (<a href=\"#X37ebab54cd7728df5f234d1188798b3f89e6038\"><span data-custom-style=\"Hyperlink\">Figure\u00a020<\/span><\/a>). At this point groundwater is flowing in some unique direction. We can represent the flow direction and magnitude with a vector (i.e., an arrow), and define a three\u2011dimensional <em>coordinate system<\/em> (the [latex]x,\\ y[\/latex], and [latex]z[\/latex] arrows) to provide a frame of reference. Any vector within a coordinate system can be resolved into <em>component vectors<\/em> in the principal ([latex]x,\\ y[\/latex], and [latex]z[\/latex]) directions. What this means is that the three component vectors add up to the groundwater flow vector.\r\n\r\n<\/div>\r\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\r\n\r\n<img style=\"display: block; margin-left: auto; margin-right: auto; ;width: 5.83333in; height: 1.83115in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file22-2.png\" alt=\"A diagram of a graph Description automatically generated\" \/>\r\n<div class=\"bc-figcaption figcaption\">\r\n<div data-custom-style=\"Image Caption\">\r\n<div class=\"figcaption-text\">\r\n\r\n<span id=\"X37ebab54cd7728df5f234d1188798b3f89e6038\" class=\"anchor\"><strong>Figure\u00a020\u00a0\u2011\u00a0<\/strong>Vector resolved into component vectors in three dimensions.<\/span>\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nEach one of the component vectors represents the flow (i.e., [latex]q[\/latex]) in each of the principal directions. There is a [latex]q[\/latex] in the x direction ([latex]q_{x}[\/latex]), in the y direction ([latex]q_{y}[\/latex]), and in the z direction ([latex]q_{z}[\/latex]). When we have a derivative that has multiple independent variables (i.e., the three directions), we use a partial differential equation, which uses the symbol [latex]\\partial[\/latex] instead of [latex]d[\/latex] to indicate the change in magnitude of a variable with respect to the independent variable in each direction. We can then express that flow vector using partial differentials as [latex]\\left( \\frac{{\\partial q}_{x}}{\\partial x} + \\frac{{\\partial q}_{y}}{\\partial y} + \\frac{{\\partial q}_{z}}{\\partial z} \\right).[\/latex]\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nGenerally, groundwater flow varies with time and we can equate the flow vector to the change in water stored in the system over an infinitesimally small time (i.e., a derivative of water level, h, with time, t) as expressed by the partial differential equation shown as Equation\u00a0<a href=\"#X7fb5e44e4ca99e717e50f51785f28daac264448\"><span data-custom-style=\"Hyperlink\">(34)<\/span><\/a>.\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 8%;\" \/> <col style=\"width: 64%;\" \/> <col style=\"width: 11%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]\\frac{{\\partial q}_{x}}{\\partial x} + \\frac{{\\partial q}_{y}}{\\partial y} + \\frac{{\\partial q}_{z}}{\\partial z} = S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/th>\r\n<th style=\"text-align: right;\"><span id=\"X7fb5e44e4ca99e717e50f51785f28daac264448\" class=\"anchor\">(34)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nwhere (parameter dimensions are dark green font with mass as M, length as L, time as T):\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 12%;\" \/> <col style=\"width: 5%;\" \/> <col style=\"width: 66%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th>[latex]\\mathbf{S}_{\\mathbf{s}}[\/latex]<\/th>\r\n<th style=\"text-align: center;\">=<\/th>\r\n<th>the <em>specific storage<\/em> (L<sup>-1<\/sup>), which is a coefficient that relates changes in hydraulic head to changes in volume of water stored in a unit volume of aquifer. This term is defined in more detail in Chapter 5.<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nDarcy\u2019s law (expressed as the specific discharge) is shown again here in Equation\u00a0<a href=\"#Xa4929b30e780544ce9431da595f4a13a3293bdd\"><span data-custom-style=\"Hyperlink\">(35)<\/span><\/a>.\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 8%;\" \/> <col style=\"width: 64%;\" \/> <col style=\"width: 11%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]q = - K\\frac{dh}{dl}[\/latex]<\/th>\r\n<th style=\"text-align: right;\"><span id=\"Xa4929b30e780544ce9431da595f4a13a3293bdd\" class=\"anchor\">(35)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nSo, we can substitute it for [latex]q[\/latex], to obtain Equation\u00a0<a href=\"#X81fad6ff2f9c70b4db719104d797170ac146ad3\"><span data-custom-style=\"Hyperlink\">(36)<\/span><\/a>.\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 8%;\" \/> <col style=\"width: 64%;\" \/> <col style=\"width: 11%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]\\frac{\\partial\\left( K_{x}\\frac{\\partial h}{\\partial x} \\right)}{\\partial x} + \\frac{\\partial\\left( K_{y}\\frac{\\partial h}{\\partial y} \\right)}{\\partial y} + \\frac{\\partial\\left( K_{z}\\frac{\\partial h}{\\partial z} \\right)}{\\partial z} = S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/th>\r\n<th style=\"text-align: right;\"><span id=\"X81fad6ff2f9c70b4db719104d797170ac146ad3\" class=\"anchor\">(36)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nThis is the primary groundwater flow equation, which describes the flow of groundwater at <em>every point<\/em> in a flow system. The left side of this equation is called the divergence of flux because it represents the amount of flux entering or leaving a point.\r\n\r\n<\/div>\r\n<\/div>\r\n<div id=\"what-does-this-mean-in-conceptual-terms\" class=\"level3\">\r\n<h3>What Does This Mean in Conceptual Terms?<\/h3>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThis is just a more complex version of the continuity equation; each term on the left is:\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\n[latex]I - O[\/latex] (<em>what goes in<\/em> minus <em>what goes out<\/em>)\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nfor each direction (x, y, and z) for a given point in the flow system, and the term on the right is:\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\n[latex]\\mathrm{\\Delta}Storage [\/latex](<em>change in storage<\/em>)\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nfor that point in the system.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nIn a 3\u2011D flow system, there are an infinite number of points, and an infinite number of equations that describe flow through each point. We cannot simply plug a few values into the equation and directly solve it; instead, we must use some indirect method to solve the equation.\r\n\r\n<\/div>\r\n<\/div>\r\n<div id=\"how-do-we-solve-these-equations\" class=\"level3\">\r\n<h3>How Do We Solve These Equations?<\/h3>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nWe use one of several methods to solve the partial differential equations that represent natural systems.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\n<em>Flow nets<\/em> are a graphical representation of a groundwater flow system that can be used to identify flow paths in an aquifer and estimate hydraulic conductivity. More information on constructing and using flow nets can be found in <a href=\"https:\/\/gw-project.org\/books\/graphical-construction-of-groundwater-flow-nets\/\"><span data-custom-style=\"Hyperlink\"><em>Graphical Construction of Groundwater Flow Nets<\/em><\/span><\/a>\u2197 by Poeter and Hsieh (2020), which is available from The Groundwater Project as a free download or can be read online.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nAnalytical equations are solutions of the groundwater flow equation for specific situations, such as a system with a simple geometry. An example of an analytical solution to the 3\u2011D groundwater equation\u2014i.e., Equation\u00a0<a href=\"#X81fad6ff2f9c70b4db719104d797170ac146ad3\"><span data-custom-style=\"Hyperlink\">(36)<\/span><\/a>\u2014is the Theis equation, which solves for radial flow to a pumping well with a specific set of assumptions. The Theis equation is discussed in more detail in Section\u00a07.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nNumerical solutions involve discretizing a system\u2014that is, reducing the infinite number of points\/equations to a finite number and simultaneously solving that set of equations. For numerical solutions, each individual equation is generally dependent on the values for the surrounding equations, often in a non-linear manner; therefore, the equations need to be solved using an iterative method that requires a computer and specialized software. These types of solutions, and the software used to develop them, are generally called <em>computer models<\/em>. Popular software packages used for developing groundwater flow models include MODFLOW and FEFLOW. Computer models are discussed in more detail in Section\u00a010.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nA more detailed discussion of the equations governing groundwater flow can be found in the seventh section of <a href=\"https:\/\/gw-project.org\/books\/hydrogeologic-properties-of-earth-materials-and-principles-of-groundwater-flow\/\"><span data-custom-style=\"Hyperlink\"><em>Hydrogeologic Properties of Earth Materials and Principles of Groundwater Flow<\/em><\/span><\/a>\u2197 by Woessner and Poeter (2020), which is available from The Groundwater Project as a free download or can be read online. In addition, a good introduction to using MODFLOW can be found in <a href=\"https:\/\/gw-project.org\/books\/getting-started-with-modflow\/\"><span data-custom-style=\"Hyperlink\"><em><u>Getting Started with MODFLOW<\/u><\/em><\/span><\/a>\u2197 by Winston (2023), which is also available as a free download from The Groundwater\u00a0Project.\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>","rendered":"<div id=\"groundwater-flow-equations\" class=\"level1\">\n<div data-custom-style=\"First Paragraph\">\n<p>Darcy\u2019s law is the basic equation that governs the one\u2011dimensional flow of water in a porous medium. So far, we have treated the gradient term as a simple quantity that can be calculated by dividing the difference in head from one point to the next by the distance between the two points (<a href=\"#X8ccfdf9def3efe9523aca1224e0b7c264c65f53\"><span data-custom-style=\"Hyperlink\">Figure\u00a017<\/span><\/a>). Also, we have assumed that wells A and B are along a line oriented with steepest head decline for the point of interest.<\/p>\n<\/div>\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.5in; height: 1.86in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file19-2.png\" alt=\"\" \/><\/p>\n<div class=\"bc-figcaption figcaption\">\n<div data-custom-style=\"Image Caption\">\n<div class=\"figcaption-text\">\n<p><span id=\"X8ccfdf9def3efe9523aca1224e0b7c264c65f53\" class=\"anchor\"><strong>Figure\u00a017\u00a0\u2011\u00a0<\/strong>Cross section showing the overall gradient between two wells.<\/span><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>It is not inappropriate to do this in practice; however, it may not be entirely accurate. The water table could look something like <a href=\"#Xb8252a415e1d32378897944ada1e65db1372435\"><span data-custom-style=\"Hyperlink\">Figure\u00a018<\/span><\/a>. In this case, the gradient is smaller near well A and larger near well B.<\/p>\n<\/div>\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.5in; height: 1.86in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file20-2.png\" alt=\"\" \/><\/p>\n<div class=\"bc-figcaption figcaption\">\n<div data-custom-style=\"Image Caption\">\n<div class=\"figcaption-text\">\n<p><span id=\"Xb8252a415e1d32378897944ada1e65db1372435\" class=\"anchor\"><strong>Figure\u00a018\u00a0\u2011\u00a0<\/strong>Cross section showing the water table between two wells.<\/span><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>The hydraulic gradient represents the slope of the water table or potentiometric surface at a point, and that slope varies from one point to the next. In practice, we can occasionally ignore the type of variations shown in <a href=\"#Xb8252a415e1d32378897944ada1e65db1372435\"><span data-custom-style=\"Hyperlink\">Figure\u00a018<\/span><\/a> and assume that they are not significant. However, it is important to have a basic understanding of the theory and mathematics underlying the hydraulic gradient so that you can decide when it is appropriate to ignore those variations and when it is important to acknowledge them.<\/p>\n<\/div>\n<div id=\"differential-equations\" class=\"level3\">\n<h3>Differential Equations<\/h3>\n<div data-custom-style=\"First Paragraph\">\n<p>Natural, dynamic, three\u2011dimensional systems (like aquifers) contain an infinite number of points, and the conditions in those systems vary from point to point. Groundwater flow is driven by a gradient. at each point. A gradient is a slope at a point (i.e., a rise of a line on a graph divided by the run), which is also known as a derivative. Derivatives are the underlying basis of differential equations. So, we use differential equations to describe and understand groundwater systems.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>A derivative is the instantaneous change in a parameter with respect to some independent variable. In Darcy\u2019s law, the hydraulic gradient is a derivative. <a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a> graphically shows the slope of the water table over progressively smaller intervals (<a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a>a through c) with <a href=\"#X4068acdddafd92e7d11d8318596b42f20b720f6\"><span data-custom-style=\"Hyperlink\">Figure\u00a019<\/span><\/a>d showing the slope at an infinitesimally small distance (i.e., at a point). We use the symbol [latex]\\mathbf{\\mathrm{\\Delta}}[\/latex] to indicate a finite amount of change in a parameter; for example, [latex]\\mathrm{\\Delta}h[\/latex] refers to a specific change in head over a given lateral distance ([latex]\\mathrm{\\Delta}l)[\/latex].<\/p>\n<\/div>\n<div data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 4.38961in; height: 8.25247in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file21-2.png\" alt=\"A diagram of a straight line Description automatically generated with medium confidence\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><span id=\"X4068acdddafd92e7d11d8318596b42f20b720f6\" class=\"anchor\"><strong>Figure\u00a019\u00a0&#8211;<\/strong>\u00a0a), b), c) Illustration of hydraulic gradient as a slope along progressively smaller finite distances ([latex]\\mathrm{\\Delta}h\/\\mathrm{\\Delta}l[\/latex]) and d) as an infinitesimal derivative at a single point ([latex]dh\/dl[\/latex]).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<div data-custom-style=\"First Paragraph\">\n<p>In a differential equation, we use the symbol d to signify the infinitesimal change in magnitude of a parameter (e.g., head) with respect to some independent variable (e.g., length along the flow path). So, for Darcy\u2019s law, we are looking at the hydraulic gradient at a single point in a single direction (along the length of the flow path), and the derivative dh\/dl represents the gradient in that direction at that single point.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Since groundwater flow (i.e., [latex]q[\/latex]) is driven by the gradient, q also varies from one point to the next. So, we can think about the change in flow with respect to length as [latex]\\frac{dq}{dl}[\/latex].<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>An aquifer, however, exists in three dimensions and contains an infinite number of points, and at each point groundwater is flowing in some direction. Let\u2019s choose a single point in our aquifer and consider how groundwater is flowing at that point (<a href=\"#X37ebab54cd7728df5f234d1188798b3f89e6038\"><span data-custom-style=\"Hyperlink\">Figure\u00a020<\/span><\/a>). At this point groundwater is flowing in some unique direction. We can represent the flow direction and magnitude with a vector (i.e., an arrow), and define a three\u2011dimensional <em>coordinate system<\/em> (the [latex]x,\\ y[\/latex], and [latex]z[\/latex] arrows) to provide a frame of reference. Any vector within a coordinate system can be resolved into <em>component vectors<\/em> in the principal ([latex]x,\\ y[\/latex], and [latex]z[\/latex]) directions. What this means is that the three component vectors add up to the groundwater flow vector.<\/p>\n<\/div>\n<div class=\"bc-figure figure\" data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 5.83333in; height: 1.83115in;\" src=\"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-content\/uploads\/sites\/42\/2026\/05\/file22-2.png\" alt=\"A diagram of a graph Description automatically generated\" \/><\/p>\n<div class=\"bc-figcaption figcaption\">\n<div data-custom-style=\"Image Caption\">\n<div class=\"figcaption-text\">\n<p><span id=\"X37ebab54cd7728df5f234d1188798b3f89e6038\" class=\"anchor\"><strong>Figure\u00a020\u00a0\u2011\u00a0<\/strong>Vector resolved into component vectors in three dimensions.<\/span><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Each one of the component vectors represents the flow (i.e., [latex]q[\/latex]) in each of the principal directions. There is a [latex]q[\/latex] in the x direction ([latex]q_{x}[\/latex]), in the y direction ([latex]q_{y}[\/latex]), and in the z direction ([latex]q_{z}[\/latex]). When we have a derivative that has multiple independent variables (i.e., the three directions), we use a partial differential equation, which uses the symbol [latex]\\partial[\/latex] instead of [latex]d[\/latex] to indicate the change in magnitude of a variable with respect to the independent variable in each direction. We can then express that flow vector using partial differentials as [latex]\\left( \\frac{{\\partial q}_{x}}{\\partial x} + \\frac{{\\partial q}_{y}}{\\partial y} + \\frac{{\\partial q}_{z}}{\\partial z} \\right).[\/latex]<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Generally, groundwater flow varies with time and we can equate the flow vector to the change in water stored in the system over an infinitesimally small time (i.e., a derivative of water level, h, with time, t) as expressed by the partial differential equation shown as Equation\u00a0<a href=\"#X7fb5e44e4ca99e717e50f51785f28daac264448\"><span data-custom-style=\"Hyperlink\">(34)<\/span><\/a>.<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 8%;\" \/>\n<col style=\"width: 64%;\" \/>\n<col style=\"width: 11%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]\\frac{{\\partial q}_{x}}{\\partial x} + \\frac{{\\partial q}_{y}}{\\partial y} + \\frac{{\\partial q}_{z}}{\\partial z} = S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/th>\n<th style=\"text-align: right;\"><span id=\"X7fb5e44e4ca99e717e50f51785f28daac264448\" class=\"anchor\">(34)<\/span><\/th>\n<\/tr>\n<\/thead>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>where (parameter dimensions are dark green font with mass as M, length as L, time as T):<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 12%;\" \/>\n<col style=\"width: 5%;\" \/>\n<col style=\"width: 66%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th>[latex]\\mathbf{S}_{\\mathbf{s}}[\/latex]<\/th>\n<th style=\"text-align: center;\">=<\/th>\n<th>the <em>specific storage<\/em> (L<sup>-1<\/sup>), which is a coefficient that relates changes in hydraulic head to changes in volume of water stored in a unit volume of aquifer. This term is defined in more detail in Chapter 5.<\/th>\n<\/tr>\n<\/thead>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>Darcy\u2019s law (expressed as the specific discharge) is shown again here in Equation\u00a0<a href=\"#Xa4929b30e780544ce9431da595f4a13a3293bdd\"><span data-custom-style=\"Hyperlink\">(35)<\/span><\/a>.<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 8%;\" \/>\n<col style=\"width: 64%;\" \/>\n<col style=\"width: 11%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]q = - K\\frac{dh}{dl}[\/latex]<\/th>\n<th style=\"text-align: right;\"><span id=\"Xa4929b30e780544ce9431da595f4a13a3293bdd\" class=\"anchor\">(35)<\/span><\/th>\n<\/tr>\n<\/thead>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>So, we can substitute it for [latex]q[\/latex], to obtain Equation\u00a0<a href=\"#X81fad6ff2f9c70b4db719104d797170ac146ad3\"><span data-custom-style=\"Hyperlink\">(36)<\/span><\/a>.<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 8%;\" \/>\n<col style=\"width: 64%;\" \/>\n<col style=\"width: 11%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]\\frac{\\partial\\left( K_{x}\\frac{\\partial h}{\\partial x} \\right)}{\\partial x} + \\frac{\\partial\\left( K_{y}\\frac{\\partial h}{\\partial y} \\right)}{\\partial y} + \\frac{\\partial\\left( K_{z}\\frac{\\partial h}{\\partial z} \\right)}{\\partial z} = S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/th>\n<th style=\"text-align: right;\"><span id=\"X81fad6ff2f9c70b4db719104d797170ac146ad3\" class=\"anchor\">(36)<\/span><\/th>\n<\/tr>\n<\/thead>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>This is the primary groundwater flow equation, which describes the flow of groundwater at <em>every point<\/em> in a flow system. The left side of this equation is called the divergence of flux because it represents the amount of flux entering or leaving a point.<\/p>\n<\/div>\n<\/div>\n<div id=\"what-does-this-mean-in-conceptual-terms\" class=\"level3\">\n<h3>What Does This Mean in Conceptual Terms?<\/h3>\n<div data-custom-style=\"First Paragraph\">\n<p>This is just a more complex version of the continuity equation; each term on the left is:<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>[latex]I - O[\/latex] (<em>what goes in<\/em> minus <em>what goes out<\/em>)<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>for each direction (x, y, and z) for a given point in the flow system, and the term on the right is:<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>[latex]\\mathrm{\\Delta}Storage[\/latex](<em>change in storage<\/em>)<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>for that point in the system.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>In a 3\u2011D flow system, there are an infinite number of points, and an infinite number of equations that describe flow through each point. We cannot simply plug a few values into the equation and directly solve it; instead, we must use some indirect method to solve the equation.<\/p>\n<\/div>\n<\/div>\n<div id=\"how-do-we-solve-these-equations\" class=\"level3\">\n<h3>How Do We Solve These Equations?<\/h3>\n<div data-custom-style=\"First Paragraph\">\n<p>We use one of several methods to solve the partial differential equations that represent natural systems.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p><em>Flow nets<\/em> are a graphical representation of a groundwater flow system that can be used to identify flow paths in an aquifer and estimate hydraulic conductivity. More information on constructing and using flow nets can be found in <a href=\"https:\/\/gw-project.org\/books\/graphical-construction-of-groundwater-flow-nets\/\"><span data-custom-style=\"Hyperlink\"><em>Graphical Construction of Groundwater Flow Nets<\/em><\/span><\/a>\u2197 by Poeter and Hsieh (2020), which is available from The Groundwater Project as a free download or can be read online.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Analytical equations are solutions of the groundwater flow equation for specific situations, such as a system with a simple geometry. An example of an analytical solution to the 3\u2011D groundwater equation\u2014i.e., Equation\u00a0<a href=\"#X81fad6ff2f9c70b4db719104d797170ac146ad3\"><span data-custom-style=\"Hyperlink\">(36)<\/span><\/a>\u2014is the Theis equation, which solves for radial flow to a pumping well with a specific set of assumptions. The Theis equation is discussed in more detail in Section\u00a07.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Numerical solutions involve discretizing a system\u2014that is, reducing the infinite number of points\/equations to a finite number and simultaneously solving that set of equations. For numerical solutions, each individual equation is generally dependent on the values for the surrounding equations, often in a non-linear manner; therefore, the equations need to be solved using an iterative method that requires a computer and specialized software. These types of solutions, and the software used to develop them, are generally called <em>computer models<\/em>. Popular software packages used for developing groundwater flow models include MODFLOW and FEFLOW. Computer models are discussed in more detail in Section\u00a010.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>A more detailed discussion of the equations governing groundwater flow can be found in the seventh section of <a href=\"https:\/\/gw-project.org\/books\/hydrogeologic-properties-of-earth-materials-and-principles-of-groundwater-flow\/\"><span data-custom-style=\"Hyperlink\"><em>Hydrogeologic Properties of Earth Materials and Principles of Groundwater Flow<\/em><\/span><\/a>\u2197 by Woessner and Poeter (2020), which is available from The Groundwater Project as a free download or can be read online. In addition, a good introduction to using MODFLOW can be found in <a href=\"https:\/\/gw-project.org\/books\/getting-started-with-modflow\/\"><span data-custom-style=\"Hyperlink\"><em><u>Getting Started with MODFLOW<\/u><\/em><\/span><\/a>\u2197 by Winston (2023), which is also available as a free download from The Groundwater\u00a0Project.<\/p>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"author":8,"menu_order":13,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-1136","chapter","type-chapter","status-publish","hentry"],"part":1477,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapters\/1136","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/wp\/v2\/users\/8"}],"version-history":[{"count":3,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapters\/1136\/revisions"}],"predecessor-version":[{"id":1851,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapters\/1136\/revisions\/1851"}],"part":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/parts\/1477"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapters\/1136\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/wp\/v2\/media?parent=1136"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/pressbooks\/v2\/chapter-type?post=1136"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/wp\/v2\/contributor?post=1136"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/wp-json\/wp\/v2\/license?post=1136"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}