Diffusion of HTO, Cl and I in Opalinus Clay samples from Mont Terri: Effect of confining pressure

Diffusion coefficients (T=23±2 °C) and accessible porosities for HTO, ³⁶Cl⁻ and ¹²⁵I⁻ were measured on Opalinus Clay (OPA) samples from the Mont Terri Underground Rock Laboratory (URL) using the through-diffusion technique. The direction of transport (diffusion) was perpendicular to bedding. Special cells that allowed the application of confining pressure were designed and constructed. The pressures ranged from 1 to 5 MPa, the latter value simulating the overburden at the Mont Terri URL (about 200 m). The test solution used in the experiments was a synthetic version of the Opalinus Clay pore water, which has Na⁺ and Cl⁻ as the main components (I=0.42 M).The measured values of the effective diffusion coefficients (D_e) and rock capacity factors (α) are: D_e=1.2–1.5×10⁻¹¹ m² s⁻¹ and α=0.09–0.11 for HTO, D_e=4.0–5.5×10⁻¹² m² s⁻¹ and α=0.05 for ³⁶Cl⁻ and D_e=3.2–4.6×10⁻¹² m² s⁻¹ and α=0.07–0.10 for ¹²⁵I⁻. For non-sorbing tracers (HTO, ³⁶Cl) the rock capacity factor α is equal to the diffusion-accessible porosity . The experimental results showed that pressure only had a small effect on the value of the diffusion coefficients. Increasing the pressure from 1 to 5 MPa resulted in a decrease of the diffusion coefficient of 17% for HTO, 28% for ³⁶Cl⁻ and 30% for ¹²⁵I⁻. Moreover, the diffusion coefficients for ³⁶Cl⁻ and ¹²⁵I⁻ are smaller than for HTO, which is consistent with an effect arising from anion exclusion.The diffusion coefficients of HTO and ¹²⁵I⁻ measured in this study are in good agreement with recent measurements at three other laboratories performed within the framework of a laboratory comparison exercise. The values of the diffusion-accessible porosities show a larger degree of scatter.     

Diffusion of tritiated water and Na through non-degraded hardened cement pastes

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and ²²Na⁺, considered to be conservative tracers, have been carried out in triplicates in a glove box under a controlled nitrogen atmosphere. Each experiment consisted of a through-diffusion test followed by an out-diffusion test.The experimental data were inversely modelled applying an automated Marquardt–Levenberg procedure. The analysis of the through-diffusion data allowed the extraction of values for the effective diffusion coefficients, D_e, and the rock capacity factor, α. Good agreement between measured and calculated tracer breakthrough curves was achieved using both a simple diffusion model without sorption and a diffusion/linear sorption model. The best-fit K_d-values were found to be consistent with R_d-values measured in previous batch-sorption experiments.The best-fit values from the through-diffusion tests were then used to predict the results of subsequent out-diffusion experiments. Good agreement between experimental data and predictions was achieved only for the case of linear sorption.Isotopic exchange can only partially account for both the amount of tracer taken up in the batch-sorption tests and the measured retardation in the diffusion experiments and, hence, additional mechanisms have to be invoked to explain the data.     

Modelling of silica diffusion experiments with Si in Boom Clay

A mathematical model describing the dissolution of nuclear glass directly disposed in clay combines a first-order dissolution rate law with the diffusion of dissolved silica in clay. According to this model, the main parameters describing the long-term dissolution of the glass are ηR, the product of the diffusion accessible porosity η and the retardation factor R, and the apparent diffusion coefficient D_app of dissolved silica in clay.For determining the migration parameters needed for long-term predictions, four Through-Diffusion (T-D) experiments and one percolation test have been performed on undisturbed clay cores. In the Through-Diffusion experiments, the concentration decrease after injection of ³²Si (radioactive labelled silica) was measured in the inlet compartment. At the end of the T-D experiments, the clay cores were cut in thin slices and the activity of labelled silica in each slice was determined. The measured activity profiles for these four clay cores are well reproducible.Since no labelled silica could be detected in the outlet compartments, the Through-Diffusion experiments are fitted by two In-Diffusion models: one model assuming linear and reversible sorption equilibrium and a second model taking into account sorption kinetics. Although the kinetic model provides better fits, due to the sufficiently long duration of the experiments, both models give approximately similar values for the fit parameters. The single percolation test leads to an apparent diffusion coefficient value about two to three times lower than those of the Through-Diffusion tests.Therefore, dissolved silica appears to be strongly retarded in Boom Clay. A retardation factor R between 100 and 300 was determined. The corresponding in situ distribution coefficient K_d is in the range 25–75 cm³ g⁻¹. The apparent diffusion coefficient of dissolved silica in Boom Clay is estimated between 2×10⁻¹³ and 7×10⁻¹³ m² s⁻¹. The pore diffusion coefficient is in the range from 6×10⁻¹¹ to 1×10⁻¹⁰ m² s⁻¹.     

