Long-term oxygen depletion from infiltrating groundwaters: model development and application to intra-glaciation and glaciation conditions

Processes that control the redox conditions in deep groundwaters have been studied. The understanding of such processes in a long-term perspective is important for the safety assessment of a deep geological repository for high-level nuclear waste. An oxidising environment at the depth of the repository would increase the solubility and mobility of many radionuclides, and increase the potential risk for radioactive contamination at the ground surface. Proposed repository concepts also include engineered barriers such as copper canisters, the corrosion of which increases considerably in an oxidising environment compared to prevailing reducing conditions. Swedish granitic rocks are typically relatively sparsely fractured and are best treated as a dual-porosity medium with fast flowing channels through fractures in the rock with a surrounding porous matrix, the pores of which are accessible from the fracture by diffusive transport. Highly simplified problems have been explored with the aim to gain understanding of the underlying transport processes, thermodynamics and chemical reaction kinetics. The degree of complexity is increased successively, and mechanisms and processes identified as of key importance are included in a model framework. For highly complex models, analytical expressions are not fully capable of describing the processes involved, and in such cases the solutions are obtained by numerical calculations. Deep in the rock the main source for reducing capacity is identified as reducing minerals. Such minerals are found inside the porous rock matrix and as infill particles or coatings in fractures in the rock. The model formulation also allows for different flow modes such as flow along discrete fractures in sparsely fractured rocks and along flowpaths in a fracture network. The scavenging of oxygen is exemplified for these cases as well as for more comprehensive applications, including glaciation considerations. Results show that chemical reaction kinetics control the scavenging of oxygen during a relatively short time with respect to the lifetime of the repository. For longer times the scavenging of oxygen is controlled by transport processes in the porous rock matrix. The penetration depth of oxygen along the flowpath depends largely on the hydraulic properties, which may vary significantly between different locations and situations. The results indicate that oxygen, in the absence of easily degradable organic matter, may reach long distances along a flow path during the life-time of the repository (hundreds to thousands of metres in a million years depending on e.g. hydraulic properties of the flow path and the availability of reducing capacity). However, large uncertainties regarding key input parameters exist leading to the conclusion that the results from the model must be treated with caution pending more accurate and validated data. Ongoing and planned experiments are expected to reduce these uncertainties, which are required in order to make more reliable predictions for a safety assessment of a nuclear waste repository.     

Linking a one-dimensional pesticide fate model to a three-dimensional groundwater model to simulate pollution risks of shallow and deep groundwater underlying fractured till

The one-dimensional pesticide fate model MACRO was loose-linked to the three-dimensional discrete fracture/matrix diffusion model FRAC3DVS to describe transport of the pesticide mecoprop in a fractured moraine till and local sand aquifer (5-5.5 m depth) overlying a regional limestone aquifer (16 m depth) at Havdrup, Denmark. Alternative approaches to describe the upper boundary in the groundwater model were examined. Field-scale simulations were run to compare a uniform upper boundary condition with a spatially variable upper boundary derived from Monte-Carlo simulations with MACRO. Plot-scale simulations were run to investigate the influence of the temporal resolution of the upper boundary conditions for fluxes in the groundwater model and the effects of different assumptions concerning the macropore/fracture connectivity between the two models. The influence of within-field variability of leaching on simulated mecoprop concentrations in the local aquifer was relatively small. A fully transient simulation with FRAC3DVS gave 20 times larger leaching to the regional aquifer compared to the case with steady-state water flow, assuming full connectivity with respect to macropores/fractures across the boundary between the two models. For fully transient simulations 'disconnecting' the macropores/fractures at the interface between the two models reduced leaching by a factor 24. A fully connected, transient simulation with FRAC3DVS, with spatially uniform upper boundary fluxes derived from a MACRO simulation with 'effective' parameters is therefore recommended for assessing leaching risks to the regional aquifer, at this, and similar sites.     

Fast method for simulation of radionuclide chain migration in dual porosity fracture rocks

