Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

Hydrogeochemical controls on uranium in aquifers of the Jacobsville Sandstone

The occurrence of elevated uranium (U) in sandstone aquifers was investigated in the Upper Peninsula of Michigan, focusing on aquifers of the Jacobsville Sandstone. The hydrogeochemical controls on groundwater U concentrations were characterized using a combination of water sampling and spectral gamma-ray logging of sandstone cliffs and residential water wells. 235U/238U isotope ratios were consistent with naturally occurring U. Approximately 25% of the 270 wells tested had U concentrations above the U.S. Environmental Protection Agency Maximum Contaminant Level (MCL) of 30 μg/L, with elevated U generally occurring in localized clusters. Water wells were logged to determine whether groundwater U anomalies could be explained by the heterogeneous distribution of U in the sandstone. Not all wells with relative U enrichment in the sandstone produced water with U above the MCL, indicating that the effect of U enrichment in the sandstone may be modified by other hydrogeochemical factors. Well water had high redox, indicating U is in its highly soluble (VI) valence. Equilibrium modeling indicated that aqueous U is complexed with carbonates. In general, wells with elevated U concentrations had low 235U/238U activity ratios. However, in some areas U concentrations and 235U/238U activity ratios were simultaneously high, possibly indicating differences in rock-water interactions. Limited groundwater age dating suggested that residence time may also help explain variations in well water U concentrations. Low levels of U enrichment (4 to 30 ppm) in the Jacobsville sandstone may make wells in its oxidized aquifers at risk for U concentrations above the MCL. On average, U concentrations were highest in heavy mineral and clay layers and rip up conglomerates. Uranium concentrations above 4 ppm also occurred in siltstones, sandstones and conglomerates. Uranium enrichment was likely controlled by deposition processes, sorption to clays, and groundwater flow, which was controlled by permeability variations in the sandstone. Low levels of U enrichment were found at deltaic, lacustrine and alluvial fan deposits and were not isolated to specific depositional environments.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.

Characterization of a Coal Fly Ash-Cement Slurry by the Absolute Foam Index

This dissertation established a standard foam index: the absolute foam index test. This test characterized a wide range of coal fly ash by the absolute volume of air-entraining admixture (AEA) necessary to produce a 15-second metastable foam in a coal fly ash-cement slurry in a specified time. The absolute foam index test was used to characterize fly ash samples having loss on ignition (LOI) values that ranged from 0.17 to 23.3 %wt. The absolute foam index characterized the fly ash samples by absolute volume of AEA, defined as the amount of undiluted AEA solution added to obtain a 15-minute endpoint signified by 15-second metastable foam. Results were compared from several foam index test time trials that used different initial test concentrations to reach termination at selected times. Based on the coefficient of variation (CV), a 15-minute endpoint, with limits of 12 to 18 minutes was chosen. Various initial test concentrations were used to accomplish consistent contact times and concentration gradients for the 15-minute test endpoint for the fly ash samples. A set of four standard concentrations for the absolute foam index test were defined by regression analyses and a procedure simplifying the test process. The set of standard concentrations for the absolute foam index test was determined by analyzing experimental results of 80 tests on coal fly ashes with loss on ignition (LOI) values ranging from 0.39 to 23.3 wt.%. A regression analysis informed selection of four concentrations (2, 6, 10, and 15 vol.% AEA) that are expected to accommodate fly ashes with 0.39 to 23.3 wt.% LOI, depending on the AEA type. Higher concentrations should be used for high-LOI fly ash when necessary. A procedure developed using these standard concentrations is expected to require only 1-3 trials to meet specified endpoint criteria for most fly ashes. The AEA solution concentration that achieved the metastable foam in the foam index test was compared to the AEA equilibrium concentration obtained from the direct adsorption isotherm test with the same fly ash. The results showed that the AEA concentration that satisfied the absolute foam index test was much less than the equilibrium concentration. This indicated that the absolute foam index test was not at or near equilibrium. Rather, it was a dynamic test where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms. Equilibrium isotherm equations obtained from direct isotherm tests were used to calculate the equilibrium concentrations and capacities of fly ash from 0.17 to 10.5% LOI. The results showed that the calculated fly ash capacity was much less than capacities obtained from isotherm tests that were conducted with higher initial concentrations. This indicated that the absolute foam index was not equilibrium. Rather, the test is dynamic where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms for fly ash of 0.17 to 10.5% LOI. Several batches of mortars were mixed for the same fly ash type increasing only the AEA concentration (dosage) in each subsequent batch. Mortar air test results for each batch showed for each increase in AEA concentration, air contents increased until a point where the next increase in AEA concentration resulted in no increase in air content. This was maximum air content that could be achieved by the particular mortar system; the system reached its air capacity at the saturation limit. This concentration of AEA was compared to the critical micelle concentration (CMC) for the AEA and the absolute foam index.

