A multi-flowpath model for the interpretation of immiscible displacement experiments in heterogeneous soil columns

This work focuses on the phenomenon of the immiscible two-phase flow of water and oil in saturated heterogeneous soil columns. The goal is to develop a fast and reliable method for quantifying soil heterogeneities for incorporation into the relevant capillary pressure and relative permeability functions. Such data are commonly used as input data in simulators of contaminant transport in the subsurface. Rate-controlled drainage experiments are performed on undisturbed soil columns and the transient response of the axial distribution of water saturation is determined from electrical measurements. The transient responses of the axial distribution of water saturation and total pressure drop are fitted with the multi-flowpath model (MFPM) where the pore space is regarded as a system of parallel paths of different permeability. The MFPM enables us to quantify soil heterogeneity at two scales: the micro-scale parameters describe on average the effects of pore network heterogeneities on the two-phase flow pattern; the macro-scale parameters indicate the variability of permeability at the scale of interconnected pore networks. The capillary pressure curve is consistent with that measured with mercury intrusion porosimetry over the low pressure range. The oil relative permeability increases sharply at a very low oil saturation (< 10^− 3) and tends to a high end value. The water relative permeability decreases abruptly at a low oil saturation (~ 0.1), whereas the irreducible wetting phase saturation is quite high. The foregoing characteristics of the two-phase flow properties are associated with critical (preferential) flowpaths that comprise a very small percentage of the total pore volume, control the overall hydraulic conductivity, and are consistent with the very broad range of pore-length scales usually probed in soil porous matrix.     

Influence of soil texture and tillage on herbicide transport

Two long-term no-till corn production studies, representing different soil texture, consistently showed higher leaching of atrazine [2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine] to groundwater in a silt loam soil than in a sandy loam soil. A laboratory leaching study was initiated using intact soil cores from the two sites to determine whether the soil texture could account for the observed differences. Six intact soil cores (16 cm dia by 20 cm high) were collected from a four-year old no-till corn plots at each of the two locations (ca. 25 km apart). All cores were mounted in funnels and the saturated hydraulic conductivity (Ksat) was measured. Three cores (from each soil texture) with the lowest Ksat were mixed and repacked. All cores were surface treated with 1.7 kg ai ha(-1) [ring-14C] atrazine, subjected to simulated rainfall at a constant 12 mm h(-1) intensity until nearly 3 pore volume of leachate was collected and analyzed for a total of 14C. On an average, nearly 40% more of atrazine was leached through the intact silt loam than the sandy loam soil cores. For both the intact and repacked cores, the initial atrazine leaching rates were higher in the silt loam than the sandy loam soils, indicating that macropore flow was a more prominent mechanism for atrazine leaching in the silt loam soil. A predominance of macropore flow in the silt loam soil, possibly due to greater aggregate stability, may account for the observed leaching patterns for both field and laboratory studies.     

Evaluating non-equilibrium solute transport in small soil columns

Displacement studies on leaching of bromide and two pesticides (atrazine and isoproturon) were conducted under unsaturated steady state flow conditions in 24 small undisturbed soil columns (5.7 cm in diameter and 10 cm long) each collected from two sites differing in soil structure and organic carbon content in North Germany. There were large and irregular variabilities in the characteristics of both soils, as well as in the shapes of breakthrough curves (BTCs) of different columns, including some with early breakthrough and increased tailing, qualitatively indicating the presence of preferential flow. It was estimated that one preferential flow column (PFC) at site A, and four at site B, contributed, respectively to 11% and 58% of the accumulated leached fraction and to more than 80% of the maximum observed standard deviation (SD) in the field-scale concentration and mass flux of pesticides at two sites. The bromide BTCs of two sites were analyzed with the equilibrium convection-dispersion equation (CDE) and a non-equilibrium two-region/mobile-immobile model. Transport parameters of these models for individual BTCs were determined using a curve fitting program, CXTFIT, and by the time moment method. For the CDE based equilibrium model, the mean values of retardation factor, R, considered separately for all columns, PFCs or non-preferential flow columns (NPFCs) were comparable for the two methods; significant differences were observed in the values of dispersion coefficients of two sites using the two estimation methods. It was inferred from the estimated parameters of non-equilibrium model that 5-12% of water at site A, and 12% at site B, was immobile during displacement in NPFCs. The corresponding values for PFCs of two sites were much larger, ranging from 25% to 51% by CXTFIT and from 24% to 72% by the moment method, suggesting the role of certain mechanisms other than immobile water in higher degrees of non-equilibrium in these columns. Peclet numbers in PFCs of both sites were consistently smaller than five, indicating the inadequacy of the non-equilibrium model to incorporate the effect of all forms of non-equilibrium in PFCs. Overall, the BTCs of individual NPFCs, PFCs and of field average concentration at the two sites were better reproduced with parameters obtained from CXTFIT than by the moment method. The moment method failed to capture the peak concentrations in PFCs, but tended to describe the desorption and tail branches of BTCs better than the curve fitting approach.     