In fractured rocks with a porous rock matrix such as granites, radionuclides will flow with the water in the fracture network. The nuclides will diffuse in and out the rock matrix where they can sorb and be considerably retarded compared to the water velocity. A water parcel entering the network will mix and split at the fracture intersections and parts of the original parcel will traverse a multitude of different fractures. The flowrates, velocities, sizes and apertures of the fractures can vary widely. Normally one must solve the transport equations for every fracture and use the effluent concentration as inlet condition to the next fracture and so on. It is shown that under some weakly simplified conditions it suffices to determine one single parameter group containing information on the flow wetted surface that a water parcel contacts along the entire path. It is also shown how this can be obtained. Then, solving the transport equations only once for time and location along the path gives the concentration and nuclide flux of every nuclide in the chain everywhere along a path. The same solution actually is valid for every path in the network. This dramatically reduces the computation effort. The same approach can be used for models based on streamtubes.     

Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms

Recent discovery of bomb-related ³⁶Cl at depth in fractured tuff in the unsaturated zone at the Yucca Mountain candidate high-level waste (HLW) repository site has called into question the usual modeling assumptions based on the equivalent continuum model (ECM). A dual continuum model (DCM) for simulating transient flow and transport at Yucca Mountain is developed. In order to ensure properly converged flow solutions, which are used in the transport simulation, a new flow solution convergence criteria is derived. An extensive series of simulation studies is presented which indicates that rapid movement of solute through the fractures will not occur unless there are intense episodic infiltration events. Movement of solute in the environs of the repository is enhanced if the properties of the tuff layer at the repository horizon are modified from current best-estimate values. Due to a large advective–dispersive coupling between the matrix and fractures, the matrix acts as a major buffer which inhibits rapid transport along the fractures. Consequently, fast movement of solutes through the fractures to the repository depth can only be explained if the matrix–fracture coupling term is significantly reduced from a value that would be calculated on the basis of data currently available.     

The effect of a biofilm on solute diffusion in fractured porous media

At sites in fractured rock where contamination has been exposed to the rock matrix for extended periods of time, the amount of contaminant mass residing in the matrix can be considerable. Even though it may be possible to diminish concentrations by the advection of clean water through the fracture features, back diffusion from mass held in the matrix will lead to a continuing source of contamination. In such an event, the development of a biofilm (a thin film of microbial mass) on the wall of the fractures may act to limit or prevent the back diffusion process. The objective of this preliminary study is to explore the influence imparted by the presence of a biofilm on the process of matrix diffusion. The investigation was conducted using radial diffusion cells constructed from rock core in which biofilm growth was stimulated in a central reservoir. Once biofilms were developed, forward diffusion experiments were conducted in which a conservative solute migrated from the central reservoir into the intact rock sample. Diffusion experiments were performed in a total of 11 diffusion cell pairs where biofilm growth was stimulated in one member of the pair and inhibited in the other. The effect of the presence of a biofilm on tracer diffusion was determined by comparison of the diffusion curves produced by each cell pair. A semi-analytical model that accounts for the presence of a biofilm was used to investigate the effect of the biofilm on mass transfer due to changes in the effective porosity, effective diffusion coefficient, and the depth of penetration of the biofilm into the intact rock. The results show that the biofilm acted to plug the rock matrix, rather than forming a discrete layer on the reservoir surface. The reduction in effective porosity due to the biofilm ranged from 6% to 52% with the majority of the samples in the 30% to 50% range. Based on the present results, with more efficient biofilm stimulation, it is reasonable to assume that a more complete plugging of the microcrack porosity might be possible, leaving a much thicker and efficient barrier than could be achieved via a surface biofilm.     

Tritium and 36Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and 36Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this "fast flow" in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891-989] and suggest that fast flow in fractures with minimal fracture-matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and 36Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Numerical analysis of water and solute transport in variably-saturated fractured clayey till

