Inverse modeling of multicomponent reactive transport through single and dual porosity media

Compacted bentonite is foreseen as buffer material for high-level radioactive waste in deep geological repositories because it provides hydraulic isolation, chemical stability, and radionuclide sorption. A wide range of laboratory tests were performed within the framework of FEBEX (Full-scale Engineered Barrier EXperiment) project to characterize buffer properties and develop numerical models for FEBEX bentonite. Here we present inverse single and dual-continuum multicomponent reactive transport models of a long-term permeation test performed on a 2.5 cm long sample of FEBEX bentonite. Initial saline bentonite porewater was flushed with 5.5 pore volumes of fresh granitic water. Water flux and chemical composition of effluent waters were monitored during almost 4 years. The model accounts for solute advection and diffusion and geochemical reactions such as aqueous complexation, acid-base, cation exchange, protonation/deprotonation by surface complexation and dissolution/precipitation of calcite, chalcedony and gypsum. All of these processes are assumed at local equilibrium. Similar to previous studies of bentonite porewater chemistry on batch systems which attest the relevance of protonation/deprotonation on buffering pH, our results confirm that protonation/deprotonation is a key process in maintaining a stable pH under dynamic transport conditions. Breakthrough curves of reactive species are more sensitive to initial porewater concentration than to effective diffusion coefficient. Optimum estimates of initial porewater chemistry of saturated compacted FEBEX bentonite are obtained by solving the inverse problem of multicomponent reactive transport. While the single-continuum model reproduces the trends of measured data for most chemical species, it fails to match properly the long tails of most breakthrough curves. Such limitation is overcome by resorting to a dual-continuum reactive transport model.     

Non-Fickian dispersion in a single fracture

Solute transport in fractured rocks is of major interest in many applications, from the petroleum industry to ground water management. This work focuses on the dispersion process in a transparent replica of a real single fracture. The fracture exhibits strong changes in heterogeneity, with the first half very heterogeneous and the second half fairly homogeneous. Three models have been used to interpret the tracer experiments: the classical advection-dispersion equation (ADE), the continuous time random walk (CTRW), and the stratified model. The main goals were to test these models and to study possible correlations between fitting parameters and heterogeneities. As expected, the solution derived from the ADE equation appears to be unable to model long-time tailing behavior. On the other hand, the results confirm the CTRW robustness and the coefficient beta seems well correlated to heterogeneities. Finally, the stratified model is also able to describe non-Fickian dispersion. The parameters defined by this model are correlated to the heterogeneities of the fracture.     

Identification of key parameters controlling dissolved oxygen migration and attenuation in fractured crystalline rocks

In the crystalline rocks of the Canadian Shield, geochemical conditions are currently reducing at depths of 500-1000 m. However, during future glacial periods, altered hydrologic conditions could potentially result in enhanced recharge of glacial melt water containing a relatively high concentration of dissolved oxygen (O2). It is therefore of interest to investigate the physical and geochemical processes, including naturally-occurring redox reactions, that may control O2 ingress. In this study, the reactive transport code MIN3P is used in combination with 2k factorial analyses to identify the most important parameters controlling oxygen migration and attenuation in fractured crystalline rocks. Scenarios considered are based on simplified conceptual models that include a single vertical fracture, or a fracture zone, contained within a rock matrix that extends from the ground surface to a depth of 500 m. Consistent with field observations, Fe(II)-bearing minerals are present in the fractures (i.e. chlorite) and the rock matrix (biotite and small quantities of pyrite). For the parameter ranges investigated, results indicate that for the single fracture case, the most influential factors controlling dissolved O2 ingress are flow velocity in the fracture, fracture aperture, and the biotite reaction rate in the rock matrix. The most important parameters for the fracture zone simulations are flow velocity in the individual fractures, pO2 in the recharge water, biotite reaction rate, and to a lesser degree the abundance and reactivity of chlorite in the fracture zone, and the fracture zone width. These parameters should therefore receive increased consideration during site characterization, and in the formulation of site-specific models intended to predict O2 behavior in crystalline rocks.     

Mixing-controlled biodegradation in a toluene plume--results from two-dimensional laboratory experiments

