{"id":238,"date":"2022-01-13T23:17:25","date_gmt":"2022-01-13T23:17:25","guid":{"rendered":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/delayed-compaction-of-aquitards-confining-beds\/"},"modified":"2022-03-30T20:19:41","modified_gmt":"2022-03-30T20:19:41","slug":"delayed-compaction-of-aquitards-confining-beds","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/delayed-compaction-of-aquitards-confining-beds\/","title":{"raw":"2.5  Delayed Compaction of Aquitards (Confining Beds)","rendered":"2.5  Delayed Compaction of Aquitards (Confining Beds)"},"content":{"raw":"<div class=\"delayed-compaction-of-aquitards-(confining-beds)\">\r\n<p class=\"import-Normal\">An aquitard (or confining bed) is a clayey\u2011silty low permeability formation that does not provide an appreciable quantity of groundwater to pumping wells; however, it can transmit appreciable water to adjacent aquifers. While flow in an aquifer is predominantly two\u2011dimensional (2\u2011D) and horizontal, particularly if wellbores are fully penetrating, flow in the aquitards separating the aquifers is mostly 1\u2011D and vertical. In a complex aquifer system (for example, Figure\u00a016) the role played by the intervening aquitards is important as they can represent a significant source of water to the aquifers and can contribute greatly to land subsidence as clay\/silt compressibility <em>c<\/em><sub><em>b<\/em><\/sub> is usually much larger than that of the sand\/gravel.<\/p>\r\n<p class=\"import-Normal\"><img class=\"alignnone\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/uploads\/sites\/20\/2022\/01\/image16.png\" alt=\"a)\u00a0Digital elevation model of the Emilia\u2011Romagna plain, Italy, and b) vertical cross section showing the complex multi\u2011aquifer system used to supply freshwater to the coastland.\" width=\"1461\" height=\"628\" \/><\/p>\r\n<p class=\"import-Normal figcaption-text\"><strong>Figure\u00a0<\/strong><strong>16<\/strong><strong>\u00a0<\/strong><strong>\u2011\u00a0<\/strong>a)\u00a0Digital elevation model of the Emilia\u2011Romagna plain, Italy, and b) vertical cross section along the A\u2011A' alignment shown in (a)\u00a0of the complex multi\u2011aquifer system used to supply freshwater to the coastland (modified after Teatini et al., 2006).<\/p>\r\n<p class=\"import-Normal\">Normally aquitard compaction is larger and delayed in time relative to aquifer compaction. The law that governs pore-water decline in the aquitard as a function of time and the factors controlling compaction are explained in this section. Darcy\u2019s law describing the velocity of groundwater flow in an aquitard can be written as shown in Equation\u00a014.<a id=\"eq_14\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle v_{z}=-K\\frac{\\partial h}{\\partial z}[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(14)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>v<\/em><sub><em>z<\/em><\/sub><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">apparent seepage velocity (LT<sup>-1<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>K<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">hydraulic conductivity (LT<sup>-1<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>h<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">hydraulic head = <em>z<\/em> + <em>p<\/em>\/<em>\u03b3<\/em> (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>z<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">vertical coordinate positive downward (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u2202h<\/em>\/<em>\u2202z<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">vertical hydraulic gradient (LL<sup>-1<\/sup>)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">The hydraulic conductivity is a function of the physical properties of fluid and soil as shown in Equation\u00a015.<a id=\"eq_15\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle K=k^{*}\\frac{\\gamma }{\\mu }[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(15)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>k*<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">intrinsic permeability (L<sup>2<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03b3<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">specific weight of water (ML<sup>-2<\/sup>T<sup>-2<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03bc<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">dynamic viscosity of water (ML<sup>-1<\/sup>T<sup>-1<\/sup>)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">Intrinsic permeability is dependent exclusively on the properties of the medium:<\/p>\r\n<p style=\"text-align: center;\"><em>k*<\/em> = <em>CD<\/em><sup>2<\/sup><\/p>\r\n<p class=\"import-Normal\">where:<\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>D<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">a representative length of the porous medium (for example, the average grain size) (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>C<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">appropriate parameter related to the soil type (dimensionless)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">Other more complex relationships (depending on porosity, mean pore diameter, and specific surface area) have been developed for intrinsic permeability of reactive clays, especially if salt is dissolved into the pore water (for example, Raffensperger and Ferrell Jr., 1991).<\/p>\r\n<p class=\"import-Normal\">Assume the initial conditions are in equilibrium, and all the hydrologic and geomechanical quantities presented here are incremental with respect to the initial conditions. Let\u2019s balance the weight of water in an elementary soil sample of initial length \u0394<em>z<\/em> and unitary cross\u2011sectional area (shown as 1 in the expressions below) between time <em class=\"import-Cambria\">t<\/em> and <em class=\"import-Cambria\">t<\/em>\u00a0+\u00a0\u0394<em class=\"import-Cambria\">t<\/em>:<\/p>\r\n\r\n<ul>\r\n \t<li class=\"import-Normal\">Inflow: (<em>\u03b3<\/em><em>v<\/em><sub>z<\/sub>)\u00a0(1)\u00a0(\u0394<em class=\"import-Cambria\">t<\/em>)<\/li>\r\n \t<li class=\"import-Normal\">Outflow: <em>\u03b3<\/em>(<em>v<\/em><sub>z<\/sub>\u00a0+\u00a0<em class=\"import-Cambria\">\u2202<\/em><em>v<\/em><sub>z<\/sub>\/<em class=\"import-Cambria\">\u2202<\/em><em>z<\/em><em>\u00a0<\/em>\u0394<em>z<\/em>)\u00a0(1)\u00a0(\u0394<em>t<\/em>)<\/li>\r\n \t<li class=\"import-Normal\">Weight of water expelled by the porous space contraction and the expansion of the water expressed by Equation 16 (we assume incompressible solid grains the total medium volume change coincides with the porous volume change):<a id=\"eq_16\"><\/a><\/li>\r\n<\/ul>\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">\u2212[(<em>\u03b3<\/em>\u2206(<em>\u03d5<\/em>\u2206<em>z<\/em>) 1 \u2206<em>p<\/em>) + (<em>\u03b3\u03d5\u03b2<\/em> 1 \u2206<em>p<\/em>)]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(16)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03b2<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">volumetric compressibility of water (ML<sup>-1<\/sup>T<sup>-2<\/sup>)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">In Equation\u00a016 the total geostatic stress <em>\u03c3<\/em><sub><em>c<\/em><\/sub> is assumed to be constant, so (from <a href=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/effective-intergranular-stress-and-soil-parameters#eq_3\">Equation\u00a03<\/a>):<\/p>\r\n<p class=\"import-Normal\" style=\"text-align: center;\"><span style=\"display: block; font-size: NaNpt; color: #; ; text-decoration: none;\">\u0394<em>\u03c3<\/em><sub><em>z<\/em><\/sub> = \u2212\u2206<em>p<\/em><\/span><\/p>\r\n<p class=\"import-Normal\">The change in pressure, \u0394<em>p<\/em><em>,<\/em> is negative when <em class=\"import-Cambria\">p<\/em> is reduced, as happens during groundwater pumping. Notice that \u0394(<em>\u03d5<\/em>\u0394<em>z<\/em>) is equal to <em class=\"import-Cambria\">\u2206<\/em>{[<em class=\"import-Cambria\">e<\/em>\/(1\u00a0+\u00a0<em class=\"import-Cambria\">e<\/em>)]<em class=\"import-Cambria\">\u2206z<\/em>} with <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em>\/(1\u00a0+\u00a0<em class=\"import-Cambria\">e<\/em>) constant because this is the solid part (grains) of the elementary volume (1)\u00a0<em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em> . Hence (from <a href=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/effective-intergranular-stress-and-soil-parameters#eq_5\">Equation\u00a05<\/a>) we have:<\/p>\r\n<p style=\"text-align: center;\">[latex]\\displaystyle \\frac{\\Delta e}{1+e}=c_{b}\\Delta p[\/latex]<\/p>\r\n<p class=\"import-Normal\">and therefore, we obtain:<\/p>\r\n\r\n<ul>\r\n \t<li class=\"import-Normal\">Outflow \u2013 Inflow = Weight of water expelled, i.e.<\/li>\r\n<\/ul>\r\n<p style=\"text-align: center;\">[latex]\\displaystyle \\gamma \\frac{\\partial v_{z}}{\\partial z}\\Delta z\\Delta t=-\\gamma (c_{b}+\\phi \\beta )\\Delta z\\Delta p[\/latex]<\/p>\r\n<p class=\"import-Normal\">Cancelling <em>\u03b3<\/em> and <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em> on both sides and remembering that the hydraulic head <em class=\"import-Cambria\">h<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">z<\/em>\u00a0+\u00a0<em class=\"import-Cambria\">p<\/em>\/<em class=\"import-Cambria\">\u03b3<\/em>, we know that <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">p<\/em> = <em class=\"import-Cambria\">\u03b3<\/em><em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">h<\/em>, and using Equation\u00a014 when the increment of time approaches zero <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">t<\/em><em class=\"import-Cambria\">\u00a0<\/em><em class=\"import-Cambria\">\u2192<\/em><em class=\"import-Cambria\">\u00a0<\/em>0 we obtain Equation\u00a017:<a id=\"eq_17\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle \\frac{\\partial }{\\partial z}\\left ( K\\frac{\\partial h}{\\partial z} \\right )=\\gamma (c_{b}+\\phi \\beta )\\frac{\\partial h}{\\partial t}[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(17)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">Solving Equation\u00a017, complemented with the appropriate top and bottom boundary conditions and initial conditions, provides the pressure dissipation within the aquitard, and hence the <em class=\"import-Cambria\">\u0394p<\/em> needed to compute the aquitard compaction versus time. The specific storage coefficient is defined in Equation\u00a018.<a id=\"eq_18\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\"><em>S<\/em><sub><em>s<\/em><\/sub> = <em>\u03b3<\/em>(<em>c<\/em><sub><em>b<\/em><\/sub> + <em>\u03d5\u03b2<\/em>)<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(18)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\"><em>S<\/em><sub><em>s<\/em><\/sub> represents the \u201c<em>specific elastic storage<\/em>\u201d [L<sup>\u2011<\/sup><sup>1<\/sup>] and along with the hydraulic conductivity, <em class=\"import-Cambria\">K<\/em>, defines Terzaghi\u2019s consolidation coefficient <em>c<\/em><sub><em>v<\/em><\/sub> that controls both magnitude and timing of aquitard compaction as shown in Equation\u00a019 given that <em class=\"import-Cambria\">\u03d5\u03b2\u00a0<\/em>&lt;&lt;\u00a0<em>c<\/em><sub><em>b<\/em><\/sub> for typical aquifer confining beds.<a id=\"eq_19\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle c_{v}=\\frac{K}{\\gamma (c_{b}+\\phi \\beta )}\\approx \\frac{K}{\\gamma c_{b}}[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(19)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">The initial conditions correspond to \u0394<em class=\"import-Cambria\">p<\/em>\u00a0=\u00a00 for the entire thickness, <em class=\"import-Cambria\">b<\/em>, of the aquitard while the boundary conditions are given by \u0394<em class=\"import-Cambria\">p<\/em> in the overlying and underlying aquifers. If the pressure drop \u0394<em class=\"import-Cambria\">p<\/em><sub class=\"import-Cambria\">0<\/sub> is the same at top (<em class=\"import-Cambria\">z<\/em>\u00a0=\u00a00) and bottom (<em class=\"import-Cambria\">z<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">b<\/em>), then pressure conditions are symmetrical above and below the middle of the aquitard, hence <span style=\"font-size: NaNpt; color: #; ; text-decoration: none;\">\u2202<em>p<\/em>\/\u2202<em>z<\/em> = 0 <\/span>at <em class=\"import-Cambria\">z<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">b<\/em>\/2. In this case, the solution to Equation\u00a017 can be expressed in terms of a series expansion. Writing the solution in terms of <em class=\"import-Cambria\">p<\/em> produces Equation\u00a020.<a id=\"eq_20\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle p=p_{0}-\\frac{4}{\\pi }\\Delta p_{0}\\left [ \\frac{\\pi }{4}-\\frac{\\textup{sin}(\\pi z\/b)}{\\textup{exp}[(\\pi \/b)^{2}c_{v}t]}-\\frac{1}{3}\\frac{\\textup{sin}(3\\pi z\/b)}{\\textup{exp}[(3\\pi \/b)^{2}c_{v}t]}-... \\right ][\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(20)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>z<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">vertical coordinate positive downward starting from the aquitard top (0 \u2264 <em>z<\/em> \u2264 <em>b<\/em>\/2) (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>t<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">time since the initial change in pressure at the aquitard boundary (T)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">For <em class=\"import-Cambria\">t<\/em> = 0, and setting <em class=\"import-Cambria\">\u03c0z<\/em>\/<em class=\"import-Cambria\">b<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">x<\/em>\u00a0(0\u00a0\u2264\u00a0<em class=\"import-Cambria\">x<\/em>\u00a0\u2264\u00a0<em class=\"import-Cambria\">\u03c0<\/em>\/2) Equation\u00a020 becomes Equation\u00a021.<a id=\"eq_21\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle p(z,0)=p_{0}(z)-\\frac{4}{\\pi }\\Delta p_{0}\\left [ \\frac{\\pi }{4}-\\left ( \\textup{sin}x+\\frac{1}{3}\\textup{sin}3x+\\frac{1}{5}\\textup{sin}5x+...\\right ) \\right ][\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(21)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">The content within the parentheses of Equation\u00a021 is the Fourier series development of the function <em class=\"import-Cambria\">f<\/em><em class=\"import-Cambria\">\u00a0<\/em>(<em class=\"import-Cambria\">x<\/em>):<\/p>\r\n<p style=\"text-align: center;\">[latex]\\displaystyle f(x)=\\begin{cases}\r\n&amp; \\ \\ \\pi \/4\\text{ for } 0\\leq x&lt; \\pi \\\\\r\n&amp; -\\pi \/4\\text{ for } \\pi &lt; x\\leq 2\\pi\r\n\\end{cases}[\/latex]<\/p>\r\n<p class=\"import-Normal\">Keeping in mind the range of variability of <em class=\"import-Cambria\">x\u00a0<\/em>(0\u00a0\u2264\u00a0<em class=\"import-Cambria\">x<\/em><span lang=\"pt-BR\" xml:lang=\"pt-BR\">\u00a0<\/span>\u2264\u00a0<em class=\"import-Cambria\">\u03c0<\/em>\/2) we conclude that Equation\u00a020 accurately represents the initial pore pressure at time <em class=\"import-Cambria\">t<\/em>\u00a0=\u00a00.<\/p>\r\n<p class=\"import-Normal\"><\/p>\r\n\r\n<\/div>","rendered":"<div class=\"delayed-compaction-of-aquitards-(confining-beds)\">\n<p class=\"import-Normal\">An aquitard (or confining bed) is a clayey\u2011silty low permeability formation that does not provide an appreciable quantity of groundwater to pumping wells; however, it can transmit appreciable water to adjacent aquifers. While flow in an aquifer is predominantly two\u2011dimensional (2\u2011D) and horizontal, particularly if wellbores are fully penetrating, flow in the aquitards separating the aquifers is mostly 1\u2011D and vertical. In a complex aquifer system (for example, Figure\u00a016) the role played by the intervening aquitards is important as they can represent a significant source of water to the aquifers and can contribute greatly to land subsidence as clay\/silt compressibility <em>c<\/em><sub><em>b<\/em><\/sub> is usually much larger than that of the sand\/gravel.<\/p>\n<p class=\"import-Normal\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/uploads\/sites\/20\/2022\/01\/image16.png\" alt=\"a)\u00a0Digital elevation model of the Emilia\u2011Romagna plain, Italy, and b) vertical cross section showing the complex multi\u2011aquifer system used to supply freshwater to the coastland.\" width=\"1461\" height=\"628\" \/><\/p>\n<p class=\"import-Normal figcaption-text\"><strong>Figure\u00a0<\/strong><strong>16<\/strong><strong>\u00a0<\/strong><strong>\u2011\u00a0<\/strong>a)\u00a0Digital elevation model of the Emilia\u2011Romagna plain, Italy, and b) vertical cross section along the A\u2011A&#8217; alignment shown in (a)\u00a0of the complex multi\u2011aquifer system used to supply freshwater to the coastland (modified after Teatini et al., 2006).<\/p>\n<p class=\"import-Normal\">Normally aquitard compaction is larger and delayed in time relative to aquifer compaction. The law that governs pore-water decline in the aquitard as a function of time and the factors controlling compaction are explained in this section. Darcy\u2019s law describing the velocity of groundwater flow in an aquitard can be written as shown in Equation\u00a014.<a id=\"eq_14\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-0812bbfbf36b83441485766273806964_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#118;&#95;&#123;&#122;&#125;&#61;&#45;&#75;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#104;&#125;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#122;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"37\" width=\"93\" style=\"vertical-align: -12px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(14)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>v<\/em><sub><em>z<\/em><\/sub><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">apparent seepage velocity (LT<sup>-1<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>K<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">hydraulic conductivity (LT<sup>-1<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>h<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">hydraulic head = <em>z<\/em> + <em>p<\/em>\/<em>\u03b3<\/em> (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>z<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">vertical coordinate positive downward (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u2202h<\/em>\/<em>\u2202z<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">vertical hydraulic gradient (LL<sup>-1<\/sup>)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">The hydraulic conductivity is a function of the physical properties of fluid and soil as shown in Equation\u00a015.<a id=\"eq_15\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-7cda336ad3222d7a5cd4687041cb8166_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#75;&#61;&#107;&#94;&#123;&#42;&#125;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#125;&#123;&#92;&#109;&#117;&#32;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"36\" width=\"70\" style=\"vertical-align: -16px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(15)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>k*<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">intrinsic permeability (L<sup>2<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03b3<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">specific weight of water (ML<sup>-2<\/sup>T<sup>-2<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03bc<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">dynamic viscosity of water (ML<sup>-1<\/sup>T<sup>-1<\/sup>)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">Intrinsic permeability is dependent exclusively on the properties of the medium:<\/p>\n<p style=\"text-align: center;\"><em>k*<\/em> = <em>CD<\/em><sup>2<\/sup><\/p>\n<p class=\"import-Normal\">where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>D<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">a representative length of the porous medium (for example, the average grain size) (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>C<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">appropriate parameter related to the soil type (dimensionless)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">Other more complex relationships (depending on porosity, mean pore diameter, and specific surface area) have been developed for intrinsic permeability of reactive clays, especially if salt is dissolved into the pore water (for example, Raffensperger and Ferrell Jr., 1991).<\/p>\n<p class=\"import-Normal\">Assume the initial conditions are in equilibrium, and all the hydrologic and geomechanical quantities presented here are incremental with respect to the initial conditions. Let\u2019s balance the weight of water in an elementary soil sample of initial length \u0394<em>z<\/em> and unitary cross\u2011sectional area (shown as 1 in the expressions below) between time <em class=\"import-Cambria\">t<\/em> and <em class=\"import-Cambria\">t<\/em>\u00a0+\u00a0\u0394<em class=\"import-Cambria\">t<\/em>:<\/p>\n<ul>\n<li class=\"import-Normal\">Inflow: (<em>\u03b3<\/em><em>v<\/em><sub>z<\/sub>)\u00a0(1)\u00a0(\u0394<em class=\"import-Cambria\">t<\/em>)<\/li>\n<li class=\"import-Normal\">Outflow: <em>\u03b3<\/em>(<em>v<\/em><sub>z<\/sub>\u00a0+\u00a0<em class=\"import-Cambria\">\u2202<\/em><em>v<\/em><sub>z<\/sub>\/<em class=\"import-Cambria\">\u2202<\/em><em>z<\/em><em>\u00a0<\/em>\u0394<em>z<\/em>)\u00a0(1)\u00a0(\u0394<em>t<\/em>)<\/li>\n<li class=\"import-Normal\">Weight of water expelled by the porous space contraction and the expansion of the water expressed by Equation 16 (we assume incompressible solid grains the total medium volume change coincides with the porous volume change):<a id=\"eq_16\"><\/a><\/li>\n<\/ul>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\">\u2212[(<em>\u03b3<\/em>\u2206(<em>\u03d5<\/em>\u2206<em>z<\/em>) 1 \u2206<em>p<\/em>) + (<em>\u03b3\u03d5\u03b2<\/em> 1 \u2206<em>p<\/em>)]<\/td>\n<td style=\"width: 10%; text-align: right;\">(16)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03b2<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">volumetric compressibility of water (ML<sup>-1<\/sup>T<sup>-2<\/sup>)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">In Equation\u00a016 the total geostatic stress <em>\u03c3<\/em><sub><em>c<\/em><\/sub> is assumed to be constant, so (from <a href=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/effective-intergranular-stress-and-soil-parameters#eq_3\">Equation\u00a03<\/a>):<\/p>\n<p class=\"import-Normal\" style=\"text-align: center;\"><span style=\"display: block; font-size: NaNpt; color: #; ; text-decoration: none;\">\u0394<em>\u03c3<\/em><sub><em>z<\/em><\/sub> = \u2212\u2206<em>p<\/em><\/span><\/p>\n<p class=\"import-Normal\">The change in pressure, \u0394<em>p<\/em><em>,<\/em> is negative when <em class=\"import-Cambria\">p<\/em> is reduced, as happens during groundwater pumping. Notice that \u0394(<em>\u03d5<\/em>\u0394<em>z<\/em>) is equal to <em class=\"import-Cambria\">\u2206<\/em>{[<em class=\"import-Cambria\">e<\/em>\/(1\u00a0+\u00a0<em class=\"import-Cambria\">e<\/em>)]<em class=\"import-Cambria\">\u2206z<\/em>} with <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em>\/(1\u00a0+\u00a0<em class=\"import-Cambria\">e<\/em>) constant because this is the solid part (grains) of the elementary volume (1)\u00a0<em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em> . Hence (from <a href=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/effective-intergranular-stress-and-soil-parameters#eq_5\">Equation\u00a05<\/a>) we have:<\/p>\n<p style=\"text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-a28ca588be2cc7fd273e9fe6809adbe6_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#101;&#125;&#123;&#49;&#43;&#101;&#125;&#61;&#99;&#95;&#123;&#98;&#125;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#112;\" title=\"Rendered by QuickLaTeX.com\" height=\"39\" width=\"103\" style=\"vertical-align: -14px;\" \/><\/p>\n<p class=\"import-Normal\">and therefore, we obtain:<\/p>\n<ul>\n<li class=\"import-Normal\">Outflow \u2013 Inflow = Weight of water expelled, i.e.<\/li>\n<\/ul>\n<p style=\"text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-de56374977032a36b217594e0faed036_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#118;&#95;&#123;&#122;&#125;&#125;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#122;&#125;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#122;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#116;&#61;&#45;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#40;&#99;&#95;&#123;&#98;&#125;&#43;&#92;&#112;&#104;&#105;&#32;&#92;&#98;&#101;&#116;&#97;&#32;&#41;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#122;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#112;\" title=\"Rendered by QuickLaTeX.com\" height=\"37\" width=\"253\" style=\"vertical-align: -12px;\" \/><\/p>\n<p class=\"import-Normal\">Cancelling <em>\u03b3<\/em> and <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">z<\/em> on both sides and remembering that the hydraulic head <em class=\"import-Cambria\">h<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">z<\/em>\u00a0+\u00a0<em class=\"import-Cambria\">p<\/em>\/<em class=\"import-Cambria\">\u03b3<\/em>, we know that <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">p<\/em> = <em class=\"import-Cambria\">\u03b3<\/em><em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">h<\/em>, and using Equation\u00a014 when the increment of time approaches zero <em class=\"import-Cambria\">\u2206<\/em><em class=\"import-Cambria\">t<\/em><em class=\"import-Cambria\">\u00a0<\/em><em class=\"import-Cambria\">\u2192<\/em><em class=\"import-Cambria\">\u00a0<\/em>0 we obtain Equation\u00a017:<a id=\"eq_17\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-58b0b19f166a7772ba828d8376c7daa3_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#125;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#122;&#125;&#92;&#108;&#101;&#102;&#116;&#32;&#40;&#32;&#75;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#104;&#125;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#122;&#125;&#32;&#92;&#114;&#105;&#103;&#104;&#116;&#32;&#41;&#61;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#40;&#99;&#95;&#123;&#98;&#125;&#43;&#92;&#112;&#104;&#105;&#32;&#92;&#98;&#101;&#116;&#97;&#32;&#41;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#104;&#125;&#123;&#92;&#112;&#97;&#114;&#116;&#105;&#97;&#108;&#32;&#116;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"43\" width=\"219\" style=\"vertical-align: -17px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(17)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">Solving Equation\u00a017, complemented with the appropriate top and bottom boundary conditions and initial conditions, provides the pressure dissipation within the aquitard, and hence the <em class=\"import-Cambria\">\u0394p<\/em> needed to compute the aquitard compaction versus time. The specific storage coefficient is defined in Equation\u00a018.<a id=\"eq_18\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><em>S<\/em><sub><em>s<\/em><\/sub> = <em>\u03b3<\/em>(<em>c<\/em><sub><em>b<\/em><\/sub> + <em>\u03d5\u03b2<\/em>)<\/td>\n<td style=\"width: 10%; text-align: right;\">(18)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\"><em>S<\/em><sub><em>s<\/em><\/sub> represents the \u201c<em>specific elastic storage<\/em>\u201d [L<sup>\u2011<\/sup><sup>1<\/sup>] and along with the hydraulic conductivity, <em class=\"import-Cambria\">K<\/em>, defines Terzaghi\u2019s consolidation coefficient <em>c<\/em><sub><em>v<\/em><\/sub> that controls both magnitude and timing of aquitard compaction as shown in Equation\u00a019 given that <em class=\"import-Cambria\">\u03d5\u03b2\u00a0<\/em>&lt;&lt;\u00a0<em>c<\/em><sub><em>b<\/em><\/sub> for typical aquifer confining beds.<a id=\"eq_19\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-7d65e359d34e0b743344ada6b3928941_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#99;&#95;&#123;&#118;&#125;&#61;&#92;&#102;&#114;&#97;&#99;&#123;&#75;&#125;&#123;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#40;&#99;&#95;&#123;&#98;&#125;&#43;&#92;&#112;&#104;&#105;&#32;&#92;&#98;&#101;&#116;&#97;&#32;&#41;&#125;&#92;&#97;&#112;&#112;&#114;&#111;&#120;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#75;&#125;&#123;&#92;&#103;&#97;&#109;&#109;&#97;&#32;&#99;&#95;&#123;&#98;&#125;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"41\" width=\"175\" style=\"vertical-align: -17px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(19)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">The initial conditions correspond to \u0394<em class=\"import-Cambria\">p<\/em>\u00a0=\u00a00 for the entire thickness, <em class=\"import-Cambria\">b<\/em>, of the aquitard while the boundary conditions are given by \u0394<em class=\"import-Cambria\">p<\/em> in the overlying and underlying aquifers. If the pressure drop \u0394<em class=\"import-Cambria\">p<\/em><sub class=\"import-Cambria\">0<\/sub> is the same at top (<em class=\"import-Cambria\">z<\/em>\u00a0=\u00a00) and bottom (<em class=\"import-Cambria\">z<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">b<\/em>), then pressure conditions are symmetrical above and below the middle of the aquitard, hence <span style=\"font-size: NaNpt; color: #; ; text-decoration: none;\">\u2202<em>p<\/em>\/\u2202<em>z<\/em> = 0 <\/span>at <em class=\"import-Cambria\">z<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">b<\/em>\/2. In this case, the solution to Equation\u00a017 can be expressed in terms of a series expansion. Writing the solution in terms of <em class=\"import-Cambria\">p<\/em> produces Equation\u00a020.<a id=\"eq_20\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-dcb49b2cc6fe5a30c6b0ceed53bd18a6_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#112;&#61;&#112;&#95;&#123;&#48;&#125;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#52;&#125;&#123;&#92;&#112;&#105;&#32;&#125;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#112;&#95;&#123;&#48;&#125;&#92;&#108;&#101;&#102;&#116;&#32;&#91;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#105;&#32;&#125;&#123;&#52;&#125;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#115;&#105;&#110;&#125;&#40;&#92;&#112;&#105;&#32;&#122;&#47;&#98;&#41;&#125;&#123;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#101;&#120;&#112;&#125;&#91;&#40;&#92;&#112;&#105;&#32;&#47;&#98;&#41;&#94;&#123;&#50;&#125;&#99;&#95;&#123;&#118;&#125;&#116;&#93;&#125;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#51;&#125;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#115;&#105;&#110;&#125;&#40;&#51;&#92;&#112;&#105;&#32;&#122;&#47;&#98;&#41;&#125;&#123;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#101;&#120;&#112;&#125;&#91;&#40;&#51;&#92;&#112;&#105;&#32;&#47;&#98;&#41;&#94;&#123;&#50;&#125;&#99;&#95;&#123;&#118;&#125;&#116;&#93;&#125;&#45;&#46;&#46;&#46;&#32;&#92;&#114;&#105;&#103;&#104;&#116;&#32;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"43\" width=\"475\" style=\"vertical-align: -17px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(20)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>z<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">vertical coordinate positive downward starting from the aquitard top (0 \u2264 <em>z<\/em> \u2264 <em>b<\/em>\/2) (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>t<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">time since the initial change in pressure at the aquitard boundary (T)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">For <em class=\"import-Cambria\">t<\/em> = 0, and setting <em class=\"import-Cambria\">\u03c0z<\/em>\/<em class=\"import-Cambria\">b<\/em>\u00a0=\u00a0<em class=\"import-Cambria\">x<\/em>\u00a0(0\u00a0\u2264\u00a0<em class=\"import-Cambria\">x<\/em>\u00a0\u2264\u00a0<em class=\"import-Cambria\">\u03c0<\/em>\/2) Equation\u00a020 becomes Equation\u00a021.<a id=\"eq_21\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-8f7e4ebf0a3bb45051f307117fbe48da_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#112;&#40;&#122;&#44;&#48;&#41;&#61;&#112;&#95;&#123;&#48;&#125;&#40;&#122;&#41;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#52;&#125;&#123;&#92;&#112;&#105;&#32;&#125;&#92;&#68;&#101;&#108;&#116;&#97;&#32;&#112;&#95;&#123;&#48;&#125;&#92;&#108;&#101;&#102;&#116;&#32;&#91;&#32;&#92;&#102;&#114;&#97;&#99;&#123;&#92;&#112;&#105;&#32;&#125;&#123;&#52;&#125;&#45;&#92;&#108;&#101;&#102;&#116;&#32;&#40;&#32;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#115;&#105;&#110;&#125;&#120;&#43;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#51;&#125;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#115;&#105;&#110;&#125;&#51;&#120;&#43;&#92;&#102;&#114;&#97;&#99;&#123;&#49;&#125;&#123;&#53;&#125;&#92;&#116;&#101;&#120;&#116;&#117;&#112;&#123;&#115;&#105;&#110;&#125;&#53;&#120;&#43;&#46;&#46;&#46;&#92;&#114;&#105;&#103;&#104;&#116;&#32;&#41;&#32;&#92;&#114;&#105;&#103;&#104;&#116;&#32;&#93;\" title=\"Rendered by QuickLaTeX.com\" height=\"43\" width=\"479\" style=\"vertical-align: -17px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(21)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">The content within the parentheses of Equation\u00a021 is the Fourier series development of the function <em class=\"import-Cambria\">f<\/em><em class=\"import-Cambria\">\u00a0<\/em>(<em class=\"import-Cambria\">x<\/em>):<\/p>\n<p style=\"text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-content\/ql-cache\/quicklatex.com-797c007d2dc90922d7b55bd48a06e65d_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#102;&#40;&#120;&#41;&#61;&#92;&#98;&#101;&#103;&#105;&#110;&#123;&#99;&#97;&#115;&#101;&#115;&#125; &#38;&#32;&#92;&#32;&#92;&#32;&#92;&#112;&#105;&#32;&#47;&#52;&#92;&#116;&#101;&#120;&#116;&#123;&#32;&#102;&#111;&#114;&#32;&#125;&#32;&#48;&#92;&#108;&#101;&#113;&#32;&#120;&#60;&#32;&#92;&#112;&#105;&#32;&#92;&#92; &#38;&#32;&#45;&#92;&#112;&#105;&#32;&#47;&#52;&#92;&#116;&#101;&#120;&#116;&#123;&#32;&#102;&#111;&#114;&#32;&#125;&#32;&#92;&#112;&#105;&#32;&#60;&#32;&#120;&#92;&#108;&#101;&#113;&#32;&#50;&#92;&#112;&#105; &#92;&#101;&#110;&#100;&#123;&#99;&#97;&#115;&#101;&#115;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"54\" width=\"254\" style=\"vertical-align: -23px;\" \/><\/p>\n<p class=\"import-Normal\">Keeping in mind the range of variability of <em class=\"import-Cambria\">x\u00a0<\/em>(0\u00a0\u2264\u00a0<em class=\"import-Cambria\">x<\/em><span lang=\"pt-BR\" xml:lang=\"pt-BR\">\u00a0<\/span>\u2264\u00a0<em class=\"import-Cambria\">\u03c0<\/em>\/2) we conclude that Equation\u00a020 accurately represents the initial pore pressure at time <em class=\"import-Cambria\">t<\/em>\u00a0=\u00a00.<\/p>\n<p class=\"import-Normal\">\n<\/div>\n","protected":false},"author":1,"menu_order":10,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-238","chapter","type-chapter","status-publish","hentry"],"part":121,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/238","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/users\/1"}],"version-history":[{"count":8,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/238\/revisions"}],"predecessor-version":[{"id":487,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/238\/revisions\/487"}],"part":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/parts\/121"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/238\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/media?parent=238"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapter-type?post=238"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/contributor?post=238"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/license?post=238"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}