This study numerically investigates the influence of initial water content and rain intensities on the preferential migration of two fluorescent tracers, Acid Yellow 7 (AY7) and Sulforhodamine B (SB), through variably-saturated fractured clayey till. The simulations are based on the numerical model HydroGeoSphere, which solves 3D variably-saturated flow and solute transport in discretely-fractured porous media. Using detailed knowledge of the matrix, fracture, and biopore properties, the numerical model is calibrated and validated against experimental high-resolution tracer images/data collected under dry and wet soil conditions and for three different rain events. The model could reproduce reasonably well the observed preferential migration of AY7 and SB through the fractured till, although it did not capture the exact depth of migration and the negligible impact of the dead-end biopores in a near-saturated matrix. A sensitivity analysis suggests fast flow mechanisms and dynamic surface coating in the biopores, and the presence of a plough pan in the till.     

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).     

Biogenic fluxes of halomethanes from Irish peatland ecosystems

Irish peatland ecosystems have been shown to be important sources of low molecular weight halocarbons. Emission of CH₃Br, CH₃Cl, CH₃I and CHCl₃ was recorded from all peatland sites monitored, with minor flux of other halocarbons at certain sites. Fluxes were found to be highly linked to incident light, with strong diurnal cycles recorded at all open peatland sites. Estimates of halomethane emissions, particularly from coastal peatland and conifer plantation forest floor sites, suggests that these ecosystems may make a significant contribution to the global budgets of several important halocarbons. Global annual fluxes of 4.7 (0.1–151.9), 0.9 (0.1–3.3), 5.5 (0.9–43.4), and 1.4 (0.1–12.8) Gg yr⁻¹ for CHCl₃, CH₃Br, CH₃Cl, and CH₃I, respectively, were determined for peatland ecosystems.     

Tritium and Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and ³⁶Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this “fast flow” in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891–989] and suggest that fast flow in fractures with minimal fracture–matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and ³⁶Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Concentrations of chlorine, bromine and iodine in Japanese rivers

Concentrations of halogens (Cl, Br and I) in 30 Japanese rivers were measured by ion chromatography and inductively coupled plasma mass spectrometry to understand their behavior in the terrestrial environment. Concentrations of Cl, Br and I in each river, obtained at 10 sampling points from the upper stream to the river mouth, tended to increase near the river mouth. The ranges of geometric means of Cl, Br and I in each river were 1.0–19.4 mg l⁻¹, 2.5–67.9 μg l⁻¹, and 0.18–8.34 μg l⁻¹, respectively. To compare halogen behavior, the concentration ratios, Br/Cl and I/Cl, were calculated. The Br/Cl range was (2.3–7.8) × 10⁻³ (geometric mean: 3.74 × 10⁻³), and it was nearly constant except for the Yoneshiro river. It was estimated that 60–80% of total Br in the middle to lower parts of this river was the excess Br. The Br chemical form in all the rivers is generally considered to be Br⁻. The I/Cl ratios had different trends in rivers flowing into the Japan Sea and Pacific Ocean, possibly due to the different geological features in the river catchments.     

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21 m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3 mx4 m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten alpha parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns.     

Hydraulic characterization for steam enhanced remediation conducted in fractured rock

To explore the viability of Steam Enhanced Remediation (SER) in fractured rock a small-scale steam injection and water/vapour extraction pilot study was conducted at the former Loring Air Force Base in northern Maine, USA. A detailed well testing program was undertaken to assist in the design of the injection and extraction well array, and to assess the possibility of off-site heat and contaminant migration. A structurally complex limestone having low matrix porosity and a sparse distribution of fractures underlies the study site. To characterize the groundwater and steam flow pathways, single-well slug tests and more than 100 pulse interference tests were conducted. The results of the well testing indicate that the study site is dominated by steeply dipping bedding plane fractures that are interconnected only between some wells in the injection/extraction array. The SER system was designed to take advantage of interconnected fractures located at depth in the eastern end of the site. An array of 29 wells located in an area of 60 by 40 m was used for steam injection and water/vapour extraction. The migration of heat was monitored in several wells using thermistor arrays having a 1.5 m vertical spacing. Temperature measurements obtained during and after the 3 month steam injection period showed that heat migration generally occurred along those fracture features identified by the pulse interference testing. Based on these results, it is concluded that the pulse interference tests were valuable in assisting with the design of the injection/extraction well geometry and in predicting the migration pathways of the hot water associated with the steam injection. The pulse interference test method should also prove useful in support of any other remedial method dependant on the fracture network for delivery of remedial fluid or extraction of contaminants.     

