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.     

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.     

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).     

An analytical model for the interpretation of pulse injection experiments performed for testing the spatial variability of clay formations

This paper presents an analytical model to describe pulse injection experiments. This model solves the advection-diffusion equation while taking into account back diffusion from the clay core to the inlet and from the outlet to the clay core. In most analytical models, back diffusion is neglected. For sufficiently high Péclet numbers, this is a good approximation. However, in experiments where the Péclet number is low, back diffusion is important and must be taken into account. An additional advantage of the present model is that both concentration and flux are conserved at the inlet and at the outlet of the clay core. This model is used to fit pulse injection experiments with iodide and tritiated water (HTO) in clay cores. The (new) model is required for fitting the experimental results since in clay layers advection is very slow leading to a low Péclet number. The experiments are performed on clay cores taken from different depths from the Boom Clay and the Ypres Clay layer under the site of the nuclear power plant of Doel (Belgium). The quality of all fits is excellent and the obtained parameter values are coherent. For HTO, the fitted value for the diffusion accessible porosity is consistent with measurements of the water content in Ypres Clay cores. In both types of clays, the apparent diffusion coefficient at zero flow is between 10(-10) and 2 x 10(-10) m(2)/s for iodide and between 2 x 10(-10) and 3 x 10(-10) m(2)/s for HTO. The dispersion length is in the order of 10(-3) m. The average value for the diffusion accessible porosity is between 0.35 and 0.4 for HTO and between 0.2 and 0.25 for iodide.     

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.     

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.     

A technique for estimating one-dimensional diffusion coefficients in low-permeability sedimentary rock using X-ray radiography: comparison with through-diffusion measurements

The measurement of diffusive properties of low-permeability rocks is of interest to the nuclear power industry, which is considering the option of deep geologic repositories for management of radioactive waste. We present a simple, non-destructive, constant source in-diffusion method for estimating one-dimensional pore diffusion coefficients (D(p)) in geologic materials based on X-ray radiography. Changes in X-ray absorption coefficient (Deltamicro) are used to quantify changes in relative concentration (C/C(0)) of an X-ray attenuating iodide tracer as the tracer solution diffuses through the rock pores. Estimated values of D(p) are then obtained by fitting an analytical solution to the measured concentration profiles over time. Measurements on samples before and after saturation with iodide can also be used to determine iodide-accessible porosity (phi(I)). To evaluate the radiography method, results were compared with traditional steady-state through-diffusion measurements on two rock types: shale and limestone. Values of D(p) of (4.8+/-2.5)x10(-11) m(2).s(-1) (mean+/-standard deviation) were measured for samples of Queenston Formation shale and (2.6+/-1.0)x10(-11) m(2).s(-1) for samples of Cobourg Formation limestone using the radiography method. The range of results for each rock type agree well with D(p) values of (4.6+/-2.0)x10(-11) m(2).s(-1) for shale and (3.5+/-1.8)x10(-11) m(2).s(-1) for limestone, calculated from through-diffusion experiments on adjacent rock samples. Low porosity (0.01 to 0.03) and heterogeneous distribution of porosity in the Cobourg Formation may be responsible for the slightly poorer agreement between radiography and through-diffusion results for limestones. Mean values of phi(I) for shales (0.060) and limestones (0.028) were close to mean porosity measurements made on bulk samples by the independent water loss technique (0.062 and 0.020 for shales and limestones, respectively). Radiography measurements offer the advantage of time-saving for diffusion experiments because the experiment does not require steady-state conditions and also allows for visualization of the small-scale heterogeneities in diffusive properties within rocks at the mm to cm scale.     

Modeling the Reduction of Vapor Phase Emissions from Surface Soils Due to Soil Matrix Effects: Porosity/Tortuosity Concepts

ABSTRACT

A major route for transport of volatile organic compounds within porous media is vapor phase diffusion. The diffusion rate through a porous medium is less than that through free-air due to the decreased cross-sectional area available for gas movement and the increased path length due to pore tortuosity. Numerous empirical expressions have been published that relate the diffusion coefficient in porous media to the diffusion coefficient in free-air (unobstructed gas phase). Published measurements of relative diffusivity and air-filled porosity were combined into a database. Empirical expressions available in the literature, including the popular Millington-Quirk equation, were evaluated along with a fourth-degree polynomial expression developed by the authors to determine the best type of equation to predict relative diffusivity as a function of air-filled porosity over the domain of values for porosity ranging from 0.071 to 1 for different types of materials. Mean square deviations were used as the statistical test to compare equations. The polynomial expression developed in this project produced a significantly different effective diffusion coefficient (1.3 x 10^-6 m²/sec) compared to values of 9.2 x 10^-6 m²/sec and 3.1 x 10^-6 m²/ sec predicted by forms of the Millington-Quirk equation for a specific case.     

