Iron and organo-bentonite for the reduction and sorption of trichloroethylene

Hybrid barriers using dechlorination and immobilization were studied to remove trichloroethylene (TCE) in this study. Hybrid barriers of iron filings and organo (hexadecyltrimethylammonium, HDTMA)-bentonite were simulated in columns to assess the performance of the hybrid barriers. TCE reduction rate for the mixture of zero valent iron (ZVI) and HDTMA-bentonite was approximately seven times higher than that for ZVI, suggesting the reduction of TCE was accelerated when HDTMA-bentonite was mixed with ZVI. For the column of two separate layers of iron and HDTMA-bentonite, TCE reduction rate was nearly similar to that for ZVI alone, but the partition coefficient (Kd) was 4.5 times higher than that for ZVI only. TCE was immobilized in the first layer with HDTMA-bentonite due to sorption, and then dechlorinated in the second layer with iron filings due to reduction. The HDTMA-bentonite and minimally-desorbed HDTMA from the organo-bentonite are believed to contribute the increase in TCE concentration on iron surface so that more TCE could be available for reduction. Therefore, the incorporation of HDTMA-bentonite into ZVI not only can effectively retard the transport of chlorinated organic contaminants from landfill leachate or oil shock in subsurface environment, also can expedite the reduction rate of TCE.     

Prediction of mass-transfer coefficient for solute transport in porous media

Several previously reported laboratory studies related to transport of solutes through packed columns were utilized to develop predictive relationships for mass-transfer rate coefficient. The data were classified into two groups: those obtained under rate-limited mass transfer between mobile and immobile water regions (physical nonequilibrium conditions), and those derived from rate-limited mass transfer between instantaneous and slow sorption sites (sorption nonequilibrium conditions). The mass-transfer coefficient in all these studies was obtained by fitting breakthrough data to a transport model employing a first-order rate limitations with a "constant" mass-transfer coefficient, independent of flow conditions. This study demonstrated that the mass-transfer coefficient in these models is dependent on system parameters including pore-water velocity, length-scale, retardation coefficient, and particle or aggregate size. Predictive relationships were developed, through regression analysis, relating mass-transfer coefficient to residence time. The developed relationships adequately estimated previously reported field mass-transfer values. Successful simulations of field desorption data reported by Bahr [J. Contam. Hydrol. 4 (1989) 205] further demonstrate the potential applicability of the developed relationships.     

Effect of cell physicochemical characteristics and motility on bacterial transport in groundwater

The influence of physicochemical characteristics and motility on bacterial transport in groundwater were examined in flow-through columns. Four strains of bacteria isolated from a crystalline rock groundwater system were investigated, with carboxylate-modified and amidine-modified latex microspheres and bromide as reference tracers. The bacterial isolates included a gram-positive rod (ML1), a gram-negative motile rod (ML2), a nonmotile mutant of ML2 (ML2m), and a gram-positive coccoid (ML3). Experiments were repeated at two flow velocities, in a glass column packed with glass beads, and in another packed with iron-oxyhydroxide coated glass beads. Bacteria breakthrough curves were interpreted using a transport equation that incorporates a sorption model from microscopic observation of bacterial deposition in flow-cell experiments. The model predicts that bacterial desorption rate will decrease exponentially with the amount of time the cell is attached to the solid surface. Desorption kinetics appeared to influence transport at the lower flow rate, but were not discernable at the higher flow rate. Iron-oxyhydroxide coatings had a lower-than-expected effect on bacterial breakthrough and no effect on the microsphere recovery in the column experiments. Cell wall type and shape also had minor effects on breakthrough. Motility tended to increase the adsorption rate, and decrease the desorption rate. The transport model predicts that at field scale, desorption rate kinetics may be important to the prediction of bacteria transport rates.     

In situ modification of an aquifer material by a cationic surfactant to enhance retardation of organic contaminants

