scholarly journals SOLUTE TRANSPORT MODELLING IN LOW-PERMEABILITY HOMOGENEOUS AND SATURATED SOIL MEDIA

2021 ◽  
Vol 36 (2) ◽  
pp. 25-32
Author(s):  
Muhammad Zaheer ◽  
Hadayat Ullah ◽  
Saad Ahmed Mashwani ◽  
Ehsan ul Haq ◽  
Syed Husnain Ali Shah ◽  
...  

Fickian and non-Fickian behaviors were often detected for contaminant transport activity owed to the preferential flow and heterogeneity of soil media. Therefore, using diverse methods to measure such composite solute transport in soil media has become an important research topic for solute transport modeling in soil media. In this article, the continuous-time random walk (CTRW) model was applied to illustrate the relative concentration of transport in low-permeability homogeneous and saturated soil media. The solute transport development was also demonstrated with the convection-dispersion equation (CDE) and Two Region Model (TRM) for comparison. CXTFIT 2.1 software was used for CDE and TRM, and CTRW Matlab Toolbox v.3.1 for the CTRW simulation of the breakthrough curve. It was found that higher values of determination coefficient (R2) and lower values of root mean square error (RMSE) concerning the best fits of CDE, TRM, and CTRW. It was found that in the comparison of CDE, TRM, and CTRW, we tend to use CTRW to describe the transport behavior well because there are prevailing Fickian and non-Fickian transport. The CTRW gives better fitting results to the breakthrough curves (BTCs) when β has an increasing pattern towards 2.00. In this study, the variation of parameters in three methods was investigated and results showed that the CTRW modeling approach is more effective to determine non-reactive contaminants concentration in low-permeability soil media at small depths.

2015 ◽  
Vol 12 (1) ◽  
pp. 67-92
Author(s):  
R. A. Crane ◽  
M. O. Cuthbert ◽  
W. Timms

Abstract. We present an interrupted-flow centrifugation technique to characterise preferential flow in low permeability media. The method entails a minimum of three phases: centrifuge induced flow, no flow and centrifuge induced flow, which may be repeated several times in order to most effectively characterise multi-rate mass transfer behaviour. In addition, the method enables accurate simulation of relevant in situ total stress conditions during flow by selecting an appropriate centrifugal force level. We demonstrate the utility of the technique for characterising the hydraulic properties of smectite clay dominated core samples. All samples exhibited a non-Fickian tracer breakthrough (early tracer arrival), combined with a decrease in tracer concentration immediately after each period of interrupted-flow. This is indicative of dual (or multi) porosity behaviour, with solute migration predominately via advection during induced flow, and via molecular diffusion (between the preferential flow network(s) and the low hydraulic conductivity domain) during interrupted-flow. Tracer breakthrough curves were simulated using a bespoke dual porosity model with excellent agreement between the data and model output (Nash–Sutcliffe model efficiency coefficient was >0.97 for all samples). In combination interrupted-flow centrifuge experiments and dual porosity transport modelling are shown to be a powerful method to characterise preferential flow in low permeability media.


2015 ◽  
Vol 19 (9) ◽  
pp. 3991-4000 ◽  
Author(s):  
R. A. Crane ◽  
M. O. Cuthbert ◽  
W. Timms

Abstract. We present an interrupted-flow centrifugation technique to characterise preferential flow in low permeability media. The method entails a minimum of three phases: centrifuge-induced flow, no flow and centrifuge-induced flow, which may be repeated several times in order to most effectively characterise multi-rate mass transfer behaviour. In addition, the method enables accurate simulation of relevant in situ total stress conditions during flow by selecting an appropriate centrifugal force. We demonstrate the utility of the technique for characterising the hydraulic properties of smectite-clay-dominated core samples. All core samples exhibited a non-Fickian tracer breakthrough (early tracer arrival), combined with a decrease in tracer concentration immediately after each period of interrupted flow. This is indicative of dual (or multi-)porosity behaviour, with solute migration predominately via advection during induced flow, and via molecular diffusion (between the preferential flow network(s) and the low hydraulic conductivity domain) during interrupted flow. Tracer breakthrough curves were simulated using a bespoke dual porosity model with excellent agreement between the data and model output (Nash–Sutcliffe model efficiency coefficient was > 0.97 for all samples). In combination, interrupted-flow centrifuge experiments and dual porosity transport modelling are shown to be a powerful method to characterise preferential flow in low permeability media.


Water ◽  
2020 ◽  
Vol 12 (7) ◽  
pp. 2049 ◽  
Author(s):  
Yingtao Hu ◽  
Wenjie Xu ◽  
Liangtong Zhan ◽  
Zuyang Ye ◽  
Yunmin Chen

