Modelling of the migration of trace elements along groundwater flowpaths by using a steady state approach application to the site at El Berrocal (Spain)

The applications of geochemical modelling to natural water systems mostly rely upon the equilibrium assumption. This is in principle justified for deep groundwater systems with long water residence times.In addition, trace element behaviour in these natural systems has been normally modelled by taking into account the precipitation and dissolution of individual trace element phases.Recent geochemical modelling work related to natural analogue systems would indicate that: (1) the equilibrium condition is restricted to a limited number of components in long water residence times and (2) the behaviour of trace metals is very much connected to the major component cycling in the system.There is a need to develop our geochemical modelling capabilities to take these two facts into account in order to be able to predict the behaviour of trace components (radionuclides) in a geological repository.In this work we report on the successful application of steady-state kinetics in conjunction with co-dissolution/co-precipitation approaches to model trace element geochemistry in the natural analogue system at El Berrocal. The evolution of major components of the system (Ca (II), Al (III), CO₃²⁻ and Si) as well the trace elements investigated (U, Ba and Mn) is well reproduced by using the coupling between steady-state kinetics and co-dissolution/co-precipitation approaches.     

Uranium transport around the reactor zone at Bangombé and Okélobondo (Oklo): examples of hydrogeological and geochemical model integration and data evaluation

The sites at Bangombé and Okélobondo (Oklo) in Gabon provide a unique opportunity to study the behaviour of products from natural nuclear reactions in the vicinity of reactor zones which were active around two billion years ago. The Commission of the European Communities initiated the Oklo Natural Analogue Programme. One of the principal aims was to study indications of present time migration of elements from the reactor zones under ambient conditions. The hydrogeological and hydrochemical data from the Oklo sites were modelled in order to better understand the geochemical behaviour of radionuclides in the natural system, by using independent models and by comparing the modelling outcome. Two modelling approaches were used: M3 code (hydrochemical mixing and mass balance model), developed by the Swedish Nuclear Fuel and Waste Management Company (SKB) and HYTEC (reactive transport model) developed by Ecole des Mines de Paris.Two different reactor zones were studied: Bangombé, a shallow site, the reactor being at 11 m depth, and OK84 at Okélobondo, situated at about 450 m depth, more comparable with a real repository location. This allowed the validation of modelling tools in two different sedimentary environments: one shallow, with a more homogeneous layering situated in an area of meteoric alteration, and the other offering the opportunity to study radionuclide migration from the reaction zone over a distance of 450 m through very heterogeneous sedimentary layers.The modeling results indicate that the chemical reactions retarding radionuclide transport are very different at the two sites. At Bangombé, the decomposition of organic material consumes oxygen and at Okélobondo the oxygen is consumed by inorganic reactions resulting, in both cases, in uranium retardation. Both modelling approaches (statistic with M3 code and deterministic with HYTEC code) could describe this situation.The goal of this exercise is to test codes which can help to describe and understand the processes taking place at the sites, validate the models with in situ data, and thus build confidence in the tools used for future site characterization. Ultimately, this allows identifying and selecting processes and parameters that can be used as input into repository performance assessment calculations and modelling exercises.     

Conceptual and Numerical Modeling of Radionuclide Transport and Retention in Near-Surface Systems

Scenarios of barrier failure and radionuclide release to the near-surface environment are important to consider within performance and safety assessments of repositories for nuclear waste. A geological repository for spent nuclear fuel is planned at Forsmark, Sweden. Conceptual and numerical reactive transport models were developed in order to assess the retention capacity of the Quaternary till and clay deposits for selected radionuclides, in the event of an activity release from the repository. The elements considered were carbon (C), chlorine (Cl), cesium (Cs), iodine (I), molybdenum (Mo), niobium (Nb), nickel (Ni), radium (Ra), selenium (Se), strontium (Sr), technetium (Tc), thorium (Th), and uranium (U). According to the numerical predictions, the repository-derived nuclides that would be most significantly retained are Th, Ni, and Cs, mainly through sorption onto clays, followed by U, C, Sr, and Ra, trapped by sorption and/or incorporation into mineral phases.     

Uranium transport around the reactor zone at Bangombé and Okélobondo (Oklo): examples of hydrogeological and geochemical model integration and data evaluation