Sorption of hexadecyltrimethylammonium chloride (HDTMA), a cationic surfactant, on aquifer material from Columbus AFB, Mississippi, U.S.A., was examined. Transport studies using flow-through columns and a box model aquifer showed that an almost stationary zone of HDTMA-modified aquifer material could be produced in situ without a significant decrease in hydraulic conductivity.Perchloroethylene (PCE) and naphthalene sorption isotherms on the HDTMA-modified aquifer material were linear, and sorption coefficients were increased by over two orders of magnitude relative to the unmodified material. The retardation of PCE by insitu emplaced HDTMA zones within a column was examined. Agreement between batch- and column-derived sorption coefficients and breakthrough curve symmetry indicates that local equilibrium was attained. Significant retardation of a naphthalene plume by an in situ emplaced surfactant zone was demonstrated in the box model aquifer system.The experimental results indicate that it is feasible to create in situ a sorbent zone within an aquifer using cationic surfactants. In most situations, the sorbent zone concept needs to be coupled with contaminant degradation processes for sorbent emplacement to be a practical tool in the remediation of groundwater contamination sites. Sorbent zones may be of benefit in the engineering of suitable environments for microbial or abiotic degradation reactions and by providing time slow reactions to occur.     

Sorption and reduction of tetrachloroethylene with zero valent iron and amphiphilic molecules

Effects of surfactants and natural organic matter (NOM) on the sorption and reduction of tetrachloroethylene (PCE) with zero valent iron (ZVI) were examined in this study. PCE reduction by ZVI depended on the ionic type of the surfactants. The removal of PCE and production of TCE with non-ionic Triton X-100 and cationic hexadecyltrimethyl-ammonium (HDTMA) at one-half and two times the critical micelle concentration (CMC) were 1.2-1.8 times higher than without surfactants because of the enhanced PCE partitioning and surface concentration by the sorbed surfactants. When anionic sodium dodecyl benzene sulfonate (SDDBS) at one-half and two times CMC and NOM at 20 mg l(-1) and 50 mg l(-1) concentrations were used, the removal of PCE doubled and TCE production decreased. In the presence of SDDBS, TCE production by ZVI was lower than with HDTMA and Triton X-100 while PCE removal was higher than with the other surfactants.     

Simultaneous removal of NO and SO2 by high-temperature fluidized zero-valent iron processes

A new approach to simultaneously remove nitrogen monoxide (NO) and sulfur dioxide (SO2) by zero valent iron (ZVI) was investigated. Three different parameters, temperature, flux, and ZVI dosage, were tested in fluidized ZVI column studies containing 500 ppmv of NO and SO2, respectively. Under the ZVI dosage of 0.5 g at flux of 0.6 L/cm2 x min for temperature 573 K, there is neither NO nor SO2 reduction. For 623 K and 673 K, complete removal for NO and > 90% removal for SO2 were achieved. For temperatures of 723 K and 773 K, 100% removal was achieved for both NO and SO2. The amounts of NO or SO2 reduction (as milligrams of NO or SO2 per gram ZVI) increased as temperature increased, and linearities were observed with both correlation coefficients > 0.97. Compared with NO, SO2 had earlier breakthrough because of a slower diffusion rate and less reactivity but higher mass reduction because of a higher molecular weight for SO2 (64 g/mol for SO2 and 30 g/mol for NO). At same temperature, both NO and SO2 reductions (as milligrams of NO or SO2 per gram of ZVI) were constant regardless of either flux or ZVI dosage variation, but breakthrough time was affected by both flux and ZVI dosage. A parameter weight of ZVI/flux (W/F) was developed to represent these two parameters at the same time to assess the breakthrough time of NO and SO2. Higher breakthrough time was achieved for higher W/F value. Moreover, interestingly, longer breakthrough time and more NO and SO2 mass reduction were achieved for combined NO and SO2 than individual NO or SO2 treated by ZVI, and both oxidation and reduction reactions occurred instead of a reduction reaction only. Chemical reactions among ZVI/NO, ZVI/ SO2, and ZVI/NO/SO2 were also proposed and verified by X-ray diffraction analyses.     

Strontium migration in a crystalline medium: effects of the presence of bentonite colloids

