Simultaneous estimation of effective and apparent diffusion coefficients in compacted bentonite

Effective diffusion coefficients (D(e)) are usually measured by means of "through-diffusion" experiments in which steady state is reached, and the "time-lag" methods are used to estimate the apparent diffusion coefficient (D(a)). For sorbing radionuclides (as caesium), the time needed to reach steady-state conditions is very large, and the precision in D(a) determinations is not satisfactory. In this paper, a method that allows determining simultaneously effective and apparent diffusion coefficients in compacted bentonite without reaching steady-state conditions is described. Basically, this method consists of an "in-diffusion" experiment in which the concentration profile in the bentonite sample is used to estimate D(a), and the temporal evolution of the solute concentration in the reservoir is used to estimate D(e). This method has several advantages over the typical "through-diffusion" experiments, in particular: (a) experiment duration is significantly shorter, (b) D(a) values are measured with greater precision and (c) it is not necessary to maintain a constant solute concentration in the reservoir. This new method has been used to estimate the effective and apparent diffusion coefficients for caesium in FEBEX bentonite and in order to validate it, the results have been compared with results previously obtained with standard methods. Experimental results have been satisfactorily modelled using a simple model of diffusion in porewater and the measured value of D(e)(Cs) is very similar to D(e)(HTO) in the same bentonite. There is no evidence of "surface diffusion" in FEBEX bentonite for caesium.     

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.     

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.     

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.     

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.     

A field experiment with variable-suction multi-compartment samplers to measure the spatio-temporal distribution of solute leaching in an agricultural soil

Solutes spread out in time and space as they move downwards from the soil surface with infiltrating water. Solute monitoring in the field is often limited to observations of resident concentrations, while flux concentrations govern the movement of solutes in soils. A recently developed multi-compartment sampler is capable of measuring fluxes at a high spatial resolution with minimal disturbance of the local pressure head field. The objective of this paper is to use this sampler to quantify the spatial and temporal variation of solute leaching below the root zone in an agricultural field under natural rainfall in winter and spring. We placed two samplers at 31 and 25 cm depth in an agricultural field, leaving the soil above undisturbed. Each sampler contained 100 separate cells of 31 × 31 mm. Water fluxes were measured every 5 min for each cell. We monitored leaching of a chloride pulse under natural rainfall by frequently extracting the collected leachate while leaving the samplers buried in situ. This experiment was followed by a dye tracer experiment. This setting yielded information that widely surpassed the information that can be provided by separate anionic and dye tracer trials, and solute transport monitoring by coring or suction cups. The detailed information provided by the samplers showed that percolation at the sampling depth started much faster (approximately 3 h after the start of rainfall) in initially wet soil (pressure head above − 65 cm) than in drier soil (more than 14 h at pressure heads below − 80 cm). At any time, 25% of the drainage passed through 5–6% of the sampled area, reflecting the effect of heterogeneity on the flow paths. The amount of solute carried by individual cells varied over four orders of magnitude. The lateral concentration differences were limited though. This suggests a convective–dispersive regime despite the short vertical travel distance. On the other hand, the dilution index indicates a slight tendency towards stochastic–convective transport at this depth. There was no evidence in the observed drainage patterns and dye stained profiles of significant disturbance of the flow field by the samplers.     

A simple dermal absorption model: Derivation and application

This paper deals with the derivation of a QSAR for the estimation of:

• the skin permeation coefficient from aqueous solutions in cm h⁻¹,

• the stratum corneum/water partition coefficient.

These QSARs enable the estimation of:

• the aqueous permeation coefficient in cm h⁻¹,

• the maximum dermal absorption in mg cm⁻² h⁻¹ from a saturated aqueous solution at steady state,

• the lag time in hours (h),

• the diffusivity of a substance in the stratum corneum in cm² h⁻¹

By using the independent variables:

• the log(octanol/water partition coefficient),

• the molecular weight,

• the water solubility.

The estimated maximum dermal absorption and the lag time were compared with some recent measured data of substances, which were not used for developing the QSARs. The estimates were generally in the same order of magnitude as the measured absorption and lag time. These QSARs are recommended for risk assessment of chemicals in the scope of the European REACH legislation.     

Exposure of acid aerosol for schoolchildren in metropolitan Taipei

