Groundwater chemistry of the Okélobondo uraninite deposit area (Oklo, Gabon): two-dimensional reactive transport modelling

The stability of uranium-bearing minerals in natural environments is of interest to evaluate the feasibility of radioactive waste repositories. The uraninite bodies, UO2(s), in the Oklo district (Gabon) are the result of a natural fission process, which took place 1970 Ma ago. These deposits can be regarded as natural analogues for spent fuel. One of the uraninite bodies, the Okélobondo deposit, is located at a depth of 300 m. Groundwater samples from boreholes located at shallow depths (100-200 m) show neutral to basic pH, anoxic conditions (Eh = 0.10 to -0.05 V) and are saturated with respect to uraninite. In contrast, deeper samples collected in the vicinity of the ore body are oxidising (Eh = 0.32-0.47 V), slightly basic (pH = 7.0-8.5) and undersaturated with respect to uraninite. These oxidising conditions at depth, if present under repository conditions, may affect the stability of uranium oxide. In order to improve our understanding of the observed site geochemistry, the available information on the lithology and groundwater flow was integrated in a reactive transport model. The chemical composition and the pH-Eh values of the water sampled above and in the western side of the Okélobondo deposit can be explained by the interaction of meteoric recharge with pelites, dolomites and sandstones. The dissolution of Fe(II)-silicates and the oxidation of the Fe(II)-aqueous species maintained the pH-Eh distribution along the Fe(2+)-Fe(OH)3(am) equilibrium, with the result that uraninite does not dissolve. This may explain the lower uranium content in the water samples from pelites and dolomites above the Okélobondo deposit. The high Mn/Fe ratio and the high pH-Eh values of the water sampled at depth, close to the Okélobondo deposit, suggest a control by the Mn(2+)-MnOOH(s) equilibrium. This control is attributed to the dissolution of a large rhodochrosite, MnCO3(s), and manganite, MnOOH(s) deposit in the recharge area on the eastern side.     

Factors influencing natural occurrence of fluoride-rich groundwaters: a case study in the southeastern part of the Korean Peninsula

Factors regulating the occurrence of fluoride in groundwater were investigated using natural isotope tracers and geochemical data in the southeastern part of Korea where about 10% of the total public water-supply wells (n=422) inspected in this study had fluoride levels exceeding the drinking water limit of Korea (1.5 mg/l). The F-rich public wells are mostly distributed along the major faults, especially in the terrain of the F-rich granitic rocks. The stable isotope analysis results provide substantial information for the relative ages of groundwaters. It is revealed that the F-rich groundwaters are deeply circulating paleogroundwaters and occur along the faults due to upward flow along the fault plane. While reacting with granitic rocks for a prolonged period, the F concentrations of groundwater are continuously enriched even after the groundwater reaches an equilibrium state with respect to fluorite (CaF2) due to removal of Ca by precipitation of calcite (CaCO3). These observations reflect that rock chemistry, groundwater age, well depth, and geologic structure are the important factors controlling the occurrence of high F groundwaters. However, high F groundwaters are rarely observed in the fault zones where the associated fractures are widely developed. Isotopic signature provides an evidence for deep penetration of recently recharged groundwater into the wide fault zone, indicating that the hydrologic condition of the fault is also an important factor controlling the occurrence of high F groundwaters.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

A natural analogue of high-pH cement pore waters from the Maqarin area of northern Jordan: Comparison of predicted and observed trace-element chemistry of uranium and selenium

Current design concepts for low-/intermediate-level radioactive waste disposal in many countries involve emplacement underground in a cementitious repository. The highly alkaline groundwaters at Maqarin, Jordan, are a good analogue for the cementitious pore waters that will be present within such a repository. A geochemical modelling study of these groundwaters has been carried out in order to test the applicability of equilibrium models in geochemical programs and their associated thermodynamic databases in such hyperalkaline conditions. This was achieved by comparison of elemental solubilities and speciations predicted by the programs with observations in the natural system. Five organisations took part in the study: AEA Technology, U.K.; Chalmers University of Technology, Sweden; MBT Tecnología Ambiental, Spain; Nagra, Switzerland; and SKB, Sweden. The modelling study was coordinated by the University of Berne.The results of the study showed good agreement between the predictions of the programs employed. Comparison of the observed solids with those predicted by the models has allowed limited validation of the databases. The results for U and Se are presented here.     