The diffusion and sorption of volatile organic compounds through kaolinitic clayey soils

Laboratory experiments to estimate the effective molecular diffusion coefficient (D(e)) and sorption coefficient (K(d)) for volatile organic compounds through natural clayey soils were conducted using diffusion testing apparatus. The compounds tested were methyl ethyl ketone (MEK), toluene and trichloroethylene (TCE). The D(e) and K(d) values were determined by a curve fitting procedure. The compound losses, and the effects of porous disks used in the apparatus were significant. The transport of MEK was faster than that of TCE and toluene because of the lower sorption to the soils. The D(e) values of all the compounds were of the order of 10(-10) m(2)/s and smaller than the diffusion coefficient in pure aqueous solution at infinite dilution (D(0)), due to the tortuosity of the samples. The effects of the sample thickness on the parameter determination were not significant. Comparison to the K(d) values estimated from batch sorption tests and from organic carbon content (f(oc))-based predictions showed that the diffusion test results were intermediate between those from the other two methods. The diffusion tests use compacted soil samples and should be more relevant to in situ conditions, but the reliability of the tests is affected by large compound losses that cause uncertainties in their interpretation. It is recommended that more than one method be used to assess K(d) values.     

In situ anion diffusion experiments using radiotracers

Diffusion experiments in compacted bentonite have been carried out in situ using the borehole laboratory CHEMLAB. The "ordinary" anion iodide and the redox-sensitive pertechnetate ion have been investigated. In spite of strongly reducing groundwater conditions, technetium was found to diffuse mostly unreduced as TcO4-, although in some spots in the compacted clay, the activity was significantly higher, which may be explained by reduction of some TcO4- by iron-containing minerals in the bentonite. The measured concentration profiles in the clay cannot be accommodated by assuming one single diffusion process. The experimental data are modeled assuming two diffusion paths, intralamellar diffusion and diffusion in external water. The apparent diffusivity for the intralamellar diffusion was found to be 8.6 x 10(-11) m2 s(-1) for iodide with a capacity factor of 0.1, while the apparent diffusivity for the diffusion in external water was found to be 5 x 10(-14) m2 s(-1) with alpha=2.26. The corresponding values for Tc were found to be Da= 6 x 10(-11) m2 s(-1), alpha=0.1 and Da= 1 x 10(-13) m2 s(-1), alpha=0.46, respectively. The diffusion constants and capacity factors obtained in this study are in accordance with data from laboratory experiments.     

Theoretical consideration of methane emission from sediments

A stationary theory of gas emission from sedimentary (active) layers of wetlands is developed. The theory takes into account methane generation in a sedimentary layer (W1) and its depth dependence, (W1(z)), the solubility, determined by Henry's constant (K1), and the mobility of methane molecules set by the methane diffusion coefficient (D1). The exponential dependence of methane generation rate decay with depth is considered in more detail. The penetration of atmospheric nitrogen into the active layer is also taken into consideration. It is shown that the value of diffusion methane flux from sedimentary layers is proportional to square root(K1D1P0W10) where P0 characterizes the atmospheric pressure and W10 characterizes the maximum generation rate. Coefficients relating the diffusion methane flux to the square root(K1D1P0W10) value are calculated for the different depth dependencies of methane generation rate. The values of these coefficients are not much different from unity for most real cases.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

Bio-oxidation of airborne volatile organic compounds in an activated sludge aeration tank

