Modelling the solute transport under nonequilibrium conditions on the basis of mass transfer equations

A solute transport model that describes nonequilibrium adsorption in soil/groundwater systems by mass transfer equations for film and intraparticle diffusion is presented. The model is useful in cases where breakthrough curve spreading cannot be explained by dispersion only. To evaluate its validity, the model was applied to several data sets from column experiments. The validity was also proved by a comparison with an analytical solution for the limiting case of predominating dispersion. Furthermore, a sensitivity analysis was performed to illustrate the influence of different process and sorption parameters (pore water velocity, intraparticle mass transfer coefficient, isotherm nonlinearity) on the shape of the calculated breakthrough curves. The application of the proposed model is discussed in comparison to the widely used dispersed flow/local equilibrium model, and a relationship between both models, which is based on a lumped parameter approach, is shown.     

Modelling of atrazine transport in the presence of surfactants

Laboratory experiments were conducted to examine the effect of detergents on transport of atrazine [2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine] through loam and sandy loam soils under saturation conditions. The Convection Dispersion Equation (CDE) was used to model and quantify the effects of detergents on atrazine model parameters: the retardation factor (R), pore velocity (v) and dispersion coefficient (D). The transport parameters were estimated using moment technique and partition coefficient obtained from batch experiments and compared with best-fitted parameters, R and D, keeping pore velocity constant. Results indicated the CDE model was not successful in predicting atrazine transport in the presence of surfactants at high concentrations. In the case of anionic surfactant with Elora loam, the average predicted R and D from moment technique of 3.4 and 11.1 cm2/h, respectively were significantly different than fitted parameters (R = 39 and D = 227 cm2/h). The poor performance of CDE in the presence of surfactants results from physiochemical changes in herbicide solubility and retention to the soil matrix rather than changes in soil hydraulic properties since the predicted pore water velocities from moment technique were similar to those measured during leaching experiments. Nevertheless, BTC analysis with CDE showed that land application of anionic surfactant (sulphonic) significantly increased R and D and decrease v for both soils. Addition of sulphonic increased R of atrazine by 12 and 26 folds for loam and sandy loam soils, respectively. On the other hand non-ionic surfactants seemed to decrease R, especially in sandy loam soil, thus facilitating atrazine leaching through soil. Non-equilibrium conditions seemed to govern atrazine transport in the presence of surfactants; double peaks in breakthrough curves were observed, indicating a need for mathematical models to account for such phenomena. Atrazine dispersion and tailing seemed to be higher through Elora loam compared to Caledon sandy loam due to higher aggregation of the Elora soil.     

Irrigation model of bleached Kraft mill wastewater through volcanic soil as a pollutants attenuation process

An irrigation process through volcanic soil columns was evaluated for bleached Kraft mill effluent pollutants retention. The system was designed to remove color and phenolic compounds and a simple kinetic model for determining the global mass transfer coefficient and the adsorption rate constant was used. The results clearly indicate that the global mass transfer coefficient values (K(c)a) and the adsorption rate constants are higher for the irrigation processes onto acidified soil. This means that the pretreatment of washing the volcanic soil with an acid solution has a positive effect on the adsorption rate for both pollutant groups. The enhanced adsorption capacity is partially explained by the activation of the metal oxides present in the soil matrix during the acid washing process. Increasing the flow rate from 1.5 to 2.5 ml/min yielded higher (K(c)a) values and adsorption rate constants for both pollutant groups. For instance, regarding color adsorption onto acidified soil, there is an increment of 43% in the (K(c)a) value for the experiment with a flow rate of 2.5 ml/min. Increasing the porosity of the column from 0.55 to 0.59, yielded a decrease in the (K(c)a) values for color and phenolic compounds adsorption processes. Onto natural soil for example, these decreases reached 21% and 24%, respectively. Therefore, the (K(c)a) value is dependent on both the liquid-phase velocity (external resistance) and the soil fraction in the column (internal resistance); making forced convection and diffusion to be the main transport mechanisms involved in the adsorption process. Analyzing the adsorption rate constants (K(c)a)/m, phenolic compounds and color adsorption rates onto acidified soil of 2.25 x 10(-6) and 2.62 x 10(-6) l/mg min were achieved for experiment 1. These adsorption rates are comparable with other adsorption systems and adsorbent materials.     

污泥活性炭固定床吸附甲苯