Identification of highly brominated analogues of Q1 in marine mammals

Three novel halogenated organic compounds (HOCs) have been identified in the blubber of marine mammals from coastal New England with the molecular formulae C(9)H(3)N(2)Br(6)Cl, C(9)H(3)N(2)Br(7), and C(9)H(4)N(2)Br(5)Cl. They were identified using high and low resolution gas chromatography mass spectrometry (GCMS) and appear to be highly brominated analogues of Q1, a heptachlorinated HOC suspected to be naturally produced. These compounds were found in Atlantic white sided dolphin (Lagenorhynchus acutus), bottlenose dolphin (Tursiops truncatus), common dolphin (Delphinus delphis), Risso's dolphin (Grampus griseus), harbor porpoise (Phocoena phocoena), beluga whale (Delphinapterus leucas), fin whale (Balaenoptera physalus), grey seal (Halichoerus grypus), harp seal (Phoca groenlandica) and a potential food source (Loligo pealei) with concentrations as high as 2.7 microg/g (lipid weight). The regiospecificity of C(9)H(3)N(2)Br(6)Cl is suggestive of a biogenic origin. Debromination of C(9)H(3)N(2)Br(6)Cl may be significant in the formation of C(9)H(4)N(2)Br(5)Cl.     

Sense or no-sense of the sum parameter for water soluble "adsorbable organic halogens" (AOX) and "absorbed organic halogens" (AOX-S18) for the assessment of organohalogens in sludges and sediments

"AOX" is the abbreviation of the sum parameter for water soluble "adsorbable organic halogens" in which 'A' stands for adsorbable, 'O' for organic and 'X' for the halogens chlorine, bromine and iodine.After the introduction of the AOX in 1976, this parameter has been correctly used for "real" AOX constituents (DDT and its metabolites, PCBs, etc.) but also misused for non-adsorbable adsorbed OX-compounds, mostly high molecular organohalogens in plants and even to inorganic compounds being neither organic nor adsorbable.The question of natural "Adsorbable Organic Halogens" (AOX) formed by living organisms and/or during natural abiogenic processes has been definitively solved by the known existence of already more than 3650 organohalogen compounds, amongst them the highly reactive, cancerogenic vinyl chloride (VC).The extension of the AOX to AOX-S18 for Sludges and Sediments, in which A stands for adsorbed (not for adsorbable) is questionable. It includes the most important water insoluble technical organochlorine product: polyvinyl chloride, PVC.In addition to organic halogens it also includes inorganic, mineralogenic halides, incorporated mainly in the crystal lattice of fine grained phyllosilicates, the typical clay minerals (kaolinite, montmorillonite, illite and chlorite) which are main constituents of sediments and sedimentary rocks representing the major part of the sedimentary cover of the earth.Other phyllosilicates, biotite and muscovite, major constituents of granites and many metamorphic rocks (gneiss and mica schist) will also contribute to the AOX-S18 especially in soils as result of weathering processes.Since chlorine is incorporated into the mineral structure and, as a consequence, not soluble by the nitric acid analytical step (pH 0.5) of the S18 determination, it will account to the AOX-S18 in the final charcoal combustion step at temperatures >950 degrees C.After heavy rainfalls sewage sludge composition is strongly influenced by mineralogenic components derived from the erosion of fine grained sediments or soils. Assuming 50% geogenic particles with a mean Cl concentration of 103 mg/kg (as in shales and clays) the mineralogenic Cl-content could add about 50 mg/kg to the organic AOX in sewage sludge.The occurrence of insoluble and non-adsorbable PVC in sewage sludge exhibits the same problems as the mineralogenic constituents: a detection as AOX-S18 is possible when the final high temperature analytical step is applied.Plants as major sources of organohalogens have never been doubted.Only recently [Science 295 (2002) 985] based on the determination of the form of Cl with near-edge X-ray adsorption fine structure (NEXAFS) spectroscopy and extended X-ray adsorption showed the variations in the inorganic Cl(-) and organo-Cl compounds with increasing humification of plant leaves from "fresh leaves--senescent leaves on plants--senescent leaves on soil--powdered top soil--isolated soil humus". His finding of exclusively inorganic Cl(-) in the starting material (fresh leaves) is controverse to our earlier results indicating the presence of ionic inorganic Cl together with water insoluble absorbed organohalogens (AOX-S18) in eight different macrophytes of both terrestrial and marine environments.Our research on AOX in interstitial water of anaerobic limnic sediments has led to the role of bromine playing in the diagenesis of the organic matter of sediments. In sediments of Lake Constance Br(-) concentrations in lake water at the sediment water interface increased from <0.01 to 0.25 mg/l in the pore water at 77 cm sediment depth.In the Neckar River a Br concentrations of 0.02 mg/l at the water/sediment interface increasing to 0.74 mg/l in pore water in 85 cm depth was found. Here a parallel development could be found with ammonium concentration and alkalinity. The very high positive correlation ammonium:bromide and bromide:alkalinity leads to the conclusion, that bromine, originally a high molecular constituent of the organic matter, is released as bromide during an early dehalogenation stage of diagenesis.The mlusion, that bromine, originally a high molecular constituent of the organic matter, is released as bromide during an early dehalogenation stage of diagenesis.The main general reason to discard the AOX sum parameter as a whole lies in the fact, that adsorbable halogenated substances cannot a priori be categorized as natural/anthropogenic, biotic/abiotic, harmful/harmless. If applied to sludges and sediments, adsorbed organohalogens are not water soluble and therefore not adsorbable, and mineralogenic halogens (X) are neither organic nor adsorbable, and therefore by definition no AOX.     