An activated sludge aeration tank (40 x 40 x 300 cm, width x length x height) with a set of 2-mm orifice air spargers was used to treat gas-borne volatile organic compounds (VOCs; toluene, p-xylene, and dichloromethane) in air streams. The effects of liquid depth (Z), aeration intensity (G/A), the overall mass-transfer rate of oxygen in clean water (KLaO2), the Henry's law constant of the tested VOC (H), and the influent gaseous VOC concentration (C0) on the efficiency of removal of VOCs were examined and compared with a literature-cited model. Results show that the measured VOC removal efficiencies and those predicted by the model were comparable at a G/A of 3.75-11.25 m3/m2 hr and C0 of approximately 1000-6000 mg/m3. Experimental data also indicated that the designed gas treatment reactor with KLaO2 = 5-15 hr(-l) could achieve > 85% removal of VOCs with H = 0.24-0.25 at an aerated liquid depth of 1 m and > 95% removal of dichloromethane with H = 0.13 at a 1-m liquid depth.     

Sorptive behavior of sorgoleone in ultisol in two solvent systems and determination of its lipophilicity

Sorgoleone (SGL) exuded by sorghum roots inhibits the development of some weeds. Due to its high hydrophobicity, it is expected that SGL presents low soil mobility and limited allelopathic activity in the field. This work aims to evaluate the sorptivity of sorgoleone in octanol-water and in soil under two solvent systems. The two solvent systems were methanol:water (60:40) (MeOH:H2O) and pure methanol (MeOH). These two solvent systems promote different conditions for SGL solubility. Treatments were arranged in a 2 x 6 factorial (solvent systems x equilibrium concentrations in the solution (EC)). For each solvent, the sorption was achieved by shaking 500 mg of soil with 10 ml of 0, 5, 10, 15, 25, 40, and 60 mg L-1 of SGL solution, during 24 h. After centrifugation, the supernatant was filtered and the SGL concentration was determined by high performance liquid chromatography (HPLC). Data of sorbed amount of SGL were submitted to variance analysis, using a hierarchic factorial model. The data of sorbed amount (x/m) and equilibrium concentration (C) were fitted to the linear (x/m = a + KdC) and to the Freundlich (x/m = KfC1/n) models. The isotherm obtained for the MeOH:H2O system presented linear shape, whereas for the MeOH system a two subsequent linear isotherm was fitted. Sorgoleone is a highly hydrophobic compound, presenting a log Kow of 6.1. The sorption of sorgoleone to the soil was very high. The organic environment stimulated the sorgoleone sorption to the soil.     

Fate and transport of monoterpenes through soils. Part I. Prediction of temperature dependent soil fate model input-parameters

Monoterpenes are C10H(n)O(n') compounds of natural origin and are potentially environmentally safe substitutes for traditional pesticides. Still, an assessment of their environmental behaviour is required. As a first step in a theoretical study focussing on monoterpenes applied as pesticides to terrestrial environments, soil fate model input-parameters were determined for 20 monoterpenes with widely different structural characteristics. Input-parameters are the water solubility (S(W)), vapour pressure (P), n-octanol-water partition coefficient (K(OW)), atmospheric air and bulk water diffusion coefficients (D(A)air and D(W)water), first order biodegradation rate constants (k), and their temperature dependence. Values for these parameters were estimated or taken from previous experimental work. The quality of the estimations was discussed by focussing on their statistics and by comparison with available experimental data. From these properties, the air-water partition coefficient (K(AW), Henry's Law constant), the interface-water partition coefficient (K(IW)) and the organic matter-water partition coefficient (K(OM)) could be estimated with varying levels of accuracy. In general, little experimental data turned out to be available on biodegradation rate constants and on the temperature dependence of physico-chemical parameters.     

Sandstones of unexpectedly high diffusibility