The sites at Bangombé and Okélobondo (Oklo) in Gabon provide a unique opportunity to study the behaviour of products from natural nuclear reactions in the vicinity of reactor zones which were active around two billion years ago. The Commission of the European Communities initiated the Oklo Natural Analogue Programme. One of the principal aims was to study indications of present time migration of elements from the reactor zones under ambient conditions. The hydrogeological and hydrochemical data from the Oklo sites were modelled in order to better understand the geochemical behaviour of radionuclides in the natural system, by using independent models and by comparing the modelling outcome. Two modelling approaches were used: M3 code (hydrochemical mixing and mass balance model), developed by the Swedish Nuclear Fuel and Waste Management Company (SKB) and HYTEC (reactive transport model) developed by Ecole des Mines de Paris. Two different reactor zones were studied: Bangombé, a shallow site, the reactor being at 11 m depth, and OK84 at Okélobondo, situated at about 450 m depth, more comparable with a real repository location. This allowed the validation of modelling tools in two different sedimentary environments: one shallow, with a more homogeneous layering situated in an area of meteoric alteration, and the other offering the opportunity to study radionuclide migration from the reaction zone over a distance of 450 m through very heterogeneous sedimentary layers. The modeling results indicate that the chemical reactions retarding radionuclide transport are very different at the two sites. At Bangombé, the decomposition of organic material consumes oxygen and at Okélobondo the oxygen is consumed by inorganic reactions resulting, in both cases, in uranium retardation. Both modelling approaches (statistic with M3 code and deterministic with HYTEC code) could describe this situation. The goal of this exercise is to test codes which can help to describe and understand the processes taking place at the sites, validate the models with in situ data, and thus build confidence in the tools used for future site characterization. Ultimately, this allows identifying and selecting processes and parameters that can be used as input into repository performance assessment calculations and modelling exercises.     

Present state and future directions of modeling of geochemistry in hydrogeological systems

A first step towards understanding and controlling the fate and dissemination of radioactive waste is to create a concise and comprehensive theoretical framework for the rather non-linear processes involved--hence, the need for geochemical models. Two classes of geochemical models are commonly used, i.e., static and hydrodynamic models. In contrast to static models, hydrodynamic models combine geochemical reactions with hydrogeological processes such as ground-water flow, diffusion and dispersion. In this review, we examine the present state of geochemical models in terms of included processes, thermodynamic databases, missing phenomena, numerical behavior and performance. It is shown that over the past decade, significant progress has been made with respect to modeling of geochemistry in hydrodynamic systems: this is illustrated by describing several applications. Finally, we focus on the perspectives of geochemical modeling in the assessment of the safety of nuclear waste disposal.     

Quantifying uranium complexation by groundwater dissolved organic carbon using asymmetrical flow field-flow fractionation

The long-term mobility of actinides in groundwaters is important for siting nuclear waste facilities and managing waste-rock piles at uranium mines. Dissolved organic carbon (DOC) may influence the mobility of uranium, but few field-based studies have been undertaken to examine this in typical groundwaters. In addition, few techniques are available to isolate DOC and directly quantify the metals complexed to it. Determination of U-organic matter association constants from analysis of field-collected samples compliments laboratory measurements, and these constants are needed for accurate transport calculations. The partitioning of U to DOC in a clay-rich aquitard was investigated in 10 groundwater samples collected between 2 and 30 m depths at one test site. A positive correlation was observed between the DOC (4-132 mg/L) and U concentrations (20-603 microg/L). The association of U and DOC was examined directly using on-line coupling of Asymmetrical Flow Field-Flow Fractionation (AsFlFFF) with UV absorbance (UVA) and inductively coupled plasma-mass spectrometer (ICP-MS) detectors. This method has the advantages of utilizing very small sample volumes (20-50 microL) as well as giving molecular weight information on U-organic matter complexes. AsFlFFF-UVA results showed that 47-98% of the DOC (4-136 mg C/L) was recovered in the AsFlFFF analysis, of which 25-64% occurred in the resolvable peak. This peak corresponded to a weight-average molecular weight of about 900-1400 Daltons (Da). In all cases, AsFlFFF-ICP-MS suggested that相似文献   

Humans and Ecosystems Over the Coming Millennia: Overview of a Biosphere Assessment of Radioactive Waste Disposal in Sweden

