Analyses of locally measured bromide breakthrough curves from a natural gradient tracer experiment at Krauthausen

Solute travel time distributions were derived from breakthrough curves (BTCs) of bromide concentrations, which were measured during a large-scale tracer experiment in a quaternary fluviatile aquifer at Krauthausen. Travel time distributions to a specific point in the aquifer were derived from locally measured BTCs, using averaged absolute concentrations ?abs(x1,t), normalized concentrations ?norm(x1,t), and velocity-weighted normalized concentrations ?vw(x1,t). The travel time distributions were characterized in terms of equivalent convective-dispersive transport parameters: the equivalent solute velocity and equivalent dispersivity. Parameters were derived from BTCs using moment analyses and least-squares fits of the 1-D convection-dispersion equation (CDE). Both local and averaged BTCs showed pronounced tailing which was not well described by the 1-D CDE and which indicates the presence of macroscopic regions with low velocities in the aquifer. Therefore, dispersivities derived from CDE fits were significantly smaller than those derived from time moments. The BTCs of ?abs(x1,t) were dominated by only a few local BTCs with high concentrations and were less representative for the travel time distribution than BTCs of averaged normalized concentrations. Dispersivities derived from ?norm(x1,t) and ?vw(x1,t) were very similar. Finally, estimates of dispersivities and vertical correlation length of lnK, gamma 3, from BTCs were in agreement with a first-order estimate of the dispersivity and gamma 3 based on grain size data and flow meter measurements.     

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.     

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.     

Characterizing pesticide sorption and degradation in microscale biopurification systems using column displacement experiments

Biopurification systems treating pesticide contaminated water are very efficient, however they operate as a black box. Processes inside the system are not yet characterized. To optimize the performance, knowledge of degradation and retention processes needs to be generated. Therefore, displacement experiments were carried out for four pesticides (isoproturon, bentazone, metalaxyl, linuron) in columns containing different organic mixtures. Bromide, isoproturon and bentazone breakthrough curves (BTCs) were well described using the convection-dispersion equation (CDE) and a first-order degradation kinetic approach. Metalaxyl and linuron BTCs were well described using the CDE model expanded with Monod-type kinetics. Freundlich sorption, first-order degradation and Monod kinetics coefficients were fitted to the BTCs. Fitted values of the distribution coefficient K_f,column were much lower than those determined from batch experiments. Based on mobility, pesticides were ranked as: bentazone > metalaxyl - isoproturon > linuron. Based on degradability, pesticides were ranked as: linuron > metalaxyl - isoproturon > bentazone.     

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.     

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.     

Use of tracer tests to investigate changes in flow and transport properties due to bioclogging of porous media

Tracer tests were conducted in three laboratory columns to study changes in the hydraulic properties of a porous medium due to bioclogging. About 30 breakthrough curves (BTCs) for each column were obtained. The BTCs were analyzed using analytical equilibrium and dual-porosity models, and estimates of the hydrodynamic dispersion and mass transfer coefficients were obtained by curve fitting. The change in transport properties developed in three stages: an initial phase (I) with no significant changes in transport properties, phase II with growth of biomass near the inlet of the columns causing changes in dispersivity, and phase III with added growth of micro-colonies deeper in the columns causing mass transfer of solutes from the water phase to the biophase. Tracer transport changed from being uniform to more non-uniform with increase in mass transfer of the tracer between the mobile phase and the immobile biomass. An increase in the bulk dispersivity value of up to one order of magnitude was observed. Numerical simulations suggest that local dispersivity values may be as much as 40 times higher in the more severe clogged areas inside the column. The bulk hydraulic conductivities of the columns decreased by up to three orders of magnitude. The hydraulic conductivity and dispersivity parameters were almost recovered after disinfection of the columns. Different models relating the changes of the hydraulic conductivity to the changes in the mobile porosity due to bioclogging were reviewed, and the micro-colony relation of Thullner et al. [Thullner, M., Zeyer, J., Kinzelbach, W., 2002. Influence of microbial growth on hydraulic properties of pore networks, Transport in Porous Media, 49, 99-122.] was found to best describe the relation between the bulk hydraulic parameters.     

Second-order modeling of arsenite transport in soils