Measurements have been made of diffusion coefficients (D(i)=-mass flux/concentration gradient) using a double reservoir, steady-state method with two tracers, CaBr(2) and amino-G-acid, on intact samples of Triassic red-bed sandstone from northwest England. Diffusibility (D'=D(i)/diffusion coefficient in water) averages 0.124, ranging between 0.075 and 0.215 (porosity 0.1 to 0.24), very similar for the two tracers. Implied tortuosities (actual path length/straight line length) average 1.21 (range 1.06 to 1.47), with constrictivities close to 1. In comparison with limited red-bed sandstone data from elsewhere, these D' values are up to 4 times greater, and tortuosity correspondingly lower. Re-interpretation of formation factor data from previous studies on shallow sandstone samples also from northwest England confirms that diffusibility is significantly higher in these sandstones than others from similar palaeoenvironment/stratigraphic units. The lower tortuosities appear to result from the relatively high permeability, open fabric of the rock, properties likely to be present in shallow sandstone systems used for water supply. It is concluded that diffusion rates may, in some shallow freshwater-containing continental sandstone systems, be significantly greater than is implied by estimates of sandstone diffusibility current in the literature.     

Qualitative and quantitative effects of ozone and/or sulfur dioxide on field-grown potato plants

Solanum tuberosum L. cv Norchip plants were grown in open-top chambers in the summer of 1986. Plants were treated with charcoal-filtered air, nonfiltered air, or nonfiltered air supplemented with 33, 66, or 99% of the ambient ozone (O3) concentrations from 1000 to 2000 h eastern daylight time daily. In addition, plants received charcoal-filtered air plus 0, 0.15 (393 microg m(-3)), 0.34 (891 microg m(-3)), or 0.61 (1598 microg m(-3)) ppm sulfur dioxide (SO2) from 0900 to 1200 h once every 14 d for a total of four treatments. Ozone induced a linear reduction in number and weight of Grade One (> 6.35-cm diameter) potato tubers and in total weight of tubers. Ozone also induced linear reductions in the percentage of dry matter of tubers and linear decreases in glucose and fructose content of Grade One tubers. Sulfur dioxide induced a stimulation and then decline of the number, percentage of dry matter, and sucrose content of Grade One tubers. The SO2 response best fit a quadratic curve. No O3 x SO2 interactions were detected for any of the yield or quality functions measured.     

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.     

Experimental and modeling studies on sorption and diffusion of radium in bentonite

The sorption and desorption behavior of radium on bentonite and purified smectite was investigated as a function of pH, ionic strength and liquid to solid ratio by batch experiments. The distribution coefficients (Kd) were in the range of 10(2) to > 10(4) ml g-1 and depended on ionic strength and pH. Most of sorbed Ra was desorbed by 1 M KCl. The results for purified smectite indicated that Ra sorption is dominated by ion exchange at layer sites of smectite, and surface complexation at edge sites may increase Ra sorption at higher pH region. Reaction parameters between Ra and smectite were determined based on an interaction model between smectite and groundwater. The reaction parameters were then used to explain the results of bentonite by considering dissolution and precipitation of minerals and soluble impurities. The dependencies of experimental Kd values on pH, ionic strength and liquid to solid ratio were qualitatively explained by the model. The modeling result for bentonite indicated that sorption of Ra on bentonite is dominated by ion exchange with smectite. The observed pH dependency was caused by changes of Ca concentration arising from dissolution and precipitation of calcite. Diffusion behavior of Ra in bentonite was also investigated as a function of dry density and ionic strength. The apparent diffusion coefficients (Da) obtained in compacted bentonite were in the range of 1.1 x 10(-11) to 2.2 x 10(-12) m2 s-1 and decreased with increasing in dry density and ionic strength. The Kd values obtained by measured effective diffusion coefficient (De) and modeled De were consistent with those by the sorption model in a deviation within one order of magnitude.     

Numerical simulation of an array of fences in Saemangeum reclaimed land

This paper discusses about the quantitative effect of windbreak fences on wind velocity in the reclaimed land at Saemangeum in South Korea. Windbreak fences were constructed in the reclaimed land to reduce the wind velocity to prevent the generation and diffusion of dust. However, up to this time, no in-depth studies were conducted to quantitatively measure the effect of the windbreak fences on wind velocity thus an optimum windbreak structure is not yet determined. Using CFD simulations, the effects of fence porosity, fence height, and the distance between the adjacent fences were investigated. A wind tunnel experiment was initially conducted and data gathered were used to develop the CFD models. From the experiments and CFD simulations, the overall percentage difference of the measured velocities was 7.20% which is generally acceptable to establishing the reliability of the CFD models. The reduction effect on wind velocity was measured in between the adjacent fences up to a height of 0.6 m from the ground surface. In terms of porosity ( = 0, 0.2, 0.4, 0.6), 0.2 was found to be the optimum value. Conversely, the effect of fence height (0.6, 0.8 and 1.0 m) showed no significant difference; therefore, 0.6 m height is recommended. In addition, the reduction effect of distance between the adjacent fences (2, 4 and 6 m) on wind velocity having a 0.2 porosity has decreased to about 75% regardless of the distance. In the case of the reclaimed land in Saemangeum, a decrease of 75% can prevent the generation and diffusion of dusts. However, the source of dusts is very large. Therefore, constructing an array of windbreak with 6 m distance between them is deemed necessary.