Modelling of the reactive transport for rock salt-brine in geological repository systems based on improved thermodynamic database (Invited)

NASA Astrophysics Data System (ADS)

Müller, W.; Alkan, H.; Xie, M.; Moog, H.; Sonnenthal, E. L.

2009-12-01

The release and migration of toxic contaminants from the disposed wastes is one of the main issues in long-term safety assessment of geological repositories. In the engineered and geological barriers around the nuclear waste emplacements chemical interactions between the components of the system may affect the isolation properties considerably. As the chemical issues change the transport properties in the near and far field of a nuclear repository, modelling of the transport should also take the chemistry into account. The reactive transport modelling consists of two main components: a code that combines the possible chemical reactions with thermo-hydrogeological processes interactively and a thermodynamic databank supporting the required parameters for the calculation of the chemical reactions. In the last decade many thermo-hydrogeological codes were upgraded to include the modelling of the chemical processes. TOUGHREACT is one of these codes. This is an extension of the well known simulator TOUGH2 for modelling geoprocesses. The code is developed by LBNL (Lawrence Berkeley National Laboratory, Univ. of California) for the simulation of the multi-phase transport of gas and liquid in porous media including heat transfer. After the release of its first version in 1998, this code has been applied and improved many times in conjunction with considerations for nuclear waste emplacement. A recent version has been extended to calculate ion activities in concentrated salt solutions applying the Pitzer model. In TOUGHREACT, the incorporated equation of state module ECO2N is applied as the EOS module for non-isothermal multiphase flow in a fluid system of H2O-NaCl-CO2. The partitioning of H2O and CO2 between liquid and gas phases is modelled as a function of temperature, pressure, and salinity. This module is applicable for waste repositories being expected to generate or having originally CO2 in the fluid system. The enhanced TOUGHREACT uses an EQ3/6-formatted database for both Pitzer ion-interaction parameters and thermodynamic equilibrium constants. The reliability of the parameters is as important as the accuracy of the modelling tool. For this purpose the project THEREDA (www.thereda.de)was set up. The project aims at a comprehensive and internally consistent thermodynamic reference database for geochemical modelling of near and far-field processes occurring in repositories for radioactive wastes in various host rock formations. In the framework of the project all data necessary to perform thermodynamic equilibrium calculations for elevated temperature in the system of oceanic salts are under revision, and it is expected that related data will be available for download by 2010-03. In this paper the geochemical issues that can play an essential role for the transport of radioactive contaminants within and around waste repositories are discussed. Some generic calculations are given to illustrate the geochemical interactions and their probable effects on the transport properties around HLW emplacements and on CO2 generating and/or containing repository systems.

Predictions of hydrothermal alteration within near-ridge oceanic crust from coordinated geochemical and fluid flow models

USGS Publications Warehouse

Wetzel, L.R.; Raffensperger, Jeff P.; Shock, E.L.

2001-01-01

Coordinated geochemical and hydrological calculations guide our understanding of the composition, fluid flow patterns, and thermal structure of near-ridge oceanic crust. The case study presented here illustrates geochemical and thermal changes taking place as oceanic crust ages from 0.2 to 1.0 Myr. Using a finite element code, we model fluid flow and heat transport through the upper few hundred meters of an abyssal hill created at an intermediate spreading rate. We use a reaction path model with a customized database to calculate equilibrium fluid compositions and mineral assemblages of basalt and seawater at 500 bars and temperatures ranging from 150 to 400??C. In one scenario, reaction path calculations suggest that volume increases on the order of 10% may occur within portions of the basaltic basement. If this change in volume occurred, it would be sufficient to fill all primary porosity in some locations, effectively sealing off portions of the oceanic crust. Thermal profiles resulting from fluid flow simulations indicate that volume changes along this possible reaction path occur primarily within the first 0.4 Myr of crustal aging. ?? 2001 Elsevier Science B.V. All rights reserved.

Interface COMSOL-PHREEQC (iCP), an efficient numerical framework for the solution of coupled multiphysics and geochemistry

NASA Astrophysics Data System (ADS)

Nardi, Albert; Idiart, Andrés; Trinchero, Paolo; de Vries, Luis Manuel; Molinero, Jorge

2014-08-01

This paper presents the development, verification and application of an efficient interface, denoted as iCP, which couples two standalone simulation programs: the general purpose Finite Element framework COMSOL Multiphysics® and the geochemical simulator PHREEQC. The main goal of the interface is to maximize the synergies between the aforementioned codes, providing a numerical platform that can efficiently simulate a wide number of multiphysics problems coupled with geochemistry. iCP is written in Java and uses the IPhreeqc C++ dynamic library and the COMSOL Java-API. Given the large computational requirements of the aforementioned coupled models, special emphasis has been placed on numerical robustness and efficiency. To this end, the geochemical reactions are solved in parallel by balancing the computational load over multiple threads. First, a benchmark exercise is used to test the reliability of iCP regarding flow and reactive transport. Then, a large scale thermo-hydro-chemical (THC) problem is solved to show the code capabilities. The results of the verification exercise are successfully compared with those obtained using PHREEQC and the application case demonstrates the scalability of a large scale model, at least up to 32 threads.

Geomechanical/Geochemical Modeling Studies Conducted within theInternational DECOVALEX Project

DOE Office of Scientific and Technical Information (OSTI.GOV)

Birkholzer, J.T.; Rutqvist, J.; Sonnenthal, E.L.

2005-10-19

The DECOVALEX project is an international cooperative project initiated by SKI, the Swedish Nuclear Power Inspectorate, with participation of about 10 international organizations. The general goal of this project is to encourage multidisciplinary interactive and cooperative research on modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in geologic formations in support of the performance assessment for underground storage of radioactive waste. One of the research tasks, initiated in 2004 by the U.S. Department of Energy (DOE), addresses the long-term impact of geomechanical and geochemical processes on the flow conditions near waste emplacement tunnels. Within this task, four international research teams conduct predictive analysismore » of the coupled processes in two generic repositories, using multiple approaches and different computer codes. Below, we give an overview of the research task and report its current status.« less

Geomechanical/ Geochemical Modeling Studies onducted Within the International DECOVALEX Project

DOE Office of Scientific and Technical Information (OSTI.GOV)

J.T. Birkholzer; J. Rutqvist; E.L. Sonnenthal

2006-02-01

The DECOVALEX project is an international cooperative project initiated by SKI, the Swedish Nuclear Power Inspectorate, with participation of about 10 international organizations. The general goal of this project is to encourage multidisciplinary interactive and cooperative research on modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in geologic formations in support of the performance assessment for underground storage of radioactive waste. One of the research tasks, initiated in 2004 by the U.S. Department of Energy (DOE), addresses the long-term impact of geomechanical and geochemical processes on the flow conditions near waste emplacement tunnels. Within this task, four international research teams conduct predictive analysismore » of the coupled processes in two generic repositories, using multiple approaches and different computer codes. Below, we give an overview of the research task and report its current status.« less

PFLOTRAN: Reactive Flow & Transport Code for Use on Laptops to Leadership-Class Supercomputers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hammond, Glenn E.; Lichtner, Peter C.; Lu, Chuan

PFLOTRAN, a next-generation reactive flow and transport code for modeling subsurface processes, has been designed from the ground up to run efficiently on machines ranging from leadership-class supercomputers to laptops. Based on an object-oriented design, the code is easily extensible to incorporate additional processes. It can interface seamlessly with Fortran 9X, C and C++ codes. Domain decomposition parallelism is employed, with the PETSc parallel framework used to manage parallel solvers, data structures and communication. Features of the code include a modular input file, implementation of high-performance I/O using parallel HDF5, ability to perform multiple realization simulations with multiple processors permore » realization in a seamless manner, and multiple modes for multiphase flow and multicomponent geochemical transport. Chemical reactions currently implemented in the code include homogeneous aqueous complexing reactions and heterogeneous mineral precipitation/dissolution, ion exchange, surface complexation and a multirate kinetic sorption model. PFLOTRAN has demonstrated petascale performance using 2{sup 17} processor cores with over 2 billion degrees of freedom. Accomplishments achieved to date include applications to the Hanford 300 Area and modeling CO{sub 2} sequestration in deep geologic formations.« less

Data Collection Handbook to Support Modeling Impacts of Radioactive Material in Soil and Building Structures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yu, Charley; Kamboj, Sunita; Wang, Cheng

2015-09-01

This handbook is an update of the 1993 version of the Data Collection Handbook and the Radionuclide Transfer Factors Report to support modeling the impact of radioactive material in soil. Many new parameters have been added to the RESRAD Family of Codes, and new measurement methodologies are available. A detailed review of available parameter databases was conducted in preparation of this new handbook. This handbook is a companion document to the user manuals when using the RESRAD (onsite) and RESRAD-OFFSITE code. It can also be used for RESRAD-BUILD code because some of the building-related parameters are included in this handbook.more » The RESRAD (onsite) has been developed for implementing U.S. Department of Energy Residual Radioactive Material Guidelines. Hydrogeological, meteorological, geochemical, geometrical (size, area, depth), crops and livestock, human intake, source characteristic, and building characteristic parameters are used in the RESRAD (onsite) code. The RESRAD-OFFSITE code is an extension of the RESRAD (onsite) code and can also model the transport of radionuclides to locations outside the footprint of the primary contamination. This handbook discusses parameter definitions, typical ranges, variations, and measurement methodologies. It also provides references for sources of additional information. Although this handbook was developed primarily to support the application of RESRAD Family of Codes, the discussions and values are valid for use of other pathway analysis models and codes.« less

A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development

USGS Publications Warehouse

Abrams , Robert H.; Loague, Keith

2000-01-01

This paper, the first of two parts [see Abrams and Loague, this issue], takes the compartmentalized approach for the geochemical evolution of redox zones presented by Abrams et al. [1998] and embeds it within a solute transport framework. In this paper the compartmentalized approach is generalized to facilitate the description of its incorporation into a solute transport simulator. An equivalent formulation is developed which removes any discontinuities that may occur when switching compartments. Rate‐limited redox reactions are modeled with a modified Monod relationship that allows either the organic substrate or the electron acceptor to be the rate‐limiting reactant. Thermodynamic constraints are used to inhibit lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium on the lower‐energy reactions. The procedure used allows any redox reaction to be simulated as being kinetically limited or thermodynamically limited, depending on local geochemical conditions. Empirical reaction inhibition methods are not needed. The sequential iteration approach (SIA), a technique which allows the number of solute transport equations to be reduced, is adopted to solve the coupled geochemical/solute transport problem. When the compartmentalized approach is embedded within the SIA, with the total analytical concentration of each component as the dependent variable in the transport equation, it is possible to reduce the number of transport equations even further than with the unmodified SIA. A one‐dimensional, coupled geochemical/solute transport simulation is presented in which redox zones evolve dynamically in time and space. The compartmentalized solute transport (COMPTRAN) model described in this paper enables the development of redox zones to be simulated under both kinetic and thermodynamic constraints. The modular design of COMPTRAN facilitates the use of many different, preexisting solute transport and geochemical codes. The companion paper [Abrams and Loague, this issue] presents examples of the application of COMPTRAN to field‐scale problems.

Reactive transport modeling in the subsurface environment with OGS-IPhreeqc

NASA Astrophysics Data System (ADS)

He, Wenkui; Beyer, Christof; Fleckenstein, Jan; Jang, Eunseon; Kalbacher, Thomas; Naumov, Dimitri; Shao, Haibing; Wang, Wenqing; Kolditz, Olaf

2015-04-01

Worldwide, sustainable water resource management becomes an increasingly challenging task due to the growth of population and extensive applications of fertilizer in agriculture. Moreover, climate change causes further stresses to both water quantity and quality. Reactive transport modeling in the coupled soil-aquifer system is a viable approach to assess the impacts of different land use and groundwater exploitation scenarios on the water resources. However, the application of this approach is usually limited in spatial scale and to simplified geochemical systems due to the huge computational expense involved. Such computational expense is not only caused by solving the high non-linearity of the initial boundary value problems of water flow in the unsaturated zone numerically with rather fine spatial and temporal discretization for the correct mass balance and numerical stability, but also by the intensive computational task of quantifying geochemical reactions. In the present study, a flexible and efficient tool for large scale reactive transport modeling in variably saturated porous media and its applications are presented. The open source scientific software OpenGeoSys (OGS) is coupled with the IPhreeqc module of the geochemical solver PHREEQC. The new coupling approach makes full use of advantages from both codes: OGS provides a flexible choice of different numerical approaches for simulation of water flow in the vadose zone such as the pressure-based or mixed forms of Richards equation; whereas the IPhreeqc module leads to a simplification of data storage and its communication with OGS, which greatly facilitates the coupling and code updating. Moreover, a parallelization scheme with MPI (Message Passing Interface) is applied, in which the computational task of water flow and mass transport is partitioned through domain decomposition, whereas the efficient parallelization of geochemical reactions is achieved by smart allocation of computational workload over multiple compute nodes. The plausibility of the new coupling is verified by several benchmark tests. In addition, the efficiency of the new coupling approach is demonstrated by its application in a large scale scenario, in which the environmental fate of pesticides in a complex soil-aquifer system is studied.

Reactive transport modeling in variably saturated porous media with OGS-IPhreeqc

NASA Astrophysics Data System (ADS)

He, W.; Beyer, C.; Fleckenstein, J. H.; Jang, E.; Kalbacher, T.; Shao, H.; Wang, W.; Kolditz, O.

2014-12-01

Worldwide, sustainable water resource management becomes an increasingly challenging task due to the growth of population and extensive applications of fertilizer in agriculture. Moreover, climate change causes further stresses to both water quantity and quality. Reactive transport modeling in the coupled soil-aquifer system is a viable approach to assess the impacts of different land use and groundwater exploitation scenarios on the water resources. However, the application of this approach is usually limited in spatial scale and to simplified geochemical systems due to the huge computational expense involved. Such computational expense is not only caused by solving the high non-linearity of the initial boundary value problems of water flow in the unsaturated zone numerically with rather fine spatial and temporal discretization for the correct mass balance and numerical stability, but also by the intensive computational task of quantifying geochemical reactions. In the present study, a flexible and efficient tool for large scale reactive transport modeling in variably saturated porous media and its applications are presented. The open source scientific software OpenGeoSys (OGS) is coupled with the IPhreeqc module of the geochemical solver PHREEQC. The new coupling approach makes full use of advantages from both codes: OGS provides a flexible choice of different numerical approaches for simulation of water flow in the vadose zone such as the pressure-based or mixed forms of Richards equation; whereas the IPhreeqc module leads to a simplification of data storage and its communication with OGS, which greatly facilitates the coupling and code updating. Moreover, a parallelization scheme with MPI (Message Passing Interface) is applied, in which the computational task of water flow and mass transport is partitioned through domain decomposition, whereas the efficient parallelization of geochemical reactions is achieved by smart allocation of computational workload over multiple compute nodes. The plausibility of the new coupling is verified by several benchmark tests. In addition, the efficiency of the new coupling approach is demonstrated by its application in a large scale scenario, in which the environmental fate of pesticides in a complex soil-aquifer system is studied.

Solid phase studies and geochemical modelling of low-cost permeable reactive barriers.

PubMed

Bartzas, Georgios; Komnitsas, Kostas

2010-11-15

A continuous column experiment was carried out under dynamic flow conditions in order to study the efficiency of low-cost permeable reactive barriers (PRBs) to remove several inorganic contaminants from acidic solutions. A 50:50 w/w waste iron/sand mixture was used as candidate reactive media in order to activate precipitation and promote sorption and reduction-oxidation mechanisms. Solid phase studies of the exhausted reactive products after column shutdown, using scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR) and X-ray diffraction (XRD), confirmed that the principal Fe corrosion products identified in the reactive zone are amorphous iron (hydr)oxides (maghemite/magnetite and goethite), intermediate products (sulfate green rust), and amorphous metal sulfides such as amFeS and/or mackinawite. Geochemical modelling of the metal removal processes, including interactions between reactive media, heavy metal ions and sulfates, and interpretation of the ionic profiles was also carried out by using the speciation/mass transfer computer code PHREEQC-2 and the WATEQ4F database. Mineralogical characterization studies as well as geochemical modelling calculations also indicate that the effect of sulfate and silica sand on the efficiency of the reactive zone should be considered carefully during design and operation of low-cost field PRBs. Copyright © 2010 Elsevier B.V. All rights reserved.

Analysis of a Complex Faulted CO 2 Reservoir Using a Three-dimensional Hydro-geochemical-Mechanical Approach

DOE PAGES

Nguyen, Ba Nghiep; Hou, Zhangshuan; Bacon, Diana H.; ...

2017-08-18

This work applies a three-dimensional (3D) multiscale approach recently developed to analyze a complex CO 2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults. The approach couples the STOMP-CO2-R code for flow and reactive transport modeling to the ABAQUS ® finite element package for geomechanical analysis. The objective is to examine the coupled hydro-geochemical-mechanical impact on the risk of hydraulic fracture and fault slip in a complex and representative CO 2 reservoir that contains two nearly parallel faults. STOMP-CO2-R/ABAQUS ® coupled analyses of this reservoir are performed assuming extensional and compressional stress regimesmore » to predict evolutions of fluid pressure, stress and strain distributions as well as potential fault failure and leakage of CO 2 along the fault damage zones. The tendency for the faults to slip and pressure margin to fracture are examined in terms of stress regime, mineral composition, crack distributions in the fault damage zones and geomechanical properties. Here, this model in combination with a detailed description of the faults helps assess the coupled hydro-geochemical-mechanical effect.« less

Analysis of a Complex Faulted CO 2 Reservoir Using a Three-dimensional Hydro-geochemical-Mechanical Approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nguyen, Ba Nghiep; Hou, Zhangshuan; Bacon, Diana H.

This work applies a three-dimensional (3D) multiscale approach recently developed to analyze a complex CO 2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults. The approach couples the STOMP-CO2-R code for flow and reactive transport modeling to the ABAQUS ® finite element package for geomechanical analysis. The objective is to examine the coupled hydro-geochemical-mechanical impact on the risk of hydraulic fracture and fault slip in a complex and representative CO 2 reservoir that contains two nearly parallel faults. STOMP-CO2-R/ABAQUS ® coupled analyses of this reservoir are performed assuming extensional and compressional stress regimesmore » to predict evolutions of fluid pressure, stress and strain distributions as well as potential fault failure and leakage of CO 2 along the fault damage zones. The tendency for the faults to slip and pressure margin to fracture are examined in terms of stress regime, mineral composition, crack distributions in the fault damage zones and geomechanical properties. Here, this model in combination with a detailed description of the faults helps assess the coupled hydro-geochemical-mechanical effect.« less

Nonstationary porosity evolution in mixing zone in coastal carbonate aquifer using an alternative modeling approach.

PubMed

Laabidi, Ezzeddine; Bouhlila, Rachida

2015-07-01

In the last few decades, hydrogeochemical problems have benefited from the strong interest in numerical modeling. One of the most recognized hydrogeochemical problems is the dissolution of the calcite in the mixing zone below limestone coastal aquifer. In many works, this problem has been modeled using a coupling algorithm between a density-dependent flow model and a geochemical model. A related difficulty is that, because of the high nonlinearity of the coupled set of equations, high computational effort is needed. During calcite dissolution, an increase in permeability can be identified, which can induce an increase in the penetration of the seawater into the aquifer. The majority of the previous studies used a fully coupled reactive transport model in order to model such problem. Romanov and Dreybrodt (J Hydrol 329:661-673, 2006) have used an alternative approach to quantify the porosity evolution in mixing zone below coastal carbonate aquifer at steady state. This approach is based on the analytic solution presented by Phillips (1991) in his book Flow and Reactions in Permeable Rock, which shows that it is possible to decouple the complex set of equation. This equation is proportional to the square of the salinity gradient, which can be calculated using a density driven flow code and to the reaction rate that can be calculated using a geochemical code. In this work, this equation is used in nonstationary step-by-step regime. At each time step, the quantity of the dissolved calcite is quantified, the change of porosity is calculated, and the permeability is updated. The reaction rate, which is the second derivate of the calcium equilibrium concentration in the equation, is calculated using the PHREEQC code (Parkhurst and Apello 1999). This result is used in GEODENS (Bouhlila 1999; Bouhlila and Laabidi 2008) to calculate change of the porosity after calculating the salinity gradient. For the next time step, the same protocol is used but using the updated porosity and permeability distributions.

Modelling of the physico-chemical behaviour of clay minerals with a thermo-kinetic model taking into account particles morphology in compacted material.

NASA Astrophysics Data System (ADS)

Sali, D.; Fritz, B.; Clément, C.; Michau, N.

2003-04-01

Modelling of fluid-mineral interactions is largely used in Earth Sciences studies to better understand the involved physicochemical processes and their long-term effect on the materials behaviour. Numerical models simplify the processes but try to preserve their main characteristics. Therefore the modelling results strongly depend on the data quality describing initial physicochemical conditions for rock materials, fluids and gases, and on the realistic way of processes representations. The current geo-chemical models do not well take into account rock porosity and permeability and the particle morphology of clay minerals. In compacted materials like those considered as barriers in waste repositories, low permeability rocks like mudstones or compacted powders will be used : they contain mainly fine particles and the geochemical models used for predicting their interactions with fluids tend to misjudge their surface areas, which are fundamental parameters in kinetic modelling. The purpose of this study was to improve how to take into account the particles morphology in the thermo-kinetic code KINDIS and the reactive transport code KIRMAT. A new function was integrated in these codes, considering the reaction surface area as a volume depending parameter and the calculated evolution of the mass balance in the system was coupled with the evolution of reactive surface areas. We made application exercises for numerical validation of these new versions of the codes and the results were compared with those of the pre-existing thermo-kinetic code KINDIS. Several points are highlighted. Taking into account reactive surface area evolution during simulation modifies the predicted mass transfers related to fluid-minerals interactions. Different secondary mineral phases are also observed during modelling. The evolution of the reactive surface parameter helps to solve the competition effects between different phases present in the system which are all able to fix the chemical elements mobilised by the water-minerals interaction processes. To validate our model we simulated the compacted bentonite (MX80) studied for engineered barriers for radioactive waste confinement and mainly composed of Na-Ca-montmorillonite. The study of particles morphology and reactive surfaces evolutions reveals that aqueous ions have a complex behaviour, especially when competitions between various mineral phases occur. In that case, our model predicts a preferential precipitation of finest particles, favouring smectites instead of zeolites. This work is a part of a PhD Thesis supported by Andra, the French Radioactive Waste Management Agency.

Geochemical evolution of groundwater in the Culebra dolomite near the Waste Isolation Pilot Plant, southeastern New Mexico, USA

USGS Publications Warehouse

Siegel, M.D.; Anderholm, S.

1994-01-01

The Culebra Dolomite Member of the Rustler Formation, a thin (10 m) fractured dolomite aquifer, lies approximately 450 m above the repository horizon of the Waste Isolation Pilot Plant (WIPP) in southeastern New Mexico, USA. Salinities of water in the Culebra range roughly from 10,000 to 200,000 mg/L within the WIPP site. A proposed model for the post-Pleistocene hydrochemical evolution of the Culebra tentatively identifies the major sources and sinks for many of the groundwater solutes. Reaction-path simulations with the PHRQPITZ code suggest that the Culebra dolomite is a partial chemical equilibrium system whose composition is controlled by an irreversible process (dissolution of evaporites) and equilibrium with gypsum and calcite. Net geochemical reactions along postulated modern flow paths, calculated with the NETPATH code, include dissolution of halite, carbonate and evaporite salts, and ion exchange. R-mode principal component analysis revealed correlations among the concentrations of Si, Mg, pH, Li, and B that are consistent with several clay-water reactions. The results of the geochemical calculations and mineralogical data are consistent with the following hydrochemical model: 1. (1) solutes are added to the Culebra by dissolution of evaporite minerals 2. (2) the solubilities of gypsum and calcite increase as the salinity increases; these minerals dissolve as chemical equilibrium is maintained between them and the groundwater 3. (3) equilibrium is not maintained between the waters and dolomite; sufficient Mg is added to the waters by dissolution of accessory carnallite or polyhalite such that the degree of dolomite supersaturation increases with ionic strength 4. (4) clays within the fractures and rock matrix exert some control on the distribution of Li, B, Mg, and Si via sorption, ion exchange, and dissolution. ?? 1994.

OpenGeoSys-GEMS: Hybrid parallelization of a reactive transport code with MPI and threads

NASA Astrophysics Data System (ADS)

Kosakowski, G.; Kulik, D. A.; Shao, H.

2012-04-01

