Dioxins,chlorophenols and other chlorinated organic pollutants in colloidal and water fractions of groundwater from a contaminated sawmill site

BACKGROUND, AIM, AND SCOPE: The distribution of chlorinated organic contaminants in groundwater and the importance of colloids were studied in groundwater from a sawmill site contaminated by chlorophenol preservatives. MATERIALS AND METHODS: The groundwater was fractionated into three different size ranges: (1) >0.7 mum, (2) 0.4-0.7 mum and (3) 0.2-0.4 mum and the filtered water phase. The concentrations of chlorophenols (CP), chlorinated phenoxy phenols (PCPP), chlorinated diphenyl ethers (PCDE), chlorinated dibenzofurans (PCDF) and chlorinated dibenzo-p-dioxins (PCDD) were determined in each fraction. The colloids were characterised regarding the chemical composition using X-ray photoelectron spectroscopy (XPS). RESULTS: Chlorophenols were mostly found in the water fraction and PCDD/Fs were found almost exclusively in the particulate fractions. For example, the filtered water phase contained 2,100 mug l(-1) and 0.72 ng l(-1) for CPs and PCDD/Fs, respectively, and the particulate fractions contained 27 mug l(-1) and 32 ng l(-1) for CPs and PCDD/Fs, respectively. XPS evaluation of the particulate phases showed no correlation between the surface chemistry of the particle properties and the distribution of chlorinated compounds. DISCUSSION: The results suggest that groundwater transport of CPs, PCPPs, PCDEs and PCDD/Fs may occur from contaminated sawmill sites and that the colloid-facilitated transport, especially of PCDD/Fs, is substantial. The results correlated well with previous studies of compounds sorbed to dissolved organic carbon, which indicate that dissolved and colloidal organic carbon facilitated the transport of PCDEs, PCDFs and PCDDs particularly. CONCLUSIONS: Several classes of chlorinated compounds were readily detected in the groundwater samples. Due to the differences in their physicochemical properties, CPs, PCPPs, PCDEs and PCDD/Fs vary in their partitioning between colloidal fractions and the filtered groundwater. The proportion of the bound fraction increased with an increasing hydrophobicity of the chlorinated compounds. The groundwater transport of colloid-associated pollutants from the site may be significant. RECOMMENDATIONS AND PERSPECTIVES: The results imply that colloidal particles <0.7 mum are freely mobile in groundwater from this site. The groundwater transport of colloid-associated pollutants may be significant. However, the extent of the problem is not yet known and, thus, further research is needed to evaluate the impact of colloidal transport of hydrophobic organic contaminants. In Sweden alone, 400 to 500 sawmill sites are estimated to be contaminated with PCDD/Fs as a result of the former use of CP-based wood preservatives. The widespread use of CP mixtures for a variety of applications, including wood preservation, indicates that potential colloidal transport will be an issue of concern in many countries.     

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.     

特大型城市客运交通碳排放与减排对策研究

本文基于对现有城市交通碳排放测算方法的比较分析,以上海市为例,采用IPCC"自下而上"法对特大型城市客运交通CO2排放进行了测算,结果显示:轨道交通是碳排放效率最高的客运方式,出租车最低;客运交通CO2排放总量增长迅速,且碳源结构发生了较大变化;近年客运交通CO2排放增量主要来自私人载客汽车,同时公务交通在客运交通碳排放中始终占较大比重。由此本文认为,控制客运交通碳排放的关键在于对以私人载客汽车和单位载客汽车为主的个体交通的管理和控制,形成以公共交通为主的交通结构。在此基础上,为了将控制碳排放纳入到城市交通政策目标中去,本文就主要城市交通政策对客运交通碳排放产生的影响进行了深入分析,并得出结论:以往的交通供给、需求管理政策对于抑制客运交通碳排放增长的作用有限;而就目前城市空间发展政策的实施效果而言,也不利于降低居民出行的碳排放水平。文章最后分别从交通供给、需求管理以及城市空间角度给出了控制客运交通碳排放的对策。     

Bacteria transport and deposition under unsaturated conditions: the role of the matrix grain size and the bacteria surface protein