The effects of bentonite colloids on strontium migration in fractured crystalline medium were investigated. We analyzed first the transport behaviour of bentonite colloids alone at different flow rates; then we compared the transport behaviour of strontium as solute and of strontium previously adsorbed onto stable bentonite colloids at a water velocity of approximately 7.1·10(-6)m/s-224m/yr. Experiments with bentonite colloids alone showed that - at the lowest water flow rate used in our experiments (7.1·10(-6)m/s) - approximately 70% of the initially injected colloids were retained in the fracture. Nevertheless, the mobile colloidal fraction, moved through the fracture without retardation, at any flow rate. Bentonite colloids deposited over the fracture surface were identified during post-mortem analyses. The breakthrough curve of strontium as a solute, presented a retardation factor, R(f)~6, in agreement with its sorption onto the granite fracture surface. The breakthrough curve of strontium in the presence of bentonite colloids was much more complex, suggesting additional contributions of colloids to strontium transport. A very small fraction of strontium adsorbed on mobile colloids moved un-retarded (R(f)=1) and this fraction was much lower than the expected, considering the quantity of strontium initially adsorbed onto colloids (90%). This behaviour suggests the hypothesis of strontium sorption reversibility from colloids. On the other hand, bentonite colloids retained within the granite fracture played a major role, contributing to a slower strontium transport in comparison with strontium as a solute. This was shown by a clear peak in the breakthrough curve corresponding to a retardation factor of approximately 20.     

Oxyanion flux characterization using passive flux meters: development and field testing of surfactant-modified granular activated carbon

We report here on the extension of Passive Flux Meter (PFM) applications for measuring fluxes of oxyanions in groundwater, and present results for laboratory and field studies. Granular activated carbon, with and without impregnated silver (GAC and SI-GAC, respectively), was modified with a cationic surfactant, hexadecyltrimethylammonium (HDTMA), to enhance the anion exchange capacity (AEC). Langmuir isotherm sorption maxima for oxyanions measured in batch experiments were in the following order: perchlorate>chromate>selenate, consistent with their selectivity. Linear sorption isotherms for several alcohols suggest that surfactant modification of GAC and SI-GAC reduced (approximately 30-45%) sorption of alcohols by GAC. Water and oxyanion fluxes (perchlorate and chromate) measured by deploying PFMs packed with surfactant-modified GAC (SM-GAC) or surfactant-modified, silver-impregnated GAC (SM-SI-GAC) in laboratory flow chambers were in close agreement with the imposed fluxes. The use of SM-SI-GAC as a PFM sorbent was evaluated at a field site with perchlorate contamination of a shallow unconfined aquifer. PFMs packed with SM-SI-GAC were deployed in three existing monitoring wells with a perchlorate concentration range of approximately 2.5 to 190 mg/L. PFM-measured, depth-averaged, groundwater fluxes ranged from 1.8 to 7.6 cm/day, while depth-averaged perchlorate fluxes varied from 0.22 to 1.7 g/m2/day. Groundwater and perchlorate flux distributions measured in two PFM deployments closely matched each other. Depth-averaged Darcy fluxes measured with PFMs were in line with an estimate from a borehole dilution test, but much smaller than those based on hydraulic conductivity and head gradients; this is likely due to flow divergence caused by well-screen clogging. Flux-averaged perchlorate concentrations measured with PFM deployments matched concentrations in groundwater samples taken from one well, but not in two other wells, pointing to the need for additional field testing. Use of the surfactant-modified GACs for measuring fluxes of other anions of environmental interest is discussed.     

Competitive sorption and diffusion of chromate and sulphate in a flow system with goethite in gel beads

Column experiments and model simulations were employed to evaluate the processes involved in multicomponent solute transport in a system with heterogeneous flow. Column experiments were performed with goethite embedded in polyacrylamide gel beads. The gel forms an immobile water region that can be accessed by diffusion. A two-region transport model with diffusion into spheres was combined with a surface complexation model to predict reactive transport in the goethite-gel bead system. Chromate and sulphate breakthrough curves were measured in a set of transport experiments, along with corresponding changes in the pH of the effluent. Sorption and transport of sulphate and chromate in separate columns were predicted from independently measured sorption parameters. The model overestimated the pH changes in the effluent, possibly because of proton buffering by the polyacrylamide gel. The effect of competitive sorption on transport was examined in experiments with both anions present. The model predicted the effect of competition very well in a system initially equilibrated with sulphate, followed by infiltration with chromate. However, when sulphate was infiltrated after equilibration with chromate, chromate desorption and sulphate adsorption were clearly overestimated by the transport model. The exchange between the more strongly bound chromate and the sulphate added subsequently may be too slow to cause a substantial chromate peak in the effluent. This suggests that the local equilibrium assumption was not applicable in this case.     