A Benzaldehyde Derivative as a Chelating Ligand: Helical Manganese(II) Coordination Polymers Assembling into a Porous Solid

A molecular, porous crystalline material constructed from neutral helical coordination polymers incorporating manganese(II) ions and two types of bridging ligands, namely the deprotonated form of 2-hydroxy-5-methoxy-3-nitrobenzaldehyde (HL) and isobutyrate (iB−), has been obtained and structurally characterized. Structural analysis reveals that within the coordination polymer each benzaldehyde derivative ligates two manganese ions in 6-membered chelating rings, and the isobutyrate ligands cooperatively chelate either two or three manganese ions. The solid state assembly of the resulting polymeric chains of formula [Mn4(L)2(iB)6]n (1), described in the polar space group R3c, is associated with tubular channels occupied by MeCN solvent molecules (1·xMeCN; x ≤ 9). TGA profiles and PXRD measurements demonstrate that the crystallinity of the solid remains intact in its fully desolvated form, and its stability and crystallinity are ensured up to a temperature of 190 °C. Gas adsorption properties of desolvated crystals were probed, but no remarkable sorption capacity of N2 and only a limited one for CO2 could be observed. Magnetic susceptibility data reveal an antiferromagnetic type of coupling between adjacent manganese(II) ions along the helical chains with energy parameters J1 = −5.9(6) cm−1 and J2 = −1.8(9) cm−1.

Pelagic molybdenum concentration anomalies and the impact of sediment resuspension on the molybdenum budget in two tidal systems of the North Sea

The seasonal dynamics of molybdenum (Mo) were studied in the water column of two tidal basins of the German Wadden Sea (Sylt-Rømø and Spiekeroog) between 2007 and 2011. In contrast to its conservative behaviour in the open ocean, both, losses of more than 50% of the usual concentration level of Mo in seawater and enrichments up to 20% were observed repeatedly in the water column of the study areas. During early summer, Mo removal by adsorption on algae-derived organic matter (e.g. after Phaeocystis blooms) is postulated to be a possible mechanism. Mo bound to organic aggregates is likely transferred to the surface sediment where microbial decomposition enriches Mo in the pore water. First δ98/95Mo data of the study area disclose residual Mo in the open water column being isotopically heavier than MOMo (Mean Ocean Molybdenum) during a negative Mo concentration anomaly, whereas suspended particulate matter shows distinctly lighter values. Based on field observations a Mo isotope enrichment factor of ε = −0.3‰ has been determined which was used to argue against sorption on metal oxide surfaces. It is suggested here that isotope fractionation is caused by biological activity and association to organic matter. Pelagic Mo concentration anomalies exceeding the theoretical salinity-based concentration level, on the other hand, cannot be explained by replenishment via North Sea waters alone and require a supply of excess Mo. Laboratory experiments with natural anoxic tidal flat sediments and modelled sediment displacement during storm events suggest fast and effective Mo release during the resuspension of anoxic sediments in oxic seawater as an important process for a recycling of sedimentary sulphide bound Mo into the water column.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Transport of volatile chlorinated hydrocarbons in unsaturated aggregated media

Transport of volatile hydrocarbons in soils is largely controlled by interactions of vapours with the liquid and solid phase. Sorption on solids of gaseous or dissolved comPounds may be important. Since the contact time between a chemical and a specific sorption site can be rather short, kinetic or mass-transfer resistance effects may be relevant. An existing mathematical model describing advection and diffusion in the gas phase and diffusional transport from the gaseous phase into an intra-aggregate water phase is modified to include linear kinetic sorption on ps-solid and water-solid interfaces. The model accounts for kinetic mass transfer between all three phases in a soil. The solution of the Laplace-transformed equations is inverted numerically. We performed transient column experiments with 1,1,2-Trichloroethane, Trichloroethylene, and Tetrachloroethylene using air-dry solid and water-saturated porous glass beads. The breakthrough curves were calculated based on independently estimated parameters. The model calculations agree well with experimental data. The different transport behaviour of the three compounds in our system primarily depends on Henry's constants.