采用ZnCl2化学活化法制备的污泥活性炭用于固定床吸附甲苯实验,研究了吸附剂床层厚度、甲苯初始浓度、气体线速等操作条件对其吸附性能的影响,并建立了污泥活性炭吸附甲苯的传质模型,计算不同操作条件下的传质区长度和物质总传质系数。结果表明,当床层厚度从3 cm增加到7 cm,传质区长度基本不变,物质总传质系数增加了66.66%;当甲苯初始浓度增大3.36倍,传质区长度增加了17.19%,而物质总传质系数降低了22.05%;在实验流速范围内,气体线速变为原来的2倍,物质总传质系数增加了65.27%,说明污泥活性炭对甲苯吸附属于外扩散控制。     

Transport of silver nanoparticles (AgNPs) in soil

The effect of soil properties on the transport of silver nanoparticles (AgNPs) was studied in a set of laboratory column experiments, using different combinations of size fractions of a Mediterranean sandy clay soil. The AgNPs with average size of ∼30 nm yielded a stable suspension in water with zeta potential of −39 mV. Early breakthrough of AgNPs in soil was observed in column transport experiments. AgNPs were found to have high mobility in soil with outlet relative concentrations ranging from 30% to 70%, depending on experimental conditions. AgNP mobility through the column decreased when the fraction of smaller soil aggregates was larger. The early breakthrough pattern was not observed for AgNPs in pure quartz columns nor for bromide tracer in soil columns, suggesting that early breakthrough is related to the nature of AgNP transport in natural soils. Micro-CT and image analysis used to investigate structural features of the soil, suggest that soil aggregate size strongly affects AgNP transport in natural soil. The retention of AgNPs in the soil column was reduced when humic acid was added to the leaching solution, while a lower flow rate (Darcy velocity of 0.17 cm/min versus 0.66 cm/min) resulted in higher retention of AgNPs in the soil. When soil residual chloride was exchanged by nitrate prior to column experiments, significantly improved mobility of AgNPs was observed in the soil column. These findings point to the importance of AgNP-soil chemical interactions as a retention mechanism, and demonstrate the need to employ natural soils rather than glass beads or quartz in representative experimental investigations.     

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.     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

Heavy metals removal in fixed-bed column by the macro fungus Pycnoporus sanguineus

The ability of Pycnoporus sanguineus to adsorb heavy metals from aqueous solution was investigated in fixed-bed column studies. The experiments were conducted to study the effect of important design parameters such as column bed height, flow rate and initial concentration of solution. The breakthrough profiles were obtained in these studies. A mathematical model based on external mass transfer and pore diffusion was used for the prediction of mass transfer coefficient and effective diffusivity of metals in macro-fungi bed. Experimental breakthrough profiles were compared with the simulated breakthrough profiles obtained from the mathematical model. Bed Depth Service Time (BDST) model was used to analyse the experimental data and evaluated the performance of biosorption column. The BDST model parameters needed for the design of biosorption columns were evaluated for lead, copper and cadmium removal in the column. The columns were regenerated by eluting the metal ions using 0.1 M hydrochloric acid solution after the adsorption studies. The columns were subjected to repeated cycles of adsorption of same metal ions and desorption to evaluate the removal efficiency after adsorption-desorption.     

Development of a continuously stirred flow cell for investigating sorption mass transfer

This paper introduces a new reversible-flow design for a continuously stirred reactor used to study sorption mass transfer in soil and solvent systems. The stirred reactor has potential advantages over conventional packed column or batch reactors because it isolates intraparticle sorption rate limitations from advective-dispersive transport, yet allows changes to flux through the reactor for analysis of sorption kinetics under dynamic conditions. Previously, stirred reactors have often failed due to clogging of sediment on the effluent frit. The reverse-flow backwashing design allows longer life and higher confidence in maintaining mixed conditions than previous designs. Mass transfer 'rate coefficients estimated from stirred and column experiments are compared; both techniques produced results consistent with a published correlation. The data also show that fitted sorption mass transfer coefficients can be strongly dependent on the choice of equilibrium partition coefficient (i.e. batch or first-moment derived values), and that the conventional two-site sorption kinetics model fails to accurately predict sorption mass transfer in the presence of changing solvent velocity through the reactor.     

Factors affecting mass transfer limited biodegradation in saturated porous media