Stability of titania nanoparticles in soil suspensions and transport in saturated homogeneous soil columns

The stability of TiO₂ nanoparticles in soil suspensions and their transport behavior through saturated homogeneous soil columns were studied. The results showed that TiO₂ could remain suspended in soil suspensions even after settling for 10 days. The suspended TiO₂ contents in soil suspensions after 24 h were positively correlated with the dissolved organic carbon and clay content of the soils, but were negatively correlated with ionic strength, pH and zeta potential. In soils containing soil particles of relatively large diameters and lower solution ionic strengths, a significant portion of the TiO₂ (18.8-83.0%) readily passed through the soils columns, while TiO₂ was significantly retained by soils with higher clay contents and salinity. TiO₂ aggregate sizes in the column outflow significantly increased after passing through the soil columns. The estimated transport distances of TiO₂ in some soils ranged from 41.3 to 370 cm, indicating potential environmental risk of TiO₂ nanoparticles to deep soil layers.     

Modelling tritium and phosphorus transport by preferential flow in structured soil

Subsurface solute transport through structured soil is studied by model interpretation of experimental breakthrough curves from tritium and phosphorus tracer tests in three intact soil monoliths. Similar geochemical conditions, with nearly neutral pH, were maintained in all the experiments. Observed transport differences for the same tracer are thus mainly due to differences in the physical transport process between the different monoliths. The modelling is based on a probabilistic Lagrangian approach that decouples physical and chemical mass transfer and transformation processes from pure and stochastic advection. Thereby, it enables explicit quantification of the physical transport process through preferential flow paths, honouring all independently available experimental information. Modelling of the tritium breakthrough curves yields a probability density function of non-reactive solute travel time that is coupled with a reaction model for linear, non-equilibrium sorption–desorption to describe the phosphorus transport. The tritium model results indicate that significant preferential flow occurs in all the experimental soil monoliths, ranging from 60–100% of the total water flow moving through only 25–40% of the total water content. In agreement with the fact that geochemical conditions were similar in all experiments, phosphorus model results yield consistent first-order kinetic parameter values for the sorption–desorption process in two of the three soil monoliths; phosphorus transport through the third monolith cannot be modelled because the apparent mean transport rate of phosphorus is anomalously rapid relative to the non-adsorptive tritium transport. The occurrence of preferential flow alters the whole shape of the phosphorus breakthrough curve, not least the peak mass flux and concentration values, and increases the transported phosphorus mass by 2–3 times relative to the estimated mass transport without preferential flow in the two modelled monoliths.     

Assessing tungsten transport in the vadose zone: from dissolution studies to soil columns

This study investigates the dissolution, sorption, leachability, and plant uptake of tungsten and alloying metals from canister round munitions in the presence of model, well characterized soils. The source of tungsten was canister round munitions, composed mainly of tungsten (95%) with iron and nickel making up the remaining fraction. Three soils were chosen for the lysimeter studies while four model soils were selected for the adsorption studies. Lysimeter soils were representatives of the typical range of soils across the continental USA; muck-peat, clay-loamy and sandy-quartzose soil. Adsorption equilibrium data on the four model soils were modeled with Langmuir and linear isotherms and the model parameters were obtained. The adsorption affinity of soils for tungsten follows the order: Pahokee peat > kaolinite > montmorillonite > illite. A canister round munition dissolution study was also performed. After 24 d, the measured dissolved concentrations were: 61.97, 3.56, 15.83 mg L⁻¹ for tungsten, iron and nickel, respectively. Lysimeter transport studies show muck peat and sandy quartzose soils having higher tungsten concentration, up to 150 mg kg⁻¹ in the upper layers of the lysimeters and a sharp decline with depth suggesting strong retardation processes along the soil profile. The concentrations of tungsten, iron and nickel in soil lysimeter effluents were very low in terms of posing any environmental concern; although no regulatory limits have been established for tungsten in natural waters. The substantial uptake of tungsten and nickel by ryegrass after 120 d of exposure to soils containing canister round munition suggests the possibility of tungsten and nickel entering the food chain.     