Metropolitan Taipei, which is located in the subtropical area, is characterized by high population and automobile densities. For convenience, most primary schools are located near major roads. This study explores the exposure of acid aerosols for schoolchildren in areas in Taipei with different traffic densities. Acid aerosols were collected by using a honeycomb denuder filter pack sampling system (HDS). Experimental results indicated that the air pollutants were significantly correlated with traffic densities. The ambient air NO₂, SO₂, HNO₃, NO₃⁻, SO₄²⁻, and aerosol acidity concentrations were 31.3 ppb, 4.7 ppb, 1.3 ppb, 1.9 μg m⁻³, 18.5 μg m⁻³, and 49.5 nmol m⁻³ in high traffic density areas, and 6.1 ppb, 1.8 ppb, 0.9 ppb, 0.7 μg m⁻³, 8.8 μg m⁻³ and 14.7 nmol m⁻³ in low traffic density areas. The exposure levels of acid aerosols for schoolchildren would be higher than the measurements because the sampling height was 5 m above the ground. The SO₂ levels were low (0.13–8.03 ppb) in the metropolitan Taipei. However, the SO₄²⁻ concentrations were relatively high, and might be attributed to natural emissions of sulfur-rich geothermal sources. The seasonal variations of acid aerosol concentrations were also observed. The high levels of acidic particles in spring time may be attributed to the Asian dust storm and low height of the mixture layer. We conclude that automobile contributed not only the primary pollutants but also the secondary acid aerosols through the photochemical reaction. Schoolchildren were exposed to twice the acid aerosol concentrations in high traffic density areas compared to those in low traffic density areas. The incidence of allergic rhinitis of schoolchildren in the high traffic density areas was the highest in spring time. Accompanied by high temperature variation and high levels of air pollution in spring, the health risk of schoolchildren had been observed.     

Volcanic sulphate and arctic dust plumes over the North Atlantic Ocean

High time resolution aerosol mass spectrometry measurements were conducted during a field campaign at Mace Head Research Station, Ireland, in June 2007. Observations on one particular day of the campaign clearly indicated advection of aerosol from volcanoes and desert plains in Iceland which could be traced with NOAA Hysplit air mass back trajectories and satellite images. In conjunction with this event, elevated levels of sulphate and light absorbing particles were encountered at Mace Head. While sulphate concentration was continuously increasing, nitrate levels remained low indicating no significant contribution from anthropogenic pollutants. Sulphate concentration increased about 3.8 μg m⁻³ in comparison with the background conditions. Corresponding sulphur flux from volcanic emissions was estimated to about 0.3 TgS yr⁻¹, suggesting that a large amount of sulphur released from Icelandic volcanoes may be distributed over distances larger than 1000 km. Overall, our results corroborate that transport of volcanogenic sulphate and dust particles can significantly change the chemical composition, size distribution, and optical properties of aerosol over the North Atlantic Ocean and should be considered accordingly by regional climate models.     

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.     

Salting-out phenomenon and 1-octanol/water partition coefficient of metalaxyl pesticide

In this paper, we present the effect of inorganic cations such as Na⁺, K⁺, Ca²⁺, Mg²⁺ on the salting-out phenomenon of metalaxyl from pure water to aqueous salt solutions. Moreover the 1-octanol/water partition coefficient in pure water is presented. To accomplish this, aqueous solubility of metalaxyl was determined in pure water, in different salt solution (NaCl, KCl, CaCl₂ and MgCl₂), and at different concentration level ranging from 0.01 to 1.5 M. The 1-octanol/water partition coefficient was determined using the static shake-flask method. Solubility was determined using dynamic saturation method for pure water in the range of 298.15-325.15 K and at 298.15 K for different salt solutions. The solubility value in pure water for studied interval was found constant (m = 3.118 × 10⁻² mol kg⁻¹).Solubility values were used to calculate the standard molar Gibbs free energy of dissolution (Δ_solG°) and transfer (Δ_trG°) at 298.15 K. The values of Δ_trG° from pure to all studied aqueous salt solutions did not exceed 2 kJ mol⁻¹, the value of Δ_solG° of dissolution is 18.5 ±0.72 kJ mol⁻¹. The 1-octanol/water partition coefficient in pure water log K_o/w is equal to 1.69. The obtained results confirm the classification of the neutral metalaxyl as a slightly hydrophobic molecule.     

Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

An axisymmetric diffusion experiment for the determination of diffusion and sorption coefficients of rock samples