Effective diffusivity of the uranyl ion in a granite from Inada, Ibaraki, Japan

The effective diffusivity of uranium(VI) in Inada granite has been determined by through-diffusion. Experiments were performed at room temperature (20–25°C) in a 0.1 mol 1⁻¹ KCl solution where uranium is present predominantly as the poorly sorbing UO₂²⁺. An effective diffusivity (D_e) of (3.6 ± 1.6) × 10⁻¹⁴ m² s⁻¹ was obtained, close to that for uranine (nonsorbing organic tracer), but one order of magnitude lower than those obtained for Sr²⁺ and NpO₂⁺, and two orders of magnitude lower than that obtained for I⁻. According to well established theory, a proportional relationship exists between D_e and the diffusivity in the bulk of the solution (D_v). The effective diffusivity obtained in granite was not proportional to D_v. This agrees with results obtained for effective diffusivity in a Swedish granite. The ratio D_e/D_v was found to be not constant but increased with D_e or D_v. This result suggests a limit to the application of the theory.     

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.     

Consistent interpretation of the results of through-, out-diffusion and tracer profile analysis for trace anion diffusion in compacted montmorillonite

Literature data for anion diffusion in compacted swelling clays contain systematic inconsistencies when the results of through-diffusion tests are compared with those of out-diffusion or tracer profile analysis. In the present work we investigated whether these inconsistencies can be explained by taking into account heterogeneities in the compacted samples; in particular increased porosities at the clay boundaries. Based on the combined results of out-diffusion, tracer profile analysis and the spatial distribution of the electrolyte anion in the clay, we conclude that the inconsistencies can indeed be resolved by taking into account a heterogeneous distribution of the total and the anion-accessible porosity. This, by definition, leads to a position dependence of the effective diffusion coefficient. Neglecting these effects results in a rather subordinate systematic error in the determination of effective diffusion coefficients of anions from through-diffusion tests with clay thicknesses in the centimetre range. However, stronger errors in terms of absolute values and conceptual interpretation may be introduced in out-diffusion tests and profile analyses of the diffused tracer. We recommend that anion diffusion tests should be accompanied by measurements of the total and anion-accessible porosity as a function of position in the direction of diffusion.     

The assessment of electromigration as a new technique to study diffusion of radionuclides in clayey soils

For the performance assessment study of a geological disposal of High Level Waste (HLW) in clayey formations, migration studies are essential. For low permeability soils (clays), classical diffusion studies take a very long time. In order to reduce the experimental time, we propose an electrical field as driving force to accelerate the migration of ionic species. This paper reports the assessment of the electromigration technique as a powerful new and fast technique for migration studies. The apparent molecular diffusion coefficient can be derived by two independent methods using the migration parameters obtained from an electromigration experiment, namely the apparent dispersion coefficient and the apparent convection velocity. First, it can be calculated from the velocity of the migrating species by the Einstein relation. But, corrections are necessary for electroosmotic flow. The apparent electroosmotic mobility is experimentally determined as 2.2·10⁻⁹ m²/Vs. Second, it can be calculated from the relation between the apparent dispersion coefficient and the total apparent convection velocity. But it is necessary to know the dispersion length of the medium. The dispersion length for Boom Clay is experimentally determined as 8·10⁻⁵ m. Because of the serious reduction in time, it becomes possible to run series of experiments at different electrical fields to obtain averaged values for the apparent molecular diffusion coefficient according to the two methods. Experiments at different electrical fields have another advantage: the intercept of the linear relationship between the total apparent convection velocity and the apparent dispersion coefficient gives the apparent molecular diffusion coefficient. The apparent molecular diffusion coefficients obtained for ⁸⁵Sr, ¹³¹I and HTO are respectively 0.8·10⁻¹¹, 15·10⁻¹¹, and 24·10⁻¹¹ m²/s. These values are confirmed by pure diffusion experiments. The excellent agreement with the apparent molecular diffusion coefficients obtained by classical diffusion tests clearly demonstrates the feasibility of the electromigration technique for the determination of diffusion coefficients.     

