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.     

Grid lysimeter study of steady state chloride transport in two Spodosol types using TDR and wick samplers.

Solute transport in soils is affected by soil layering and soil-specific morphological properties. We studied solute transport in two sandy Spodosols: a dry Spodosol developed under oxidizing conditions of relatively deep groundwater and a wet Spodosol under periodically reducing conditions above a shallow groundwater table. The wet Spodosol is characterized by a diffuse and heterogeneous humus-B-horizon (i.e., Spodic horizon), whereas the dry Spodosol has a sharp Spodic horizon. Drainage fluxes were moderately variable with a coefficient of variation (CV) of 25% in the wet Spodosol and 17% in the dry Spodosol. Solute transport in 1-m-long and 0.8-m-diameter soil columns was investigated using spatial averages of solute concentrations measured by a network of 36 Time Domain Reflectometry (TDR) probes. In the dry Spodosol, solute transport evolves from stochastic-convective to convective-dispersive at a depth of 0.25 m, coinciding with the depth of the Spodic horizon. Chloride breakthrough at the bottom of the soil columns was adequately well predicted by a convection-dispersion model. In the wet Spodosol, solute transport was heterogeneous over the entire depth of the column. Chloride breakthrough at 1 m depth was predicted best using a stochastic-convective transport model. The TDR sampling volume of 36 probes was too small to capture the heterogeneous flow and concomitant transport in the wet Spodosol.     

Transport modeling of copper and cadmium with linear and nonlinear retardation factors

The predictive accuracy of using the one-dimensional advection–dispersion equation to evaluate the fate and transport of solute in a soil column is usually dependent on the proper determination of chemical retardation factors. Typically, the distribution coefficient (K_d) obtained by fitting the linear sorption isotherm has been extensively used to consider general geochemical reactions on solute transport in a low-concentration range. However, the linear distribution coefficient cannot be adequately utilized to describe the solute fate at a higher concentration level. This study employed the nonlinear equilibrium-controlled sorption parameters to determine the retardation factor used in column leaching experiments. Copper and cadmium transportation in a lateritic silty-clay soil column was examined. Through the explicit finite-difference calculations with a third-order total-variation-diminishing (TVD) numerical solution scheme, all results of the theoretical copper and cadmium breakthrough curves (BTCs) simulated by using the Freundlich nonlinear retardation factors revealed good agreement with the experimental observations.     

One-dimensional simulation of solute transfer in saturated-unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous-heterogeneous, saturated-unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection-dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term. The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

An interpretation of potential scale dependence of the effective matrix diffusion coefficient

Matrix diffusion is an important process for solute transport in fractured rock, and the matrix diffusion coefficient is a key parameter for describing this process. Previous studies have indicated that the effective matrix diffusion coefficient values, obtained from a large number of field tracer tests, are enhanced in comparison with local values and may increase with test scale. In this study, we have performed numerical experiments to investigate potential mechanisms behind possible scale-dependent behavior. The focus of the experiments is on solute transport in flow paths having geometries consistent with percolation theories and characterized by multiple local flow loops formed mainly by small-scale fractures. The water velocity distribution through a flow path was determined using discrete fracture network flow simulations, and solute transport was calculated using a previously derived impulse-response function and a particle-tracking scheme. Values for effective (or up-scaled) transport parameters were obtained by matching breakthrough curves from numerical experiments with an analytical solution for solute transport along a single fracture. Results indicate that a combination of local flow loops and the associated matrix diffusion process, together with scaling properties in flow path geometry, seems to be an important mechanism causing the observed scale dependence of the effective matrix diffusion coefficient (at a range of scales).     

A numerical model for transient-hysteretic flow and solute transport in unsaturated porous media