This is an overview of the strategy used to describe the effects of a potential release from a radioactive waste repository on human exposure and future environments. It introduces a special issue of AMBIO, in which 13 articles show ways of understanding and characterizing the future. The study relies mainly on research performed in the context of a recent safety report concerning a repository for spent nuclear fuel in Sweden (the so-called SR-Site project). The development of a good understanding of on-site processes and acquisition of site-specific data facilitated the development of new approaches for assessment of surface ecosystems. A systematic and scientifically coherent methodology utilizes the understanding of the current spatial and temporal dynamics as an analog for future conditions. We conclude that future ecosystem can be inferred from a few variables and that this multidisciplinary approach is relevant in a much wider context than radioactive waste.     

Ecohydrological Responses to Diversion of Groundwater: Case Study of a Deep-Rock Repository for Spent Nuclear Fuel in Sweden

Planning and license applications concerning groundwater diversion in areas containing water-dependent or water-favored habitats must take into account both hydrological effects and associated ecological consequences. There is at present no established methodology to assess such ecohydrological responses. Thus, this paper describes a new stepwise methodology to assess ecohydrological responses to groundwater diversion from, e.g., water-drained pits, shafts, tunnels, and caverns in rock below the groundwater table. The methodology is illustrated using the planned deep-rock repository for spent nuclear fuel at Forsmark in central Sweden as a case study, offering access to a unique hydrological and ecological dataset. The case study demonstrates that results of ecohydrological assessments can provide useful inputs to planning of monitoring programs and mitigation measures in infrastructure projects. As a result of the assessment, artificial water supply to wetlands is planned in order to preserve biological diversity, nature values, and vulnerable species.     

Comparison of solubility databases used for HLW performance assessment

Elemental solubility limits are employed in many performance assessment studies of the disposal of high-level radioactive waste. A comparison of solubility databases used in a number of recent assessments from Canada, Finland, Japan, Sweden, Switzerland and the U.S.A. is presented. Despite similarities in the chemical conditions under which such solubilities are defined, there are very large ranges in the values selected for most elements (several orders of magnitude). In some cases, differences in solubilities can be attributed to differences in the thermodynamic databases used but, due to generally incomplete documentation, both of such databases and the criteria for data selection, detailed analysis of the reason for such discrepancies is precluded. A key factor, however, is the fairly subjective choice of the solid phase considered to be “solubility-limiting”. Procedures are recommended which would increase the transparency of solubility database compilation in future analyses.     

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.     

Modelling geochemical and microbial consumption of dissolved oxygen after backfilling a high level radiactive waste repository

Dissolved oxygen (DO) left in the voids of buffer and backfill materials of a deep geological high level radioactive waste (HLW) repository could cause canister corrosion. Available data from laboratory and in situ experiments indicate that microbes play a substantial role in controlling redox conditions near a HLW repository. This paper presents the application of a coupled hydro-bio-geochemical model to evaluate geochemical and microbial consumption of DO in bentonite porewater after backfilling of a HLW repository designed according to the Swedish reference concept. In addition to geochemical reactions, the model accounts for dissolved organic carbon (DOC) respiration and methane oxidation. Parameters for microbial processes were derived from calibration of the REX in situ experiment carried out at the Asp? underground laboratory. The role of geochemical and microbial processes in consuming DO is evaluated for several scenarios. Numerical results show that both geochemical and microbial processes are relevant for DO consumption. However, the time needed to consume the DO trapped in the bentonite buffer decreases dramatically from several hundreds of years when only geochemical processes are considered to a few weeks when both geochemical reactions and microbially-mediated DOC respiration and methane oxidation are taken into account simultaneously.     

Biodynamic modelling and the prediction of accumulated trace metal concentrations in the polychaete Arenicola marina

The use of biodynamic models to understand metal uptake directly from sediments by deposit-feeding organisms still represents a special challenge. In this study, accumulated concentrations of Cd, Zn and Ag predicted by biodynamic modelling in the lugworm Arenicola marina have been compared to measured concentrations in field populations in several UK estuaries. The biodynamic model predicted accumulated field Cd concentrations remarkably accurately, and predicted bioaccumulated Ag concentrations were in the range of those measured in lugworms collected from the field. For Zn the model showed less but still good comparability, accurately predicting Zn bioaccumulation in A. marina at high sediment concentrations but underestimating accumulated Zn in the worms from sites with low and intermediate levels of Zn sediment contamination. Therefore, it appears that the physiological parameters experimentally derived for A. marina are applicable to the conditions encountered in these environments and that the assumptions made in the model are plausible.     