Application of the method of temporal moments to interpret solute transport with sorption and degradation

In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data.     

Influence of modified soils on the removal of diesel fuel oil from water and the growth of oil degradation micro-organism

Three soils were modified with two kinds of cationic surfactants in order to increase their sorptive capabilities for organic contaminants. Sorption of diesel fuel oil in water by these modified soils had been investigated. Modified soils can effectively sorb diesel fuel oil from water. The sorption capability of modified soils is: HDTMA-black soil > HDTMA-yellow brown soil > HDTMA-red soil > TMA-black soil > TMA-yellow brown soil > TMA-red soil. Sorption of diesel fuel oil by natural soils and HDTMA modified soils is via partition, the sorption isotherms can be expressed by Henry equation, and logK(SOM) is 2.42-2.80, logK(HDTMA) is 3.37-3.60. Sorption isotherms of TMA modified soils can be expressed by Langmuir equation, the saturation sorption capacities are 1150 (TMA-black soil), 750 (TMA-yellow-brown soil), 171 mg/kg (TMA-red soil), respectively. A diesel fuel oil degradation micro-organism (Pseudomonas sp.) was isolated in the lab. To test the influence of the modified soils on the micro-organism, various growth curves of Pseudomonas in different conditions were drawn. Pseudomonas can grow very well with natural soils and TMA modified soils. The acclimation period of Pseudomonas is reduced. As to HDTMA modified soils, HDTMA loading amount is very important. When HDTMA loading amount is no higher than 0.5 CEC, the micro-organism can grow very well after a long acclimation period.     

Simultaneous sorption of lead and chlorobenzene by organobentonite

Clays or organoclays have been used as a barrier to prevent the transport of hazardous contaminants in landfills. However, clays are known to effectively sorb mostly inorganic contaminants, while organoclays are mainly used for organic contaminants. Since the organoclays are basically clay particles modified with cationic surfactants, there might exist an optimal coverage of cationic surfactant on the clay particles to sorb both inorganic and organic contaminants. In order to determine the optimal mass of cationic surfactants on the bentonites, sodium bentonites were treated with various ratios of hexadecyltrimethylammonium (HDTMA) to bentonites. Chlorobenzene and lead were selected as representative contaminants. When either chlorobenzene or lead exists as a single contaminant, chlorobenzene sorption increased with increasing HDTMA to bentonite ratios, and lead sorption decreased with increasing HDTMA to bentonite ratios. Sorption of chlorobenzene was a function of HDTMA coverage on the bentonites, while lead sorption was much more influenced by the initial lead concentration rather than the mass of HDTMA added to the bentonites.     

Movement of metolachlor and terbuthylazine in core and packed soil columns

Movement of metolachlor and terbuthylazine including a bromide tracer was studied in core and packed soil columns in PVC pipes (80 mm diameter, 15 mm depth) with two German soil types viz: silt loam and loamy silt. The breakthrough curves (BTCs) for bromide indicated some preferential flow of water both under conventional tillage (CN) and no-till (NT) simulation with silt loam soil. The herbicides leached to a greater extent in NT columns than in CN columns. Leaching was higher in loamy silt soil than in silt loam soil under CN conditions. This result is in agreement with the higher sorption capacity of silt loam having higher organic carbon compared to loamy silt having low organic carbon. Adsorption strength of the herbicides did not affect their breakthrough time, but was reflected in the slope and maximum height of the BTCs. The BTCs showed the expected inverse relationship between leaching and adsorption with greater mobility of the weakly-sorbed metolachlor than the more strongly sorbed terbuthylazine. Maximum amounts of the applied herbicides were recovered from the top soil layer in intact columns. Metolachlor was more mobile in packed columns than in core columns.     

Transport and transformation of sulfadiazine in soil columns packed with a silty loam and a loamy sand