A two-dimensional flow and transport model was developed for simulating transient water flow and nonreactive solute transport in heterogeneous, unsaturated porous media containing air and water. The model is composed of a unique combination of robust and accurate numerical algorithms for solving the Richards', Darcy flux, and advection-dispersion equations. The mixed form of Richards' equation is solved using a finite-element formulation and a modified Picard iteration scheme. Mass lumping is employed to improve solution convergence and stability behavior. The flow algorithm accounts for hysteresis in the pressure head-water content relationship. Darcy fluxes are approximated with a Galerkin and Petrov-Galerkin finite-element method developed for random heterogeneous porous media. The transport equation is solved using an Eulerian-Lagrangian method. A multi-step, fourth-order Runge-Kutta, reverse particle tracking technique and a quadratic-linear interpolation scheme are shown to be superior for determining the advective concentration. A Galerkin finite-element method is used for approximating the dispersive flux. The unsaturated flow and transport model was applied to a variety of rigorous problems and was found to produce accurate, mass-conserving solutions when compared to analytical solutions and published numerical results.     

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.     

Non-Fickian dispersion in a single fracture

Solute transport in fractured rocks is of major interest in many applications, from the petroleum industry to ground water management. This work focuses on the dispersion process in a transparent replica of a real single fracture. The fracture exhibits strong changes in heterogeneity, with the first half very heterogeneous and the second half fairly homogeneous. Three models have been used to interpret the tracer experiments: the classical advection-dispersion equation (ADE), the continuous time random walk (CTRW), and the stratified model. The main goals were to test these models and to study possible correlations between fitting parameters and heterogeneities. As expected, the solution derived from the ADE equation appears to be unable to model long-time tailing behavior. On the other hand, the results confirm the CTRW robustness and the coefficient beta seems well correlated to heterogeneities. Finally, the stratified model is also able to describe non-Fickian dispersion. The parameters defined by this model are correlated to the heterogeneities of the fracture.     

Bacteriophage transport through a fining-upwards sedimentary sequence: laboratory experiments and simulation

A column containing four concentric layers of progressively finer-grained glass beads (graded column) was used to study the transport of the bacteriophage T7 in water flowing parallel to layering through a fining-upwards (FU) sedimentary structure. By passing a pulse of T7, and a conservative solute tracer upwards through a column packed with a single bead size (uniform column), the capacity of each bead type to attenuate the bacteriophage was determined. Solute and bacteriophage responses were modelled using an analytical solution to the advection-dispersion equation, with first-order kinetic deposition simulating bacteriophage attenuation. Resulting deposition constants for different flow velocities indicated that filtration theory-determined values differed from experimentally determined values by less than 10%. In contrast, the responses of solute and bacteriophage tracers passing upwards through graded columns could not be reproduced with a single analytical solution. However, a flux-weighted summation of four one-dimensional advective-dispersive analytical terms approximated solute breakthrough curves. The prolonged tailing observed in the resulting curve resembled that typically generated from field-based tracer test data, reflecting the potential importance of textural heterogeneity in the transport of dissolved substances in groundwater. Moreover, bacteriophage deposition terms, determined from filtration theory, reproduced the T7 breakthrough curve once desorption and inactivation on grain surfaces were incorporated. To evaluate the effect of FU sequences on mass transport processes in more detail, bacteriophage passage through sequences resembling those sampled from a FU bed in a fluvioglacial gravel pit were carried out using an analogous approach to that employed in the laboratory. Both solute and bacteriophage breakthrough responses resembled those generated from field-based test data and in the graded column experiments. Comparisons with the results of simulations using averaged hydraulic conductivities show that simulations employing averaged parameters overestimate bacteriophage travel times and underestimate masses recovered and peak concentrations.     

The effects of matrix diffusion on radionuclide migration in rock column experiments

Rock column experiments were performed to examine the effects of matrix diffusion and hydrodynamic dispersion on the migration of radionuclides at the laboratory scale. Tritiated water and chloride transportation was studied in intact mica gneiss and in altered more porous tonalite columns with narrow flow channels. The column diffusion properties were estimated prior to water flow experiments using the gas diffusion method with helium as the tracer gas. The numerical compartment model for advection and dispersion, with and without matrix diffusion, was used to interpret the tracer transport in the columns. Matrix diffusion behavior was also distinguished from dominating hydrodynamic dispersion in rock column experiments at the slowest water flow rates.     

