Modeling colloid transport for performance assessment

The natural system is expected to contribute to isolation at the proposed high-level nuclear waste (HLW) geologic repository at Yucca Mountain, NV (YM). In developing performance assessment (PA) computer models to simulate long-term behavior at YM, colloidal transport of radionuclides has been proposed as a critical factor because of the possible reduced interaction with the geologic media. Site-specific information on the chemistry and natural colloid concentration of saturated zone groundwaters in the vicinity of YM is combined with a surface complexation sorption model to evaluate the impact of natural colloids on calculated retardation factors (RF) for several radioelements of concern in PA. Inclusion of colloids into the conceptual model can reduce the calculated effective retardation significantly. Strongly sorbed radionuclides such as americium and thorium are most affected by pseudocolloid formation and transport, with a potential reduction in RF of several orders of magnitude. Radioelements that are less strongly sorbed under YM conditions, such as uranium and neptunium, are not affected significantly by colloid transport, and transport of plutonium in the valence state is only moderately enhanced. Model results showed no increase in the peak mean annual total effective dose equivalent (TEDE) within a compliance period of 10,000 years, although this is strongly dependent on container life in the base case scenario. At longer times, simulated container failures increase and the TEDE from the colloidal models increased by a factor of 60 from the base case. By using mechanistic models and sensitivity analyses to determine what parameters and transport processes affect the TEDE, colloidal transport in future versions of the TPA code can be represented more accurately.     

Laboratory migration experiments with radionuclides and natural colloids in a granite fracture

Natural colloids in groundwater could facilitate radionuclide transport, provided the colloids are mobile, are present in sufficient concentrations and can adsorb radionuclides. This paper describes the results of a laboratory migration study carried out with combinations of radionuclides and natural colloids within a fracture in a large granite block to experimentally determine the impact of colloids on radionuclide transport. The 85Sr used in this study is an example of a moderately sorbing radionuclide, while the 241Am is typical of a strongly sorbed radionuclide with very low solubility. The natural colloids used in this study were isolated from granite groundwater from Atomic Energy of Canada (AECL) Underground Research Laboratory (URL), and consisted of mostly 1-10 nm organic colloids, along with lesser amounts of 10-450 nm colloids (organics and aluminosilicates). The measured coefficients for radionuclide sorption onto these colloids were between 3 x 10(2) and 1 x 10(3) ml/g for 85Sr, and between 7 x 10(4) and 7 x 10(5) mg/l for 241Am. The 85Sr sorption on the natural colloids appeared to be reversible. Migration experiments in the granite block were carried out by establishing a flow field between two boreholes (out of a total of nine) intersecting a main horizontal fracture. These experiments showed that dissolved 85Sr behaved as a moderately sorbing tracer, while dissolved 241Am was completely adsorbed by the fracture surfaces and showed no evidence of transport. However, when natural colloids were injected together with dissolved 241Am, a small amount of 241Am transport was observed, demonstrating the ability of natural colloids to facilitate the transport of radionuclides with low solubility. Natural colloids had only a minor effect on the transport of 85Sr. In a separate experiment to test the effect of higher colloid concentrations on 85Sr migration, synthetic colloids were produced from Avonlea bentonite. The introduction of a relatively high concentration of bentonite colloids actually reduced 85Sr transport because, compared to natural colloids, the bentonite colloids were less mobile and they sorbed 85Sr more strongly.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

The kinetic stability of colloid-associated plutonium: settling characteristics and species transformation

The colloid-associated plutonium-239, as the dominant species of Pu in natural environment, was formed through sorption of Pu onto in situ colloids. In case of chemical perturbations present in pore water, the fate and transport of Pu would be therefore impacted by changes in sorption affinity of Pu for the colloid surfaces. The present study reveals that colloidal ²³⁹Pu exhibited the kinetic stability in two respects. First, in situ colloids isolated from the vadose zone sediments at Lop Nor, when in contact with solutions of high ion concentrations or low pH, were significantly aggregated and then exhibited fast settling. Kinetics settling characteristics were described by the parameters, including settling index and characteristic time. Second, Pu dissociation from colloid surfaces occurred immediately after the introduction of Na⁺. However, the dissolved species was still unstable and had the potential for re-association with the fraction of colloids that had not settled out from the suspensions due to small size and then remained in suspension. This implies that Pu sorption sites on initial colloids were changed to the sites of suspended colloids.     