OpenGeoSys-GEMS is a generic purpose reactive transport code based on the operator splitting approach. The code couples the Finite-Element groundwater flow and multi-species transport modules of the OpenGeoSys (OGS) project (http://www.ufz.de/index.php?en=18345) with the GEM-Selektor research package to model thermodynamic equilibrium of aquatic (geo)chemical systems utilizing the Gibbs Energy Minimization approach (http://gems.web.psi.ch/). The combination of OGS and the GEM-Selektor kernel (GEMS3K) is highly flexible due to the object-oriented modular code structures and the well defined (memory based) data exchange modules. Like other reactive transport codes, the practical applicability of OGS-GEMS is often hampered by the long calculation time and large memory requirements. • For realistic geochemical systems which might include dozens of mineral phases and several (non-ideal) solid solutions the time needed to solve the chemical system with GEMS3K may increase exceptionally. • The codes are coupled in a sequential non-iterative loop. In order to keep the accuracy, the time step size is restricted. In combination with a fine spatial discretization the time step size may become very small which increases calculation times drastically even for small 1D problems. • The current version of OGS is not optimized for memory use and the MPI version of OGS does not distribute data between nodes. Even for moderately small 2D problems the number of MPI processes that fit into memory of up-to-date workstations or HPC hardware is limited. One strategy to overcome the above mentioned restrictions of OGS-GEMS is to parallelize the coupled code. For OGS a parallelized version already exists. It is based on a domain decomposition method implemented with MPI and provides a parallel solver for fluid and mass transport processes. In the coupled code, after solving fluid flow and solute transport, geochemical calculations are done in form of a central loop over all finite element nodes with calls to GEMS3K and consecutive calculations of changed material parameters. In a first step the existing MPI implementation was utilized to parallelize this loop. Calculations were split between the MPI processes and afterwards data was synchronized by using MPI communication routines. Furthermore, multi-threaded calculation of the loop was implemented with help of the boost thread library (http://www.boost.org). This implementation provides a flexible environment to distribute calculations between several threads. For each MPI process at least one and up to several dozens of worker threads are spawned. These threads do not replicate the complete OGS-GEM data structure and use only a limited amount of memory. Calculation of the central geochemical loop is shared between all threads. Synchronization between the threads is done by barrier commands. The overall number of local threads times MPI processes should match the number of available computing nodes. The combination of multi-threading and MPI provides an effective and flexible environment to speed up OGS-GEMS calculations while limiting the required memory use. Test calculations on different hardware show that for certain types of applications tremendous speedups are possible.

Hydro-geochemical modeling of the spatial and the temporal geochemical variations of the granitic Strengbach catchment springs (Vosges massif, France)

NASA Astrophysics Data System (ADS)

Ackerer, Julien; Chabaux, François; Lucas, Yann; Pierret, Marie Claire; Viville, Daniel; Fritz, Bertrand; Clement, Alain; Beaulieu, Emilie; Negrel, Philippe

2017-04-01

Regular analysis of the major element concentrations in waters from springs emerging on the Strengbach catchment is made for more than 20 years (OHGE, Observatoire Hydro-Géochimique de l'Environnement). These data confirm the spatial variability of geochemical characteristics of the Strengbach springs linked, at least partly, to the lithological variability of the substratum (Pierret et al., 2014). The data also indicate that at the first order, the geochemical fluxes exported from each spring are mainly linked to the spring discharges, without significant variations of the relationships linking these two parameters between 1990 and 2010. There is also no observation of significant variations for the dissolved silica and for most of the cationic concentrations with time. Only a significant decrease of the Ca concentrations is observed for the Strengbach springs from 1990 to 2010. Numerical simulations, performed with the KIRMAT hydro-geochemical code, show that such a decrease can be considered as the response in the "bedrock" of the water-rock interactions to the variations of the soil solution chemical compositions recorded over the last 20 years, marked by a significant increase of pH and decrease of Ca concentrations. In particular, the modeling results show that the Ca concentration decrease is controlled by the couple apatite/clays, and that significant modifications of the apatite dissolution rate and clay compositions occurred between 1990 and 2010. This study shows that the temporal evolution of the Strengbach spring chemistry cannot be explained by the only variations of the clay mineral compositions, i.e. a modification of the chemical composition of the precipitated clays or a modification of the ionic exchange capacity of the clay minerals, but that it is definitely the interrelations between the apatite and the clay minerals that are involved.

An interactive code (NETPATH) for modeling NET geochemical reactions along a flow PATH, version 2.0

USGS Publications Warehouse

Plummer, Niel; Prestemon, Eric C.; Parkhurst, David L.

1994-01-01

NETPATH is an interactive Fortran 77 computer program used to interpret net geochemical mass-balance reactions between an initial and final water along a hydrologic flow path. Alternatively, NETPATH computes the mixing proportions of two to five initial waters and net geochemical reactions that can account for the observed composition of a final water. The program utilizes previously defined chemical and isotopic data for waters from a hydrochemical system. For a set of mineral and (or) gas phases hypothesized to be the reactive phases in the system, NETPATH calculates the mass transfers in every possible combination of the selected phases that accounts for the observed changes in the selected chemical and (or) isotopic compositions observed along the flow path. The calculations are of use in interpreting geochemical reactions, mixing proportions, evaporation and (or) dilution of waters, and mineral mass transfer in the chemical and isotopic evolution of natural and environmental waters. Rayleigh distillation calculations are applied to each mass-balance model that satisfies the constraints to predict carbon, sulfur, nitrogen, and strontium isotopic compositions at the end point, including radiocarbon dating. DB is an interactive Fortran 77 computer program used to enter analytical data into NETPATH, and calculate the distribution of species in aqueous solution. This report describes the types of problems that can be solved, the methods used to solve problems, and the features available in the program to facilitate these solutions. Examples are presented to demonstrate most of the applications and features of NETPATH. The codes DB and NETPATH can be executed in the UNIX or DOS1 environment. This report replaces U.S. Geological Survey Water-Resources Investigations Report 91-4078, by Plummer and others, which described the original release of NETPATH, version 1.0 (dated December, 1991), and documents revisions and enhancements that are included in version 2.0. 1 The use of trade, brand or product names in this report is for identification purposes only and does not constitute endorsement by the U.S. Geological Survey.

Modelling of Pesticide Transport During An Injection Experiment In A Physical and Geochemical Heterogeneous Aquifer

NASA Astrophysics Data System (ADS)

Hojberg, A. L.; Engesgaard, P.; Bjerg, P. L.

The fate of selected pesticides under natural groundwater conditions was studied by natural gradient short and long term injection experiments in a shallow uncon- fined aerobic aquifer. Bentazone, DNOC, MCPP, dichlorprop, isoproturon, and BAM (dichlobenil metabolite) were injected in aqueous solution with bromide as a nonre- active tracer. The Bromide and pesticide plumes were sampled during the initial 25 m of migration in a dense monitoring net of multilevel samplers. The aquifer was physical and geochemical heterogeneous, which affected transport of several of the pesticides. A 3D reactive transport code was developed including one- and two-site linear/nonlinear equilibrium/nonequilibrium sorption and first-order as well as single Monod degradation kinetic coupled to microbial growth. Model simulations demon- strated that microbial growth was likely supported by the phenoxy acids MCPP and dichlorprop, while degradation of DNOC was adequately described by first-order degradation with no initial lag time. An observed vertical increase in pH was observed at the site and implemented in the transport code. The numerical analysis indicated that degradation of the three degradable pesticides may have been affected by vertical pH variations. Spatial variability in observed DNOC sorption was similarly suspected to be an effect of varying pH. pH dependency on DNOC sorption was confirmed by the model recognized by a match to observed breakthrough at the individual sampling points, when pH variation was included in the simulations.

Analyzing legacy U.S. Geological Survey geochemical databases using GIS: applications for a national mineral resource assessment

USGS Publications Warehouse

Yager, Douglas B.; Hofstra, Albert H.; Granitto, Matthew

2012-01-01

This report emphasizes geographic information system analysis and the display of data stored in the legacy U.S. Geological Survey National Geochemical Database for use in mineral resource investigations. Geochemical analyses of soils, stream sediments, and rocks that are archived in the National Geochemical Database provide an extensive data source for investigating geochemical anomalies. A study area in the Egan Range of east-central Nevada was used to develop a geographic information system analysis methodology for two different geochemical datasets involving detailed (Bureau of Land Management Wilderness) and reconnaissance-scale (National Uranium Resource Evaluation) investigations. ArcGIS was used to analyze and thematically map geochemical information at point locations. Watershed-boundary datasets served as a geographic reference to relate potentially anomalous sample sites with hydrologic unit codes at varying scales. The National Hydrography Dataset was analyzed with Hydrography Event Management and ArcGIS Utility Network Analyst tools to delineate potential sediment-sample provenance along a stream network. These tools can be used to track potential upstream-sediment-contributing areas to a sample site. This methodology identifies geochemically anomalous sample sites, watersheds, and streams that could help focus mineral resource investigations in the field.

A 1-D evolutionary model for icy satellites, applied to Enceladus

NASA Astrophysics Data System (ADS)

Prialnik, Dina; Malamud, Uri

2015-11-01

A 1-D long-term evolution code for icy satellites is presented, which couples multiple processes: water migration, geochemical reactions, water and silicate phase transitions, crystallization, compaction by self-gravity, and ablation. The code takes into account various energy sources: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy, and insolation. It includes heat transport by conduction, convection, and advection.The code is applied to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, and adopting a homogeneous initial structure. Assuming that the satellite has been losing water continually along its evolution, it follows that it was formed as a more massive, more ice-rich and more porous object, and gradually transformed into its present day state, due to sustained tidal heating. Several initial compositions and evolution scenarios are considered, and the evolution is simulated for the age of the Solar System. The results corresponding to the present configuration are confronted with the available observational constraints. The present configuration is shown to be differentiated into a pure icy mantle, several tens of km thick, overlying a rocky core, composed of dehydrated rock in the central part and hydrated rock in the outer part. Such a differentiated structure is obtained not only for Enceladus, but for other medium size ice-rich bodies as well.Predictions for Enceladus are a higher rock/ice mass ratio than previously assumed, and a thinner ice mantle, compatible with recent estimates based on gravity field measurements. Although, obviously, the 1-D model cannot be used to explain local phenomena, it sheds light on the internal structure invoked in explanations of localized features and activities.

Radiological performance assessment for the E-Area Vaults Disposal Facility. Appendices A through M

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cook, J.R.

1994-04-15

These document contains appendices A-M for the performance assessment. They are A: details of models and assumptions, B: computer codes, C: data tabulation, D: geochemical interactions, E: hydrogeology of the Savannah River Site, F: software QA plans, G: completeness review guide, H: performance assessment peer review panel recommendations, I: suspect soil performance analysis, J: sensitivity/uncertainty analysis, K: vault degradation study, L: description of naval reactor waste disposal, M: porflow input file. (GHH)

Methods for Estimating Physicochemical Properties of Inorganic Chemicals of Environmental Concern. Volume 1

DTIC Science & Technology

1987-03-01

available. Kincaid et al. (221 provide an overview of the geochemical code history (up to 1983) in which they group models into four major families...plant. The *Ixample is adapted from a study by Morel ,’t al. 1291 for Los Angeles County. It is prv.ented in four ;ections: I ’ Problem statement, ý2...presents environmentally important priperty dJata for several groups of elements or compcui’ds. The data include, for Pxrimpie, precal- culated

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nguyen, Ba Nghiep; Hou, Zhangshuan; Bacon, Diana H.

This paper applies a multiscale hydro-geochemical-mechanical approach to analyze faulted CO 2 reservoirs using the STOMP-CO 2-R code that is coupled to the ABAQUS® finite element package. STOMP-CO 2-R models the reactive transport of CO 2 causing mineral volume fraction changes that are captured by an Eshelby-Mori-Tanka model implemented in ABAQUS®. A three-dimensional (3D) STOMP-CO 2-R model for a reservoir containing an inclined fault was built to analyze a formation containing a reaction network with 5 minerals: albite, anorthite, calcite, kaolinite and quartz. A 3D finite element mesh that exactly maps the STOMP-CO 2-R grid is developed for coupled hydro-geochemical-mechanicalmore » analyses. The model contains alternating sandstone and shale layers. The impact of reactive transport of CO 2 on the geomechanical properties of reservoir rocks and seals are studied in terms of mineral composition changes that affect their geomechanical responses. Simulations assuming extensional and compressional stress regimes with and without coupled geochemistry are performed to study the stress regime effect on the risk of hydraulic fracture. The tendency for the fault to slip is examined in terms of stress regime, geomechanical and geochemical-mechanical effects as well as fault inclination. The results show that mineralogical changes due to long-term injection of CO 2 reduce the permeability and elastic modulus of the reservoir, leading to increased risk of hydraulic fracture in the injection location and at the caprock seal immediately above the injection zone. Fault slip is not predicted to occur. However, fault inclination and stress regime have an important impact on the slip tendency factor.« less

A Cloud Based Framework For Monitoring And Predicting Subsurface System Behaviour

NASA Astrophysics Data System (ADS)

Versteeg, R. J.; Rodzianko, A.; Johnson, D. V.; Soltanian, M. R.; Dwivedi, D.; Dafflon, B.; Tran, A. P.; Versteeg, O. J.

2015-12-01

Subsurface system behavior is driven and controlled by the interplay of physical, chemical, and biological processes which occur at multiple temporal and spatial scales. Capabilities to monitor, understand and predict this behavior in an effective and timely manner are needed for both scientific purposes and for effective subsurface system management. Such capabilities require three elements: Models, Data and an enabling cyberinfrastructure, which allow users to use these models and data in an effective manner. Under a DOE Office of Science funded STTR award Subsurface Insights and LBNL have designed and implemented a cloud based predictive assimilation framework (PAF) which automatically ingests, controls quality and stores heterogeneous physical and chemical subsurface data and processes these data using different inversion and modeling codes to provide information on the current state and evolution of subsurface systems. PAF is implemented as a modular cloud based software application with five components: (1) data acquisition, (2) data management, (3) data assimilation and processing, (4) visualization and result delivery and (5) orchestration. Serverside PAF uses ZF2 (a PHP web application framework) and Python and both open source (ODM2) and in house developed data models. Clientside PAF uses CSS and JS to allow for interactive data visualization and analysis. Client side modularity (which allows for a responsive interface) of the system is achieved by implementing each core capability of PAF (such as data visualization, user configuration and control, electrical geophysical monitoring and email/SMS alerts on data streams) as a SPA (Single Page Application). One of the recent enhancements is the full integration of a number of flow and mass transport and parameter estimation codes (e.g., MODFLOW, MT3DMS, PHT3D, TOUGH, PFLOTRAN) in this framework. This integration allows for autonomous and user controlled modeling of hydrological and geochemical processes. In our presentation we will discuss our software architecture and present the results of using these codes and the overall developed performance of our framework using hydrological, geochemical and geophysical data from the LBNL SFA2 Rifle field site.

A computer program incorporating Pitzer's equations for calculation of geochemical reactions in brines

USGS Publications Warehouse

Plummer, Niel; Parkhurst, D.L.; Fleming, G.W.; Dunkle, S.A.

1988-01-01

The program named PHRQPITZ is a computer code capable of making geochemical calculations in brines and other electrolyte solutions to high concentrations using the Pitzer virial-coefficient approach for activity-coefficient corrections. Reaction-modeling capabilities include calculation of (1) aqueous speciation and mineral-saturation index, (2) mineral solubility, (3) mixing and titration of aqueous solutions, (4) irreversible reactions and mineral water mass transfer, and (5) reaction path. The computed results for each aqueous solution include the osmotic coefficient, water activity , mineral saturation indices, mean activity coefficients, total activity coefficients, and scale-dependent values of pH, individual-ion activities and individual-ion activity coeffients , and scale-dependent values of pH, individual-ion activities and individual-ion activity coefficients. A data base of Pitzer interaction parameters is provided at 25 C for the system: Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O, and extended to include largely untested literature data for Fe(II), Mn(II), Sr, Ba, Li, and Br with provision for calculations at temperatures other than 25C. An extensive literature review of published Pitzer interaction parameters for many inorganic salts is given. Also described is an interactive input code for PHRQPITZ called PITZINPT. (USGS)

Decreasing Kd uncertainties through the application of thermodynamic sorption models.

PubMed

Domènech, Cristina; García, David; Pękala, Marek

2015-09-15

Radionuclide retardation processes during transport are expected to play an important role in the safety assessment of subsurface disposal facilities for radioactive waste. The linear distribution coefficient (Kd) is often used to represent radionuclide retention, because analytical solutions to the classic advection-diffusion-retardation equation under simple boundary conditions are readily obtainable, and because numerical implementation of this approach is relatively straightforward. For these reasons, the Kd approach lends itself to probabilistic calculations required by Performance Assessment (PA) calculations. However, it is widely recognised that Kd values derived from laboratory experiments generally have a narrow field of validity, and that the uncertainty of the Kd outside this field increases significantly. Mechanistic multicomponent geochemical simulators can be used to calculate Kd values under a wide range of conditions. This approach is powerful and flexible, but requires expert knowledge on the part of the user. The work presented in this paper aims to develop a simplified approach of estimating Kd values whose level of accuracy would be comparable with those obtained by fully-fledged geochemical simulators. The proposed approach consists of deriving simplified algebraic expressions by combining relevant mass action equations. This approach was applied to three distinct geochemical systems involving surface complexation and ion-exchange processes. Within bounds imposed by model simplifications, the presented approach allows radionuclide Kd values to be estimated as a function of key system-controlling parameters, such as the pH and mineralogy. This approach could be used by PA professionals to assess the impact of key geochemical parameters on the variability of radionuclide Kd values. Moreover, the presented approach could be relatively easily implemented in existing codes to represent the influence of temporal and spatial changes in geochemistry on Kd values. Copyright © 2015 Elsevier B.V. All rights reserved.

Modeling porosity reductions caused by mineral fouling in continuous-wall permeable reactive barriers.

PubMed

Li, Lin; Benson, Craig H; Lawson, Elizabeth M

2006-02-01

A study was conducted to assess key factors to include when modeling porosity reductions caused by mineral fouling in permeable reactive barriers (PRBs) containing granular zero valent iron. The public domain codes MODFLOW and RT3D were used and a geochemical algorithm was developed for RT3D to simulate geochemical reactions occurring in PRBs. Results of simulations conducted with the model show that the largest porosity reductions occur between the entrance and mid-plane of the PRB as a result of precipitation of carbonate minerals and that smaller porosity reductions occur between the mid-plane and exit face due to precipitation of ferrous hydroxide. These findings are consistent with field and laboratory observations, as well as modeling predictions made by others. Parametric studies were conducted to identify the most important variables to include in a model evaluating porosity reduction. These studies showed that three minerals (CaCO3, FeCO3, and Fe(OH)2 (am)) account for more than 99% of the porosity reductions that were predicted. The porosity reduction is sensitive to influent concentrations of HCO3-, Ca2+, CO3(2-), and dissolved oxygen, the anaerobic iron corrosion rate, and the rates of CaCO3 and FeCO3 formation. The predictions also show that porosity reductions in PRBs can be spatially variable and mineral forming ions penetrate deeper into the PRB as a result of flow heterogeneities, which reflects the balance between the rate of mass transport and geochemical reaction rates. Level of aquifer heterogeneity and the contrast in hydraulic conductivity between the aquifer and PRB are the most important hydraulic variables affecting porosity reduction. Spatial continuity of aquifer hydraulic conductivity is less significant.

Implementation of agronomical and geochemical modules into a 3D groundwater code for assessing nitrate storage and transport through unconfined Chalk aquifer

NASA Astrophysics Data System (ADS)

Picot-Colbeaux, Géraldine; Devau, Nicolas; Thiéry, Dominique; Pettenati, Marie; Surdyk, Nicolas; Parmentier, Marc; Amraoui, Nadia; Crastes de Paulet, François; André, Laurent

2016-04-01

Chalk aquifer is the main water resource for domestic water supply in many parts in northern France. In same basin, groundwater is frequently affected by quality problems concerning nitrates. Often close to or above the drinking water standards, nitrate concentration in groundwater is mainly due to historical agriculture practices, combined with leakage and aquifer recharge through the vadose zone. The complexity of processes occurring into such an environment leads to take into account a lot of knowledge on agronomy, geochemistry and hydrogeology in order to understand, model and predict the spatiotemporal evolution of nitrate content and provide a decision support tool for the water producers and stakeholders. To succeed in this challenge, conceptual and numerical models representing accurately the Chalk aquifer specificity need to be developed. A multidisciplinary approach is developed to simulate storage and transport from the ground surface until groundwater. This involves a new agronomic module "NITRATE" (NItrogen TRansfer for Arable soil to groundwaTEr), a soil-crop model allowing to calculate nitrogen mass balance in arable soil, and the "PHREEQC" numerical code for geochemical calculations, both coupled with the 3D transient groundwater numerical code "MARTHE". Otherwise, new development achieved on MARTHE code allows the use of dual porosity and permeability calculations needed in the fissured Chalk aquifer context. This method concerning the integration of existing multi-disciplinary tools is a real challenge to reduce the number of parameters by selecting the relevant equations and simplifying the equations without altering the signal. The robustness and the validity of these numerical developments are tested step by step with several simulations constrained by climate forcing, land use and nitrogen inputs over several decades. In the first time, simulations are performed in a 1D vertical unsaturated soil column for representing experimental nitrates vertical soil profiles (0-30m depth experimental measurements in Somme region). In the second time, this approach is used to simulate with a 3D model a drinking water catchment area in order to compared nitrate contents time series calculated and measured in the domestic water pumping well since 1995 (field in northern France - Avre Basin region). This numerical tool will help the decision-making in all activities in relation with water uses.

SilMush: A procedure for modeling of the geochemical evolution of silicic magmas and granitic rocks

NASA Astrophysics Data System (ADS)

Hertogen, Jan; Mareels, Joyce

2016-07-01

A boundary layer crystallization modeling program is presented that specifically addresses the chemical fractionation in silicic magma systems and the solidification of plutonic bodies. The model is a Langmuir (1989) type approach and does not invoke crystal settling in high-viscosity silicic melts. The primary aim is to model a granitic rock as a congealed crystal-liquid mush, and to integrate major element and trace element modeling. The procedure allows for some exploratory investigation of the exsolution of H2O-fluids and of the fluid/melt partitioning of trace elements. The procedure is implemented as a collection of subroutines for the MS Excel spreadsheet environment and is coded in the Visual Basic for Applications (VBA) language. To increase the flexibility of the modeling, the procedure is based on discrete numeric process simulation rather than on solution of continuous differential equations. The program is applied to a study of the geochemical variation within and among three granitic units (Senones, Natzwiller, Kagenfels) from the Variscan Northern Vosges Massif, France. The three units cover the compositional range from monzogranite, over syenogranite to alkali-feldspar granite. An extensive set of new major element and trace element data is presented. Special attention is paid to the essential role of accessory minerals in the fractionation of the Rare Earth Elements. The crystallization model is able to reproduce the essential major and trace element variation trends in the data sets of the three separate granitic plutons. The Kagenfels alkali-feldspar leucogranite couples very limited variation in major element composition to a considerable and complex variation of trace elements. The modeling results can serve as a guide for the reconstruction of the emplacement sequence of petrographically distinct units. Although the modeling procedure essentially deals with geochemical fractionation within a single pluton, the modeling results bring up a number of questions about the petrogenetic relationships among parental magmas of nearly coeval granitic units emplaced in close proximity.

Physical model of a fumarolic system inferred from a high-resolution 3-D Resistivity image of Solfatara volcano

NASA Astrophysics Data System (ADS)

Gresse, Marceau; Vandemeulebrouck, Jean; Byrdina, Svetlana; Chiodini, Giovanni; Rinaldi, Antonio Pio; Johnson, Timothy C.; Ricci, Tullio; Petrillo, Zaccaria; Vilardo, Giuseppe; Lebourg, Thomas; Mangiacapra, Annarita

2017-04-01

Solfatara crater, located inside the Phlegrean Fields caldera, is showing a significant unrest activity since 10 years with a increase of ground deformation, degassing and heating. Electrical Resistivity Imaging was performed between 2012 and 2016 with the purpose of improving our knowledge of the shallow hydrothermal system. The complete dataset includes 43,432 D-C measurements inverted using the E4D code. This 3-D inversion was compared with the mappings of surface temperature, diffuse soil CO2 flux and self-potential in order to better constrain the interpretation of the observed resistivity structure in terms of lithological contrasts and hydrothermal signatures. For the first time, we highlighted in 3-D the main geological units: Monte Olibano lava dome and Solfatara crypto-dome appear as two relatively resistive bodies (50-100 Ω.m). Furthermore, the resistivity model clearly revealed the contrasting geometry of the hydrothermal circulation in the Solfatara crater. A channel-like conductive structure (7 Ω.m) represents the condensate that flows from the main fumarolic area down to the liquid-dominated Fangaia mud pool. This interpretation is consistent with the negative Self-Potential anomaly and with the surface observations. We imaged at a metric-resolution the two main fumaroles, Bocca Grande and Bocca Nuova, that have the following geochemical characteristics. Bocca Grande vent: 162°C, ˜150 t of CO2 released per day with a mass ratio CO2/H20 = 0.4 and Bocca Nuova vent: 148°C, ˜50 t of CO2 released per day with a mass ratio CO2/H20 = 0.45. The differences between these geochemical characteristics could lead one to believe that they are fed by two distinct sources at depth. On the contrary, our resistivity model shows that the two fumarolic vents are directly connected to a common resistive body (30-50 Ω.m) at a depth of 50 meters. This structure likely represents a single gas reservoir feeding the two fumaroles. Its depth corresponds indeed to a steam source at a pressure of 6 bar and at a temperature of least 165 °C. The geophysical images combined with the geochemical data allowed us to build up a multiphase fluid flow model of the Bocca Grande and and Bocca Nuova fumaroles using the TOUGH 2 code. Our results show that the distinct resistivity structure, temperature, and water content of the both fumaroles are due to the particular geometry of the condensate flow that intersects and contaminates the Bocca Nuova but not the Bocca Grande fumarole. These results indicate the necessity to combine geophysical and geochemical approaches in order to better apprehend the structure complexity and the dynamics of fumaroles and hydrothermal systems.

SEDIMENT GEOCHEMICAL MODEL

EPA Science Inventory

Until recently, sediment geochemical models (diagenetic models) have been only able to explain sedimentary flux and concentration profiles for a few simplified geochemical cycles (e.g., nitrogen, carbon and sulfur). However with advances in numerical methods, increased accuracy ...

Solid phase evolution in the Biosphere 2 hillslope experiment as predicted by modeling of hydrologic and geochemical fluxes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dontsova, K.; Steefel, C.I.; Desilets, S.

2009-07-15

A reactive transport geochemical modeling study was conducted to help predict the mineral transformations occurring over a ten year time-scale that are expected to impact soil hydraulic properties in the Biosphere 2 (B2) synthetic hillslope experiment. The modeling sought to predict the rate and extent of weathering of a granular basalt (selected for hillslope construction) as a function of climatic drivers, and to assess the feedback effects of such weathering processes on the hydraulic properties of the hillslope. Flow vectors were imported from HYDRUS into a reactive transport code, CrunchFlow2007, which was then used to model mineral weathering coupled tomore » reactive solute transport. Associated particle size evolution was translated into changes in saturated hydraulic conductivity using Rosetta software. We found that flow characteristics, including velocity and saturation, strongly influenced the predicted extent of incongruent mineral weathering and neo-phase precipitation on the hillslope. Results were also highly sensitive to specific surface areas of the soil media, consistent with surface reaction controls on dissolution. Effects of fluid flow on weathering resulted in significant differences in the prediction of soil particle size distributions, which should feedback to alter hillslope hydraulic conductivities.« less

Modeling the hydrothermal circulation and the hydrogen production at the Rainbow site with Cast3M

NASA Astrophysics Data System (ADS)

Perez, F.; Mügler, C.; Charlou, J.; Jean-baptiste, P.

2012-12-01

On the Mid-Atlantic Ridge, the Rainbow venting site is described as an ultramafic-hosted active hydrothermal site and releases high fluxes of methane and hydrogen [1, 2]. This behavior has first been interpreted as the result of serpentinization processes. But geochemical reactions involving olivine and plagioclase assemblages, and leading to chlorite, tremolite, talc and magnetite assemblages, could contribute to the observed characteristics of the exiting fluid [2]. The predominance of one of these geochemical reactions or their coexistence strongly depend on the hydrothermal fluid circulation. We developed and validated a 2D/3D numerical model using a Finite Volume method to simulate heat driven fluid flows in the framework of the Cast3M code [3, 4]. We also developed a numerical model for hydrogen production and transport that is based on experimental studies of the serpentinization processes [5-6]. This geochemical model takes into account the exothermic and water-consuming behavior of the serpentinization reaction and it can be coupled to our thermo-hydrogeological model. Our simulations provide temperatures, mass fluxes and venting surface areas very close to those estimated in-situ [7]. We showed that a single-path model [8] was necessary to simulate high values such as the in-situ measured temperatures and estimated water mass fluxes of the Rainbow site [7]. This single-path model will be used to model the production and transport of hydrogen at the Rainbow hydrothermal site. References [1]Charlou et al. (2010) AGU Monograph series. [2]Seyfried et al. (2011) Geochim. Cosmochim. Acta 75, 1574-1593. [3]http://www-cast3m.cea.fr. [4]Martin & Fyfe (1970) Chem. Geol. 6, 185-202. [5] Marcaillou et al. (2011) Earth and Planet. Sci. Lett. 303, 281-290. [6]Malvoisin et al. (2012) JGR, 117, B01104. [7]Perez et al. (2012) submited to Computational Geosciences. [8]Lowell & Germanovich (2004) AGU, Washington DC, USA.

A combined PHREEQC-2/parallel fracture model for the simulation of laminar/non-laminar flow and contaminant transport with reactions

NASA Astrophysics Data System (ADS)

Masciopinto, Costantino; Volpe, Angela; Palmiotta, Domenico; Cherubini, Claudia

2010-09-01

A combination of a parallel fracture model with the PHREEQC-2 geochemical model was developed to simulate sequential flow and chemical transport with reactions in fractured media where both laminar and turbulent flows occur. The integration of non-laminar flow resistances in one model produced relevant effects on water flow velocities, thus improving model prediction capabilities on contaminant transport. The proposed conceptual model consists of 3D rock-blocks, separated by horizontal bedding plane fractures with variable apertures. Particle tracking solved the transport equations for conservative compounds and provided input for PHREEQC-2. For each cluster of contaminant pathways, PHREEQC-2 determined the concentration for mass-transfer, sorption/desorption, ion exchange, mineral dissolution/precipitation and biodegradation, under kinetically controlled reactive processes of equilibrated chemical species. Field tests have been performed for the code verification. As an example, the combined model has been applied to a contaminated fractured aquifer of southern Italy in order to simulate the phenol transport. The code correctly fitted the field available data and also predicted a possible rapid depletion of phenols as a result of an increased biodegradation rate induced by a simulated artificial injection of nitrates, upgradient to the sources.

A Thermal Model for the Differentiation of Asteroid 4 Vesta, Based on Radiogenic and Collisional Heating

NASA Astrophysics Data System (ADS)

Ghosh, Amitabha

A finite element code has been developed to study the thermal history of asteroid 4 Vesta. This is the first attempt to model the thermal history of a differentiated asteroid, from accretion through core and crust formation and subsequent cooling until geochemical closure is attained. Previous thermal models were simpler formulations aimed at explaining metamorphism and aqueous alteration in unmelted asteroids. The results of the simulation are consistent with chronological measurements of cumulate and noncumulate eucrites, meteorites belonging to the HED suite, believed to have been derived from 4 Vesta. The work solves major problems with the hypothesis of heating by decay of 26Al, an extinct radionuclide, believed to be a plausible heat source in the early solar system. The simulation draws a model chronology of Vesta and predicts the time interval of accretion at 2.85 Myrs, the absolute times (with respect to CAI formation) of core formation at 4.58 Myrs, crust formation at 6.58 Myrs and geochemical closure on Vesta at ~100 Myrs. It is concluded that neither collisional heating nor heating due to the radioactive decay of 60Fe caused any perceptible difference in the whole-body thermal history of Vesta. Further, the thermal model suggested that the olivine-rich spot observed on Vesta may not be excavated mantle material, but may be unmelted near-surface material that escaped the asteroid's differentiation history.

Modelling of chemical degradation of blended cement-based materials by leaching cycles with Callovo-Oxfordian porewater

NASA Astrophysics Data System (ADS)

Olmeda, Javier; Henocq, Pierre; Giffaut, Eric; Grivé, Mireia

2017-06-01

The present work describes a thermodynamic model based on pore water replacement cycles to simulate the chemical evolution of blended cement (BFS + FA) by interaction with external Callovo-Oxfordian (COx) pore water. In the framework of the radioactive waste management, the characterization of the radionuclide behaviour (solubility/speciation, adsorption) in cementitious materials needs to be done for several chemical degradation states (I to IV). In particular, in the context of the deep geological radioactive waste disposal project (Cigéo), cement-based materials will be chemically evolved with time in contact with the host-rock (COx formation). The objective of this study is to provide an equilibrium solution composition for each degradation state for a CEM-V cement-based material to support the adsorption and diffusion experiments reproducing any state of degradation. Calculations have been performed at 25 °C using the geochemical code PhreeqC and an up-to-date thermodynamic database (ThermoChimie v.9.0.b) coupled to SIT approach for ionic strength correction. The model replicates experimental data with accuracy. The approach followed in this study eases the analysis of the chemical evolution in both aqueous and solid phase to obtain a fast assessment of the geochemical effects associated to an external water intrusion of variable composition on concrete structures.

COTHERM: Modelling fluid-rock interactions in Icelandic geothermal systems

NASA Astrophysics Data System (ADS)

Thien, Bruno; Kosakowski, Georg; Kulik, Dmitrii

2014-05-01

Mineralogical alteration of reservoir rocks, driven by fluid circulation in natural or enhanced geothermal systems, is likely to influence the long-term performance of geothermal power generation. A key factor is the change of porosity due to dissolution of primary minerals and precipitation of secondary phases. Porosity changes will affect fluid circulation and solute transport, which, in turn, influence mineralogical alteration. This study is part of the Sinergia COTHERM project (COmbined hydrological, geochemical and geophysical modeling of geotTHERMal systems) that is an integrative research project aimed at improving our understanding of the sub-surface processes in magmatically-driven natural geothermal systems. We model the mineralogical and porosity evolution of Icelandic geothermal systems with 1D and 2D reactive transport models. These geothermal systems are typically high enthalphy systems where a magmatic pluton is located at a few kilometers depth. The shallow plutons increase the geothermal gradient and trigger the circulation of hydrothermal waters with a steam cap forming at shallow depth. We investigate two contrasting geothermal systems: Krafla, for which the water recharge consists of meteoritic water; and Reykjanes, for which the water recharge mainly consists of seawater. The initial rock composition is a fresh basalt. We use the GEM-Selektor geochemical modeling package [1] for calculation of kinetically controlled mineral equilibria between the rock and the ingression water. We consider basalt minerals dissolution kinetics according to Palandri & Kharaka [2]. Reactive surface areas are assumed to be geometric surface areas, and are corrected using a spherical-particle surface/mass relationship. For secondary minerals, we consider the partial equilibrium assuming that the primary mineral dissolution is slow, and the secondary mineral precipitation is fast. Comparison of our modeling results with the mineralogical assemblages observed in the field by Gudmundsson & Arnorsson [3] and by Icelandic partners of the COTHERM project suggests that the concept of partial equilibrium with instantaneous precipitation of secondary minerals is not sufficient to satisfactorily describe the experimental data. Considering kinetic controls also for secondary minerals appears as indispensable to properly describe the geothermal system evolution using a reactive transport modelling approach [4]. [1] Kulik D.A., Wagner T., Dmytrieva S.V., Kosakowski G., Hingerl F.F., Chudnenko K.V., Berner U., 2013. GEM-Selektor geochemical modeling package: revised algorithm and GEMS3K numerical kernel for coupled simulation codes. Computational Geosciences 17, 1-24. http://gems.web.psi.ch. [2] Palandri, J.L., Kharaka, Y.K., 2004. A compilation of rate parameters of water-mineral interaction kinetics for application to geochemical modelling. U.S.Geological Survey, Menlo Park, CA, pp. 1-64. [3] Gudmundsson B.T., Arnorsson S., 2005. Secondary mineral-fluid equilibria in the Krafla and Namafjall geothermal systems, Iceland. Applied Geochememistry 20, 1607-1625. [4] Kosakowski, G., & Watanabe, N., 2013. OpenGeoSys-Gem: A numerical tool for calculating geochemical and porosity changes in saturated and partially saturated media. Physics and Chemistry of the Earth, Parts A/B/C. doi:10.1016/j.pce.2013.11.008

Generic reactive transport codes as flexible tools to integrate soil organic matter degradation models with water, transport and geochemistry in soils

NASA Astrophysics Data System (ADS)

Jacques, Diederik; Gérard, Fréderic; Mayer, Uli; Simunek, Jirka; Leterme, Bertrand

2016-04-01

A large number of organic matter degradation, CO2 transport and dissolved organic matter models have been developed during the last decades. However, organic matter degradation models are in many cases strictly hard-coded in terms of organic pools, degradation kinetics and dependency on environmental variables. The scientific input of the model user is typically limited to the adjustment of input parameters. In addition, the coupling with geochemical soil processes including aqueous speciation, pH-dependent sorption and colloid-facilitated transport are not incorporated in many of these models, strongly limiting the scope of their application. Furthermore, the most comprehensive organic matter degradation models are combined with simplified representations of flow and transport processes in the soil system. We illustrate the capability of generic reactive transport codes to overcome these shortcomings. The formulations of reactive transport codes include a physics-based continuum representation of flow and transport processes, while biogeochemical reactions can be described as equilibrium processes constrained by thermodynamic principles and/or kinetic reaction networks. The flexibility of these type of codes allows for straight-forward extension of reaction networks, permits the inclusion of new model components (e.g.: organic matter pools, rate equations, parameter dependency on environmental conditions) and in such a way facilitates an application-tailored implementation of organic matter degradation models and related processes. A numerical benchmark involving two reactive transport codes (HPx and MIN3P) demonstrates how the process-based simulation of transient variably saturated water flow (Richards equation), solute transport (advection-dispersion equation), heat transfer and diffusion in the gas phase can be combined with a flexible implementation of a soil organic matter degradation model. The benchmark includes the production of leachable organic matter and inorganic carbon in the aqueous and gaseous phases, as well as different decomposition functions with first-order, linear dependence or nonlinear dependence on a biomass pool. In addition, we show how processes such as local bioturbation (bio-diffusion) can be included implicitly through a Fickian formulation of transport of soil organic matter. Coupling soil organic matter models with generic and flexible reactive transport codes offers a valuable tool to enhance insights into coupled physico-chemical processes at different scales within the scope of C-biogeochemical cycles, possibly linked with other chemical elements such as plant nutrients and pollutants.

Geochemical Modeling of Carbon Sequestration, MMV, and EOR in the Illinois Basin

USGS Publications Warehouse

Berger, P.M.; Roy, W.R.; Mehnert, E.

2009-01-01

The Illinois State Geologic Survey is conducting several ongoing CO2 sequestration projects that require geochemical models to gain an understanding of the processes occurring in the subsurface. The ISGS has collected brine and freshwater samples associated with an enhanced oil recovery project in the Loudon oil field. Geochemical modeling allows us to understand reactions with carbonate and silicate minerals in the reservoir, and the effects they have had on brine composition. For the Illinois Basin Decatur project, geochemical models should allow predictions of the reactions that will take place before CO2 injection begins. ?? 2009 Elsevier Ltd. All rights reserved.

Lagrangian simulation of mixing and reactions in complex geochemical systems

NASA Astrophysics Data System (ADS)

Engdahl, Nicholas B.; Benson, David A.; Bolster, Diogo

2017-04-01

Simulations of detailed geochemical systems have traditionally been restricted to Eulerian reactive transport algorithms. This note introduces a Lagrangian method for modeling multicomponent reaction systems. The approach uses standard random walk-based methods for the particle motion steps but allows the particles to interact with each other by exchanging mass of their various chemical species. The colocation density of each particle pair is used to calculate the mass transfer rate, which creates a local disequilibrium that is then relaxed back toward equilibrium using the reaction engine PhreeqcRM. The mass exchange is the only step where the particles interact and the remaining transport and reaction steps are entirely independent for each particle. Several validation examples are presented, which reproduce well-known analytical solutions. These are followed by two demonstration examples of a competitive decay chain and an acid-mine drainage system. The source code, entitled Complex Reaction on Particles (CRP), and files needed to run these examples are hosted openly on GitHub (https://github.com/nbengdahl/CRP), so as to enable interested readers to readily apply this approach with minimal modifications.

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

PubMed

Yang, Changbing; Samper, Javier; Molinero, Jorge; Bonilla, Mercedes

2007-08-15

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.

Prediction of the effects of soil-based countermeasures on soil solution chemistry of soils contaminated with radiocesium using the hydrogeochemical code PHREEQC.

PubMed

Hormann, Volker; Kirchner, Gerald

2002-04-22

For agriculturally used areas, which are contaminated by the debris from a nuclear accident, the use of chemical amendmends (e.g. potassium chloride and lime) is among the most common soil-based countermeasures. These countermeasures are intended to reduce the plant uptake of radionuclides (mainly 137Cs and 90Sr) by competitive inhibition by chemically similar ions. So far, the impacts of countermeasures on soil solution composition - and thus, their effectiveness - have almost exclusively been established experimentally, since they depend on mineral composition and chemical characteristics of the soil affected. In this study, which focuses on caesium contamination, the well-established code PHREEQC was used as a geochemical model to calculate the changes in the ionic compositions of soil solutions, which result from the application of potassium or ammonium in batch equilibrium experiments. The simple ion exchange model used by PHREEQC was improved by taking into account selective sorption of Cs+, NH4+ and K+ by clay minerals. Calculations were performed with three different initial soil solution compositions, corresponding to particular soil types (loam, sand, peat). For loamy and sandy soils, our calculational results agree well with experimental data reported by Nisbet (Effectiveness of soil-based countermeasures six months and one year after contamination of five diverse soil types with caesium-134 and strontium-90. Contract Report NRPB-M546, National Radiation Protection Board, Chilton, 1995.). For peat, discrepancies were found indicating that for organic soils a reliable set of exchange constants of the relevant cations still has to be determined experimentally. For cesium, however, these discrepancies almost disappeared if selective sites were assumed to be inaccessible. Additionally, results of sensitivity analyses are presented by which the influence of the main soil parameters on Cs+ concentrations in solution after soil treatment has been systematically studied. It is shown that calculating the impacts of soil-based chemical countermeasures on soil solution chemistry using geochemical codes such as PHREEQC offers an attractive alternative to establishing these impacts by often time-consuming and site-specific experiments.

Development of a Model of Geophysical and Geochemical Controls on Abiotic Carbon Cycling on Earth-Like Planets

NASA Astrophysics Data System (ADS)

Neveu, M.; Felton, R.; Domagal-Goldman, S. D.; Desch, S. J.; Arney, G. N.

2017-12-01

About 20 Earth-sized planets (0.6-1.6 Earth masses and radii) have now been discovered beyond our solar system [1]. Although such planets are prime targets in the upcoming search for atmospheric biosignatures, their composition, geology, and climate are essentially unconstrained. Yet, developing an understanding of how these factors influence planetary evolution through time and space is essential to establishing abiotic backgrounds against which any deviations can provide evidence for biological activity. To this end, we are building coupled geophysical-geochemical models of abiotic carbon cycling on such planets. Our models are controlled by atmospheric factors such as temperature and composition, and compute interior inputs to atmospheric species. They account for crustal weathering, ocean-atmosphere equilibria, and exchange with the deep interior as a function of planet composition and size (and, eventually, age).Planets in other solar systems differ from the Earth not only in their bulk physical properties, but also likely in their bulk chemical composition [2], which influences key parameters such as the vigor of mantle convection and the near-surface redox state. Therefore, simulating how variations in such parameters affect carbon cycling requires us to simulate the above processes from first principles, rather than by using arbitrary parameterizations derived from observations as is often done with models of carbon cycling on Earth [3] or extrapolations thereof [4]. As a first step, we have developed a kinetic model of crustal weathering using the PHREEQC code [5] and kinetic data from [6]. We will present the ability of such a model to replicate Earth's carbon cycle using, for the time being, parameterizations for surface-interior-atmosphere exchange processes such as volcanism (e.g., [7]).[1] exoplanet.eu, 7/28/2017.[2] Young et al. (2014) Astrobiology 14, 603-626.[3] Lerman & Wu (2008) Kinetics of Global Geochemical Cycles. In Kinetics of Water-Rock Interaction (Brantley et al., eds.), Springer, New York.[4] Edson et al. (2012) Astrobiology 12, 562-571.[5] Parkhurst & Appelo (2013) USGS Techniques and Methods 6-A43.[6] Palandri & Kharaka (2008) USGS Report 2004-1068.[7] Kite et al. (2009) ApJ 700, 1732-1749.

Hiereachical Bayesian Model for Combining Geochemical and Geophysical Data for Environmental Applications Software

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Jinsong

2013-05-01

Development of a hierarchical Bayesian model to estimate the spatiotemporal distribution of aqueous geochemical parameters associated with in-situ bioremediation using surface spectral induced polarization (SIP) data and borehole geochemical measurements collected during a bioremediation experiment at a uranium-contaminated site near Rifle, Colorado. The SIP data are first inverted for Cole-Cole parameters including chargeability, time constant, resistivity at the DC frequency and dependence factor, at each pixel of two-dimensional grids using a previously developed stochastic method. Correlations between the inverted Cole-Cole parameters and the wellbore-based groundwater chemistry measurements indicative of key metabolic processes within the aquifer (e.g. ferrous iron, sulfate, uranium)more » were established and used as a basis for petrophysical model development. The developed Bayesian model consists of three levels of statistical sub-models: 1) data model, providing links between geochemical and geophysical attributes, 2) process model, describing the spatial and temporal variability of geochemical properties in the subsurface system, and 3) parameter model, describing prior distributions of various parameters and initial conditions. The unknown parameters are estimated using Markov chain Monte Carlo methods. By combining the temporally distributed geochemical data with the spatially distributed geophysical data, we obtain the spatio-temporal distribution of ferrous iron, sulfate and sulfide, and their associated uncertainity information. The obtained results can be used to assess the efficacy of the bioremediation treatment over space and time and to constrain reactive transport models.« less

A Spatially Constrained Multi-autoencoder Approach for Multivariate Geochemical Anomaly Recognition

NASA Astrophysics Data System (ADS)

Lirong, C.; Qingfeng, G.; Renguang, Z.; Yihui, X.

2017-12-01

Separating and recognizing geochemical anomalies from the geochemical background is one of the key tasks in geochemical exploration. Many methods have been developed, such as calculating the mean ±2 standard deviation, and fractal/multifractal models. In recent years, deep autoencoder, a deep learning approach, have been used for multivariate geochemical anomaly recognition. While being able to deal with the non-normal distributions of geochemical concentrations and the non-linear relationships among them, this self-supervised learning method does not take into account the spatial heterogeneity of geochemical background and the uncertainty induced by the randomly initialized weights of neurons, leading to ineffective recognition of weak anomalies. In this paper, we introduce a spatially constrained multi-autoencoder (SCMA) approach for multivariate geochemical anomaly recognition, which includes two steps: spatial partitioning and anomaly score computation. The first step divides the study area into multiple sub-regions to segregate the geochemical background, by grouping the geochemical samples through K-means clustering, spatial filtering, and spatial constraining rules. In the second step, for each sub-region, a group of autoencoder neural networks are constructed with an identical structure but different initial weights on neurons. Each autoencoder is trained using the geochemical samples within the corresponding sub-region to learn the sub-regional geochemical background. The best autoencoder of a group is chosen as the final model for the corresponding sub-region. The anomaly score at each location can then be calculated as the euclidean distance between the observed concentrations and reconstructed concentrations of geochemical elements.The experiments using the geochemical data and Fe deposits in the southwestern Fujian province of China showed that our SCMA approach greatly improved the recognition of weak anomalies, achieving the AUC of 0.89, compared with the AUC of 0.77 using a single deep autoencoder approach.

Contaminant source identification using semi-supervised machine learning

NASA Astrophysics Data System (ADS)

Vesselinov, Velimir V.; Alexandrov, Boian S.; O'Malley, Daniel

2018-05-01

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. In this paper, we propose a new contaminant source identification approach that performs decomposition of the observation mixtures based on Non-negative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the unknown number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. NMFk is tested on synthetic and real-world site data. The NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).