Microbial degradation rates in the subsurface are not only limited by the physiological capacity of the organisms, but also by inefficient supply of nutrients to the microbes. Although mass transfer limitation of biodegradation in the subsurface has been postulated for years, experimental evidence is still scarce. In the column experiments described here, diffusive transport of 4-nitroanisole from the bulk solution to cells of Rhodococcus opacus strain AS2 immobilized on glass beads or sand appeared to be responsible for the slow transformation rates observed. Assuming steady state, we applied a coupled transformation/transport equation to these data (Best equation) and apparent bead-related mass transfer coefficients were found to increase in proportion to the surface area covered with bacteria. This implies that mass transfer coefficients for individual cells remained constant. In an idealized oligotrophic environment where cells are only loosely clustered and do not shield each other, we would therefore expect biodegradation rates to be independent from the longitudinal distribution of the total biomass along a given flow path. Moreover, apparent mass transfer coefficients increased with the grain size of the column fillings, but did not change upon varying the flow rate. With a limiting external transport step, overall transformation fluxes do not become saturated at concentrations as low as predicted for Michaelis-Menten-type kinetics. Mass transfer limitation thus offers a justification for the common assumption that biodegradation rates in the subsurface follow first order kinetics in a wide concentration range.     

Evidence of one-dimensional scale-dependent fractional advection-dispersion

A semi-analytical inverse method and the corresponding program FADEMain for parameter estimation of the fractional advection-dispersion equation (FADE) were developed in this paper. We have analyzed Huang et al.'s [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.] laboratory experimental data of conservative solute transport in 12.5-m long homogeneous and heterogeneous soil columns to test the non-Fickian dispersion theory of FADE. The dispersion coefficient was calculated by fitting the analytical solution of FADE to the measured data at different transport scales. We found that the dispersion coefficient increased exponentially with transport scale for the homogeneous column, whereas it increased with transport scale in a power law function for the heterogeneous column. The scale effect of the dispersion coefficient in the heterogeneous soil was much more significant comparing to that in the homogeneous soil. The increasing rate of dispersion coefficient versus transport distance was smaller for FADE than that for the advection-dispersion equation (ADE). Finite difference numerical approximations of the scale-dependent FADE were established to interpret the experimental results. The numerical solutions were found to be adequate for predicting scale-dependent transport in the homogeneous column, while the prediction for the heterogeneous column was less satisfactory.     

Airflow dispersion in unsaturated soil

Dispersion data is abundant for water flow in the saturated zone but is lacking for airflow in unsaturated soil. However, for remediation processes such as soil vapour extraction, characterization of airflow dispersion is necessary for improved modelling and prediction capabilities. Accordingly, gas-phase tracer experiments were conducted in five soils ranging from uniform sand to clay at air-dried and wetted conditions. The disturbed soils were placed in one-dimensional stainless steel columns, with sulfur hexafluoride used as the inert tracer. The tested interstitial velocities were typical of those present in the vicinity of a soil vapour extraction well, while wetting varied according to the water-holding capacity of the soils. Results gave dispersivities that varied between 0.42 and 2.6 cm, which are typical of values in the literature. In air-dried soils, dispersion was found to increase with the pore size variability of the soil. For wetted soils, particle shape was an important factor at low water contents, while at high water contents, the proportion of macroporous space filled with water was important. The relative importance of diffusion decreased with increasing interstitial velocity and water content and was, in general, found to be minor compared to mechanical mixing across all conditions studied.     

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.     

Nonequilibrium sorption of phenols onto geosorbents: the impact of pH on intraparticle mass transfer

While the pH effect on sorption equilibrium of weak acids on natural sorbents was investigated in a number of studies, less is known about the pH dependence of sorption kinetics. This paper investigates the impact of pH on sorption kinetics during the transport of some selected phenols through a sandy aquifer material. Breakthrough curves measured in column experiments were analyzed using a mass transfer based nonequilibrium model designated as dispersed flow, film and particle diffusion model (DF-FPDM). In this model, the rate limiting intraparticle diffusion is characterized by the mass transfer coefficient, kSaV, which can be determined from breakthrough curves by curve fitting. The experimental results indicate that the kSaV is pH-dependent and inversely correlated with the pH-dependent distribution coefficient, K(d,app). Regression equations are presented that may be used to estimate approximate values of intraparticle mass transfer coefficients on the basis of experimentally determined or LFER predicted distribution coefficients.     