Colloid‐facilitated transport of metolachlor through intact soil columns

Abstract

This study evaluated the role of water dispersible colloids with diverse physicochemical and mineralogical characteristics in facilitating the transport of metolachlor through macropores of intact soil columns. The soil columns represented upper solum horizons of an Alfisol in the Bluegrass region of Kentucky. Three different colloid suspensions tagged with metolachlor [2‐chloro‐N‐(2‐ethyl‐6‐methylphenyl)‐N‐(2‐methoxy‐l‐methylethyl)acetamide] were introduced at a constant flux into undisturbed soil columns. The eluents were collected and analyzed periodically for colloid and metolachlor concentrations. Colloid recovery in the eluents ranged from 54 to 90 %. The presence of colloids enhanced the transport of metolachlor by 22 to 70 % depending on the colloid type and mobility. Colloids with higher pH, organic carbon, cation exchange capacity (CEC), total exchangeable bases (TEB), surface area (SA), and electrophoretic mobility (EM), showed better mobility, greater affinity for interaction with the herbicide and, thus, greater potential to co‐transport metolachlor. In contrast, increased level of kaolinite, Fe, and Al inhibited metolachlor adsorption and transport. In spite of the increased transportability of metolachlor by the presence of soil colloids, the colloid bound herbicide portion accounted for a very small part of the observed increase. This suggests that surface site exclusion mechanisms and preferential sorption induced by the presence of colloids are more important than ion exchange phenomena in promoting herbicide mobility in subsurface environments.     

Field-scale colloid migration experiments in a granite fracture

An understanding of particle migration in fractured rock, required to assess the potential for colloid-facilitated transport of radionuclides, can best be evaluated when the results of laboratory experiments are demonstrated in the field. Field-scale migration experiments with silica colloids were carried out at AECL's Underground Research Laboratory (URL), located in southern Manitoba, to develop the methodology for large-scale migration experiments and to determine whether colloid transport is possible over distances up to 17 m. In addition, these experiments were designed to evaluate the effects of flow rate and flow path geometry, and to determine whether colloid tracers could be used to provide additional information on subsurface transport to that provided by conservative tracers alone. The colloid migration studies were carried out as part of AECL's Transport Properties in Highly Fractured Rock Experiment, the objective of which was to develop and demonstrate methods for evaluating the solute transport characteristics of zones of highly fractured rock. The experiments were carried out within fracture zone 2 as two-well recirculating, two-well non-recirculating, and convergent flow tests, using injection rates of 5 and 10 1 min⁻¹. Silica colloids with a 20 nm size were used because they are potentially mobile due to their stability, small size and negative surface charge. The shapes of elution profiles for colloids and conservative tracers were similar, demonstrating that colloids can migrate over distances of 17 m. The local region of drawdown towards the URL shaft affected colloid migration and, to a lesser extent, conservative tracer migration within the flow field established by the two-well tracer tests. These results indicate that stable colloids, with sizes as small as 20 nm, have different migration properties from dissolved conservative tracers. - 1997 Atomic Energy of Canada.     

Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

Numerical modeling of humic colloid borne americium (III) migration in column experiments using the transport/speciation code K1D and the KICAM model

