Diffusion of HTO, 36Cl- and 125I- in Opalinus Clay samples from Mont Terri. Effect of confining pressure

Diffusion coefficients (T=23 +/- 2 degrees C) and accessible porosities for HTO, 36Cl(-) and 125I(-) 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 (alpha) are: D(e)=1.2-1.5 x 10(-11) m(2) s(-1) and alpha=0.09-0.11 for HTO, D(e)=4.0-5.5 x 10(-12) m(2) s(-1) and alpha=0.05 for 36Cl(-) and D(e)=3.2-4.6 x 10(-12) m(2) s(-1) and alpha=0.07-0.10 for 125I(-). For non-sorbing tracers (HTO, 36Cl) the rock capacity factor alpha is equal to the diffusion-accessible porosity epsilon. 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 approximately 17% for HTO, approximately 28% for 36Cl(-) and approximately 30% for 125I(-). Moreover, the diffusion coefficients for 36Cl(-) and 125I(-) are smaller than for HTO, which is consistent with an effect arising from anion exclusion. The diffusion coefficients of HTO and 125I(-) 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.     

Through-electromigration: a new method of investigating pore connectivity and obtaining formation factors

The retardation of radionuclides and other contaminants in fractured crystalline rock is strongly associated with the diffusive properties of the rock matrix. At present, the scientific community is divided concerning the question of long-range pore connectivity in intrusive igneous rock. This paper presents a fast new method, called the through-electromigration method, of obtaining formation factors and investigating pore connectivity. The method involves the migration of an ionic tracer through a rock sample with an electrical potential gradient as the main driving force. The method is analogous to the through-diffusion method but the experimental time is reduced by orders of magnitude. This enables investigations of pore connectivity, as measurements can be made on longer samples. In a preliminary investigation, the new method is compared to the traditional through-diffusion method as well as to rock resistivity methods. The diffusive properties of nine granitic rock samples from Laxemar in Sweden, ranging from 15 to 121 mm in length, have been investigated and the results are compared.     

Effect of temperature on the containment properties of argillaceous rocks: The case study of Callovo-Oxfordian claystones

Heat generated by high level radioactive wastes could alter the performance of a clay repository. It was intended to investigate the effect of such a thermal period on the diffusive properties of Callovo-Oxfordian claystones. Thus, through-diffusion experiments with HTO, Cl-36, Na-22 and Cs-137 were performed before, during and after stages of heating at 80°C that lasted for up to one year. A special attention was paid to limit the occurrence of any chemical disturbance. Therefore (i) the temperature was raised to 80°C, then progressively brought back to 21°C, thanks to three intermediate temperature stages, and (ii) specific synthetic solutions were used for each temperature, chemistry of which being close to the equilibrium state, especially with respect to the carbonate and sulphate minerals. It was found that experiments carried out at 80°C showed a clear increase of the effective diffusion coefficient values for the four tracers with respect to those obtained at 21°C (by a factor of 3 for HTO and Cl-36, 5 for Na-22 and 2 for Cs-137). On the other hand, the porosity and rock capacity values did not exhibit any significant discrepancy between 21°C and 80°C, indicating no observable damage of both the pore conducing network and the sorption properties of clay minerals. The Stokes-Einstein relationship, based on the temperature dependency of the viscosity of bulk water, could be used to describe the temperature dependence of the diffusion of HTO and Cl-36 but failed to describe the diffusive evolution of the two sorbing cations, Na-22 and Cs-137. Furthermore, experiments performed after the thermal period led to diffusive properties well matching those obtained before heating. All these results suggest that at the lab scale the heating of rock samples would not alter the claystone containment properties.     

Matrix diffusion coefficients in volcanic rocks at the Nevada test site: influence of matrix porosity, matrix permeability, and fracture coating minerals