Characteristics of nitrate horizontal transport in a paddy field of the Tai Lake region,China

The characteristics of nitrate horizontal transport in a major paddy soil, Wu Shan soil in the Tai Lake region, were studied. The concentration of nitrate during horizontal movement decreased with the increasing in distance from the tracer source, the change following a logarithmic function. The concentration of the nitrate was strongly correlated with the soil moisture content, as an exponential function. The horizontal transport velocity of nitrate was significantly correlated with the distance of the tracer source as power function. Therefore, the velocity of nitrate horizontal transport was controlled by the concentration gradient of nitrate, and soil water potential gradient from beginning to the 20 cm mark in the horizontal column. However, the velocity of nitrate horizontal was stable beyond 20 cm, where it was controlled by soil matric potential.     

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.     

Street canyon ventilation and atmospheric turbulence

Operational models for pollutant dispersion in urban areas require an estimate of the turbulent transfer between the street canyons and the overlying atmospheric flow. To date, the mechanisms that govern this process remain poorly understood. We have studied the mass exchange between a street canyon and the atmospheric flow above it by means of wind tunnel experiments. Fluid velocities were measured with a Particle Image Velocimetry system and passive scalar concentrations were measured using a Flame Ionisation Detector. The mass-transfer velocity between the canyon and the external flow has been estimated by measuring the cavity wash-out time. A two-box model, used to estimate the transfer velocity for varying dynamical conditions of the external flow, has been used to interpret the experimental data. This study sheds new light on the mechanisms which drive the ventilation of a street canyon and illustrates the influence of the external turbulence on the transfer process.     

Transverse vertical dispersion in groundwater and the capillary fringe

Transverse dispersion is the most relevant process in mass transfer of contaminants across the capillary fringe (both directions), dilution of contaminants, and mixing of electron acceptors and electron donors in biodegrading groundwater plumes. This paper gives an overview on literature values of transverse vertical dispersivities alpha(tv) measured at different flow velocities and compares them to results from well-controlled laboratory-tank experiments on mass transfer of trichloroethene (TCE) across the capillary fringe. The measured values of transverse vertical dispersion in the capillary fringe region were larger than in fully saturated media, which is credited to enhanced tortuosity of the flow paths due to entrapped air within the capillary fringe. In all cases, the values observed for alpha(tv) were < 1 mm. The new measurements and the literature values indicate that alpha(tv) apparently declines with increasing flow velocity. The latter is attributed to incomplete diffusive mixing at the pore scale (pore throats). A simple conceptual model, based on the mean square displacement and the pore size accounting for only partial diffusive mixing at increasing flow velocities, shows very good agreement with measured and published data.     

Choice of dispersion coefficients in reactive transport calculations on smoothed fields

Local dispersion dominates the mixing of compounds that are introduced separately into the subsurface and do not partition into any other than the aqueous phase. Thus, reactions between these compounds are controlled by dispersive mixing if they are not limited by kinetics. I quantify longitudinal dispersive mixing by the longitudinal effective dispersion coefficient of a conservative tracer introduced by a point-like injection [Water Resour. Res. 36 (12) (2000) 3591-3604]. In the upscaling of mixing-controlled reactive transport, I apply the mean velocity and the effective dispersion coefficient to the macroscopic transport calculations, whereas the reactive parameters on the macro-scale are identical to those on the local scale. The applicability of the approach is demonstrated for the transport of compounds undergoing a second-order irreversible bimolecular reaction. Ten realizations of a two-dimensional heterogeneous log-conductivity field are considered. Using the effective dispersion parameters, the overall mass balance is met in the ensemble average, whereas solute spreading is underestimated. I assess the lack of spreading by the difference between the expected macrodispersion and effective dispersion coefficients. I extend the approach to simulations on log-conductivity fields obtained by kriging of regularly spaced conductivity measurements. These fields contain the large-scale features of the true fields but do not resolve the small-scale variability. For the calculations on the kriged fields, the corresponding conditional covariance is substituted into the analytical expressions of effective dispersion, yielding a correction effective dispersion coefficient. The comparison between simulations on the fully resolved fields and on the kriged fields indicates that the approach is valid for wide plumes meeting the ergodicity condition. The high variability of mixing on small scales unresolved by kriging, however, leads to severe uncertainty when mixing-controlled reactions are predicted for narrow plumes.