Kinetics of OH-initiated oxidation of oxygenated organic compounds in the aqueous phase: new rate constants, structure–activity relationships and atmospheric implications

The kinetics of OH oxidation of several organic compounds of atmospheric relevance were measured in the aqueous phase. Relative kinetics were performed using various organic references and OH sources. After validation of the protocol, temperature-dependent rate constants for the reactions of OH radical with ethyl ter-butyl ether (, E_a/R=580 (±560) K), n-butyl acetate ( (±0.4)×10⁹ M⁻¹ s⁻¹, E_a/R=1000 (±200) K), acetone ( (±0.05)×10⁹ M⁻¹ s⁻¹, E_a/R=1400 (±500) K), methyl ethyl ketone (, E_a/R=1200 (±200) K), methyl iso-butyl ketone (, E_a/R=1200 (±300) K) and methylglyoxal (, E_a/R=1100 (±300) K) were determined. A non-Arrhenius behavior was found for phenol, in good agreement with the contribution of an OH addition to the mechanism, which also includes H-abstraction by OH radicals. Global rate constants of acetaldehyde, propionaldehyde, butyraldehyde and valeraldehyde were studied at 298 K only, as these compounds partly hydrate in the aqueous phase. All the obtained data (except those of phenol) complemented by literature data were used to investigate three methods to estimate rate constants for H-abstraction reactions of OH radicals in aqueous solutions when measured data were not available: Evans-Polanyi-type correlations, comparisons with gas-phase data, structure activity relationships (SAR). The results show that the SAR method is promising; however, the data set is currently too small to extend this method to temperatures other than 298 K. The atmospheric impact of aqueous phase OH oxidation of water-soluble organic compounds is discussed with the determination of their global atmospheric lifetimes, taking into account both gas- and aqueous-phase reactivities. The results show that atmospheric droplets can act as powerful photoreactors to eliminate soluble organic compounds from the atmosphere.     

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish Äspö Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing ²²Na⁺, ⁸⁵Sr²⁺, ⁴⁷Ca²⁺and more strongly sorbing ⁸⁶Rb⁺, ¹³³Ba²⁺, ¹³⁷Cs⁺.Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity.The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations.After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted K_d values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning K_ds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.     

Measurement of total OH reactivity by laser-induced pump and probe technique—comprehensive observations in the urban atmosphere of Tokyo

Total OH reactivity was observed by use of the laser-induced pump and probe technique, and the urban air quality in Tokyo was diagnosed comprehensively. The concentrations of NO_x, CO, O₃, non-methane hydrocarbons (NMHCs) and oxygenated volatile organic compounds (OVOCs) were observed simultaneously. The observations were conducted in July and August 2003, and in January, February, May, and November 2004. Generally, the observed OH reactivity was higher than the calculated values derived using the observed concentrations of the trace species. The differences between the observed and calculated values in summer, spring, and autumn were approximately 30%. However, the difference in winter was smaller than those in the other seasons. In addition, while the differences observed in summer, spring, and autumn correlated with the total reactivity of the OVOCs (Σ_i k_OVOCi[OVOC_i](s⁻¹), k_i is rate constant of its compounds with OH), the correlations were not confirmed in the case of winter because atmospheric oxidation was less active and OVOCs levels were low in winter. These results suggest that the secondary products of the photochemical reactions in the atmosphere would be a missing sink for the OH loss process in the urban area.     

Influence of cosolvents on quinoline sorption by subsurface materials and clays

Quinoline sorption was measured on Na-saturated subsurface materials, a natural clay isolate, and montmorillonite, and was dominated by exchange of the quinolinium ion. In water/methanol and water/acetone mixtures, quinoline sorption on the subsoils and clays was lower than from water. In cosolvent, sorption followed the ionization fraction indicating the continued predominance of ion exchange. The reduction in quinoline sorption by cosolvent was similar for all the subsoils and clays indicating commonality in the surface-solute-solvent interaction. Conditional equilibrium constants (_cK_ex) for quinoline exchange on the subsoils in water/methanol mixtures decreased log-linearly with mole percent cosolvent up to 20% methanol. This decrease closely followed the increase in quinoline solubility in the cosolvent mixtures. Acetone caused greater reduction in sorption than methanol, at comparable mole percent, in accordance with its lower dielectric constant and enhanced solvating power. A generalized thermodynamic approach based on the concept of transfer activity coefficients was developed to account for the cosolvent effect on _cK_ex, and was successfully applied to the quinoline sorption data. The thermodynamic analysis suggested that enhanced solvation of the organic cation in the bulk solvent and desolvation of Na⁺ at the charged surface predominate the cosolvent effect.     