Contaminant source identification using semi-supervised machine learning

DOE PAGES

Vesselinov, Velimir Valentinov; Alexandrov, Boian S.; O’Malley, Dan

2017-11-08

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may needmore » to be simulated in these models which further complicates the analyses. In this paper, we propose a new contaminant source identification approach that performs decomposition of the observation mixtures based on Non-negative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the unknown number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. NMFk is tested on synthetic and real-world site data. Finally, the NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).« less

Contaminant source identification using semi-supervised machine learning

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vesselinov, Velimir Valentinov; Alexandrov, Boian S.; O’Malley, Dan

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may needmore » to be simulated in these models which further complicates the analyses. In this paper, we propose a new contaminant source identification approach that performs decomposition of the observation mixtures based on Non-negative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the unknown number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. NMFk is tested on synthetic and real-world site data. Finally, the NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).« less

A modified procedure for mixture-model clustering of regional geochemical data

USGS Publications Warehouse

Ellefsen, Karl J.; Smith, David B.; Horton, John D.

2014-01-01

A modified procedure is proposed for mixture-model clustering of regional-scale geochemical data. The key modification is the robust principal component transformation of the isometric log-ratio transforms of the element concentrations. This principal component transformation and the associated dimension reduction are applied before the data are clustered. The principal advantage of this modification is that it significantly improves the stability of the clustering. The principal disadvantage is that it requires subjective selection of the number of clusters and the number of principal components. To evaluate the efficacy of this modified procedure, it is applied to soil geochemical data that comprise 959 samples from the state of Colorado (USA) for which the concentrations of 44 elements are measured. The distributions of element concentrations that are derived from the mixture model and from the field samples are similar, indicating that the mixture model is a suitable representation of the transformed geochemical data. Each cluster and the associated distributions of the element concentrations are related to specific geologic and anthropogenic features. In this way, mixture model clustering facilitates interpretation of the regional geochemical data.

Comparison of U-spatial statistics and C-A fractal models for delineating anomaly patterns of porphyry-type Cu geochemical signatures in the Varzaghan district, NW Iran

NASA Astrophysics Data System (ADS)

Ghezelbash, Reza; Maghsoudi, Abbas

2018-05-01

The delineation of populations of stream sediment geochemical data is a crucial task in regional exploration surveys. In this contribution, uni-element stream sediment geochemical data of Cu, Au, Mo, and Bi have been subjected to two reliable anomaly-background separation methods, namely, the concentration-area (C-A) fractal and the U-spatial statistics methods to separate geochemical anomalies related to porphyry-type Cu mineralization in northwest Iran. The quantitative comparison of the delineated geochemical populations using the modified success-rate curves revealed the superiority of the U-spatial statistics method over the fractal model. Moreover, geochemical maps of investigated elements revealed strongly positive correlations between strong anomalies and Oligocene-Miocene intrusions in the study area. Therefore, follow-up exploration programs should focus on these areas.

Use of the Water, Energy, and Biogeochemical Model (WEBMOD) to Simulate Water Quality at Five U.S. Geological Survey Research Watersheds

NASA Astrophysics Data System (ADS)

Webb, R. M.; Leavesley, G. H.; Shanley, J. B.; Peters, N. E.; Aulenbach, B. T.; Blum, A. E.; Campbell, D. H.; Clow, D. W.; Mast, M. A.; Stallard, R. F.; Larsen, M. C.; Troester, J. W.; Walker, J. F.; White, A. F.

2003-12-01

The Water, Energy, and Biogeochemical Model (WEBMOD) was developed as an aid to compare and contrast basic hydrologic and biogeochemical processes active in the diverse hydroclimatic regions represented by the five U.S. Geological Survey (USGS) Water, Energy, and Biogeochemical Budget (WEBB) sites: Loch Vale, Colorado; Trout Lake, Wisconsin; Sleepers River, Vermont; Panola Mountain, Georgia; and Luquillo Experimental Forest, Puerto Rico. WEBMOD simulates solute concentrations for vegetation canopy, snow pack, impermeable ground, leaf litter, unsaturated and saturated soil zones, riparian zones and streams using selected process modules coupled within the USGS Modular Modeling System (MMS). Source codes for the MMS hydrologic modules include the USGS Precipitation Runoff Modeling System, the National Weather Service Hydro-17 snow model, and TOPMODEL. The hydrologic modules distribute precipitation and temperature to predict evapotranspiration, snow accumulation, snow melt, and streamflow. Streamflow generation mechanisms include infiltration excess, saturated overland flow, preferential lateral flow, and base flow. Input precipitation chemistry, and fluxes and residence times predicted by the hydrologic modules are input into the geochemical module where solute concentrations are computed for a series of discrete well-mixed reservoirs using calls to the geochemical engine PHREEQC. WEBMOD was used to better understand variations in water quality observed at the WEBB sites from October 1991 through September 1997. Initial calibrations were completed by fitting the simulated hydrographs with those measured at the watershed outlets. Model performance was then refined by comparing the predicted export of conservative chemical tracers such as chloride, with those measured at the watershed outlets. The model succeeded in duplicating the temporal variability of net exports of major ions from the watersheds.

User’s guide for GcClust—An R package for clustering of regional geochemical data

USGS Publications Warehouse

Ellefsen, Karl J.; Smith, David B.

2016-04-08

GcClust is a software package developed by the U.S. Geological Survey for statistical clustering of regional geochemical data, and similar data such as regional mineralogical data. Functions within the software package are written in the R statistical programming language. These functions, their documentation, and a copy of the user’s guide are bundled together in R’s unit of sharable code, which is called a “package.” The user’s guide includes step-by-step instructions showing how the functions are used to cluster data and to evaluate the clustering results. These functions are demonstrated in this report using test data, which are included in the package.

Three-dimensional Numerical Models of the Cocos-northern Nazca Slab Gap

NASA Astrophysics Data System (ADS)

Jadamec, M.; Fischer, K. M.

2012-12-01

In contrast to anisotropy beneath the middle of oceanic plates, seismic observations in subduction zones often indicate mantle flow patterns that are not easily explained by simple coupling of the subducting and overriding plates to the mantle. For example, in the Costa Rica-Nicaragua subduction zone local S shear wave splitting measurements combined with geochemical data indicate trench parallel flow in the mantle wedge with flow rates of 6.3-19 cm/yr, which is on order of or may be up to twice the subducting plate velocity. We construct geographically referenced high-resolution three-dimensional (3D) geodynamic models of the Cocos-northern Nazca subduction system to investigate what is driving the northwest directed, and apparently rapid, trench-parallel flow in the mantle wedge beneath Costa Rica-Nicaragua. We use the SlabGenerator code to construct a 3D plate configuration that is used as input to the community mantle convection code, CitcomCU. Models are run on over 400 CPUs on XSEDE, with a mesh resolution of up to 3 km at the plate boundary. Seismicity and seismic tomography delineate the shape and depth of the Cocos and northern Nazca slabs. The subducting plate thermal structure is based on a plate cooling model and ages from the seafloor age grid. Overriding plate thickness is constrained by the ages from the sea floor age grid where available and the depth to the lithosphere-asthenosphere boundary from the greatest negative gradient in absolute shear wave velocity. The geodynamic models test the relative controls of the change in the dip of the Cocos plate and the slab gap between the Cocos and northern Nazca plates in driving the mantle flow beneath Central America. The models also investigate the effect of a non-Newtonian rheology in dynamically generating a low viscosity mantle wedge and how this controls mantle flow rates. To what extent the Cocos-northern Nazca slab gap channelizes mantle flow between Central and South America has direct application to geochemical and geologic studies of the region. In addition, 3D geodynamic models of this kind can further test the hypothesis of rapid mantle flow in subduction zones as a global process and the non-Newtonian rheology as a mechanism for decoupling the mantle from lithospheric plate motion.

Controls on the geochemical evolution of Prairie Pothole Region lakes and wetlands over decadal time scales

USGS Publications Warehouse

Goldhaber, Martin B.; Mills, Christopher T.; Mushet, David M.; McCleskey, R. Blaine; Rover, Jennifer

2016-01-01

One hundred sixty-seven Prairie Pothole lakes, ponds and wetlands (largely lakes) previously analyzed chemically during the late 1960’s and early to mid-1970’s were resampled and reanalyzed in 2011–2012. The two sampling periods differed climatically. The earlier sampling took place during normal to slightly dry conditions, whereas the latter occurred during and immediately following exceptionally wet conditions. As reported previously in Mushet et al. (2015), the dominant effect was expansion of the area of these lakes and dilution of their major ions. However, within that context, there were significant differences in the evolutionary pathways of major ions. To establish these pathways, we employed the inverse modeling computer code NetpathXL. This code takes the initial and final lake composition and, using mass balance constrained by the composition of diluting waters, and input and output of phases, calculates plausible geochemical evolution pathways. Despite the fact that in most cases major ions decreased, a subset of the lakes had an increase in SO42−. This distinction is significant because SO42− is the dominant anion in a majority of Prairie Pothole Region wetlands and lakes. For lakes with decreasing SO42−, the proportion of original lake water required for mass balance was subordinate to rainwater and/or overland flow. In contrast, lakes with increasing SO42− between the two sampling episodes tended to be dominated by original lake water. This suite of lakes tended to be smaller and have lower initial SO42−concentrations such that inputs of sulfur from dissolution of the minerals gypsum or pyrite had a significant impact on the final sulfur concentration given the lower dilution factors. Thus, our study provides context for how Prairie Pothole Region water bodies evolve geochemically as climate changes. Because wetland geochemistry in turn controls the ecology of these water bodies, this research contributes to the prediction of the impact of climate change on this important complex of ecosystems.

Impact of the volume of gaseous phase in closed reactors on ANC results and modelling

NASA Astrophysics Data System (ADS)

Drapeau, Clémentine; Delolme, Cécile; Lassabatere, Laurent; Blanc, Denise

2016-04-01

The understanding of the geochemical behavior of polluted solid materials is often challenging and requires huge expenses of time and money. Nevertheless, given the increasing amounts of polluted solid materials and related risks for the environment, it is more and more crucial to understand the leaching of majors and trace metals elements from these matrices. In the designs of methods to quantify pollutant solubilization, the combination of experimental procedures with modeling approaches has recently gained attention. Among usual methods, some rely on the association of ANC and geochemical modeling. ANC experiments - Acid Neutralization Capacity - consists in adding known quantities of acid or base to a mixture of water and contaminated solid materials at a given liquid / solid ratio in closed reactors. Reactors are agitated for 48h and then pH, conductivity, redox potential, carbon, majors and heavy metal solubilized are quantified. However, in most cases, the amounts of matrix and water do not reach the total volume of reactors, leaving some space for air (gaseous phase). Despite this fact, no clear indication is given in standard procedures about the effect of this gaseous phase. Even worse, the gaseous phase is never accounted for when exploiting or modeling ANC data. The gaseous phase may exchange CO2 with the solution, which may, in turn, impact both pH and element release. This study lies within the most general framework for the use of geochemical modeling for the prediction of ANC results for the case of pure phases to real phase assemblages. In this study, we focus on the effect of the gaseous phase on ANC experiments on different mineral phases through geochemical modeling. To do so, we use PHREEQC code to model the evolution of pH and element release (including majors and heavy metals) when several matrices are put in contact with acid or base. We model the following scenarios for the gaseous phase: no gas, contact with the atmosphere (open system) and real reactors conditions (semi-closed systems). The solid phases tested are pure phases (calcite, sulfides, etc.) and synthetic assemblages mimicking complex polluted matrices. The modeling clearly shows that the systems are sensitive to the opening to the atmosphere. If the open system and the system with no gas are entirely different, "real" reactors also differ significantly from the other systems. Apparently, the presence of the gaseous phase in reactors greatly impacts pH and element release. This parameter should be accounted for in ANC experimental procedures and modeling. In addition to this numerical study, experimental results, previously obtained for urban polluted sediments, are analyzed in lights of the findings of the numerical study. This step allows us to strengthen conclusions and to pinpoint at the necessity to account for the gaseous phase when performing and modeling ANC experiments.

Uncertainty quantification and experimental design based on unsupervised machine learning identification of contaminant sources and groundwater types using hydrogeochemical data

NASA Astrophysics Data System (ADS)

Vesselinov, V. V.

2017-12-01

Identification of the original groundwater types present in geochemical mixtures observed in an aquifer is a challenging but very important task. Frequently, some of the groundwater types are related to different infiltration and/or contamination sources associated with various geochemical signatures and origins. The characterization of groundwater mixing processes typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical species. Numerous geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. As a result, these types of model analyses are typically extremely challenging. Here, we demonstrate a new contaminant source identification approach that performs decomposition of the observation mixtures based on Nonnegative Matrix Factorization (NMF) method for Blind Source Separation (BSS), coupled with a custom semi-supervised clustering algorithm. Our methodology, called NMFk, is capable of identifying (a) the number of groundwater types and (b) the original geochemical concentration of the contaminant sources from measured geochemical mixtures with unknown mixing ratios without any additional site information. We also demonstrate how NMFk can be extended to perform uncertainty quantification and experimental design related to real-world site characterization. The NMFk algorithm works with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios). The NMFk algorithm has been extensively tested on synthetic datasets; NMFk analyses have been actively performed on real-world data collected at the Los Alamos National Laboratory (LANL) groundwater sites related to Chromium and RDX contamination.

Motivation, description, and summary status of geomechanical andgeochemical modeling studies in Task D of the InternationalDECOVALEX-THMC Project

DOE Office of Scientific and Technical Information (OSTI.GOV)

Birkholzer, J.T.; Barr, D.; Rutqvist, J.

2005-11-15

The DECOVALEX project is an international cooperativeproject initiated by SKI, the Swedish Nuclear Power Inspectorate, withparticipation of about 10 international organizations. The general goalof this project is to encourage multidisciplinary interactive andcooperative research on modelling coupledthermo-hydro-mechanical-chemical (THMC) processes in geologic formationsin support of the performance assessment for underground storage ofradioactive waste. One of the research tasks, initiated in 2004 by theU.S. Department of Energy (DOE), addresses the long-term impact ofgeomechanical and geochemical processes on the flow conditions near wasteemplacement tunnels. Within this task, four international research teamsconduct predictive analysis of the coupled processes in two genericrepositories, using multiple approaches andmore » different computer codes.Below, we give an overview of the research task and report its currentstatus.« less

MODELING MONOMETHYLMERCURY AND TRIBUTYLTIN SPECIATION WITH EPA'S GEOCHEMICAL SPECIATION MODEL MINTEQA2

EPA Science Inventory

Given the complexity of the various, simultaneous (and competing) equilibrium reactions governing the speciation of ionic species in aquatic systems, EPA has developed and distributed the geochemical speciation model MINTEQA2 (Brown and Allison, 1987, Allison et al., 1991; Hydrog...

Programming PHREEQC calculations with C++ and Python a comparative study

USGS Publications Warehouse

Charlton, Scott R.; Parkhurst, David L.; Muller, Mike

2011-01-01

The new IPhreeqc module provides an application programming interface (API) to facilitate coupling of other codes with the U.S. Geological Survey geochemical model PHREEQC. Traditionally, loose coupling of PHREEQC with other applications required methods to create PHREEQC input files, start external PHREEQC processes, and process PHREEQC output files. IPhreeqc eliminates most of this effort by providing direct access to PHREEQC capabilities through a component object model (COM), a library, or a dynamically linked library (DLL). Input and calculations can be specified through internally programmed strings, and all data exchange between an application and the module can occur in computer memory. This study compares simulations programmed in C++ and Python that are tightly coupled with IPhreeqc modules to the traditional simulations that are loosely coupled to PHREEQC. The study compares performance, quantifies effort, and evaluates lines of code and the complexity of the design. The comparisons show that IPhreeqc offers a more powerful and simpler approach for incorporating PHREEQC calculations into transport models and other applications that need to perform PHREEQC calculations. The IPhreeqc module facilitates the design of coupled applications and significantly reduces run times. Even a moderate knowledge of one of the supported programming languages allows more efficient use of PHREEQC than the traditional loosely coupled approach.

Field-Integrated Studies of Long-Term Sustainability of Chromium Bioreduction at Hanford 100H Site

DOE Office of Scientific and Technical Information (OSTI.GOV)

Long, Philip E.

2006-06-01

The objectives of the project are to investigate coupled hydraulic, geochemical, and microbial conditions, and to determine the critical biogeochemical parameters necessary to maximize the extent of Cr(VI) bioreduction and minimize Cr(III) reoxidation in groundwater. Specific goals of the project are as follows: (1) Field testing and monitoring of Cr(VI) bioreduction in ground water and its transformation into insoluble species of Cr(III) at the Hanford 100H site, to develop the optimal strategy of water sampling for chemical, microbial, stable isotope analyses, and noninvasive geophysical monitoring; (2) Bench-scale flow and transport investigations using columns of undisturbed sediments to obtain diffusion andmore » kinetic parameters needed for the development of a numerical model, predictions of Cr(VI) bioreduction, and potential of Cr(III) reoxidation; and (3) Development of a multiphase, multi-component 3D reactive transport model and a code, TOUGHREACT-BIO, to predict coupled biogeochemical-hydrological processes associated with bioremediation, and to calibrate and validate the developed code based on the results of bench-scale and field-scale Cr(VI) biostimulation experiments in ground water at the Hanford Site.« less

GEOCHEMICAL MODELING OF ARSENIC SPECIATION AND MOBILIZATION: IMPLICATIONS FOR BIOREMEDIATION

EPA Science Inventory

Geochemical modeling techniques were used to examine the biogeochemical linkages between Fe, S, and As in shallow alluvial aquifers. We modeled: 1) the adsorption and desorption of As on the surface of hydrous ferric oxides (HFO’s) in stream beds under aerobic conditions; 2) red...

Adaptive Multiscale Modeling of Geochemical Impacts on Fracture Evolution

NASA Astrophysics Data System (ADS)

Molins, S.; Trebotich, D.; Steefel, C. I.; Deng, H.

2016-12-01

Understanding fracture evolution is essential for many subsurface energy applications, including subsurface storage, shale gas production, fracking, CO2 sequestration, and geothermal energy extraction. Geochemical processes in particular play a significant role in the evolution of fractures through dissolution-driven widening, fines migration, and/or fracture sealing due to precipitation. One obstacle to understanding and exploiting geochemical fracture evolution is that it is a multiscale process. However, current geochemical modeling of fractures cannot capture this multi-scale nature of geochemical and mechanical impacts on fracture evolution, and is limited to either a continuum or pore-scale representation. Conventional continuum-scale models treat fractures as preferential flow paths, with their permeability evolving as a function (often, a cubic law) of the fracture aperture. This approach has the limitation that it oversimplifies flow within the fracture in its omission of pore scale effects while also assuming well-mixed conditions. More recently, pore-scale models along with advanced characterization techniques have allowed for accurate simulations of flow and reactive transport within the pore space (Molins et al., 2014, 2015). However, these models, even with high performance computing, are currently limited in their ability to treat tractable domain sizes (Steefel et al., 2013). Thus, there is a critical need to develop an adaptive modeling capability that can account for separate properties and processes, emergent and otherwise, in the fracture and the rock matrix at different spatial scales. Here we present an adaptive modeling capability that treats geochemical impacts on fracture evolution within a single multiscale framework. Model development makes use of the high performance simulation capability, Chombo-Crunch, leveraged by high resolution characterization and experiments. The modeling framework is based on the adaptive capability in Chombo which not only enables mesh refinement, but also refinement of the model-pore scale or continuum Darcy scale-in a dynamic way such that the appropriate model is used only when and where it is needed. Explicit flux matching provides coupling betwen the scales.

Uranium plume persistence impacted by hydrologic and geochemical heterogeneity in the groundwater and river water interaction zone of Hanford site

NASA Astrophysics Data System (ADS)

Chen, X.; Zachara, J. M.; Vermeul, V. R.; Freshley, M.; Hammond, G. E.

2015-12-01

The behavior of a persistent uranium plume in an extended groundwater- river water (GW-SW) interaction zone at the DOE Hanford site is dominantly controlled by river stage fluctuations in the adjacent Columbia River. The plume behavior is further complicated by substantial heterogeneity in physical and geochemical properties of the host aquifer sediments. Multi-scale field and laboratory experiments and reactive transport modeling were integrated to understand the complex plume behavior influenced by highly variable hydrologic and geochemical conditions in time and space. In this presentation we (1) describe multiple data sets from field-scale uranium adsorption and desorption experiments performed at our experimental well-field, (2) develop a reactive transport model that incorporates hydrologic and geochemical heterogeneities characterized from multi-scale and multi-type datasets and a surface complexation reaction network based on laboratory studies, and (3) compare the modeling and observation results to provide insights on how to refine the conceptual model and reduce prediction uncertainties. The experimental results revealed significant spatial variability in uranium adsorption/desorption behavior, while modeling demonstrated that ambient hydrologic and geochemical conditions and heterogeneities in sediment physical and chemical properties both contributed to complex plume behavior and its persistence. Our analysis provides important insights into the characterization, understanding, modeling, and remediation of groundwater contaminant plumes influenced by surface water and groundwater interactions.

VIRUS TRANSPORT IN PHYSICALLY AND GEOCHEMICALLY HETEROGENEOUS SUBSURFACE POROUS MEDIA. (R826179)

EPA Science Inventory

A two-dimensional model for virus transport in physically and geochemically heterogeneous subsurface porous media is presented. The model involves solution of the advection–dispersion equation, which additionally considers virus inactivation in the solution, as well as ...

INEEL Subregional Conceptual Model Report Volume 2: Summary of Existing Knowledge of Geochemical Influences on the Fate and Transport of Contaminants in the Subsurface at the INEEL

DOE Office of Scientific and Technical Information (OSTI.GOV)

Paul L. Wichlacz; Robert C. Starr; Brennon Orr

2003-09-01

This document summarizes previous descriptions of geochemical system conceptual models for the vadose zone and groundwater zone (aquifer) beneath the Idaho National Engineering and Environmental Laboratory (INEEL). The primary focus is on groundwater because contaminants derived from wastes disposed at INEEL are present in groundwater, groundwater provides a pathway for potential migration to receptors, and because geochemical characteristics in and processes in the aquifer can substantially affect the movement, attenuation, and toxicity of contaminants. The secondary emphasis is perched water bodies in the vadose zone. Perched water eventually reaches the regional groundwater system, and thus processes that affect contaminants inmore » the perched water bodies are important relative to the migration of contaminants into groundwater. Similarly, processes that affect solutes during transport from nearsurface disposal facilities downward through the vadose zone to the aquifer are relevant. Sediments in the vadose zone can affect both water and solute transport by restricting the downward migration of water sufficiently that a perched water body forms, and by retarding solute migration via ion exchange. Geochemical conceptual models have been prepared by a variety of researchers for different purposes. They have been published in documents prepared by INEEL contractors, the United States Geological Survey (USGS), academic researchers, and others. The documents themselves are INEEL and USGS reports, and articles in technical journals. The documents reviewed were selected from citation lists generated by searching the INEEL Technical Library, the INEEL Environmental Restoration Optical Imaging System, and the ISI Web of Science databases. The citation lists were generated using the keywords ground water, groundwater, chemistry, geochemistry, contaminant, INEL, INEEL, and Idaho. In addition, a list of USGS documents that pertain to the INEEL was obtained and manually searched. The documents that appeared to be the most pertinent were selected from further review. These documents are tabulated in the citation list. This report summarizes existing geochemical conceptual models, but does not attempt to generate a new conceptual model or select the ''right'' model. This document is organized as follows. Geochemical models are described in general in Section 2. Geochemical processes that control the transport and fate of contaminants introduced into groundwater are described in Section 3. The natural geochemistry of the Eastern Snake River Plain Aquifer (SRPA) is described in Section 4. The effect of waste disposal on the INEEL subsurface is described in Section 5. The geochemical behavior of the major contaminants is described in Section 6. Section 7 describes the site-specific geochemical models developed for various INEEL facilities.« less

A Watched Ocean World Never Boils: Inspecting the Geochemical Impact on Ocean Worlds from Their Thermal Evolution

NASA Astrophysics Data System (ADS)

Spiers, E. M.; Schmidt, B. E.

2018-05-01

I aim to acquire better understanding of coupled thermal evolution and geochemical fluxes of an ocean world through a box model. A box model divides the system into plainer elements with realistically-solvable, dynamic equations.

Geochemical Characterization Using Geophysical Data and Markov Chain Monte Carlo Methods

NASA Astrophysics Data System (ADS)

Chen, J.; Hubbard, S.; Rubin, Y.; Murray, C.; Roden, E.; Majer, E.

2002-12-01

Although the spatial distribution of geochemical parameters is extremely important for many subsurface remediation approaches, traditional characterization of those parameters is invasive and laborious, and thus is rarely performed sufficiently to describe natural hydrogeological variability at the field-scale. This study is an effort to jointly use multiple sources of information, including noninvasive geophysical data, for geochemical characterization of the saturated and anaerobic portion of the DOE South Oyster Bacterial Transport Site in Virginia. Our data set includes hydrogeological and geochemical measurements from five boreholes and ground-penetrating radar (GPR) and seismic tomographic data along two profiles that traverse the boreholes. The primary geochemical parameters are the concentrations of extractable ferrous iron Fe(II) and ferric iron Fe(III). Since iron-reducing bacteria can reduce Fe(III) to Fe(II) under certain conditions, information about the spatial distributions of Fe(II) and Fe(III) may indicate both where microbial iron reduction has occurred and in which zone it is likely to occur in the future. In addition, as geochemical heterogeneity influences bacterial transport and activity, estimates of the geochemical parameters provide important input to numerical flow and contaminant transport models geared toward bioremediation. Motivated by our previous research, which demonstrated that crosshole geophysical data could be very useful for estimating hydrogeological parameters, we hypothesize in this study that geochemical and geophysical parameters may be linked through their mutual dependence on hydrogeological parameters such as lithofacies. We attempt to estimate geochemical parameters using both hydrogeological and geophysical measurements in a Bayesian framework. Within the two-dimensional study domain (12m x 6m vertical cross section divided into 0.25m x 0.25m pixels), geochemical and hydrogeological parameters were considered as data if they were available from direct measurements or as variables otherwise. To estimate the geochemical parameters, we first assigned a prior model for each variable and a likelihood model for each type of data, which together define posterior probability distributions for each variable on the domain. Since the posterior probability distribution may involve hundreds of variables, we used a Markov Chain Monte Carlo (MCMC) method to explore each variable by generating and subsequently evaluating hundreds of realizations. Results from this case study showed that although geophysical attributes are not necessarily directly related to geochemical parameters, geophysical data could be very useful for providing accurate and high-resolution information about geochemical parameter distribution through their joint and indirect connections with hydrogeological properties such as lithofacies. This case study also demonstrated that MCMC methods were particularly useful for geochemical parameter estimation using geophysical data because they allow incorporation into the procedure of spatial correlation information, measurement errors, and cross correlations among different types of parameters.

Delineation of geochemical anomalies based on stream sediment data utilizing fractal modeling and staged factor analysis

NASA Astrophysics Data System (ADS)

Afzal, Peyman; Mirzaei, Misagh; Yousefi, Mahyar; Adib, Ahmad; Khalajmasoumi, Masoumeh; Zarifi, Afshar Zia; Foster, Patrick; Yasrebi, Amir Bijan

2016-07-01

Recognition of significant geochemical signatures and separation of geochemical anomalies from background are critical issues in interpretation of stream sediment data to define exploration targets. In this paper, we used staged factor analysis in conjunction with the concentration-number (C-N) fractal model to generate exploration targets for prospecting Cr and Fe mineralization in Balvard area, SE Iran. The results show coexistence of derived multi-element geochemical signatures of the deposit-type sought and ultramafic-mafic rocks in the NE and northern parts of the study area indicating significant chromite and iron ore prospects. In this regard, application of staged factor analysis and fractal modeling resulted in recognition of significant multi-element signatures that have a high spatial association with host lithological units of the deposit-type sought, and therefore, the generated targets are reliable for further prospecting of the deposit in the study area.

Subsurface Reactive Transport Modelling of the Lateritic Ni mineralization in New Caledonia: A coupled Thermo-Hydro-Geochemical Approach

NASA Astrophysics Data System (ADS)

Myagkiy, Andrey; Golfier, Fabrice; Truche, Laurent; Cathelineau, Michel

2017-04-01

This research proposes a subsurface reactive geochemical transport modelling of the development of a nickel laterite profile in New Caledonia over the past few million years. Such a regolith formation from ultramafic bedrock was not yet modelled and gives new profound insights into the Ni vertical mobility, its retention processes in a soil profile and relative enrichment, that are still poorly studied. The downward progression of the lateritization front is allowed by the leaching of the soluble elements (Si, Mg and Ni) through drainage system, represented by porous column of peridotite. Particular emphasis is placed on the detailed understanding of Ni redistribution as a function of time and depth triggered by Ni-bearing silicate precipitation (i.e. garnierite) and by sorption or recrystallization process with goethite. Current work consists of the following models: i) 1-D calculations that are done at 25oC with the code PHREEQC associated with the llnl thermodynamic database and ii) 2-D model that handles coupled thermo-hydro-chemical processes and is calculated on the interface Comsol-Phreeqc (iCP, Nardi et al., 2014). The impact of i) fluid flow in fractures and ii) recharge rate along with iii) hydraulic and iv) geothermal gradients are considered here. While the first model gives profound insights into the vertical mobility of metals upon the formation of laterite (Myagkiy et al, submitted), the latter one additionally allows to describe heterogeneities of mineralizing distributions due to the influence of preferential pathways (fractures), convective flows and lateral transfers. Our long-term 1-D simulations (10 Ma) clearly demonstrate that the Ni enrichment and thickening of iron-rich zone are governed by the vertical progression of the pH front. At the same time 2-D modelling shows reactivation of Ni from oxide zone and it subsequent redistribution and concentration in saprolite. Such a model appears to be of importance in attempt of explanation Ni mineralization processes, revealing the main keys to understanding the trace elements mobility in ultramafic environment. Myagkiy A, Truche L, Cathelineau M, Golfier F. "Revealing the conditions of Ni mineralization in laterite profile of New Caledonia: insights from reactive geochemical transport modelling" Preprint submitted to Chemical Geology (September 28, 2016). Nardi A, Idiart A, Trinchero P, de Vries LM, and Molinero J. "Interface COMSOL-PHREEQC (iCP), an efficient numerical framework for the solution of coupled multiphysics and geochemistry."Computers & Geosciences 69 (2014): 10-21.

Computer finds ore

NASA Astrophysics Data System (ADS)

Bell, Peter M.

Artificial intelligence techniques are being used for the first time to evaluate geophysical, geochemical, and geologic data and theory in order to locate ore deposits. After several years of development, an intelligent computer code has been formulated and applied to the Mount Tolman area in Washington state. In a project funded by the United States Geological Survey and the National Science Foundation a set of computer programs, under the general title Prospector, was used successfully to locate a previously unknown ore-grade porphyry molybdenum deposit in the vicinity of Mount Tolman (Science, Sept. 3, 1982).The general area of the deposit had been known to contain exposures of porphyry mineralization. Between 1964 and 1978, exploration surveys had been run by the Bear Creek Mining Company, and later exploration was done in the area by the Amax Corporation. Some of the geophysical data and geochemical and other prospecting surveys were incorporated into the programs, and mine exploration specialists contributed to a set of rules for Prospector. The rules were encoded as ‘inference networks’ to form the ‘expert system’ on which the artificial intelligence codes were based. The molybdenum ore deposit discovered by the test is large, located subsurface, and has an areal extent of more than 18 km2.

Surficial geologic maps along the riparian zone of the Animas River and its headwater tributaries, Silverton to Durango, Colorado, with upper Animas River watershed gradient profiles

USGS Publications Warehouse

Blair, R.W.; Yager, D.B.; Church, S.E.

2002-01-01

This product consists of Adobe Acrobat .PDF format documents for 10 surficial geologic strip maps along the Animas River watershed from its major headwater tributaries, south to Durango, Colorado. The Animas River originates in the San Juan Mountains north of the historic mining town of Silverton, Colorado. The surficial geologic maps identify surficial deposits, such as flood-plain and terrace gravels, alluvial fans, glacial till, talus, colluvium, landslides, and bogs. Sixteen primary units were mapped that included human-related deposits and structures, eight alluvial, four colluvial, one glacial, travertine deposits, and undifferentiated bedrock. Each of the surficial geologic strip maps has .PDF links to surficial geology photographs, which enable the user to take a virtual tour of these deposits. Geochemical data collected from mapped surficial deposits that pre- and postdate mining activity have aided in determining the geochemical baseline in the watershed. Several photographs with their corresponding geochemical baseline profiles are accessible through .PDF links from several of the maps. A single coverage for all surficial deposits mapped is included as an ArcInfo shape file as an Arc Export format .e00 file. A gradient map for major headwater tributary streams to the Animas River is also included. The gradient map has stream segments that are color-coded based on relative variations in slope and .PDF format links to each stream gradient profile. Stream gradients were derived from U.S. Geological Survey 10-m digital elevation model data. This project was accomplished in support of the U.S. Geological Survey's Abandoned Mine Lands Initiative in the San Juan Mountains, Colorado.

Geochemical Reaction Mechanism Discovery from Molecular Simulation

DOE PAGES

Stack, Andrew G.; Kent, Paul R. C.

2014-11-10

Methods to explore reactions using computer simulation are becoming increasingly quantitative, versatile, and robust. In this review, a rationale for how molecular simulation can help build better geochemical kinetics models is first given. We summarize some common methods that geochemists use to simulate reaction mechanisms, specifically classical molecular dynamics and quantum chemical methods and discuss their strengths and weaknesses. Useful tools such as umbrella sampling and metadynamics that enable one to explore reactions are discussed. Several case studies wherein geochemists have used these tools to understand reaction mechanisms are presented, including water exchange and sorption on aqueous species and mineralmore » surfaces, surface charging, crystal growth and dissolution, and electron transfer. The impact that molecular simulation has had on our understanding of geochemical reactivity are highlighted in each case. In the future, it is anticipated that molecular simulation of geochemical reaction mechanisms will become more commonplace as a tool to validate and interpret experimental data, and provide a check on the plausibility of geochemical kinetic models.« less

Simulating Donnan equilibria based on the Nernst-Planck equation

NASA Astrophysics Data System (ADS)

Gimmi, Thomas; Alt-Epping, Peter

2018-07-01