Leaching of 2,4-D from a soil in the presence of beta-cyclodextrin: laboratory columns experiments

This study reports on the effect of the presence of beta-cyclodextrin (beta-CD) on the adsorption and mobility of the pesticide 2,4-dichlorophenoxyacetic acid (2,4-D) through soil columns. The previous application of beta-CD to the soil produced a retarded leaching of 2,4-D through the soil column, due probably to herbicide adsorption on the soil through beta-CD adsorbed. However, the application of beta-CD solution to the soil column where 2,4-D had been previously adsorbed, led to the complete desorption of the herbicide, due to the formation of water-soluble 1:1 inclusion complexes between 2,4-D and beta-CD. Beta-CD can be viewed as a microscopic organic-phase extractant. It can be an advantage to remove from soil pesticides which are able to form inclusion complexes with cyclodextrins, making them possible candidates for use in in situ remediation efforts.     

Characterization of preferential flow in undisturbed, structured soil columns using a vertical TDR probe

Rapid movement of agricultural chemicals through soil to groundwater via preferential flow pathways is one cause of water contamination. Previous studies have shown that time domain reflectometry (TDR) could be used to characterize solute transport in soil. However, previous studies have only scarcely addressed preferential flow. This study presents an extended application of TDR for determining preferential flow properties. A TDR method was tested in carefully controlled laboratory experiments using 20-cm long and 12-cm diameter undisturbed, structured soil columns. The method used a vertically installed TDR probe and a short pulse of tracer application to obtain residual mass (RM) breakthrough curves (BTC). The RM BTC obtained from TDR were used to estimate mobile/immobile model (MIM) parameters that were compared to the parameter estimates from effluent data. A conventional inverse curve fitting method (CXTFIT) was used to estimate parameters. The TDR-determined parameters were then used to generate calculated effluent BTC for comparison with observed effluent BTC for the same soil columns. Time moments of the calculated and observed BTC were calculated to quantitatively evaluate the calculated BTC. Overall, the RM BTC obtained from TDR were similar to the RM BTC obtained from effluent data. The TDR-determined parameters corresponded well to the parameters obtained from the effluent data, although they were not within the 95% confidence intervals. Correlation coefficients between the parameters obtained from TDR and from effluent data for the immobile water fraction (theta im/theta), mass exchange coefficient (alpha), and dispersion coefficient (Dm) were 0.95, 0.95, and 0.99, respectively. For three of the four soil cores, theta im/theta ranged from 0.42 to 0.82, indicating considerable preferential flow. The TDR-calculated effluent BTC also were similar to the observed effluent BTC having an average coefficient of determination of 0.94. Time moments obtained from calculated BTC were representative of those obtained from observed BTC. The vertical TDR probe method was simple and minimally destructive and provided representative preferential flow properties that enabled the characterization of solute transport in soil.     

Simulation model for biotransformation of xenobiotics and chemotaxis in soil columns

In this paper a model is presented which can be used to simulate the behaviour of xenobiotic chemicals in soil columns with respect to their physical and chemical properties. Terms describing biological transformation of xenobiotics are also included in the model. It incorporates microbial growth following Monod kinetics and a chemotactic response of the transforming bacteria towards the xenobiotic substrate. The model appeared to yield good simulations of an experiment by Van der Meer et al. (1987) who investigated the degradation of 1,2-dichlorobenzene in soil columns inoculated with Pseudomonas sp. strain P51. The behaviour and the fate of 1,4-dichlorobenzene as found by Kuhn et al. (1985) can also be simulated using this model but their results were also adequately simulated using a simple second order model. The results generated by the model correspond to kinetic parameters obtained in other studies. It is concluded that the model is a useful tool for the investigation of the activity of bacteria degrading xenobiotics in soil columns, provided that the microbial parameters can be determined in independent experiments, and that the active microbial mass in the soil can be measured.     