The humic colloid borne Am(III) transport was investigated in column experiments for Gorleben groundwater/sand systems. It was found that the interaction of Am with humic colloids is kinetically controlled, which strongly influences the migration behavior of Am(III). These kinetic effects have to be taken into account for transport/speciation modeling. The kinetically controlled availability model (KICAM) was developed to describe actinide sorption and transport in laboratory batch and column experiments. Application of the KICAM requires a chemical transport/speciation code, which simultaneously models both kinetically controlled processes and equilibrium reactions. Therefore, the code K1D was developed as a flexible research code that allows the inclusion of kinetic data in addition to transport features and chemical equilibrium. This paper presents the verification of K1D and its application to model column experiments investigating unimpeded humic colloid borne Am migration. Parmeters for reactive transport simulations were determined for a Gorleben groundwater system of high humic colloid concentration (GoHy 2227). A single set of parameters was used to model a series of column experiments. Model results correspond well to experimental data for the unretarded humic borne Am breakthrough.     

Actinide geochemistry: From the molecular level to the real system

Geochemical processes leading to either mobilization or retention of radionuclides in an aquifer system are significantly influenced by their interaction with rock, sediment and colloid surfaces. Therefore, a sound safety assessment of nuclear waste disposal requires the elucidation and quantification of those processes. State-of-the-art analytical techniques as e.g. laser- and X-ray spectroscopy are increasingly applied to study solid–liquid interface reactions to obtain molecular level speciation insight.We have studied the sorption of trivalent lanthanides and actinides onto aluminium oxides, hydroxides and purified clay minerals by the time-resolved laser fluorescence spectroscopy and X-ray-absorption spectroscopy. Chemical constitution and structure of surface bound actinides are proposed based on spectroscopic information. Open questions still remain with regard to the exact nature of mineral surface ligands and the mineral/water interface. Similarities of spectroscopic data obtained for M(III) sorbed onto γ-alumina, and clay minerals suggest the formation of very comparable inner-sphere surface complexes such as S–O–An(III)(OH)_x^(2 − x)(H₂O)_5 − x at pH > 5. Those speciation data are found consistent with those predicted by surface complexation modelling. The applicability of data obtained for pure mineral phases to actinide sorption onto heterogeneously composed natural clay rock is examined by experiments and by geochemical modelling. Good agreement of experiment and model calculations is found for U(VI) and trivalent actinide/lanthanide sorption to natural clay rock. The agreement of spectroscopy, geochemical modelling and batch experiments with natural rock samples and purified minerals increases the reliability in model predictions.The assessment of colloid borne actinide migration observed in various laboratory and field studies calls for detailed information on actinide–colloid interaction. Kinetic stabilization of colloid bound actinides can be due to inclusion into inorganic colloid matrix or by macromolecular rearrangement in case of organic, humic/fulvic like colloids. Only a combination of spectroscopy, microscopy and classical batch sorption experiments can help to elucidate the actinide–colloid interaction mechanisms and thus contribute to the assessment of colloids for radionuclide migration.     

Transport of a non-volatile contaminant in unsaturated porous media in the presence of colloids

Stable colloidal particles can travel long distances in subsurface environments and carry particle-reactive contaminants with them to locations further than predicted by the conventional advective-dispersive transport equation. When such carriers exist in a saturated porous medium, the system can be idealized as consisting of three phases: an aqueous phase, a carrier phase, and a stationary solid matrix phase. However, when colloids are present in an unsaturated porous medium, the system representation should include one more phase, i.e. the air phase. In the work reported, a mathematical model was developed to describe the transport and fate of the colloidal particles and a non-volatile contaminant in unsaturated porous media. The model is based on mass balance equations in a four-phase porous medium. Colloid mass transfer mechanisms among aqueous, solid matrix, and air phases, and contaminant mass transfer between aqueous and colloid phases are represented by kinetic expressions. Governing equations are non-dimensionalized and solved to investigate colloid and contaminant transport in an unsaturated porous medium. A sensitivity analysis of the transport model was utilized to assess the effects of several parameters on model behavior. The colloid transport model matches successfully with experimental data of Wan and Wilson. The presence of air-water interface retards the colloid transport significantly counterbalancing the facilitating effect of colloids. However, the retardation of contaminant transport by colloids is highly dependent on the properties of the contaminant and the colloidal surface.     