Monobromo and higher brominated congeners of the marine halogenated natural product 2,3,3',4,4',5,5'-heptachloro-1'-methyl-1,2'-bipyrrole (Q1)

The marine halogenated natural product 2,3,3',4,4',5,5'-heptachloro-1'-methyl-1,2'-bipyrrole (Q1) is widely distributed in the environment. In this study, we screened samples which have previously been found to contain remarkably high residues of Q1 (blubber of marine mammals from Australia, samples from Antarctica, human milk from the Faroe Island) for the additional presence of mixed chlorinated and brominated congeners. Using GC/ECNI-MS, all samples tested were positive and many contained four out of five possible bromohexachloro congeners (BrCl6-MBPs), five out of 14 possible dibromopentachloro congeners (Br2Cl5-MBPs), five of 21 possible tribromotetrachloro-congeners (Br3Cl4-MBPs), as well as several higher brominated congeners. About 20 heptahalo congeners of Q1 are described for the first time in the scientific literature. Isomers eluted within about one minute, respectively. Hence it is possible, that the peak clusters identified may be composed of more, co-eluting congeners. Similarities in the GC/ECNI-MS mass spectra with polychlorinated biphenyls (PCBs) were addressed. We also suggest an acronym system similar to that in use for polychlorinated biphenyls that may simplify the use of this substance class in scientific papers. In the samples from Australia, BrCl6-MBPs and Br2Cl5-MBPs amounted for 7-27.5% and 0.4-4.2% of Q1, respectively whereas Br3Cl4-MBPs and higher brominated MBPs were found in the range of <1% of Q1 or less.     

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.     

Mass removal of chlorinated ethenes from rough-walled fractures using permanganate

In situ chemical oxidation (ISCO) employing permanganate is an emerging technology that has been successful at enhancing mass removal from DNAPL source zones in unconsolidated media at the pilot-scale. The focus of this study was to evaluate the applicability of flushing a permanganate solution across two single vertical fractures in a laboratory environment to remove free phase DNAPL. The fracture experiments were designed to represent a portion of a larger fractured aquifer system impacted by a near-surface DNAPL spill over a shallow fractured rock aquifer. Each fracture was characterized by hydraulic and tracer tests, and the aperture field for one of the fractures was mapped using a co-ordinate measurement machine. Following DNAPL emplacement, a series of water and permanganate flushes were performed. To support observations from the fracture experiments, a set of batch experiments was conducted. The data from both fracture experiments showed that the post-oxidation effluent concentration was not impacted by the oxidant flush; however, changes in the aperture distribution, flow field, and flow rate were observed. These changes resulted in a significant decrease to the mass loading from the fractures, and were attributed to the build-up of oxidation by-products (manganese oxides and carbon dioxide) within the fracture which was corroborated by the batch experiment data and visual examination of the walls of one fracture. These results provide insight into the potential impact that a permanganate solution and oxidation by-products can have on the aperture distribution within a fracture and on DNAPL mass transfer rates. A permanganate flush or injection completed within a fractured rock aquifer may lead to the development of an insoluble product adjacent to the DNAPL which results in the reduction or complete elimination of advective regions near the DNAPL and reduces mass transfer rates. This outcome would have significant implications on the plume generating potential of the remaining DNAPL.     