Understanding ion transport through clays and clay membranes is important for many geochemical and environmental applications. Ion transport is affected by electrostatic forces exerted by charged clay surfaces. Anions are partly excluded from pore water near these surfaces, whereas cations are enriched. Such effects can be modeled by the Donnan approach. Here we introduce a new, comparatively simple way to represent Donnan equilibria in transport simulations. We include charged surfaces as immobile ions in the balance equation and calculate coupled transport of all components, including the immobile charges, with the Nernst-Planck equation. This results in an additional diffusion potential that influences ion transport, leading to Donnan ion distributions while maintaining local charge balance. The validity of our new approach was demonstrated by comparing Nernst-Planck simulations using the reactive transport code Flotran with analytical solutions available for simple Donnan systems. Attention has to be paid to the numerical evaluation of the electrochemical migration term in the Nernst-Planck equation to obtain correct results for asymmetric electrolytes. Sensitivity simulations demonstrate the influence of various Donnan model parameters on simulated anion accessible porosities. It is furthermore shown that the salt diffusion coefficient in a Donnan pore depends on local concentrations, in contrast to the aqueous salt diffusion coefficient. Our approach can be easily implemented into other transport codes. It is versatile and facilitates, for instance, assessing the implications of different activity models for the Donnan porosity.

Integrated Disposal Facility FY 2012 Glass Testing Summary Report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pierce, Eric M.; Kerisit, Sebastien N.; Krogstad, Eirik J.

2013-03-29

PNNL is conducting work to provide the technical basis for estimating radionuclide release from the engineered portion of the disposal facility for Hanford immobilized low-activity waste (ILAW). Before the ILAW can be disposed, DOE must conduct a performance assessment (PA) for the Integrated Disposal Facility (IDF) that describes the long-term impacts of the disposal facility on public health and environmental resources. As part of the ILAW glass testing program, PNNL is implementing a strategy, consisting of experimentation and modeling, to provide the technical basis for estimating radionuclide release from the glass waste form in support of future IDF PAs. Keymore » activities in FY12 include upgrading the STOMP/eSTOMP codes to do near-field modeling, geochemical modeling of PCT tests to determine the reaction network to be used in the STOMP codes, conducting PUF tests on selected glasses to simulate and accelerate glass weathering, developing a Monte Carlo simulation tool to predict the characteristics of the weathered glass reaction layer as a function of glass composition, and characterizing glasses and soil samples exhumed from an 8-year lysimeter test. The purpose of this report is to summarize the progress made in fiscal year (FY) 2012 and the first quarter of FY 2013 toward implementing the strategy with the goal of developing an understanding of the long-term corrosion behavior of LAW glasses.« less

A review of single-sample-based models and other approaches for radiocarbon dating of dissolved inorganic carbon in groundwater

USGS Publications Warehouse

Han, L. F; Plummer, Niel

2016-01-01

Numerous methods have been proposed to estimate the pre-nuclear-detonation 14C content of dissolved inorganic carbon (DIC) recharged to groundwater that has been corrected/adjusted for geochemical processes in the absence of radioactive decay (14C0) - a quantity that is essential for estimation of radiocarbon age of DIC in groundwater. The models/approaches most commonly used are grouped as follows: (1) single-sample-based models, (2) a statistical approach based on the observed (curved) relationship between 14C and δ13C data for the aquifer, and (3) the geochemical mass-balance approach that constructs adjustment models accounting for all the geochemical reactions known to occur along a groundwater flow path. This review discusses first the geochemical processes behind each of the single-sample-based models, followed by discussions of the statistical approach and the geochemical mass-balance approach. Finally, the applications, advantages and limitations of the three groups of models/approaches are discussed.The single-sample-based models constitute the prevailing use of 14C data in hydrogeology and hydrological studies. This is in part because the models are applied to an individual water sample to estimate the 14C age, therefore the measurement data are easily available. These models have been shown to provide realistic radiocarbon ages in many studies. However, they usually are limited to simple carbonate aquifers and selection of model may have significant effects on 14C0 often resulting in a wide range of estimates of 14C ages.Of the single-sample-based models, four are recommended for the estimation of 14C0 of DIC in groundwater: Pearson's model, (Ingerson and Pearson, 1964; Pearson and White, 1967), Han & Plummer's model (Han and Plummer, 2013), the IAEA model (Gonfiantini, 1972; Salem et al., 1980), and Oeschger's model (Geyh, 2000). These four models include all processes considered in single-sample-based models, and can be used in different ranges of 13C values.In contrast to the single-sample-based models, the extended Gonfiantini & Zuppi model (Gonfiantini and Zuppi, 2003; Han et al., 2014) is a statistical approach. This approach can be used to estimate 14C ages when a curved relationship between the 14C and 13C values of the DIC data is observed. In addition to estimation of groundwater ages, the relationship between 14C and δ13C data can be used to interpret hydrogeological characteristics of the aquifer, e.g. estimating apparent rates of geochemical reactions and revealing the complexity of the geochemical environment, and identify samples that are not affected by the same set of reactions/processes as the rest of the dataset. The investigated water samples may have a wide range of ages, and for waters with very low values of 14C, the model based on statistics may give more reliable age estimates than those obtained from single-sample-based models. In the extended Gonfiantini & Zuppi model, a representative system-wide value of the initial 14C content is derived from the 14C and δ13C data of DIC and can differ from that used in single-sample-based models. Therefore, the extended Gonfiantini & Zuppi model usually avoids the effect of modern water components which might retain ‘bomb’ pulse signatures.The geochemical mass-balance approach constructs an adjustment model that accounts for all the geochemical reactions known to occur along an aquifer flow path (Plummer et al., 1983; Wigley et al., 1978; Plummer et al., 1994; Plummer and Glynn, 2013), and includes, in addition to DIC, dissolved organic carbon (DOC) and methane (CH4). If sufficient chemical, mineralogical and isotopic data are available, the geochemical mass-balance method can yield the most accurate estimates of the adjusted radiocarbon age. The main limitation of this approach is that complete information is necessary on chemical, mineralogical and isotopic data and these data are often limited.Failure to recognize the limitations and underlying assumptions on which the various models and approaches are based can result in a wide range of estimates of 14C0 and limit the usefulness of radiocarbon as a dating tool for groundwater. In each of the three generalized approaches (single-sample-based models, statistical approach, and geochemical mass-balance approach), successful application depends on scrutiny of the isotopic (14C and 13C) and chemical data to conceptualize the reactions and processes that affect the 14C content of DIC in aquifers. The recently developed graphical analysis method is shown to aid in determining which approach is most appropriate for the isotopic and chemical data from a groundwater system.

Geochemical Models of Water-Quality Changes During Aquifer Storage Recovery (ASR) Cycle Tests, Phase 1: Geochemical Models Using Existing Data

DTIC Science & Technology

2006-09-01

Richardson, in review). Figure 1 shows the lithostratigraphic setting for Eocene through Miocene strata, and the occurrence of hydrostratigraphic units of...basal Haw- thorn unit lies unconformably on lithologies informally called “ Eocene limestones,” which consist of Suwannee Limestone, Ocala Limestone

Performance of a Zerovalent Iron Reactive Barrier for the Treatment of Arsenic in Groundwater: Part 2. Geochemical Modeling and Solid Phase Studies

EPA Science Inventory

Arsenic uptake processes were evaluated in a zerovalent iron reactive barrier installed at a lead smelting facility using geochemical modeling, solid-phase analysis, and X-ray absorption spectroscopy techniques. Aqueous speciation of arsenic plays a key role in directing arsenic...

A kinetics database and scripts for PHREEQC

NASA Astrophysics Data System (ADS)

Hu, B.; Zhang, Y.; Teng, Y.; Zhu, C.

2017-12-01

Kinetics of geochemical reactions has been increasingly used in numerical models to simulate coupled flow, mass transport, and chemical reactions. However, the kinetic data are scattered in the literature. To assemble a kinetic dataset for a modeling project is an intimidating task for most. In order to facilitate the application of kinetics in geochemical modeling, we assembled kinetics parameters into a database for the geochemical simulation program, PHREEQC (version 3.0). Kinetics data were collected from the literature. Our database includes kinetic data for over 70 minerals. The rate equations are also programmed into scripts with the Basic language. Using the new kinetic database, we simulated reaction path during the albite dissolution process using various rate equations in the literature. The simulation results with three different rate equations gave difference reaction paths at different time scale. Another application involves a coupled reactive transport model simulating the advancement of an acid plume in an acid mine drainage site associated with Bear Creek Uranium tailings pond. Geochemical reactions including calcite, gypsum, and illite were simulated with PHREEQC using the new kinetic database. The simulation results successfully demonstrated the utility of new kinetic database.

A new statistical model to find bedrock, a prequel to geochemical mass balance

NASA Astrophysics Data System (ADS)

Fisher, B.; Rendahl, A. K.; Aufdenkampe, A. K.; Yoo, K.

2016-12-01

We present a new statistical model to assess weathering trends in deep weathering profiles. The Weathering Trends (WT) model is presented as an extension of the geochemical mass balance model (Brimhall & Dietrich, 1987), and is available as an open-source R library on GitHub (https://github.com/AaronRendahl/WeatheringTrends). WT uses element concentration data to determine the depth to fresh bedrock by assessing the maximum extent of weathering for all elements and the model applies confidence intervals on the depth to bedrock. WT models near-surface features and the shape of the weathering profile using a log transformation of data to capture the magnitude of changes that are relevant to geochemical kinetics and thermodynamics. The WT model offers a new, enhanced opportunity to characterize and understand biogeochemical weathering in heterogeneous rock types. We apply the model to two 21-meter drill cores in the Laurels Schist bedrock in the Christina River Basin Critical Zone Observatory in the Pennsylvania Piedmont. The Laurels Schist had inconclusive weathering indicators prior to development and application of WT model. The model differentiated between rock variability and weathering to delineate the maximum extent of weathering at 12.3 (CI 95% [9.2, 21.3]) meters in Ridge Well 1 and 7.2 (CI 95% [4.3, 13.0]) meters in Interfluve Well 2. The modeled extent to weathering is decoupled from the water table at the ridge, but coincides with the water table at the interfluve. These depths were applied as the parent material for the geochemical mass balance for the Laurels Schist. We test statistical approaches to assess the variability and correlation of immobile elements to facilitate the selection of the best immobile element for use in both models. We apply the model to other published data where the geochemical mass balance was applied, to demonstrate how the WT model provides additional information about weathering depth and weathering trends.

A COMSOL-GEMS interface for modeling coupled reactive-transport geochemical processes

NASA Astrophysics Data System (ADS)

Azad, Vahid Jafari; Li, Chang; Verba, Circe; Ideker, Jason H.; Isgor, O. Burkan

2016-07-01

An interface was developed between COMSOL MultiphysicsTM finite element analysis software and (geo)chemical modeling platform, GEMS, for the reactive-transport modeling of (geo)chemical processes in variably saturated porous media. The two standalone software packages are managed from the interface that uses a non-iterative operator splitting technique to couple the transport (COMSOL) and reaction (GEMS) processes. The interface allows modeling media with complex chemistry (e.g. cement) using GEMS thermodynamic database formats. Benchmark comparisons show that the developed interface can be used to predict a variety of reactive-transport processes accurately. The full functionality of the interface was demonstrated to model transport processes, governed by extended Nernst-Plank equation, in Class H Portland cement samples in high pressure and temperature autoclaves simulating systems that are used to store captured carbon dioxide (CO2) in geological reservoirs.

Geoneutrinos and Heat Production in the Earth: Constraints and Implications

ScienceCinema

McDonough, Bill

2017-12-29

Recent results from antineutrino (geoneutrino) studies at KamLAND are coincident with geochemical models of Th and U in the Earth.  KamLAND and Borexino detectors are on line, thus uncertainties in counting statistics will be reduced as data are accumulated.  The SNO+ detector, situated in the middle of the North American plate will come on line in ~3 yrs and will be best suited to yield a precise estimate of the continental contribution to the Earth’s Th & U budget.  The distribution of heat producing elements in the Earth drives convection and plate tectonics.  Geochemical models posit that ~40% of the heat producing elements are in the continental crust, with the remainder in the mantle.  Although models of core formation allow for the incorporation of heat producing elements, the core contribution of radiogenic heating is considered to be negligible.  Most parameterized convection models for the Earth require significant amounts of radiogenic heating of the Earth, a factor of two greater than geochemical models predict.  The initial KamLAND results challenge these geophysical models and support geochemical models calling for a significant contribution from secular cooling of the mantle.

Manual hierarchical clustering of regional geochemical data using a Bayesian finite mixture model

USGS Publications Warehouse

Ellefsen, Karl J.; Smith, David

2016-01-01

Interpretation of regional scale, multivariate geochemical data is aided by a statistical technique called “clustering.” We investigate a particular clustering procedure by applying it to geochemical data collected in the State of Colorado, United States of America. The clustering procedure partitions the field samples for the entire survey area into two clusters. The field samples in each cluster are partitioned again to create two subclusters, and so on. This manual procedure generates a hierarchy of clusters, and the different levels of the hierarchy show geochemical and geological processes occurring at different spatial scales. Although there are many different clustering methods, we use Bayesian finite mixture modeling with two probability distributions, which yields two clusters. The model parameters are estimated with Hamiltonian Monte Carlo sampling of the posterior probability density function, which usually has multiple modes. Each mode has its own set of model parameters; each set is checked to ensure that it is consistent both with the data and with independent geologic knowledge. The set of model parameters that is most consistent with the independent geologic knowledge is selected for detailed interpretation and partitioning of the field samples.

Geochemistry and the Understanding of Groundwater Systems

NASA Astrophysics Data System (ADS)

Glynn, P. D.; Plummer, L. N.; Weissmann, G. S.; Stute, M.

2009-12-01

Geochemical techniques and concepts have made major contributions to the understanding of groundwater systems. Advances continue to be made through (1) development of measurement and characterization techniques, (2) improvements in computer technology, networks and numerical modeling, (3) investigation of coupled geologic, hydrologic, geochemical and biologic processes, and (4) scaling of individual observations, processes or subsystem models into larger coherent model frameworks. Many applications benefit from progress in these areas, such as: (1) understanding paleoenvironments, in particular paleoclimate, through the use of groundwater archives, (2) assessing the sustainability (recharge and depletion) of groundwater resources, and (3) their vulnerability to contamination, (4) evaluating the capacity and consequences of subsurface waste isolation (e.g. geologic carbon sequestration, nuclear and chemical waste disposal), (5) assessing the potential for mitigation/transformation of anthropogenic contaminants in groundwater systems, and (6) understanding the effect of groundwater lag times in ecosystem-scale responses to natural events, land-use changes, human impacts, and remediation efforts. Obtaining “representative” groundwater samples is difficult and progress in obtaining “representative” samples, or interpreting them, requires new techniques in characterizing groundwater system heterogeneity. Better characterization and simulation of groundwater system heterogeneity (both physical and geochemical) is critical to interpreting the meaning of groundwater “ages”; to understanding and predicting groundwater flow, solute transport, and geochemical evolution; and to quantifying groundwater recharge and discharge processes. Research advances will also come from greater use and progress (1) in the application of environmental tracers to ground water dating and in the analysis of new geochemical tracers (e.g. compound specific isotopic analyses, noble gas isotopes, analyses of natural organic tracers), (2) in inverse geochemical and hydrological modeling, (3) in the understanding and simulation of coupled biological, geological, geochemical and hydrological processes, and (4) in the description and quantification of processes occurring at the boundaries of groundwater systems (e.g. unsaturated zone processes, groundwater/surface water interactions, impacts of changing geomorphology and vegetation). Improvements are needed in the integration of widely diverse information. Better techniques are needed to construct coherent conceptual frameworks from individual observations, simulated or reconstructed information, process models, and intermediate scale models. Iterating between data collection, interpretation, and the application of forward, inverse, and statistical modeling tools is likely to provide progress in this area. Quantifying groundwater system processes by using an open-system thermodynamic approach in a common mass- and energy-flow framework will also facilitate comparison and understanding of diverse processes.

Probabilistic, sediment-geochemical parameterisation of the groundwater compartment of the Netherlands for spatially distributed, reactive transport modelling

NASA Astrophysics Data System (ADS)

Janssen, Gijs; Gunnink, Jan; van Vliet, Marielle; Goldberg, Tanya; Griffioen, Jasper

2017-04-01

Pollution of groundwater aquifers with contaminants as nitrate is a common problem. Reactive transport models are useful to predict the fate of such contaminants and to characterise the efficiency of mitigating or preventive measures. Parameterisation of a groundwater transport model on reaction capacity is a necessary step during building the model. Two Dutch, national programs are combined to establish a methodology for building a probabilistic model on reaction capacity of the groundwater compartment at the national scale: the Geological Survey program and the NHI Netherlands Hydrological Instrument program. Reaction capacity is considered as a series of geochemical characteristics that control acid/base condition, redox condition and sorption capacity. Five primary reaction capacity variables are characterised: 1. pyrite, 2. non-pyrite, reactive iron (oxides, siderite and glauconite), 3. clay fraction, 4. organic matter and 5. Ca-carbonate. Important reaction capacity variables that are determined by more than one solid compound are also deduced: 1. potential reduction capacity (PRC) by pyrite and organic matter, 2. cation-exchange capacity (CEC) by organic matter and clay content, 3. carbonate buffering upon pyrite oxidation (CPBO) by carbonate and pyrite. Statistical properties of these variables are established based on c. 16,000 sediment geochemical analyses. The first tens of meters are characterised based on 25 regions using combinations of lithological class and geological formation as strata. Because of both less data and more geochemical uniformity, the deeper subsurface is characterised in a similar way based on 3 regions. The statistical data is used as input in an algoritm that probabilistically calculates the reaction capacity per grid cell. First, the cumulative frequency distribution (cfd) functions are calculated from the statistical data for the geochemical strata. Second, all voxel cells are classified into the geochemical strata. Third, the cfd functions are used to put random reaction capacity variables into the hydrological voxel model. Here, the distribution can be conditioned on two variables. Two important variables are clay content and depth. The first is valid because more dense data is available for clay content than for geochemical variables as pyrite and probabilistic, lithological models are also built at TNO Geological Survey. The second is important to account for locally different depths at which the redox cline between NO3-rich and Fe(II)-rich groundwater occurs within the first tens of meters of the subsurface. An extensive data-set of groundwater quality analyses is used to derive criteria for depth variability of the redox cline. The result is a unique algoritm in order to obtain heterogeneous geochemical reaction capacity models of the entire groundwater compartment of the Netherlands.

Conceptual models in exploration geochemistry-The Basin and Range Province of the Western United States and Northern Mexico

USGS Publications Warehouse

Lovering, T.G.; McCarthy, J.H.

1978-01-01

This summary of geochemical exploration in the Basin and Range Province is another in the series of reviews of geochemical-exploration applications covering a large region; this series began in 1975 with a summary for the Canadian Cordillera and Canadian Shield, and was followed in 1976 by a similar summary for Scandinavia (Norden). Rather than adhering strictly to the type of conceptual models applied in those papers, we have made use of generalized landscape geochemistry models related to the nature of concealment of ore deposits. This study is part of a continuing effort to examine and evaluate geochemical-exploration practices in different areas of the world. Twenty case histories of the application of geochemical exploration in both district and regional settings illustrate recent developments in techniques and approaches. Along with other published reports these case histories, exemplifying generalized models of concealed deposits, provide data used to evaluate geochemical-exploration programs and specific sample media. Because blind deposits are increasingly sought in the Basin and Range Province, the use of new sample media or anomaly-enhancement techniques is a necessity. Analysis of vapors or gases emanating from blind deposits is a promising new technique. Certain fractions of stream sediments show anomalies that are weak or not detected in conventional minus 80-mesh fractions. Multi-element analysis of mineralized bedrock may show zoning patterns that indicate depth or direction of ore. Examples of the application of these and other, more conventional methods are indicated in the case histories. The final section of this paper contains a brief evaluation of the applications of all types of sample media to geochemical exploration in the arid environment of the Basin and Range Province. ?? 1978.

Monitoring CO2 sequestration into deep saline aquifer and associated salt intrusion using coupled multiphase flow modeling and time-lapse electrical resistivity tomography

NASA Astrophysics Data System (ADS)

Lu, C.; Zhang, C.; Huang, H.; Johnson, T.

2012-12-01

Geological sequestration of carbon dioxide (CO2) into the subsurface has been considered as one solution to reduce greenhouse emission to the atmosphere. Successful sequestration process requires efficient and adequate monitoring of injected fluids as they migrate into the aquifer to evaluate flow path, leakage, and geochemical interactions between CO2 and geologic media. In this synthetic field scale study, we have integrated 3D multiphase flow modeling code PFLOTRAN with 3D time-laps electrical resistivity tomography (ERT) to gain insight into the supercritical (SC) CO2 plumes movement in the deep saline aquifer and associated brine intrusion into shallower fresh water aquifer. A parallel ERT forward and inverse modeling package was introduced, and related algorithms are briefly described. The capabilities and limitations of ERT in monitoring CO2 migration are assessed by comparing the results from PFLOTRAN simulations with the ERT inversion results. In general, our study shows the ERT inversion results compare well with PFLOTRAN with reasonable discrepancies, indicating that the ERT can capture the actual CO2 plume dynamics and brine intrusion. Detailed comparisons on the location, size and volume of CO2 plume show the ERT method underestimated area review and overestimated total plume volume in the predictions of SC CO2 movements. These comparisons also show the ERT method constantly overestimate salt intrusion area and underestimated total solute amount in the predictions of brine filtration. Our study shows that together with other geochemical and geophysical methods, ERT is a potentially useful monitoring tool in detecting the SC CO2 and formation fluid migrations.

Hydrogeology, ground-water quality, and source of ground water causing water-quality changes in the Davis well field at Memphis, Tennessee

USGS Publications Warehouse

Parks, William S.; Mirecki, June E.; Kingsbury, James A.

1995-01-01

NETPATH geochemical model code was used to mix waters from the alluvial aquifer with water from the Memphis aquifer using chloride as a conservative tracer. The resulting models indicated that a mixture containing 3 percent alluvial aquifer water mixed with 97 percent unaffected Memphis aquifer water would produce the chloride concentration measured in water from the Memphis aquifer well most affected by water-quality changes. NETPATH also was used to calculate mixing percentages of alluvial and Memphis aquifer Abstract waters based on changes in the concentrations of selected dissolved major inorganic and trace element constituents that define the dominant reactions that occur during mixing. These models indicated that a mixture containing 18 percent alluvial aquifer water and 82 percent unaffected Memphis aquifer water would produce the major constituent and trace element concentrations measured in water from the Memphis aquifer well most affected by water-quality changes. However, these model simulations predicted higher dissolved methane concentrations than were measured in water samples from the Memphis aquifer wells.

A scaling relationship for impact-induced melt volume

NASA Astrophysics Data System (ADS)

Nakajima, M.; Rubie, D. C.; Melosh, H., IV; Jacobson, S. A.; Golabek, G.; Nimmo, F.; Morbidelli, A.

2016-12-01

During the late stages of planetary accretion, protoplanets experience a number of giant impacts and extensive mantle melting. The impactor's core sinks through the molten part of the target mantle (magma ocean) and experiences metal-silicate partitioning (e.g., Stevenson, 1990). For understanding the chemical evolution of the planetary mantle and core, we need to determine the impact-induced melt volume because the partitioning strongly depends on the ranges of the pressures and temperatures within the magma ocean. Previous studies have investigated the effects of small impacts (i.e. impact cratering) on melt volume, but those for giant impacts are not well understood yet. Here, we perform giant impact simulations to derive a scaling law for melt volume as a function of impact velocity, impact angle, and impactor-to-target mass ratio. We use two different numerical codes, namely smoothed particle hydrodynamics we developed (SPH, a particle method) and the code iSALE (a grid-based method) to compare their outcomes. Our simulations show that these two codes generally agree as long as the same equation of state is used. We also find that some of the previous studies developed for small impacts (e.g., Abramov et al., 2012) overestimate giant impact melt volume by orders of magnitudes partly because these models do not consider self-gravity of the impacting bodies. Therefore, these models may not be extrapolated to large impacts. Our simulations also show that melt volume can be scaled by the total mass of the system. In this presentation, we further discuss geochemical implications for giant impacts on planets, including Earth and Mars.

Geochemical Data Package for Performance Assessment Calculations Related to the Savannah River Site

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kaplan, Daniel I.

The Savannah River Site (SRS) disposes of low-level radioactive waste (LLW) and stabilizes high-level radioactive waste (HLW) tanks in the subsurface environment. Calculations used to establish the radiological limits of these facilities are referred to as Performance Assessments (PA), Special Analyses (SA), and Composite Analyses (CA). The objective of this document is to revise existing geochemical input values used for these calculations. This work builds on earlier compilations of geochemical data (2007, 2010), referred to a geochemical data packages. This work is being conducted as part of the on-going maintenance program of the SRS PA programs that periodically updates calculationsmore » and data packages when new information becomes available. Because application of values without full understanding of their original purpose may lead to misuse, this document also provides the geochemical conceptual model, the approach used for selecting the values, the justification for selecting data, and the assumptions made to assure that the conceptual and numerical geochemical models are reasonably conservative (i.e., bias the recommended input values to reflect conditions that will tend to predict the maximum risk to the hypothetical recipient). This document provides 1088 input parameters for geochemical parameters describing transport processes for 64 elements (>740 radioisotopes) potentially occurring within eight subsurface disposal or tank closure areas: Slit Trenches (ST), Engineered Trenches (ET), Low Activity Waste Vault (LAWV), Intermediate Level (ILV) Vaults, Naval Reactor Component Disposal Areas (NRCDA), Components-in-Grout (CIG) Trenches, Saltstone Facility, and Closed Liquid Waste Tanks. The geochemical parameters described here are the distribution coefficient, Kd value, apparent solubility concentration, k s value, and the cementitious leachate impact factor.« less

Reactive transport modeling at uranium in situ recovery sites: uncertainties in uranium sorption on iron hydroxides

USGS Publications Warehouse

Johnson, Raymond H.; Tutu, Hlanganani; Brown, Adrian; Figueroa, Linda; Wolkersdorfer, Christian

2013-01-01

Geochemical changes that can occur down gradient from uranium in situ recovery (ISR) sites are important for various stakeholders to understand when evaluating potential effects on surrounding groundwater quality. If down gradient solid-phase material consists of sandstone with iron hydroxide coatings (no pyrite or organic carbon), sorption of uranium on iron hydroxides can control uranium mobility. Using one-dimensional reactive transport models with PHREEQC, two different geochemical databases, and various geochemical parameters, the uncertainties in uranium sorption on iron hydroxides are evaluated, because these oxidized zones create a greater risk for future uranium transport than fully reduced zones where uranium generally precipitates.

Mercury Slovenian soils: High, medium and low sample density geochemical maps

NASA Astrophysics Data System (ADS)

Gosar, Mateja; Šajn, Robert; Teršič, Tamara

2017-04-01

Regional geochemical survey was conducted in whole territory of Slovenia (20273 km2). High, medium and low sample density surveys were compared. High sample density represented the regional geochemical data set supplemented by local high-density sampling data (irregular grid, n=2835). Medium-density soil sampling was performed in a 5 x 5 km grid (n=817) and low-density geochemical survey was conducted in a sampling grid 25 x 25 km (n=54). Mercury distribution in Slovenian soils was determined with models of mercury distribution in soil using all three data sets. A distinct Hg anomaly in western part of Slovenia is evident on all three models. It is a consequence of 500-years of mining and ore processing in the second largest mercury mine in the world, the Idrija mine. The determined mercury concentrations revealed an important difference between the western and the eastern parts of the country. For the medium scale geochemical mapping is the median value (0.151 mg /kg) for western Slovenia almost 2-fold higher than the median value (0.083 mg/kg) in eastern Slovenia. Besides the Hg median for the western part of Slovenia exceeds the Hg median for European soil by a factor of 4 (Gosar et al., 2016). Comparing these sample density surveys, it was shown that high sampling density allows the identification and characterization of anthropogenic influences on a local scale, while medium- and low-density sampling reveal general trends in the mercury spatial distribution, but are not appropriate for identifying local contamination in industrial regions and urban areas. The resolution of the pattern generated is the best when the high-density survey on a regional scale is supplemented with the geochemical data of the high-density surveys on a local scale. References: Gosar, M, Šajn, R, Teršič, T. Distribution pattern of mercury in the Slovenian soil: geochemical mapping based on multiple geochemical datasets. Journal of geochemical exploration, 2016, 167/38-48.

Reactive transport and mass balance modeling of the Stimson sedimentary formation and altered fracture zones constrain diagenetic conditions at Gale crater, Mars

NASA Astrophysics Data System (ADS)

Hausrath, E. M.; Ming, D. W.; Peretyazhko, T. S.; Rampe, E. B.

2018-06-01

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosity in Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, τ, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past martian climate and environments.

A geochemical and geophysical reappraisal to the significance of the recent unrest at Campi Flegrei caldera (Southern Italy)

NASA Astrophysics Data System (ADS)

Moretti, Roberto; De Natale, Giuseppe; Troise, Claudia

2017-04-01

Volcanic unrest at calderas involve complex interaction between magma, hydrothermal fluids and crustal stress and strain. Campi Flegrei caldera (CFc), located in the Naples (Italy) area and characterised by the highest volcanic risk on Earth for the extreme urbanisation, undergoes unrest phenomena involving several meters of uplift and intense shallow micro-seismicity since several decades. Despite unrest episodes display in the last decade only moderate ground deformation and seismicity, current interpretations of geochemical data point to a highly pressurized hydrothermal system. We show that at CFc, the usual assumption of vapour-liquid coexistence in the fumarole plumes leads to largely overestimated hydrothermal pressures and, accordingly, interpretations of elevated unrest. By relaxing unconstrained geochemical assumptions, we infer an alternative model yielding better agreement between geophysical and geochemical observations. The model reconciles discrepancies between what observed 1) for two decades since the 1982-84 large unrest, when shallow magma was supplying heat and fluids to the hydrothermal system, and 2) in the last decade. Compared to the 1980's unrest, the post-2005 phenomena are characterized by much lower aquifers overpressure and magmatic involvement, as indicated by geophysical data and despite large changes in geochemical indicators. Our interpretation points out a model in which shallow sills, intruded during 1969-1984, have completely cooled, so that fumarole emissions are affected now by deeper, CO2-richer, magmatic gases producing a relatively modest heating and overpressure of the hydrothermal system. Our results do have important implications on the short-term eruption hazard assessment and on the best strategies for monitoring and interpreting geochemical data.

Experimental validation of Swy-2 clay standard's PHREEQC model

NASA Astrophysics Data System (ADS)

Szabó, Zsuzsanna; Hegyfalvi, Csaba; Freiler, Ágnes; Udvardi, Beatrix; Kónya, Péter; Székely, Edit; Falus, György

2017-04-01

One of the challenges of the present century is to limit the greenhouse gas emissions for the mitigation of climate change which is possible for example by a transitional technology, CCS (Carbon Capture and Storage) and, among others, by the increase of nuclear proportion in the energy mix. Clay minerals are considered to be responsible for the low permeability and sealing capacity of caprocks sealing off stored CO2 and they are also the main constituents of bentonite in high level radioactive waste disposal facilities. The understanding of clay behaviour in these deep geological environments is possible through laboratory batch experiments of well-known standards and coupled geochemical models. Such experimentally validated models are scarce even though they allow deriving more precise long-term predictions of mineral reactions and rock and bentonite degradation underground and, therefore, ensuring the safety of the above technologies and increase their public acceptance. This ongoing work aims to create a kinetic geochemical model of Na-montmorillonite standard Swy-2 in the widely used PHREEQC code, supported by solution and mineral composition results from batch experiments. Several four days experiments have been carried out in 1:35 rock:water ratio at atmospheric conditions, and with inert and CO2 supercritical phase at 100 bar and 80 ⁰C relevant for the potential Hungarian CO2 reservoir complex. Solution samples have been taken during and after experiments and their compositions were measured by ICP-OES. The treated solid phase has been analysed by XRD and ATR-FTIR and compared to in-parallel measured references (dried Swy-2). Kinetic geochemical modelling of the experimental conditions has been performed by PHREEQC version 3 using equations and kinetic rate parameters from the USGS report of Palandri and Kharaka (2004). The visualization of experimental and numerous modelling results has been automatized by R. Experiments and models show very fast reactions under the studied conditions and increased reactivity in presence of scCO2. A model sensitivity analysis has pointed out that the continuously changing solution composition results cannot be described by the change of the uncertain reactive surface area of mineral phases in the model and still several orders of magnitude different ion-concentrations are predicted. However, by considering the clay standard's cation exchange capacity divided proportionally among interlayer cations of Na-montmorillonite, the measured variation can be described on an order of magnitude level. It is furthermore indicated that not only the interlayer cations take part in this process but a minor proportion of other, structural ions as well, differently in the reference and scCO2 environments. Experimental methodological aspects of the work, such as solution sampling, solid sample post-experimental treatment, solution and solid sample analysis sensitivity, expected experimental by-products etc. are also to be addressed.

Inverse Modeling of Water-Rock-CO2 Batch Experiments: Potential Impacts on Groundwater Resources at Carbon Sequestration Sites.

PubMed

Yang, Changbing; Dai, Zhenxue; Romanak, Katherine D; Hovorka, Susan D; Treviño, Ramón H

2014-01-01

This study developed a multicomponent geochemical model to interpret responses of water chemistry to introduction of CO2 into six water-rock batches with sedimentary samples collected from representative potable aquifers in the Gulf Coast area. The model simulated CO2 dissolution in groundwater, aqueous complexation, mineral reactions (dissolution/precipitation), and surface complexation on clay mineral surfaces. An inverse method was used to estimate mineral surface area, the key parameter for describing kinetic mineral reactions. Modeling results suggested that reductions in groundwater pH were more significant in the carbonate-poor aquifers than in the carbonate-rich aquifers, resulting in potential groundwater acidification. Modeled concentrations of major ions showed overall increasing trends, depending on mineralogy of the sediments, especially carbonate content. The geochemical model confirmed that mobilization of trace metals was caused likely by mineral dissolution and surface complexation on clay mineral surfaces. Although dissolved inorganic carbon and pH may be used as indicative parameters in potable aquifers, selection of geochemical parameters for CO2 leakage detection is site-specific and a stepwise procedure may be followed. A combined study of the geochemical models with the laboratory batch experiments improves our understanding of the mechanisms that dominate responses of water chemistry to CO2 leakage and also provides a frame of reference for designing monitoring strategy in potable aquifers.

SubductionGenerator: A program to build three-dimensional plate configurations

NASA Astrophysics Data System (ADS)

Jadamec, M. A.; Kreylos, O.; Billen, M. I.; Turcotte, D. L.; Knepley, M.

2016-12-01

Geologic, geochemical, and geophysical data from subduction zones indicate that a two-dimensional paradigm for plate tectonic boundaries is no longer adequate to explain the observations. Many open source software packages exist to simulate the viscous flow of the Earth, such as the dynamics of subduction. However, there are few open source programs that generate the three-dimensional model input. We present an open source software program, SubductionGenerator, that constructs the three-dimensional initial thermal structure and plate boundary structure. A 3D model mesh and tectonic configuration are constructed based on a user specified model domain, slab surface, seafloor age grid file, and shear zone surface. The initial 3D thermal structure for the plates and mantle within the model domain is then constructed using a series of libraries within the code that use a half-space cooling model, plate cooling model, and smoothing functions. The code maps the initial 3D thermal structure and the 3D plate interface onto the mesh nodes using a series of libraries including a k-d tree to increase efficiency. In this way, complicated geometries and multiple plates with variable thickness can be built onto a multi-resolution finite element mesh with a 3D thermal structure and 3D isotropic shear zones oriented at any angle with respect to the grid. SubductionGenerator is aimed at model set-ups more representative of the earth, which can be particularly challenging to construct. Examples include subduction zones where the physical attributes vary in space, such as slab dip and temperature, and overriding plate temperature and thickness. Thus, the program can been used to construct initial tectonic configurations for triple junctions and plate boundary corners.