Evaluating equilibrium and non-equilibrium transport of bromide and isoproturon in disturbed and undisturbed soil columns

In this study, displacement experiments of isoproturon were conducted in disturbed and undisturbed columns of a silty clay loam soil under similar rainfall intensities. Solute transport occurred under saturated conditions in the undisturbed soil and under unsaturated conditions in the sieved soil because of a greater bulk density of the compacted undisturbed soil compared to the sieved soil. The objective of this work was to determine transport characteristics of isoproturon relative to bromide tracer. Triplicate column experiments were performed with sieved (structure partially destroyed to simulate conventional tillage) and undisturbed (structure preserved) soils. Bromide experimental breakthrough curves were analyzed using convective-dispersive and dual-permeability (DP) models (HYDRUS-1D). Isoproturon breakthrough curves (BTCs) were analyzed using the DP model that considered either chemical equilibrium or non-equilibrium transport. The DP model described the bromide elution curves of the sieved soil columns well, whereas it overestimated the tailing of the bromide BTCs of the undisturbed soil columns. A higher degree of physical non-equilibrium was found in the undisturbed soil, where 56% of total water was contained in the slow-flow matrix, compared to 26% in the sieved soil. Isoproturon BTCs were best described in both sieved and undisturbed soil columns using the DP model combined with the chemical non-equilibrium. Higher degradation rates were obtained in the transport experiments than in batch studies, for both soils. This was likely caused by hysteresis in sorption of isoproturon. However, it cannot be ruled out that higher degradation rates were due, at least in part, to the adopted first-order model. Results showed that for similar rainfall intensity, physical and chemical non-equilibrium were greater in the saturated undisturbed soil than in the unsaturated sieved soil. Results also suggested faster transport of isoproturon in the undisturbed soil due to higher preferential flow and lower fraction of equilibrium sorption sites.     

Coupled production and transport of selenium vapor in unsaturated soil: evaluation by experiments and numerical simulation

Volatilization of selenium (Se) from soil to the atmosphere involves several sequential chemical reactions that form volatile Se species, followed by transport of the gaseous Se through the soil. This paper describes a numerical model that simulates the chemical and physical processes governing the production and transport of Se vapor in unsaturated soil. The model couples the four Se species involved in the production of Se vapor through chemical reactions, and allows each to migrate through the soil by advection, liquid or vapor diffusion depending on its affinity for the dissolved or vapor phase. The coupled transformations and transport of the four Se species, i.e., selenate, selenite, elemental and organic Se, and Se vapor, were calculated based on the Crank-Nicolson finite difference method. The model was used to analyze fluxes of Se vapor measured from a soil amended with inorganic Se in the form of selenate and covered with unamended clean soil of various thicknesses. Evolution of Se vapor from the soil was very fast, with measurable amounts of Se detected within 24 h. The peak of Se volatilization, detected at the 6th day, reached 3.31 Se microgram/day for the uncovered soil, but was reduced to near the detection limit (0.05 microgram/day) in the presence of a 8- or 16-cm clean soil cover. With two reaction rate coefficients fitted to the data, the model described Se volatilization very well. The estimated rate coefficient of Se methylation was unexpectedly high, with a value of 0.167/day. The net volatilization of Se, however, was severely inhibited by the fast demethylation, i.e., the reverse reaction which converted volatile Se species back into nonvolatile forms. As a result, Se vapor only penetrated a few centimeters in the soil. The demethylation rate coefficient, assessed by independent transport experiments using dimethyl selenide, was estimated as 186.8/day, corresponding to a half-life of only 5.3 min for Se vapor. Results of this study indicated that rapid demethylation of Se vapor during its diffusive transport through a soil is probably an important limiting factor in the volatilization of Se under natural conditions.     

Modeling solute transport in one-dimensional homogeneous and heterogeneous soil columns with continuous time random walk