The influence of contact area, caused by normal deformation, on fluid flow and solute transport through three-dimensional (3D) rock fractures is investigated. Fracture surfaces with different Hurst exponents (H) were generated numerically using the modified successive random addition (SRA) method. By applying deformations normal to the fracture surface (Δu), a series of fracture models with different aperture distributions and contact area ratios (c) were simulated. The results show that the contact area between the two fracture surfaces increases and more void spaces are reduced as deformation (Δu) increases. The streamlines in the rough-walled fractures show that the contact areas result in preferential flow paths and fingering type transport. The non-Fickian characteristics of the “early arrival” and “long tail” in all of the breakthrough curves (BTCs) for fractures with different deformation (Δu) and Hurst parameters (H) were determined. The solute concentration distribution index (CDI), which quantifies the uniformity of the concentration distribution within the fracture, decreases exponential as deformation (Δu) and/or contact area ratios (c) increase, indicating that increased contact area can result in a larger delay rate of mass exchange between the immobile zone around the contact areas and the main flow channel, thus, resulting in a longer time for the solute to fill the entire fracture. The BTCs were analyzed using the continuous time random walk (CTRW) inverse model. The inverse modeling results show that the dispersion exponent β decreases from 1.92 to 0.81 as c increases and H decreases, suggesting that the increase in contact area and fracture surfaces enhance the magnitude of the non-Fickian transport.


2014 ◽  
Vol 139 (6) ◽  
pp. 634-641 ◽  
Author(s):  
Tyler C. Hoskins ◽  
James S. Owen ◽  
Jeb S. Fields ◽  
James E. Altland ◽  
Zachary M. Easton ◽  
...  

An understanding of how dissolved mineral nutrient ions (solutes) move through pine bark substrates during the application of irrigation water is vital to better understand nutrient transport and leaching from containerized crops during an irrigation event. However, current theories on solute transport processes in soilless systems are largely based on research in mineral soils and thus do not necessarily explain solute transport in soilless substrates. A study was conducted to characterize solute transport through a 9 pine bark:1 sand (by volume) substrate by developing and analyzing breakthrough curves (BTCs). Columns filled with pine bark substrate were subjected to the application of a nutrient solution (tracer) and deionized water under saturated and unsaturated conditions. Effluent drained from the columns during these applications was collected and analyzed to determine the effluent concentration (C) of the bulk ions in solution through electrical conductivity (EC) and nitrate (NO3–), phosphate, and potassium (K+) concentrations. The BTCs were developed by plotting C relative to the concentration of the input solution (Co) (i.e., relative concentration = C/Co) as a function of the cumulative effluent volume. Solutes broke through the column earlier (i.e., with less cumulative effluent) and the transition from C/Co = 0 to 1 occurred more abruptly under unsaturated than saturated conditions. Movement of the anion, NO3–, through the substrate was observed to occur more quickly than the cation K+. Throughout the experiment, 37% of the applied K+ was retained by the pine bark. The adsorption of K+ to pine bark cation exchange sites displaced calcium (Ca2+) and magnesium (Mg2+), of which the combined equivalent charge accounted for 43.1% of the retained K+. These results demonstrate the relative ease that negatively charged fertilizer ions could move through a pine bark substrate while solution is actively flowing through substrate pores such as during irrigation events. This approach to evaluating solute transport may be used in horticultural research to better understand how mineral nutrients move through and subsequently leach from soilless substrates during irrigation. Expanding this knowledge base may lead to the refinement of production practices that improve nutrient and water use efficiency in container nurseries.


2005 ◽  
Vol 69 (2) ◽  
pp. 291-300 ◽  
Author(s):  
Young-Jin Kim ◽  
Christophe J. G. Darnault ◽  
Nathan O. Bailey ◽  
J.-Yves Parlange ◽  
Tammo S. Steenhuis

2021 ◽  
pp. 126397
Author(s):  
Jiaxu Zhang ◽  
Chunhui Lu ◽  
Chengji Shen ◽  
Chenming Zhang ◽  
Jun Kong ◽  
...  

2009 ◽  
Vol 13 (6) ◽  
pp. 935-944 ◽  
Author(s):  
A. E. Anderson ◽  
M. Weiler ◽  
Y. Alila ◽  
R. O. Hudson

Abstract. Preferential flow paths have been found to be important for runoff generation, solute transport, and slope stability in many areas around the world. Although many studies have identified the particular characteristics of individual features and measured the runoff generation and solute transport within hillslopes, very few studies have determined how individual features are hydraulically connected at a hillslope scale. In this study, we used dye staining and excavation to determine the morphology and spatial pattern of a preferential flow network over a large scale (30 m). We explore the feasibility of extending small-scale dye staining techniques to the hillslope scale. We determine the lateral preferential flow paths that are active during the steady-state flow conditions and their interaction with the surrounding soil matrix. We also calculate the velocities of the flow through each cross-section of the hillslope and compare them to hillslope scale applied tracer measurements. Finally, we investigate the relationship between the contributing area and the characteristics of the preferential flow paths. The experiment revealed that larger contributing areas coincided with highly developed and hydraulically connected preferential flow paths that had flow with little interaction with the surrounding soil matrix. We found evidence of subsurface erosion and deposition of soil and organic material laterally and vertically within the soil. These results are important because they add to the understanding of the runoff generation, solute transport, and slope stability of preferential flow-dominated hillslopes.


Sign in / Sign up

Export Citation Format

Share Document