A Quantitative Geochemical Target for Modeling the Formation of the Earth and Moon

NASA Technical Reports Server (NTRS)

Boyce, Jeremy W.; Barnes, Jessica J.; McCubbin, Francis M.

2017-01-01

The past decade has been one of geochemical, isotopic, and computational advances that are bringing the laboratory measurements and computational modeling neighborhoods of the Earth-Moon community to ever closer proximity. We are now however in the position to become even better neighbors: modelers can generate testable hypthotheses for geochemists; and geochemists can provide quantitive targets for modelers. Here we present a robust example of the latter based on Cl isotope measurements of mare basalts.

A reactive transport model for mercury fate in soil--application to different anthropogenic pollution sources.

PubMed

Leterme, Bertrand; Blanc, Philippe; Jacques, Diederik

2014-11-01

Soil systems are a common receptor of anthropogenic mercury (Hg) contamination. Soils play an important role in the containment or dispersion of pollution to surface water, groundwater or the atmosphere. A one-dimensional model for simulating Hg fate and transport for variably saturated and transient flow conditions is presented. The model is developed using the HP1 code, which couples HYDRUS-1D for the water flow and solute transport to PHREEQC for geochemical reactions. The main processes included are Hg aqueous speciation and complexation, sorption to soil organic matter, dissolution of cinnabar and liquid Hg, and Hg reduction and volatilization. Processes such as atmospheric wet and dry deposition, vegetation litter fall and uptake are neglected because they are less relevant in the case of high Hg concentrations resulting from anthropogenic activities. A test case is presented, assuming a hypothetical sandy soil profile and a simulation time frame of 50 years of daily atmospheric inputs. Mercury fate and transport are simulated for three different sources of Hg (cinnabar, residual liquid mercury or aqueous mercuric chloride), as well as for combinations of these sources. Results are presented and discussed with focus on Hg volatilization to the atmosphere, Hg leaching at the bottom of the soil profile and the remaining Hg in or below the initially contaminated soil layer. In the test case, Hg volatilization was negligible because the reduction of Hg(2+) to Hg(0) was inhibited by the low concentration of dissolved Hg. Hg leaching was mainly caused by complexation of Hg(2+) with thiol groups of dissolved organic matter, because in the geochemical model used, this reaction only had a higher equilibrium constant than the sorption reactions. Immobilization of Hg in the initially polluted horizon was enhanced by Hg(2+) sorption onto humic and fulvic acids (which are more abundant than thiols). Potential benefits of the model for risk management and remediation of contaminated sites are discussed.

The use of multifractal modelling for targeting resources from soil and stream geochemistry data: the case of the Variscan basement of the Iberian Peninsula

NASA Astrophysics Data System (ADS)

Gonçalves, Mario; Mateus, Antonio

2016-04-01

The safeguarding of access/use of many critical raw materials for Society requires that much of previously dismissed areas for exploration must be re-evaluated with new criteria in which the significance of "anomaly" should not be treated independently of the geochemical signals of the ore-forming processes and how the different chemical elements are interrelated. For much of the previous decade, several multifractal methods were methodically being refined as automatic tools to analyze and detect geochemical anomalies. These included the early concentration-area method (Cheng et al., 1994), singularity mapping (Cheng, 2007), and spectrum-area (Cheng et al., 2000), which has been recently combined with the bi-dimensional empirical mode decomposition (Xu et al., 2016) as a tool to separate different contributing sources of an otherwise complex geochemical pattern. We propose yet another approach, the use of geochemical indexes, which links to the geological and ore-forming processes known to define a given region in order to assess much of these numerical approaches. Therefore, we picked several areas from the Variscan basement in Portugal, with different geologic and metallogentic contexts, some of them previously analyzed with multifractal methods (Gonçalves et al., 2001; Jesus et al., 2013) and a multi-element geochemical campaign on which to test the different multifractal methods combined with the geochemical indexes, as an advantageous alternative to principal component mapping, for example. Some preliminary essays with stochastic models similar to those reported in Gonçalves (2001) and Agterberg (2007), with different overprinted pulses are presented as well. Acknowledgments: This is a contribution from UID/GEO/50019/2013 - Instituto Dom Luiz, supported by FCT. Agterberg, 2007, Math. Geol., 39, 1. Cheng et al, 1994, J. Geochem. Explor., 51, 109. Cheng et al., 2000, Nat. Resour. Res, 9, 43. Cheng, 2007, Ore Geol. Rev., 32, 314. Gonçalves, 2001, Math. Geol., 33, 41. Gonçalves et al., 2001, J. Geochem. Explor., 72, 91. Jesus et al., 2013, J. Geochem. Explor., 126-127, 23. Xu et al., 2016, J. Geochem. Explor., in press

Compositional data analysis and geochemical modeling of CO2-water-rock interactions in three provinces of Korea.

PubMed

Kim, Seong Hee; Choi, Byoung-Young; Lee, Gyemin; Yun, Seong-Taek; Kim, Soon-Oh

2017-12-20

The CO 2 -rich spring water (CSW) occurring naturally in three provinces, Kangwon (KW), Chungbuk (CB), and Gyeongbuk (GB) of South Korea was classified based on its hydrochemical properties using compositional data analysis. Additionally, the geochemical evolution pathways of various CSW were simulated via equilibrium phase modeling (EPM) incorporated in the PHREEQC code. Most of the CSW in the study areas grouped into the Ca-HCO 3 water type, but some samples from the KW area were classified as Na-HCO 3 water. Interaction with anorthite is likely to be more important than interaction with carbonate minerals for the hydrochemical properties of the CSW in the three areas, indicating that the CSW originated from interactions among magmatic CO 2 , deep groundwater, and bedrock-forming minerals. Based on the simulation results of PHREEQC EPM, the formation temperatures of the CSW within each area were estimated as 77.8 and 150 °C for the Ca-HCO 3 and Na-HCO 3 types of CSW, respectively, in the KW area; 138.9 °C for the CB CSW; and 93.0 °C for the GB CSW. Additionally, the mixing ratios between simulated carbonate water and shallow groundwater were adjusted to 1:9-9:1 for the CSW of the GB area and the Ca-HCO 3 -type CSW of the KW area, indicating that these CSWs were more affected by carbonate water than by shallow groundwater. On the other hand, mixing ratios of 1:9-5:5 and 1:9-3:7 were found for the Na-HCO 3 -type CSW of the KW area and for the CSW of the CB area, respectively, suggesting a relatively small contribution of carbonate water to these CSWs. This study proposes a systematic, but relatively simple, methodology to simulate the formation of carbonate water in deep environments and the geochemical evolution of CSW. Moreover, the proposed methodology could be applied to predict the behavior of CO 2 after its geological storage and to estimate the stability and security of geologically stored CO 2 .

Generic Procedure for Coupling the PHREEQC Geochemical Modeling Framework with Flow and Solute Transport Simulators

NASA Astrophysics Data System (ADS)

Wissmeier, L. C.; Barry, D. A.

2009-12-01

Computer simulations of water availability and quality play an important role in state-of-the-art water resources management. However, many of the most utilized software programs focus either on physical flow and transport phenomena (e.g., MODFLOW, MT3DMS, FEFLOW, HYDRUS) or on geochemical reactions (e.g., MINTEQ, PHREEQC, CHESS, ORCHESTRA). In recent years, several couplings between both genres of programs evolved in order to consider interactions between flow and biogeochemical reactivity (e.g., HP1, PHWAT). Software coupling procedures can be categorized as ‘close couplings’, where programs pass information via the memory stack at runtime, and ‘remote couplings’, where the information is exchanged at each time step via input/output files. The former generally involves modifications of software codes and therefore expert programming skills are required. We present a generic recipe for remotely coupling the PHREEQC geochemical modeling framework and flow and solute transport (FST) simulators. The iterative scheme relies on operator splitting with continuous re-initialization of PHREEQC and the FST of choice at each time step. Since PHREEQC calculates the geochemistry of aqueous solutions in contact with soil minerals, the procedure is primarily designed for couplings to FST’s for liquid phase flow in natural environments. It requires the accessibility of initial conditions and numerical parameters such as time and space discretization in the input text file for the FST and control of the FST via commands to the operating system (batch on Windows; bash/shell on Unix/Linux). The coupling procedure is based on PHREEQC’s capability to save the state of a simulation with all solid, liquid and gaseous species as a PHREEQC input file by making use of the dump file option in the TRANSPORT keyword. The output from one reaction calculation step is therefore reused as input for the following reaction step where changes in element amounts due to advection/dispersion are introduced as irreversible reactions. An example for the coupling of PHREEQC and MATLAB for the solution of unsaturated flow and transport is provided.

Surface water quality assessment of the Upper Illinois River basin in Illinois, Indiana, and Wisconsin : geochemical data for fine-fraction streambed sediment from high- and low-order streams, 1987

USGS Publications Warehouse

Colman, John A.; Sanzolone, R.F.

1991-01-01

Geochemical data are presented from a synoptic survey of 46 elements in fine-fraction streambed sediments of the Upper Illinois River Basin during the fall of 1987. The survey was a component study of the Illinois pilot project of the U.S. Geological Survey's National Water-Quality Assessment program. Most of the sampling sites were randomly chosen--135 on main stems of rivers and 238 on first- and second-order streams. In addition, 196 samples were collected for quality-assurance and special-study purposes. The report includes element concentration data and summary-statistics tables of percentiles, nested analysis of variance, and correlation coefficients. All concentration data are included in tabular form and can be selected by map reference number, latitude and longitude, or remark code indicating purpose for collecting sample.

Segmentation of singularity maps in the context of soil porosity

NASA Astrophysics Data System (ADS)

Martin-Sotoca, Juan J.; Saa-Requejo, Antonio; Grau, Juan; Tarquis, Ana M.

2016-04-01

Geochemical exploration have found with increasingly interests and benefits of using fractal (power-law) models to characterize geochemical distribution, including concentration-area (C-A) model (Cheng et al., 1994; Cheng, 2012) and concentration-volume (C-V) model (Afzal et al., 2011) just to name a few examples. These methods are based on the singularity maps of a measure that at each point define areas with self-similar properties that are shown in power-law relationships in Concentration-Area plots (C-A method). The C-A method together with the singularity map ("Singularity-CA" method) define thresholds that can be applied to segment the map. Recently, the "Singularity-CA" method has been applied to binarize 2D grayscale Computed Tomography (CT) soil images (Martin-Sotoca et al, 2015). Unlike image segmentation based on global thresholding methods, the "Singularity-CA" method allows to quantify the local scaling property of the grayscale value map in the space domain and determinate the intensity of local singularities. It can be used as a high-pass-filter technique to enhance high frequency patterns usually regarded as anomalies when applied to maps. In this work we will put special attention on how to select the singularity thresholds in the C-A plot to segment the image. We will compare two methods: 1) cross point of linear regressions and 2) Wavelets Transform Modulus Maxima (WTMM) singularity function detection. REFERENCES Cheng, Q., Agterberg, F. P. and Ballantyne, S. B. (1994). The separation of geochemical anomalies from background by fractal methods. Journal of Geochemical Exploration, 51, 109-130. Cheng, Q. (2012). Singularity theory and methods for mapping geochemical anomalies caused by buried sources and for predicting undiscovered mineral deposits in covered areas. Journal of Geochemical Exploration, 122, 55-70. Afzal, P., Fadakar Alghalandis, Y., Khakzad, A., Moarefvand, P. and Rashidnejad Omran, N. (2011) Delineation of mineralization zones in porphyry Cu deposits by fractal concentration-volume modeling. Journal of Geochemical Exploration, 108, 220-232. Martín-Sotoca, J. J., Tarquis, A. M., Saa-Requejo, A. and Grau, J. B. (2015). Pore detection in Computed Tomography (CT) soil images through singularity map analysis. Oral Presentation in PedoFract VIII Congress (June, La Coruña - Spain).

A geochemical examination of humidity cell tests

USGS Publications Warehouse

Maest, Ann; Nordstrom, D. Kirk

2017-01-01

Humidity cell tests (HCTs) are long-term (20 to >300 weeks) leach tests that are considered by some to be the among the most reliable geochemical characterization methods for estimating the leachate quality of mined materials. A number of modifications have been added to the original HCT method, but the interpretation of test results varies widely. We suggest that the HCTs represent an underutilized source of geochemical data, with a year-long test generating approximately 2500 individual chemical data points. The HCT concentration peaks and valleys can be thought of as a “chromatogram” of reactions that may occur in the field, whereby peaks in concentrations are associated with different geochemical processes, including sulfate salt dissolution, sulfide oxidation, and dissolution of rock-forming minerals, some of which can neutralize acid. Some of these reactions occur simultaneously, some do not, and geochemical modeling can be used to help distinguish the dominant processes. Our detailed examination, including speciation and inverse modeling, of HCTs from three projects with different geology and mineralization shows that rapid sulfide oxidation dominates over a limited period of time that starts between 40 and 200 weeks of testing. The applicability of laboratory tests results to predicting field leachate concentrations, loads, or rates of reaction has not been adequately demonstrated, although early flush releases and rapid sulfide oxidation rates in HCTs should have some relevance to field conditions. Knowledge of possible maximum solute concentrations is needed to design effective treatment and mitigation approaches. Early flush and maximum sulfide oxidation results from HCTs should be retained and used in environmental models. Factors that complicate the use of HCTs include: sample representation, time for microbial oxidizers to grow, sample storage before testing, geochemical reactions that add or remove constituents, and the HCT results chosen for use in modeling the environmental performance at mine sites. Improved guidance is needed for more consistent interpretation and use of HCT results that rely on identifying: the geochemical processes; the mineralogy, including secondary mineralogy; the available surface area for reactions; and the influence of hydrologic processes on leachate concentrations in runoff, streams, and groundwater.

Geochemical Constraints on the Size of the Moon-Forming Giant Impact

NASA Astrophysics Data System (ADS)

Piet, Hélène; Badro, James; Gillet, Philippe

2017-12-01

Recent models involving the Moon-forming giant impact hypothesis have managed to reproduce the striking isotopic similarity between the two bodies, albeit using two extreme models: one involves a high-energy small impactor that makes the Moon out of Earth's proto-mantle; the other supposes a gigantic collision between two half-Earths creating the Earth-Moon system from both bodies. Here we modeled the geochemical influence of the giant impact on Earth's mantle and found that impactors larger than 15% of Earth mass result in mantles always violating the present-day concentrations of four refractory moderately siderophile trace elements (Ni, Co, Cr, and V). In the aftermath of the impact, our models cannot further discriminate between a fully and a partially molten bulk silicate Earth. Then, the preservation of primordial geochemical reservoirs predating the Moon remains the sole argument against a fully molten mantle after the Moon-forming impact.

Selected Geochemical Data for Modeling Near-Surface Processes in Mineral Systems

USGS Publications Warehouse

Giles, Stuart A.; Granitto, Matthew; Eppinger, Robert G.

2009-01-01

The database herein was initiated, designed, and populated to collect and integrate geochemical, geologic, and mineral deposit data in an organized manner to facilitate geoenvironmental mineral deposit modeling. The Microsoft Access database contains data on a variety of mineral deposit types that have variable environmental effects when exposed at the ground surface by mining or natural processes. The data tables describe quantitative and qualitative geochemical analyses determined by 134 analytical laboratory and field methods for over 11,000 heavy-mineral concentrate, rock, sediment, soil, vegetation, and water samples. The database also provides geographic information on geology, climate, ecoregion, and site contamination levels for over 3,000 field sites in North America.

Obtaining higher-accuracy estimates of water-rich rocks and water-poor sand dunes on Mars in active neutron experiments

NASA Astrophysics Data System (ADS)

Gabriel, T. S. J.; Hardgrove, C.; Litvak, M. L.; Nowicki, S.; Mitrofanov, I. G.; Boynton, W. V.; Fedosov, F.; Golovin, D.; Jun, I.; Mischna, M.; Tate, C. G.; Moersch, J.; Harshman, K.; Kozyrev, A.; Malakhov, A. V.; Mokrousov, M.; Nikiforov, S.; Sanin, A. B.; Vostrukhin, A.; Thompson, L. M.

2017-12-01

The Dynamic Albedo of Neutrons (DAN) experiment on the Mars Science Laboratory Curiosity Rover delivers high-energy (14.1 MeV) pulses of neutrons into the surface when operating in "active" mode. Neutrons are moderated in the subsurface and return to two detectors to provide a time-of-flight profile in 64 time-bins in epithermal and thermal energy ranges. Results are compared to simulations of the experiment in the Monte Carlo N-Particle Transport Code where several aspects are modeled including the DAN detectors, neutron source, rover components, and underlying rock. Models can be improved by increasing the fidelity of the rock geochemistry as informed by instruments including the Alpha Particle X-Ray Spectrometer (APXS). Furthermore, increasing the fidelity of the rock morphology in models is enabled by the suite of imaging instruments on the rover.To rapidly interpret DAN data a set of pre-simulated generic rock density and bulk geochemistry models are compared to several DAN active observations. While, to first order, this methodology provides an indication of significant geochemical changes in the subsurface, higher-fidelity models should be used to provide accurate constraints on water content, depth of geologic layers, or abundance of neutron absorbers. For example, in high-silicon, low-iron rocks observed along the rover's traverse, generic models can differ by several wt%H2O from models that use APXS measurements of nearby drill samples. Accurate measurements of high-silicon targets are necessary in outlining the extent of aqueous alteration and hydrothermal activity in Gale Crater. Additionally, we find that for DAN active experiments over sand dunes best-fit models can differ by greater than 0.5 wt%HO when the upper layer density is reduced by 0.6 g/cm3 to account for the low-bulk density of sand. In areas where the rock geochemistry differs little from generic models the difference in results is expectedly less disparate. We report refined wt%HO values for high-silicon, aqueously-altered rock and comparatively dry sand dunes along the rover traverse. We also outline the methodology for providing accurate geochemical and morphological constraints using DAN active measurements.

Modeling Anaerobic Soil Organic Carbon Decomposition in Arctic Polygon Tundra: Insights into Soil Geochemical Influences on Carbon Mineralization: Modeling Archive

DOE Data Explorer

Zheng, Jianqiu; Thornton, Peter; Painter, Scott; Gu, Baohua; Wullschleger, Stan; Graham, David

2018-06-13

This anaerobic carbon decomposition model is developed with explicit representation of fermentation, methanogenesis and iron reduction by combining three well-known modeling approaches developed in different disciplines. A pool-based model to represent upstream carbon transformations and replenishment of DOC pool, a thermodynamically-based model to calculate rate kinetics and biomass growth for methanogenesis and Fe(III) reduction, and a humic ion-binding model for aqueous phase speciation and pH calculation are implemented into the open source geochemical model PHREEQC (V3.0). Installation of PHREEQC is required to run this model.

Strategies to predict metal mobility in surficial mining environments

USGS Publications Warehouse

Smith, Kathleen S.

2007-01-01

This report presents some strategies to predict metal mobility at mining sites. These strategies are based on chemical, physical, and geochemical information about metals and their interactions with the environment. An overview of conceptual models, metal sources, and relative mobility of metals under different geochemical conditions is presented, followed by a discussion of some important physical and chemical properties of metals that affect their mobility, bioavailability, and toxicity. The physical and chemical properties lead into a discussion of the importance of the chemical speciation of metals. Finally, environmental and geochemical processes and geochemical barriers that affect metal speciation are discussed. Some additional concepts and applications are briefly presented at the end of this report.

Singularity analysis based on wavelet transform of fractal measures for identifying geochemical anomaly in mineral exploration

NASA Astrophysics Data System (ADS)

Chen, Guoxiong; Cheng, Qiuming

2016-02-01

Multi-resolution and scale-invariance have been increasingly recognized as two closely related intrinsic properties endowed in geofields such as geochemical and geophysical anomalies, and they are commonly investigated by using multiscale- and scaling-analysis methods. In this paper, the wavelet-based multiscale decomposition (WMD) method was proposed to investigate the multiscale natures of geochemical pattern from large scale to small scale. In the light of the wavelet transformation of fractal measures, we demonstrated that the wavelet approximation operator provides a generalization of box-counting method for scaling analysis of geochemical patterns. Specifically, the approximation coefficient acts as the generalized density-value in density-area fractal modeling of singular geochemical distributions. Accordingly, we presented a novel local singularity analysis (LSA) using the WMD algorithm which extends the conventional moving averaging to a kernel-based operator for implementing LSA. Finally, the novel LSA was validated using a case study dealing with geochemical data (Fe2O3) in stream sediments for mineral exploration in Inner Mongolia, China. In comparison with the LSA implemented using the moving averaging method the novel LSA using WMD identified improved weak geochemical anomalies associated with mineralization in covered area.

Geochemical Constraints on the Size of the Moon — Forming Giant Impact

NASA Astrophysics Data System (ADS)

Piet, H.; Badro, J.; Gillet, P.

2018-05-01

We use the partitioning of siderophile trace elements to model the geochemical influence of the Moon-forming giant impact on Earth’s mantle during core formation. We find the size of the impactor to be 15% of Earth mass or smaller.

Chemical mixing model studies of lunar orbital geochemical data - Apollo 16 and 17 highlands compositions

NASA Technical Reports Server (NTRS)

Spudis, P. D.; Hawke, B. R.

1982-01-01

Chemical mixing model studies of lunar geochemical data for the central and Taurus-Littrow lunar highlands were performed utilizing pristine highland rock types as end member compositions. The central highlands show considerable diversity in composition; anorthosite is the principal rock type in the Apollo 16/Descartes region, while norite predominates in the highlands west of the landing site. This change in crustal composition is coincident with a major color boundary seen in earth-based multispectral data and probably represents the presence of distinct geochemical provinces within the central highlands. The Taurus-Littrow highlands are dominated by norite; anorthosite is far less abundant than in the central highlands. This suggests that the impact target for the Serenitatis basin was different than that of the Nectaris basin and further strengthens the hypothesis that the lunar highlands are petrologically heterogeneous on a regional basis. It is suggested that the lunar highlands should be viewed in terms of geochemical provinces that have undergone distinct and complex igneous and impact histories.

Establishing nursery estuary otolith geochemical tags for Sea Bass (Dicentrarchus labrax): Is temporal stability estuary dependent?

NASA Astrophysics Data System (ADS)

Ryan, Diarmuid; Wögerbauer, Ciara; Roche, William

2016-12-01

The ability to determine connectivity between juveniles in nursery estuaries and adult populations is an important tool for fisheries management. Otoliths of juvenile fish contain geochemical tags, which reflect the variation in estuarine elemental chemistry, and allow discrimination of their natal and/or nursery estuaries. These tags can be used to investigate connectivity patterns between juveniles and adults. However, inter-annual variability of geochemical tags may limit the accuracy of nursery origin determinations. Otolith elemental composition was used to assign a single cohort of 0-group sea bass Dicentrarchus labrax to their nursery estuary thus establishing an initial baseline for stocks in waters around Ireland. Using a standard LDFA model, high classification accuracies to nursery sites (80-88%) were obtained. Temporal stability of otolith geochemical tags was also investigated to assess if annual sampling is required for connectivity studies. Geochemical tag stability was found to be strongly estuary dependent.

Energy-Constrained Recharge, Assimilation, and Fractional Crystallization (EC-RAχFC): A Visual Basic computer code for calculating trace element and isotope variations of open-system magmatic systems

NASA Astrophysics Data System (ADS)

Bohrson, Wendy A.; Spera, Frank J.

2007-11-01

Volcanic and plutonic rocks provide abundant evidence for complex processes that occur in magma storage and transport systems. The fingerprint of these processes, which include fractional crystallization, assimilation, and magma recharge, is captured in petrologic and geochemical characteristics of suites of cogenetic rocks. Quantitatively evaluating the relative contributions of each process requires integration of mass, species, and energy constraints, applied in a self-consistent way. The energy-constrained model Energy-Constrained Recharge, Assimilation, and Fractional Crystallization (EC-RaχFC) tracks the trace element and isotopic evolution of a magmatic system (melt + solids) undergoing simultaneous fractional crystallization, recharge, and assimilation. Mass, thermal, and compositional (trace element and isotope) output is provided for melt in the magma body, cumulates, enclaves, and anatectic (i.e., country rock) melt. Theory of the EC computational method has been presented by Spera and Bohrson (2001, 2002, 2004), and applications to natural systems have been elucidated by Bohrson and Spera (2001, 2003) and Fowler et al. (2004). The purpose of this contribution is to make the final version of the EC-RAχFC computer code available and to provide instructions for code implementation, description of input and output parameters, and estimates of typical values for some input parameters. A brief discussion highlights measures by which the user may evaluate the quality of the output and also provides some guidelines for implementing nonlinear productivity functions. The EC-RAχFC computer code is written in Visual Basic, the programming language of Excel. The code therefore launches in Excel and is compatible with both PC and MAC platforms. The code is available on the authors' Web sites http://magma.geol.ucsb.edu/and http://www.geology.cwu.edu/ecrafc) as well as in the auxiliary material.

Lead Isotopes in Olivine-Phyric Shergottite Tissint: Implications for the Geochemical Evolution of the Shergottite Source Mantle

NASA Technical Reports Server (NTRS)

Moriwaki, R.; Usui, T.; Simon, J. I.; Jones, J. H.; Yokoyama, T.

2015-01-01

Geochemically-depleted shergottites are basaltic rocks derived from a martian mantle source reservoir. Geochemical evolution of the martian mantle has been investigated mainly based on the Rb-Sr, Sm-Nd, and Lu-Hf isotope systematics of the shergottites [1]. Although potentially informative, U-Th- Pb isotope systematics have been limited because of difficulties in interpreting the analyses of depleted meteorite samples that are more susceptible to the effects of near-surface processes and terrestrial contamination. This study conducts a 5-step sequential acid leaching experiment of the first witnessed fall of the geochemically-depleted olivinephyric shergottite Tissint to minimize the effect of low temperature distrubence. Trace element analyses of the Tissint acid residue (mostly pyroxene) indicate that Pb isotope compositions of the residue do not contain either a martian surface or terrestrial component, but represent the Tissint magma source [2]. The residue has relatively unradiogenic initial Pb isotopic compositions (e.g., 206Pb/204Pb = 10.8136) that fall within the Pb isotope space of other geochemically-depleted shergottites. An initial µ-value (238U/204Pb = 1.5) of Tissint at the time of crystallization (472 Ma [3]) is similar to a time-integrated mu- value (1.72 at 472 Ma) of the Tissint source mantle calculated based on the two-stage mantle evolution model [1]. On the other hand, the other geochemically-depleted shergottites (e.g., QUE 94201 [4]) have initial µ-values of their parental magmas distinctly lower than those of their modeled source mantle. These results suggest that only Tissint potentially reflects the geochemical signature of the shergottite mantle source that originated from cumulates of the martian magma ocean

Calculation of individual isotope equilibrium constants for implementation in geochemical models

USGS Publications Warehouse

Thorstenson, Donald C.; Parkhurst, David L.

2002-01-01

Theory is derived from the work of Urey to calculate equilibrium constants commonly used in geochemical equilibrium and reaction-transport models for reactions of individual isotopic species. Urey showed that equilibrium constants of isotope exchange reactions for molecules that contain two or more atoms of the same element in equivalent positions are related to isotope fractionation factors by , where is n the number of atoms exchanged. This relation is extended to include species containing multiple isotopes, for example and , and to include the effects of nonideality. The equilibrium constants of the isotope exchange reactions provide a basis for calculating the individual isotope equilibrium constants for the geochemical modeling reactions. The temperature dependence of the individual isotope equilibrium constants can be calculated from the temperature dependence of the fractionation factors. Equilibrium constants are calculated for all species that can be formed from and selected species containing , in the molecules and the ion pairs with where the subscripts g, aq, l, and s refer to gas, aqueous, liquid, and solid, respectively. These equilibrium constants are used in the geochemical model PHREEQC to produce an equilibrium and reaction-transport model that includes these isotopic species. Methods are presented for calculation of the individual isotope equilibrium constants for the asymmetric bicarbonate ion. An example calculates the equilibrium of multiple isotopes among multiple species and phases.

A sorption model for alkalis in cement-based materials - Correlations with solubility and electrokinetic properties

NASA Astrophysics Data System (ADS)

Henocq, Pierre

2017-06-01

In cement-based materials, radionuclide uptake is mainly controlled by calcium silicate hydrates (C-S-H). This work presents an approach for defining a unique set of parameters of a surface complexation model describing the sorption behavior of alkali ions on the C-S-H surface. Alkali sorption processes are modeled using the CD-MUSIC function integrated in the Phreeqc V.3.0.6 geochemical code. Parameterization of the model was performed based on (1) retention, (2) zeta potential, and (3) solubility experimental data from the literature. This paper shows an application of this model to sodium ions. It was shown that retention, i.e. surface interactions, and solubility are closely related, and a consistent sorption model for radionuclides in cement-based materials requires a coupled surface interaction/chemical equilibrium model. In case of C-S-H with low calcium-to-silicon ratios, sorption of sodium ions on the C-S-H surface strongly influences the chemical equilibrium of the C-S-H + NaCl system by significantly increasing the aqueous calcium concentration. The close relationship between sorption and chemical equilibrium was successfully illustrated by modeling the effect of the solid-to-liquid ratio on the calcium content in solution in the case of C-S-H + NaCl systems.

Introducing 3D U-statistic method for separating anomaly from background in exploration geochemical data with associated software development

NASA Astrophysics Data System (ADS)

Ghannadpour, Seyyed Saeed; Hezarkhani, Ardeshir

2016-03-01

The U-statistic method is one of the most important structural methods to separate the anomaly from the background. It considers the location of samples and carries out the statistical analysis of the data without judging from a geochemical point of view and tries to separate subpopulations and determine anomalous areas. In the present study, to use U-statistic method in three-dimensional (3D) condition, U-statistic is applied on the grade of two ideal test examples, by considering sample Z values (elevation). So far, this is the first time that this method has been applied on a 3D condition. To evaluate the performance of 3D U-statistic method and in order to compare U-statistic with one non-structural method, the method of threshold assessment based on median and standard deviation (MSD method) is applied on the two example tests. Results show that the samples indicated by U-statistic method as anomalous are more regular and involve less dispersion than those indicated by the MSD method. So that, according to the location of anomalous samples, denser areas of them can be determined as promising zones. Moreover, results show that at a threshold of U = 0, the total error of misclassification for U-statistic method is much smaller than the total error of criteria of bar {x}+n× s. Finally, 3D model of two test examples for separating anomaly from background using 3D U-statistic method is provided. The source code for a software program, which was developed in the MATLAB programming language in order to perform the calculations of the 3D U-spatial statistic method, is additionally provided. This software is compatible with all the geochemical varieties and can be used in similar exploration projects.

Surrogate model approach for improving the performance of reactive transport simulations

NASA Astrophysics Data System (ADS)

Jatnieks, Janis; De Lucia, Marco; Sips, Mike; Dransch, Doris

2016-04-01

Reactive transport models can serve a large number of important geoscientific applications involving underground resources in industry and scientific research. It is common for simulation of reactive transport to consist of at least two coupled simulation models. First is a hydrodynamics simulator that is responsible for simulating the flow of groundwaters and transport of solutes. Hydrodynamics simulators are well established technology and can be very efficient. When hydrodynamics simulations are performed without coupled geochemistry, their spatial geometries can span millions of elements even when running on desktop workstations. Second is a geochemical simulation model that is coupled to the hydrodynamics simulator. Geochemical simulation models are much more computationally costly. This is a problem that makes reactive transport simulations spanning millions of spatial elements very difficult to achieve. To address this problem we propose to replace the coupled geochemical simulation model with a surrogate model. A surrogate is a statistical model created to include only the necessary subset of simulator complexity for a particular scenario. To demonstrate the viability of such an approach we tested it on a popular reactive transport benchmark problem that involves 1D Calcite transport. This is a published benchmark problem (Kolditz, 2012) for simulation models and for this reason we use it to test the surrogate model approach. To do this we tried a number of statistical models available through the caret and DiceEval packages for R, to be used as surrogate models. These were trained on randomly sampled subset of the input-output data from the geochemical simulation model used in the original reactive transport simulation. For validation we use the surrogate model to predict the simulator output using the part of sampled input data that was not used for training the statistical model. For this scenario we find that the multivariate adaptive regression splines (MARS) method provides the best trade-off between speed and accuracy. This proof-of-concept forms an essential step towards building an interactive visual analytics system to enable user-driven systematic creation of geochemical surrogate models. Such a system shall enable reactive transport simulations with unprecedented spatial and temporal detail to become possible. References: Kolditz, O., Görke, U.J., Shao, H. and Wang, W., 2012. Thermo-hydro-mechanical-chemical processes in porous media: benchmarks and examples (Vol. 86). Springer Science & Business Media.

Geochemistry of groundwater in the Beaver and Camas Creek drainage basins, eastern Idaho

USGS Publications Warehouse

Rattray, Gordon W.; Ginsbach, Michael L.

2014-01-01

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, is studying the fate and transport of waste solutes in the eastern Snake River Plain (ESRP) aquifer at the Idaho National Laboratory (INL) in eastern Idaho. This effort requires an understanding of the natural and anthropogenic geochemistry of groundwater at the INL and of the important physical and chemical processes controlling the geochemistry. In this study, the USGS applied geochemical modeling to investigate the geochemistry of groundwater in the Beaver and Camas Creek drainage basins, which provide groundwater recharge to the ESRP aquifer underlying the northeastern part of the INL. Data used in this study include petrology and mineralogy from 2 sediment and 3 rock samples, and water-quality analyses from 4 surface-water and 18 groundwater samples. The mineralogy of the sediment and rock samples was analyzed with X-ray diffraction, and the mineralogy and petrology of the rock samples were examined in thin sections. The water samples were analyzed for field parameters, major ions, silica, nutrients, dissolved organic carbon, trace elements, tritium, and the stable isotope ratios of hydrogen, oxygen, carbon, sulfur, and nitrogen. Groundwater geochemistry was influenced by reactions with rocks of the geologic terranes—carbonate rocks, rhyolite, basalt, evaporite deposits, and sediment comprised of all of these rocks. Agricultural practices near and south of Dubois and application of road anti-icing liquids on U.S. Interstate Highway 15 were likely sources of nitrate, chloride, calcium, and magnesium to groundwater. Groundwater geochemistry was successfully modeled in the alluvial aquifer in Camas Meadows and the ESRP fractured basalt aquifer using the geochemical modeling code PHREEQC. The primary geochemical processes appear to be precipitation or dissolution of calcite and dissolution of silicate minerals. Dissolution of evaporite minerals, associated with Pleistocene Lake Terreton, is an important contributor of solutes in the Mud Lake-Dubois area. Oxidation-reduction reactions are important influences on the chemistry of groundwater at Camas Meadows and the Camas National Wildlife Refuge. In addition, mixing of different groundwaters or surface water with groundwater appears to be an important physical process influencing groundwater geochemistry in much of the study area, and evaporation may be an important physical process influencing the groundwater geochemistry of the Camas National Wildlife Refuge. The mass-balance modeling results from this study provide an explanation of the natural geochemistry of groundwater in the ESRP aquifer northeast of the INL, and thus provide a starting point for evaluating the natural and anthropogenic geochemistry of groundwater at the INL.

Final Report: Molecular mechanisms and kinetics of microbial anaerobic nitrate-dependent U(IV) and Fe(II) oxidation

DOE Office of Scientific and Technical Information (OSTI.GOV)

O'Day, Peggy A.; Asta, Maria P.; Kanematsu, Masakazu

2015-02-27