In this paper, we used the continuous time random walk (CTRW) framework to characterize the transport process in 1250-cm long one-dimensional homogenous and heterogeneous soil columns at the experiments conducted by Huang et al. [Huang, K., Toride, N., van Genuchten, M.Th., 1995. Experimental investigation of solute transport in large, homogeneous and heterogeneous, saturated soil columns. Trans. Porous Media. 18, 283-302]. The transport process was also simulated by using the advection-dispersion equation (ADE) and the spatial fractional advection-dispersion equation (FADE) for comparison. In the homogeneous soil column, the non-Fickian behavior is found at the distances less than 1000cm with beta values larger than 1.60, but less than 2, and Fickian form transport is obtained at distances larger than 1000cm with beta values larger than 2. In the heterogeneous soil column, we found the most anomalous behavior at distances from 200cm to 700cm with beta values ranging from 0.894 to 0.958, and non-Fickian transport process is observed at distances larger than 800cm with beta values in the range between 1 and 1.3. More significant non-Fickian behavior is found for transport in the heterogeneous soil column than that in the homogeneous soil column. The CTRW fits to the breakthrough curves (BTCs) have lower values of root mean square error (RMSE) and higher values of determination coefficient (r(2)), with respect to the fits of ADE and FADE. The CTRW model also is better captures the full evolution of BTCs, and especially their tails.     

Monitoring transport of selected pesticides and phenols in soil columns by high performance liquid chromatography

Abstract

Soil columns were used to study pesticides and phenols transport under rapid infiltration land treatment conditions. An analytical procedure is described for the quantitative determination of atrazine, diuron, carbofuran, phenol, 2,4‐dinitrophenol, 2,4‐dimethylphenol, and 2,4‐dichlorophenol in soil and wastewater. Recoveries of all analytes were greater than 90%. The method detection limits for all analytes were ≤0.03 μg/ml (s/n=4) in wastewater and ≤ 0.1 μg/ml (s/n=5) in soil.     

Leaching of benzimidazole antiparasitics in soil columns and in soil columns amended with sheep excreta

Environmental Science and Pollution Research - Benzimidazoles are anthelmintics frequently used in sheep farming due to the high susceptibility of these animals to parasitic diseases. Sheep excreta...     

Phosphorus fate and transport in soil columns loaded intermittently with influent of high phosphorus concentrations.

In this study, several columns of different lengths were filled with composite soils sampled from the field at corresponding depths and then loaded intermittently with influent of a high phosphorus concentration to evaluate phosphorus fate and transport in soil. The results indicate that the height of the mass transfer zone, solvent pore velocity, and soil's life expectancy for phosphorus removal increased with depth, while the retained phosphorus per kilogram of soil and the linear adsorption equilibrium coefficient, R, decreased with depth. An equation was developed to link liquid-phase phosphorus with solvent traveling time and soil depth. The results of X-ray diffraction and washout tests indicate that calcium-phosphorus precipitation and/or crystal growth occurred in the columns. The new protocol is useful for evaluation of phosphorus fate and transport in other subsurface systems, because it allows flexible adjustments in hydraulic loadings, feed solution, and sampling schemes.     

Dual-domain solute transfer and transport processes: evaluation in batch and transport experiments

Soil macropore networks establish a dual-domain transport scenario in which water and solutes are preferentially channeled through soil macropores while slowly diffusing into and out of the bulk soil matrix. The influence of macropore networks on intra-ped solute diffusion and preferential transport in a soil typical of subsurface-drained croplands in the Midwestern United States was studied in batch- and column-scale experiments. In the batch diffusion studies with soil aggregates, the estimated diffusion radius (length) of the soil aggregates corresponded to the half-spacing of the aggregate fissures, suggesting that the intra-ped fissures reduced the diffusion impedance and preferentially allowed solutes to diffuse into the soil matrix. In the column-scale solute transport experiments, the average diffusion radius (estimated from HYDRUS-2D simulations and a first-order diffusive transfer term) was nearly double that of the batch-scale study. This increase may be attributed to a loss of pore continuity and a compounding of the small diffusion impedance through macropores at the larger scale. The column-scale solute transport experiments also suggest that two preferential networks exist in the soil. At and near soil saturation, a primary network of large macropores (possibly root channels and earthworm burrows) dominate advective transport, causing a high degree of physical and sorption nonequilibrium and simultaneous breakthrough of a nonreactive (bromide) and a reactive (alachlor) solute. As the saturation level decreases, the primary network drains, while transport through smaller macropores (possibly intra-ped features) continues, resulting in a reduced degree of nonequilibrium and separation in the breakthrough curves of bromide and alachlor.