{"id":1150,"date":"2026-05-12T21:46:24","date_gmt":"2026-05-12T21:46:24","guid":{"rendered":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/chapter\/steady-state-versus-transient-groundwater-flow\/"},"modified":"2026-05-13T16:22:36","modified_gmt":"2026-05-13T16:22:36","slug":"steady-state-versus-transient-groundwater-flow","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/basic-hydrogeology-an-introduction-to-the-fundamentals-of-groundwater-science\/chapter\/steady-state-versus-transient-groundwater-flow\/","title":{"raw":"6.2 Steady State versus Transient Groundwater Flow","rendered":"6.2 Steady State versus Transient Groundwater Flow"},"content":{"raw":"
\r\n
\r\n

Steady\u2011State Flow<\/h3>\r\n
\r\n\r\nConsider the continuity equation; otherwise known as the law of mass conservation<\/em>, or the continuity principle<\/em>:\r\n\r\n<\/div>\r\n
\r\n
[latex]I - O[\/latex] (what goes in minus what goes out<\/em>)<\/blockquote>\r\n<\/div>\r\n
\r\n
[latex]= \\mathrm{\\Delta}Storage [\/latex](change in storage in the system<\/em>) \/<\/strong> [latex]\\mathrm{\\Delta}t [\/latex](change in time<\/em>).<\/blockquote>\r\n<\/div>\r\n
\r\n\r\nThis principle states that mass is conserved in a system, and any change of mass entering a system must be balanced by a corresponding change in either the mass leaving the system, or the mass stored in the system. In other words, if the mass going in equals the mass going out, then the mass in storage will not change, and [latex]\\mathrm{\\Delta}S\/\\mathrm{\\Delta}t[\/latex] will equal zero. In the case where [latex]I[\/latex] is equal to [latex]O[\/latex], and [latex]\\mathrm{\\Delta}Storage\/\\mathrm{\\Delta}t[\/latex] equals zero, we say that the system is at steady state. If the system is at steady state, heads at any given point in the aquifer do not change through time, the [latex]\\partial h[\/latex] on the right side of the equation equals zero, as shown by Equation\u00a0(46)<\/span><\/a> and Equation\u00a0(47)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n\r\n\r\n
<\/th>\r\n[latex]\\frac{{\\partial q}_{x}}{\\partial x} + \\frac{{\\partial q}_{y}}{\\partial y} + \\frac{{\\partial q}_{z}}{\\partial z} = 0[\/latex]<\/th>\r\n(46)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n
<\/td>\r\n[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} = 0[\/latex]<\/td>\r\n(47)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n\r\nIf we assume the aquifer is homogeneous\u2014that is, the permeability at one point is the same as at any other point\u2014then the hydraulic conductivity ([latex]K[\/latex]) will not change with respect to location, and we can remove [latex]K_{x},\\ K_{y},[\/latex] and [latex]K_{z}[\/latex] from inside the differentials as shown in Equation\u00a0(48)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n
<\/th>\r\n[latex]K_{x}\\frac{\\partial\\left( \\frac{\\partial h}{\\partial x} \\right)}{\\partial x} + K_{y}\\frac{\\partial\\left( \\frac{\\partial h}{\\partial y} \\right)}{\\partial y} + K_{z}\\frac{\\partial\\left( \\frac{\\partial h}{\\partial z} \\right)}{\\partial z} = 0[\/latex]<\/th>\r\n(48)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n
\r\n\r\nIf we also assume that the permeability of our aquifer is isotropic<\/em>, then the hydraulic conductivity ([latex]K[\/latex]) at any point in the aquifer will be the same in all directions, and [latex]K_{x}\\ = \\ K_{y}\\ = \\ K_{z}\\ = \\ K[\/latex]. In that case, we can factor out [latex]K[\/latex] and rewrite it as shown by Equation\u00a0(49)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n
<\/th>\r\n[latex]K\\left( \\frac{\\partial\\left( \\frac{\\partial h}{\\partial x} \\right)}{\\partial x} + \\frac{\\partial\\left( \\frac{\\partial h}{\\partial y} \\right)}{\\partial y} + \\frac{\\partial\\left( \\frac{\\partial h}{\\partial z} \\right)}{\\partial z} \\right) = 0[\/latex]<\/th>\r\n(49)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n
\r\n\r\nAlso, we can use a simpler notation for the derivatives and rewrite it as shown by Equation\u00a0(50)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n
<\/th>\r\n[latex]K\\left( \\frac{\\partial^{2}h}{{\\partial x}^{2}} + \\frac{\\partial^{2}h}{{\\partial y}^{2}} + \\frac{\\partial^{2}h}{{\\partial z}^{2}} \\right) = 0[\/latex]<\/th>\r\n(50)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n
\r\n\r\nEach one of the derivatives in the parentheses is a 2nd order partial differential. This simply means that it isn\u2019t only that head changes with distance, but the change in head with distance changes with distance (or, in other words, the gradient changes with distance). If we then divide both sides by [latex]K[\/latex], the equation reduces to Equation\u00a0(51)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n
<\/th>\r\n[latex]\\frac{\\partial^{2}h}{{\\partial x}^{2}} + \\frac{\\partial^{2}h}{{\\partial y}^{2}} + \\frac{\\partial^{2}h}{{\\partial z}^{2}} = 0[\/latex]<\/th>\r\n(51)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<\/table>\r\n
\r\n\r\nEquation\u00a0(51)<\/span><\/a> is called the Laplace equation<\/em>. This equation describes the steady\u2011state flow of groundwater (i.e., no changes in storage) in a homogeneous and isotropic aquifer. As we can see, the assumptions of homogeneity and isotropy greatly simplify the equation.\r\n\r\n<\/div>\r\n<\/div>\r\n
\r\n

Transient Flow<\/h3>\r\n
\r\n\r\nSteady state means that the amount going in equals the amount going out. In groundwater systems, this means that the amount recharging the system\u2014infiltration of rain, loss from streams, leakage from other formations\u2014is equal to the amount discharging the system\u2014spring flow, discharge to streams, pumping by wells. In field settings, recharge is never exactly equal to discharge\u2014although, in some cases, they may be pretty close over the long term, such that it is often appropriate to assume a system is in a steady state when dealing with regional flow systems. That is discussed more in subsequent sections of this book. Aquifers are in a constant state of dynamic change; recharge varies, discharge varies, the volume in storage varies, and head varies. If it rains more in one year than it usually does, recharge increases and heads increase. If we pump a water from a well, we remove water from storage and heads decline. These scenarios depict transient flow conditions\u2014conditions where heads and the volume of water in storage are changing through time.\r\n\r\n<\/div>\r\n
\r\n\r\nThe key point to consider here is time. Continuously pumping water from a well will result in ongoing changes in head, and we as hydrogeologists are interested in knowing things like how much we can pump from a well in a day, how fast heads will decline when we pump, and how long it will take for heads to return to their original level if we stop pumping. The right side of the equation essentially represents this aspect of aquifers, so the right side is zero for steady state flow conditions.\r\n\r\n<\/div>\r\n
\r\n\r\nWhen flow varies with time, conditions are transient and the right side of Equation (50)<\/span><\/a> is replaced with Equation\u00a0(52)<\/span><\/a> which describes the change in storage in the aquifer with time, becoming the transient groundwater flow equation shown in Equation (53)<\/span><\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n\r\n\r\n
<\/th>\r\n[latex]\\frac{\\mathrm{\\Delta}Storage}{t} \\approx S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/th>\r\n(52)<\/span><\/th>\r\n<\/tr>\r\n<\/thead>\r\n
<\/td>\r\n[latex]K\\left( \\frac{\\partial^{2}h}{{\\partial x}^{2}} + \\frac{\\partial^{2}h}{{\\partial y}^{2}} + \\frac{\\partial^{2}h}{{\\partial z}^{2}} \\right) = S_{s}\\frac{\\partial h}{\\partial t}[\/latex]<\/td>\r\n(53)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n\r\nIt is useful to remember that specific storage is a coefficient, just like [latex]K[\/latex] is a coefficient. A coefficient simply relates two quantities to one another:\r\n\r\n<\/div>\r\n