Factors influencing natural occurrence of fluoride-rich groundwaters: a case study in the southeastern part of the Korean Peninsula

Factors regulating the occurrence of fluoride in groundwater were investigated using natural isotope tracers and geochemical data in the southeastern part of Korea where about 10% of the total public water-supply wells (n=422) inspected in this study had fluoride levels exceeding the drinking water limit of Korea (1.5 mg/l). The F-rich public wells are mostly distributed along the major faults, especially in the terrain of the F-rich granitic rocks. The stable isotope analysis results provide substantial information for the relative ages of groundwaters. It is revealed that the F-rich groundwaters are deeply circulating paleogroundwaters and occur along the faults due to upward flow along the fault plane. While reacting with granitic rocks for a prolonged period, the F concentrations of groundwater are continuously enriched even after the groundwater reaches an equilibrium state with respect to fluorite (CaF2) due to removal of Ca by precipitation of calcite (CaCO3). These observations reflect that rock chemistry, groundwater age, well depth, and geologic structure are the important factors controlling the occurrence of high F groundwaters. However, high F groundwaters are rarely observed in the fault zones where the associated fractures are widely developed. Isotopic signature provides an evidence for deep penetration of recently recharged groundwater into the wide fault zone, indicating that the hydrologic condition of the fault is also an important factor controlling the occurrence of high F groundwaters.     

Analysis of water flow paths: methodology and example calculations for a potential geological repository in Sweden

Safety assessment related to the siting of a geological repository for spent nuclear fuel deep in the bedrock requires identification of potential flow paths and the associated travel times for radionuclides originating at repository depth. Using the Laxemar candidate site in Sweden as a case study, this paper describes modeling methodology, data integration, and the resulting water flow models, focusing on the Quaternary deposits and the upper 150 m of the bedrock. Example simulations identify flow paths to groundwater discharge areas and flow paths in the surface system. The majority of the simulated groundwater flow paths end up in the main surface waters and along the coastline, even though the particles used to trace the flow paths are introduced with a uniform spatial distribution at a relatively shallow depth. The calculated groundwater travel time, determining the time available for decay and retention of radionuclides, is on average longer to the coastal bays than to other biosphere objects at the site. Further, it is demonstrated how GIS-based modeling can be used to limit the number of surface flow paths that need to be characterized for safety assessment. Based on the results, the paper discusses an approach for coupling the present models to a model for groundwater flow in the deep bedrock.     

Modeling the long-term transport and accumulation of radionuclides in the landscape for derivation of dose conversion factors

To evaluate the radiological impact of potential releases to the biosphere from a geological repository for spent nuclear fuel, it is necessary to assess the long-term dynamics of the distribution of radionuclides in the environment. In this paper, we propose an approach for making prognoses of the distribution and fluxes of radionuclides released from the geosphere, in discharges of contaminated groundwater, to an evolving landscape. The biosphere changes during the temperate part (spanning approximately 20,000 years) of an interglacial period are handled by building biosphere models for the projected succession of situations. Radionuclide transport in the landscape is modeled dynamically with a series of interconnected radioecological models of those ecosystem types (sea, lake, running water, mire, agricultural land and forest) that occur at present, and are projected to occur in the future, in a candidate area for a geological repository in Sweden. The transformation between ecosystems is modeled as discrete events occurring every thousand years by substituting one model by another. Examples of predictions of the radionuclide distribution in the landscape are presented for several scenarios with discharge locations varying in time and space. The article also outlines an approach for estimating the exposure of man resulting from all possible reasonable uses of a potentially contaminated landscape, which was used for derivation of Landscape Dose Factors.     

Modeling Emissions from Elevated Cylindrical Bins

Abstract

This paper demonstrates how wind tunnel modeling data that accurately describe plume characteristics near an unconventional emission source can be used to improve the near-field downwind plume profiles predicted by conventional air pollution dispersion models. The study considers a vertical, cylindrical-shaped, elevated bin similar to large product storage bins that can be found at many industrial plant sites. Two dispersion models are considered: the U.S. Environmental Protection Agency's ISC2(ST) model and the Ontario Ministry of the Environment and Energy's GAS model. The wind tunnel study showed that plume behavior was contrary to what was predicted using conventional dispersion models such as ISC2(ST) and GAS and default values of input parameters. The wind tunnel data were used to develop a protocol for correcting the dispersion models inputs, resulting in a substantial improvement in the accuracy of the dispersion estimates.