In this project, we combined molecular genetic, spectroscopic, and microscopic techniques with kinetic and reactive transport studies to describe and quantify biotic and abiotic mechanisms underlying anaerobic, nitrate-dependent U(IV) and Fe(II) oxidation, which influences the long-term efficacy of in situ reductive immobilization of uranium at DOE sites. In these studies, Thiobacillus denitrificans, an autotrophic bacterium that catalyzes anaerobic U(IV) and Fe(II) oxidation, was used to examine coupled oxidation-reduction processes under either biotic (enzymatic) or abiotic conditions in batch and column experiments with biogenically produced UIVO2(s). Synthesis and quantitative analysis of coupled chemical and transport processes were done with the reactivemore » transport modeling code Crunchflow. Research focused on identifying the primary redox proteins that catalyze metal oxidation, environmental factors that influence protein expression, and molecular-scale geochemical factors that control the rates of biotic and abiotic oxidation.« less

Review of simulation techniques for Aquifer Thermal Energy Storage (ATES)

NASA Astrophysics Data System (ADS)

Mercer, J. W.; Faust, C. R.; Miller, W. J.; Pearson, F. J., Jr.

1981-03-01

The analysis of aquifer thermal energy storage (ATES) systems rely on the results from mathematical and geochemical models. Therefore, the state-of-the-art models relevant to ATES were reviewed and evaluated. These models describe important processes active in ATES including ground-water flow, heat transport (heat flow), solute transport (movement of contaminants), and geochemical reactions. In general, available models of the saturated ground-water environment are adequate to address most concerns associated with ATES; that is, design, operation, and environmental assessment. In those cases where models are not adequate, development should be preceded by efforts to identify significant physical phenomena and relate model parameters to measurable quantities.

Modeling biogechemical reactive transport in a fracture zone

DOE Office of Scientific and Technical Information (OSTI.GOV)

Molinero, Jorge; Samper, Javier; Yang, Chan Bing, and Zhang, Guoxiang

2005-01-14

A coupled model of groundwater flow, reactive solute transport and microbial processes for a fracture zone of the Aspo site at Sweden is presented. This is the model of the so-called Redox Zone Experiment aimed at evaluating the effects of tunnel construction on the geochemical conditions prevailing in a fracture granite. It is found that a model accounting for microbially-mediated geochemical processes is able to reproduce the unexpected measured increasing trends of dissolved sulfate and bicarbonate. The model is also useful for testing hypotheses regarding the role of microbial processes and evaluating the sensitivity of model results to changes inmore » biochemical parameters.« less

Strontium Isotopes and the Reconstruction of the Chaco Regional System: Evaluating Uncertainty with Bayesian Mixing Models

PubMed Central

Drake, Brandon Lee; Wills, Wirt H.; Hamilton, Marian I.; Dorshow, Wetherbee

2014-01-01

Strontium isotope sourcing has become a common and useful method for assigning sources to archaeological artifacts. In Chaco Canyon, an Ancestral Pueblo regional center in New Mexico, previous studies using these methods have suggested that significant portion of maize and wood originate in the Chuska Mountains region, 75 km to the East. In the present manuscript, these results were tested using both frequentist methods (to determine if geochemical sources can truly be differentiated) and Bayesian methods (to address uncertainty in geochemical source attribution). It was found that Chaco Canyon and the Chuska Mountain region are not easily distinguishable based on radiogenic strontium isotope values. The strontium profiles of many geochemical sources in the region overlap, making it difficult to definitively identify any one particular geochemical source for the canyon's pre-historic maize. Bayesian mixing models support the argument that some spruce and fir wood originated in the San Mateo Mountains, but that this cannot explain all 87Sr/86Sr values in Chaco timber. Overall radiogenic strontium isotope data do not clearly identify a single major geochemical source for maize, ponderosa, and most spruce/fir timber. As such, the degree to which Chaco Canyon relied upon outside support for both food and construction material is still ambiguous. PMID:24854352

IRON PRECIPITATION AND ARSENIC ATTENUATION - ASSESSMENT OF ARSENIC NATURAL ATTENUATION OF THE SUBSURFACE USING A GEOCHEMICAL MODEL (PHREEQC)

EPA Science Inventory

Laboratory experiments show that amorphous and poorly crystallized ferric iron hydroxides have much greater capacity to attenuate arsenic compared to clays and other aluminosilicate minerals. Studies (e.g., Lin and Qvarfort, 1996) showed that a sudden change in geochemical condit...

Improved Geothermometry Through Multivariate Reaction-path Modeling and Evaluation of Geomicrobiological Influences on Geochemical Temperature Indicators: Final Report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mattson, Earl; Smith, Robert; Fujita, Yoshiko

2015-03-01

The project was aimed at demonstrating that the geothermometric predictions can be improved through the application of multi-element reaction path modeling that accounts for lithologic and tectonic settings, while also accounting for biological influences on geochemical temperature indicators. The limited utilization of chemical signatures by individual traditional geothermometer in the development of reservoir temperature estimates may have been constraining their reliability for evaluation of potential geothermal resources. This project, however, was intended to build a geothermometry tool which can integrate multi-component reaction path modeling with process-optimization capability that can be applied to dilute, low-temperature water samples to consistently predict reservoirmore » temperature within ±30 °C. The project was also intended to evaluate the extent to which microbiological processes can modulate the geochemical signals in some thermal waters and influence the geothermometric predictions.« less

Interactive client side data visualization with d3.js

NASA Astrophysics Data System (ADS)

Rodzianko, A.; Versteeg, R.; Johnson, D. V.; Soltanian, M. R.; Versteeg, O. J.; Girouard, M.

2015-12-01

Geoscience data associated with near surface research and operational sites is increasingly voluminous and heterogeneous (both in terms of providers and data types - e.g. geochemical, hydrological, geophysical, modeling data, of varying spatiotemporal characteristics). Such data allows scientists to investigate fundamental hydrological and geochemical processes relevant to agriculture, water resources and climate change. For scientists to easily share, model and interpret such data requires novel tools with capabilities for interactive data visualization. Under sponsorship of the US Department of Energy, Subsurface Insights is developing the Predictive Assimilative Framework (PAF): a cloud based subsurface monitoring platform which can manage, process and visualize large heterogeneous datasets. Over the last year we transitioned our visualization method from a server side approach (in which images and animations were generated using Jfreechart and Visit) to a client side one that utilizes the D3 Javascript library. Datasets are retrieved using web service calls to the server, returned as JSON objects and visualized within the browser. Users can interactively explore primary and secondary datasets from various field locations. Our current capabilities include interactive data contouring and heterogeneous time series data visualization. While this approach is very powerful and not necessarily unique, special attention needs to be paid to latency and responsiveness issues as well as to issues as cross browser code compatibility so that users have an identical, fluid and frustration-free experience across different computational platforms. We gratefully acknowledge support from the US Department of Energy under SBIR Award DOE DE-SC0009732, the use of data from the Lawrence Berkeley National Laboratory (LBNL) Sustainable Systems SFA Rifle field site and collaboration with LBNL SFA scientists.

Management and assimilation of diverse, distributed watershed datasets

NASA Astrophysics Data System (ADS)

Varadharajan, C.; Faybishenko, B.; Versteeg, R.; Agarwal, D.; Hubbard, S. S.; Hendrix, V.

2016-12-01

The U.S. Department of Energy's (DOE) Watershed Function Scientific Focus Area (SFA) seeks to determine how perturbations to mountainous watersheds (e.g., floods, drought, early snowmelt) impact the downstream delivery of water, nutrients, carbon, and metals over seasonal to decadal timescales. We are building a software platform that enables integration of diverse and disparate field, laboratory, and simulation datasets, of various types including hydrological, geological, meteorological, geophysical, geochemical, ecological and genomic datasets across a range of spatial and temporal scales within the Rifle floodplain and the East River watershed, Colorado. We are using agile data management and assimilation approaches, to enable web-based integration of heterogeneous, multi-scale dataSensor-based observations of water-level, vadose zone and groundwater temperature, water quality, meteorology as well as biogeochemical analyses of soil and groundwater samples have been curated and archived in federated databases. Quality Assurance and Quality Control (QA/QC) are performed on priority datasets needed for on-going scientific analyses, and hydrological and geochemical modeling. Automated QA/QC methods are used to identify and flag issues in the datasets. Data integration is achieved via a brokering service that dynamically integrates data from distributed databases via web services, based on user queries. The integrated results are presented to users in a portal that enables intuitive search, interactive visualization and download of integrated datasets. The concepts, approaches and codes being used are shared across various data science components of various large DOE-funded projects such as the Watershed Function SFA, Next Generation Ecosystem Experiment (NGEE) Tropics, Ameriflux/FLUXNET, and Advanced Simulation Capability for Environmental Management (ASCEM), and together contribute towards DOE's cyberinfrastructure for data management and model-data integration.

Uranium resource assessment through statistical analysis of exploration geochemical and other data. Final report. [Codes EVAL, SURE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koch, G.S. Jr.; Howarth, R.J.; Schuenemeyer, J.H.

1981-02-01

We have developed a procedure that can help quadrangle evaluators to systematically summarize and use hydrogeochemical and stream sediment reconnaissance (HSSR) and occurrence data. Although we have not provided an independent estimate of uranium endowment, we have devised a methodology that will provide this independent estimate when additional calibration is done by enlarging the study area. Our statistical model for evaluation (system EVAL) ranks uranium endowment for each quadrangle. Because using this model requires experience in geology, statistics, and data analysis, we have also devised a simplified model, presented in the package SURE, a System for Uranium Resource Evaluation. Wemore » have developed and tested these models for the four quadrangles in southern Colorado that comprise the study area; to investigate their generality, the models should be applied to other quandrangles. Once they are calibrated with accepted uranium endowments for several well-known quadrangles, the models can be used to give independent estimates for less-known quadrangles. The point-oriented models structure the objective comparison of the quandrangles on the bases of: (1) Anomalies (a) derived from stream sediments, (b) derived from waters (stream, well, pond, etc.), (2) Geology (a) source rocks, as defined by the evaluator, (b) host rocks, as defined by the evaluator, and (3) Aerial radiometric anomalies.« less

Produced water re-injection in a non-fresh water aquifer with geochemical reaction, hydrodynamic molecular dispersion and adsorption kinetics controlling: model development and numerical simulation

NASA Astrophysics Data System (ADS)

Obe, Ibidapo; Fashanu, T. A.; Idialu, Peter O.; Akintola, Tope O.; Abhulimen, Kingsley E.

2017-06-01

An improved produced water reinjection (PWRI) model that incorporates filtration, geochemical reaction, molecular transport, and mass adsorption kinetics was developed to predict cake deposition and injectivity performance in hydrocarbon aquifers in Nigeria oil fields. Thus, the improved PWRI model considered contributions of geochemical reaction, adsorption kinetics, and hydrodynamic molecular dispersion mechanism to alter the injectivity and deposition of suspended solids on aquifer wall resulting in cake formation in pores during PWRI and transport of active constituents in hydrocarbon reservoirs. The injectivity decline and cake deposition for specific case studies of hydrocarbon aquifers in Nigeria oil fields were characterized with respect to its well geometry, lithology, and calibrations data and simulated in COMSOL multiphysics software environment. The PWRI model was validated by comparisons to assessments of previous field studies based on data and results supplied by operator and regulator. The results of simulation showed that PWRI performance was altered because of temporal variations and declinations of permeability, injectivity, and cake precipitation, which were observed to be dependent on active adsorption and geochemical reaction kinetics coupled with filtration scheme and molecular dispersion. From the observed results and findings, transition time t r to cake nucleation and growth were dependent on aquifer constituents, well capacity, filtration coefficients, particle-to-grain size ratio, water quality, and more importantly, particle-to-grain adsorption kinetics. Thus, the results showed that injectivity decline and permeability damage were direct contributions of geochemical reaction, hydrodynamic molecular diffusion, and adsorption kinetics to the internal filtration mechanism, which are largely dependent on the initial conditions of concentration of active constituents of produced water and aquifer capacity.

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

USGS Publications Warehouse

Charlton, Scott R.; Parkhurst, David L.

2011-01-01

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

USGS Publications Warehouse

Charlton, S.R.; Parkhurst, D.L.

2011-01-01

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server-for example, Excel??, Visual Basic??, Python, or MATLAB??. PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations. ?? 2011.

Updated Reference Model for Heat Generation in the Lithosphere

NASA Astrophysics Data System (ADS)

Wipperfurth, S. A.; Sramek, O.; Roskovec, B.; Mantovani, F.; McDonough, W. F.

2017-12-01

Models integrating geophysics and geochemistry allow for characterization of the Earth's heat budget and geochemical evolution. Global lithospheric geophysical models are now constrained by surface and body wave data and are classified into several unique tectonic types. Global lithospheric geochemical models have evolved from petrological characterization of layers to a combination of petrologic and seismic constraints. Because of these advances regarding our knowledge of the lithosphere, it is necessary to create an updated chemical and physical reference model. We are developing a global lithospheric reference model based on LITHO1.0 (segmented into 1°lon x 1°lat x 9-layers) and seismological-geochemical relationships. Uncertainty assignments and correlations are assessed for its physical attributes, including layer thickness, Vp and Vs, and density. This approach yields uncertainties for the masses of the crust and lithospheric mantle. Heat producing element abundances (HPE: U, Th, and K) are ascribed to each volume element. These chemical attributes are based upon the composition of subducting sediment (sediment layers), composition of surface rocks (upper crust), a combination of petrologic and seismic correlations (middle and lower crust), and a compilation of xenolith data (lithospheric mantle). The HPE abundances are correlated within each voxel, but not vertically between layers. Efforts to provide correlation of abundances horizontally between each voxel are discussed. These models are used further to critically evaluate the bulk lithosphere heat production in the continents and the oceans. Cross-checks between our model and results from: 1) heat flux (Artemieva, 2006; Davies, 2013; Cammarano and Guerri, 2017), 2) gravity (Reguzzoni and Sampietro, 2015), and 3) geochemical and petrological models (Rudnick and Gao, 2014; Hacker et al. 2015) are performed.

Long-term climate change and the geochemical cycle of carbon

NASA Technical Reports Server (NTRS)

Marshall, Hal G.; Walker, James C. G.; Kuhn, William R.

1988-01-01

The response of the coupled climate-geochemical system to changes in paleography is examined in terms of the biogeochemical carbon cycle. The simple, zonally averaged energy balance climate model combined with a geochemical carbon cycle model, which was developed to study climate changes, is described. The effects of latitudinal distributions of the continents on the carbon cycle are investigated, and the global silicate weathering rate as a function of latitude is measured. It is observed that a concentration of land area at high altitudes results in a high CO2 partial pressure and a high global average temperature, and for land at low latitudes a cold globe and ice are detected. It is noted that the CO2 greenhouse feedback effect is potentially strong and has a stabilizing effect on the climate system.

Conference on Early Mars: Geologic and Hydrologic Evolution, Physical and Chemical Environments, and the Implications for Life

NASA Technical Reports Server (NTRS)

Clifford, S. M. (Editor); Treiman, A. H. (Editor); Newsom, H. E. (Editor); Farmer, J. D. (Editor)

1997-01-01

Topics considered include: Geology alteration and life in an extreme environment; developing a chemical code to identify magnetic biominerals; effect of impacts on early Martin geologic evolution; spectroscopic identification of minerals in Hematite-bearing soils and sediments; exopaleontology and the search for a Fossil record on Mars; geochemical evolution of the crust of Mars; geological evolution of the early earth;solar-wind-induced erosion of the Mars atmosphere. Also included geological evolution of the crust of Mars.

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

PubMed Central

2011-01-01

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface. PMID:21943229

Mineralization of Basalts in the CO 2-H 2O-H 2S System

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schaef, Herbert T.; McGrail, B. Peter; Owen, Antionette T.

2013-05-10

Basalt samples representing five different formations were immersed in water equilibrated with supercritical carbon dioxide containing 1% hydrogen sulfide (H2S) at reservoir conditions (100 bar, 90°C) for up to 3.5 years. Surface coatings in the form of pyrite and metal cation substituted carbonates were identified as reaction products associated with all five basalts. In some cases, high pressure tests contained excess H2S, which produced the most corroded basalts and largest amount of secondary products. In comparison, tests containing limited amounts of H2S appeared least reacted with significantly less concentrations of reaction products. In all cases, pyrite appeared to precede carbonation,more » and in some instances, was observed in the absence of carbonation such as in cracks, fractures, and within the porous glassy mesostasis. Armoring reactions from pyrite surface coatings observed in earlier shorter duration tests were found to be temporary with carbonate mineralization observed with all the basalts tested in these long duration experiments. Geochemical simulations conducted with the geochemical code EQ3/6 accurately predicted early pyrite precipitation followed by formation of carbonates. Reactivity with H2S was correlated with measured Fe(II)/Fe(III) ratios in the basalts with more facile pyrite formation occurring with basalts containing more Fe(III) phases. These experimental and modeling results confirm potential for long term sequestration of acid gas mixtures in continental flood basalt formations.« less

BLT-EC (Breach, Leach and Transport-Equilibrium Chemistry) data input guide. A computer model for simulating release and coupled geochemical transport of contaminants from a subsurface disposal facility

DOE Office of Scientific and Technical Information (OSTI.GOV)

MacKinnon, R.J.; Sullivan, T.M.; Kinsey, R.R.

1997-05-01

The BLT-EC computer code has been developed, implemented, and tested. BLT-EC is a two-dimensional finite element computer code capable of simulating the time-dependent release and reactive transport of aqueous phase species in a subsurface soil system. BLT-EC contains models to simulate the processes (container degradation, waste-form performance, transport, chemical reactions, and radioactive production and decay) most relevant to estimating the release and transport of contaminants from a subsurface disposal system. Water flow is provided through tabular input or auxiliary files. Container degradation considers localized failure due to pitting corrosion and general failure due to uniform surface degradation processes. Waste-form performancemore » considers release to be limited by one of four mechanisms: rinse with partitioning, diffusion, uniform surface degradation, and solubility. Transport considers the processes of advection, dispersion, diffusion, chemical reaction, radioactive production and decay, and sources (waste form releases). Chemical reactions accounted for include complexation, sorption, dissolution-precipitation, oxidation-reduction, and ion exchange. Radioactive production and decay in the waste form is simulated. To improve the usefulness of BLT-EC, a pre-processor, ECIN, which assists in the creation of chemistry input files, and a post-processor, BLTPLOT, which provides a visual display of the data have been developed. BLT-EC also includes an extensive database of thermodynamic data that is also accessible to ECIN. This document reviews the models implemented in BLT-EC and serves as a guide to creating input files and applying BLT-EC.« less

Experimental Study of Cement - Sandstone/Shale - Brine - CO2 Interactions

PubMed Central

2011-01-01

Background Reactive-transport simulation is a tool that is being used to estimate long-term trapping of CO2, and wellbore and cap rock integrity for geologic CO2 storage. We reacted end member components of a heterolithic sandstone and shale unit that forms the upper section of the In Salah Gas Project carbon storage reservoir in Krechba, Algeria with supercritical CO2, brine, and with/without cement at reservoir conditions to develop experimentally constrained geochemical models for use in reactive transport simulations. Results We observe marked changes in solution composition when CO2 reacted with cement, sandstone, and shale components at reservoir conditions. The geochemical model for the reaction of sandstone and shale with CO2 and brine is a simple one in which albite, chlorite, illite and carbonate minerals partially dissolve and boehmite, smectite, and amorphous silica precipitate. The geochemical model for the wellbore environment is also fairly simple, in which alkaline cements and rock react with CO2-rich brines to form an Fe containing calcite, amorphous silica, smectite and boehmite or amorphous Al(OH)3. Conclusions Our research shows that relatively simple geochemical models can describe the dominant reactions that are likely to occur when CO2 is stored in deep saline aquifers sealed with overlying shale cap rocks, as well as the dominant reactions for cement carbonation at the wellbore interface. PMID:22078161

Evaluation of the transport matrix method for simulation of ocean biogeochemical tracers

NASA Astrophysics Data System (ADS)

Kvale, Karin F.; Khatiwala, Samar; Dietze, Heiner; Kriest, Iris; Oschlies, Andreas

2017-06-01

Conventional integration of Earth system and ocean models can accrue considerable computational expenses, particularly for marine biogeochemical applications. Offline numerical schemes in which only the biogeochemical tracers are time stepped and transported using a pre-computed circulation field can substantially reduce the burden and are thus an attractive alternative. One such scheme is the transport matrix method (TMM), which represents tracer transport as a sequence of sparse matrix-vector products that can be performed efficiently on distributed-memory computers. While the TMM has been used for a variety of geochemical and biogeochemical studies, to date the resulting solutions have not been comprehensively assessed against their online counterparts. Here, we present a detailed comparison of the two. It is based on simulations of the state-of-the-art biogeochemical sub-model embedded within the widely used coarse-resolution University of Victoria Earth System Climate Model (UVic ESCM). The default, non-linear advection scheme was first replaced with a linear, third-order upwind-biased advection scheme to satisfy the linearity requirement of the TMM. Transport matrices were extracted from an equilibrium run of the physical model and subsequently used to integrate the biogeochemical model offline to equilibrium. The identical biogeochemical model was also run online. Our simulations show that offline integration introduces some bias to biogeochemical quantities through the omission of the polar filtering used in UVic ESCM and in the offline application of time-dependent forcing fields, with high latitudes showing the largest differences with respect to the online model. Differences in other regions and in the seasonality of nutrients and phytoplankton distributions are found to be relatively minor, giving confidence that the TMM is a reliable tool for offline integration of complex biogeochemical models. Moreover, while UVic ESCM is a serial code, the TMM can be run on a parallel machine with no change to the underlying biogeochemical code, thus providing orders of magnitude speed-up over the online model.

High-temperature life without photosynthesis as a model for Mars

NASA Technical Reports Server (NTRS)

Shock, E. L.

1997-01-01

Discoveries in biology and developments in geochemistry over the past two decades have lead to a radical revision of concepts relating to the upper temperature at which life thrives, the genetic relationships among all life on Earth, links between organic and inorganic compounds in geologic processes, and the geochemical supply of metabolic energy. It is now apparent that given a source of geochemical energy, in the form of a mixture of compounds that is far from thermodynamic equilibrium, microorganisms can take advantage of the energy and thrive without the need for photosynthesis as a means of primary productivity. This means that life can exist in the subsurface of a planet such as Mars without necessarily exhibiting a surface expression. Theoretical calculations quantify the geochemically provided metabolic energy available to hyperthermophilic organisms in submarine hydrothermal systems on the Earth, and help to explain the enormous biological productivity of these systems. Efforts to place these models in the context of the early Earth reveal that substantial geochemical energy would have been available and that organic synthesis would have been thermodynamically favored as hydrothermal fluids mix with seawater.

Geochemical study of groundwater at Sandia National Laboratories/New Mexico and Kirtland Air Force Base

DOE Office of Scientific and Technical Information (OSTI.GOV)

NONE

The U.S. Department of Energy (DOE) Grand Junction Projects Office (GJPO) and its contractor, Rust Geotech, support the Kirtland Area Office by assisting Sandia National Laboratories/New Mexico (Sandia/NM) with remedial action, remedial design, and technical support of its Environmental Restoration Program. To aid in determining groundwater origins and flow paths, the GJPO was tasked to provide interpretation of groundwater geochemical data. The purpose of this investigation was to describe and analyze the groundwater geochemistry of the Sandia/NM Kirtland Air Force Base (KAFB). Interpretations of groundwater origins are made by using these data and the results of {open_quotes}mass balance{close_quotes} and {open_quotes}reactionmore » path{close_quote} modeling. Additional maps and plots were compiled to more fully comprehend the geochemical distributions. A more complete set of these data representations are provided in the appendices. Previous interpretations of groundwater-flow paths that were based on well-head, geologic, and geochemical data are presented in various reports and were used as the basis for developing the models presented in this investigation.« less

Porosity Development in a Coastal Setting: A Reactive Transport Model to Assess the Influence of Heterogeneity of Hydrological, Geochemical and Lithological Conditions

NASA Astrophysics Data System (ADS)

Maqueda, A.; Renard, P.; Cornaton, F. J.

2014-12-01

Coastal karst networks are formed by mineral dissolution, mainly calcite, in the freshwater-saltwater mixing zone. The problem has been approached first by studying the kinetics of calcite dissolution and then coupling ion-pairing software with flow and mass transport models. Porosity development models require high computational power. A workaround to reduce computational complexity is to assume the calcite dissolution reaction is relatively fast, thus equilibrium chemistry can be used to model it (Sanford & Konikow, 1989). Later developments allowed the full coupling of kinetics and transport in a model. However kinetics effects of calcite dissolution were found negligible under the single set of assumed hydrological and geochemical boundary conditions. A model is implemented with the coupling of FeFlow software as the flow & transport module and PHREEQC4FEFLOW (Wissmeier, 2013) ion-pairing module. The model is used to assess the influence of heterogeneities in hydrological, geochemical and lithological boundary conditions on porosity evolution. The hydrologic conditions present in the karst aquifer of Quintana Roo coast in Mexico are used as a guide for generating inputs for simulations.

Characteristics and habitat of deep vs. shallow slow slip events

NASA Astrophysics Data System (ADS)

Wipperfurth, S. A.; Sramek, O.; Roskovec, B.; Mantovani, F.; McDonough, W. F.

2016-12-01

Models integrating geophysics and geochemistry allow for characterization of the Earth's heat budget and geochemical evolution. Global lithospheric geophysical models are now constrained by surface and body wave data and are classified into several unique tectonic types. Global lithospheric geochemical models have evolved from petrological characterization of layers to a combination of petrologic and seismic constraints. Because of these advances regarding our knowledge of the lithosphere, it is necessary to create an updated chemical and physical reference model. We are developing a global lithospheric reference model based on LITHO1.0 (segmented into 1°lon x 1°lat x 9-layers) and seismological-geochemical relationships. Uncertainty assignments and correlations are assessed for its physical attributes, including layer thickness, Vp and Vs, and density. This approach yields uncertainties for the masses of the crust and lithospheric mantle. Heat producing element abundances (HPE: U, Th, and K) are ascribed to each volume element. These chemical attributes are based upon the composition of subducting sediment (sediment layers), composition of surface rocks (upper crust), a combination of petrologic and seismic correlations (middle and lower crust), and a compilation of xenolith data (lithospheric mantle). The HPE abundances are correlated within each voxel, but not vertically between layers. Efforts to provide correlation of abundances horizontally between each voxel are discussed. These models are used further to critically evaluate the bulk lithosphere heat production in the continents and the oceans. Cross-checks between our model and results from: 1) heat flux (Artemieva, 2006; Davies, 2013; Cammarano and Guerri, 2017), 2) gravity (Reguzzoni and Sampietro, 2015), and 3) geochemical and petrological models (Rudnick and Gao, 2014; Hacker et al. 2015) are performed.

Matlab Geochemistry: An open source geochemistry solver based on MRST

NASA Astrophysics Data System (ADS)

McNeece, C. J.; Raynaud, X.; Nilsen, H.; Hesse, M. A.

2017-12-01

The study of geological systems often requires the solution of complex geochemical relations. To address this need we present an open source geochemical solver based on the Matlab Reservoir Simulation Toolbox (MRST) developed by SINTEF. The implementation supports non-isothermal multicomponent aqueous complexation, surface complexation, ion exchange, and dissolution/precipitation reactions. The suite of tools available in MRST allows for rapid model development, in particular the incorporation of geochemical calculations into transport simulations of multiple phases, complex domain geometry and geomechanics. Different numerical schemes and additional physics can be easily incorporated into the existing tools through the object-oriented framework employed by MRST. The solver leverages the automatic differentiation tools available in MRST to solve arbitrarily complex geochemical systems with any choice of species or element concentration as input. Four mathematical approaches enable the solver to be quite robust: 1) the choice of chemical elements as the basis components makes all entries in the composition matrix positive thus preserving convexity, 2) a log variable transformation is used which transfers the nonlinearity to the convex composition matrix, 3) a priori bounds on variables are calculated from the structure of the problem, constraining Netwon's path and 4) an initial guess is calculated implicitly by sequentially adding model complexity. As a benchmark we compare the model to experimental and semi-analytic solutions of the coupled salinity-acidity transport system. Together with the reservoir simulation capabilities of MRST the solver offers a promising tool for geochemical simulations in reservoir domains for applications in a diversity of fields from enhanced oil recovery to radionuclide storage.

Fractionation and risk assessment of Fe and Mn in surface sediments from coastal sites of Sonora, Mexico (Gulf of California).

PubMed

Jara-Marini, Martín E; García-Camarena, Raúl; Gómez-Álvarez, Agustín; García-Rico, Leticia

2015-07-01

The aim of this study was to evaluate Fe and Mn distribution in geochemical fractions of the surface sediment of four oyster culture sites in the Sonora coast, Mexico. A selective fractionation scheme to obtain five fractions was adapted for the microwave system. Surface sediments were analyzed for carbonates, organic matter contents, and Fe and Mn in geochemical fractions. The bulk concentrations of Fe ranged from 10,506 to 21,918 mg/kg (dry weight, dry wt), and the bulk concentrations of Mn ranged from 185.1 to 315.9 mg/kg (dry wt) in sediments, which was low and considered as non-polluted in all of the sites. The fractionation study indicated that the major geochemical phases for the metals were the residual, as well as the Fe and Mn oxide fractions. The concentrations of metals in the geochemical fractions had the following order: residual > Fe and Mn oxides > organic matter > carbonates > interchangeable. Most of the Fe and Mn were linked to the residual fraction. Among non-residual fractions, high percentages of Fe and Mn were linked to Fe and Mn oxides. The enrichment factors (EFs) for the two metals were similar in the four studied coasts, and the levels of Fe and Mn are interpreted as non-enrichment (EF < 1) because the metals concentrations were within the baseline concentrations. According to the environmental risk assessment codes, Fe and Mn posed no risk and low risk, respectively. Although the concentrations of Fe and Mn were linked to the residual fraction, the levels in non-residual fractions may significantly result in the transference of other metals, depending on several physico-chemical and biological factors.

Mineral dissolution and precipitation during CO 2 injection at the Frio-I Brine Pilot: Geochemical modeling and uncertainty analysis

DOE PAGES

Ilgen, A. G.; Cygan, R. T.

2015-12-07

During the Frio-I Brine Pilot CO 2 injection experiment in 2004, distinct geochemical changes in response to the injection of 1600 tons of CO 2 were recorded in samples collected from the monitoring well. Previous geochemical modeling studies have considered dissolution of calcite and iron oxyhydroxides, or release of adsorbed iron, as the most likely sources of the increased ion concentrations. We explore in this modeling study possible alternative sources of the increasing calcium and iron, based on the data from the detailed petrographic characterization of the Upper Frio Formation “C”. Particularly, we evaluate whether dissolution of pyrite andmore » oligoclase (anorthite component) can account for the observed geochemical changes. Due to kinetic limitations, dissolution of pyrite and anorthite cannot account for the increased iron and calcium concentrations on the time scale of the field test (10 days). However, dissolution of these minerals is contributing to carbonate and clay mineral precipitation on the longer time scales (1000 years). The one-dimensional reactive transport model predicts carbonate minerals, dolomite and ankerite, as well as clay minerals kaolinite, nontronite and montmorillonite, will precipitate in the Frio Formation “C” sandstone as the system progresses towards chemical equilibrium during a 1000-year period. Cumulative uncertainties associated with using different thermodynamic databases, activity correction models (Pitzer vs. B-dot), and extrapolating to reservoir temperature, are manifested in the difference in the predicted mineral phases. Furthermore, these models are consistent with regards to the total volume of mineral precipitation and porosity values which are predicted to within 0.002%.« less

Advances in the hydrogeochemistry and microbiology of acid mine waters

USGS Publications Warehouse

Nordstrom, D. Kirk

2000-01-01

The last decade has witnessed a plethora of research related to the hydrogeochemistry and microbiology of acid mine waters and associated tailings and waste-rock waters. Numerous books, reviews, technical papers, and proceedings have been published that examine the complex bio-geochemical process of sulfide mineral oxidation, develop and apply geochemical models to site characterization, and characterize the microbial ecology of these environments. This review summarizes many of these recent works, and provides references for those investigating this field. Comparisons of measured versus calculated Eh and measured versus calculated pH for water samples from several field sites demonstrate the reliability of some current geochemical models for aqueous speciation and mass balances. Geochemical models are not, however, used to predict accurately time-dependent processes but to improve our understanding of these systems and to constrain possible processes that contribute to actual or potential water quality issues. Microbiological studies are demonstrating that there is much we have yet to learn about the types of different microorganisms and their function and ecology in mine-waste environments. A broad diversity of green algae, bacteria, archaea, yeasts, and fungi are encountered in acid mine waters, and a better understanding of their ecology and function may potentially enhance remediation possibilities as well as our understanding of the evolution of life.

Approach for environmental baseline water sampling

USGS Publications Warehouse

Smith, K.S.

2011-01-01

Samples collected during the exploration phase of mining represent baseline conditions at the site. As such, they can be very important in forecasting potential environmental impacts should mining proceed, and can become measurements against which future changes are compared. Constituents in stream water draining mined and mineralized areas tend to be geochemically, spatially, and temporally variable, which presents challenges in collecting both exploration and baseline water-quality samples. Because short-term (daily) variations can complicate long-term trends, it is important to consider recent findings concerning geochemical variability of stream-water constituents at short-term timescales in designing sampling plans. Also, adequate water-quality information is key to forecasting potential ecological impacts from mining. Therefore, it is useful to collect baseline water samples adequate tor geochemical and toxicological modeling. This requires complete chemical analyses of dissolved constituents that include major and minor chemical elements as well as physicochemical properties (including pH, specific conductance, dissolved oxygen) and dissolved organic carbon. Applying chemical-equilibrium and appropriate toxicological models to water-quality information leads to an understanding of the speciation, transport, sequestration, bioavailability, and aquatic toxicity of potential contaminants. Insights gained from geochemical and toxicological modeling of water-quality data can be used to design appropriate mitigation and for economic planning for future mining activities.

Merging Hydrologic, Geochemical, and Geophysical Approaches to Understand the Regolith Architecture of a Deeply Weathered Piedmont Critical Zone

NASA Astrophysics Data System (ADS)

Cosans, C.; Moore, J.; Harman, C. J.

2017-12-01

Located in the deeply weathered Piedmont in Maryland, Pond Branch has a rich legacy of hydrological and geochemical research dating back to the first geochemical mass balance study published in 1970. More recently, geophysical investigations including seismic and electrical resistivity tomography have characterized the subsurface at Pond Branch and contributed to new hypotheses about critical zone evolution. Heterogeneity in electrical resistivity in the shallow subsurface may suggest disparate flow paths for recharge, with some regions with low hydraulic conductivity generating perched flow, while other hillslope sections recharge to the much deeper regolith boundary. These shallow and deep flow paths are hypothesized to be somewhat hydrologically and chemically connected, with the spatially and temporally discontinuous connections resulting in different hydraulic responses to recharge and different concentrations of weathering solutes. To test this hypothesis, we combined modeling and field approaches. We modeled weathering solutes along the hypothesized flow paths using PFLOTRAN. We measured hydrologic gradients in the hillslopes and riparian zone using piezometer water levels. We collected geochemical data including major ions and silica. Weathering solute concentrations were measured directly in the precipitation, hillslope springs, and the riparian zone for comparison to modeled concentration values. End member mixing methods were used to determine contributions of precipitation, hillslopes, and riparian zone to the stream. Combining geophysical, geochemical, and hydrological methods may offer insights into the source of stream water and controls on chemical weathering. Previous hypotheses that Piedmont critical zone architecture results from a balance of erosion, soil, and weathering front advance rates cannot account for the inverted regolith structure observed through seismic investigations at Pond Branch. Recent alternative hypotheses including weathering along tectonically-induced fractures and weathering front advance have been proposed, but additional data are needed to test them. Developing a thorough, nuanced understanding of the geochemical and hydrological behavior of Pond Branch may help test and refine hypotheses for Piedmont critical zone evolution.

Modelling of clay diagenesis using a combined approach of crystalchemistry and thermochemistry: a case study in the smectite illitization.

NASA Astrophysics Data System (ADS)