Diffusion anisotropy is a critical property in predicting migration of substances in sedimentary formations with very low permeability. The diffusion anisotropy of sedimentary rocks has been evaluated mainly from laboratory diffusion experiments, in which the directional diffusivities are separately estimated by through-diffusion experiments using different rock samples, or concurrently by in-diffusion experiments in which only the tracer profile in a rock block is measured. To estimate the diffusion anisotropy from a single rock sample, this study proposes an axisymmetric diffusion test, in which tracer diffuses between a cylindrical rock sample and a surrounding solution reservoir. The tracer diffusion between the sample and reservoir can be monitored from the reservoir tracer concentrations, and the tracer profile could also be obtained after dismantling the sample. Semi-analytical solutions are derived for tracer concentrations in both the reservoir and sample, accounting for an anisotropic diffusion tensor of rank two as well as the dilution effects from sampling and replacement of reservoir solution. The transient and steady-state analyses were examined experimentally and numerically for different experimental configurations, but without the need for tracer profiling. These experimental configurations are tested for in- and out-diffusion experiments using Koetoi and Wakkanai mudstones and Shirahama sandstone, and are scrutinized by a numerical approach to identify favorable conditions for parameter estimation. The analysis reveals the difficulty in estimating diffusion anisotropy; test configurations are proposed for enhanced identifiability of diffusion anisotropy. Moreover, it is demonstrated that the axisymmetric diffusion test is efficient in obtaining the sorption parameter from both steady-state and transient data, and in determining the effective diffusion coefficient if isotropic diffusion is assumed. Moreover, measuring reservoir concentrations in an axisymmetric diffusion experiment coupled with tracer profiling may be a promising approach to estimate of diffusion anisotropy of sedimentary rocks.     

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.     

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.     

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.     

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.     

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

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

C2–C8 NMHCs over the Eastern Mediterranean: Seasonal variation and impact on regional oxidation chemistry

More than 2500 measurements of C₂–C₈ non-methane hydrocarbons (NMHCs) have been conducted at Finokalia sampling station on the island of Crete over a thirty-month period (September 2003–February 2006), to investigate the factors controlling NMHC levels and estimate their role in the oxidizing capacity of the Eastern Mediterranean atmosphere. Atmospheric concentrations of NMHCs range from below the detection limit (5 pptv) to a few ppbv and present a hydroxyl radical (OH) driven seasonal pattern with lower values during summer. The diel variability was also influenced by the reaction of the NMHC with the OH radical, exhibiting a nighttime maximum and a midday or early afternoon minimum. Long-lived compounds demonstrate higher concentrations under the influence of the northern sector (European continent), indicating that besides chemistry, transport significantly contributes to NMHCs levels in the area. Based on the observed NMHCs diurnal cycles, mean OH radical levels of 3.5 × 10⁶ molecules cm⁻³ have been derived for May–October period.     

Ionic contribution to the self-potential signals associated with a redox front

In contaminant plumes or in the case of ore bodies, a source current density is produced at depth in response to the presence of a gradient of the redox potential. Two charge carriers can exist in such a medium: electrons and ions. Two contributions to the source current density are associated with these charge carriers (i) the gradient of the chemical potential of the ionic species and (ii) the gradient of the chemical potential of the electrons (i.e., the gradient of the redox potential). We ran a set of experiments in which a geobattery is generated using electrolysis reactions of a pore water solution containing iron. A DC power supply is used to impose a difference of electrical potential of 3 V between a working platinum electrode (anode) and an auxiliary platinum electrode (cathode). Both electrodes inserted into a tank filled with a well-calibrated sand infiltrated by a (0.01 mol L^− 1 KCl + 0.0035 mol L^− 1 FeSO₄) solution. After the direct current is turned off, we follow the pH, the redox potential, and the self-potential at several time intervals. The self-potential anomalies amount to a few tens of millivolts after the current is turned off and decreases over time. After several days, all the redox-active compounds produced initially by the electrolysis reactions are consumed through chemical reactions and the self-potential anomalies fall to zero. The resulting self-potential anomalies are shown to be much weaker than the self-potential anomalies observed in the presence of an electronic conductor in the laboratory or in the field. In the presence of a biotic or an abiotic electronic conductor, the self-potential anomalies can amount to a few hundred millivolts. These observations point out indirectly the potential role of bacteria forming biofilms in the transfer of electrons through sharp redox potential gradient in contaminant plumes that are rich in organic matter.