Concerning the transport of the veterinary antibiotic sulfadiazine (SDZ) little is known about its possible degradation during transport. Also its sorption behaviour is not yet completely understood. We investigated the transport of SDZ in soil columns with a special emphasis on the detection of transformation products in the outflow of the soil columns and on modelling of the concentration distribution in the soil columns afterwards. We used disturbed soil columns near saturation, packed with a loamy sand and a silty loam. SDZ was applied as a 0.57 mg L(-1) solution at a constant flow rate of 0.25 cm h(-1) for 68 h. Breakthrough curves (BTC) of SDZ and its transformation products 4-(2-iminopyrimidin-1(2H)-yl)aniline and 4-hydroxy-SDZ were measured for both soils. For the silty loam we additionally measured a BTC for an unknown transformation product which we only detected in the outflow samples of this soil. After the leaching experiments the (14)C-concentration was quantified in different layers of the soil columns. The transformation rates were low with mean SDZ mass fractions in the outflow samples of 95% for the loamy sand compared to 97% for the silty loam. The formation of 4-(2-iminopyrimidin-1(2H)-yl)aniline appears to be light dependent and did probably not occur in the soils, but afterwards. In the soil columns most of the (14)C was found near the soil surface. The BTCs in both soils were described well by a model with one reversible (kinetic) and one irreversible sorption site. Sorption kinetics played a more prominent role than sorption capacity. The prediction of the (14)C -concentration profiles was improved by applying two empirical models other than first order to predict irreversible sorption, but also these models were not able to describe the (14)C concentration profiles correctly. Irreversible sorption of sulfadiazine still is not well understood.     

Determination of bromacil transport as a function of water and carbon content in soils

This study was conducted to determine the significance of bromacil transport as a function of water and carbon content in soils and to explore the implications of neglecting sorption when making assessments of travel time of bromacil through the vadose zone. Equilibrium batch sorption tests were performed for loamy sand and sandy soil added with four different levels of powdered activated carbon (PAC) content (0, 0.01, 0.05, and 0.1%). Column experiments were also conducted at various water and carbon contents under steady-state flow conditions. The first set of column experiments was conducted in loamy sand containing 1.5% organic carbon under three different water contents (0.23, 0.32, and 0.41) to measure breakthrough curves (BTCs) of bromide and bromacil injected as a square pulse. In the second set of column experiments, BTCs of bromide and bromacil injected as a front were measured in saturated sandy columns at the four different PAC levels given above. Column breakthrough data were analyzed with both equilibrium and nonequilibrium (two-site) convection-dispersion equation (CDE) models to determine transport and sorption parameters under various water and carbon contents. Analysis with batch data indicated that neglect of the partition-related term in the calculation of solute velocity may lead to erroneous estimation of travel time of bromacil, i.e. an overestimation of the solute velocity by a factor of R. The column experiments showed that arrival time of the bromacil peak was larger than that of the bromide peak in soils, indicating that transport of bromacil was retarded relative to bromide in the observed conditions. Extent of bromacil retardation (R) increased with decreasing water content and increasing PAC content, supporting the importance of retardation in the estimation of travel time of bromacil even at small amounts of organic carbon for soils with lower water content.     

AMMONIA REMOVAL FROM COMPOST LACHATE USING ZEOLITE. II. A STUDY USING CONTINUOUS FLOW PACKED COLUMNS

Bench-scale packed zeolite columns were set up and operated to investigate the continuous removal of ammonium ions from compost leachate. The effects of hydraulic retention time (HRT), and particle size of the zeolite on the ammonia adsorption capacity were studied. For both the coarse particle and the powdered zeolite columns, higher ammonia removal efficiencies were achieved with longer HRT (i.e., lower influent flow rate) tests. At the same HRT, ammonia removal efficiencies from tests with powdered zeolite were generally 20% higher than tests with the coarse particle zeolite. A HRT of 6 hours was found appropriate for efficient ammonia removal, and an operating capacity of 1.31 mg N/g zeolite was obtained. Over 98% of the ammonia input from the influent was consistently removed for over 5 bed volumes (BV) of compost leachate flowing through the zeolite column. Zeolite proved to have a great potential as a medium for ammonia removal in treating composting leachate.     

Ammonia removal from compost leachate using zeolite. II. A study using continuous flow packed columns