Estimation of apparent rate coefficients for phenanthrene and pentachlorophenol interacting with sediments

To gain information on organic pollutants in water-sediment systems, a compartment model was applied to study the sorption course of phenanthrene and pentachlorophenol (PCP) in sediments. The model described the time-dependent interaction of phenanthrene and PCP with operationally defined reversible and irreversible (or slowly reversible) sediment fractions. The interactions between these fractions were described using first order differential equations. By fitting the models to the experimental data, apparent rate constants were obtained using numerical optimization software. The model optimizations showed that the amount of reversible phase increased rapidly in the first 10 d with the sorption time, then decreased after 10 d, while the amount of irreversible phase increased in the total sorption course. That suggested the mass transport between reversible phase and irreversible phase. The extraction efficiency with hot methanol ranged from 36% to 103% for phenanthrene and from 65% to 101% for PCP, with the trend of decreasing with sorption time.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Analysis of colloid and tracer breakthrough curves

We consider the dispersion and elution of colloids and dissolved nonsorbing tracers within saturated heterogeneous porous media. Since flow path geometry in natural systems is often ill-characterized macroscopic (mean) flow rates and dispersion tensors are utilized in order to account for the sub-model scale microscopic fluctuations in media structure (and the consequent hydrodynamic profile). Even for tracer migration and dispersal this issue is far from settled.Here we consider how colloid and tracer migration phenomena can be treated consistently. Theoretical calculations for model flow geometries yield two quantitative predictions for the transport of free (not yet captured) colloids with reference to a non-sorbing dissolved tracer within the same medium: the average migration velocity of the free colloids is higher than that of the tracer; and that the ratio of the equivalent hydrodynamic dispersion rates of colloids and tracer is dependent only upon properties of the colloids and the porous medium, it is independent of pathlengths and fluid flux, once length scales are large enough.The first of these is well known, since even in simple flow paths free colloids must stay more centre stream. The second, if validated suggests how solute and colloid dispersion may be dealt with consistently in macroscopic migration models. This is crucial since dispersion is usually ill-characterized and unaddressed by the experimental literature. In this paper we present evidence based upon an existing Drigg field injection test for the validity of these predictions.We show that starting from experimental data the fitted dispersion rates of both colloids and non-sorbing tracers increase with the measured elution rates (obeying slightly different rules for tracers and colloids); and that the ratio of colloid and nonsorbing tracer elution rates, and the ratio of colloid and nonsorbing tracer dispersion rates may be dependent upon properties of the colloids and the medium (not the flow regime).It is important to realize that even for unretarded species, an earlier peak in the breakthrough curve does not necessarily correspond to a faster mean elution rate, or vice versa. But rather that a colloid may elute faster but disperse less than an equivalent tracer. Hence its peak may be retarded compared to that of the tracer, even assuming no retardation. Hence one must consider a combination of mean elution rate and mean dispersion rate, and not rely on “peak times” to corroborate chromatographic effects. The importance of this lies in the fact that these processes are not independent and yet upscale differently. Thus realistic estimates of effective colloid dispersion rates should be upscaled in a way consistent with that adopted for tracers within the same system.     

The concentration, apparent molecular weight and chemical reactivity of silica from groundwater in southern Nevada

Sorption of radionuclides, metals and organic compounds to colloidal particles has been suggested to increase the mobility of these pollutants in groundwater. Because silicates and alumino-silicates can be important components of groundwater colloids, we have conducted a study to characterize the nature of silica in various springs and wells in Southern Nevada and to determine the extent that silica may be associated with colloidal particles that can participate in pollutant transport. The total silica content was measured using inductively coupled plasma emission spectroscopy (ICP). In addition, reactive silica was measured using the silica molybdate colorimetric technique. The apparent molecular weight of the silica was investigated using split-flow-lateral-transport-thin-cell-fractionation (SPLITT) which can readily distinguish between colloidal and low molecular weight associations. This study indicates that silica does not tend to form stable inorganic colloids in Southern Nevada groundwaters but exists as low molecular weight species. However, water from one of the test facilities on the Nevada Test Site (NTS) did contain stable siliceous colloids that could have important implications for the modeling the transport of radionuclides at this site.     