Use of carbonised beet pulp carbon for removal of Remazol Turquoise Blue-G 133 from aqueous solution

Carbonised beet pulp (BPC) produced from agricultural solid waste by-product in sugar industry was used as adsorbent for the removal of Remazol Turquoise Blue-G 133 (RTB-G 133) dye in this study. The kinetics and equilibrium of sorption process were investigated with respect to pH, temperature and initial dye concentration. Adsorption studies with real textile wastewater were also performed. The results showed that adsorption was a strongly pH-dependent process, and optimum pH was determined as 1.0. The maximum dye adsorption capacity was obtained as 47.0 mg g^?1at the temperature of 25 °C at this pH value. The Freundlich and Langmuir adsorption models were used for describing the adsorption equilibrium data of the dye, and isotherm constants were evaluated depending on sorption temperature. Equilibrium data of RTB-G 133 sorption fitted very well to the Freundlich isotherm. Mass transfer and kinetic models were applied to the experimental data to examine the mechanisms of adsorption and potential rate-controlling steps. It was found that both external mass transfer and intra-particle diffusion played an important role in the adsorption mechanisms of dye and adsorption kinetics followed the pseudo second-order type kinetic model. The thermodynamic analysis indicated that the sorption process was exothermic and spontaneous in nature.     

Oxidative dissolution of uraninite precipitated on Navajo sandstone

Column and batch experiments were conducted with sandstone and ground water samples to investigate oxidation of uraninite precipitated by microbially mediated reduction of U(VI), a contaminant in ground water beneath a uranium mill tailings site near Tuba City, AZ, USA. Uraninite precipitated together with mackinawite (FeS_0.9) because Fe(III) from the sandstone and sulfate, another contaminant in the water were reduced together with U(VI). After completion of U(VI) reduction, experiments were conducted to find out whether uraninite is protected by mackinawite against reoxidation. Uncontaminated ground water from the same site, containing 7 mg/l of dissolved oxygen, was passed through the columns or mixed with sandstone in batch experiments. The results showed that small masses of uraninite, 0.1 μg/g of sandstone, are protected by mackinawite from reoxidation. Uraninite masses on the order of 0.1 μg/g correspond to U(VI) concentrations of 0.5 mg/l, typically encountered in uranium contaminated ground waters. Mackinawite is an effective buffer and is formed in sufficient quantity to provide long-term protection of uraninite. Uranium concentrations in ground water passed through the columns are too low (4 μg/l) to distinguish between dissolution and oxidative dissolution of uraninite. However, batch experiments showed that uraninite oxidation takes place.     

Evaluating mobilization and transport of arsenic in sediments and groundwaters of Aquia aquifer, Maryland, USA