Effects of pH and ionic strength on sulfamethoxazole and ciprofloxacin transport in saturated porous media

Many antibiotics regarded as emerging contaminants have been frequently detected in soils and groundwater; however, their transport behaviors in soils remain largely unknown. This study examined the transport of two antibiotics, sulfamethoxazole (SMZ) and ciprofloxacin (CIP), in saturated porous media. Laboratory columns packed with quartz sand was used to test the effects of solution pH and ionic strength (IS) on their retention and transport. The results showed that these two antibiotics behaved differently in the saturated sand columns. In general, SMZ manifested a much higher mobility than CIP for all experimental conditions tested. Almost all SMZ transported through the columns within one pore volume in deionized water (i.e., pH=5.6, IS=0), but no CIP was detected in the effluents under the same condition after extended column flushing. Perturbations in solution pH (5.6 and 9.5) and IS (0 and 0.1M) showed no effect on SMZ transport in the saturated columns. When pH increased to 9.5, however, ~93% of CIP was eluted from the sand columns. Increase of IS from 0 to 0.1M also slightly changed the distribution of adsorbed CIP within the sand column at pH 5.6, but still no CIP was detected in the effluents. A mathematical model based on advection-dispersion equation coupled with equilibrium and kinetic reactions successfully simulated the transport of the antibiotics in water-saturated porous media with R(2)=0.99.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Simulation and analysis of solute transport in 2D fracture/pipe networks: the SOLFRAC program

The Time Domain Random Walk (TDRW) method has been recently developed by Delay and Bodin [Delay, F. and Bodin, J., 2001. Time domain random walk method to simulate transport by advection-dispersion and matrix diffusion in fracture networks. Geophys. Res. Lett., 28(21): 4051-4054.] and Bodin et al. [Bodin, J., Porel, G. and Delay, F., 2003c. Simulation of solute transport in discrete fracture networks using the time domain random walk method. Earth Planet. Sci. Lett., 6566: 1-8.] for simulating solute transport in discrete fracture networks. It is assumed that the fracture network can reasonably be represented by a network of interconnected one-dimensional pipes (i.e. flow channels). Processes accounted for are: (1) advection and hydrodynamic dispersion in the channels, (2) matrix diffusion, (3) diffusion into stagnant zones within the fracture planes, (4) sorption reactions onto the fracture walls and in the matrix, (5) linear decay, and (6) mass sharing at fracture intersections. The TDRW method is handy and very efficient in terms of computation costs since it allows for the one-step calculation of the particle residence time in each bond of the network. This method has been programmed in C++, and efforts have been made to develop an efficient and user-friendly software, called SOLFRAC. This program is freely downloadable at the URL (labo.univ-poitiers.fr/hydrasa/intranet/telechargement.htm). It calculates solute transport into 2D pipe networks, while considering different types of injections and different concepts of local dispersion within each flow channel. Post-simulation analyses are also available, such as the mean velocity or the macroscopic dispersion at the scale of the entire network. The program may be used to evaluate how a given transport mechanism influences the macroscopic transport behaviour of fracture networks. It may also be used, as is the case, e.g., with analytical solutions, to interpret laboratory or field tracer test experiments performed in single fractures.     

Field-scale effective matrix diffusion coefficient for fractured rock: results from literature survey