Diffusion cell experiments were conducted to measure nonsorbing solute matrix diffusion coefficients in forty-seven different volcanic rock matrix samples from eight different locations (with multiple depth intervals represented at several locations) at the Nevada Test Site. The solutes used in the experiments included bromide, iodide, pentafluorobenzoate (PFBA), and tritiated water ((3)HHO). The porosity and saturated permeability of most of the diffusion cell samples were measured to evaluate the correlation of these two variables with tracer matrix diffusion coefficients divided by the free-water diffusion coefficient (D(m)/D*). To investigate the influence of fracture coating minerals on matrix diffusion, ten of the diffusion cells represented paired samples from the same depth interval in which one sample contained a fracture surface with mineral coatings and the other sample consisted of only pure matrix. The log of (D(m)/D*) was found to be positively correlated with both the matrix porosity and the log of matrix permeability. A multiple linear regression analysis indicated that both parameters contributed significantly to the regression at the 95% confidence level. However, the log of the matrix diffusion coefficient was more highly-correlated with the log of matrix permeability than with matrix porosity, which suggests that matrix diffusion coefficients, like matrix permeabilities, have a greater dependence on the interconnectedness of matrix porosity than on the matrix porosity itself. The regression equation for the volcanic rocks was found to provide satisfactory predictions of log(D(m)/D*) for other types of rocks with similar ranges of matrix porosity and permeability as the volcanic rocks, but it did a poorer job predicting log(D(m)/D*) for rocks with lower porosities and/or permeabilities. The presence of mineral coatings on fracture walls did not appear to have a significant effect on matrix diffusion in the ten paired diffusion cell experiments.     

Tracer diffusion coefficients in sedimentary rocks: correlation to porosity and hydraulic conductivity

Matrix diffusion is an important transport process in geologic materials of low hydraulic conductivity. For predicting the fate and transport of contaminants, a detailed understanding of the diffusion processes in natural porous media is essential. In this study, diffusive tracer transport (iodide) was investigated in a variety of geologically different limestone and sandstone rocks. Porosity, structural and mineralogical composition, hydraulic conductivity, and other rock properties were determined. The effective diffusion coefficients were measured using the time-lag method. The results of the diffusion experiments indicate that there is a close relationship between total porosity and the effective diffusion coefficient of a rock (analogous to Archie's Law). Consequently, the tortousity factor can be expressed as a function of total porosity. The relationship fits best for thicker samples (> 1.0 cm) with high porosities (> 20%), because of the reduced influence of heterogeneity in larger samples. In general, these correlations appear to be a simple way to determine tortuosity and the effective diffusion coefficient from easy to determine rock porosity values.     

Direct simulation of heterogeneous diffusion and inversion procedure applied to an out-diffusion experiment. Test case of Palmottu granite

An out-diffusion laboratory experiment using a non-reactive tracer was fitted using the Time Domain Diffusion (TDD) method. This rapid particle tracking method allows simulation of the heterogeneous diffusion based on pore-scale images and local values of diffusivities. The superimposed porosity and mineral 2D maps act as computation grids to condition diffusion pathways. We focused on a Palmottu granite sample, in which the connected pore space has a composite microstructure with cracks linking microporous minerals and is above the percolation threshold. Three main results were achieved: (i) When compared to the fitting obtained with one coefficient (best mean square residual R = 1.6 x 10(-2)), diffusion is shown to be suitably characterised with two coefficients related to cracks and microporous minerals (best R = 6.5 x 10(-4)), (ii) rather than imposing a local apparent diffusion coefficient D(a) independent of the local porosity Phi, a best fit is obtained by applying Archie's relationship D(a) = D(0) x G with G = Phi(m) to each pixel of the calculation grids (G is the geometry factor, D(0) is the diffusion coefficient in free fluid, and m is Archie's exponent), and (iii) the order of magnitude of the fitted diffusion coefficient or Archie's exponents (m=0 for microcracks and m=1.82 for microporous minerals) is physically realistic.     

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

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and 22Na(+), 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, alpha. 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.     

Determination of spatially-resolved porosity, tracer distributions and diffusion coefficients in porous media using MRI measurements and numerical simulations

This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials.     

Contribution of the diffusive exchange method to the characterization of pore-water in consolidated argillaceous rocks

Consolidated argillaceous rocks, potential host rocks for nuclear waste disposal, display particular properties (low water content, very small pore size) that require specific methods, especially for the characterization of their pore-water. The diffusive exchange method was applied to the claystone of Tournemire by taking into account the effect of swelling or micro-cracks induced by drilling. The good agreement of diffusion parameters with those obtained from the through-diffusion method showed that the extent of this disturbance was low. This method led to stable isotope content values in pore-water in good agreement with those measured in the fracture fluids located in the vicinity, validating the method. On the contrary, the vacuum distillation technique gave a systematic depletion in heavy isotopes of water, probably induced by an incomplete extraction of pore-water. For the chloride tracer, the determination of the accessible porosity on each sample allowed the estimate of chloride concentration of pore-water to be improved, by taking into account all the petrophysical heterogeneities.     