Groundwater and sediment samples were collected along a flow path in the Aquia aquifer (Paleocene), Maryland in order to examine and study the factors influencing "evolution" of arsenic (As) in these groundwaters. The Aquia crops out near Washington, DC, where it is unconfined, and extends approximately 90 km down dip to the south and east towards and beneath the Chesapeake Bay. The studied flow path was chosen owing to (i) the number of accessible wells, (ii) differences in total dissolved As concentrations in groundwaters from some of the sampled wells, which reach values >/=667 nmol kg(-1) or >/=50 ppb, and (iii) the distinct difference in total dissolved As concentrations in Aquia groundwaters between the northern and southern portions of the study area. In groundwater samples, in situ separation of inorganic As species [As(III) and As(V)] were performed by using anion exchange chromatography. Subsequently, As concentrations were determined by inductively coupled plasma mass spectrometry. In situ measurements of Fe concentrations and speciation, dissolved S(-II) concentrations, pH, alkalinity, and oxidation-reduction potential (Eh) were determined to establish the oxidation-reduction conditions and solution chemistry along the flow path. Concentrations of As in 12 analyzed groundwater samples range from approximately 0.75 to 1 072 nmol kg(-1), and As(III) concentrations ranging from 0.24 to 980 nmol kg(-1) appears to be the dominant form of As in solution. 50% of the studied wells yielded groundwaters with concentrations that exceed the US EPA's Maximum Contaminant Level for As in drinking water of 133 nmol kg(-1) or 10 ppb. In order to examine the solid phase speciation of As within the aquifer sediments, we collected a number of Aquia sediment samples from a drill core that was archived at the Maryland Geological Survey. These sediment samples were evaluated using a previously established sequential extractions procedure. Solid phase As concentrations range between 973 and 2,012 nmol kg(-1). Additionally, petrographic, X-Ray diffraction and diffuse reflectance spectroscopy analyses of the Aquia sediments reveal presence of glauconite, and smectite along with goethite and hematite within the samples. Here, we present the possible mechanisms responsible for the elevated As concentrations in the studied groundwaters of the Aquia aquifer.     

Investigation of surfactant-enhanced dissolution of entrapped nonaqueous phase liquid chemicals in a two-dimensional groundwater flow field.

Because of their low solubility, waste chemicals in the form of nonaqueous phase liquids (NAPLs) that are entrapped in subsurface formations act as long-term sources of groundwater contamination. In the design of remediation schemes that use surfactants, it is necessary to estimate the mass transfer rate coefficients under multi-dimensional flow fields that exit at field sites. In this study, we investigate mass transfer under a two-dimensional flow field to obtain an understanding of the basic mechanisms of surfactant-enhanced dissolution and to quantify the mass transfer rates. Enhanced dissolution experiments in a two-dimensional test cell were conducted to measure rates of mass depletion from entrapped NAPLs to a flowing aqueous phase containing a surfactant. In situ measurement of transient saturation changes using a gamma attenuation system revealed dissolution patterns that are affected by the dimensionality of the groundwater flow field. Numerical modeling of local flow fields that changed with time, due to depletion of NAPL sources, enabled the examination of the basic mechanisms of NAPL dissolution in complex groundwater systems. Through nonlinear regression analysis, mass transfer rates were correlated to porous media properties, NAPL saturation and aqueous phase velocity. Results from the experiments and numerical analyses were used to identify deficiencies in existing methods of analysis that uses assumptions of one-dimensional flow, homogeneity of aquifer properties, local equilibrium and idealized transient mass transfer.     

Hydrogeological and biogeochemical constrains of arsenic mobilization in shallow aquifers from the Hetao basin, Inner Mongolia

Little is known about the importance of drainage/irrigation channels and biogeochemical processes in arsenic distribution of shallow groundwaters from the Hetao basin. This investigation shows that although As concentrations are primarily dependent on reducing conditions, evaporation increases As concentration in the centre of palaeo-lake sedimentation. Near drainage channels, groundwater As concentrations are the lowest in suboxic-weakly reducing conditions. Results demonstrate that both drainage and irrigation channels produce oxygen-rich water that recharges shallow groundwaters and therefore immobilize As. Groundwater As concentration increases with a progressive decrease in redox potential along the flow path in an alluvial fan. A negative correlation between SO₄²⁻ concentrations and δ³⁴S values indicates that bacterial reduction of SO₄²⁻ occurs in reducing aquifers. Due to high concentrations of Fe (>0.5 mg L⁻¹), reductive dissolution of Fe oxides is believed to cause As release from aquifer sediments. Target aquifers for safe drinking water resources are available in alluvial fans and near irrigation channels.     

Biosorption of Cr(VI) and Zn(II) ions from aqueous solution onto the solid biodiesel waste residue: mechanistic, kinetic and thermodynamic studies

