Identifiability of diffusion and sorption parameters from in situ diffusion experiments by using simultaneously tracer dilution and claystone data

In situ diffusion experiments are performed in geological formations at underground research laboratories to overcome the limitations of laboratory diffusion experiments and investigate scale effects. Tracer concentrations are monitored at the injection interval during the experiment (dilution data) and measured from host rock samples around the injection interval at the end of the experiment (overcoring data). Diffusion and sorption parameters are derived from the inverse numerical modeling of the measured tracer data. The identifiability and the uncertainties of tritium and Na-22(+) diffusion and sorption parameters are studied here by synthetic experiments having the same characteristics as the in situ diffusion and retention (DR) experiment performed on Opalinus Clay. Contrary to previous identifiability analyses of in situ diffusion experiments, which used either dilution or overcoring data at approximate locations, our analysis of the parameter identifiability relies simultaneously on dilution and overcoring data, accounts for the actual position of the overcoring samples in the claystone, uses realistic values of the standard deviation of the measurement errors, relies on model identification criteria to select the most appropriate hypothesis about the existence of a borehole disturbed zone and addresses the effect of errors in the location of the sampling profiles. The simultaneous use of dilution and overcoring data provides accurate parameter estimates in the presence of measurement errors, allows the identification of the right hypothesis about the borehole disturbed zone and diminishes other model uncertainties such as those caused by errors in the volume of the circulation system and the effective diffusion coefficient of the filter. The proper interpretation of the experiment requires the right hypothesis about the borehole disturbed zone. A wrong assumption leads to large estimation errors. The use of model identification criteria helps in the selection of the best model. Small errors in the depth of the overcoring samples lead to large parameter estimation errors. Therefore, attention should be paid to minimize the errors in positioning the depth of the samples. The results of the identifiability analysis do not depend on the particular realization of random numbers. (C) 2012 Elsevier B.V. All rights reserved.

An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

Tracer profiles across argillaceous formations: A tool to constrain transport processes.

The spatial distributions of non-reactive natural tracers (anions, stable water isotopes, noble gases) in pore water of clay-rich formations were studied at nine sites. Regular curved profiles were identified in most cases. Transport modeling considering diffusion, advection and available constraints on the paleo-hydrogeological evolution indicates generally that diffusion alone can explain the observations, whereas a marked advective component would distort the profiles and so is not consistent with the data.

Tracer applications of noble gas radionuclides in the geosciences

Noble gas radionuclides, including 81Kr (t1/2 = 229,000 years), 85Kr (t1/2 = 10.8 years), and 39Ar (t1/2 = 269 years), possess nearly ideal chemical and physical properties for studies of earth and environmental processes. Recent advances in Atom Trap Trace Analysis (ATTA), a laser-based atom counting method, have enabled routine measurements of the radiokrypton isotopes, as well as the demonstration of the ability to measure 39Ar in environmental samples. Here we provide an overview of the ATTA technique, and a survey of recent progress made in several laboratories worldwide.We review the application of noble gas radionuclides in the geosciences and discuss how ATTA can help advance these fields, specifically: determination of groundwater residence times using 81Kr, 85Kr, and 39Ar; dating old glacial ice using 81Kr; and an 39Ar survey of the main water masses of the oceans, to study circulation pathways and estimate mean residence times. Other scientific questions involving a deeper circulation of fluids in the Earth's crust and mantle are also within the scope of future applications. We conclude that the geoscience community would greatly benefit from an ATTA facility dedicated to this field, with instrumentation for routine measurements, as well as for research on further development of ATTA methods.

A multi-tracer study of groundwater origin and transit-time in the aquifers of the Venice region (Italy)

Located in the northeastern region of Italy, the Venetian Plain (VP) is a sedimentary basin containing an extensively exploited groundwater system. The northern part is characterised by a large undifferentiated phreatic aquifer constituted by coarse grain alluvial deposits and recharged by local rainfalls and discharges from the rivers Brenta and Piave. The southern plain is characterised by a series of aquitards and sandy aquifers forming a well-defined artesian multi-aquifer system. In order to determine origins, transit times and mixing proportions of different components in groundwater (GW), a multi tracer study (H, He/He, C, CFC, SF, Kr, Ar, Sr/Sr, O, H, cations, and anions) has been carried out in VP between the rivers Brenta and Piave. The geochemical pattern of GW allows a distinction of the different water origins in the system, in particular based on View the MathML source HCO3-,SO42-,Ca/Mg,NO3-, O, H. A radiogenic Sr signature clearly marks GW originated from the Brenta and Tertiary catchments. End-member analysis and geochemical modelling highlight the existence of a mixing process involving waters recharged from the Brenta and Piave rivers, from the phreatic aquifer and from another GW reservoirs characterised by very low mineralization. Noble gas excesses in respect to atmospheric equilibrium occur in all samples, particularly in the deeper aquifers of the Piave river, but also in phreatic water of the undifferentiated aquifers. He–H ages in the phreatic aquifer and in the shallower level of the multi-aquifer system indicate recharge times in the years 1970–2008. The progression of H–He ages with the distance from the recharge areas together with initial tritium concentration (H + Hetrit) imply an infiltration rate of about 1 km/y and the absence of older components in these GW. SF and Kr data corroborate these conclusions. H − He ages in the deeper artesian aquifers suggest a dilution process with older, tritium free waters. C Fontes–Garnier model ages of the old GW components range from 1 to 12 ka, yielding an apparent GW velocity of about 1–10 m/y. Increase of radiogenic He follows the progression of C ages. Ar, radiogenic He and C tracers yield model-dependent age-ranges in overall good agreement once diffusion of C from aquitards, GW dispersion, lithogenic Ar production, and He production-rate heterogeneities are taken into account. The rate of radiogenic He increase with time, deduced by comparison with C model ages, is however very low compared to other studies. Comparison with C and C data obtained 40 years ago on the same aquifer system shows that exploitation of GW caused a significant loss of the old groundwater reservoir during this time.

Development of a specific tracer for metabolic imaging of alveolar echinococcosis: A preclinical study

Positron emission tomography (PET)-computed tomography (CT) using [18F]-fluorodeoxyglucose (FDG) (FDG-PET/CT) is a valuable method for initial staging and follow up of patients with alveolar echinococcosis (AE). However, the cells responsible for FDG uptake have not been clearly identified. The main goal of our study was to evaluate the uptake of PET tracers by the cells involved in the host-parasite reaction around AE lesions as the first step to develop a specific PET tracer that would allow direct assessment of parasite viability in AE. Candidate molecules ([18F]-fluorotyrosine (FET), [18F]-fluorothymidine (FLT), and [18F]-fluorometylcholine (FMC), were compared to FDG by in vitro studies on human leukocytes and parasite vesicles. Our results confirmed that FDG was mainly consumed by immune cells and showed that FLT was the best candidate tracer for parasite metabolism. Indeed, parasite cells exhibited high uptake of FLT. We also performed PET/CT scans in mice infected intraperitoneally with E. multilocularis metacestodes. PET images showed no FDG or FLT uptake in parasitic lesions. This preliminary study assessed the metabolic activity of human leukocytes and AE cells using radiolabeling. Future studies could develop a specific PET tracer for AE lesions to improve lesion detection and echinococcosis treatment in patients. Our results demonstrated that a new animal model is needed for preclinical PET imaging to better mimic human hepatic and/or periparasitic metabolism.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.