An evaluation of contaminant migration patterns at two waste disposal sites on fractured porous media in terms of the equivalent porous medium (EPM) model

Contamination has occurred many non-indurated and bedrock systems wherein the groundwater flows almost exclusively through a network of connected, open fractures. The matrix surrounding the fractures often possesses porosity which allows contaminant diffusion into the matrix. If the diffusion rates are fast relative to the fracture groundwater velocity, transport effects may be predicted by considering the system to be an equivalent porous medium (EPM). The rapidity with which fracture/immobile-matrix equilibrium is established will be determined in part by the: fracture aperture (2b); interfracture spacing (2B); porosity in the immobile matrix (θ_im); and the matrix diffusion coefficient (D′). Two systems which are characterized by very different values of the above parameters have been studied by our laboratories. At Alkali Lake, Oregon, the EPM approach describes contaminant transport well. At Bayview Park, Ontario, the EPM approach is not appropriate. Several features of the two sites are compared to illustrate the different nature of these two sites. These features include: (1) natural characteristics of the groundwater systems; (2) contaminant distributions; (3) observed transport; and (4) computed fracture/immobile-matrix diffusion times.     

The use of clays to sequestrate organic pollutants. Leaching experiments

Leaching experiments are performed from clay-pollutant systems in order to evaluate the capability of clays to sequestrate organic pollutants from wastewaters. Reference kaolinite KGa-1b, montmorrillonite SWy-2 and reference soil BCR^®-700 are the sorbent materials. 2,4,6-trichloroaniline (2,4,6-TCA) and 4-chlorophenol (4-CP) are the typical pollutants, sorbed at amounts of 10.0 mg g⁻¹ and 5.8 mg g⁻¹ on SWy-2 and 7.3 mg g⁻¹ and 2.2 mg g⁻¹ on KGa-1b, respectively.The leaching agents are ultrapure water and model solutions of acid rain and surface waters that simulate meteoric leaching. 1.0 mM HNO₃, 1.0 mM H₂SO₄ solutions and a methanol/water 50/50 (v/v) mixture simulate leaching agents of industrial source.The results are compared and the preferential capability of the clays to sequestrate the more lipophilic 2,4,6-TCA is evidenced.The bond interactions are discussed and explained through preferential adsorption reactions. For montmorrillonite also a simultaneous intercalation in the phyllosilicate interlayer is proposed.     

Dissolution behavior of representative elements from red mud (RM) by leaching with titanium white waste acid (TWWA)

Titanium white waste acid (TWWA) was used to dissolve the representative elements from red mud (RM) to achieve the goal of “treating waste with waste.” The leaching parameters on the leaching efficiency of Na, Sc, and Al were investigated, in which the analysis of XRD and SEM–EDS on RM and leaching residue was performed. The leaching kinetics of Na, Sc, and Al was studied with unreacted shrinking core model (USCM). The results show that the dealkalization efficiency was close to 100%, and the leaching efficiency of Sc and Al was 82% and 75%, respectively. Cancrinite was dissolved from RM, and then the elements such as Na, Al, and Ca reacted with H₂SO₄ of TWWA. Na existed in the leaching liquor in the form of ions. Ca reacted with sulfuric acid to form anhydrite, which existed in the leaching residue. The particles of RM became smaller and dispersed with each other by acid leaching. The leaching apparent activation energy of Na, Sc, and Al was 4.947 kJ/mol, 6.361 kJ/mol, and 31.666 kJ/mol, respectively. The leaching kinetic equation of Na, Sc, and Al was 1???(1???a)^2/3?=?0.084·exp[??595.05/T]·t by external diffusion, 1???2a/3???(1???a)^2/3?=?0.021·exp[??765.16/T]·t by internal diffusion, and ln(1???a)/3?+?(1???a)^?2/3 – 1?=?67.12·exp[??3808.8/T]·t by joint action, respectively.

     

Heavy-metal ion complexation by particulate matter in the leachate of solid waste: a multi-method approach