Geloni, Claudio; Previde Massara, Elisabetta; Di Paola, Eleonora; Ortenzi, Andrea; Gherardi, Fabrizio; Blanc, Philippe

2017-04-01

Diagenetic transformations occurring in clayey and arenaceous sediments is investigated in a number of hydrocarbon reservoirs with an integrated approach that combines mineralogical analysis, crystalchemistry, estimation of thermochemical parameters of clay minerals, and geochemical modelling. Because of the extremely variable crystalchemistry of clays, especially in the smectite - illite compositional range, the estimation of thermochemical parameters of site-specific clay-rich rocks is crucial to investigate water-rock equilibria and to predict mineralogical evolutionary patterns at the clay-sandstone interface. The task of estimating the thermochemical properties of clay minerals and predicting diagenetic reactions in natural reservoirs is accomplished through the implementation of an informatized, procedure (IP) that consists of: (i) laboratory analysis of smectite, illite and mixed layers (I/S) for the determination of their textural characteristics and chemical composition; (ii) estimation of the thermodynamic and structural parameters (enthalpy, entropy, and free energy of formation, thermal capacity, molar volume, molar weight) with a MS Excel tool (XLS) specifically developed at the French Bureau of Geological and Mining Researches (BRGM); (iii) usage of the SUPCRT (Johnson et al., 1992) software package (thereinafter, SSP) to derive log K values to be incorporated in thermodynamic databases of the standard geochemical codes; (iv) check of the consistency of the stability domains calculated with these log K values with relevant predominance diagrams; (v) final application of geochemical and reactive transport models to investigate the reactive mechanisms under different thermal conditions (40-150°C). All the simulations consider pore waters having roughly the same chemical composition of reservoir pore waters, and are performed with The Geochemist Workbench (Bethke and Yeakel, 2015), PHREEQC (Parkhurst, 1999) and TOUGHREACT (Xu, 2006). The overall procedure benefits from: (i) (minor) improvements of the I/O structure of the SSP; (ii) the development of a suite of python scripts to automate the steps needed to augment the thermodynamic database by integrating the external information provided by potential users with the XLS tool and the SSP; (iii) the creation of specific outputs to allow for more convenient handling and inspection of computed parameters of the thermodynamic database. A case study focused on non-isothermal smectite-illite transformation is presented to show the capability of our numerical models to account for clay compaction under 1D geometry conditions. This model considers fluid flow driven by the compaction of a clay layer, and chemistry-fluid flow mutual feedback with the underlying sandstone during the advancement of the diagenesis. Due to this complex interaction, as a result of the smectite-illite transformation in the clays, significant quartz cementation affects the sandstone adjacent to the compacting clay.

Lead Isotope Compositions of Acid Residues from Olivine-Phyric Shergottite Tissint: Implications for Heterogeneous Shergottite Source Reservoirs

NASA Technical Reports Server (NTRS)

Moriwaki, R.; Usui, T.; Yokoyama, T.; Simon, J. I.; Jones, J. H.

2015-01-01

Geochemical studies of shergottites suggest that their parental magmas reflect mixtures between at least two distinct geochemical source reservoirs, producing correlations between radiogenic isotope compositions and trace element abundances. These correlations have been interpreted as indicating the presence of a reduced, incompatible element- depleted reservoir and an oxidized, incompatible- element-enriched reservoir. The former is clearly a depleted mantle source, but there is ongoing debate regarding the origin of the enriched reservoir. Two contrasting models have been proposed regarding the location and mixing process of the two geochemical source reservoirs: (1) assimilation of oxidized crust by mantle derived, reduced magmas, or (2) mixing of two distinct mantle reservoirs during melting. The former requires the ancient Martian crust to be the enriched source (crustal assimilation), whereas the latter requires isolation of a long-lived enriched mantle domain that probably originated from residual melts formed during solidification of a magma ocean (heterogeneous mantle model). This study conducts Pb isotope and trace element concentration analyses of sequential acid-leaching fractions (leachates and the final residues) from the geochemically depleted olivine-phyric shergottite Tissint. The results suggest that the Tissint magma is not isotopically uniform and sampled at least two geochemical source reservoirs, implying that either crustal assimilation or magma mixing would have played a role in the Tissint petrogenesis.

A combined geodynamical-geochemical modelling approach to investigating the Lu-Hf isotopic evolution of the terrestrial mantle and crust

NASA Astrophysics Data System (ADS)

Jones, R.; Van Keken, P. E.; Hauri, E.; Vervoort, J. D.; Ballentine, C. J.

2017-12-01

The chemical and isotopic evolution of the Earth's mantle is largely influenced by the formation of oceanic and continental crust at spreading ridges and through arc volcanism, and the subsequent recycling of this crust back into the mantle via subduction. In this study we use a combined geodynamical-geochemical modelling approach to investigate the Lu-Hf isotopic evolution of the terrestrial mantle and crust. We utilise the geodynamic mantle convection model developed by Brandenburg et al., 2008. This model satisfies the geophysical constraints of oceanic heat flow and average plate velocities, as well as geochemical observations such as 40Ar in the atmosphere. It has also been shown to reproduce the observed geochemical distributions in multiple isotope systems (U-Th-Pb, Rb-Sr, Sm-Nd, and Re-Os) that define the DMM, HIMU and EM1 mantle endmembers. We go on to extend this application to investigate the Lu-Hf isotope system, specifically in combination with Sm-Nd. The model has been updated to include a self-consistent reorganisation of the plates with regions of up-/down-wellings. The model is initiated at 4.55 Ga, assumes continental crust is produced from 4 Ga and that a transition from `dry' to `wet' subduction occurs at 3 Ga. The results of the geodynamic model suggest that the ƐHf composition and evolution of the upper mantle can be generated through the extraction and recycling of oceanic crust, which creates an enriched and radiogenic reservoir at the core-mantle boundary. The formation of continental crust, which is extracted at each time-step from the oceanic crust to imitate subduction zone processes, and the recycling of this continental crust as sediments, plays a lesser role. Depending on the selected partition coefficients DMM, FOZO and HIMU mantle endmember compositions are also produced via the simple extraction and recycling of oceanic crust. The formation of continental crust produces spread in the ƐNd vs. ƐHf array and extends the model values into the HIMU region of the terrestrial array. We go on to use this geodynamic-geochemical model to investigate different models of continental growth, by observing the effects on the coupled crustal-mantle reservoirs. Brandenburg, J.P., Hauri, E.H., van Keken, P.E., Ballentine, C.J., 2008. Earth and Planetary Science Letters 276, 1-13.

Investigating ebullition in a sand column using dissolved gas analysis and reactive transport modeling

USGS Publications Warehouse

Amos, Richard T.; Mayer, K. Ulrich

2006-01-01

Ebullition of gas bubbles through saturated sediments can enhance the migration of gases through the subsurface, affect the rate of biogeochemical processes, and potentially enhance the emission of important greenhouse gases to the atmosphere. To better understand the parameters controlling ebullition, methanogenic conditions were produced in a column experiment and ebullition through the column was monitored and quantified through dissolved gas analysis and reactive transport modeling. Dissolved gas analysis showed rapid transport of CH4 vertically through the column at rates several times faster than the bromide tracer and the more soluble gas CO2, indicating that ebullition was the main transport mechanism for CH4. An empirically derived formulation describing ebullition was integrated into the reactive transport code MIN3P allowing this process to be investigated on the REV scale in a complex geochemical framework. The simulations provided insights into the parameters controlling ebullition and show that, over the duration of the experiment, 36% of the CH4 and 19% of the CO2 produced were transported to the top of the column through ebullition.

Coupling ANIMO and MT3DMS for 3D regional-scale modeling of nutrient transport in soil and groundwater

NASA Astrophysics Data System (ADS)

Janssen, G.; Del Val Alonso, L.; Groenendijk, P.; Griffioen, J.

2012-12-01

We developed an on-line coupling between the 1D/quasi-2D nutrient transport model ANIMO and the 3D groundwater transport model code MT3DMS. ANIMO is a detailed, process-oriented model code for the simulation of nitrate leaching to groundwater, N- and P-loads on surface waters and emissions of greenhouse gasses. It is the leading nutrient fate and transport code in the Netherlands where it is used primarily for the evaluation of fertilization related legislation. In addition, the code is applied frequently in international research projects. MT3DMS is probably the most commonly used groundwater solute transport package worldwide. The on-line model coupling ANIMO-MT3DMS combines the state-of-the-art descriptions of the biogeochemical cycles in ANIMO with the advantages of using a 3D approach for the transport through the saturated domain. These advantages include accounting for regional lateral transport, considering groundwater-surface water interactions more explicitly, and the possibility of using MODFLOW to obtain the flow fields. An additional merit of the on-line coupling concept is that it preserves feedbacks between the saturated and unsaturated zone. We tested ANIMO-MT3DMS by simulating nutrient transport for the period 1970-2007 in a Dutch agricultural polder catchment covering an area of 118 km2. The transient groundwater flow field had a temporal resolution of one day and was calculated with MODFLOW-MetaSWAP. The horizontal resolution of the model grid was 100x100m and consisted of 25 layers of varying thickness. To keep computation times manageable, we prepared MT3DMS for parallel computing, which in itself is a relevant development for a large community of groundwater transport modelers. For the parameterization of the soil, we applied a standard classification approach, representing the area by 60 units with unique combinations of soil type, land use and geohydrological setting. For the geochemical parameterization of the deeper subsurface, however, we applied a novel geostatistical technique, which allocates reactivity parameters to the grid cells by sampling from these parameters' cumulative frequency distribution (CDF) functions. These CDF functions are derived for each relevant geohydrological unit present in the model domain, from datasets of groundwater and sediment analyses. The nutrient loads on the surface water system and the nutrient concentrations in groundwater, simulated by the transport model, are in fair agreement with field measurements. The experience with the test model constitutes a proof-of-concept, justifying further developments towards application of ANIMO-MT3DMS in actual regional decision-making processes.

Reactive alteration of Mt. Simon sandstone during CO2-rich brine injection: A coupled experimental and modeling study

NASA Astrophysics Data System (ADS)

Dávila Ordoñez, M. G.; Zahasky, C.; Crandall, D.; Druhan, J. L.

2017-12-01

Thus far, one million metric tons of CO2 have been injected into the lower Mt. Simon formation as part of the Decatur CO2 Capture and Storage Project. Micro-seismic events were observed within the CO2 plume both during and after pressurization associated with the primary injection. The Mt. Simon reservoir rock consists of 76.5 wt.% quartz, 2.1 wt.% calcite, 17.3 wt.% K-feldspar, 1.1 wt.% chlorite, 0.7 wt.% illite and lesser extents of siderite, kaolinite, dolomite and marcasite, and is thus anticipated to become geochemically altered by exposure to acidified CO2-rich brine. However, the extent to which the geochemical reactivity contributes to structural weakening is unknown. To explore relationships between the principle geochemical reactions, evolution of fluid transport properties and physical alteration, we performed a series of flow-through experiments using Mt. Simon core (5 cm diameter, ranging from 4.3 - 8.6 cm length) and fluids representative of acidified reservoir brine. Experiments were operated under P = 1450 bar, Pconfining = 1900 - 3000 bar and T = 53 ºC conditions, and flow rates varied from 0.08 to 5.00 mL h-1 over a period of 166 h. A 2D reactive transport code (Crunch-Tope) was used to simulate these experiments, constrained by measured time series aqueous concentrations of Ca, Mg, S, Si, K and Fe and pH during the CO2-rich brine interaction. The model domain was divided into 30 nodes in x at a spacing of 0.12 cm, and 40 nodes in y at a spacing of 0.22 cm, and initial permeability measured for the core was specified and allowed to evolve over the course of the simulation using measured flow rate as a constraint. All relevant kinetic and thermodynamic reaction parameters were obtained from the literature. Solute time series from both experiments and simulations indicated that the acidified brine introduced continuously into the column promoted dissolution of K-feldspar, chloride, illite, pyrite and calcite, and the precipitation of Ca-, Fe- and Si -bearing secondary phases, resulting in a net porosity increase at the inlet. Despite this opening of the inlet pore space, permeability decreased over the length of the column (kfinal/kinitial = 0.76), thus altering local resistance to fluid phase pressure gradients.

Early Earth plume-lid tectonics: A high-resolution 3D numerical modelling approach

NASA Astrophysics Data System (ADS)

Fischer, R.; Gerya, T.

2016-10-01

Geological-geochemical evidence point towards higher mantle potential temperature and a different type of tectonics (global plume-lid tectonics) in the early Earth (>3.2 Ga) compared to the present day (global plate tectonics). In order to investigate tectono-magmatic processes associated with plume-lid tectonics and crustal growth under hotter mantle temperature conditions, we conduct a series of 3D high-resolution magmatic-thermomechanical models with the finite-difference code I3ELVIS. No external plate tectonic forces are applied to isolate 3D effects of various plume-lithosphere and crust-mantle interactions. Results of the numerical experiments show two distinct phases in coupled crust-mantle evolution: (1) a longer (80-100 Myr) and relatively quiet 'growth phase' which is marked by growth of crust and lithosphere, followed by (2) a short (∼20 Myr) and catastrophic 'removal phase', where unstable parts of the crust and mantle lithosphere are removed by eclogitic dripping and later delamination. This modelling suggests that the early Earth plume-lid tectonic regime followed a pattern of episodic growth and removal also called episodic overturn with a periodicity of ∼100 Myr.

Testing and Resilience of the Impact Origin of the Moon

NASA Technical Reports Server (NTRS)

Righter, K.; Canup, R. M.

2016-01-01

The leading hypothesis for the origin of the Moon is the giant impact model, which grew out of the post-Apollo science community. The hypothesis was able to explain the high E-M system angular momentum, the small lunar core, and consistent with the idea that the early Moon melted substantially. The standard hypothesis requires that the Moon be made entirely from the impactor, strangely at odds with the nearly identical oxygen isotopic composition of the Earth and Moon, compositions that might be expected to be different if Moon came from a distinct impactor. Subsequent geochemical research has highlighted the similarity of both geochemical and isotopic composition of the Earth and Moon, and measured small but significant amounts of volatiles in lunar glassy materials, both of which are seemingly at odds with the standard giant impact model. Here we focus on key geochemical measurements and spacecraft observations that have prompted a healthy re-evaluation of the giant impact model, provide an overview of physical models that are either newly proposed or slightly revised from previous ideas, to explain the new datasets.

Assessment of groundwater quality: a fusion of geochemical and geophysical information via Bayesian neural networks.

PubMed

Maiti, Saumen; Erram, V C; Gupta, Gautam; Tiwari, Ram Krishna; Kulkarni, U D; Sangpal, R R

2013-04-01

Deplorable quality of groundwater arising from saltwater intrusion, natural leaching and anthropogenic activities is one of the major concerns for the society. Assessment of groundwater quality is, therefore, a primary objective of scientific research. Here, we propose an artificial neural network-based method set in a Bayesian neural network (BNN) framework and employ it to assess groundwater quality. The approach is based on analyzing 36 water samples and inverting up to 85 Schlumberger vertical electrical sounding data. We constructed a priori model by suitably parameterizing geochemical and geophysical data collected from the western part of India. The posterior model (post-inversion) was estimated using the BNN learning procedure and global hybrid Monte Carlo/Markov Chain Monte Carlo optimization scheme. By suitable parameterization of geochemical and geophysical parameters, we simulated 1,500 training samples, out of which 50 % samples were used for training and remaining 50 % were used for validation and testing. We show that the trained model is able to classify validation and test samples with 85 % and 80 % accuracy respectively. Based on cross-correlation analysis and Gibb's diagram of geochemical attributes, the groundwater qualities of the study area were classified into following three categories: "Very good", "Good", and "Unsuitable". The BNN model-based results suggest that groundwater quality falls mostly in the range of "Good" to "Very good" except for some places near the Arabian Sea. The new modeling results powered by uncertainty and statistical analyses would provide useful constrain, which could be utilized in monitoring and assessment of the groundwater quality.

Multivariate analysis of ATR-FTIR spectra for assessment of oil shale organic geochemical properties

USGS Publications Warehouse

Washburn, Kathryn E.; Birdwell, Justin E.

2013-01-01

In this study, attenuated total reflectance (ATR) Fourier transform infrared spectroscopy (FTIR) was coupled with partial least squares regression (PLSR) analysis to relate spectral data to parameters from total organic carbon (TOC) analysis and programmed pyrolysis to assess the feasibility of developing predictive models to estimate important organic geochemical parameters. The advantage of ATR-FTIR over traditional analytical methods is that source rocks can be analyzed in the laboratory or field in seconds, facilitating more rapid and thorough screening than would be possible using other tools. ATR-FTIR spectra, TOC concentrations and Rock–Eval parameters were measured for a set of oil shales from deposits around the world and several pyrolyzed oil shale samples. PLSR models were developed to predict the measured geochemical parameters from infrared spectra. Application of the resulting models to a set of test spectra excluded from the training set generated accurate predictions of TOC and most Rock–Eval parameters. The critical region of the infrared spectrum for assessing S1, S2, Hydrogen Index and TOC consisted of aliphatic organic moieties (2800–3000 cm−1) and the models generated a better correlation with measured values of TOC and S2 than did integrated aliphatic peak areas. The results suggest that combining ATR-FTIR with PLSR is a reliable approach for estimating useful geochemical parameters of oil shales that is faster and requires less sample preparation than current screening methods.

PhreeqcRM: A reaction module for transport simulators based on the geochemical model PHREEQC

USGS Publications Warehouse

Parkhurst, David L.; Wissmeier, Laurin

2015-01-01

PhreeqcRM is a geochemical reaction module designed specifically to perform equilibrium and kinetic reaction calculations for reactive transport simulators that use an operator-splitting approach. The basic function of the reaction module is to take component concentrations from the model cells of the transport simulator, run geochemical reactions, and return updated component concentrations to the transport simulator. If multicomponent diffusion is modeled (e.g., Nernst–Planck equation), then aqueous species concentrations can be used instead of component concentrations. The reaction capabilities are a complete implementation of the reaction capabilities of PHREEQC. In each cell, the reaction module maintains the composition of all of the reactants, which may include minerals, exchangers, surface complexers, gas phases, solid solutions, and user-defined kinetic reactants.PhreeqcRM assigns initial and boundary conditions for model cells based on standard PHREEQC input definitions (files or strings) of chemical compositions of solutions and reactants. Additional PhreeqcRM capabilities include methods to eliminate reaction calculations for inactive parts of a model domain, transfer concentrations and other model properties, and retrieve selected results. The module demonstrates good scalability for parallel processing by using multiprocessing with MPI (message passing interface) on distributed memory systems, and limited scalability using multithreading with OpenMP on shared memory systems. PhreeqcRM is written in C++, but interfaces allow methods to be called from C or Fortran. By using the PhreeqcRM reaction module, an existing multicomponent transport simulator can be extended to simulate a wide range of geochemical reactions. Results of the implementation of PhreeqcRM as the reaction engine for transport simulators PHAST and FEFLOW are shown by using an analytical solution and the reactive transport benchmark of MoMaS.

PHREEQCI; a graphical user interface for the geochemical computer program PHREEQC

USGS Publications Warehouse

Charlton, Scott R.; Macklin, Clifford L.; Parkhurst, David L.

1997-01-01

PhreeqcI is a Windows-based graphical user interface for the geochemical computer program PHREEQC. PhreeqcI provides the capability to generate and edit input data files, run simulations, and view text files containing simulation results, all within the framework of a single interface. PHREEQC is a multipurpose geochemical program that can perform speciation, inverse, reaction-path, and 1D advective reaction-transport modeling. Interactive access to all of the capabilities of PHREEQC is available with PhreeqcI. The interface is written in Visual Basic and will run on personal computers under the Windows(3.1), Windows95, and WindowsNT operating systems.

A graphical method to evaluate predominant geochemical processes occurring in groundwater systems for radiocarbon dating

USGS Publications Warehouse

Han, Liang-Feng; Plummer, Niel; Aggarwal, Pradeep

2012-01-01

A graphical method is described for identifying geochemical reactions needed in the interpretation of radiocarbon age in groundwater systems. Graphs are constructed by plotting the measured 14C, δ13C, and concentration of dissolved inorganic carbon and are interpreted according to specific criteria to recognize water samples that are consistent with a wide range of processes, including geochemical reactions, carbon isotopic exchange, 14C decay, and mixing of waters. The graphs are used to provide a qualitative estimate of radiocarbon age, to deduce the hydrochemical complexity of a groundwater system, and to compare samples from different groundwater systems. Graphs of chemical and isotopic data from a series of previously-published groundwater studies are used to demonstrate the utility of the approach. Ultimately, the information derived from the graphs is used to improve geochemical models for adjustment of radiocarbon ages in groundwater systems.

GIS Methodic and New Database for Magmatic Rocks. Application for Atlantic Oceanic Magmatism.

NASA Astrophysics Data System (ADS)

Asavin, A. M.

2001-12-01

There are several geochemical Databases in INTERNET available now. There one of the main peculiarities of stored geochemical information is geographical coordinates of each samples in those Databases. As rule the software of this Database use spatial information only for users interface search procedures. In the other side, GIS-software (Geographical Information System software),for example ARC/INFO software which using for creation and analyzing special geological, geochemical and geophysical e-map, have been deeply involved with geographical coordinates for of samples. We join peculiarities GIS systems and relational geochemical Database from special software. Our geochemical information system created in Vernadsky Geological State Museum and institute of Geochemistry and Analytical Chemistry from Moscow. Now we tested system with data of geochemistry oceanic rock from Atlantic and Pacific oceans, about 10000 chemical analysis. GIS information content consist from e-map covers Wold Globes. Parts of these maps are Atlantic ocean covers gravica map (with grid 2''), oceanic bottom hot stream, altimeteric maps, seismic activity, tectonic map and geological map. Combination of this information content makes possible created new geochemical maps and combination of spatial analysis and numerical geochemical modeling of volcanic process in ocean segment. Now we tested information system on thick client technology. Interface between GIS system Arc/View and Database resides in special multiply SQL-queries sequence. The result of the above gueries were simple DBF-file with geographical coordinates. This file act at the instant of creation geochemical and other special e-map from oceanic region. We used more complex method for geophysical data. From ARC\\View we created grid cover for polygon spatial geophysical information.

A quantitative speciation model for the adsorption of organic pollutants on activated carbon.

PubMed

Grivé, M; García, D; Domènech, C; Richard, L; Rojo, I; Martínez, X; Rovira, M

2013-01-01

Granular activated carbon (GAC) is commonly used as adsorbent in water treatment plants given its high capacity for retaining organic pollutants in aqueous phase. The current knowledge on GAC behaviour is essentially empirical, and no quantitative description of the chemical relationships between GAC surface groups and pollutants has been proposed. In this paper, we describe a quantitative model for the adsorption of atrazine onto GAC surface. The model is based on results of potentiometric titrations and three types of adsorption experiments which have been carried out in order to determine the nature and distribution of the functional groups on the GAC surface, and evaluate the adsorption characteristics of GAC towards atrazine. Potentiometric titrations have indicated the existence of at least two different families of chemical groups on the GAC surface, including phenolic- and benzoic-type surface groups. Adsorption experiments with atrazine have been satisfactorily modelled with the geochemical code PhreeqC, assuming that atrazine is sorbed onto the GAC surface in equilibrium (log Ks = 5.1 ± 0.5). Independent thermodynamic calculations suggest a possible adsorption of atrazine on a benzoic derivative. The present work opens a new approach for improving the adsorption capabilities of GAC towards organic pollutants by modifying its chemical properties.

Leach test of cladding removal waste grout using Hanford groundwater

DOE Office of Scientific and Technical Information (OSTI.GOV)

Serne, R.J.; Martin, W.J.; Legore, V.L.

1995-09-01

This report describes laboratory experiments performed during 1986-1990 designed to produce empirical leach rate data for cladding removal waste (CRW) grout. At the completion of the laboratory work, funding was not available for report completion, and only now during final grout closeout activities is the report published. The leach rates serve as inputs to computer codes used in assessing the potential risk from the migration of waste species from disposed grout. This report discusses chemical analyses conducted on samples of CRW grout, and the results of geochemical computer code calculations that help identify mechanisms involved in the leaching process. Themore » semi-infinite solid diffusion model was selected as the most representative model for describing leaching of grouts. The use of this model with empirically derived leach constants yields conservative predictions of waste release rates, provided no significant changes occur in the grout leach processes over long time periods. The test methods included three types of leach tests--the American Nuclear Society (ANS) 16.1 intermittent solution exchange test, a static leach test, and a once-through flow column test. The synthetic CRW used in the tests was prepared in five batches using simulated liquid waste spiked with several radionuclides: iodine ({sup 125}I), carbon ({sup 14}C), technetium ({sup 99}Tc), cesium ({sup 137}Cs), strontium ({sup 85}Sr), americium ({sup 241}Am), and plutonium ({sup 238}Pu). The grout was formed by mixing the simulated liquid waste with dry blend containing Type I and Type II Portland cement, class F fly ash, Indian Red Pottery clay, and calcium hydroxide. The mixture was allowed to set and cure at room temperature in closed containers for at least 46 days before it was tested.« less

ab initio MD simulations of geomaterials with ~1000 atoms

NASA Astrophysics Data System (ADS)

Martin, G. B.; Kirtman, B.; Spera, F. J.

2009-12-01

In the last two decades, ab initio studies of materials using Density Functional Theory (DFT) have increased exponentially in popularity. DFT codes are now used routinely to simulate properties of geomaterials--mainly silicates and geochemically important metals such as Fe. These materials are ubiquitous in the Earth’s mantle and core and in terrestrial exoplanets. Because of computational limitations, most First Principles Molecular Dynamics (FPMD) calculations are done on systems of only ~100 atoms for a few picoseconds. While this approach can be useful for calculating physical quantities related to crystal structure, vibrational frequency, and other lattice-scale properties (especially in crystals), it is statistically marginal for duplicating physical properties of the liquid state like transport and structure. In MD simulations in the NEV ensemble, temperature (T), and pressure (P) fluctuations scale as N-1/2; small particle number (N) systems are therefore characterized by greater statistical state point location uncertainty than large N systems. Previous studies have used codes such as VASP where CPU time increases with N2, making calculations with N much greater than 100 impractical. SIESTA (Soler, et al. 2002) is a DFT code that enables electronic structure and MD computations on larger systems (N~103) by making some approximations, such as localized numerical orbitals, that would be useful in modeling some properties of geomaterials. Here we test the applicability of SIESTA to simulate geosilicates, both hydrous and anhydrous, in the solid and liquid state. We have used SIESTA for lattice calculations of brucite, Mg(OH)2, that compare very well to experiment and calculations using CRYSTAL, another DFT code. Good agreement between more classical DFT calculations and SIESTA is needed to justify study of geosilicates using SIESTA across a range of pressures and temperatures relevant to the Earth’s interior. Thus, it is useful to adjust parameters in SIESTA in accordance with calculations from CRYSTAL as a check on feasibility. Results are reported here that suggest SIESTA may indeed be useful to model silicate liquids at very high T and P.

A coupled THC model of the FEBEX in situ test with bentonite swelling and chemical and thermal osmosis.

PubMed

Zheng, Liange; Samper, Javier; Montenegro, Luis

2011-09-25

The performance assessment of a geological repository for radioactive waste requires quantifying the geochemical evolution of the bentonite engineered barrier. This barrier will be exposed to coupled thermal (T), hydrodynamic (H), mechanical (M) and chemical (C) processes. This paper presents a coupled THC model of the FEBEX (Full-scale Engineered Barrier EXperiment) in situ test which accounts for bentonite swelling and chemical and thermal osmosis. Model results attest the relevance of thermal osmosis and bentonite swelling for the geochemical evolution of the bentonite barrier while chemical osmosis is found to be almost irrelevant. The model has been tested with data collected after the dismantling of heater 1 of the in situ test. The model reproduces reasonably well the measured temperature, relative humidity, water content and inferred geochemical data. However, it fails to mimic the solute concentrations at the heater-bentonite and bentonite-granite interfaces because the model does not account for the volume change of bentonite, the CO(2)(g) degassing and the transport of vapor from the bentonite into the granite. The inferred HCO(3)(-) and pH data cannot be explained solely by solute transport, calcite dissolution and protonation/deprotonation by surface complexation, suggesting that such data may be affected also by other reactions. Published by Elsevier B.V.

The oceanic islands - Azores. [geological, geophysical and geochemical features

NASA Technical Reports Server (NTRS)

Ridley, W. I.; Watkins, N. D.; Macfarlane, D. J.

1974-01-01

A presentation is made of the known geological, geophysical, and geochemical data on the Azores. The regional setting of the islands is described; under the geological heading, surface geology and petrochemistry are discussed; and paleomagnetism, marine magnetic surveys, gravity, seismology, and heat flow are treated in the geophysics category. A model for the origin of the Azores is constructed on the basis of these observations.

Benchmark Problems of the Geothermal Technologies Office Code Comparison Study

DOE Office of Scientific and Technical Information (OSTI.GOV)

White, Mark D.; Podgorney, Robert; Kelkar, Sharad M.

A diverse suite of numerical simulators is currently being applied to predict or understand the performance of enhanced geothermal systems (EGS). To build confidence and identify critical development needs for these analytical tools, the United States Department of Energy, Geothermal Technologies Office has sponsored a Code Comparison Study (GTO-CCS), with participants from universities, industry, and national laboratories. A principal objective for the study was to create a community forum for improvement and verification of numerical simulators for EGS modeling. Teams participating in the study were those representing U.S. national laboratories, universities, and industries, and each team brought unique numerical simulationmore » capabilities to bear on the problems. Two classes of problems were developed during the study, benchmark problems and challenge problems. The benchmark problems were structured to test the ability of the collection of numerical simulators to solve various combinations of coupled thermal, hydrologic, geomechanical, and geochemical processes. This class of problems was strictly defined in terms of properties, driving forces, initial conditions, and boundary conditions. Study participants submitted solutions to problems for which their simulation tools were deemed capable or nearly capable. Some participating codes were originally developed for EGS applications whereas some others were designed for different applications but can simulate processes similar to those in EGS. Solution submissions from both were encouraged. In some cases, participants made small incremental changes to their numerical simulation codes to address specific elements of the problem, and in other cases participants submitted solutions with existing simulation tools, acknowledging the limitations of the code. The challenge problems were based on the enhanced geothermal systems research conducted at Fenton Hill, near Los Alamos, New Mexico, between 1974 and 1995. The problems involved two phases of research, stimulation, development, and circulation in two separate reservoirs. The challenge problems had specific questions to be answered via numerical simulation in three topical areas: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Whereas the benchmark class of problems were designed to test capabilities for modeling coupled processes under strictly specified conditions, the stated objective for the challenge class of problems was to demonstrate what new understanding of the Fenton Hill experiments could be realized via the application of modern numerical simulation tools by recognized expert practitioners.« less

Integrating surrogate models into subsurface simulation framework allows computation of complex reactive transport scenarios

NASA Astrophysics Data System (ADS)

De Lucia, Marco; Kempka, Thomas; Jatnieks, Janis; Kühn, Michael

2017-04-01

Reactive transport simulations - where geochemical reactions are coupled with hydrodynamic transport of reactants - are extremely time consuming and suffer from significant numerical issues. Given the high uncertainties inherently associated with the geochemical models, which also constitute the major computational bottleneck, such requirements may seem inappropriate and probably constitute the main limitation for their wide application. A promising way to ease and speed-up such coupled simulations is achievable employing statistical surrogates instead of "full-physics" geochemical models [1]. Data-driven surrogates are reduced models obtained on a set of pre-calculated "full physics" simulations, capturing their principal features while being extremely fast to compute. Model reduction of course comes at price of a precision loss; however, this appears justified in presence of large uncertainties regarding the parametrization of geochemical processes. This contribution illustrates the integration of surrogates into the flexible simulation framework currently being developed by the authors' research group [2]. The high level language of choice for obtaining and dealing with surrogate models is R, which profits from state-of-the-art methods for statistical analysis of large simulations ensembles. A stand-alone advective mass transport module was furthermore developed in order to add such capability to any multiphase finite volume hydrodynamic simulator within the simulation framework. We present 2D and 3D case studies benchmarking the performance of surrogates and "full physics" chemistry in scenarios pertaining the assessment of geological subsurface utilization. [1] Jatnieks, J., De Lucia, M., Dransch, D., Sips, M.: "Data-driven surrogate model approach for improving the performance of reactive transport simulations.", Energy Procedia 97, 2016, p. 447-453. [2] Kempka, T., Nakaten, B., De Lucia, M., Nakaten, N., Otto, C., Pohl, M., Chabab [Tillner], E., Kühn, M.: "Flexible Simulation Framework to Couple Processes in Complex 3D Models for Subsurface Utilization Assessment.", Energy Procedia, 97, 2016 p. 494-501.

Interactive Visual Analytics Approch for Exploration of Geochemical Model Simulations with Different Parameter Sets

NASA Astrophysics Data System (ADS)

Jatnieks, Janis; De Lucia, Marco; Sips, Mike; Dransch, Doris

2015-04-01

Many geoscience applications can benefit from testing many combinations of input parameters for geochemical simulation models. It is, however, a challenge to screen the input and output data from the model to identify the significant relationships between input parameters and output variables. For addressing this problem we propose a Visual Analytics approach that has been developed in an ongoing collaboration between computer science and geoscience researchers. Our Visual Analytics approach uses visualization methods of hierarchical horizontal axis, multi-factor stacked bar charts and interactive semi-automated filtering for input and output data together with automatic sensitivity analysis. This guides the users towards significant relationships. We implement our approach as an interactive data exploration tool. It is designed with flexibility in mind, so that a diverse set of tasks such as inverse modeling, sensitivity analysis and model parameter refinement can be supported. Here we demonstrate the capabilities of our approach by two examples for gas storage applications. For the first example our Visual Analytics approach enabled the analyst to observe how the element concentrations change around previously established baselines in response to thousands of different combinations of mineral phases. This supported combinatorial inverse modeling for interpreting observations about the chemical composition of the formation fluids at the Ketzin pilot site for CO2 storage. The results indicate that, within the experimental error range, the formation fluid cannot be considered at local thermodynamical equilibrium with the mineral assemblage of the reservoir rock. This is a valuable insight from the predictive geochemical modeling for the Ketzin site. For the second example our approach supports sensitivity analysis for a reaction involving the reductive dissolution of pyrite with formation of pyrrothite in presence of gaseous hydrogen. We determine that this reaction is thermodynamically favorable under a broad range of conditions. This includes low temperatures and absence of microbial catalysators. Our approach has potential for use in other applications that involve exploration of relationships in geochemical simulation model data.

Geochemical modeling of magma mixing and magma reservoir volumes during early episodes of Kīlauea Volcano's Pu`u `Ō`ō eruption

NASA Astrophysics Data System (ADS)

Shamberger, Patrick J.; Garcia, Michael O.

2007-02-01

Geochemical modeling of magma mixing allows for evaluation of volumes of magma storage reservoirs and magma plumbing configurations. A new analytical expression is derived for a simple two-component box-mixing model describing the proportions of mixing components in erupted lavas as a function of time. Four versions of this model are applied to a mixing trend spanning episodes 3 31 of Kilauea Volcano’s Puu Oo eruption, each testing different constraints on magma reservoir input and output fluxes. Unknown parameters (e.g., magma reservoir influx rate, initial reservoir volume) are optimized for each model using a non-linear least squares technique to fit model trends to geochemical time-series data. The modeled mixing trend closely reproduces the observed compositional trend. The two models that match measured lava effusion rates have constant magma input and output fluxes and suggest a large pre-mixing magma reservoir (46±2 and 49±1 million m3), with little or no volume change over time. This volume is much larger than a previous estimate for the shallow, dike-shaped magma reservoir under the Puu Oo vent, which grew from ˜3 to ˜10 12 million m3. These volumetric differences are interpreted as indicating that mixing occurred first in a larger, deeper reservoir before the magma was injected into the overlying smaller reservoir.