Diffusive transport through compacted Na- and Ca-bentonite

The effect of exchangeable cation — Na⁺ and Ca ²⁺ — on the diffusive transport of I⁻, Sr ²⁺ and ³H (as HTO) in compacted bentonite was examined using a through-diffusion method. Total intrinsic diffusion coefficients, D_i, were determined from the steady-state flux of the diffusants through the clays, and apparent diffusion coefficients, D_a, were obtained from the time lag technique. The clays were compacted to a dry bulk density of 1.3 Mg/m³, and Na-bentonite was saturated with a solution of 100 mol NaCl/m3 and Ca-bentonite with one of 50 mol CaCl₂/m³. The D_i values for all diffusants are 2 to 6 times higher in the Ca- than Na-clay. We attribute this to the larger quasicrystal, or particle, size of Ca- compared to Na-bentonite. Hence, Ca-bentonite has a greater proportion of relatively large pores; this was confirmed by Hg intrusion porosimetry. This means the diffusion pathways in Ca-bentonite are less tortuous than those in Na-bentonite. Moreover, in some cases the effective porosity, or the porosity available for diffusive transport, may be greater in Ca-bentonite. The D_a values are inversely proportional to the distribution coefficients of the diffusants with the clays.     

Diffusive flux and pore anisotropy in sedimentary rocks

Diffusion of dissolved contaminants into or from bedrock matrices can have a substantial impact on both the extent and longevity of dissolved contaminant plumes. For layered rocks, bedding orientation can have a significant impact on diffusion. A series of laboratory experiments was performed on minimally disturbed bedrock cores to measure the diffusive flux both parallel and normal to mineral bedding of four different anisotropic sedimentary rocks. Measured effective diffusion coefficients ranged from 4.9×10(-8) to 6.5×10(-7)cm(2)/s. Effective diffusion coefficients differed by as great as 10-folds when comparing diffusion normal versus parallel to bedding. Differences in the effective diffusion coefficients corresponded to differences in the "apparent" porosity in the orientation of diffusion (determined by determining the fraction of pore cross-sectional area measured using scanning electron microscopy), with the difference in apparent porosity between normal and parallel bedding orientations differing by greater than 2-folds for two of the rocks studied. Existing empirical models failed to provide accurate predictions of the effective diffusion coefficient in either bedding orientation for all four rock types studied, indicating that substantial uncertainty exists when attempting to predict diffusive flux through sedimentary rocks containing mineral bedding. A modified model based on the apparent porosity of the rocks provided a reasonable prediction of the experimental diffusion data.     

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.     

In situ diffusion experiment in granite: phase I

A program of in situ experiments, supported by laboratory studies, was initiated to study diffusion in sparsely fractured rock (SFR), with a goal of developing an understanding of diffusion processes within intact crystalline rock. Phase I of the in situ diffusion experiment was started in 1996, with the purpose of developing a methodology for estimating diffusion parameter values. Four in situ diffusion experiments, using a conservative iodide tracer, were performed in highly stressed SFR at a depth of 450 m in the Underground Research Laboratory (URL). The experiments, performed over a 2 year period, yielded rock permeability estimates of 2 x 10(-21) m(2) and effective diffusion coefficients varying from 2.1 x 10(-14) to 1.9 x 10(-13) m(2)/s, which were estimated using the MOTIF code. The in situ diffusion profiles reveal a characteristic "dog leg" pattern, with iodide concentrations decreasing rapidly within a centimeter of the open borehole wall. It is hypothesized that this is an artifact of local stress redistribution and creation of a zone of increased constrictivity close to the borehole wall. A comparison of estimated in situ and laboratory diffusivities and permeabilities provides evidence that the physical properties of rock samples removed from high-stress regimes change. As a result of the lessons learnt during Phase I, a Phase II in situ program has been initiated to improve our general understanding of diffusion in SFR.     

Characterization of controlled-release KMnO4 (CRP) barrier system for groundwater remediation: a pilot-scale flow-tank study