Matrix diffusion is an important mechanism for solute transport in fractured rock. We recently conducted a literature survey on the effective matrix diffusion coefficient, D_m^e, a key parameter for describing matrix diffusion processes at the field scale. Forty field tracer tests at 15 fractured geologic sites were surveyed and selected for the study, based on data availability and quality. Field-scale D_m^e values were calculated, either directly using data reported in the literature, or by reanalyzing the corresponding field tracer tests. The reanalysis was conducted for the selected tracer tests using analytic or semi-analytic solutions for tracer transport in linear, radial, or interwell flow fields. Surveyed data show that the scale factor of the effective matrix diffusion coefficient (defined as the ratio of D_m^e to the lab-scale matrix diffusion coefficient, D_m, of the same tracer) is generally larger than one, indicating that the effective matrix diffusion coefficient in the field is comparatively larger than the matrix diffusion coefficient at the rock-core scale. This larger value can be attributed to the many mass-transfer processes at different scales in naturally heterogeneous, fractured rock systems.Furthermore, we observed a moderate, on average trend toward systematic increase in the scale factor with observation scale. This trend suggests that the effective matrix diffusion coefficient is likely to be statistically scale-dependent. The scale-factor value ranges from 0.5 to 884 for observation scales from 5 to 2000 m. At a given scale, the scale factor varies by two orders of magnitude, reflecting the influence of differing degrees of fractured rock heterogeneity at different geologic sites. In addition, the surveyed data indicate that field-scale longitudinal dispersivity generally increases with observation scale, which is consistent with previous studies. The scale-dependent field-scale matrix diffusion coefficient (and dispersivity) may have significant implications for assessing long-term, large-scale radionuclide and contaminant transport events in fractured rock, both for nuclear waste disposal and contaminant remediation.     

Analytical model of solute transport by unsteady unsaturated gravitational infiltration

Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils.     

Effects of pore volume-transmissivity correlation on transport phenomena

The relevant velocity that describes transport phenomena in a porous medium is the pore velocity. For this reason, one needs not only to describe the variability of transmissivity, which fully determines the Darcy velocity field for given source terms and boundary conditions, but also any variability of the pore volume. We demonstrate that hydraulically equivalent media with exactly the same transmissivity field can produce dramatic differences in the displacement of a solute if they have different pore volume distributions. In particular, we demonstrate that correlation between pore volume and transmissivity leads to a much smoother and more homogeneous solute distribution. This was observed in a laboratory experiment performed in artificial fractures made of two plexiglass plates into which a space-dependent aperture distribution was milled. Using visualization by a light transmission technique, we observe that the solute behaviour is much smoother and more regular after the fractures are filled with glass powder, which plays the role of a homogeneous fault gouge material. This is due to a perfect correlation between pore volume and transmissivity that causes pore velocity to be not directly dependent on the transmissivity, but only indirectly through the hydraulic gradient, which is a much smoother function due to the diffusive behaviour of the flow equation acting as a filter. This smoothing property of the pore volume-transmissivity correlation is also supported by numerical simulations of tracer tests in a dipole flow field. Three different conceptual models are used: an empty fracture, a rough-walled fracture filled with a homogeneous material and a parallel-plate fracture with a heterogeneous fault gouge. All three models are hydraulically equivalent, yet they have a different pore volume distribution. Even if piezometric heads and specific flow rates are exactly the same at any point of the domain, the transport process differs dramatically. These differences make it important to discriminate in situ among different conceptual models in order to simulate correctly the transport phenomena. For this reason, we study the solute breakthrough and recovery curves at the extraction wells. Our numerical case studies show that discrimination on the basis of such data might be impossible except under very favourable conditions, i.e. the integral scale of the transmissivity field has to be known and small compared to the dipole size. If the latter conditions are satisfied, discrimination between the rough-walled fracture filled with a homogeneous material and the other two models becomes possible, whereas the parallel-plate fracture with a heterogeneous fault gouge and the empty fracture still show identifiability problems. The latter may be solved by inspection of aperture and pressure testing.     

Semi-analytical solution of two-dimensional sharp interface LNAPL transport models

In this paper, we present semi-analytical solutions for two-dimensional equations governing transport of Light Non-Aqueous Phase Liquids (LNAPL) in unconfined aquifers. The proposed model is based on sharp interface displacement and steady groundwater flow assumptions, where both the water–LNAPL interface and the LNAPL–air interface are represented as sharp interfaces. In the case of steady groundwater flow, these equations can be reduced to a two-dimensional nonlinear solute transport equation, with the LNAPL thickness in the free product lens being the primary unknown variable. The linearized form of this solute transport equation falls into the category of two-dimensional transport equation with time-dependent dispersion coefficients. This equation can be solved analytically for an infinite domain region. In this paper, the general form of the analytical solution for the transport equation, as well as the solutions for some specific cases are presented. To demonstrate the utility of the proposed solution, numerical results obtained for two example problems are discussed and presented comparatively with a finite-element solution and other more restrictive solutions available in the literature. Although the solutions discussed in this paper have some simplifying assumptions, such as sharp-interfaces between fluid phases, steady groundwater flow and homogeneous aquifer properties, the semi-analytical solutions presented in this study may be used effectively as bench mark solutions in evaluating LNAPL migration in the subsurface. These solutions are simple and cost effective to implement and may be used in the calibration of other more complex numerical solutions that can be found in the literature.     