Predicting Biological Information Flow in a Model Oxygen Minimum Zone

NASA Astrophysics Data System (ADS)

Louca, S.; Hawley, A. K.; Katsev, S.; Beltran, M. T.; Bhatia, M. P.; Michiels, C.; Capelle, D.; Lavik, G.; Doebeli, M.; Crowe, S.; Hallam, S. J.

2016-02-01

Microbial activity drives marine biochemical fluxes and nutrient cycling at global scales. Geochemical measurements as well as molecular techniques such as metagenomics, metatranscriptomics and metaproteomics provide great insight into microbial activity. However, an integration of molecular and geochemical data into mechanistic biogeochemical models is still lacking. Recent work suggests that microbial metabolic pathways are, at the ecosystem level, strongly shaped by stoichiometric and energetic constraints. Hence, models rooted in fluxes of matter and energy may yield a holistic understanding of biogeochemistry. Furthermore, such pathway-centric models would allow a direct consolidation with meta'omic data. Here we present a pathway-centric biogeochemical model for the seasonal oxygen minimum zone in Saanich Inlet, a fjord off the coast of Vancouver Island. The model considers key dissimilatory nitrogen and sulfur fluxes, as well as the population dynamics of the genes that mediate them. By assuming a direct translation of biocatalyzed energy fluxes to biosynthesis rates, we make predictions about the distribution and activity of the corresponding genes. A comparison of the model to molecular measurements indicates that the model explains observed DNA, RNA, protein and cell depth profiles. This suggests that microbial activity in marine ecosystems such as oxygen minimum zones is well described by DNA abundance, which, in conjunction with geochemical constraints, determines pathway expression and process rates. Our work further demonstrates how meta'omic data can be mechanistically linked to environmental redox conditions and biogeochemical processes.

Modeling background radiation using geochemical data: A case study in and around Cameron, Arizona.

PubMed

Marsac, Kara E; Burnley, Pamela C; Adcock, Christopher T; Haber, Daniel A; Malchow, Russell L; Hausrath, Elisabeth M

2016-12-01

This study compares high resolution forward models of natural gamma-ray background with that measured by high resolution aerial gamma-ray surveys. The ability to predict variations in natural background radiation levels should prove useful for those engaged in measuring anthropogenic contributions to background radiation for the purpose of emergency response and homeland security operations. The forward models are based on geologic maps and remote sensing multi-spectral imagery combined with two different sources of data: 1) bedrock geochemical data (uranium, potassium and thorium concentrations) collected from national databases, the scientific literature and private companies, and 2) the low spatial resolution NURE (National Uranium Resource Evaluation) aerial gamma-ray survey. The study area near Cameron, Arizona, is located in an arid region with minimal vegetation and, due to the presence of abandoned uranium mines, was the subject of a previous high resolution gamma-ray survey. We found that, in general, geologic map units form a good basis for predicting the geographic distribution of the gamma-ray background. Predictions of background gamma-radiation levels based on bedrock geochemical analyses were not as successful as those based on the NURE aerial survey data sorted by geologic unit. The less successful result of the bedrock geochemical model is most likely due to a number of factors including the need to take into account the evolution of soil geochemistry during chemical weathering and the influence of aeolian addition. Refinements to the forward models were made using ASTER visualizations to create subunits of similar exposure rate within the Chinle Formation, which contains multiple lithologies and by grouping alluvial units by drainage basin rather than age. Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.

Modeling background radiation using geochemical data: A case study in and around Cameron, Arizona

DOE Office of Scientific and Technical Information (OSTI.GOV)

Marsac, Kara E.; Burnley, Pamela C.; Adcock, Christopher T.

Here, this study compares high-resolution forward models of natural gamma-ray background with that measured by high resolution aerial gamma-ray surveys. The ability to predict variations in natural background radiation levels should prove useful for those engaged in measuring anthropogenic contributions to background radiation for the purpose of emergency response and homeland security operations. The forward models are based on geologic maps and remote sensing multi-spectral imagery combined with two different sources of data: 1) bedrock geochemical data (uranium, potassium and thorium concentrations) collected from national databases, the scientific literature and private companies, and 2) the low spatial resolution NURE (Nationalmore » Uranium Resource Evaluation) aerial gamma-ray survey. The study area near Cameron, Arizona, is located in an arid region with minimal vegetation and, due to the presence of abandoned uranium mines, was the subject of a previous high resolution gamma-ray survey. We found that, in general, geologic map units form a good basis for predicting the geographic distribution of the gamma-ray background. Predictions of background gamma-radiation levels based on bedrock geochemical analyses were not as successful as those based on the NURE aerial survey data sorted by geologic unit. The less successful result of the bedrock geochemical model is most likely due to a number of factors including the need to take into account the evolution of soil geochemistry during chemical weathering and the influence of aeolian addition. Refinements to the forward models were made using ASTER visualizations to create subunits of similar exposure rate within the Chinle Formation, which contains multiple lithologies and by grouping alluvial units by drainage basin rather than age.« less

Modeling background radiation using geochemical data: A case study in and around Cameron, Arizona

DOE PAGES

Marsac, Kara E.; Burnley, Pamela C.; Adcock, Christopher T.; ...

2016-09-16

Here, this study compares high-resolution forward models of natural gamma-ray background with that measured by high resolution aerial gamma-ray surveys. The ability to predict variations in natural background radiation levels should prove useful for those engaged in measuring anthropogenic contributions to background radiation for the purpose of emergency response and homeland security operations. The forward models are based on geologic maps and remote sensing multi-spectral imagery combined with two different sources of data: 1) bedrock geochemical data (uranium, potassium and thorium concentrations) collected from national databases, the scientific literature and private companies, and 2) the low spatial resolution NURE (Nationalmore » Uranium Resource Evaluation) aerial gamma-ray survey. The study area near Cameron, Arizona, is located in an arid region with minimal vegetation and, due to the presence of abandoned uranium mines, was the subject of a previous high resolution gamma-ray survey. We found that, in general, geologic map units form a good basis for predicting the geographic distribution of the gamma-ray background. Predictions of background gamma-radiation levels based on bedrock geochemical analyses were not as successful as those based on the NURE aerial survey data sorted by geologic unit. The less successful result of the bedrock geochemical model is most likely due to a number of factors including the need to take into account the evolution of soil geochemistry during chemical weathering and the influence of aeolian addition. Refinements to the forward models were made using ASTER visualizations to create subunits of similar exposure rate within the Chinle Formation, which contains multiple lithologies and by grouping alluvial units by drainage basin rather than age.« less

How to build stable geochemical reservoirs on Mars?

NASA Astrophysics Data System (ADS)

Plesa, Ana-Catalina; Tosi, Nicola; Breuer, Doris

2014-05-01

To explain the complex thermo-chemical processes needed for the formation of distinct and stable geochemical reservoirs early in the thermo-chemical evolution of Mars, most geochemical studies argue that fractional crystallization of a global magma ocean may reproduce the isotopic characteristic of the SNCs [1, 2]. However, geodynamical models show that such scenario is difficult to reconcile with other observations like late volcanic activity and crustal density values as obtained from gravity and topography modelling [3, 4]. The stable density gradient, which establishes after the mantle overturn has completed, inhibits thermal convection. Albeit capable to provide stable reservoirs, this scenario suggests a conductive mantle after the overturn which on the one hand fails to sample deep regions of the mantle and on the other hand is clearly at odds with the volcanic history of Mars. This is best explained by assuming a convective mantle and partial melting as the principal agents responsible for the generation and evolution of Martian volcanism. Therefore, in this work an alternative scenario for the formation of early stable geochemical reservoirs is presented similar to the model of [5]. We investigate the influence of partial melting on mantle dynamics, crustal formation, and volcanic outgassing of a one-plate planet using a 2D mantle convection code. When melt is extracted to form crust, the mantle material left behind is more buoyant than its parent material and depleted in radioactive heat sources. The extracted heat-producing elements are then enriched in the crust, which also has an insulating effect due to its lower thermal conductivity compared to the mantle. In addition, partial melting can influence the mantle rheology through the dehydration (water depletion) of the mantle material by volcanic outgassing. As a consequence, the viscosity of water-depleted regions increases more than two orders of magnitude compared to water-saturated rocks resulting in slower cooling rates. The most important parameter influencing the thermo-chemical evolution is the assumed density difference between the primitive and the depleted mantle material (i.e., between peridotite and harzburgite). With small or negligible values of compositional buoyancy, crustal formation including crustal delamination is very efficient, also resulting in efficient processing and degassing of the mantle. The entire convecting mantle below the stagnant lid depletes continuously with time. In contrast, with increasing compositional buoyancy, crustal formation and mantle degassing are strongly suppressed although partial melting is substantially prolonged in the thermal evolution. The crust shows strong lateral variations in thickness, and crustal delamination is reduced and occurs only locally. Furthermore, two to four different mantle reservoirs can form depending on the initial temperature distribution [6]. Some of these reservoirs can be sustained during the entire evolution whereas others change with time - a scenario possibly valid for Mars as it may explain the isotope characteristic of the Martian meteorites. References: [1] Elkins-Tanton et al., 2005, EPSL; [2] Debaille et al., 2009, Nature; [3] Tosi et al., 2013, JGR; [4] Plesa et al., submitted to EPSL; [5] Ogawa and Yanagisawa 2011, JGR; [6] Plesa and Breuer, 2013, PSS.

MetPetDB: A database for metamorphic geochemistry

NASA Astrophysics Data System (ADS)

Spear, Frank S.; Hallett, Benjamin; Pyle, Joseph M.; Adalı, Sibel; Szymanski, Boleslaw K.; Waters, Anthony; Linder, Zak; Pearce, Shawn O.; Fyffe, Matthew; Goldfarb, Dennis; Glickenhouse, Nickolas; Buletti, Heather

2009-12-01

We present a data model for the initial implementation of MetPetDB, a geochemical database specific to metamorphic rock samples. The database is designed around the concept of preservation of spatial relationships, at all scales, of chemical analyses and their textural setting. Objects in the database (samples) represent physical rock samples; each sample may contain one or more subsamples with associated geochemical and image data. Samples, subsamples, geochemical data, and images are described with attributes (some required, some optional); these attributes also serve as search delimiters. All data in the database are classified as published (i.e., archived or published data), public or private. Public and published data may be freely searched and downloaded. All private data is owned; permission to view, edit, download and otherwise manipulate private data may be granted only by the data owner; all such editing operations are recorded by the database to create a data version log. The sharing of data permissions among a group of collaborators researching a common sample is done by the sample owner through the project manager. User interaction with MetPetDB is hosted by a web-based platform based upon the Java servlet application programming interface, with the PostgreSQL relational database. The database web portal includes modules that allow the user to interact with the database: registered users may save and download public and published data, upload private data, create projects, and assign permission levels to project collaborators. An Image Viewer module provides for spatial integration of image and geochemical data. A toolkit consisting of plotting and geochemical calculation software for data analysis and a mobile application for viewing the public and published data is being developed. Future issues to address include population of the database, integration with other geochemical databases, development of the analysis toolkit, creation of data models for derivative data, and building a community-wide user base. It is believed that this and other geochemical databases will enable more productive collaborations, generate more efficient research efforts, and foster new developments in basic research in the field of solid earth geochemistry.

Community-Based Development of Standards for Geochemical and Geochronological Data

NASA Astrophysics Data System (ADS)

Lehnert, K. A.; Walker, D.; Vinay, S.; Djapic, B.; Ash, J.; Falk, B.

2007-12-01

The Geoinformatics for Geochemistry (GfG) Program (www.geoinfogeochem.org) and the EarthChem project (www.earthchem.org) aim to maximize the application of geochemical data in Geoscience research and education by building a new advanced data infrastructure for geochemistry that facilitates the compilation, communication, serving, and visualization of geochemical data and their integration with the broad Geoscience data set. Building this new data infrastructure poses substantial challenges that are primarily cultural in nature, and require broad community involvement in the development and implementation of standards for data reporting (e.g., metadata for analytical procedures, data quality, and analyzed samples), data publication, and data citation to achieve broad acceptance and use. Working closely with the science community, with professional societies, and with editors and publishers, recommendations for standards for the reporting of geochemical and geochronological data in publications and to data repositories have been established, which are now under consideration for adoption in journal and agency policies. The recommended standards are aligned with the GfG and EarthChem data models as well as the EarthChem XML schema for geochemical data. Through partnerships with other national and international data management efforts in geochemistry and in the broader marine and terrestrial geosciences, GfG and EarthChem seek to integrate their development of geochemical metadata standards, data format, and semantics with relevant existing and emerging standards and ensure compatibility and compliance.

Using multifractal modeling as a standard tool in geochemical exploration for predicting mineralized areas

NASA Astrophysics Data System (ADS)

Gonçalves, Mario A.

2015-04-01

It has been 20 years since the pioneering work of Cheng et al (1994) first proposed a quantitative relationship for the areas enclosing concentration values of an element above given thresholds and their distribution in the field, known as concentration-area (CA) method, which is based in multifractal theory. The method allows the definition of geochemical anomalies in wide set of geological backgrounds but it took nearly 15 years before it became a widely used methodology for mineral exploration. The method was also extended to 1D and 3D data sets. It is worth noting the variety of methods that spanned from the theory of fractals. Building on previous models, including multiplicative cascades and size-grade relationships, increasing evidence points to the powerful tools of fractal theory to describe and model ore deposit distribution and formation. However, while much of these approaches become complex and not easy to use, the CA method is remarkable for its utter simplicity and disarming results obtained when confronted with the geological reality in the field. This is most valued by companies and professionals undertaking geochemical exploration surveys for the characterization or refining of potential ore targets or known mineralized areas. Several approaches have combined the CA method with geostatistic modeling and simulation and other established statistical techniques in order to enhance anomalous threshold identification. Examples are not restricted to geochemical exploration alone, other applications being studies on environmental change. Some of these examples will be addressed as they have been applied to different regions in the world, but particular emphasis will be put on geochemical exploration surveys in different geotectonic units of the Variscan basement in the Iberian Peninsula. These include quartz-vein gold mineralization in Northern Portugal and several surveys for base metals over two wide areas, which served to re-evaluate much of the scattered geochemical data sets that have been accumulating for decades of mining exploration in Southern Portugal. The studied zones include: the tectonic controlled quartz-vein Au-Sb mineralizations, the gabbroic and ultramafic complex of the southern border of the Ossa-Morena Zone, and the rocks belonging to the World-class massive sulfide province, the Iberian Pyrite Belt (IPB). The methodology used the CA method but also variogram analysis and modelling to outline and classify different sets of mineral deposits before confirmation in the field. This diversity of geologic contexts serves to show how effective and powerful the CA method can be, since it not only enhances already known mineralizations, it allowed the screening and identification of several new mineralized spots that have been previously overlooked. This has been of particularly economic importance because a major re-analysis of data and new exploration campaigns are currently under way for the next years in the IPB, with the potential for opening a new paradigm in the exploration for massive sulfide deposits in the region. Cheng et al, 1994, J. Geochem. Explor., 51, 109.

Chemical modeling constraints on Martian surface mineralogies formed in an early, warm, wet climate, and speculations on the occurrence of phosphate minerals in the Martian regolith

NASA Technical Reports Server (NTRS)

Plumlee, Geoffrey S.; Ridley, W. Ian; Debraal, Jeffrey D.

1992-01-01

This is one in a series of reports summarizing our chemical modeling studies of water-rock-gas interactions at the martian surface through time. The purpose of these studies is to place constraints on possible mineralogies formed at the martian surface and to model the geochemical implications of martian surficial processes proposed by previous researchers. Plumlee and Ridley summarize geochemical processes that may have occurred as a result of inferred volcano- and impact-driven hydrothermal activity on Mars. DeBraal et al. model the geochemical aspects of water-rock interactions and water evaporation near 0 C, as a prelude to future calculations that will model sub-0 C brine-rock-clathrate interactions under the current martian climate. In this report, we discuss reaction path calculations that model chemical processes that may have occurred at the martian surface in a postulated early, warm, wet climate. We assume a temperature of 25 C in all our calculations. Processes we model here include (1) the reaction of rainwater under various ambient CO2 and O2 pressures with basaltic rocks at the martian surface, (2) the formation of acid rain by volcanic gases such as HCl and SO2, (3) the reactions of acid rain with basaltic surficial materials, and (4) evaporation of waters resulting from rainwater-basalt interactions.

Hydrologic and geochemical data assimilation at the Hanford 300 Area

NASA Astrophysics Data System (ADS)

Chen, X.; Hammond, G. E.; Murray, C. J.; Zachara, J. M.

2012-12-01

In modeling the uranium migration within the Integrated Field Research Challenge (IFRC) site at the Hanford 300 Area, uncertainties arise from both hydrologic and geochemical sources. The hydrologic uncertainty includes the transient flow boundary conditions induced by dynamic variations in Columbia River stage and the underlying heterogeneous hydraulic conductivity field, while the geochemical uncertainty is a result of limited knowledge of the geochemical reaction processes and parameters, as well as heterogeneity in uranium source terms. In this work, multiple types of data, including the results from constant-injection tests, borehole flowmeter profiling, and conservative tracer tests, are sequentially assimilated across scales within a Bayesian framework to reduce the hydrologic uncertainty. The hydrologic data assimilation is then followed by geochemical data assimilation, where the goal is to infer the heterogeneous distribution of uranium sources using uranium breakthrough curves from a desorption test that took place at high spring water table. We demonstrate in our study that Ensemble-based data assimilation techniques (Ensemble Kalman filter and smoother) are efficient in integrating multiple types of data sequentially for uncertainty reduction. The computational demand is managed by using the multi-realization capability within the parallel PFLOTRAN simulator.

Geochemical evolution of groundwater in the Mud Lake area, eastern Idaho, USA

USGS Publications Warehouse

Rattray, Gordon W.

2015-01-01

Groundwater with elevated dissolved-solids concentrations—containing large concentrations of chloride, sodium, sulfate, and calcium—is present in the Mud Lake area of Eastern Idaho. The source of these solutes is unknown; however, an understanding of the geochemical sources and processes controlling their presence in groundwater in the Mud Lake area is needed to better understand the geochemical sources and processes controlling the water quality of groundwater at the Idaho National Laboratory. The geochemical sources and processes controlling the water quality of groundwater in the Mud Lake area were determined by investigating the geology, hydrology, land use, and groundwater geochemistry in the Mud Lake area, proposing sources for solutes, and testing the proposed sources through geochemical modeling with PHREEQC. Modeling indicated that sources of water to the eastern Snake River Plain aquifer were groundwater from the Beaverhead Mountains and the Camas Creek drainage basin; surface water from Medicine Lodge and Camas Creeks, Mud Lake, and irrigation water; and upward flow of geothermal water from beneath the aquifer. Mixing of groundwater with surface water or other groundwater occurred throughout the aquifer. Carbonate reactions, silicate weathering, and dissolution of evaporite minerals and fertilizer explain most of the changes in chemistry in the aquifer. Redox reactions, cation exchange, and evaporation were locally important. The source of large concentrations of chloride, sodium, sulfate, and calcium was evaporite deposits in the unsaturated zone associated with Pleistocene Lake Terreton. Large amounts of chloride, sodium, sulfate, and calcium are added to groundwater from irrigation water infiltrating through lake bed sediments containing evaporite deposits and the resultant dissolution of gypsum, halite, sylvite, and bischofite.

Benchmarking Geant4 for simulating galactic cosmic ray interactions within planetary bodies

DOE PAGES

Mesick, K. E.; Feldman, W. C.; Coupland, D. D. S.; ...

2018-06-20

Galactic cosmic rays undergo complex nuclear interactions with nuclei within planetary bodies that have little to no atmosphere. Radiation transport simulations are a key tool used in understanding the neutron and gamma-ray albedo coming from these interactions and tracing these signals back to geochemical composition of the target. In this paper, we study the validity of the code Geant4 for simulating such interactions by comparing simulation results to data from the Apollo 17 Lunar Neutron Probe Experiment. Different assumptions regarding the physics are explored to demonstrate how these impact the Geant4 simulation results. In general, all of the Geant4 resultsmore » over-predict the data, however, certain physics lists perform better than others. Finally, in addition, we show that results from the radiation transport code MCNP6 are similar to those obtained using Geant4.« less

Benchmarking Geant4 for simulating galactic cosmic ray interactions within planetary bodies

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mesick, K. E.; Feldman, W. C.; Coupland, D. D. S.

Galactic cosmic rays undergo complex nuclear interactions with nuclei within planetary bodies that have little to no atmosphere. Radiation transport simulations are a key tool used in understanding the neutron and gamma-ray albedo coming from these interactions and tracing these signals back to geochemical composition of the target. In this paper, we study the validity of the code Geant4 for simulating such interactions by comparing simulation results to data from the Apollo 17 Lunar Neutron Probe Experiment. Different assumptions regarding the physics are explored to demonstrate how these impact the Geant4 simulation results. In general, all of the Geant4 resultsmore » over-predict the data, however, certain physics lists perform better than others. Finally, in addition, we show that results from the radiation transport code MCNP6 are similar to those obtained using Geant4.« less

Data Qualification Report For: Thermodynamic Data File, DATA0.YMP.R0 For Geochemical Code, EQ3/6 

DOE Office of Scientific and Technical Information (OSTI.GOV)

P.L. Cloke

The objective of this work is to evaluate the adequacy of chemical thermodynamic data provided by Lawrence Livermore National Laboratory (LLNL) as DataO.ymp.ROA in response to an input request submitted under AP-3.14Q. This request specified that chemical thermodynamic data available in the file, Data0.com.R2, be updated, improved, and augmented for use in geochemical modeling used in Process Model Reports (PMRs) for Engineered Barrier Systems, Waste Form, Waste Package, Unsaturated Zone, and Near Field Environment, as well as for Performance Assessment. The data are qualified in the temperature range 0 to 100 C. Several Data Tracking Numbers (DTNs) associated with Analysis/Modelmore » Reports (AMR) addressing various aspects of the post-closure chemical behavior of the waste package and the Engineered Barrier System that rely on EQ316 outputs to which these data are used as input, are Principal Factor affecting. This qualification activity was accomplished in accordance with the AP-SIII.2Q using the Technical Assessment method. A development plan, TDP-EBS-MD-000044, was prepared in accordance with AP-2.13Q and approved by the Responsible Manager. In addition, a Process Control Evaluation was performed in accordance with AP-SV.1Q. The qualification method, selected in accordance with AP-SIII.2Q, was Technical Assessment. The rationale for this approach is that the data in File Data0.com.R2 are considered Handbook data and therefore do not themselves require qualification. Only changes to Data0.com.R2 required qualification. A new file has been produced which contains the database Data0.ymp.R0, which is recommended for qualification as a result of this action. Data0.ymp.R0 will supersede Data0.com.R2 for all Yucca Mountain Project (YMP) activities.« less

Geochemical Fate and Transport of Diphenhydramine and Cetirizine in Soil

NASA Astrophysics Data System (ADS)

Wireman, R.; Rutherford, C. J.; Vulava, V. M.; Cory, W. C.

2015-12-01

Pharmaceuticals compounds presence in natural soils and water around the world has become a growing concern. These compounds are being discharged into the environment through treated wastewater or municipal sludge applications. The main goal of this study is determine their geochemical fate in natural soils. In this study we investigated sorption and transport behavior of diphenhydramine (DPH) and cetirizine (CTZ) in natural soils. These two commonly-used antihistamines are complex aromatic hydrocarbons with polar functional groups. Two clean acidic soils (pH~4.5) were used for these studies - an A-horizon soil that had higher organic matter content (OM, 7.6%) and a B-horizon soil that had lower OM (1.6%), but higher clay content (5.1%). Sorption isotherms were measured using batch reactor experiments. Data indicated that sorption was nonlinear and that it was stronger in clay-rich soils. The pKa's of DPH and CTZ are 8.98 and 8.27 respectively, i.e., these compounds are predominantly in cationic form at soil pH. In these forms, they preferentially sorb to negatively charged mineral surfaces (e.g., clay) present in the soils. Soil clay mineral characterization indicated that kaolinite was the dominant clay mineral present along with small amount of montmorillonite. The nonlinear sorption isotherms were fitted with Freundlich model. Transport behavior of both compounds was measured using glass chromatography columns. As expected both DPH and CTZ were strongly retained in the clay-rich soil as compared with OM-rich soil. The asymmetrical shape of the breakthrough curves indicated that there were likely two separate sorption sites in the soil, each with different reaction rates with each compound. A two-region advection-dispersion transport code was used to model the transport breakthrough curves. There was no evidence of transformation or degradation of the compounds during our sorption and transport studies.

An archean suture zone in the Tobacco Root Mountains? (1984) Evolution of Archean Continental Crust, SW Montana (1985)

NASA Technical Reports Server (NTRS)

Mogk, D. W.; Kain, L.

1985-01-01

The Lake Plateau area of the Beartooth Mountains, Montana were mapped and geochemically sampled. The allochthonous nature of the Stillwater Complex was interpreted as a Cordilleran-style continental margin. The metamorphic and tectonic history of the Beartooth Mountains was addressed. The Archean geology of the Spanish Peaks area, northern Madison Range was addressed. A voluminous granulite terrain of supracrustal origin was identified, as well as a heretofore unknown Archean batholithic complex. Mapping, petrologic, and geochemical investigations of the Blacktail Mountains, on the western margin of the Wyoming Province, are completed. Mapping at a scale of 1:24000 in the Archean rocks of the Gravelly Range is near completion. This sequence is dominantly of stable-platform origin. Samples were collected for geothermometric/barometric analysis and for U-Pb zircon age dating. The analyses provide the basis for additional geochemical and geochronologic studies. A model for the tectonic and geochemical evolution of the Archean basement of SW Montana is presented.

A coupled THC model of the FEBEX in situ test with bentonite swelling and chemical and thermal osmosis

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zheng, L.; Samper, J.; Montenegro, L.

The performance assessment of a geological repository for radioactive waste requires quantifying the geochemical evolution of the bentonite engineered barrier. This barrier will be exposed to coupled thermal (T), hydrodynamic (H), mechanical (M) and chemical (C) processes. This paper presents a coupled THC model of the FEBEX (Full-scale Engineered Barrier EXperiment) in situ test which accounts for bentonite swelling and chemical and thermal osmosis. Model results attest the relevance of thermal osmosis and bentonite swelling for the geochemical evolution of the bentonite barrier while chemical osmosis is found to be almost irrelevant. The model has been tested with data collectedmore » after the dismantling of heater 1 of the in situ test. The model reproduces reasonably well the measured temperature, relative humidity, water content and inferred geochemical data. However, it fails to mimic the solute concentrations at the heater-bentonite and bentonite-granite interfaces because the model does not account for the volume change of bentonite, the CO{sub 2}(g) degassing and the transport of vapor from the bentonite into the granite. The inferred HCO{sub 3}{sup -} and pH data cannot be explained solely by solute transport, calcite dissolution and protonation/deprotonation by surface complexation, suggesting that such data may be affected also by other reactions.« less

Geochemical modeling of reactions and partitioning of trace metals and radionuclides during titration of contaminated acidic sediments.

PubMed

Zhang, Fan; Luo, Wensui; Parker, Jack C; Spalding, Brian P; Brooks, Scott C; Watson, David B; Jardine, Philip M; Gu, Baohua

2008-11-01

Many geochemical reactions that control aqueous metal concentrations are directly affected by solution pH. However, changes in solution pH are strongly buffered by various aqueous phase and solid phase precipitation/dissolution and adsorption/desorption reactions. The ability to predict acid-base behavior of the soil-solution system is thus critical to predict metal transport under variable pH conditions. This studywas undertaken to develop a practical generic geochemical modeling approach to predict aqueous and solid phase concentrations of metals and anions during conditions of acid or base additions. The method of Spalding and Spalding was utilized to model soil buffer capacity and pH-dependent cation exchange capacity by treating aquifer solids as a polyprotic acid. To simulate the dynamic and pH-dependent anion exchange capacity, the aquifer solids were simultaneously treated as a polyprotic base controlled by mineral precipitation/ dissolution reactions. An equilibrium reaction model that describes aqueous complexation, precipitation, sorption and soil buffering with pH-dependent ion exchange was developed using HydroGeoChem v5.0 (HGC5). Comparison of model results with experimental titration data of pH, Al, Ca, Mg, Sr, Mn, Ni, Co, and SO4(2-) for contaminated sediments indicated close agreement suggesting that the model could potentially be used to predictthe acid-base behavior of the sediment-solution system under variable pH conditions.

Model simulations of a field experiment on cation exchange-affected multicomponent solute transport in a sandy aquifer

NASA Astrophysics Data System (ADS)

Bjerg, Poul L.; Ammentorp, Hans C.; Christensen, Thomas H.

1993-04-01

A large-scale and long-term field experiment on cation exchange in a sandy aquifer has been modelled by a three-dimensional geochemical transport model. The geochemical model includes cation-exchange processes using a Gaines-Thomas expression, the closed carbonate system and the effects of ionic strength. Information on geology, hydrogeology and the transient conservative solute transport behaviour was obtained from a dispersion study in the same aquifer. The geochemical input parameters were carefully examined. CEC and selectivity coefficients were determined on the actual aquifer material by batch experiments and by the composition of the cations on the exchange complex. Potassium showed a non-ideal exchange behaviour with KCa selectivity coefficients indicating dependency on equivalent fraction and K + concentration in the aqueous phase. The model simulations over a distance of 35 m and a period of 250 days described accurately the observed attenuation of Na and the expelled amounts of Ca and Mg. Also, model predictions of plateau zones, formed by interaction with the background groundwater, in general agreed satisfactorily with the observations. Transport of K was simulated over a period of 800 days due to a substantially attenuation in the aquifer. The observed and the predicted breakthrough curves showed a reasonable accordance taking the duration of the experiment into account. However, some discrepancies were observed probably caused by the revealed non-ideal exchange behaviour of K +.

Geochemical modeling of iron, sulfur, oxygen and carbon in a coastal plain aquifer

NASA Astrophysics Data System (ADS)

Brown, C. J.; Schoonen, M. A. A.; Candela, J. L.

2000-11-01

Fe(III) reduction in the Magothy aquifer of Long Island, NY, results in high dissolved-iron concentrations that degrade water quality. Geochemical modeling was used to constrain iron-related geochemical processes and redox zonation along a flow path. The observed increase in dissolved inorganic carbon is consistent with the oxidation of sedimentary organic matter coupled to the reduction of O 2 and SO 42- in the aerobic zone, and to the reduction of SO 42- in the anaerobic zone; estimated rates of CO 2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO 2 production calculated from dissolved inorganic carbon mass transfer (2.55×10 -4 to 48.6×10 -4 mmol l -1 yr-1) generally were comparable to the calculated rates of CO 2 production by the combined reduction of O 2, Fe(III) and SO 42- (1.31×10 -4 to 15×10 -4 mmol l -1 yr-1). The overall increase in SO 42- concentrations along the flow path, together with the results of mass-balance calculations, and variations in δ34S values along the flow path indicate that SO 42- loss through microbial reduction is exceeded by SO 42- gain through diffusion from sediments and through the oxidation of FeS 2. Geochemical and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO 42--rich zones have been depleted by microbial reduction and resulted in localized SO 42--reducing zones in which the formation of iron disulfides decreases dissolved iron concentrations. These localized zones of SO 42- reduction, which are important for assessing zones of low dissolved iron for water-supply development, could be overlooked by aquifer studies that rely only on groundwater data from well-water samples for geochemical modeling.

An array processing system for lunar geochemical and geophysical data

NASA Technical Reports Server (NTRS)

Eliason, E. M.; Soderblom, L. A.

1977-01-01

A computerized array processing system has been developed to reduce, analyze, display, and correlate a large number of orbital and earth-based geochemical, geophysical, and geological measurements of the moon on a global scale. The system supports the activities of a consortium of about 30 lunar scientists involved in data synthesis studies. The system was modeled after standard digital image-processing techniques but differs in that processing is performed with floating point precision rather than integer precision. Because of flexibility in floating-point image processing, a series of techniques that are impossible or cumbersome in conventional integer processing were developed to perform optimum interpolation and smoothing of data. Recently color maps of about 25 lunar geophysical and geochemical variables have been generated.

The geochemical evolution of riparian ground water in a forested piedmont catchment

USGS Publications Warehouse

Burns, Douglas A.; Plummer, Niel; McDonnell, Jeffrey J.; Busenberg, Eurybiades; Casile, Gerolamo C.; Kendall, Carol; Hooper, Richard P.; Freer, James E.; Peters, Norman E.; Beven, Keith; Schlosser, Peter

2003-01-01

The principal weathering reactions and their rates in riparian ground water were determined at the Panola Mountain Research Watershed (PMRW) near Atlanta, Georgia. Concentrations of major solutes were measured in ground water samples from 19 shallow wells completed in the riparian (saprolite) aquifer and in one borehole completed in granite, and the apparent age of each sample was calculated from chloroflourocarbons and tritium/helium-3 data. Concentrations of SiO2, Na+, and Ca2+ generally increased downvalley and were highest in the borehole near the watershed outlet. Strong positive correlations were found between the concentrations of these solutes and the apparent age of ground water that was modern (zero to one year) in the headwaters, six to seven years midway down the valley, and 26 to 27 years in the borehole, located ∼500 m downstream from the headwaters. Mass-balance modeling of chemical evolution showed that the downstream changes in ground water chemistry could be largely explained by weathering of plagioclase to kaolinite, with possible contributions from weathering of K-feldspar, biotite, hornblende, and calcite. The in situ rates of weathering reactions were estimated by combining the ground water age dates with geochemical mass-balance modeling results. The weathering rate was highest for plagioclase (∼6.4 μmol/L/year), but could not be easily compared with most other published results for feldspar weathering at PMRW and elsewhere because the mineral-surface area to which ground water was exposed during geochemical evolution could not be estimated. However, a preliminary estimate of the mineral-surface area that would have contacted the ground water to provide the observed solute concentrations suggests that the plagioclase weathering rate calculated in this study is similar to the rate calculated in a previous study at PMRW, and three to four orders of magnitude slower than those published in previous laboratory studies of feldspar weathering. An accurate model of the geochemical evolution of riparian ground water is necessary to accurately model the geochemical evolution of stream water at PMRW.

A geochemical and geophysical reappraisal to the significance of the recent unrest at Campi Flegrei caldera (Southern Italy)

NASA Astrophysics Data System (ADS)

Moretti, Roberto; De Natale, Giuseppe; Troise, Claudia

2017-03-01

Volcanic unrest at calderas involves complex interaction between magma, hydrothermal fluids, and crustal stress and strain. Campi Flegrei caldera (CFc), located in the Naples (Italy) area and characterized by the highest volcanic risk on Earth for the extreme urbanization, undergoes unrest phenomena involving several meters of uplift and intense shallow microseismicity since several decades. Despite unrest episodes display in the last decade only moderate ground deformation and seismicity, current interpretations of geochemical data point to a highly pressurized hydrothermal system. We show that at CFc, the usual assumption of vapor-liquid coexistence in the fumarole plumes leads to largely overestimated hydrothermal pressures and, accordingly, interpretations of elevated unrest. By relaxing unconstrained geochemical assumptions, we infer an alternative model yielding better agreement between geophysical and geochemical observations. The model reconciles discrepancies between what observed (1) for two decades since the 1982-1984 large unrest, when shallow magma was supplying heat and fluids to the hydrothermal system, and (2) in the last decade. Compared to the 1980's unrest, the post-2005 phenomena are characterized by much lower aquifers overpressure and magmatic involvement, as indicated by geophysical data and despite large changes in geochemical indicators. Our interpretation points out a model in which shallow sills, intruded during 1969-1984, have completely cooled, so that fumarole emissions are affected now by deeper, CO2-richer, magmatic gases producing the modest heating and overpressure of the hydrothermal system. Our results have important implications on the short-term eruption hazard assessment and on the best strategies for monitoring and interpreting geochemical data.