Anisotropic diffusion at the field scale in a 4-year multi-tracer diffusion and retention experiment – I: Insights from the experimental data

Abstract Claystones are considered worldwide as barrier materials for nuclear waste repositories. In the Mont Terri underground research laboratory (URL), a nearly 4-year diffusion and retention (DR) experiment has been performed in Opalinus Clay. It aimed at (1) obtaining data at larger space and time scales than in laboratory experiments and (2) under relevant in situ conditions with respect to pore water chemistry and mechanical stress, (3) quantifying the anisotropy of in situ diffusion, and (4) exploring possible effects of a borehole-disturbed zone. The experiment included two tracer injection intervals in a borehole perpendicular to bedding, through which traced artificial pore water (APW) was circulated, and a pressure monitoring interval. The APW was spiked with neutral tracers (HTO, HDO, H2O-18), anions (Br, I, SeO4), and cations (Na-22, Ba-133, Sr-85, Cs-137, Co-60, Eu-152, stable Cs, and stable Eu). Most tracers were added at the beginning, some were added at a later stage. The hydraulic pressure in the injection intervals was adjusted according to the measured value in the pressure monitoring interval to ensure transport by diffusion only. Concentration time-series in the APW within the borehole intervals were obtained, as well as 2D concentration distributions in the rock at the end of the experiment after overcoring and subsampling which resulted in �250 samples and �1300 analyses. As expected, HTO diffused the furthest into the rock, followed by the anions (Br, I, SeO4) and by the cationic sorbing tracers (Na-22, Ba-133, Cs, Cs-137, Co-60, Eu-152). The diffusion of SeO4 was slower than that of Br or I, approximately proportional to the ratio of their diffusion coefficients in water. Ba-133 diffused only into �0.1 m during the �4 a. Stable Cs, added at a higher concentration than Cs-137, diffused further into the rock than Cs-137, consistent with a non-linear sorption behavior. The rock properties (e.g., water contents) were rather homogeneous at the centimeter scale, with no evidence of a borehole-disturbed zone. In situ anisotropy ratios for diffusion, derived for the first time directly from field data, are larger for HTO and Na-22 (�5) than for anions (�3�4 for Br and I). The lower ionic strength of the pore water at this location (�0.22 M) as compared to locations of earlier experiments in the Mont Terri URL (�0.39 M) had no notable effect on the anion accessible pore fraction for Cl, Br, and I: the value of 0.55 is within the range of earlier data. Detailed transport simulations involving different codes will be presented in a companion paper.

In-situ diffusion of HTO, 22Na+, Cs+ and I- in Opalinus Clay at the Mont Terri underground rock laboratory

The diffusion properties of the Opalinus Clay were studied in the underground research laboratory at Mont Terri (Canton Jura, Switzerland) and the results were compared with diffusion data measured in the laboratory on small-scale samples. The diffusion of HTO, Na-22(+), Cs+ and I- were investigated for a period of 10 months. The diffusion equipment used in the field experiment was designed in such a way that a solution of tracers was circulated through a sintered metal screen placed at the end of a borehole drilled in the formation. The concentration decrease caused by the diffusion of tracers into the rock could be followed with time and allowed first estimations of the effective diffusion coefficient. After 10 months, the diffusion zone was over-cored and the tracer profiles measured. From these profiles, effective diffusion coefficients and rock capacity factors Could be extracted by applying a two-dimensional transport model including diffusion and sorption. The simulations were done with the reactive transport code CRUNCH. In addition, results obtained from through-diffusion experiments oil small-sized samples with HTO, Cl-36(-) and Na-22(+) are presented and compared with the in situ data. In all cases. excellent agreement between the two data sets exists. Results for Cs+ indicated five times higher diffusion rates relative to HTO. Corresponding laboratory diffusion measurements are still lacking. However. our Cs+ data are in qualitative agreement wish through-diffusion data for Callovo-Oxfordian argillite rock samples. which also indicate significantly higher effective diffusivities for Cs+ relative to HTO.

