Transport of simazine in unsaturated sandy soil and predictions of its leaching under hypothetical field conditions

The potential contamination of groundwater by herbicides is often controlled by processes in the vadose zone, through which herbicides travel before entering groundwater. In the vadose zone, both physical and chemical processes affect the fate and transport of herbicides, therefore it is important to represent these processes by mathematical models to predict contaminant movement. To simulate the movement of simazine, a herbicide commonly used in Chilean vineyards, batch and miscible displacement column experiments were performed on a disturbed sandy soil to quantify the primary parameters and processes of simazine transport. Chloride (Cl(-)) was used as a non-reactive tracer, and simazine as the reactive tracer. The Hydrus-1D model was used to estimate the parameters by inversion from the breakthrough curves of the columns and to evaluate the potential groundwater contamination in a sandy soil from the Casablanca Valley, Chile. The two-site, chemical non-equilibrium model was observed to best represent the experimental results of the miscible displacement experiments in laboratory soil columns. Predictions of transport under hypothetical field conditions using the same soil from the column experiments were made for 40 years by applying herbicide during the first 20 years, and then halting the application and considering different rates of groundwater recharge. For recharge rates smaller than 84 mm year(-1), the predicted concentration of simazine at a depth of 1 m is below the U.S. EPA's maximum contaminant levels (4 microg L(-1)). After eight years of application at a groundwater recharge rate of 180 mm year(-1) (approximately 50% of the annual rainfall), simazine was found to reach the groundwater (located at 1 m depth) at a higher concentration (more than 40 microg L(-1)) than the existing guidelines in the USA and Europe.     

Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.     

Contaminant degradation in physically and chemically heterogeneous aquifers

This paper examines the importance of the correlation between hydraulic conductivity (K) and degradation rate constant (k) during the transport of reactive contaminants in heterogeneous aquifers. We simulated reactive transport in an ensemble of two-dimensional heterogeneous aquifers. Two sets of transport simulations were conducted: one in which a perfect positive correlation was assumed between ln(K) and ln(k), and one in which a perfect negative correlation was assumed. We found that the sign of the correlation has important consequences for the contaminant transport. Qualitatively, a negative correlation leads to significantly more pronounced "fingering" of the contaminant plume than does a positive correlation, with potentially important consequences for downgradient receptors. Quantitatively, the expected behavior (as quantified by the contaminant mass remaining in the aquifer) is statistically different between the positive and negative cases: on average, more contaminant mass persists when K and k are negatively correlated. Also, the negative correlation leads to more variability between realizations of the ensemble, whereas a positive correlation induces relatively little variability between realizations. We discuss the implications of these findings for the management of contaminated aquifers.     

How Plants Cope with Foreign Compounds. Translocation of xenobiotic glutathione conjugates in roots of barley (Hordeum vulgare) (9 pp)