In this present study, the biosorption of Cr(VI) and Zn(II) ions from synthetic aqueous solution on defatted J atropha oil cake (DJOC) was investigated. The effect of various process parameters such as the initial pH, adsorbent dosage, initial metal ion concentration and contact time has been studied in batch-stirred experiments. Maximum removal of Cr(VI) and Zn(II) ions in aqueous solution was observed at pH 2.0 and pH. 5.0, respectively. The removal efficiency of Cr(VI) and Zn(II) ions from the aqueous solution was found to be 72.56 and 79.81 %, respectively, for initial metal ion concentration of 500 mg/L at 6 g/L dosage concentration. The biosorbent was characterized by Fourier transform infrared, scanning electron microscopy and zero point charge. Equilibrium data were fitted to the Langmuir, Freundlich, Temkin and Dubinin–Radushkevich isotherm models and the best fit is found to be with the Freundlich isotherm for both Cr(VI) and Zn(II) metal ions. The kinetic data obtained at different metal ion concentration have been analysed using the pseudo-first-order, pseudo-second-order and intraparticle diffusion models and were found to follow the pseudo-second-order kinetic model. The values of mass transfer diffusion coefficients (D _e) were determined by Boyd model and compared with literature values. Various thermodynamic parameters, such as ΔG°, ΔH° and ΔS°, were analysed using the equilibrium constant values (K _e) obtained from experimental data at different temperatures. The results showed that biosorption of Cr(VI) and Zn(II) ions onto the DJOC system is more spontaneous and exothermic in nature. The results indicate that DJOC was shown to be a promising adsorbent for the removal of Cr(VI) and Zn(II) ions from aqueous solution.     

2,4-二氯苯酚在土壤中的吸附及其批实验与柱实验的分配系数比较

为了研究2,4-二氯苯酚在土壤中的吸附及比较其批实验与柱实验的分配系数Kd,开展了2,4-二氯苯酚的批实验（不同液固比条件下）和柱实验。通过分析结果可知,在批实验中,不同液固比条件下2,4-二氯苯酚达到平衡的时间类似,都在60～70 h,吸附动力学曲线符合伪二级动力学方程,吸附规律是：液固比越大,平衡吸附量增大,反应速率常数K2减小,初始吸附速率常数减小;Kd随液固比增大而降低,范围在2.91～2.12 L/kg。柱实验结果表明,2,4-二氯苯酚的贯穿曲线可以很好地用化学非平衡模型来拟合,通过模型拟合得到的Kd值要低于批实验的结果。该研究对表征2,4-二氯苯酚在环境中的行为、预测其对土壤和地下水的污染及其治理提供了依据。     

Equilibration time scales of organic aerosol inside thermodenuders: Evaporation kinetics versus thermodynamics

The interpretation of thermodenuder (TD) data often relies on the assumption that thermodynamic equilibrium is reached inside the instrument. We modeled the evaporation of three organic aerosol types (adipic acid, α-pinene SOA and aged OA) inside a thermodenuder with a mass transfer model, and calculated equilibration time scales for these systems at realistic conditions. The equilibrium times varied from less than a second to several hours, decreasing with increasing aerosol concentrations, decreasing particle sizes, decreasing volatilities and increasing mass accommodation coefficients. The results indicate that generally TDs measure particle evaporation rates rather than equilibria, and time-dependent modeling of the evaporation is usually needed to interpret the data. Measurements at varying residence times and temperatures, on the other hand, are desirable to investigate the equilibration of the studied aerosol and decouple the kinetic effects from the effects caused by the thermodynamic properties of the aerosol. Organic aerosol is likely to be further from equilibrium under typical field conditions compared with laboratory data. When determining the aerosol properties from TD data, assuming incorrectly equilibrium results in under-prediction of the vaporization enthalpy of the evaporating species. Similar under-estimation is predicted if multicomponent aerosols are approximated with single-component properties.     

Soil column experiments used as a means to assess transport,sorption, and biodegradation of pesticides in groundwater