Strontium migration in a crystalline medium: effects of the presence of bentonite colloids

The effects of bentonite colloids on strontium migration in fractured crystalline medium were investigated. We analyzed first the transport behaviour of bentonite colloids alone at different flow rates; then we compared the transport behaviour of strontium as solute and of strontium previously adsorbed onto stable bentonite colloids at a water velocity of approximately 7.1·10(-6)m/s-224m/yr. Experiments with bentonite colloids alone showed that - at the lowest water flow rate used in our experiments (7.1·10(-6)m/s) - approximately 70% of the initially injected colloids were retained in the fracture. Nevertheless, the mobile colloidal fraction, moved through the fracture without retardation, at any flow rate. Bentonite colloids deposited over the fracture surface were identified during post-mortem analyses. The breakthrough curve of strontium as a solute, presented a retardation factor, R(f)~6, in agreement with its sorption onto the granite fracture surface. The breakthrough curve of strontium in the presence of bentonite colloids was much more complex, suggesting additional contributions of colloids to strontium transport. A very small fraction of strontium adsorbed on mobile colloids moved un-retarded (R(f)=1) and this fraction was much lower than the expected, considering the quantity of strontium initially adsorbed onto colloids (90%). This behaviour suggests the hypothesis of strontium sorption reversibility from colloids. On the other hand, bentonite colloids retained within the granite fracture played a major role, contributing to a slower strontium transport in comparison with strontium as a solute. This was shown by a clear peak in the breakthrough curve corresponding to a retardation factor of approximately 20.     

Trace elements in atmospheric particulate matter over a coal burning power production area of western Macedonia, Greece

Total suspended particle (TSP) concentrations were determined in the Eordea basin (western Macedonia, Greece), an area with intensive lignite burning for power generation. The study was conducted over a one-year period (November 2000–November 2001) at 10 sites located at variable distances from the power plants. Ambient TSP samples were analyzed for 27 major, minor and trace elements. Annual means of TSP concentrations ranged between 47 ± 33 μg m⁻³ and 110 ± 50 μg m⁻³ at 9 out of the 10 sites. Only the site closest to the power stations and the lignite conveyor belts exhibited annual TSP levels (210 ± 97 μg m⁻³) exceeding the European standard (150 μg m⁻³, 80/779/EEC). Concentrations of TSP and almost all elemental components exhibited significant spatial variations; however, the elemental profiles of TSP were quite similar among all sites suggesting that they are affected by similar source types. At all sites, statistical analysis indicated insignificant (P < 0.05) seasonal variation for TSP concentrations. Some elements (Cl, As, Pb, Br, Se, S, Cd) exhibited significantly higher concentrations at certain sites during the cold period suggesting more intense emissions from traffic, domestic heating and other combustion sources. On the contrary, concentrations significantly higher in the warm period were found at other sites mainly for crustal elements (Ti, Mn, K, P, Cr, etc.) suggesting stronger influence from soil resuspension and/or fly ash in the warm months. The most enriched elements against local soil or road dust were S, Cl, Cu, As, Se, Br, Cd and Pb, whereas negligible enrichment was found for Ti, Mn, Mg, Al, Si, P, Cr. At most sites, highest concentrations of TSP and elemental components were associated with low- to moderate-speed winds favoring accumulation of emissions from local sources. Influences from the power generation were likely at those sites located closest to the power plants and mining activities.