Background, Aim and Scope Numerous herbicides and xenobiotic organic pollutants are detoxified in plants to glutathione conjugates. Following this enzyme catalyzed reaction, xenobiotic GS-conjugates are thought to be compartmentalized in the vacuole of plant cells. In the present study, evidence is presented for long range transport of these conjugates in plants, rather than storage in the vacuole. To our knowledge this is the first report about the unidirectional long range transport of xenobiotic conjugates in plants and the exudation of a glutathione conjugate from the root tips. This could mean that plants possess an excretion system for unwanted compounds which give them similar advantages as animals. Materials and Methods: Barley plants (Hordeum vulgare L. cv. Cherie) were grown in Petri dishes soaked with tap water in the greenhouse. - Fluorescence Microscopy. Monobromo- and Monochlorobimane, two model xenobiotics that are conjugated rapidly in plant cells with glutathione, hereby forming fluorescent metabolites, were used as markers for our experiments. Their transport in the root could be followed sensitively with very good temporal and spatial resolution. Roots of barley seedlings were cut under water and the end at which xenobiotics were applied was fixed in an aperture with a thin latex foil and transferred into a drop of water on a cover slide. The cover slide was fixed in a measuring chamber on the stage of an inverse fluorescence microscope (Zeiss Axiovert 100). - Spectrometric enzyme assay. Glutathione S-transferase (GST) activity was determined in the protein extracts following established methods. Aliquots of the enzyme extract were incubated with 1-chloro-2,4-dinitrobenzene (CDNB), or monochlorobimane. Controls lacking enzyme or GSH were measured. - Pitman chamber experiments. Ten days old barley plants or detached roots were inserted into special incubation chambers, either complete with tips or decapitated, as well as 10 days old barley plants without root tips. Compartment A was filled with a transport medium and GSH conjugate or L-cysteine conjugate. Compartments B and C contained sugar free media. Samples were taken from the root tip containing compartment C and the amount of conjugate transported was determined spectro-photometrically. Results: The transport in roots is unidirectional towards the root tips and leads to exsudation of the conjugates at rates between 20 and 200 nmol min-1. The microscopic studies have been complemented by transport studies in small root chambers and spectroscopic quantification of dinitrobenzene-conjugates. The latter experiments confirm the microscopic studies. Furthermore it was shown that glutathione conjugates are transported at higher rates than cysteine conjugates, despite of their higher molecular weights. This observation points to the existence of glutathione specific carriers and a specific role of glutathione in the root. Discussion: It can be assumed that long distance transport of glutathione conjugates within the plant proceeds like GSH or amino acid transport in both, phloem and xylem. The high velocity of this translocation of the GS-X is indicative of an active transport. For free glutathione, a rapid transport-system is essential because an accumulation of GSH in the root tip inhibits further uptake of sulfur. Taking into account that all described MRP transporters and also the GSH plasmalemma ATPases have side activities for glutathione derivatives and conjugates, co-transport of these xenobiotic metabolites seems credible. - On the other hand, when GS-B was applied to the root tips from the outside, no significant uptake was observed. Thus it can be concluded that only those conjugates can be transported in the xylem which are formed inside the root apex. Having left the root once, there seems to be no return into the root vessels, probably because of a lack of inward directed transporters. Conclusions: Plants seem to possess the capability to store glutathione conjugates in the vacuole, but under certain conditions, these metabolites might also undergo long range transport, predominantly into the plant root. The transport seems dependent on specific carriers and is unidirectional, this means that xenobiotic conjugates from the rhizosphere are not taken up again. The exudation of xenobiotic metabolites offers an opportunity to avoid the accumulation of such compounds in the plant. Recommendations and Perspectives: The role of glutathione and glutathione related metabolites in the rhizosphere has not been studied in any detail, and only scattered data are available on interactions between the plant root and rhizosphere bacteria that encounter such conjugates. The final fate of these compounds in the root zone has also not been addressed so far. It will be interesting to study effects of the exuded metabolites on the biology of rhizosphere bacteria and fungi.     

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

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

An approach to optimise nutrient management in environmental sanitation systems despite limited data

The material flow analysis method can be used to assess the impact of environmental sanitation systems on resource consumption and environmental pollution. However, given the limited access to reliable data, application of this data-intensive method in developing countries may be difficult. This paper presents an approach allowing to develop material flow models despite limited data availability. Application of an iterative procedure is of key importance: model parameter values should first be assessed on the basis of a literature review and by eliciting expert judgement. If model outputs are not plausible, sensitive input parameters should be reassessed more accurately. Moreover, model parameters can be expressed as probability distributions and variable uncertainty estimated by using Monte Carlo simulation. The impact of environmental sanitation systems on the phosphorus load discharged into surface water in Hanoi, Vietnam, is simulated by applying the proposed approach.     

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.     

Elimination of Methanol and Ethanol by Biofiltration: An Experimental Study

This study presents the experimental results obtained during long-term operation of two biofilters treating two alcohols: methanol and ethanol. The biofilters used for this purpose were previously packed with a compost material made from tobacco processing residues. The alcohols concentrations tested lay between 0.40 and 3.20 g/m³ for methanol, and 0.55 and 5.05 g/m³ for ethanol. The empty bed residence time in each biofilter was 60 s. Biofilter inlet loads of less than 190 and 300 g/m³/h for the methanol and ethanol additions respectively, were thereafter evaluated. In addition, the concentrations of nutrient nitrogen were also varied, from 0.1 to 2.0 g-N/l and from 0.3 to 11.3 g-N/l for the ethanol and the methanol, respectively. The results thus obtained have made it possible to select the optimal nitrogen concentrations which, for the cases examined, turn out to be 0.3 g-N/l and between 2 and 3.8 g-N/l for the ethanol and methanol substrates, respectively. The maximum elimination capacities obtained in this study were 82 and 150 g/m³/h, respectively for the methanol and ethanol cases. It was therefore concluded that, for a readily biodegradable compound such as ethanol, the nitrogen requirement is substantially lower than that needed for the methanol degradation, the latter appearing to be more difficult to degrade biologically under similar operating conditions. The production rate of the co-product carbon dioxide during methanol and ethanol biofiltration was also investigated. Also, a good correlation was found to exist between the temperature and the conversion achieved in the biofilter.     

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.     

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.