The metal ion binding characteristics of particulate matter obtained from column experiments on the anaerobic digestion of solid waste were studied using a titrimetric approach. The experimental set-up allowed us to study the dynamics of particle bound ligand concentrations during digestion processes typically found in landfills.We developed a continuous titration method by simultaneously using a Cd-sensitive and pH electrode and combining metal and acid/base titrations. This technique allows for a more precise determination of pK_a-log K_M pairs for each ligand than metal titrations alone. The results were compared with titration methods using differential pulse anodic stripping voltammetry (DPASV) and atomic absorption spectroscopy (AAS) with longer equilibration times in order to further characterize ligand properties such as reaction kinetics, the electrochemical lability of the respective complex during DPASV, the distinction between metal adsorption to particulate matter and metal complexation by soluble ligands adhered to particles, reversibility of the binding process by competition studies, and resistance against purging with nitrogen gas.The properties of seven major metal binding ligands were identified and assignments to the most likely functional groups were made. The most important ligand properties are for ligand A: pK_a ≈ 9.2, log K_cd ≈ 7.0 fast reaction kinetics (mercapto groups); ligand B: pK_a = 4.8, log K_Cd ≈ 6.0, slow reaction kinetics (chelates with 3 or 4 carboxylic groups); ligand C: pK_a ≈ 6.0, log K_Cd ≈ 13.0, irreversible metal binding at basic pH-values (uptake inside bacterial cells); ligand D: pK_a = 7.7, log K_Cd = 4.0, runs parallel to N content of particulate matter with digestion time (primary amines neighboring oxo groups); ligand E: pK_a ≈ 12.0, log K_Cd = 9.0, runs parallel to P content of particulate matter (phosphate); ligand F:pK_a > 9.0, log K_{Cdf = pKa + 0.4}, runs parallel to N content of particulate matter (primary amines neighboring SH groups); and ligand G: pK_a ≤ 4.8, log K_Pb ≈ 4.3, strong Pb²⁺ ligand, even at low pH-values.Metal ions were found to be irreversibly bound by ligand C at low heavy-metal concentratins, whereas at higher concentrations the binding is reversible and can be predicted using the mass of the digestion process (methanogenic phase). All other ligands have their concentration maximum in the transition phase between acetogenic and methanogenic phase.     

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak ‘starter paths’ formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N₂ gas pressure gradient (ΔPN₂) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock.     

Partitioning of organochlorine pesticides from water to polyethylene passive samplers

The mass transfer rates and equilibrium partitioning behaviour of 14 diverse organochlorine pesticides (OCP) between water and polyethylene (PE) passive samplers, cut from custom made PE sheets and commercial polyethylene plastic bags, were quantified. Overall mass transfer coefficients, k_O, estimated PE membrane diffusion coefficients, D_PE, and PE-water partitioning coefficients, K_PE-water, are reported. In addition, the partitioning of three polycyclic aromatic hydrocarbons (PAHs) from water to PE is quantified and compared with literature values. K_PE-water values agreed mostly within a factor of two for both passive samplers and also with literature values for the reference PAHs. As PE is expected to exhibit similar sorption behaviour to long-chain alkanes, PE-water partitioning coefficients were compared to hexadecane-water partitioning coefficients estimated with the SPARC online calculator, COSMOtherm and a polyparameter linear free energy relationship based on the Abraham approach. The best correlation for all compounds tested was with COSMOtherm estimated hexadecane-water partitioning coefficients.     

Impact of initial and boundary conditions on preferential flow

Preferential flow in soil is approached by a water-content wave, WCW, that proceeds downward from the ground surface. WCWs were obtained from sprinkler experiments with infiltration rates varying from 5 to 40 mm h^− 1. TDR-probes and tensiometers measured volumetric water contents θ(z,t) at seven depths, and capillary heads, h(z,t) at six depths in a column of an undisturbed soil. The wave is characterized by the velocity of the wetting front, c_W, the amplitude, w_S, and the final water content, θ. We tested with uni-variate and bi-variate linear regressions the impacts of initial volumetric water contents, θ_ini, and input rates, q_S, on c_W, w_S and θ.The test showed that θ_ini influenced θ and w_S and q_S effected c_W. The expected proportionality of w_S ≈ qs^1/3 was weak and c_W ≈ qs^2/3 was strong.     

Flux density and breakthrough times for water and tracer in a spatially variable, compacted clay soil

Flux density values, computed from observed infiltration and outflow measurements at 184 locations in a 0.3-m-thick, 9m × 23m layer of compacted clay subsoil, are compared to effective flux density values that are based on breakthrough time distributions for water and Br⁻ tracer over the same area. Results suggest that both water and tracer move at similar rates, but considerably faster than expected, on the basis of flux density alone, and that only a small fraction of the total pore space is involved in active transport. The ramifications of these findings are explored against the background of effective porosity, degree of compaction, and observed changes in bulk density with time.