Unsaturated (80% water saturated) packed column experiments were conducted to investigate the influence of grain size distribution and bacteria surface macromolecules on bacteria (Rhodococcus rhodochrous) transport and deposition mechanisms. Three sizes of silica sands were used in these transport experiments, and their median grain sizes were 607, 567, and 330 microm. The amount of retained bacteria increased with decreasing sand size, and most of the deposited bacteria were found adjacent to the column inlet. The deposition profiles were not consistent with predictions based on classical filtration theory. The experimental data could be accurately characterized using a mathematical model that accounted for first-order attachment, detachment, and time and depth-dependent straining processes. Visual observations of the bacteria deposition as well as mathematical modelling indicated that straining was the dominant mechanism of deposition in these sands (78-99.6% of the deposited bacteria), which may have been enhanced due to the tendency of this bacterium to form aggregates. An additional unsaturated experiment was conducted to better deduce the role of bacteria surface macromolecules on attachment and straining processes. In this case, the bacteria surface was treated using a proteolitic enzyme. This technique was assessed by examining the Fourier-transform infrared spectrum and hydrophobicity of untreated and enzyme treated cells. Both of these analytical procedures demonstrated that this enzymatic treatment removed the surface proteins and/or associated macromolecules. Transport and modelling studies conducted with the enzyme treated bacteria, revealed a decrease in attachment, but that straining was not significantly affected by this treatment.     

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.     

Semi-analytical solution of two-dimensional sharp interface LNAPL transport models

In this paper, we present semi-analytical solutions for two-dimensional equations governing transport of Light Non-Aqueous Phase Liquids (LNAPL) in unconfined aquifers. The proposed model is based on sharp interface displacement and steady groundwater flow assumptions, where both the water–LNAPL interface and the LNAPL–air interface are represented as sharp interfaces. In the case of steady groundwater flow, these equations can be reduced to a two-dimensional nonlinear solute transport equation, with the LNAPL thickness in the free product lens being the primary unknown variable. The linearized form of this solute transport equation falls into the category of two-dimensional transport equation with time-dependent dispersion coefficients. This equation can be solved analytically for an infinite domain region. In this paper, the general form of the analytical solution for the transport equation, as well as the solutions for some specific cases are presented. To demonstrate the utility of the proposed solution, numerical results obtained for two example problems are discussed and presented comparatively with a finite-element solution and other more restrictive solutions available in the literature. Although the solutions discussed in this paper have some simplifying assumptions, such as sharp-interfaces between fluid phases, steady groundwater flow and homogeneous aquifer properties, the semi-analytical solutions presented in this study may be used effectively as bench mark solutions in evaluating LNAPL migration in the subsurface. These solutions are simple and cost effective to implement and may be used in the calibration of other more complex numerical solutions that can be found in the literature.     

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.     

Dairy washwater treatment using a horizontal flow biofilm system

In Ireland, dairy farmyard washwater commonly comprises farmyard run-off and dairy parlour washings. Land-spreading is the most widely used method for treating this wastewater. However, this method can be labour intensive and can cause, in some cases, the degradation of surface and ground waters, mainly due to nitrogen contamination. In this study, a horizontal flow biofilm reactor (HFBR) with step-feed was constructed and tested in the laboratory, to remove organic carbon and nitrogen from a agricultural strength synthetic washwater (SWW). The HFBR had an average top plan surface area (TPSA) of 0.1002 m(2) and consisted of a stack of 45 polystyrene horizontal sheets--15 sheets embedded with 25 mm deep frustums above 30 sheets with 10 mm deep frustums. The frustums acted as miniature reservoirs. The sheets were alternately offset to allow the wastewater to flow horizontally along each sheet and vertically from sheet to sheet down through the reactor. Biofilms developed on the sheets and treated the wastewater. During the 212-d study, the total hydraulic loading rate based on the TPSA of the sheets was 35 l m(-2) d(-1). SWW was pumped for 10 min each hour, in a step feed arrangement at a rate of 23.33 l m(-2) d(-1) on to the top sheet during Phases 1 and 2, and 11.67 l m(-2) d(-1) onto Sheet 16 during Phase 1 (days 1-92) and onto Sheet 30 during Phase 2 (days 93-212). The substrate loading rate during Phases 1 and 2 was 94.8 g total chemical oxygen demand (COD) m(-2) d(-1) and 10.5 g total nitrogen (TN) m(-2) d(-1), based on the TPSA. At steady state in Phase 2, the unit achieved excellent carbon removal of 99.7% 5-day biochemical oxygen demand (BOD(5)) and 96.7% total COD, equivalent to TPSA removal rates of 67.5 g BOD(5)m(-2)d(-1) and 91.7 g COD m(-2) d(-1). The nitrogen removal percentages were 98.3% total ammonium-nitrogen (NH(4)-N(t)) and 72.8% TN, which equated to TPSA removal rates of 4.8 g NH(4)-N(t) m(-2) d(-1) and 7.6g TN m(-2) d(-1). No sloughing of solids or clogging of media occurred during the study. The unit was simple to construct and operate, with little maintenance.     

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.     

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.