Diffusion of HTO, Br-, I-, Cs+, Sr-85(2+) and Co-60(2+) in a clay formation: Results and modelling from an in situ experiment in Opalinus Clay

The migration of radioactive and chemical contaminants in clay materials and argillaceous host rocks is characterised by diffusion and retention processes. Valuable information on such processes can be gained by combining diffusion studies at laboratory scale with field migration tests. In this work, the outcome of a multi-tracer in situ migration test performed in the Opalinus Clay formation in the Mont Terri underground rock laboratory (Switzerland) is presented. Thus, 1.16 x 10(5) Bq/L of HTO, 3.96 x 10(3) Bq/L of Sr-85, 6.29 x 10(2) Bq/L of Co-60, 2.01 x 10(-3) mol/L Cs, 9.10 x 10(-4) mol/L I and 1.04 x 10(-3) mol/L Br were injected into the borehole. The decrease of the radioisotope concentrations in the borehole was monitored using in situ gamma-spectrometry. The other tracers were analyzed with state-of-the-art laboratory procedures after sampling of small water aliquots from the reservoir. The diffusion experiment was carried out over a period of one year after which the interval section was overcored and analyzed. Based on the experimental data from the tracer evolution in the borehole and the tracer profiles in the rock, the diffusion of tracers was modelled with the numerical code CRUNCH. The results obtained for HTO (H-3), I- and Br- confirm previous lab and in situ diffusion data. Anionic fluxes into the formation were smaller compared to HTO because of anion exclusion effects. The migration of the cations Sr-85(2+), Cs+ and Co-60(2+) was found to be governed by both diffusion and sorption processes. For Sr-85(2+), the slightly higher diffusivity relative to HTO and the low sorption value are consistent with laboratory diffusion measurements on small-scale samples. In the case of Cs+, the numerically deduced high diffusivity and the Freundlich-type sorption behaviour is also supported by ongoing laboratory data. For Co, no laboratory diffusion data were yet available for comparison; however, the modelled data suggests that Co-60(2+) sorption was weaker than would be expected from available batch sorption data. Overall, the results demonstrate the feasibility of the experimental setup for obtaining high-quality diffusion data for conservative and sorbing tracers. (C) 2007 Elsevier Ltd. All rights reserved.

Long-term diffusion experiment at Mont Terri: first results from field and laboratory data

The diffusion of radionuclides is an important safety aspect for nuclear waste disposal in argillaceous host rocks. A long-term diffusion experiment, termed DI-A, is being carried out at the Mont Terri Rock Laboratory in the Opalinus Clay formation. The aim of this experiment is the understanding of the migration and sorption behaviour of cationic and anionic species in consolidated clays. This study reports on the experimental layout and the first results obtained from the DI-A experiment, which include the investigation of HTO, Na-22(+), Cs+, and I- migration during a period of 1 year by analysing these tracers in the water circulating in the borehole. In addition, results obtained from through-diffusion experiments on small-sized samples with HTO, I-, and Cl-36(-) are presented. The decrease of tracer concentrations in the borehole is fastest for Cs+, followed by Na-22(+), HTO, and finally I-. The chemical composition of the artificial pore water in the borehole shows very little variation with time, thus indicating almost no chemical disturbance around the borehole. Through-diffusion experiments in the laboratory that were performed parallel to the bedding plane with two different methods yielded effective diffusion coefficients for HTO of 4-5 X 10(-11) m(2) s(-1) and significantly lower ones for anions Cl- and I- (0.7-1.6 X 10(-11) m(2) s(-1)). The results indicate the importance of anion exclusion effects arising from the negatively charged clay surfaces. Furthermore, they demonstrate the anisotropic diffusion properties of the clay formation with significantly increased diffusion rates parallel to bedding relative to the perpendicular direction. The tracer data of the in situ experiment were successfully described with 2D diffusion models using diffusion and sorption parameters obtained from the above mentioned and other laboratory studies. The modelling results indicate that HTO and I- diffused with no retardation. The retardation of Na+ and Cs+ could be described by empirical sorption expressions from previously derived batch sorption (Cs+) or diffusion (Na+) experiments. Overall, the obtained results demonstrate the feasibility of the technical concept to study the diffusion of nonsorbing and sorbing tracers in consolidated clays. (C) 2004 Elsevier B.V. All rights reserved.