Analytical solutions for reactive transport under an infiltration-redistribution cycle

Transport of reactive solute in unsaturated soils under an infiltration-redistribution cycle is investigated. The study is based on the model of vertical flow and transport in the unsaturated zone proposed by Indelman et al. [J. Contam. Hydrol. 32 (1998) 77], and generalizes it by accounting for linear nonequilibrium kinetics. An exact analytical solution is derived for an irreversible desorption reaction. The transport of solute obeying linear kinetics is modeled by assuming equilibrium during the redistribution stage. The model which accounts for nonequilibrium during the infiltration and assumes equilibrium at the redistribution stage is termed partial equilibrium infiltration-redistribution model (PEIRM). It allows to derive approximate closed form solutions for transport in one-dimensional homogeneous soils. These solutions are further applied to computing the field-scale concentration by adopting the Dagan and Bresler [Soil Sci. Soc. Am. J. 43 (1979) 461] column model. The effect of soil heterogeneity on the solute spread is investigated by modeling the hydraulic saturated conductivity as a random function of horizontal coordinates. The quality of the PEIRM is illustrated by calculating the critical values of the Damk?hler number which provide the achievable accuracy in estimating the solute mass in the mobile phase. The distinguishing feature of transport during the infiltration-redistribution cycle as compared to that of infiltration only is the finite depth of solute penetration. For irreversible desorption, the maximum solute penetration W/theta(r) is determined by the amount of applied water W and the residual water content theta(r). For sorption-desorption kinetics, the maximum depth of penetration z(r)(e, infinity ) also depends on the ratio between the rate of application and the column-saturated conductivity. It is shown that z(r)(e, infinity ) is bounded between the depths W/(theta(r)+K(d)) and W/theta(r) corresponding to the maximum solute penetration for equilibrium transport and for irreversible desorption, respectively. This feature of solute penetration explains the unusual phenomena of plume contraction after an initial period of spreading [Lessoff, S.C., Indelman, P., Dagan, G., 2002. Solute transport in infiltration-redistribution cycles in heterogeneous soils. In Raats, P.A.C., Smiles, D.,Warrick, A.W. (Eds), Environmental Mechanics: Water, Mass and Energy Transport in the Biosphere. American Geophysical Union, pp. 133-144]. Unlike transport under equilibrium conditions, when the solute is completely concentrated at the front, the solute under nonequilibrium conditions is spread out behind the front. Heterogeneity leads to additional spreading of the plume.     

Influence of methanol on retention of hydrophobic organic chemicals in soil leaching column chromatography

The influence of methanol in methanol-water mixed eluents on the capacity factor (k'), an important parameter which could depict leaching potential of hydrophobic organic chemicals (HOCs) in soil leaching column chromatography (SLCC), was investigated. Two reference soils, GSE 17201 obtained from Bayer Landwirtschaftszentrum, Monheim, Germany and SP 14696 from LUFA, Spencer, Germany, were used as packing materials in soil columns, and isocratic elution with methanol-water mixtures at different volume fractions of methanol (phi) were tested. Short-term exposure of the column (packed with the GSE 17201 soil) to the eluents increased solute retention by a certain (23% log-unit) degree evaluated through a correlation with the retention on the same soil column but unpreconditioned by methanol-containing eluents. Long-term exposure of soil columns to the eluents did not influence the solute retention. A log-linear equation, log k' = log k'(w) - S(phi), could well and generally describe the retention of HOCs in SLCC. For the compounds of homologous series, logk'(w) had good linear relationship with S, indicating the hydrophobic partition mechanism existing in the retention process.