Influence of the colloid type on the transfer of 60Co and 85Sr in silica sand column under varying physicochemical conditions

The influence of two types of colloids (natural organic matter, NOM), a colloid with high affinity for radionuclides (RN(s)), and hydrophilic synthetic latex (SHL), a colloid with low affinity for RN(s) on the transfer of (60)Co and (85)Sr in a silica sand column was studied under different physicochemical conditions: pH (4.9), ionic strength (10(-3) M and 10(-2) M), concentration of colloids (100 mg l(-1), 10 mg l(-1)), flow velocity (12.4 cm h(-1) and 3.7 cm h(-1)), water saturation of the column (100% and 70%). In the absence of colloids, the transfer of (60)Co and (85)Sr was retarded compared to the transfer of the conservative tracer. In the presence of colloids and according to the specific physicochemical conditions, an acceleration or retardation of (60)Co and (85)Sr transfer was observed compared to their transfer in the absence of colloids. Our results evidenced that any colloids even with low reactivity could significantly modify the RN transfer. However, the extent to which the transfer was influenced differs according to the colloid type; the NOM exhibiting higher impact than SHL. Batch experiments helped in interpreting of the interactions between the colloids, RN(s) and solid phase observed in column.     

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.     

A novel two-dimensional model for colloid transport in physically and geochemically heterogeneous porous media

A two-dimensional model for colloid transport in geochemically and physically heterogeneous porous media is presented. The model considers patchwise geochemical heterogeneity, which is suitable to describe the chemical variability of many surficial aquifers with ferric oxyhydroxide-coated porous matrix, as well as spatial variability of hydraulic conductivity, which results in heterogeneous flow field. The model is comprised of a transient fluid flow equation, a transient colloid transport equation, and an equation for the dynamics of colloid deposition and release. Numerical simulations were carried out with the model to investigate the colloid transport behavior in layered and randomly heterogeneous porous media. Results demonstrate that physical and geochemical heterogeneities markedly affect the colloid transport behavior. Layered physical or geochemical heterogeneity can result in distinct preferential flow paths of colloidal particles. Furthermore, the combined effect of layered physical and geochemical heterogeneity may result in enhanced or reduced preferential flow of colloids. Random distribution of physical heterogeneity (hydraulic conductivity) results in a random flow field and an irregularly distributed colloid concentration profile in the porous medium. Contrary to random physical heterogeneity, the effect of random patchwise geochemical heterogeneity on colloid transport behavior is not significant. It is mostly the mean value of geochemical heterogeneity rather than its distribution that governs the colloid transport behavior.     

Comparative kinetic desorption of 60Co, 85Sr and 134Cs from a contaminated natural silica sand column: influence of varying physicochemical conditions and dissolved organic matter

In order to determine the mechanisms of the retention of 60Co, 85Sr and 134Cs in natural silica sand columns, desorption experiments were performed by changes of pH and ionic strength and by injection of natural organic matter (NOM). Injection of KCl (0.1 M) resulted in a high release of 60Co (60-100%) and 85Sr (72-100%) but a smaller release of 134Cs (31-66%). Only limited release of 60Co (66%) and 85Sr (71%) and no release of 134Cs were observed by injection of NOM. The different percentages of desorption were related to the chemical characteristics of the organic colloids previously retained in columns before the desorption step. The results evidenced different sorption processes on energetically heterogeneous surface sites. According to the initial conditions, the binding of the radionuclides to the solid phase resulted from weak and easily reversible sorption processes to strong association probably by inner sphere complexes. The rather weak release of 134Cs by KCl was attributed to the strong retention of 134Cs by clay coatings on the natural silica sand surfaces.