Soil column experiments are used to investigate the fate of three pesticides of high, intermediate, and low solubility in groundwater: N- phosphonomethyl glycine (glyphosate); O,O-diethyl-S-[(ethylthio)methyl]phosphorodithioate (phorate); (2,4-dichlorophenoxy)acetic acid (2,4-D). Feed solutions are prepared by adding each pesticide (100 mg/L glyphosate, 50 μ g/L phorate, 50 mg/L 2,4-D) along with conservative tracer, KBr, in synthetic groundwater. The concentration of the pesticides in effluents is detected by ion chromatography (glyphosate, 2,4-D) and GC-FID (phorate). The Br^? breakthrough curves are employed to estimate the dispersion coefficient and mean pore velocity in each column. Solute transport and reactive models accounting for equilibrium/non-equilibrium sorption and biodegradation are coupled with inverse modeling numerical codes to estimate the kinetic parameters for all pesticides.     

Processes controlling metal transport and retention as metal-contaminated groundwaters efflux through estuarine sediments

Factors affecting the transport and retention of Cd, Cr, Cu, Ni, Pb and Zn in acidic groundwaters as they pass through estuarine sediments were investigated using column experiments. Acidic groundwaters caused the rapid dissolution of iron sulfide (AVS) and other iron and manganese phases from sediments that are important for metal binding and buffering. Metal breakthrough to overlying water occurred in the order of Ni>Zn>Cd>Cu>Cr/Pb. Metal transport increased as the sediment permeability increased, reflecting the low resistance to flow caused by larger sand-sized particles and the decreased abundance of metal adsorption sites on these materials. Metal mobility increased as the groundwater pH decreased, as flow rate or metal concentrations increased, and as the exposure duration increased. Groundwater Cr and Pb were promptly attenuated by the sediments, the mobility of Cu was low and decreased rapidly as sediment pH increased above 4.5, while Cd, Ni and Zn were the most easily transported to the surface sediments and released to the overlying waters. For groundwaters of pH 3, metal migration velocities through sandy sediments were generally 0.5-2% (Cr, Pb), 1-6% (Cu) and 4-13% (Cd, Ni, Zn) of the total groundwater velocity (9-700 m/yr). The oxidative precipitation of Fe(II) and Mn(II) in the groundwaters did not affect metal mobility through the sediments. The results indicated that the efflux of acidic and metal-contaminated groundwater through estuarine sediments would affect organisms resident in sandy sediments more greatly than organisms resident in fine-grained, silty, sediments.     

Long-term geochemical evolution of the near field repository: insights from reactive transport modelling and experimental evidences

The KBS-3 underground nuclear waste repository concept designed by the Swedish Nuclear Fuel and Waste Management Co. (SKB) includes a bentonite buffer barrier surrounding the copper canisters and the iron insert where spent nuclear fuel will be placed. Bentonite is also part of the backfill material used to seal the access and deposition tunnels of the repository. The bentonite barrier has three main safety functions: to ensure the physical stability of the canister, to retard the intrusion of groundwater to the canisters, and in case of canister failure, to retard the migration of radionuclides to the geosphere. Laboratory experiments (< 10 years long) have provided evidence of the control exerted by accessory minerals and clay surfaces on the pore water chemistry. The evolution of the pore water chemistry will be a primordial factor on the long-term stability of the bentonite barrier, which is a key issue in the safety assessments of the KBS-3 concept.In this work we aim to study the long-term geochemical evolution of bentonite and its pore water in the evolving geochemical environment due to climate change. In order to do this, reactive transport simulations are used to predict the interaction between groundwater and bentonite which is simulated following two different pathways: (1) groundwater flow through the backfill in the deposition tunnels, eventually reaching the top of the deposition hole, and (2) direct connection between groundwater and bentonite rings through fractures in the granite crosscutting the deposition hole. The influence of changes in climate has been tested using three different waters interacting with the bentonite: present-day groundwater, water derived from ice melting, and deep-seated brine. Two commercial bentonites have been considered as buffer material, MX-80 and Deponit CA-N, and one natural clay (Friedland type) for the backfill. They show differences in the composition of the exchangeable cations and in the accessory mineral content. Results from the simulations indicate that pore water chemistry is controlled by the equilibrium with the accessory minerals, especially carbonates. pH is buffered by precipitation/dissolution of calcite and dolomite, when present. The equilibrium of these minerals is deeply influenced by gypsum dissolution and cation exchange reactions in the smectite interlayer. If carbonate minerals are initially absent in bentonite, pH is then controlled by surface acidity reactions in the hydroxyl groups at the edge sites of the clay fraction, although its buffering capacity is not as strong as the equilibrium with carbonate minerals. The redox capacity of the bentonite pore water system is mainly controlled by Fe(II)-bearing minerals (pyrite and siderite). Changes in the groundwater composition lead to variations in the cation exchange occupancy, and dissolution–precipitation of carbonate minerals and gypsum. The most significant changes in the evolution of the system are predicted when ice-melting water, which is highly diluted and alkaline, enters into the system. In this case, the dissolution of carbonate minerals is enhanced, increasing pH in the bentonite pore water. Moreover, a rapid change in the population of exchange sites in the smectite is expected due to the replacement of Na for Ca.     