Ar-40/Ar-39 Evidence for Middle Proterozoic (1300-1500 Ma) Slow Cooling of the Southern Black Hills, South Dakota, Midcontinent, North America: Implications for Early Proterozoic P-T Evolution and Posttectonic Magmatism

Ar-40/Ar-39 total gas and plateau dates from muscovite and biotite in the southern Black Hills, South Dakota, provide evidence for a period of Middle Proterozoic slow cooling. Early Proterozoic (1600-1650 Ma) mica dates were obtained from metasedimentary rocks located in a synformal structure between the Harney Peak and Bear Mountain domes and also south of Bear Mountain. Metamorphic rocks from the dome areas and undeformed samples of the similar to 1710 Ma Harney Peak Granite (HPG) yield Middle Proterozoic mica dates (similar to 1270-1500 Ma). Two samples collected between the synform and Bear Mountain dome yield intermediate total gas mica dates of similar to 1550 Ma. We suggest two end-member interpretations to explain the map pattern of cooling ages: (1) subhorizontal slow cooling of an area which exhibits variation in mica Ar retention intervals or (2) mild folding of a Middle Proterozoic (similar to 1500 Ma) similar to 300 degrees C isotherm. According to the second interpretation, the preservation of older dates between the domes may reflect reactivation of a preexisting synformal structure (and downwarping of relatively cold rocks) during a period of approximately east-west contraction and slow uplift during the Middle Proterozoic. The mica data, together with hornblende data from the Black Hills published elsewhere, indicate that the ambient country-rock temperature at the 3-4 kbar depth of emplacement of the HPG was between 350 degrees C and 500 degrees C, suggesting that the average upper crustal geothermal gradient was 25 degrees-40 degrees C/km prior to intrusion. The thermochronologic data suggest HPG emplacement was followed by a similar to 200 m.y. period of stability and tectonic quiescence with little uplift. We propose that crust thickened during the Early Proterozoic was uplifted and erosionally(?) thinned prior to similar to 1710 Ma and that the HPG magma was emplaced into isostatically stable crust of relatively normal thickness. We speculate that uplift and crustal thinning prior to HPG intrusion was the result of differential thinning of the subcrustal lithosphere beneath the Black Hills. If so, this process would have also caused an increase in mantle heat flux across the Moho and triggered vapor-absent melting of biotite to produce the HPG magma. This scenario for posttectonic granite generation is supported, in part, by the fact that in the whole of the Black Hills, the HPG is spatially associated with the deepest exposed Early Proterozoic country rock.

Polychlorinated Biphenyls in Glaciers. 2. Model Results of Deposition and Incorporation Processes

In previous work, Alpine glaciers have been identified as a secondary source of persistent organic pollutants (POPs). However, detailed understanding of the processes organic chemicals undergo in a glacial system was missing. Here, we present results from a chemical fate model describing deposition and incorporation of polychlorinated biphenyls (PCBs) into an Alpine glacier (Fiescherhorn, Switzerland) and an Arctic glacier (Lomonosovfonna, Norway). To understand PCB fate and dynamics, we investigate the interaction of deposition, sorption to ice and particles in the atmosphere and within the glacier, revolatilization, diffusion and degradation, and discuss the effects of these processes on the fate of individual PCB congeners. The model is able to reproduce measured absolute concentrations in the two glaciers for most PCB congeners. While the model generally predicts concentration profiles peaking in the 1970s, in the measurements, this behavior can only be seen for higher-chlorinated PCB congeners on Fiescherhorn glacier. We suspect seasonal melt processes are disturbing the concentration profiles of the lower-chlorinated PCB congeners. While a lower-chlorinated PCB congener is mainly deposited by dry deposition and almost completely revolatilized after deposition, a higher-chlorinated PCB congener is predominantly transferred to the glacier surface by wet deposition and then is incorporated into the glacier ice. The incorporated amounts of PCBs are higher on the Alpine glacier than on the Arctic glacier due to the higher precipitation rate and aerosol particle concentration on the former. Future studies should include the effects of seasonal melt processes, calculate the quantities of PCBs incorporated into the entire glacier surface, and estimate the quantity of chemicals released from glaciers to determine the importance of glaciers as a secondary source of organic chemicals to remote aquatic ecosystems.