Rate limited processes including kinetic adsorption-desorption can greatly impact the fate and behavior of toxic arsenic compounds in heterogeneous soils. In this study, miscible displacement column experiments were carried out to investigate the extent of reactivity during transport of arsenite in soils. Arsenite breakthrough curves (BTCs) of Olivier and Windsor soils exhibited strong retardation with diffusive effluent fronts followed by slow release or tailing during leaching. Such behavior is indicative of the dominance of kinetic retention reactions for arsenite transport in the soil columns. Sharp decrease or increase in arsenite concentration in response to flow interruptions (stop-flow) further verified that non-equilibrium conditions are dominant. After some 40-60 pore volumes of continued leaching, 30-70% of the applied arsenite was retained by the soil in the columns. Furthermore, continued arsenite slow release for months was evident by the high levels of residual arsenite concentrations observed during leaching. In contrast, arsenite transport in a reference sand material exhibited no retention where complete mass recovery in the effluent solution was attained. A second-order model (SOM) which accounts for equilibrium, reversible, and irreversible retention mechanisms was utilized to describe arsenite transport results from the soil columns. Based on inverse and predictive modeling results, the SOM model successfully depicted arsenite BTCs from several soil columns. Based on inverse and predictive modeling results, a second-order model which accounts for kinetic reversible and irreversible reactions is recommended for describing arsenite transport in soils.     

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.     

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.     

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.     

Effects of degree of water saturation on dispersivity and immobile water in sandy soil columns

Three natural nonaggregated soil samples, with similar grain-size distributions, have been used to determine the dispersive behavior of porous media under steady, saturated and unsaturated flow conditions. Tritium was used as a tracer and was found to have no sorption on the solid matrix. Generated breakthrough curves (BTCs) for the unsaturated experiments were symmetrical with no evidence of tailing. The unsaturated experiments for two of the soils were adequately described by considering all the water in the pore volume as mobile. However, about 10% of the pore water, independent of the degree of saturation, was found to be immobile in the case of the third soil during unsaturated flow. For this soil, there was no mass transfer between the two water regions, indicating that the immobile water is essentially isolated from the flowing water fraction. For all three soils, dispersivity under unsaturated conditions was found to be higher, independent of the degree of water saturation, than the value determined for the saturated experiments. This is inconsistent with what would be expected from the simple bundle-of-capillary-tubes model and does not agree well with a more sophisticated conceptualization of the porous medium. The data, however, clearly indicate a wider range in pore-water velocities when these soils are desaturated.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

Modelling field-scale cadmium transport below the root zone of a sewage sludge amended soil in an arid region in Central Iran

Addition of trace metals such as cadmium to soils in metal-rich sewage sludge may result in contamination of soil and groundwater. This study addresses the plot-scale transport of Cd derived from sewage sludge in a layered clay soil in an arid region of central Iran. Sewage sludge was enriched by Cd at rates of 38 and 80 mg kg(-1) and applied to experimental soil plots using a complete random block design with three replicates. Cadmium concentration was measured as a function of depth after 185 and 617 days. HYDRUS-1D and MACRO codes were calibrated for Cd transport in the site treated with 80 mg kg(-1) sewage sludge. Model parameters were estimated by inverse modelling using the SUFI-2 procedure. The site treated with 38 mg kg(-1) cadmium was used to test the calibrated models. Both convection-dispersion equation (CDE) and non-equilibrium CDE in HYDRUS-1D produced reasonable calibration results. However, the estimated Freundlich sorption constants were significantly smaller than those measured in a batch study. A site tracer experiment revealed the existence of substantial macropore flow. For this reason we applied MACRO to account for this process. The calibration and test results with MACRO were as good as those obtained by HYDRUS-1D with the difference that adsorption constants were much closer to the measured ones. This indicates that in HYDRUS-1D, the adsorption parameters were underestimated in order to allow a deeper transport of Cd which had actually occurred due to macropore flow. A 20-year simulation scenario depicting the long-term effect of sludge application indicated small risk of groundwater contamination. However, high concentration of Cd near the soil surface raises a concern about the crop Cd uptake which should be further investigated.     

Transport and reduction of nitrate in clayey till underneath forest and arable land