Various abiotic and biotic processes such as sorption, dilution, and degradation are known to affect the fate of organic contaminants, such as petroleum hydrocarbons in saturated porous media. Reactive transport modeling of such plumes indicates that the biodegradation of organic pollutants is, in many cases, controlled by mixing and therefore occurs locally at the plume's fringes, where electron donors and electron-acceptors mix. Herein, we aim to test whether this hypothesis can be verified by experimental results obtained from aerobic and anaerobic degradation experiments in two-dimensional sediment microcosms. Toluene was selected as a model compound for oxidizable contaminants. The two-dimensional microcosm was filled with quartz sand and operated under controlled flow conditions simulating a contaminant plume in otherwise uncontaminated groundwater. Aerobic degradation of toluene by Pseudomonas putida mt-2 reduced a continuous 8.7 mg L(-1) toluene concentration by 35% over a transport distance of 78 cm in 15.5 h. In comparison, under similar conditions Aromatoleum aromaticum strain EbN1 degraded 98% of the toluene infiltrated using nitrate (68.5+/-6.2 mg L(-1)) as electron acceptor. A major part of the biodegradation activity was located at the plume fringes and the slope of the electron-acceptor gradient was steeper during periods of active biodegradation. The distribution of toluene and the significant overlap of nitrate at the plume's fringe indicate that biokinetic and/or microscale transport processes may constitute additional limiting factors. Experimental data is corroborated with results from a reactive transport model using double Monod kinetics. The outcome of the study shows that in order to simulate degradation in contaminant plumes, detailed data sets are required to test the applicability of models. These will have to deal with the incorporation of existing parameters coding for substrate conversion kinetics and microbial growth.     

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.     

危险品道路运输泄漏引发水污染事故的定量风险评估方法研究

介绍了由道路交通事故引发的地表水污染事故风险定量分析方法,重点分析了突发水污染事故源发生概率和危害后果表征方法.以评价水域内道路危险品运输泄漏事故发生概率为水污染事故源发生概率.污染事故后果危害水平的表征:1)敏感区(点)最高浓度; 2)敏感区(点)超标总历时;3)水质超标的范围(如超标河流长度).评估方法为陆源造成的突发水污染事故的预防管理与应急计划的制订实施提供了科学依据.     

Comparison of Microscopic and Macroscopic Modeling Approaches for Subsurface Contaminant Transport

The practice of contaminant transport and remediation has shown significant progress in recent years. However, despite the significant progress made, remediation efforts are often delayed by extremely long breakthrough curve tails that render efforts to bring the level of contaminants below the regulatory standards inefficient. One hypothesis is that these long tails are due to the reservoir-like slow diffusive processes in soil micropore zones. This study compares the effects of micropores at macroscopic and microscopic levels and establishes a link between these approaches for validation and calibration purposes. The link between macroscopic and microscopic levels is established through comparisons and testing of the two models while incorporating appropriate scale and boundary effects. Despite the differences in conceptual approaches and simulation time, the two approaches rendered meaningful results. The link helps forecast the effects of micropore zone transport processes in the subsurface efficiently and thus allows development of numerical tools that could contribute towards more efficient remediation design.     

Downward Migration of Dense Bottom Currents

The focus of this paper is on the dynamics of a dense current flowing along the continental slope, and the frictionally induced downward motion it experiences. In particular, the movement of the lateral boundaries where the isopycnals meet the bottom are considered. The current is taken to be wide compared to the Rossby radius, which is in accordance with observations and makes the dynamics quasigeostrophic. The time development of the plume thickness is studied, using three different parameterisations for the bottom friction.Independently of the choice of parameterisation, the following results are obtained: In the central part of the plume friction acts as a diffusive process to minimise the curvature of the upper surface of the plume. At the upper edge the plume quickly approaches a state with small slope, i.e. small geostrophic velocity and small frictionally induced downward flow. At the lower edge a tongue of dense water shoots out creating a widening layer. The thickness of the migrating layer is approximately equal to the boundary layer depth and its downward speed is comparable to the along-slope geostrophic velocity. The downslope end of the migrating layer may form a steep front which requires some precautions in the numerical procedure.     

A Global Model Tracking Water, Nitrogen, and Land Inputs and Virtual Transfers from Industrialized Meat Production and Trade

Rising populations and incomes throughout the world have boosted meat demand by over 75% in the last 20years, intensifying pressures on production systems and the natural resources to which they are linked. As a growing proportion of global meat production is traded, the environmental impacts of production become increasingly separated from where the meat is consumed. In this paper, we quantify the use of three important resources associated with industrial livestock production and trade—water, land, and nitrogen—using a country-specific model that combines trade, agronomic, biogeochemical, and hydrological data. Our model focuses on pigs and chickens, as these animals are raised predominantly in intensive systems using concentrated, compound feeds. The results describe the geographical patterns of environmental resource use due to meat production, trade, and consumption. We show that US feed, animal, and meat destined for export require almost as much nitrogen and land, and 20% more water, than products destined for domestic consumption. Model results also demonstrate that among various production factors, improvements in crop yields and animal feed conversion efficiencies result in the most significant reductions in environmental harm. By explicitly tracking the externalities of meat production, we hope to bolster suppliers’ accountability and provide better information to meat consumers.

Climate change, urban air problems and transport policy in the European Union

This paper surveys some recent studies on conventional air pollution and climate change in the transport sector in Europe. Fuel efficiency standards, car emission standards and transport pricing instruments are analysed from an economic perspective taking into account environmental and economic efficiency objectives. This revised version was published online in July 2006 with corrections to the Cover Date.