Bench-scale packed zeolite columns were set up and operated to investigate the continuous removal of ammonium ions from compost leachate. The effects of hydraulic retention time (HRT), and particle size of the zeolite on the ammonia adsorption capacity were studied. For both the coarse particle and the powdered zeolite columns, higher ammonia removal efficiencies were achieved with longer HRT (i.e., lower influent flow rate) tests. At the same HRT, ammonia removal efficiencies from tests with powdered zeolite were generally 20% higher than tests with the coarse particle zeolite. A HRT of 6 hours was found appropriate for efficient ammonia removal, and an operating capacity of 1.31 mg N/g zeolite was obtained. Over 98% of the ammonia input from the influent was consistently removed for over 5 bed volumes (BV) of compost leachate flowing through the zeolite column. Zeolite proved to have a great potential as a medium for ammonia removal in treating composting leachate.     

Effect of pore-water velocity on chemical nonequilibrium transport of Cd,Zn, and Pb in alluvial gravel columns

This paper investigates the effects of pore-water velocity on chemical nonequilibrium during transport of Cd, Zn, and Pb through alluvial gravel columns. Three pore-water velocities ranging from 3 to 60 m/day were applied to triplicate columns for each metal. Model results for the symmetric breakthrough curves (BTCs) of tritium (3H2O) data suggest that physical nonequilibrium components were absent in the uniformly packed columns used in these studies. As a result, values of pore-water velocity and dispersion coefficient were estimated from fitting 3H2O BTCs to an equilibrium model. The BTCs of metals display long tailing, indicating presence of chemical nonequilibrium in the system, which was further supported by the decreased metal concentrations during flow interruption. The BTCs of the metals were analysed using a two-site model, and transport parameters were derived using the CXTFIT curve-fitting program. The model results indicate that the partitioning coefficient (beta), forward rate (k1), and backward rate (k2) are positively correlated with pore-water velocity (V); while the retardation factor (R), mass transfer coefficient ((omega), and ratio of k1/k2 are inversely correlated with V. There is no apparent relationship between the fraction of exchange sites at equilibrium (f) and V. The influence of Von k2 is much greater than on R, beta, omega, and k1. A one-order-of-magnitude change in V would cause a two-order-of-magnitude change in k2 while resulting in only a one order-of-magnitude change in R, beta, omega, and k1. The forward rates for the metals are found to be two to three orders-of-magnitude greater than the corresponding backward rate. However, the difference between the two rates reduces with increasing pore-water velocity. Model results also suggest that Cd and Zn behave similarly, while Pb is much more strongly sorbed. At input concentrations of about 4 mg/l and pore-water velocities of 3-60 m/day in the groundwater within alluvial gravel, this study suggests retardation factors of 26-289 for Cd, 24-255 for Zn, and 322-6377 for Pb.     

Continuous time random walk model better describes the tailing of atrazine transport in soil

Contaminant transport in soils is complicated and involves some physical and chemical nonequilibrium processes. In this research, the soil column displacement experiments of Cl⁻ and atrazine under different flow velocities were carried out. The data sets of Cl⁻ transport in sandy loam fitted to the convection dispersion equation (CDE) and the two-region model (TRM) indicated that the effects of physical nonequilibrium process produced by immobile water on the breakthrough curves (BTCs) of Cl⁻ and atrazine transport through the repacking soil columns were negligible. The two-site model (TSM) and the continuous time random walk (CTRW) were also used to fit atrazine transport behavior at the flow rate of 19.86 cm h⁻¹. The CTRW convincingly captured the full evolution of atrazine BTC in the soil column, especially for the part of long tailing. However, the TSM failed to characterize the tailing of atrazine BTC in the soil column. The calculated fraction of equilibrium sorption sites, F, ranging from 0.78 to 0.80 for all flow rates suggested the contribution of nonequilibrium sorption sites to the asymmetry of atrazine BTCs. Furthermore, the data sets for the flow rates of 6.68 cm h⁻¹ and 32.81 cm h⁻¹ were predicted by the TSM and the CTRW. As to the flow rate of 6.68 cm h⁻¹, the CTRW predicted the entire BTC of atrazine transport better than the TSM did. For the flow rate of 32.81 cm h⁻¹, the CTRW characterized the late part of the tail better, while the TSM failed to predict the tailings of atrazine BTC.