Soil column experiments used as a means to assess transport, sorption, and biodegradation of pesticides in groundwater

Soil column experiments are used to investigate the fate of three pesticides of high, intermediate, and low solubility in groundwater: N- phosphonomethyl glycine (glyphosate); O,O-diethyl-S-[(ethylthio)methyl]phosphorodithioate (phorate); (2,4-dichlorophenoxy)acetic acid (2,4-D). Feed solutions are prepared by adding each pesticide (100 mg/L glyphosate, 50 micro g/L phorate, 50 mg/L 2,4-D) along with conservative tracer, KBr, in synthetic groundwater. The concentration of the pesticides in effluents is detected by ion chromatography (glyphosate, 2,4-D) and GC-FID (phorate). The Br(-) breakthrough curves are employed to estimate the dispersion coefficient and mean pore velocity in each column. Solute transport and reactive models accounting for equilibrium/non-equilibrium sorption and biodegradation are coupled with inverse modeling numerical codes to estimate the kinetic parameters for all pesticides.     

Modelling the solute transport under nonequilibrium conditions on the basis of mass transfer equations

A solute transport model that describes nonequilibrium adsorption in soil/groundwater systems by mass transfer equations for film and intraparticle diffusion is presented. The model is useful in cases where breakthrough curve spreading cannot be explained by dispersion only. To evaluate its validity, the model was applied to several data sets from column experiments. The validity was also proved by a comparison with an analytical solution for the limiting case of predominating dispersion. Furthermore, a sensitivity analysis was performed to illustrate the influence of different process and sorption parameters (pore water velocity, intraparticle mass transfer coefficient, isotherm nonlinearity) on the shape of the calculated breakthrough curves. The application of the proposed model is discussed in comparison to the widely used dispersed flow/local equilibrium model, and a relationship between both models, which is based on a lumped parameter approach, is shown.     

Intercomparison of reactive transport models applied to UO2 oxidative dissolution and uranium migration

Oxidative dissolution of uranium dioxide (UO₂) and the subsequent migration of uranium in a subsurface environment and an underground waste disposal have been simulated with reactive transport models. In these systems, hydrogeological and chemical processes are closely entangled and their interdependency has been analyzed in detail, notably with respect to redox reactions, kinetics of mineralogical evolution and hydrodynamic migration of species of interest.Different codes, where among CASTEM, CHEMTRAP and HYTEC, have been used as an intercomparison and verification exercise. Although the agreement between codes is satisfactory, it is shown that the discretization method of the transport equation (i.e. finite elements (FE) versus mixed-hybrid FE and finite differences) and the sequential coupling scheme may lead to systematic discrepancies.