{"id":1678,"date":"2026-02-14T04:24:04","date_gmt":"2026-02-14T04:24:04","guid":{"rendered":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/ch044-xhtml\/"},"modified":"2026-08-01T10:47:31","modified_gmt":"2026-08-01T10:47:31","slug":"hantush-equation-for-a-leaky-confined-system-with-water-released-from-confining-bed-storage","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/hantush-equation-for-a-leaky-confined-system-with-water-released-from-confining-bed-storage\/","title":{"raw":"9.3 Hantush Equation for a Leaky Confined System with Water Released from Confining Bed Storage","rendered":"9.3 Hantush Equation for a Leaky Confined System with Water Released from Confining Bed Storage"},"content":{"raw":"<div id=\"hantush-equation-for-a-leaky-confined-system-with-water-released-from-confining-bed-storage\" class=\"level1\">\r\n<div data-custom-style=\"First Paragraph\">\r\n<div data-custom-style=\"First Paragraph\">\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nWhen pumping a confined aquifer with a permeable confining bed contributing water to the pumped aquifer, the reduction in head in the aquifer creates a hydraulic gradient in the confining bed toward the aquifer. As heads in the confining bed are lowered, water is released from storage and flows to the confined aquifer. If the confining bed is relativity thin, the gradient change reaches the top of the bed relatively quickly and then stabilizes such that release of water from storage within the confining bed ceases. If there is an aquifer on the other side of the confining bed, it provides a source of water that flows through the confining bed to the pumped aquifer. Hantush (1960) developed a set of analytical solutions that account for leakage of water from aquitard storage into a confined aquifer as it is pumped. His work also addresses leakage into a confined system from vertical flow though confining beds as described by Hantush\u2011Jacob and illustrated in <a href=\"\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/formulation-of-equations-to-address-leaky-confined-conditions\/#X57b0920e91bd230bfa3e910c632b20097e16b1b\">Figure 41<\/a>. Conceptual models of such leakage are illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>.\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: 5.4956in; height: 4.37662in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file53-1.jpg\" alt=\"\" \/>\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=\"X6b6c2ec720b85923f014beffc01d41c6be08513\" class=\"anchor\"><strong>Figure 49 \u2011<\/strong> Conceptual models that can be represented by the Hantush equation (1960).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"Block Text\">\r\n<blockquote>a) Drawdown in a confined aquifer that receives water from storage in one or two confining beds when the beds overlying and underlying the confining beds are not permeable.<\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"Block Text\">\r\n<blockquote>b) Drawdown in a confined aquifer that receives water from storage in overlying and\/or underlying confining beds, as well as leakage of water through the confining bed from a highly transmissive unit such as an unconfined water table aquifer.<\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"Block Text\">\r\n<blockquote>c) Drawdown in a confined aquifer receiving water via direct flow through a confining bed that has negligible storage from a highly transmissive aquifer above and\/or below the confining bed.<\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThe general analytical solution for leaky conditions given by Hantush (1960) is presented in Equation\u00a0<a href=\"#X15cd29ce3b49460451a071492929753350a452b\">(55)<\/a>. The different conditions shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a> can be represented using variations of this equation. Similar to the Theis well function for drawdown when pumping a fully confined water bearing unit, the Hantush integral includes an infinite series and can be represented using an approximation as shown in Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a>. The Hantush well function uses the parameter <em>b<\/em> to account for the relative properties of the aquifer and its overlying and underlying layers. Drawdown in a confined unit with one or two adjacent confining unit(s) that release(s) water from storage that is overlain or underlain by an impermeable unit is calculated as shown in Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a>.\r\n\r\n<\/div>\r\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 75%;\" \/> <col style=\"width: 15%;\" \/><\/colgroup>\r\n<tbody>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ \\int_{u}^{\\infty}{\\ \\frac{1}{y}\\ exp\\ \\left\\{ - y \\right\\}\\ erfc\\left( \\frac{\\frac{\\beta}{\\sqrt{u}}}{\\sqrt{y(y - u}\\ )} \\right)}dy[\/latex]<\/th>\r\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X15cd29ce3b49460451a071492929753350a452b\" class=\"anchor\">(55)<\/span><\/th>\r\n<\/tr>\r\n<tr>\r\n<td><\/td>\r\n<td style=\"text-align: center;\">[latex]H(u,\\beta) = W(u) - \\ \\frac{4\\beta}{(\\pi{u)}^{0.5}}\\ \\left\\lbrack 0.2577 + 0.6631\\exp\\left( \\frac{- u}{2}\\ \\right) \\right\\rbrack[\/latex]<\/td>\r\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\" class=\"anchor\">(56)<\/span><\/td>\r\n<\/tr>\r\n<tr>\r\n<td><\/td>\r\n<td style=\"text-align: center;\">[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ H(u,\\beta)[\/latex]<\/td>\r\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\" class=\"anchor\">(57)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nwhere:\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 14%;\" \/> <col style=\"width: 3%;\" \/> <col style=\"width: 67%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th style=\"text-align: right;\"><em>s<\/em>(<em>r<\/em>,<em>t<\/em>)<\/th>\r\n<th style=\"text-align: center;\">=<\/th>\r\n<th>drawdown at a radial distance <em>r<\/em> and time <em>t<\/em> (L)<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>radial distance to the observation well (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>t<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>time (T)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>Q<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>constant well discharge rate (L<sup>3<\/sup>T<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>transmissivity of the confined pumped unit (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>y<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>variable of integration<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\">erfc<\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>complementary error function (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>\u03b2<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>represents the relative transmissibly and storage of the pumped unit and the contributing overlying and\/or underlying units as shown in Equations\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> and <a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>u<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td><em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\">H(<em>u<\/em>,<em>\u03b2<\/em>)<\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>Hantush well function integral (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\">W(<em>u<\/em>)<\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>Theis well function (dimensionless)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nWhen storage comes from only one confining bed, <em>b<\/em> is described by Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a>. This represents the case illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a with one aquitard\u2014which can be above or below the pumped unit. In this case the other aquitard shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a is impermeable. For the setting with two confining beds, <em>b<\/em> includes terms for both aquitards and is described by Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a>.\r\n\r\n<\/div>\r\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 75%;\" \/> <col style=\"width: 15%;\" \/><\/colgroup>\r\n<tbody>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]\\beta = \\frac{r}{4\\ }\\sqrt{\\frac{K'S'}{b'TS}\\ }[\/latex]<\/th>\r\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X35b59d96334da087626420f015e0710c839455d\" class=\"anchor\">(58)<\/span><\/th>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 75%;\" \/> <col style=\"width: 15%;\" \/><\/colgroup>\r\n<tbody>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]\\beta = \\frac{r}{4\\ }\\left( \\sqrt{\\frac{K'S'}{b'TS}\\ }\\ + \\sqrt{\\frac{K\"S\"}{b^{\"}TS}} \\right)[\/latex]<\/th>\r\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X7ef8f34ca870785ea61794264c127ffbbce8681\" class=\"anchor\">(59)<\/span><\/th>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nwhere:\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 14%;\" \/> <col style=\"width: 3%;\" \/> <col style=\"width: 67%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th style=\"text-align: right;\"><em>r<\/em><\/th>\r\n<th style=\"text-align: left;\">=<\/th>\r\n<th>radial distance to the pumping well (L)<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>K\u2019<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>vertical hydraulic conductivity of the first confining unit (LT<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>K\u201d<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>vertical hydraulic conductivity of the second confining unit (LT<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S\u2019<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>storativity of the first confining unit (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S\u201d<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>storativity of the second confining unit (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>b\u2019<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>thickness of the first confining unit (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>b\u201d<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>thickness of the second confining unit (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>transmissivity of the confined unit being pumped (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\r\n<td style=\"text-align: left;\">=<\/td>\r\n<td>storativity of the confined unit being pumped (dimensionless)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nDrawdown in the confined aquifer can be calculated with Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> by using Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a> to calculate the Hantush well function (H(<em>u<\/em>,<em>\u03b2<\/em>)) or reading\/interpolating the value from a table (<a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a>) or type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>).\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: 2.20962in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file54-1.jpg\" 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=\"X32295cc7f397b2f30defdbe73d897d6eb898586\" class=\"anchor\"><strong>Figure 50 \u2011<\/strong> Tables of the Hantush well function for variables <em>u<\/em> and <em>\u03b2<\/em> (from Fetter (2001)). Type curves of H(<em>u<\/em>,<em>\u03b2<\/em>) versus 1\/<em>u<\/em> are illustrated in <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a><\/span>.\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\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: 5.8in; height: 3.52in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file55-1.jpg\" alt=\"\" \/>\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=\"X21bd8ced2a72427e089642240b99c2f3bd335eb\" class=\"anchor\"><strong>Figure 51 \u2011<\/strong> Type curves of H(<em>u<\/em>,<em><u>\u03b2<\/u><\/em>) versus 1\/<em>u<\/em> plotted on log\u2011log axes (modified from Reed (1980); Kruseman &amp; de RIdder (2000); Kasenow (2001)).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThe table of values (<a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a>) and plotted type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>) represent a subset of H(<em>u<\/em>,<em>\u03b2<\/em>) values that can be generated by substituting combinations of <em>u<\/em> and <em>\u03b2<\/em> into Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a> or by using mathematical software packages to solve the integral portion of Equation\u00a0<a href=\"#X15cd29ce3b49460451a071492929753350a452b\">(55)<\/a>. The required values are commonly interpolated from the table or curves presented in <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> and <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a> respectively.\r\n\r\n<\/div>\r\n<div id=\"using-the-hantush-equation-to-predict-drawdown-in-leaky-confined-units-with-water-released-from-confining-bed-storage\" class=\"level3\">\r\n<h3>Using the Hantush Equation to Predict Drawdown in Leaky Confined Units with Water Released from Confining Bed Storage<\/h3>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nWhen the hydraulic properties of the sequence of geologic units as illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a are known, the Hantush equation can be used to approximate drawdown at an observation point at any time <em>t<\/em> after the pumping commenced.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\n<strong>Example<\/strong>\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nA fully penetrating well pumps at 500 m<sup>3<\/sup>\/d for 0.2\u00a0d from a 30 m thick leaky confined sand rich aquifer that is overlain by a silt rich aquitard that is 20 m thick. The confining units are overlain and underlain by impermeable units (<a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a). The observation well is located 53 m from the pumping well. The transmissivity of the confined aquifer is 300 m<sup>2<\/sup>\/d and the storativity is 0.00003. The vertical hydraulic conductivity of the confining bed is 0.3 m\/d and the storativity of the confining bed is 0.00001. Calculate the drawdown in the observation well after 0.2 d of pumping.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nUsing the definition of <em>u<\/em> (<em>u<\/em>=<em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)) along with Equations <a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> and <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a>, and the information listed above, H(<em>u<\/em>,<em>b<\/em>) is obtained by first computing <em>u<\/em> and <em>\uf062<\/em>, and then using <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> to obtain the Hantush well function.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]u = \\frac{r^{2}S}{4Tt} = \\frac{{(53\\ m)}^{2}\\ 0.00003}{4\\ \\left( 300\\frac{m^{2}}{d} \\right)0.2\\ d} = 0.0004[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]\\beta = \\frac{r}{4\\ }\\sqrt{\\frac{K'S'}{b'TS}\\ } = \\frac{53\\ m}{4\\ }\\sqrt{\\frac{0.3\\frac{m}{d}(0.00001)}{20\\ m\\ \\left( 300\\frac{m^{2}}{d} \\right)\\ 0.00003}} = 0.05[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nWith <em>u<\/em> = 0.0004 and <em>b<\/em> = 0.05, H(<em>u<\/em>,<em>b<\/em>) is obtained by interpolation from <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> or <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>. With <em>u<\/em> = 0.0004, <em>b<\/em> = 0.05 then H(<em>u<\/em>,<em>b<\/em>) is about 6.52.\r\n\r\n<\/div>\r\n<table style=\"width: 37%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 8%;\" \/> <col style=\"width: 18%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th><em>u<\/em><\/th>\r\n<th><em>\u03b2<\/em><\/th>\r\n<th>H(<em>u<\/em>,<em>\u03b2<\/em>)<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td>0.0001<\/td>\r\n<td>0.05<\/td>\r\n<td>7.1771<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>0.0004<\/td>\r\n<td>0.05<\/td>\r\n<td>6.52 (interpolated)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>0.0005<\/td>\r\n<td>0.05<\/td>\r\n<td>6.3523<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Captioned Figure\">\r\n\r\n<img style=\"display: block; margin-left: auto; margin-right: auto; ;width: 5.8in; height: 4.13in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file56-1.jpg\" alt=\"\" \/>\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<strong>Figure 52 \u2011<\/strong> Hantush type curves plotted on log\u2011log scales. The dashed vertical blue line is the computed value of 1\/<em>u<\/em> (2,500) for the example problem. The red line is the interpolated value of <em>\u03b2<\/em>, about 0.05, and the black dot dashed line is the interpolated value of H(<em>u<\/em>,<em>\u03b2<\/em>), about 5.3. This value is lower than derived from the tables. The close plotting of the type curves made interpretation difficult for the designated <em>\u03b2\u00a0<\/em>value (modified from Kruseman &amp; de RIdder, 2000; Kasenow, 2001).\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThen Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> provides the value of drawdown.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]s(r,t) = \\frac{500\\frac{m^{3}}{d}}{4\\ (3.14)\\ 300\\frac{m^{2}}{d}}\\ 6.52 = 0.87\\ m[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<div id=\"hantush-curve-matching-method-to-compute-t-and-s-from-a-pumping-test-in-leaky-confined-unit-with-aquitard-storage\" class=\"level3\">\r\n<h3>Hantush Curve Matching Method to Compute T and S from a Pumping Test in Leaky Confined Unit with Aquitard Storage<\/h3>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThe Hantush equation can be used to determine <em>T<\/em> and <em>S<\/em> from time\u2011drawdown data collected in an observation well when a leaky confined aquifer is pumped and water is released from aquitard storage. The hydraulic conductivity and storage properties of adjacent confining units can also be derived. Standard curve\u2011matching techniques described in Section 8.3.1 can be used.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nAfter data from a pumping test has been collected, <em>T<\/em> and <em>S<\/em> for a leaky confined aquifer with the release of water from confining bed storage is estimated by first plotting drawdown versus time on a log\u2011log graph at the same scale as the type curves as shown in (<a href=\"#X87f36bd97f7d73b20af2133215aeb032eb87b68\">Figure 53<\/a>).\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: 4in; height: 2.44in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file57-1.jpg\" alt=\"\" \/>\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=\"X87f36bd97f7d73b20af2133215aeb032eb87b68\" class=\"anchor\"><strong>Figure 53 \u2011<\/strong> Log\u2011log plot of the response of an observation well to the pumping of a well in a leaky confined aquifer with water being released from confining bed storage.<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<div data-custom-style=\"First Paragraph\">\r\n\r\nThe data graph (<a href=\"#X87f36bd97f7d73b20af2133215aeb032eb87b68\">Figure 53<\/a>) is overlain on the type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>). Keeping the axes parallel, the field data graph and type curve graph are overlain and shifted until the curves match (<a href=\"#X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\">Figure 54<\/a>). A match point is selected from the overlapping portion of the two plots and values of H(<em>u<\/em>,<em>\u03b2<\/em>), 1\/<em>u<\/em>, <em>s<\/em>, <em>t<\/em> and <em>\u03b2<\/em> are obtained. Estimates of <em>T<\/em> and <em>S<\/em> for the leaky confined unit as well as its hydraulic conductivity, and, at early times, storativity values of the confining beds can be estimated using Equations <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a> and <a href=\"#X60f5ece6500904699333d9dcddeea65ef1620a1\">(61)<\/a> which are obtained from rearranging Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> and the definition of <em>u<\/em> (<em>u<\/em>=<em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)), respectively. Early time refers to the period after pumping begins when the confining unit is releasing water from storage. Early time is less than <em>b\u2019S\u2019\/10K\u2019<\/em>.\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: 5.28182in; height: 4.07756in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file58-1.jpg\" alt=\"\" \/>\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=\"X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\" class=\"anchor\"><strong>Figure 54 \u2011<\/strong> Curve matching of the data from an observation well with the Hantush type curves. The scales of both graphs are equal and their axes are kept parallel. The value of\u00a0<em>\u03b2<\/em> is obtained from the corresponding curve and the match point values are obtained at the location of the orange triangle (modified from Reed, 1980; Kruseman &amp; de Ridder, 2000).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 75%;\" \/> <col style=\"width: 15%;\" \/><\/colgroup>\r\n<tbody>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]T = \\frac{Q}{4\\pi s}\\ H(u,\\beta)[\/latex]<\/th>\r\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X6e72c096754d5b8d21d281aebd8e380fb388d01\" class=\"anchor\">(60)<\/span><\/th>\r\n<\/tr>\r\n<tr>\r\n<td><\/td>\r\n<td style=\"text-align: center;\">[latex]S = \\frac{u(4Tt)}{r^{2}}[\/latex]<\/td>\r\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"X60f5ece6500904699333d9dcddeea65ef1620a1\" class=\"anchor\">(61)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nwhere:\r\n\r\n<\/div>\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 14%;\" \/> <col style=\"width: 3%;\" \/> <col style=\"width: 67%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th style=\"text-align: right;\"><em>T<\/em><\/th>\r\n<th style=\"text-align: center;\">=<\/th>\r\n<th>transmissivity (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>s<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>drawdown (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>Q<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>constant pumping rate (L<sup>3<\/sup>T<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\">H(<em>u<\/em>,<em>\u03b2<\/em>)<\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>Hantush Well Function (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>u<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td><em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>) (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>radial distance to the observation well (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>t<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>time since pumping started (T)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>storativity of the confined aquifer (dimensionless)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nOnce <em>T<\/em> and <em>S<\/em> are computed for the pumped confined aquifer, estimates of the vertical hydraulic conductivity and storage properties of the confining unit can be derived as described by Schwartz and Zhang (2003). For conditions shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a, when <em>K\u201db\u201dS\u201d<\/em> = 0 then <em>K\u2019S\u2019<\/em> is calculated using Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a> which is a rearrangement of Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a>. When <em>K\u201dS\u201d<\/em>=<em>K\u2019S\u2019<\/em> (<a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a) then <em>K\u2019S\u2019<\/em> is calculated using Equation\u00a0<a href=\"#X052ae1b94024351dd201771b56fed592938d5c4\">(63)<\/a> which is a rearrangement of Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> .\r\n\r\n<\/div>\r\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 10%;\" \/> <col style=\"width: 75%;\" \/> <col style=\"width: 15%;\" \/><\/colgroup>\r\n<tbody>\r\n<tr>\r\n<th><\/th>\r\n<th style=\"text-align: center;\">[latex]K'S' = \\ \\frac{16\\ b'TS\\beta^{2}}{r^{2}}[\/latex]<\/th>\r\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X8ed3b31b003f54e3dc5adaaaa2484dc41261504\" class=\"anchor\">(62)<\/span><\/th>\r\n<\/tr>\r\n<tr>\r\n<td><\/td>\r\n<td style=\"text-align: center;\">[latex]K'S' = \\ \\frac{16\\ TS\\beta^{2}}{r^{2}}\\ \\left( \\frac{b'b\"}{\\text{b' + b\"} + 2(b'b\")^{0.5}} \\right)[\/latex]<\/td>\r\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"X052ae1b94024351dd201771b56fed592938d5c4\" class=\"anchor\">(63)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nwhere:\r\n<table style=\"width: 85%;\" data-custom-style=\"Table\"><colgroup> <col style=\"width: 14%;\" \/> <col style=\"width: 3%;\" \/> <col style=\"width: 67%;\" \/> <\/colgroup>\r\n<thead>\r\n<tr>\r\n<th style=\"text-align: right;\"><em>\u03b2<\/em><\/th>\r\n<th style=\"text-align: center;\">=<\/th>\r\n<th>defined by Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> for Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a>, and Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> for Equation\u00a0<a href=\"#X052ae1b94024351dd201771b56fed592938d5c4\">(63)<\/a><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>radial distance to the pumping well (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S\u2019<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>storativity of the first confining unit (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S\u201d<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>storativity of the second confining unit (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>storativity of the confined unit being pumped (dimensionless)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>transmissivity of the confined unit being pumped (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>K\u2019<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>vertical hydraulic conductivity of the first confining unit (LT<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>b\u2019<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>thickness of the first confining unit (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>K\u201d<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>vertical hydraulic conductivity of the second confining unit (LT<sup>\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"text-align: right;\"><em>b\u201d<\/em><\/td>\r\n<td style=\"text-align: center;\">=<\/td>\r\n<td>thickness of the second confining unit (L)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\n<strong>Example<\/strong>\r\n\r\nThe data presented in <a href=\"#X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\">Figure 54<\/a> represent drawdown in a fully\u2011penetrating observation well responding to a fully\u2011penetrating well in a confined aquifer that is pumped continuously for 1 day at 1,000 m<sup>3<\/sup>\/d. The production well and observation well are completed in a 20 m thick sandstone unit. The observation well is 100 m from the pumping well. The confining unit overlying the leaky confined unit is a sandy shale that is 10 m thick. Based on this information and the curve matching results, the <em>T<\/em> and <em>S<\/em> of the confined sandstone unit can be computed as follows using Equations <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a> and <a href=\"#X60f5ece6500904699333d9dcddeea65ef1620a1\">(61)<\/a>.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]T = \\frac{1000\\ \\frac{m^{3}}{d}}{4\\ (3.14)\\ 35\\ m}\\ (1) = 2.3\\frac{m^{2}}{d}[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]S = \\frac{1\\ (4)\\ 2.3\\frac{m^{2}}{d}0.000125\\ }{{(100\\ m)}^{2}} = \\ 0.0000001\\ or\\ 1x10^{- 7}[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n\r\nA rearrangement of Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> is used to solve for hydraulic conductivity of the confining bed using either a known or estimated value for storativity of the confining bed. Methods to estimate the confined storativity are provided in <span id=\"TextLinkToBox2\" class=\"anchor\"><a href=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/box-2-estimating-storativity-and-specific-storage-ss\/\">Box 2<\/a><\/span>. When the confining unit is fractured rock, Lohman (1972) suggests confining unit storativity can be estimated by multiplying the thickness of the unit in meters by 0.0000003\/m. In this example the confining bed storativity is estimated as 10 m times 0.0000003\/m and equals 0.000003. <em>\u03b2<\/em> from the curve match is 10. Then Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a> is rearranged to solve for K\u2019.\r\n\r\n<\/div>\r\n<div data-custom-style=\"Body Text\">\r\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\r\n\r\n[latex]K' = \\ \\frac{16\\ b'TS\\beta^{2}}{S'r^{2}} = \\ \\frac{16\\ (10\\ m)\\ 2.3\\frac{m^{2}}{d}(0.0000001)\\ ({10)}^{2}}{(0.000003)\\ {(100\\ m)}^{2}} = 0.12\\frac{m}{d}[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>","rendered":"<div id=\"hantush-equation-for-a-leaky-confined-system-with-water-released-from-confining-bed-storage\" class=\"level1\">\n<div data-custom-style=\"First Paragraph\">\n<div data-custom-style=\"First Paragraph\">\n<div data-custom-style=\"First Paragraph\">\n<p>When pumping a confined aquifer with a permeable confining bed contributing water to the pumped aquifer, the reduction in head in the aquifer creates a hydraulic gradient in the confining bed toward the aquifer. As heads in the confining bed are lowered, water is released from storage and flows to the confined aquifer. If the confining bed is relativity thin, the gradient change reaches the top of the bed relatively quickly and then stabilizes such that release of water from storage within the confining bed ceases. If there is an aquifer on the other side of the confining bed, it provides a source of water that flows through the confining bed to the pumped aquifer. Hantush (1960) developed a set of analytical solutions that account for leakage of water from aquitard storage into a confined aquifer as it is pumped. His work also addresses leakage into a confined system from vertical flow though confining beds as described by Hantush\u2011Jacob and illustrated in <a href=\"\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/formulation-of-equations-to-address-leaky-confined-conditions\/#X57b0920e91bd230bfa3e910c632b20097e16b1b\">Figure 41<\/a>. Conceptual models of such leakage are illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>.<\/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: 5.4956in; height: 4.37662in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file53-1.jpg\" alt=\"\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><span id=\"X6b6c2ec720b85923f014beffc01d41c6be08513\" class=\"anchor\"><strong>Figure 49 \u2011<\/strong> Conceptual models that can be represented by the Hantush equation (1960).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote><p>a) Drawdown in a confined aquifer that receives water from storage in one or two confining beds when the beds overlying and underlying the confining beds are not permeable.<\/p><\/blockquote>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote><p>b) Drawdown in a confined aquifer that receives water from storage in overlying and\/or underlying confining beds, as well as leakage of water through the confining bed from a highly transmissive unit such as an unconfined water table aquifer.<\/p><\/blockquote>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote><p>c) Drawdown in a confined aquifer receiving water via direct flow through a confining bed that has negligible storage from a highly transmissive aquifer above and\/or below the confining bed.<\/p><\/blockquote>\n<\/div>\n<div data-custom-style=\"First Paragraph\">\n<p>The general analytical solution for leaky conditions given by Hantush (1960) is presented in Equation\u00a0<a href=\"#X15cd29ce3b49460451a071492929753350a452b\">(55)<\/a>. The different conditions shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a> can be represented using variations of this equation. Similar to the Theis well function for drawdown when pumping a fully confined water bearing unit, the Hantush integral includes an infinite series and can be represented using an approximation as shown in Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a>. The Hantush well function uses the parameter <em>b<\/em> to account for the relative properties of the aquifer and its overlying and underlying layers. Drawdown in a confined unit with one or two adjacent confining unit(s) that release(s) water from storage that is overlain or underlain by an impermeable unit is calculated as shown in Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a>.<\/p>\n<\/div>\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 75%;\" \/>\n<col style=\"width: 15%;\" \/><\/colgroup>\n<tbody>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ \\int_{u}^{\\infty}{\\ \\frac{1}{y}\\ exp\\ \\left\\{ - y \\right\\}\\ erfc\\left( \\frac{\\frac{\\beta}{\\sqrt{u}}}{\\sqrt{y(y - u}\\ )} \\right)}dy[\/latex]<\/th>\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X15cd29ce3b49460451a071492929753350a452b\" class=\"anchor\">(55)<\/span><\/th>\n<\/tr>\n<tr>\n<td><\/td>\n<td style=\"text-align: center;\">[latex]H(u,\\beta) = W(u) - \\ \\frac{4\\beta}{(\\pi{u)}^{0.5}}\\ \\left\\lbrack 0.2577 + 0.6631\\exp\\left( \\frac{- u}{2}\\ \\right) \\right\\rbrack[\/latex]<\/td>\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\" class=\"anchor\">(56)<\/span><\/td>\n<\/tr>\n<tr>\n<td><\/td>\n<td style=\"text-align: center;\">[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ H(u,\\beta)[\/latex]<\/td>\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\" class=\"anchor\">(57)<\/span><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>where:<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 14%;\" \/>\n<col style=\"width: 3%;\" \/>\n<col style=\"width: 67%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th style=\"text-align: right;\"><em>s<\/em>(<em>r<\/em>,<em>t<\/em>)<\/th>\n<th style=\"text-align: center;\">=<\/th>\n<th>drawdown at a radial distance <em>r<\/em> and time <em>t<\/em> (L)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>radial distance to the observation well (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>t<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>time (T)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>Q<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>constant well discharge rate (L<sup>3<\/sup>T<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>transmissivity of the confined pumped unit (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>y<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>variable of integration<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\">erfc<\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>complementary error function (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>\u03b2<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>represents the relative transmissibly and storage of the pumped unit and the contributing overlying and\/or underlying units as shown in Equations\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> and <a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>u<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td><em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\">H(<em>u<\/em>,<em>\u03b2<\/em>)<\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>Hantush well function integral (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\">W(<em>u<\/em>)<\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>Theis well function (dimensionless)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>When storage comes from only one confining bed, <em>b<\/em> is described by Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a>. This represents the case illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a with one aquitard\u2014which can be above or below the pumped unit. In this case the other aquitard shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a is impermeable. For the setting with two confining beds, <em>b<\/em> includes terms for both aquitards and is described by Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a>.<\/p>\n<\/div>\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 75%;\" \/>\n<col style=\"width: 15%;\" \/><\/colgroup>\n<tbody>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]\\beta = \\frac{r}{4\\ }\\sqrt{\\frac{K'S'}{b'TS}\\ }[\/latex]<\/th>\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X35b59d96334da087626420f015e0710c839455d\" class=\"anchor\">(58)<\/span><\/th>\n<\/tr>\n<\/tbody>\n<\/table>\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 75%;\" \/>\n<col style=\"width: 15%;\" \/><\/colgroup>\n<tbody>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]\\beta = \\frac{r}{4\\ }\\left( \\sqrt{\\frac{K'S'}{b'TS}\\ }\\ + \\sqrt{\\frac{K\"S\"}{b^{\"}TS}} \\right)[\/latex]<\/th>\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X7ef8f34ca870785ea61794264c127ffbbce8681\" class=\"anchor\">(59)<\/span><\/th>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>where:<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 14%;\" \/>\n<col style=\"width: 3%;\" \/>\n<col style=\"width: 67%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th style=\"text-align: right;\"><em>r<\/em><\/th>\n<th style=\"text-align: left;\">=<\/th>\n<th>radial distance to the pumping well (L)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"text-align: right;\"><em>K\u2019<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>vertical hydraulic conductivity of the first confining unit (LT<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>K\u201d<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>vertical hydraulic conductivity of the second confining unit (LT<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S\u2019<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>storativity of the first confining unit (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S\u201d<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>storativity of the second confining unit (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>b\u2019<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>thickness of the first confining unit (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>b\u201d<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>thickness of the second confining unit (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>transmissivity of the confined unit being pumped (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\n<td style=\"text-align: left;\">=<\/td>\n<td>storativity of the confined unit being pumped (dimensionless)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>Drawdown in the confined aquifer can be calculated with Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> by using Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a> to calculate the Hantush well function (H(<em>u<\/em>,<em>\u03b2<\/em>)) or reading\/interpolating the value from a table (<a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a>) or type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>).<\/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: 2.20962in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file54-1.jpg\" alt=\"\" \/><\/p>\n<div class=\"bc-figcaption figcaption\">\n<div data-custom-style=\"Image Caption\">\n<div class=\"figcaption-text\">\n<p><span id=\"X32295cc7f397b2f30defdbe73d897d6eb898586\" class=\"anchor\"><strong>Figure 50 \u2011<\/strong> Tables of the Hantush well function for variables <em>u<\/em> and <em>\u03b2<\/em> (from Fetter (2001)). Type curves of H(<em>u<\/em>,<em>\u03b2<\/em>) versus 1\/<em>u<\/em> are illustrated in <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a><\/span>.<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 5.8in; height: 3.52in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file55-1.jpg\" alt=\"\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><span id=\"X21bd8ced2a72427e089642240b99c2f3bd335eb\" class=\"anchor\"><strong>Figure 51 \u2011<\/strong> Type curves of H(<em>u<\/em>,<em><u>\u03b2<\/u><\/em>) versus 1\/<em>u<\/em> plotted on log\u2011log axes (modified from Reed (1980); Kruseman &amp; de RIdder (2000); Kasenow (2001)).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<div data-custom-style=\"First Paragraph\">\n<p>The table of values (<a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a>) and plotted type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>) represent a subset of H(<em>u<\/em>,<em>\u03b2<\/em>) values that can be generated by substituting combinations of <em>u<\/em> and <em>\u03b2<\/em> into Equation\u00a0<a href=\"#Xed9b7003846df5cf9bdf4022b5769bb80d71c8b\">(56)<\/a> or by using mathematical software packages to solve the integral portion of Equation\u00a0<a href=\"#X15cd29ce3b49460451a071492929753350a452b\">(55)<\/a>. The required values are commonly interpolated from the table or curves presented in <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> and <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a> respectively.<\/p>\n<\/div>\n<div id=\"using-the-hantush-equation-to-predict-drawdown-in-leaky-confined-units-with-water-released-from-confining-bed-storage\" class=\"level3\">\n<h3>Using the Hantush Equation to Predict Drawdown in Leaky Confined Units with Water Released from Confining Bed Storage<\/h3>\n<div data-custom-style=\"First Paragraph\">\n<p>When the hydraulic properties of the sequence of geologic units as illustrated in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a are known, the Hantush equation can be used to approximate drawdown at an observation point at any time <em>t<\/em> after the pumping commenced.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p><strong>Example<\/strong><\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>A fully penetrating well pumps at 500 m<sup>3<\/sup>\/d for 0.2\u00a0d from a 30 m thick leaky confined sand rich aquifer that is overlain by a silt rich aquitard that is 20 m thick. The confining units are overlain and underlain by impermeable units (<a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a). The observation well is located 53 m from the pumping well. The transmissivity of the confined aquifer is 300 m<sup>2<\/sup>\/d and the storativity is 0.00003. The vertical hydraulic conductivity of the confining bed is 0.3 m\/d and the storativity of the confining bed is 0.00001. Calculate the drawdown in the observation well after 0.2 d of pumping.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>Using the definition of <em>u<\/em> (<em>u<\/em>=<em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)) along with Equations <a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> and <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a>, and the information listed above, H(<em>u<\/em>,<em>b<\/em>) is obtained by first computing <em>u<\/em> and <em>\uf062<\/em>, and then using <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> to obtain the Hantush well function.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]u = \\frac{r^{2}S}{4Tt} = \\frac{{(53\\ m)}^{2}\\ 0.00003}{4\\ \\left( 300\\frac{m^{2}}{d} \\right)0.2\\ d} = 0.0004[\/latex]<\/p>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]\\beta = \\frac{r}{4\\ }\\sqrt{\\frac{K'S'}{b'TS}\\ } = \\frac{53\\ m}{4\\ }\\sqrt{\\frac{0.3\\frac{m}{d}(0.00001)}{20\\ m\\ \\left( 300\\frac{m^{2}}{d} \\right)\\ 0.00003}} = 0.05[\/latex]<\/p>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>With <em>u<\/em> = 0.0004 and <em>b<\/em> = 0.05, H(<em>u<\/em>,<em>b<\/em>) is obtained by interpolation from <a href=\"#X32295cc7f397b2f30defdbe73d897d6eb898586\">Figure 50<\/a> or <a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>. With <em>u<\/em> = 0.0004, <em>b<\/em> = 0.05 then H(<em>u<\/em>,<em>b<\/em>) is about 6.52.<\/p>\n<\/div>\n<table style=\"width: 37%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 8%;\" \/>\n<col style=\"width: 18%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th><em>u<\/em><\/th>\n<th><em>\u03b2<\/em><\/th>\n<th>H(<em>u<\/em>,<em>\u03b2<\/em>)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td>0.0001<\/td>\n<td>0.05<\/td>\n<td>7.1771<\/td>\n<\/tr>\n<tr>\n<td>0.0004<\/td>\n<td>0.05<\/td>\n<td>6.52 (interpolated)<\/td>\n<\/tr>\n<tr>\n<td>0.0005<\/td>\n<td>0.05<\/td>\n<td>6.3523<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Captioned Figure\">\n<p><img decoding=\"async\" style=\"display: block; margin-left: auto; margin-right: auto; ;width: 5.8in; height: 4.13in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file56-1.jpg\" alt=\"\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><strong>Figure 52 \u2011<\/strong> Hantush type curves plotted on log\u2011log scales. The dashed vertical blue line is the computed value of 1\/<em>u<\/em> (2,500) for the example problem. The red line is the interpolated value of <em>\u03b2<\/em>, about 0.05, and the black dot dashed line is the interpolated value of H(<em>u<\/em>,<em>\u03b2<\/em>), about 5.3. This value is lower than derived from the tables. The close plotting of the type curves made interpretation difficult for the designated <em>\u03b2\u00a0<\/em>value (modified from Kruseman &amp; de RIdder, 2000; Kasenow, 2001).<\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<div data-custom-style=\"First Paragraph\">\n<p>Then Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> provides the value of drawdown.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]s(r,t) = \\frac{500\\frac{m^{3}}{d}}{4\\ (3.14)\\ 300\\frac{m^{2}}{d}}\\ 6.52 = 0.87\\ m[\/latex]<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div id=\"hantush-curve-matching-method-to-compute-t-and-s-from-a-pumping-test-in-leaky-confined-unit-with-aquitard-storage\" class=\"level3\">\n<h3>Hantush Curve Matching Method to Compute T and S from a Pumping Test in Leaky Confined Unit with Aquitard Storage<\/h3>\n<div data-custom-style=\"First Paragraph\">\n<p>The Hantush equation can be used to determine <em>T<\/em> and <em>S<\/em> from time\u2011drawdown data collected in an observation well when a leaky confined aquifer is pumped and water is released from aquitard storage. The hydraulic conductivity and storage properties of adjacent confining units can also be derived. Standard curve\u2011matching techniques described in Section 8.3.1 can be used.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>After data from a pumping test has been collected, <em>T<\/em> and <em>S<\/em> for a leaky confined aquifer with the release of water from confining bed storage is estimated by first plotting drawdown versus time on a log\u2011log graph at the same scale as the type curves as shown in (<a href=\"#X87f36bd97f7d73b20af2133215aeb032eb87b68\">Figure 53<\/a>).<\/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: 4in; height: 2.44in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file57-1.jpg\" alt=\"\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><span id=\"X87f36bd97f7d73b20af2133215aeb032eb87b68\" class=\"anchor\"><strong>Figure 53 \u2011<\/strong> Log\u2011log plot of the response of an observation well to the pumping of a well in a leaky confined aquifer with water being released from confining bed storage.<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<div data-custom-style=\"First Paragraph\">\n<p>The data graph (<a href=\"#X87f36bd97f7d73b20af2133215aeb032eb87b68\">Figure 53<\/a>) is overlain on the type curves (<a href=\"#X21bd8ced2a72427e089642240b99c2f3bd335eb\">Figure 51<\/a>). Keeping the axes parallel, the field data graph and type curve graph are overlain and shifted until the curves match (<a href=\"#X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\">Figure 54<\/a>). A match point is selected from the overlapping portion of the two plots and values of H(<em>u<\/em>,<em>\u03b2<\/em>), 1\/<em>u<\/em>, <em>s<\/em>, <em>t<\/em> and <em>\u03b2<\/em> are obtained. Estimates of <em>T<\/em> and <em>S<\/em> for the leaky confined unit as well as its hydraulic conductivity, and, at early times, storativity values of the confining beds can be estimated using Equations <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a> and <a href=\"#X60f5ece6500904699333d9dcddeea65ef1620a1\">(61)<\/a> which are obtained from rearranging Equation\u00a0<a href=\"#Xcd361c4ebba676e40fed3f3f0c755daeb579a9f\">(57)<\/a> and the definition of <em>u<\/em> (<em>u<\/em>=<em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>)), respectively. Early time refers to the period after pumping begins when the confining unit is releasing water from storage. Early time is less than <em>b\u2019S\u2019\/10K\u2019<\/em>.<\/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: 5.28182in; height: 4.07756in;\" src=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-content\/uploads\/sites\/40\/2026\/08\/file58-1.jpg\" alt=\"\" \/><\/p>\n<\/div>\n<div data-custom-style=\"Block Text\">\n<blockquote>\n<div class=\"figcaption-text\">\n<p><span id=\"X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\" class=\"anchor\"><strong>Figure 54 \u2011<\/strong> Curve matching of the data from an observation well with the Hantush type curves. The scales of both graphs are equal and their axes are kept parallel. The value of\u00a0<em>\u03b2<\/em> is obtained from the corresponding curve and the match point values are obtained at the location of the orange triangle (modified from Reed, 1980; Kruseman &amp; de Ridder, 2000).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 75%;\" \/>\n<col style=\"width: 15%;\" \/><\/colgroup>\n<tbody>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]T = \\frac{Q}{4\\pi s}\\ H(u,\\beta)[\/latex]<\/th>\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X6e72c096754d5b8d21d281aebd8e380fb388d01\" class=\"anchor\">(60)<\/span><\/th>\n<\/tr>\n<tr>\n<td><\/td>\n<td style=\"text-align: center;\">[latex]S = \\frac{u(4Tt)}{r^{2}}[\/latex]<\/td>\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"X60f5ece6500904699333d9dcddeea65ef1620a1\" class=\"anchor\">(61)<\/span><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>where:<\/p>\n<\/div>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 14%;\" \/>\n<col style=\"width: 3%;\" \/>\n<col style=\"width: 67%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th style=\"text-align: right;\"><em>T<\/em><\/th>\n<th style=\"text-align: center;\">=<\/th>\n<th>transmissivity (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"text-align: right;\"><em>s<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>drawdown (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>Q<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>constant pumping rate (L<sup>3<\/sup>T<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\">H(<em>u<\/em>,<em>\u03b2<\/em>)<\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>Hantush Well Function (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>u<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td><em>r<\/em><sup>2<\/sup><em>S<\/em>\/(4<em>Tt<\/em>) (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>radial distance to the observation well (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>t<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>time since pumping started (T)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>storativity of the confined aquifer (dimensionless)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>Once <em>T<\/em> and <em>S<\/em> are computed for the pumped confined aquifer, estimates of the vertical hydraulic conductivity and storage properties of the confining unit can be derived as described by Schwartz and Zhang (2003). For conditions shown in <a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a, when <em>K\u201db\u201dS\u201d<\/em> = 0 then <em>K\u2019S\u2019<\/em> is calculated using Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a> which is a rearrangement of Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a>. When <em>K\u201dS\u201d<\/em>=<em>K\u2019S\u2019<\/em> (<a href=\"#X6b6c2ec720b85923f014beffc01d41c6be08513\">Figure 49<\/a>a) then <em>K\u2019S\u2019<\/em> is calculated using Equation\u00a0<a href=\"#X052ae1b94024351dd201771b56fed592938d5c4\">(63)<\/a> which is a rearrangement of Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> .<\/p>\n<\/div>\n<table class=\"gw-eq-table\" style=\"width: 100%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 10%;\" \/>\n<col style=\"width: 75%;\" \/>\n<col style=\"width: 15%;\" \/><\/colgroup>\n<tbody>\n<tr>\n<th><\/th>\n<th style=\"text-align: center;\">[latex]K'S' = \\ \\frac{16\\ b'TS\\beta^{2}}{r^{2}}[\/latex]<\/th>\n<th style=\"text-align: right; white-space: nowrap;\"><span id=\"X8ed3b31b003f54e3dc5adaaaa2484dc41261504\" class=\"anchor\">(62)<\/span><\/th>\n<\/tr>\n<tr>\n<td><\/td>\n<td style=\"text-align: center;\">[latex]K'S' = \\ \\frac{16\\ TS\\beta^{2}}{r^{2}}\\ \\left( \\frac{b'b\"}{\\text{b' + b\"} + 2(b'b\")^{0.5}} \\right)[\/latex]<\/td>\n<td style=\"text-align: right; white-space: nowrap;\"><span id=\"X052ae1b94024351dd201771b56fed592938d5c4\" class=\"anchor\">(63)<\/span><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<div data-custom-style=\"Body Text\">\n<p>where:<\/p>\n<table style=\"width: 85%;\" data-custom-style=\"Table\">\n<colgroup>\n<col style=\"width: 14%;\" \/>\n<col style=\"width: 3%;\" \/>\n<col style=\"width: 67%;\" \/> <\/colgroup>\n<thead>\n<tr>\n<th style=\"text-align: right;\"><em>\u03b2<\/em><\/th>\n<th style=\"text-align: center;\">=<\/th>\n<th>defined by Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> for Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a>, and Equation\u00a0<a href=\"#X7ef8f34ca870785ea61794264c127ffbbce8681\">(59)<\/a> for Equation\u00a0<a href=\"#X052ae1b94024351dd201771b56fed592938d5c4\">(63)<\/a><\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td style=\"text-align: right;\"><em>r<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>radial distance to the pumping well (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S\u2019<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>storativity of the first confining unit (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S\u201d<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>storativity of the second confining unit (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>S<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>storativity of the confined unit being pumped (dimensionless)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>T<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>transmissivity of the confined unit being pumped (L<sup>2<\/sup>T<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>K\u2019<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>vertical hydraulic conductivity of the first confining unit (LT<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>b\u2019<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>thickness of the first confining unit (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>K\u201d<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>vertical hydraulic conductivity of the second confining unit (LT<sup>\u20111<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: right;\"><em>b\u201d<\/em><\/td>\n<td style=\"text-align: center;\">=<\/td>\n<td>thickness of the second confining unit (L)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p><strong>Example<\/strong><\/p>\n<p>The data presented in <a href=\"#X81e66f51c403b0a47c65b949c145b3a0b3a3ff3\">Figure 54<\/a> represent drawdown in a fully\u2011penetrating observation well responding to a fully\u2011penetrating well in a confined aquifer that is pumped continuously for 1 day at 1,000 m<sup>3<\/sup>\/d. The production well and observation well are completed in a 20 m thick sandstone unit. The observation well is 100 m from the pumping well. The confining unit overlying the leaky confined unit is a sandy shale that is 10 m thick. Based on this information and the curve matching results, the <em>T<\/em> and <em>S<\/em> of the confined sandstone unit can be computed as follows using Equations <a href=\"#X6e72c096754d5b8d21d281aebd8e380fb388d01\">(60)<\/a> and <a href=\"#X60f5ece6500904699333d9dcddeea65ef1620a1\">(61)<\/a>.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]T = \\frac{1000\\ \\frac{m^{3}}{d}}{4\\ (3.14)\\ 35\\ m}\\ (1) = 2.3\\frac{m^{2}}{d}[\/latex]<\/p>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]S = \\frac{1\\ (4)\\ 2.3\\frac{m^{2}}{d}0.000125\\ }{{(100\\ m)}^{2}} = \\ 0.0000001\\ or\\ 1x10^{- 7}[\/latex]<\/p>\n<\/div>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<p>A rearrangement of Equation\u00a0<a href=\"#X35b59d96334da087626420f015e0710c839455d\">(58)<\/a> is used to solve for hydraulic conductivity of the confining bed using either a known or estimated value for storativity of the confining bed. Methods to estimate the confined storativity are provided in <span id=\"TextLinkToBox2\" class=\"anchor\"><a href=\"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/box-2-estimating-storativity-and-specific-storage-ss\/\">Box 2<\/a><\/span>. When the confining unit is fractured rock, Lohman (1972) suggests confining unit storativity can be estimated by multiplying the thickness of the unit in meters by 0.0000003\/m. In this example the confining bed storativity is estimated as 10 m times 0.0000003\/m and equals 0.000003. <em>\u03b2<\/em> from the curve match is 10. Then Equation\u00a0<a href=\"#X8ed3b31b003f54e3dc5adaaaa2484dc41261504\">(62)<\/a> is rearranged to solve for K\u2019.<\/p>\n<\/div>\n<div data-custom-style=\"Body Text\">\n<div class=\"display-math\" data-custom-style=\"Display Equation\">\n<p>[latex]K' = \\ \\frac{16\\ b'TS\\beta^{2}}{S'r^{2}} = \\ \\frac{16\\ (10\\ m)\\ 2.3\\frac{m^{2}}{d}(0.0000001)\\ ({10)}^{2}}{(0.000003)\\ {(100\\ m)}^{2}} = 0.12\\frac{m}{d}[\/latex]<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"author":8,"menu_order":37,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-1678","chapter","type-chapter","status-publish","hentry"],"part":611,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1678","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/users\/8"}],"version-history":[{"count":8,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1678\/revisions"}],"predecessor-version":[{"id":3759,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1678\/revisions\/3759"}],"part":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/parts\/611"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1678\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/media?parent=1678"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapter-type?post=1678"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/contributor?post=1678"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/license?post=1678"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}