{"id":133,"date":"2024-12-12T23:56:31","date_gmt":"2024-12-12T23:56:31","guid":{"rendered":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/chapter\/2-7-mechanisms-of-contaminants-transport-through-soil\/"},"modified":"2026-03-16T13:55:23","modified_gmt":"2026-03-16T13:55:23","slug":"2-7-mechanisms-of-contaminants-transport-through-soil","status":"publish","type":"chapter","link":"https:\/\/libraryresources.nse.org.ng\/fundamentalsoffoundationengineering\/chapter\/2-7-mechanisms-of-contaminants-transport-through-soil\/","title":{"raw":"2.7 Mechanisms of contaminants transport through soil","rendered":"2.7 Mechanisms of contaminants transport through soil"},"content":{"raw":"When a pollutant escapes from a repository, it moves through the vadose zone, <\/em>which is the zone between the soil surface and the groundwater table (Figure 2.10<\/a>). The soil grains retain part of the pollutant, either chemically or mechanically, while part of the pollutant reaches the aquifer. Water-soluble pollutants are then carried away downstream, following the groundwater flow, which is determined by the hydraulic gradient.\n\nAs pollutants are being carried away, contamination of the aquifer spreads and substance concentration decreases, as we move from the main zone to the neutral zone (Figure 2.10<\/a>). Especially for saturated or partially saturated soils and pollutants that can be dissolved or suspended in water (light NAPLs that float are excluded), three basic mechanisms can be identified:\n

2.7.1 Advection or Hydraulic transportation<\/h1>\nThe contaminant is carried away by the groundwater though the soil voids, following the hydraulic gradient (Figure 2.10<\/a>). Assuming saturated flow, the seepage velocity will be equal to:\n\n(2.1)<\/strong> [latex]V =\\dfrac{q}{n}[\/latex]\n\n[caption id=\"attachment_131\" align=\"aligncenter\" width=\"900\"]\"Diagram Figure 2.10. Pollutant movement through the subsoil.[\/caption]\n\nwhere n<\/em> is the soil porosity (n=Vvoids<\/sub><\/em>\/V<\/em>) and q<\/em> is Darcy's water flux, equal to the volume Q<\/em> flowing though a unit cross-sectional area A<\/em>:\n\n(2.2) <\/strong>[latex]q = \\dfrac{Q}{A} = - k\\dfrac{{\\partial h}}{{\\partial x}} = ki[\/latex]\n\nIn Eq. 2.1 k<\/em> is the permeability of the aquifer, h<\/em> is the hydraulic head, i<\/em> is the hydraulic gradient and x<\/em> defines the flow direction (Figure 2.10<\/a>). Based on the above the adjective mass flux of the contaminant, defined as the mass flowing through a unit cross-sectional area per time (units: mass\/(time length2<\/sup>)), is given as:\n\n(2.3) <\/strong>[latex]{J_{advection}} = qc = nVc[\/latex]<\/em>\n\nwhere c<\/em> is the concentration of the contaminant in the liquid phase. The concentration c<\/em> in a specific location constantly varies with time, while concentration in a specific volume of water remains constant, as this volume of water flows following the hydraulic gradient.\n\n
\n\n

2.7.2 Diffusion or molecular diffusion<\/h1>\n[caption id=\"attachment_131\" align=\"aligncenter\" width=\"600\"]\"Illustration<\/a> Figure 2.11. Diffusion transportation mechanism from a microscopic and macroscopic point of view.[\/caption]\n\nContrary to advection transport, which refers to macroscopic movement, diffusion is a molecular based phenomenon. The contaminant is diffused in the groundwater, due to variation in concentration, meaning that the substance moves from areas of high to areas of low concentration. Figure 2.11<\/a> illustrates this from a microscopic and macroscopic point of view. Initially there are solute molecules (or ions) of the contaminant only on the left side of a barrier. These molecules or ions move randomly. As a result of this motion, if we remove the barrier, the solute diffuses to fill the whole container, moving from the high concentration area (left) to the low concentration area (right). Diffusion will continue until the concentration of the contaminant becomes equal everywhere. The rate at which these irregularities in concentration disappear can be quantified via Fick's first law of diffusion:\n\n(2.4) <\/strong>[latex]{J_{diffusion}} = - D\\dfrac{{\\partial c}}{{\\partial x}}[\/latex]<\/i>\n\nEq. 2.4 suggests that the diffusive mass flux of the contaminant Jdiffusion<\/sub><\/em>, which provides the mass of the contaminant flowing through a unit cross-sectional area per time (units: mass\/(time length2<\/sup>)), depends on the difference in concentration, the inverse of the distance, and a diffusion coefficient D<\/em>. This is a function of the type of cation, temperature, soil porosity and soil tortuosity. In other words, diffusion of the contaminant does not depend on groundwater flow, and can take place in still water too.\n\n
\n\n

2.7.3 Dispersion or mechanical dispersion<\/h1>\nMovement of the contaminant is due to the presence of interconnected voids of the soil skeleton, of various shapes. Thus in microscopic scale, the velocity of water moving though the voids varies from place to place (Figure 2.12<\/a>), and may diverge from the mean (macroscopic) groundwater flow i.e., the substance can be transported perpendicularly to the macroscopic groundwater flow too. Similar to the diffusion mechanism described above, the dispersive mass flux of the contaminant Jdispersion<\/sub><\/em> provides the mass of the contaminant flowing through a unit cross-sectional area per time:\n\n(2.5) <\/strong>[latex]{J_{dispersion}} = - {D_d}\\dfrac{{\\partial c}}{{\\partial x}} [\/latex]<\/i>\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0<\/em>\n\nJ<\/em>dispersion<\/em><\/sub> (units: mass\/(time length2<\/sup>)) depends again on the difference in concentration, the inverse of the distance, and a dispersion coefficient Dd<\/sub><\/em>, which is a function of the seepage velocity and of the longitudinal dispersivity.\n\nOften, these mechanisms may act antagonistically (Figure 2.13<\/a>) or simultaneously (Figure 2.14<\/a>). More specifically, the advection mechanism tends to transport the contaminant from right to left, from the area of high hydraulic potential to the area of lower hydraulic potential, following the hydraulic gradient. The diffusion mechanism on the other hand tends to transport the contaminant from the area of high concentration to the area of low concentration (Figures 2.13<\/a> and 2.14<\/a>).\n\n[caption id=\"attachment_131\" align=\"aligncenter\" width=\"600\"]\"Diagram Figure 2.12. Mechanical dispersion of contaminant through soil.[\/caption]\n\nTherefore to calculate the evolution of contaminant concentration with time we need to solve the continuity Eq. 2.6, which describes mathematically the transport of a conserved quantity (here the mass of the contaminant):\n\n(2.6)<\/strong> [latex]\\dfrac{{\\partial m}}{{\\partial t}} + \\nabla J = S [\/latex]\n\nwhere m <\/em>is the mass density (mass of contaminant per unit volume), J <\/em>is the mass flux vector considering the above mechanisms and S<\/em> quantifies the increase (S<\/em> > 0) or reduction (S<\/em> < 0) of the contaminant mass. Positive S<\/em> terms are referred to as \u201csources\u201d and negative S<\/em> terms are referred to as \u201csinks\u201d.\n\nThe gradient of the mass flux vector:\n\n(2.7)<\/strong> [latex]\\nabla J = \\dfrac{{\\partial J}}{{\\partial x}}i + \\dfrac{{\\partial J}}{{\\partial y}}j + \\dfrac{{\\partial J}}{{\\partial z}}k,{\\rm{with \\: \\:}}\\hat i = \\left[ {\\begin{array}{*{20}{c}}1\\\\0\\\\0\\end{array}} \\right]{\\rm{ }}\\hat j = \\left[ {\\begin{array}{*{20}{c}}0\\\\1\\\\0\\end{array}} \\right]{\\rm{ }}\\hat k = \\left[ {\\begin{array}{*{20}{c}}0\\\\0\\\\1\\end{array}} \\right] [\/latex]\n\nwill provide the direction towards which the mass flux vector rises faster in the three-dimensional space, and its magnitude provides a measure of how fast the mass flux rises in this direction. For example, if the mass flux vector due to advection and diffusion is J = -3x-4y <\/em>where x <\/em>is the horizontal axis and y <\/em>is the vertical axis, then:\n\n(2.8)<\/strong> [latex]\\nabla J = - 3i - 4j,{\\rm{and \\:}}\\left| J \\right| = \\sqrt {{3^2} + {4^2}} = 5\\left( {\\dfrac{{{\\rm{mass}}}}{{{\\rm{time}} \\times {{\\rm{length}}^{\\rm{2}}}}}} \\right) [\/latex]\n\n[caption id=\"attachment_131\" align=\"aligncenter\" width=\"600\"]\"Diagram Figure 2.13. Antagonistic action of two mechanisms of contaminant transport.[\/caption]\n\n[caption id=\"attachment_131\" align=\"aligncenter\" width=\"600\"]\"Diagram<\/a> Figure 2.14. Combined action of two mechanisms of contaminant transport.[\/caption]\n\nOther, non-mechanical mechanisms may also act to reduce and disperse the pollution load, and are introduced as \u201csinks\u201d in Eq. 2.6:\n