Transport and reduction of nitrate in a typically macroporous clayey till were examined at variable flow rate and nitrate flux. The experiments were carried out using saturated, large diameter (0.5 m), undisturbed soil columns (LUC), from a forest and nearby agricultural sites. Transport of nitrate was controlled by flow along the macropores (fractures and biopores) in the columns. Nitrate reduction (denitrification) determined under active flow mainly followed first order reactions with half-lives (t(1/2)) increasing with depth (1.5-3.5 m) from 7 to 35 days at the forest site and 1-7 h at the agricultural site. Nitrate reduction was likely due to microbial degradation of accumulated organic matter coupled with successive consumption of O2 and NO3- in the macropore water followed by reductive dissolution of Fe and Mn from minerals along the macropores. Concentrations of total organic carbon measured in soil samples were near identical at the two study sites and consequently not useful as indicator for the observed differences in nitrate reduction. Instead the high reduction rates at the agricultural site were positively correlated with elevated concentration of water-soluble organic carbon and nitrate-removing bacteria relative to the forest site. After high concentrations of water-soluble organic carbon in the columns from the agricultural site were leached they lost their elevated reduction rates, which, however, was successfully re-established by infiltration of new reactive organics represented by pesticides. Simulations using a calibrated discrete fracture matrix diffusion (DFMD) model could reasonably reproduce the denitrification and resulting flux of nitrate observed during variable flow rate from the columns.     

Effect of seepage conditions on chemical attenuation of arsenic by soils across an abandoned mine site

The effect of seepage velocity on the As leaching/adsorption by soils collected from abandoned mine sites was evaluated under batch equilibrium and different seepage settings. The breakthrough curves (BTCs) of As leaching from the mine soil column initially displayed the peak export and gradually leveled off over the leaching experiment. A similar As peak was observed after a flow interruption period. Adsorption by downgradient soils was clearly nonlinear, as Freundlich N was <1. In the BTCs of the layered columns, where downgradient soils were overloaded above the mine soil, the extended lag period of As appearance and lower steady-state As concentration observed for slow seepage velocity supported the idea of kinetically limited As attenuation driven by soil adsorption. The perturbation of As concentration was insignificant when the intra-column As concentration gradient was higher. The As concentration drop and time to recovery were greater for less adsorptive soil and fast seepage velocity. Desorption of As from soils retrieved from both batch adsorption and column experiment demonstrate hysteric behavior. The results of this work demonstrated that the risk of As leaching from an abandoned mine site can be greatly attenuated by intermediate downgradient soils via chemical adsorption, which tends to be kinetically limited and energetically hysteric (i.e., non-identical energy pathway).     

Using aliphatic alcohols as gaseous tracers in determination of water contents and air-water interfacial areas in unsaturated sands

A new type of gaseous tracer utilizing nontoxic aliphatic alcohols for the determination of water content and air-water interfacial area is tested on unsaturated sands of low water content. Alcohol vapors are generated at room temperature and passed through the experimental sand column. Breakthrough curves (BTCs) of these vapors are obtained by monitoring their effluent concentrations using GC-FID. The retardation factor with respect to each vapor transport process is obtained by optimizing BTCs data using the CXTFIT program in the reverse problem mode. The water content and the interfacial area are subsequently calculated from their retardation factors by both equilibrium and nonequilibrium transport models. Experimental results indicate that the pentanol tracer is feasible in the determination of water content at conditions when the degree of water saturation is low. In the determination of air-water interfacial area, decanol is selected due to its interfacial adsorption characteristics. By comparing to interfacial areas from theoretical predictions as well as other conventional tarcer methods, the ones determined from the decanol tracer tests are found to be close to the true interfacial areas when the water content is low.     

The influence of organic ligands on the retention of lead in soil

Organic acids are commonly produced and exuded by plant roots and soil microorganisms. Some of these organic compounds are effective chelating agents and have the potential to enhance metal mobility. The effect of citrate and salicylate on the leaching of lead in soil was investigated in a laboratory experiment. In short-term batch experiments, adsorption of lead to soil was slightly enhanced with increasing salicylate concentration (500-5000 microM) but decreased significantly in the presence of citrate. These observations suggested that citrate may enhance Pb leaching, but this was not observed in the column study. Soluble Pb in the presence and absence citrate or salicylate (up to 5000 microM) was added to soil columns at a moderate flow rate, but no Pb was observed to emerge from the soil in any of the soil columns. Rapid biodegradation of citrate in soil eliminated potential complexing ability. Breakthrough of Pb from soil was noted only when using small columns at high flow rates (>20 pore volumes per day). Under these conditions of physical and chemical non-equilibrium, citrate was not degraded and significantly enhanced Pb mobility. As in the batch adsorption experiments, the presence of salicylate reduced Pb leaching. Considering the extreme conditions required to induce Pb leaching, it is likely that Pb will remain relatively immobile in soil even in the presence of a strong complexing agent such as citrate.