Release and spreading of permanganate (MnO(4)(-)) in the well-based controlled-release potassium permanganate (KMnO(4)) barrier system (CRP system) was investigated by conducting column release tests, model simulations, soil oxidant demand (SOD) analyses, and pilot-scale flow-tank experiments. A large flow tank (L x W x D=8m x 4m x 3m) was constructed. Pilot-scale CRP pellets (OD x L=0.05 m x1.5m; n=110) were manufactured by mixing approximately 198 kg of KMnO(4) powders with paraffin wax and silica sands in cylindrical moulds. The CRP system (L x W x D=3m x 4m x 1.5m) comprising 110 delivery wells in three discrete barriers was constructed in the flow tank. Natural sands (organic carbon content=0.18%; SOD=3.7-11 g MnO(4)(-)kg(-1)) were used as porous media. Column release tests and model simulations indicated that the CRP system could continuously release MnO(4)(-) over several years, with slowly decreasing release rates of 2.5 kg d(-1) (day one), 109 g d(-1) (day 100), 58 g d(-1) (year one), 22 g d(-1) (year five), and 12 g d(-1) (year 10). Mean MnO(4)(-) concentrations within the CRP system ranged from 0.5 to 6 mg l(-1) during the 42 days of testing period. The continuously releasing MnO(4)(-) was gradually removed by SOD limiting the length of MnO(4)(-) zone in the porous media. These data suggested that the CRP system could create persistent and confined oxidation zone in the subsurface. Through development of advanced tools for describing agent transport and facilitating lateral agent spreading, the CRP system could provide new approach for long-term in situ treatment of contaminant plumes in groundwater.     

Pentachlorophenol association with fulvic acids from recycled wastes

The apparent water solubility of pentachlorophenol was measured at pH=6 and at 25 degrees C in pure water, aqueous solutions of three salts (NaCl, KNO(3) and CaCl(2) at 0.010, 0.10 and 1.0M) and in aqueous solutions of three fulvic acids samples extracted from a natural soil (sFA), composted sewage sludge (csFA) and composted livestock's material (lsFA). A solubility enhancement method was developed for the measurement of partition coefficients (K(oc), L/kg organic carbon). Pentachlorophenol associates strongly with the fulvic acid samples and the calculated K(oc) were the following (averages and standard deviations): (sFA) (211+/-22) x 10(2), (csFA) (253+/-26) x 10(2), (lsFA) (235+/-10) x 10(2). For comparison purposes the K(oc) for pyrene were also calculated for the three FA samples and were the following: (sFA) (119+/-10) x 10(2), (csFA) (239+/-21) x 10(2), (lsFA) (92+/-10) x 10(2). The analysis of variance (one-way ANOVA) of the effect of the type of FA sample on the solubilization of pentachlorophenol and pyrene shows that this factor causes significant differences on the aqueous solubilization of these two organic substances.     

Measurement of effective air diffusion coefficients for trichloroethene in undisturbed soil cores

In this study, we measure effective diffusion coefficients for trichloroethene in undisturbed soil samples taken from Picatinny Arsenal, New Jersey. The measured effective diffusion coefficients ranged from 0.0053 to 0.0609 cm2/s over a range of air-filled porosity of 0.23-0.49. The experimental data were compared to several previously published relations that predict diffusion coefficients as a function of air-filled porosity and porosity. A multiple linear regression analysis was developed to determine if a modification of the exponents in Millington's [Science 130 (1959) 100] relation would better fit the experimental data. The literature relations appeared to generally underpredict the effective diffusion coefficient for the soil cores studied in this work. Inclusion of a particle-size distribution parameter, d10, did not significantly improve the fit of the linear regression equation. The effective diffusion coefficient and porosity data were used to recalculate estimates of diffusive flux through the subsurface made in a previous study performed at the field site. It was determined that the method of calculation used in the previous study resulted in an underprediction of diffusive flux from the subsurface. We conclude that although Millington's [Science 130 (1959) 100] relation works well to predict effective diffusion coefficients in homogeneous soils with relatively uniform particle-size distributions, it may be inaccurate for many natural soils with heterogeneous structure and/or non-uniform particle-size distributions.     

Modeling solute diffusion in the presence of pore-scale heterogeneity: method development and an application to the Culebra dolomite member of the Rustler Formation, New Mexico, USA

Previous studies have revealed the presence of pore-scale variability in diffusivity in the Culebra (dolomite) member of the Rustler Formation, NM. In this study, eight laboratory-scale diffusion experiments on five Culebra samples were analyzed using a methodology for modeling solute diffusion through porous media in the presence of multiple matrix diffusivities, Dp. A lognormal distribution of Dp is assumed within each of the lab samples. The estimated standard deviation (sigma d) of ln(Dp) within each sample ranges from 0 to 1, with most values lying between 0.5 and 1. The variability over all samples leads to a combined sigma d in the range of 1.0-1.2, which is consistent with the distribution of independently determined formation factor measurements for similar Culebra samples. A comparison of our estimation results to other rock properties suggests that, at the lab-scale, the geometric mean of Dp increases with bulk porosity and the quantity of macroscopic features such as vugs and fractures. However, sigma d appears to be determined by variability within such macroscopic features and/or by micropore-scale heterogeneity. In addition, comparison of these experiments to those at larger spatial scales suggests that increasing sample volume results in an increase in sigma d.     

Henry's law constants for hexachlorobenzene, p,p'-DDE and components of technical chlordane and estimates of gas exchange for Lake Ontario

The Henry's law constants (HLC) for trans- and cis-chlordane (TC, CC), trans-nonachlor (TN), hexachlorobenzene (HCB) and p,p'-DDE were determined by the gas-stripping method over a temperature range of 5-35 degrees C. The HLC variation versus temperature (K) was described by logH=m/T+b. Parameters of this equation were (with standard deviations) TC: m=-1524+/-158, b=6.58+/-0.54; CC: m=-1786+/-209, b=7.42+/-0.71; TN m=-2068+/-284, b=8.44+/-0.97; HCB: m=-3013+/-174, b=11.60+/-0.59 and p,p'-DDE: m=-2043+/-240, b=8.37+/-0.82. The HLCs (Pa m3 mol(-1)) at 25 degrees C (298.15 K) were: TC=29; CC=27; TN=32; p,p'-DDE=33 and HCB=35. These HLCs values were used to calculate fugacity ratios from paired air and water data from Lake Ontario, July 1998. The resulting fugacity ratios predict that volatilization was occurring for all compounds during that month.     

Modelling of silica diffusion experiments with 32Si 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 etaR, the product of the diffusion accessible porosity eta 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 32Si (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(3) g(-1). The apparent diffusion coefficient of dissolved silica in Boom Clay is estimated between 2 x 10(-13) and 7 x 10(-13) m(2) s(-1). The pore diffusion coefficient is in the range from 6 x 10(-11) to 1 x 10(-10) m(2) s(-1).     

Persistent organochlorine pesticides in human milk samples from Australia

Human milk has been used as a surrogate for the assessment of body burden and exposure to persistent lipophilic organochlorine pesticides (OCPs). With the exception of restricted use of mirex, the use of persistent, lipophilic OCPs has been banned in Australia since the 1980s. The decline of human body burden of OCPs following their ban in many industrialised countries is well recorded worldwide from the 1970s until the 1990s though little is known on whether these trends are continuing. In this study, 157 human milk samples collected during 2002 and 2003 as well as 24 samples collected in 1993 were analysed as 20 regional pools for 17 OCPs. OCPs were detected in all pooled human milk samples from 2002/03 typically with highest concentrations of p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE) (mean+/-standard deviation; median concentration 311+/-174; 279 ng g(-1) lipid) followed by beta-hexachlorocyclohexane (beta-HCH) (80+/-173; 21 ng g(-1) lipid). Other OCPs consistently detected included dieldrin (16+/-6; 17 ng g(-1) lipid), hexachlorobenzene (HCB) (18+/-16; 14 ng g(-1) lipid), transnonachlor (11+/-5; 9 ng g(-1) lipid) and p,p'-dichlorodiphenyltrichloroethane (p,p'-DDT) (9+/-6; 7 ng g(-1) lipid). The results from this study indicated that following a substantial decline of OCP concentrations from the early 1980s to the 1990s, little decline could be observed since then. Furthermore the use of multivariate statistics indicated some regional trends with slightly higher levels of the broadly used insecticides DDT and HCH in both historic and recent samples from Melbourne, whereas, sample pools collected from mothers that lived in rural Queensland and New South Wales as well as Adelaide and Sydney showed comparatively higher levels of heptachlor and dieldrin - both of which have been used for termite treatment. These results indicate that even 20 years after the discontinuation of usage, historical use of OCPs rather than exposure via global transport of OCPs is responsible for continuous low exposure in Australia.