Structurally Controlled Geothermal Systems in the Central Cascades Arc-Backarc Regime, Oregon

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

Wannamaker, Philip E.

The goal of this project has been to analyze available magnetotelluric (MT) geophysical surveys, structural geology based on mapping and LiDAR, and fluid geochemical data, to identify high-temperature fluid upwellings, critically stressed rock volumes, and other evidence of structurally-controlled geothermal resources. Data were to be integrated to create conceptual models of volcanic-hosted geothermal resources along the Central Cascades arc segment, especially in the vicinity of Mt. Jefferson to Three Sisters. LiDAR data sets available at Oregon State University (OSU) allowed detailed structural geology modeling through forest canopy. Copious spring and well fluid chemistries, including isotopes, were modeled using Geo-T andmore » TOUGHREACT software.« less

Hydrogeological modeling constraints provided by geophysical and geochemical mapping of a chlorinated ethenes plume in northern France

NASA Astrophysics Data System (ADS)

Razafindratsima, Stephen; Guérin, Roger; Bendjoudi, Hocine; de Marsily, Ghislain

2014-09-01

A methodological approach is described which combines geophysical and geochemical data to delineate the extent of a chlorinated ethenes plume in northern France; the methodology was used to calibrate a hydrogeological model of the contaminants' migration and degradation. The existence of strong reducing conditions in some parts of the aquifer is first determined by measuring in situ the redox potential and dissolved oxygen, dissolved ferrous iron and chloride concentrations. Electrical resistivity imaging and electromagnetic mapping, using the Slingram method, are then used to determine the shape of the pollutant plume. A decreasing empirical exponential relation between measured chloride concentrations in the water and aquifer electrical resistivity is observed; the resistivity formation factor calculated at a few points also shows a major contribution of chloride concentration in the resistivity of the saturated porous medium. MODFLOW software and MT3D99 first-order parent-daughter chain reaction and the RT3D aerobic-anaerobic model for tetrachloroethene (PCE)/trichloroethene (TCE) dechlorination are finally used for a first attempt at modeling the degradation of the chlorinated ethenes. After calibration, the distribution of the chlorinated ethenes and their degradation products simulated with the model approximately reflects the mean measured values in the observation wells, confirming the data-derived image of the plume.

Predictive Modeling of Terrestrial Radiation Exposure from Geologic Materials

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

Malchow, Russell L.; Haber, Daniel University of Nevada, Las Vegas; Burnley, Pamela

2015-01-01

Aerial gamma ray surveys are important for those working in nuclear security and industry for determining locations of both anthropogenic radiological sources and natural occurrences of radionuclides. During an aerial gamma ray survey, a low flying aircraft, such as a helicopter, flies in a linear pattern across the survey area while measuring the gamma emissions with a sodium iodide (NaI) detector. Currently, if a gamma ray survey is being flown in an area, the only way to correct for geologic sources of gamma rays is to have flown the area previously. This is prohibitively expensive and would require complete nationalmore » coverage. This project’s goal is to model the geologic contribution to radiological backgrounds using published geochemical data, GIS software, remote sensing, calculations, and modeling software. K, U and Th are the three major gamma emitters in geologic material. U and Th are assumed to be in secular equilibrium with their daughter isotopes. If K, U, and Th abundance values are known for a given geologic unit the expected gamma ray exposure rate can be calculated using the Grasty equation or by modeling software. Monte Carlo N-Particle Transport software (MCNP), developed by Los Alamos National Laboratory, is modeling software designed to simulate particles and their interactions with matter. Using this software, models have been created that represent various lithologies. These simulations randomly generate gamma ray photons at energy levels expected from natural radiologic sources. The photons take a random path through the simulated geologic media and deposit their energy at the end of their track. A series of nested spheres have been created and filled with simulated atmosphere to record energy deposition. Energies deposited are binned in the same manner as the NaI detectors used during an aerial survey. These models are used in place of the simplistic Grasty equation as they take into account absorption properties of the lithology which the simplistic equation ignores.« less

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

Numerical Modelling of Speleothem and Dripwater Chemistry: Interpreting Coupled Trace Element and Isotope Proxies for Climate Reconstructions

NASA Astrophysics Data System (ADS)

Owen, R.; Day, C. C.; Henderson, G. M.

2016-12-01

Speleothem palaeoclimate records are widely used but are often difficult to interpret due to the geochemical complexity of the soil-karst-cave system. Commonly analysed proxies (e.g. δ18O, δ13C and Mg/Ca) may be affected by multiple processes along the water flow path from atmospheric moisture source through to the cave drip site. Controls on speleothem chemistry include rainfall and aerosol chemistry, bedrock chemistry, temperature, soil pCO2, the degree of open-system dissolution and prior calcite precipitation. Disentangling the effects of these controls is necessary to fully interpret speleothem palaeoclimate records. To quantify the effects of these processes, we have developed an isotope-enabled numerical model based on the geochemical modelling software PHREEQC. The model calculates dripwater chemistry and isotopes through equilibrium bedrock dissolution and subsequent iterative CO2 degassing and calcite precipitation. This approach allows forward modelling of dripwater and speleothem proxies, both chemical (e.g. Ca concentration, pH, Mg/Ca and Sr/Ca ratios) and isotopic (e.g. δ18O, δ13C, δ44Ca and radiocarbon content), in a unified framework. Potential applications of this model are varied and the model may be readily expanded to include new isotope systems or processes. Here we focus on calculated proxy co-variation due to changes in model parameters. Examples include: - The increase in Ca concentration, decrease in δ13C and increase in radiocarbon content as bedrock dissolution becomes more open-system. - Covariation between δ13C, δ44Ca and trace metal proxies (e.g. Mg/Ca) predicted by changing prior calcite precipitation. - The effect of temperature change on all proxies through the soil-karst-cave system. Separating the impact of soil and karst processes on geochemical proxies allows more quantitative reconstruction of the past environment, and greater understanding in modern cave monitoring studies.

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

Newberry Volcano EGS Demonstration Stimulation Modeling

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

Trenton T. Cladouhos, Matthew Clyne, Maisie Nichols,; Susan Petty, William L. Osborn, Laura Nofziger

2011-10-23

As a part of Phase I of the Newberry Volcano EGS Demonstration project, several data sets were collected to characterize the rock volume around the well. Fracture, fault, stress, and seismicity data has been collected by borehole televiewer, LiDAR elevation maps, and microseismic monitoring. Well logs and cuttings from the target well (NWG 55-29) and core from a nearby core hole (USGS N-2) have been analyzed to develop geothermal, geochemical, mineralogical and strength models of the rock matrix, altered zones, and fracture fillings (see Osborn et al., this volume). These characterization data sets provide inputs to models used to planmore » and predict EGS reservoir creation and productivity. One model used is AltaStim, a stochastic fracture and flow software model developed by AltaRock. The software's purpose is to model and visualize EGS stimulation scenarios and provide guidance for final planning. The process of creating an AltaStim model requires synthesis of geologic observations at the well, the modeled stress conditions, and the stimulation plan. Any geomechanical model of an EGS stimulation will require many assumptions and unknowns; thus, the model developed here should not be considered a definitive prediction, but a plausible outcome given reasonable assumptions. AltaStim is a tool for understanding the effect of known constraints, assumptions, and conceptual models on plausible outcomes.« less

National geochemical data base; PLUTO geochemical data base for the United States

USGS Publications Warehouse

Baedecker, Philip A.; Grossman, Jeffrey N.; Buttleman, Kim P.

1998-01-01

The PLUTO CD-ROM data base contains inorganic geothermal data obtained by the analytical laboratories of the Geologic Division of the U.S. Geological Survey (USGS) for the United States, including Hawaii and Alaska, in support of USGS program activities requiring chemical data. This CD-ROM was produced in accordance with the ISO 9660 standard and can be accessed by any computer system that has the appropriate software to read the ISO 9660 discs; however, the disc is intended for use in a DOS environment.

ThermoFit: A Set of Software Tools, Protocols and Schema for the Organization of Thermodynamic Data and for the Development, Maintenance, and Distribution of Internally Consistent Thermodynamic Data/Model Collections

NASA Astrophysics Data System (ADS)

Ghiorso, M. S.

2013-12-01

Internally consistent thermodynamic databases are critical resources that facilitate the calculation of heterogeneous phase equilibria and thereby support geochemical, petrological, and geodynamical modeling. These 'databases' are actually derived data/model systems that depend on a diverse suite of physical property measurements, calorimetric data, and experimental phase equilibrium brackets. In addition, such databases are calibrated with the adoption of various models for extrapolation of heat capacities and volumetric equations of state to elevated temperature and pressure conditions. Finally, these databases require specification of thermochemical models for the mixing properties of solid, liquid, and fluid solutions, which are often rooted in physical theory and, in turn, depend on additional experimental observations. The process of 'calibrating' a thermochemical database involves considerable effort and an extensive computational infrastructure. Because of these complexities, the community tends to rely on a small number of thermochemical databases, generated by a few researchers; these databases often have limited longevity and are universally difficult to maintain. ThermoFit is a software framework and user interface whose aim is to provide a modeling environment that facilitates creation, maintenance and distribution of thermodynamic data/model collections. Underlying ThermoFit are data archives of fundamental physical property, calorimetric, crystallographic, and phase equilibrium constraints that provide the essential experimental information from which thermodynamic databases are traditionally calibrated. ThermoFit standardizes schema for accessing these data archives and provides web services for data mining these collections. Beyond simple data management and interoperability, ThermoFit provides a collection of visualization and software modeling tools that streamline the model/database generation process. Most notably, ThermoFit facilitates the rapid visualization of predicted model outcomes and permits the user to modify these outcomes using tactile- or mouse-based GUI interaction, permitting real-time updates that reflect users choices, preferences, and priorities involving derived model results. This ability permits some resolution of the problem of correlated model parameters in the common situation where thermodynamic models must be calibrated from inadequate data resources. The ability also allows modeling constraints to be imposed using natural data and observations (i.e. petrologic or geochemical intuition). Once formulated, ThermoFit facilitates deployment of data/model collections by automated creation of web services. Users consume these services via web-, excel-, or desktop-clients. ThermoFit is currently under active development and not yet generally available; a limited capability prototype system has been coded for Macintosh computers and utilized to construct thermochemical models for H2O-CO2 mixed fluid saturation in silicate liquids. The longer term goal is to release ThermoFit as a web portal application client with server-based cloud computations supporting the modeling environment.

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.

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.

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.

Geologic Communications | Alaska Division of Geological & Geophysical

Science.gov Websites

improves a database for the Division's digital and map-based geological, geophysical, and geochemical data interfaces DGGS metadata and digital data distribution - Geospatial datasets published by DGGS are designed to be compatible with a broad variety of digital mapping software, to present DGGS's geospatial data

Modeling low-temperature geochemical processes: Chapter 2

USGS Publications Warehouse

Nordstrom, D. Kirk; Campbell, Kate M.

2014-01-01

This chapter provides an overview of geochemical modeling that applies to water–rock interactions under ambient conditions of temperature and pressure. Topics include modeling definitions, historical background, issues of activity coefficients, popular codes and databases, examples of modeling common types of water–rock interactions, and issues of model reliability. Examples include speciation, microbial redox kinetics and ferrous iron oxidation, calcite dissolution, pyrite oxidation, combined pyrite and calcite dissolution, dedolomitization, seawater–carbonate groundwater mixing, reactive-transport modeling in streams, modeling catchments, and evaporation of seawater. The chapter emphasizes limitations to geochemical modeling: that a proper understanding and ability to communicate model results well are as important as completing a set of useful modeling computations and that greater sophistication in model and code development is not necessarily an advancement. If the goal is to understand how a particular geochemical system behaves, it is better to collect more field data than rely on computer codes.

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.

The Possibility of Serpentinization on Enceladus

NASA Astrophysics Data System (ADS)

Boye Nissen, Silas; Taubner, Ruth-Sophie; Leitner, Johannes J.; Firneis, Maria G.

2015-04-01

Enceladus, an icy moon of Saturn, is considered a potential place for a second genesis of life in a rather easy accessible environment for future space missions. Plumes erupting material out to Saturn's E-ring have been detected. This and the latest gravity data received by NASA's Cassini probe highly suggest a regional subsurface water reservoir to exist. The Cassini mission has detected all the necessary elements for life in these plumes and few ecosystems known from the Earth could in theory exist on this moon. These ecosystems rely on a sufficient hydrogen source. On Earth, this molecule is known to be produced when the oceans interact with silicates in the mantle through a geochemical process called serpentinization. The aim of this work is first of all to determine whether serpentinization is possible on Enceladus at an assumed pressure range of 25-50 bars and a temperature range of 0-50 °C. If it is possible, then the second aim is to estimate the hydrogen production rate. Therefore, geochemical modeling with the EQ3/6 software package has been applied and thermody-namic databases have been created with the use of the DBCreate program. Due to the present lack of data on the chemical composition of the subsurface sea it is assumed to be equal to the chemical composition of the plumes. Finally, the importance of the pH value (7-9) is studied. Models based on Enceladus being a captured comet or aggregated ring material are used to specify the core composition. These are different combinations of pyroxenes and olivines or no silicates, respectively. The results show that serpentinization is ongoing in all models with silicates at these low pressure and temperature ranges. Hydrogen is only produced in the presence of an iron-rich end-member of either pyroxene (ferrosilite) or olivine (fayalite). In these cases, the dihydrogen concentration is estimated to be around 2.0 mg/kg solution. Changing the pH value did not make any differences to the results. It is discussed how the complex aqueous solution assumed (based on the composition of the plumes) affect the serpentinization process on Enceladus. In addition, it is discussed if the pressure and temperature ranges are too small to make a significant change on the dihydrogen concentration. This work shows that serpentinization seems to be possible at the low pressure and temperature values found at the water/rock boundary on Enceladus. To improve the estimate of the hydrogen production rate, measurements for the silicate ions in the aqueous solution on Enceladus are needed. We would like to thank Wolfgang Bach from the Research Group "Petrology of the Ocean Crust" (University of Bremen) for his very helpful support on applying the geochemical modeling software.

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.

Integration of bio- and geoscience data with the ODM2 standards and software ecosystem for the CZOData and BiG CZ Data projects

NASA Astrophysics Data System (ADS)

Aufdenkampe, A. K.; Mayorga, E.; Horsburgh, J. S.; Lehnert, K. A.; Zaslavsky, I.

2015-12-01

We have developed a family of solutions to the challenges of integrating diverse data from of biological and geological (BiG) disciplines for Critical Zone (CZ) science. These standards and software solutions have been developed around the new Observations Data Model version 2.0 (ODM2, http://ODM2.org), which was designed as a profile of the Open Geospatial Consortium's (OGC) Observations and Measurements (O&M) standard. The ODM2 standards and software ecosystem has at it's core an information model that balances specificity with flexibility to powerfully and equally serve the needs of multiple dataset types, from multivariate sensor-generated time series to geochemical measurements of specimen hierarchies to multi-dimensional spectral data to biodiversity observations. ODM2 has been adopted as the information model guiding the next generation of cyberinfrastructure development for the Interdisciplinary Earth Data Alliance (http://www.iedadata.org/) and the CUAHSI Water Data Center (https://www.cuahsi.org/wdc). Here we present several components of the ODM2 standards and software ecosystem that were developed specifically to help CZ scientists and their data managers to share and manage data through the national Critical Zone Observatory data integration project (CZOData, http://criticalzone.org/national/data/) and the bio integration with geo for critical zone science data project (BiG CZ Data, http://bigcz.org/). These include the ODM2 Controlled Vocabulary system (http://vocabulary.odm2.org), the YAML Observation Data Archive & exchange (YODA) File Format (https://github.com/ODM2/YODA-File) and the BiG CZ Toolbox, which will combine easy-to-install ODM2 databases (https://github.com/ODM2/ODM2) with a variety of graphical software packages for data management such as ODMTools (https://github.com/ODM2/ODMToolsPython) and the ODM2 Streaming Data Loader (https://github.com/ODM2/ODM2StreamingDataLoader).

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.

Reactive transport simulation via combination of a multiphase-capable transport code for unstructured meshes with a Gibbs energy minimization solver of geochemical equilibria

NASA Astrophysics Data System (ADS)

Fowler, S. J.; Driesner, T.; Hingerl, F. F.; Kulik, D. A.; Wagner, T.

2011-12-01

We apply a new, C++-based computational model for hydrothermal fluid-rock interaction and scale formation in geothermal reservoirs. The model couples the Complex System Modelling Platform (CSMP++) code for fluid flow in porous and fractured media (Matthai et al., 2007) with the Gibbs energy minimization numerical kernel GEMS3K of the GEM-Selektor (GEMS3) geochemical modelling package (Kulik et al., 2010) in a modular fashion. CSMP++ includes interfaces to commercial file formats, accommodating complex geometry construction using CAD (Rhinoceros) and meshing (ANSYS) software. The CSMP++ approach employs finite element-finite volume spatial discretization, implicit or explicit time discretization, and operator splitting. GEMS3K can calculate complex fluid-mineral equilibria based on a variety of equation of state and activity models. A selection of multi-electrolyte aqueous solution models, such as extended Debye-Huckel, Pitzer (Harvie et al., 1984), EUNIQUAC (Thomsen et al., 1996), and the new ELVIS model (Hingerl et al., this conference), makes it well-suited for application to a wide range of geothermal conditions. An advantage of the GEMS3K solver is simultaneous consideration of complex solid solutions (e.g., clay minerals), gases, fluids, and aqueous solutions. Each coupled simulation results in a thermodynamically-based description of the geochemical and physical state of a hydrothermal system evolving along a complex P-T-X path. The code design allows efficient, flexible incorporation of numerical and thermodynamic database improvements. We demonstrate the coupled code workflow and applicability to compositionally and physically complex natural systems relevant to enhanced geothermal systems, where temporally and spatially varying chemical interactions may take place within diverse lithologies of varying geometry. Engesgaard, P. & Kipp, K. L. (1992). Water Res. Res. 28: 2829-2843. Harvie, C. E.; Møller, N. & Weare, J. H. (1984). Geochim. Cosmochim. Acta 48: 723-751. Kulik, D. A., Wagner, T., Dmytrieva S. V, et al. (2010). GEM-Selektor home page, Paul Scherrer Institut. Available at http://gems.web.psi.ch. Matthäi, S. K., Geiger, S., Roberts, S. G., Paluszny, A., Belayneh, M., Burri, A., Mezentsev, A., Lu, H., Coumou, D., Driesner, T. & Heinrich C. A. (2007). Geol. Soc. London, Spec. Publ. 292: 405-429. Thomsen, K. Rasmussen, P. & Gani, R. (1996). Chem. Eng. Sci. 51: 3675-3683.

VS2DRTI: Simulating Heat and Reactive Solute Transport in Variably Saturated Porous Media.

PubMed

Healy, Richard W; Haile, Sosina S; Parkhurst, David L; Charlton, Scott R

2018-01-29

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid-rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult-to-use models. To address the need for a simple and easy-to-use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two-dimensional, constant-density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature-dependent cation exchange. VS2DRTI is freely available public domain software. © 2018, National Ground Water Association.

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.

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

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

Integrating Thermodynamic Models in Geodynamic Simulations: The Example of the Community Software ASPECT

NASA Astrophysics Data System (ADS)

Dannberg, J.; Heister, T.; Grove, R. R.; Gassmoeller, R.; Spiegelman, M. W.; Bangerth, W.

2017-12-01

Earth's surface shows many features whose genesis can only be understood through the interplay of geodynamic and thermodynamic models. This is particularly important in the context of melt generation and transport: Mantle convection determines the distribution of temperature and chemical composition, the melting process itself is then controlled by the thermodynamic relations and in turn influences the properties and the transport of melt. Here, we present our extension of the community geodynamics code ASPECT, which solves the equations of coupled magma/mantle dynamics, and allows to integrate different parametrizations of reactions and phase transitions: They may alternatively be implemented as simple analytical expressions, look-up tables, or computed by a thermodynamics software. As ASPECT uses a variety of numerical methods and solvers, this also gives us the opportunity to compare different approaches of modelling the melting process. In particular, we will elaborate on the spatial and temporal resolution that is required to accurately model phase transitions, and show the potential of adaptive mesh refinement when applied to melt generation and transport. We will assess the advantages and disadvantages of iterating between fluid dynamics and chemical reactions derived from thermodynamic models within each time step, or decoupling them, allowing for different time step sizes. Beyond that, we will expand on the functionality required for an interface between computational thermodynamics and fluid dynamics models from the geodynamics side. Finally, using a simple example of melting of a two-phase, two-component system, we compare different time-stepping and solver schemes in terms of accuracy and efficiency, in dependence of the time scales of fluid flow and chemical reactions relative to each other. Our software provides a framework to integrate thermodynamic models in high resolution, 3d simulations of coupled magma/mantle dynamics, and can be used as a tool to study links between physical processes and geochemical signals in the Earth.

EQ3NR, a computer program for geochemical aqueous speciation-solubility calculations: Theoretical manual, user`s guide, and related documentation (Version 7.0); Part 3

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

Wolery, T.J.

1992-09-14

EQ3NR is an aqueous solution speciation-solubility modeling code. It is part of the EQ3/6 software package for geochemical modeling. It computes the thermodynamic state of an aqueous solution by determining the distribution of chemical species, including simple ions, ion pairs, and complexes, using standard state thermodynamic data and various equations which describe the thermodynamic activity coefficients of these species. The input to the code describes the aqueous solution in terms of analytical data, including total (analytical) concentrations of dissolved components and such other parameters as the pH, pHCl, Eh, pe, and oxygen fugacity. The input may also include a desiredmore » electrical balancing adjustment and various constraints which impose equilibrium with special pure minerals, solid solution end-member components (of specified mole fractions), and gases (of specified fugacities). The code evaluates the degree of disequilibrium in terms of the saturation index (SI = 1og Q/K) and the thermodynamic affinity (A = {minus}2.303 RT log Q/K) for various reactions, such as mineral dissolution or oxidation-reduction in the aqueous solution itself. Individual values of Eh, pe, oxygen fugacity, and Ah (redox affinity) are computed for aqueous redox couples. Equilibrium fugacities are computed for gas species. The code is highly flexible in dealing with various parameters as either model inputs or outputs. The user can specify modification or substitution of equilibrium constants at run time by using options on the input file.« less

Module-oriented modeling of reactive transport with HYTEC

NASA Astrophysics Data System (ADS)

van der Lee, Jan; De Windt, Laurent; Lagneau, Vincent; Goblet, Patrick

2003-04-01

The paper introduces HYTEC, a coupled reactive transport code currently used for groundwater pollution studies, safety assessment of nuclear waste disposals, geochemical studies and interpretation of laboratory column experiments. Based on a known permeability field, HYTEC evaluates the groundwater flow paths, and simulates the migration of mobile matter (ions, organics, colloids) subject to geochemical reactions. The code forms part of a module-oriented structure which facilitates maintenance and improves coding flexibility. In particular, using the geochemical module CHESS as a common denominator for several reactive transport models significantly facilitates the development of new geochemical features which become automatically available to all models. A first example shows how the model can be used to assess migration of uranium from a sub-surface source under the effect of an oxidation front. The model also accounts for alteration of hydrodynamic parameters (local porosity, permeability) due to precipitation and dissolution of mineral phases, which potentially modifies the migration properties in general. The second example illustrates this feature.

Parametrization of environment by geochemistry of the varved clastic and bio- chemogenic lake sediments

NASA Astrophysics Data System (ADS)

Kalugin, Ivan; Darin, Andrey; Babich, Valery; Markovich, Tatiana; Meydan, Feray

2017-04-01

As it well known, recent quantitative estimations of high-resolution environmental variability are based on geochemical records in lake sediments. Naturally, annually laminated sediments (varves) are the best objects for paleoclimatic study, because they allow to investigate seasonal variability for understanding long-term environmental pattern. Also, varved sediments seem to be applied as the model for identification of element-indicators for non-laminated sediments. The XRF scanner on Synchrotron Radiation provides big geochemical dataset for next mathematic treatment, including time series construction. XRF scanning realizes rapid and non-destructive determinations more than 30 trace elements in a range of concentration from 1 up to 10000 ppm in annual layers. That makes sedimentary cores comparable with tree-rings. Geochemical and physicochemical investigation of lake sediments provides basic information to identify geochemical signals with paleoclimate. In general, sediment consists of mineral component, organics and carbonates. The proportions between these components are affected by environmental parameters, because measured element content or their combinations show correlation with meteodata on instrumental time interval. That allows applying geochemical variability to reconstruct the environmental parameters in the form of time series. The proportions between main components are controlled by temperature, atmospheric precipitation, water salinity and other external forcings. So, layered structure of lake bottom sediments and detectable elements content variability both represent a continuous record of environmental history. Element composition and it's climatic response. Bottom sediments represent conditions of physical weathering, temperate bioproductivity and aridity, which concern to mountain lakes within extra tropical zone. The numerical values of the parameters can be computed by software of physical-chemical modeling for gas+water+rock multisystems. Mineral matter responses to runoff. Mineral clastic part is correlated with x-ray density. It includes "clastic" rock-forming - Si, Al , Ti, Fe, Mg, Ca, K and trace elements such as Sr, Rb, Y, Zr, REE etc. Organic component of sediment more reflects temperature by means of productivity in the catchment and waterbody. Organophillic elements are Br, I, U and others soluble elements correlated with organic Carbon or LOI<500oC. Bio-chemogenic component is more characteristic for saline lakes, where Ca-, Mg- and Sr- carbonates precipitated in dependence of temperature, aridity and water salinity. Separate geochemical indicators are directly used for paleo- environmental evaluation. For example, elements with changing valency may be a proxy of outer conditions. Fe is strictly connected with sulfur in sulphide under anoxic conditions. And also Fe forms siderite in carbonate ion saturated, but calcium poor, water in the sedimentation system. Mn-enriched layers, crusts and nodules mark usually a long - term pauses of sedimentation in oxic systems. Mo/Mn ratio is good correlated with anoxic atmosphere. And so on. The work is supported by grants RFBR 16-05-00641, 16-05-00657, 15-55-46001.

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.

Microbial respiration and dissolution precipitation reactions of minerals: thermo-kinetics and reactive transport modelling

NASA Astrophysics Data System (ADS)

Azaroual, M. M.; Parmentier, M.; Andre, L.; Croiset, N.; Pettenati, M.; Kremer, S.

2010-12-01

Microbial processes interact closely with abiotic geochemical reactions and mineralogical transformations in several hydrogeochemical systems. Reactive transport models are aimed to analyze these complex mechanisms integrating as well as the degradation of organic matter as the redox reactions involving successive terminal electron acceptors (TEAPs) mediated by microbes through the continuum of unsaturated zone (soil) - saturated zone (aquifer). The involvement of microbial processes in reactive transport in soil and subsurface geologic greatly complicates the mastery of the major mechanisms and the numerical modelling of these systems. The introduction of kinetic constraints of redox reactions in aqueous phase requires the decoupling of equilibrium reactions and the redefinition of mass balance of chemical elements including the concept of basis species and secondary species of thermodynamic databases used in geochemical modelling tools. An integrated methodology for modelling the reactive transport has been developed and implemented to simulate the transfer of arsenic, denitrification processes and the role of metastable aqueous sulfur species with pyrite and organic matter as electron donors entities. A mechanistic rate law of microbial respiration in various geochemical environments was used to simulate reactive transport of arsenic, nitrate and organic matter combined to the generalized rate law of mineral dissolution - precipitation reactions derived from the transition state theory was used for dissolution - precipitation of silica, aluminosilicate, carbonate, oxyhydroxide, and sulphide minerals. The kinetic parameters are compiled from the literature measurements based on laboratory constrained experiments and field observations. Numerical simulations, using the geochemical software PHREEQC, were performed aiming to identify the key reactions mediated by microbes in the framework of in the first hand the concept of the unsaturated - saturated zones of an artificial recharge of deep aquifers system and in a second hand an acid mine drainage system. A large amount of data is available on the old mine site of Cheni (France). This field data on acid mine drainage are compared to a thermokinetic model including biological kinetics, precipitation-dissolution kinetics and surface complexation on ferrihydrite. The kinetic parameters are from literature and from a fitting on batch biological experiments. The integrated approach combining reaction kinetics and biogeochemical thermodynamic constraints is successfully applied to denitrification experiments in the presence of acetate and pyrite conducted in the laboratory for batch and column systems. The powerful of this coupled approach allows a fine description of the different transition species from nitrate to nitrogen. The fitted kinetic parameters established for modelling these laboratory results are thus extended to simulate the denitrification processes in a field case where organic matter and pyrite FeS2 are the electron donors and O2, NO3, Fe(OH)3, SO4 are the electron acceptors in the framework of a continuum UZ - SZ aiming to identify the stabilized redox zones of acid mine drainage. The detailed results obtained on two actual case studies will be presented.

Benchmarking a Visual-Basic based multi-component one-dimensional reactive transport modeling tool

NASA Astrophysics Data System (ADS)

Torlapati, Jagadish; Prabhakar Clement, T.

2013-01-01

We present the details of a comprehensive numerical modeling tool, RT1D, which can be used for simulating biochemical and geochemical reactive transport problems. The code can be run within the standard Microsoft EXCEL Visual Basic platform, and it does not require any additional software tools. The code can be easily adapted by others for simulating different types of laboratory-scale reactive transport experiments. We illustrate the capabilities of the tool by solving five benchmark problems with varying levels of reaction complexity. These literature-derived benchmarks are used to highlight the versatility of the code for solving a variety of practical reactive transport problems. The benchmarks are described in detail to provide a comprehensive database, which can be used by model developers to test other numerical codes. The VBA code presented in the study is a practical tool that can be used by laboratory researchers for analyzing both batch and column datasets within an EXCEL platform.

Autonomous cloud based site monitoring through hydro geophysical data assimilation, processing and result delivery

NASA Astrophysics Data System (ADS)

Versteeg, R.; Johnson, D. V.; Rodzianko, A.; Zhou, H.; Dafflon, B.; Leger, E.; de Kleine, M.

2017-12-01

Understanding of processes in the shallow subsurface requires that geophysical, biogeochemical, hydrological and remote sensing datasets are assimilated, processed and interpreted. Multiple enabling software capabilities for process understanding have been developed by the science community. These include information models (ODM2), reactive transport modeling (PFLOTRAN, Modflow, CLM, Landlab), geophysical inversion (E4D, BERT), parameter estimation (PEST, DAKOTA), visualization (ViSiT, Paraview, D3, QGIS) as well as numerous tools written in python and R for petrophysical mapping, stochastic modeling, data analysis and so on. These capabilities use data collected using sensors and analytical tools developed by multiple manufacturers which produce many different measurements. While scientists obviously leverage tools, capabilities and lessons learned from one site at other sites, the current approach to site characterization and monitoring is very labor intensive and does not scale well. Our objective is to be able to monitor many (hundreds - thousands) of sites. This requires that monitoring can be done in a near time, affordable, auditable and essentially autonomous manner. For this we have developed a modular vertically integrated cloud based software framework which was designed from the ground up for effective site and process monitoring. This software framework (PAF - Predictive Assimilation Framework) is multitenant software and provides automation of data ingestion, processing and visualization of hydrological, geochemical and geophysical (ERT/DTS) data. The core organizational element of PAF is a project/user one in which capabilities available to users are controlled by a combination of available data and access permissions. All PAF capabilities are exposed through APIs, making it easy to quickly add new components. PAF is fully integrated with newly developed autonomous electrical geophysical hardware and thus allows for automation of electrical geophysical ingestion and processing and the ability for co analysis and visualization of the raw and processed data with other data of interest (e.g. soil temperature, soil moisture, precipitation). We will demonstrate current PAF capabilities and discuss future efforts.

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

A Geochemical Reaction Model for Titration of Contaminated Soil and Groundwater at the Oak Ridge Reservation

NASA Astrophysics Data System (ADS)

Zhang, F.; Parker, J. C.; Gu, B.; Luo, W.; Brooks, S. C.; Spalding, B. P.; Jardine, P. M.; Watson, D. B.

2007-12-01

This study investigates geochemical reactions during titration of contaminated soil and groundwater at the Oak Ridge Reservation in eastern Tennessee. The soils and groundwater exhibits low pH and high concentrations of aluminum, calcium, magnesium, manganese, various trace metals such as nickel and cobalt, and radionuclides such as uranium and technetium. The mobility of many of the contaminant species diminishes with increasing pH. However, base additions to increase pH are strongly buffered by various precipitation/dissolution and adsorption/desorption reactions. The ability to predict acid-base behavior and associated geochemical effects is thus critical to evaluate remediation performance of pH manipulation strategies. This study was undertaken to develop a practical but generally applicable geochemical model to predict aqueous and solid-phase speciation during soil and groundwater titration. To model titration in the presence of aquifer solids, an approach proposed by Spalding and Spalding (2001) was utilized, which treats aquifer solids as a polyprotic acid. Previous studies have shown that Fe and Al-oxyhydroxides strongly sorb dissolved Ni, U and Tc species. In this study, since the total Fe concentration is much smaller than that of Al, only ion exchange reactions associated with Al hydroxides are considered. An equilibrium reaction model that includes aqueous complexation, precipitation, ion exchange, and soil buffering reactions was developed and implemented in the code HydroGeoChem 5.0 (HGC5). Comparison of model results with experimental titration curves for contaminated groundwater alone and for soil- water systems indicated close agreement. This study is expected to facilitate field-scale modeling of geochemical processes under conditions with highly variable pH to develop practical methods to control contaminant mobility at geochemically complex sites.

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.

Geochemical baseline distribution of harmful elements in the surface soils of Campania region.

NASA Astrophysics Data System (ADS)

Albanese, Stefano; Lima, Annamaria; Qu, Chengkai; Cicchella, Domenico; Buccianti, Antonella; De Vivo, Benedetto

2015-04-01

Environmental geochemical mapping has assumed an increasing relevance and the separation of values to discriminate between anthropogenic pollution and natural (geogenic) sources has become crucial to address environmental problems affecting the quality of life of human beings. In the last decade, a number of geochemical prospecting projects, mostly focused on surface soils (topsoils), were carried out at different scales (from regional to local) across the whole Campania region (Italy) to characterize the distribution of both harmful elements and persistent organic pollutants (POP) in the environment and to generating a valuable database to serve as reference in developing geomedical studies. During the 2014, a database reporting the distribution of 53 chemical elements in 3536 topsoil samples, collected across the whole region, was completed. The geochemical data, after necessary quality controls, were georeferenced and processed in a geochemistry dedicated GIS software named GEODAS. For each considered element a complete set of maps was generated to depict both the discrete and the spatially continuous (interpolated) distribution of elemental concentrations across the region. The interpolated maps were generated using the Multifractal Inverse Distance eighted (MIDW) algorithm. Subsequently, the S-A method, also implemented in GEODAS, was applied to MIDW maps to eliminate spatially limited anomalies from the original grid and to generate the distribution patterns of geochemical baselines for each element. For a selected group of elements geochemical data were also treated by means of a Compositional Data Analysis (CoDA) aiming at investigating the regionalised structure of the data by considering the joint behaviour of several elements constituting for each sample its whole composition. A regional environmental risk assessment was run on the basis of the regional distribution of heavy metals in soil, land use types and population. The risk assessment produced a ranking of priorities and located areas of regional territory where human health risk is more relevant and follow-up activities are required.

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

Modeling Low-temperature Geochemical Processes

NASA Astrophysics Data System (ADS)

Nordstrom, D. K.

2003-12-01

Geochemical modeling has become a popular and useful tool for a wide number of applications from research on the fundamental processes of water-rock interactions to regulatory requirements and decisions regarding permits for industrial and hazardous wastes. In low-temperature environments, generally thought of as those in the temperature range of 0-100 °C and close to atmospheric pressure (1 atm=1.01325 bar=101,325 Pa), complex hydrobiogeochemical reactions participate in an array of interconnected processes that affect us, and that, in turn, we affect. Understanding these complex processes often requires tools that are sufficiently sophisticated to portray multicomponent, multiphase chemical reactions yet transparent enough to reveal the main driving forces. Geochemical models are such tools. The major processes that they are required to model include mineral dissolution and precipitation; aqueous inorganic speciation and complexation; solute adsorption and desorption; ion exchange; oxidation-reduction; or redox; transformations; gas uptake or production; organic matter speciation and complexation; evaporation; dilution; water mixing; reaction during fluid flow; reaction involving biotic interactions; and photoreaction. These processes occur in rain, snow, fog, dry atmosphere, soils, bedrock weathering, streams, rivers, lakes, groundwaters, estuaries, brines, and diagenetic environments. Geochemical modeling attempts to understand the redistribution of elements and compounds, through anthropogenic and natural means, for a large range of scale from nanometer to global. "Aqueous geochemistry" and "environmental geochemistry" are often used interchangeably with "low-temperature geochemistry" to emphasize hydrologic or environmental objectives.Recognition of the strategy or philosophy behind the use of geochemical modeling is not often discussed or explicitly described. Plummer (1984, 1992) and Parkhurst and Plummer (1993) compare and contrast two approaches for modeling groundwater chemistry: (i) "forward modeling," which predicts water compositions from hypothesized reactions and user assumptions and (ii) "inverse modeling," which uses water, mineral, and isotopic compositions to constrain hypothesized reactions. These approaches simply reflect the amount of information one has to work with. With minimal information on a site, a modeler is forced to rely on forward modeling. Optimal information would include detailed mineralogy on drill cores or well cuttings combined with detailed water analyses at varying depths and sufficient spatial distribution to follow geochemical reactions and mixing of waters along defined flow paths. With optimal information, a modeler will depend on inverse modeling.This chapter outlines the main concepts and key developments in the field of geochemical modeling for low-temperature environments and illustrates their use with examples. It proceeds with a short discussion of what modeling is, continues with concepts and definitions commonly used, and follows with a short history of geochemical models, a discussion of databases, the codes that embody models, and recent examples of how these codes have been used in water-rock interactions. An important new stage of development seems to have been reached in this field with questions of reliability and validity of models. Future work will be obligated to document ranges of certainty and sources of uncertainty, sensitivity of models and codes to parameter errors and assumptions, propagation of errors, and delineation of the range of applicability.

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.

Distributed All-Optical Sensor to Detect dCO2 in Aqueous Environments

NASA Astrophysics Data System (ADS)

Bhatia, S.; Coelho, J.; Melo, L.; Davies, B.; Ahmed, F.; Bao, B.; Wild, P.; Risk, D. A.; Sinton, D.; Jun, M.

2012-12-01

Already a proven technology for temperature and stress measurements, an all-optical sensor to detect dCO2 is being developed for deployment in challenging environments. Optical sensors function under high pressure, do not require electronics and therefore experience no magnetic interference. They are also able to transmit signals over long distances with minimal losses. The dCO2 sensor's principal application is in measurement monitoring and verification of carbon capture and storage sites; however, it could also be useful in ocean, fresh water, and transition environments. The objective for the first phase of development was to detect a CO2 signal in laboratory tests. The developmental program incorporated experiments to detect CO2 under high pressure (1400 psi) in aqueous environments. Laboratory testing involved a custom pressure cell, off-the-shelf and custom long period gratings written in SMF125 fiber. Femptosecond laser micromachining was used to test alternative long period grating (LPG) and cutout shapes to maximize evanescent field interaction with the environment. A comprehensive program of geochemical modeling using PHREEQC 2 was used to identify the diversity of species in environments of interest that could exert confounding influence. Purchased UV-LPG responded to changes in concentration of scCO2 in brine at high pressure. Signal differences between CO2-saturated brine and pure brine were also noted under the same, high pressure conditions. Geochemical modeling software, PHREEQC 2, revealed a diversity of species in environments of interest whose concentrations varied strongly with temperature and pH. The modeling program's detailed characterization of environments informed work currently being undertaken as part of Phase 2, to develop a CO2-selective membrane to filter out measurement artifact.

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.

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.

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.

The site-scale saturated zone flow model for Yucca Mountain: Calibration of different conceptual models and their impact on flow paths

USGS Publications Warehouse

Zyvoloski, G.; Kwicklis, E.; Eddebbarh, A.-A.; Arnold, B.; Faunt, C.; Robinson, B.A.

2003-01-01

This paper presents several different conceptual models of the Large Hydraulic Gradient (LHG) region north of Yucca Mountain and describes the impact of those models on groundwater flow near the potential high-level repository site. The results are based on a numerical model of site-scale saturated zone beneath Yucca Mountain. This model is used for performance assessment predictions of radionuclide transport and to guide future data collection and modeling activities. The numerical model is calibrated by matching available water level measurements using parameter estimation techniques, along with more informal comparisons of the model to hydrologic and geochemical information. The model software (hydrologic simulation code FEHM and parameter estimation software PEST) and model setup allows for efficient calibration of multiple conceptual models. Until now, the Large Hydraulic Gradient has been simulated using a low-permeability, east-west oriented feature, even though direct evidence for this feature is lacking. In addition to this model, we investigate and calibrate three additional conceptual models of the Large Hydraulic Gradient, all of which are based on a presumed zone of hydrothermal chemical alteration north of Yucca Mountain. After examining the heads and permeabilities obtained from the calibrated models, we present particle pathways from the potential repository that record differences in the predicted groundwater flow regime. The results show that Large Hydraulic Gradient can be represented with the alternate conceptual models that include the hydrothermally altered zone. The predicted pathways are mildly sensitive to the choice of the conceptual model and more sensitive to the quality of calibration in the vicinity on the repository. These differences are most likely due to different degrees of fit of model to data, and do not represent important differences in hydrologic conditions for the different conceptual models. ?? 2002 Elsevier Science B.V. All rights reserved.

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.

Detectability of Wellbore CO2 Leakage using the Magnetotelluric Method

NASA Astrophysics Data System (ADS)

Yang, X.; Buscheck, T. A.; Mansoor, K.; Carroll, S.

2016-12-01

We assessed the effectiveness of the magnetotelluric (MT) method in detecting CO2 and brine leakage through a wellbore, which penetrates a CO2 storage reservoir, into overlying aquifers, 0 to 1720 m in depth, in support of the USDOE National Risk Assessment Partnership (NRAP) monitoring program. Synthetic datasets based on the Kimberlina site in the southern San Joaquin Basin, California were created using CO2 storage reservoir models, wellbore leakage models, and groundwater/geochemical models of the overlying aquifers. The species concentrations simulated with the groundwater/geochemical models were converted into bulk electrical conductivity (EC) distributions as the MT model input. Brine and CO2 leakage into the overlying aquifers increases ion concentrations, and thus results in an EC increase, which may be detected by the MT method. Our objective was to estimate and maximize the probability of leakage detection using the MT method. The MT method is an electromagnetic geophysical technique that images the subsurface EC distribution by measuring natural electric and magnetic fields in the frequency range from 0.01 Hz to 1 kHz with sensors on the ground surface. The ModEM software was used to predict electromagnetic responses from brine and CO2 leakage and to invert synthetic MT data for recovery of subsurface conductivity distribution. We are in the process of building 1000 simulations for ranges of permeability, leakage flux, and hydraulic gradient to study leakage detectability and to develop an optimization method to answer when, where and how an MT monitoring system should be deployed to maximize the probability of leakage detection. This work was sponsored by the USDOE Fossil Energy, National Energy Technology Laboratory, managed by Traci Rodosta and Andrea McNemar. This work was performed under the auspices of the USDOE by LLNL under contract DE-AC52-07NA27344. LLNL IM release number is LLNL-ABS-699276.

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.

NASA Tech Briefs, March 2012

NASA Technical Reports Server (NTRS)

2012-01-01

The topics include: 1) Spectral Profiler Probe for In Situ Snow Grain Size and Composition Stratigraphy; 2) Portable Fourier Transform Spectroscopy for Analysis of Surface Contamination and Quality Control; 3) In Situ Geochemical Analysis and Age Dating of Rocks Using Laser Ablation-Miniature Mass Spectrometer; 4) Physics Mining of Multi-Source Data Sets; 5) Photogrammetry Tool for Forensic Analysis; 6) Connect Global Positioning System RF Module; 7) Simple Cell Balance Circuit; 8) Miniature EVA Software Defined Radio; 9) Remotely Accessible Testbed for Software Defined Radio Development; 10) System-of-Systems Technology-Portfolio-Analysis Tool; 11) VESGEN Software for Mapping and Quantification of Vascular Regulators; 12) Constructing a Database From Multiple 2D Images for Camera Pose Estimation and Robot Localization; 13) Adaption of G-TAG Software for Validating Touch and Go Asteroid Sample Return Design Methodology; 14) 3D Visualization for Phoenix Mars Lander Science Operations; 15) RxGen General Optical Model Prescription Generator; 16) Carbon Nanotube Bonding Strength Enhancement Using Metal Wicking Process; 17) Multi-Layer Far-Infrared Component Technology; 18) Germanium Lift-Off Masks for Thin Metal Film Patterning; 19) Sealing Materials for Use in Vacuum at High Temperatures; 20) Radiation Shielding System Using a Composite of Carbon Nanotubes Loaded With Electropolymers; 21) Nano Sponges for Drug Delivery and Medicinal Applications; 22) Molecular Technique to Understand Deep Microbial Diversity; 23) Methods and Compositions Based on Culturing Microorganisms in Low Sedimental Fluid Shear Conditions; 24) Secure Peer-to-Peer Networks for Scientific Information Sharing; 25) Multiplexer/Demultiplexer Loading Tool (MDMLT); 26) High-Rate Data-Capture for an Airborne Lidar System; 27) Wavefront Sensing Analysis of Grazing Incidence Optical Systems; 28) Foam-on-Tile Damage Model; 29) Instrument Package Manipulation Through the Generation and Use of an Attenuated-Fluent Gas Fold; 30) Multicolor Detectors for Ultrasensitive Long-Wave Imaging Cameras; 31) Lunar Reconnaissance Orbiter (LRO) Command and Data Handling Flight Electronics Subsystem; and 32) Electro-Optic Segment-Segment Sensors for Radio and Optical Telescopes.

Building continental-scale 3D subsurface layers in the Digital Crust project: constrained interpolation and uncertainty estimation.

NASA Astrophysics Data System (ADS)

Yulaeva, E.; Fan, Y.; Moosdorf, N.; Richard, S. M.; Bristol, S.; Peters, S. E.; Zaslavsky, I.; Ingebritsen, S.

2015-12-01

The Digital Crust EarthCube building block creates a framework for integrating disparate 3D/4D information from multiple sources into a comprehensive model of the structure and composition of the Earth's upper crust, and to demonstrate the utility of this model in several research scenarios. One of such scenarios is estimation of various crustal properties related to fluid dynamics (e.g. permeability and porosity) at each node of any arbitrary unstructured 3D grid to support continental-scale numerical models of fluid flow and transport. Starting from Macrostrat, an existing 4D database of 33,903 chronostratigraphic units, and employing GeoDeepDive, a software system for extracting structured information from unstructured documents, we construct 3D gridded fields of sediment/rock porosity, permeability and geochemistry for large sedimentary basins of North America, which will be used to improve our understanding of large-scale fluid flow, chemical weathering rates, and geochemical fluxes into the ocean. In this talk, we discuss the methods, data gaps (particularly in geologically complex terrain), and various physical and geological constraints on interpolation and uncertainty estimation.

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.

Resolving the Multi-scale Behavior of Geochemical Weathering in the Critical Zone Using High Resolution Hydro-geochemical Models

NASA Astrophysics Data System (ADS)

Pandey, S.; Rajaram, H.

2015-12-01

This work investigates hydrologic and geochemical interactions in the Critical Zone (CZ) using high-resolution reactive transport modeling. Reactive transport models can be used to predict the response of geochemical weathering and solute fluxes in the CZ to changes in a dynamic environment, such as those pertaining to human activities and climate change in recent years. The scales of hydrology and geochemistry in the CZ range from days to eons in time and centimeters to kilometers in space. Here, we present results of a multi-dimensional, multi-scale hydro-geochemical model to investigate the role of subsurface heterogeneity on the formation of mineral weathering fronts in the CZ, which requires consideration of many of these spatio-temporal scales. The model is implemented using the reactive transport code PFLOTRAN, an open source subsurface flow and reactive transport code that utilizes parallelization over multiple processing nodes and provides a strong framework for simulating weathering in the CZ. The model is set up to simulate weathering dynamics in the mountainous catchments representative of the Colorado Front Range. Model parameters were constrained based on hydrologic, geochemical, and geophysical observations from the Boulder Creek Critical Zone Observatory (BcCZO). Simulations were performed in fractured rock systems and compared with systems of heterogeneous and homogeneous permeability fields. Tracer simulations revealed that the mean residence time of solutes was drastically accelerated as fracture density increased. In simulations that include mineral reactions, distinct signatures of transport limitations on weathering arose when discrete flow paths were included. This transport limitation was related to both advective and diffusive processes in the highly heterogeneous systems (i.e. fractured media and correlated random permeability fields with σlnk > 3). The well-known time-dependence of mineral weathering rates was found to be the most pronounced in the fractured systems, with a departure from the maximum system-averaged dissolution rate occurring after ~100 kyr followed by a gradual decrease in the reaction rate with time that persists beyond 104 kyr.

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.

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.

Hydro-geochemical modeling of subalpine urbanized area: geochemical characterization of the shallow and deep aquifers of the urban district of Como (first results).

NASA Astrophysics Data System (ADS)

Terrana, Silvia; Brunamonte, Fabio; Frascoli, Francesca; Ferrario, Maria Francesca; Michetti, Alessandro Maria; Pozzi, Andrea; Gambillara, Roberto; Binda, Gilberto

2016-04-01

One of the greatest environmental and social-economics threats is climate change. This topic, in the next few years, will have a significant impact on the availability of water resources of many regions. This is compounded by the strong anthropization of water systems that shows an intensification of conflicts for water resource exploitation. Therefore, it is necessary a sustainable manage of natural resources thorough knowledge of the hosting territories. The development of investigation and data processing methods are essential to reduce costs for the suitable use and protection of resources. Identify a sample area for testing the best approach is crucial. This research aims to find a valid methodology for the characterization, modeling and management of subalpine urban aquifers, and the urban district of Como appears perfect. The city of Como is located at the southern end of the western sector of Lake Como (N Italy). It is a coastal town, placed on a small alluvial plain, therefore in close communication with the lake. The plain is drained by two streams, which are presently artificially buried, and have an underground flow path in the urban section till the mouth. This city area, so, is suitable for this project as it is intensely urbanized, its dimensions is not too extensive and it is characterized by two aquifers very important and little known. These are a shallow aquifer and a deep aquifer, which are important not only for any water supply, but also for the stability of the ground subsidence in the city. This research is also the opportunity to work in a particular well-known area with high scientific significance; however, there is complete absence of information regarding the deep aquifer. Great importance has also the chosen and used of the more powerful open source software for this type of area, such as PHREEQC, EnvironInsite, PHREEQE etc., used for geological and geochemical data processing. The main goal of this preliminary work is the characterization of both aquifers and the comprehension of the interaction between the shallow one and deep one. The first results of geochemical and hydrological analyses in nine piezometer, seven in the shallow aquifer and two in the deep one, permit to identify a general characterization of groundwater: the waters of both the aquifers are calcium bicarbonate, with discontinuous enrichment in sulphate and silicate. The conductibility values are higher in the shallow aquifer then in the deep aquifer, this is an indication of the high flow rate of the deep one and of the higher level of vulnerability of the shallow water respect the urban pollutant. Particular is also the relatively high temperature for the deep water (15.5°C to 16.1°C). The next steps of work are an increase of geochemical analysis, with trace element, and the correlation with hydro-chemical surface water.

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.

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

On the Future of Thermochemical Databases, the Development of Solution Models and the Practical Use of Computational Thermodynamics in Volcanology, Geochemistry and Petrology: Can Innovations of Modern Data Science Democratize an Oligarchy?

NASA Astrophysics Data System (ADS)

Ghiorso, M. S.

2014-12-01

Computational thermodynamics (CT) has now become an essential tool of petrologic and geochemical research. CT is the basis for the construction of phase diagrams, the application of geothermometers and geobarometers, the equilibrium speciation of solutions, the construction of pseudosections, calculations of mass transfer between minerals, melts and fluids, and, it provides a means of estimating materials properties for the evaluation of constitutive relations in fluid dynamical simulations. The practical application of CT to Earth science problems requires data. Data on the thermochemical properties and the equation of state of relevant materials, and data on the relative stability and partitioning of chemical elements between phases as a function of temperature and pressure. These data must be evaluated and synthesized into a self consistent collection of theoretical models and model parameters that is colloquially known as a thermodynamic database. Quantitative outcomes derived from CT reply on the existence, maintenance and integrity of thermodynamic databases. Unfortunately, the community is reliant on too few such databases, developed by a small number of research groups, and mostly under circumstances where refinement and updates to the database lag behind or are unresponsive to need. Given the increasing level of reliance on CT calculations, what is required is a paradigm shift in the way thermodynamic databases are developed, maintained and disseminated. They must become community resources, with flexible and assessable software interfaces that permit easy modification, while at the same time maintaining theoretical integrity and fidelity to the underlying experimental observations. Advances in computational and data science give us the tools and resources to address this problem, allowing CT results to be obtained at the speed of thought, and permitting geochemical and petrological intuition to play a key role in model development and calibration.

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.

Geochemical and radionuclide profile of Tuzla geothermal field, Turkey.

PubMed

Baba, Alper; Deniz, Ozan; Ozcan, Hasan; Erees, Serap F; Cetiner, S Ziya

2008-10-01

Tuzla geothermal basin is situated in north-western Turkey on the Biga Peninsula, which is located at the west end of the Northern Anatolian Fault system. Soil and water samples were collected between August 2003 and June 2004 to initiate development of a geochemical profile of surface and subsurface waters in the geothermal basin and radionuclide concentrations in soils. All water samples were found to fall within Turkish Water Quality Class 4, meaning they were remarkably contaminated for any water consumption sector (industrial, human use or agricultural) based on sodium and chloride ions. Such waters could be used only after appropriate water treatment. The water samples are of the chloride type in terms of geochemical evaluation. Preliminary geochemical evidence shows that the N-S flowing part of the Tuzla River acts as a natural barrier within the basin. Heavy metal concentrations in the soil samples show slight elevations, especially those obtained from the east part of the basin where thermal springs are dominant. Geochemical calculations were carried out with PHREEQC software to determine equilibrium concentration of chemical species and saturation indices, by which it is suggested that chloride is the most important ligand to mobilize the heavy metals in the studied system. In addition, the activity concentration and gamma-absorbed dose rates of the terrestrial naturally occurring radionuclides were determined in the soil using gamma-ray spectrometry. The soil activity ranged from 42.77 to 988.66 Bq kg(-1) (averaging 138 Bq kg(-1)) for ( 238 )U, 13.27 to 106.31 Bq kg(-1) (averaging 32.42 Bq kg(-1)) for ( 232 )Th, and 99.28 to 935.36 Bq kg(-1) (averaging 515.44 Bq kg(-1)) for ( 40 )K. The highest value of ( 238 )U was found in the soil samples obtained from an area close to the hot spring.

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.

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.

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.

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.

BAT3 Analyzer: Real-Time Data Display and Interpretation Software for the Multifunction Bedrock-Aquifer Transportable Testing Tool (BAT3)

USGS Publications Warehouse

Winston, Richard B.; Shapiro, Allen M.

2007-01-01

The BAT3 Analyzer provides real-time display and interpretation of fluid pressure responses and flow rates measured during geochemical sampling, hydraulic testing, or tracer testing conducted with the Multifunction Bedrock-Aquifer Transportable Testing Tool (BAT3) (Shapiro, 2007). Real-time display of the data collected with the Multifunction BAT3 allows the user to ensure that the downhole apparatus is operating properly, and that test procedures can be modified to correct for unanticipated hydraulic responses during testing. The BAT3 Analyzer can apply calibrations to the pressure transducer and flow meter data to display physically meaningful values. Plots of the time-varying data can be formatted for a specified time interval, and either saved to files, or printed. Libraries of calibrations for the pressure transducers and flow meters can be created, updated and reloaded to facilitate the rapid set up of the software to display data collected during testing with the Multifunction BAT3. The BAT3 Analyzer also has the functionality to estimate calibrations for pressure transducers and flow meters using data collected with the Multifunction BAT3 in conjunction with corroborating check measurements. During testing with the Multifunction BAT3, and also after testing has been completed, hydraulic properties of the test interval can be estimated by comparing fluid pressure responses with model results; a variety of hydrogeologic conceptual models of the formation are available for interpreting fluid-withdrawal, fluid-injection, and slug tests.

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.

The application of remote sensing techniques to the study of ophiolites

NASA Astrophysics Data System (ADS)

Khan, Shuhab D.; Mahmood, Khalid

2008-08-01

Satellite remote sensing methods are a powerful tool for detailed geologic analysis, especially in inaccessible regions of the earth's surface. Short-wave infrared (SWIR) bands are shown to provide spectral information bearing on the lithologic, structural, and geochemical character of rock bodies such as ophiolites, allowing for a more comprehensive assessment of the lithologies present, their stratigraphic relationships, and geochemical character. Most remote sensing data are widely available for little or no cost, along with user-friendly software for non-specialists. In this paper we review common remote sensing systems and methods that allow for the discrimination of solid rock (lithologic) components of ophiolite complexes and their structural relationships. Ophiolites are enigmatic rock bodies which associated with most, if not all, plate collision sutures. Ophiolites are ideal for remote sensing given their widely recognized diversity of lithologic types and structural relationships. Accordingly, as a basis for demonstrating the utility of remote sensing techniques, we briefly review typical ophiolites in the Tethyan tectonic belt. As a case study, we apply integrated remote sensing studies of a well-studied example, the Muslim Bagh ophiolite, located in Balochistan, western Pakistan. On this basis, we attempt to demonstrate how remote sensing data can validate and reconcile existing information obtained from field studies. The lithologic and geochemical diversity of Muslim Bagh are representative of Tethyan ophiolites. Despite it's remote location it has been extensively mapped and characterized by structural and geochemical studies, and is virtually free of vegetative cover. Moreover, integrating the remote sensing data with 'ground truth' information thus offers the potential of an improved template for interpreting remote sensing data sets of other ophiolites for which little or no field information is available.

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

Mineralogical determination and geo-chemical modeling of chromium release from AOD slag: Distribution and leachability aspects.

PubMed

Li, Junguo; Liu, Bao; Zeng, Yanan; Wang, Ziming

2017-01-01

AOD (argon oxygen decarburization) slag, which is the by-product of the stainless steel refining process, is a recyclable slag because of its high content of calcium and silicon. The leaching toxicity cannot be ignored in the recycling process because the slag contains a certain amount of Cr. In this study, the mineral analysis, batch leaching tests and thermodynamic and kinetic modeling by PHREEQC combined with FactSage software were performed to explore the influence of the dissolution of primary minerals and the precipitation of secondary minerals on the elution of Cr from AOD slag. The results indicated that the main minerals in the original AOD slag are larnite, merwinite, pyroxene and periclase. Cr was dispersed in the mineral phases mentioned above. The simulation of Cr leaching controlled by Cr(III)-hydroxide corresponded better to the batch leaching tests, while the Cr leaching controlled by chromite or double control was underestimated. Increasing the L/S ratio enhances the pH of the leachate and restrains the elution of Cr from the AOD slag. Copyright © 2016 Elsevier Ltd. All rights reserved.

Automating U-Pb IDTIMS data reduction and reporting: Cyberinfrastructure meets geochronology

NASA Astrophysics Data System (ADS)

Bowring, J. F.; McLean, N.; Walker, J. D.; Ash, J. M.

2009-12-01

We demonstrate the efficacy of an interdisciplinary effort between software engineers and geochemists to produce working cyberinfrastructure for geochronology. This collaboration between CIRDLES, EARTHTIME and EarthChem has produced the software programs Tripoli and U-Pb_Redux as the cyber-backbone for the ID-TIMS community. This initiative incorporates shared isotopic tracers, data-reduction algorithms and the archiving and retrieval of data and results. The resulting system facilitates detailed inter-laboratory comparison and a new generation of cooperative science. The resolving power of geochronological data in the earth sciences is dependent on the precision and accuracy of many isotopic measurements and corrections. Recent advances in U-Pb geochronology have reinvigorated its application to problems such as precise timescale calibration, processes of crustal evolution, and early solar system dynamics. This project provides a heretofore missing common data reduction protocol, thus promoting the interpretation of precise geochronology and enabling inter-laboratory comparison. U-Pb_Redux is an open-source software program that provides end-to-end support for the analysis of uranium-lead geochronological data. The system reduces raw mass spectrometer data to U-Pb dates, allows users to interpret ages from these data, and then provides for the seamless federation of the results, coming from many labs, into a community web-accessible database using standard and open techniques. This EarthChem GeoChron database depends also on keyed references to the SESAR sample database. U-Pb_Redux currently provides interactive concordia and weighted mean plots and uncertainty contribution visualizations; it produces publication-quality concordia and weighted mean plots and customizable data tables. This initiative has achieved the goal of standardizing the data elements of a complete reduction and analysis of uranium-lead data, which are expressed using extensible markup language schema definition (XSD) artifacts. U-Pb_Redux leverages the freeware program Tripoli, which imports raw mass spectrometer data files and supports interactive review and archiving of isotopic data. Tripoli facilitates the visualization of temporal trends and scatter during measurement, statistically rigorous filtering of data and supports oxide and fractionation corrections. The Cyber Infrastructure Research and Development Lab for the Earth Sciences (CIRDLES) collaboratively integrates domain-specific software engineering with the efforts EARTHTIME and Earthchem. The EARTHTIME initiative pursues consensus-based approaches to geochemical data reduction, and the EarthChem initiative pursues the creation of data repositories for all geochemical data. CIRDLES develops software and systems for geochronology. This collaboration benefits the earth sciences by enabling geochemists to focus on their specialties using robust software that produces reliable results. This collaboration benefits software engineering by providing research opportunities to improve process methodologies used in the design and implementation of domain-specific solutions.

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.

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.

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

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

Influence of physical factors and geochemical conditions on groundwater acidification during enhanced reductive dechlorination

NASA Astrophysics Data System (ADS)

Brovelli, A.; Barry, D. A.; Robinson, C.; Gerhard, J.

2010-12-01

Enhanced reductive dehalogenation is an attractive in situ treatment technology for chlorinated contaminants. The process includes two acid-forming microbial reactions: fermentation of an organic substrate resulting in short-chain fatty acids, and dehalogenation resulting in hydrochloric acid. The accumulation of acids and the resulting drop of groundwater pH are controlled by the mass and distribution of chlorinated solvents in the source zone, type of electron donor, availability of alternative terminal electron acceptors and presence of soil mineral phases able to buffer the pH (such as carbonates). Groundwater acidification may reduce or halt microbial activity, and thus dehalogenation, significantly increasing the time and costs required to remediate the aquifer. For this reason, research in this area is gaining increasing attention. In previous work (Robinson et al., 2009 407:4560, Sci. Tot. Environ, Robinson and Barry, 2009 24:1332, Environ. Model. & Software, Brovelli et al., 2010, submitted), a detailed geochemical and groundwater flow model able to predict the pH change occurring during reductive dehalogenation was developed. The model accounts for the main processes influencing groundwater pH, including the groundwater composition, the electron donor used and soil mineral phase interactions. In this study, the model was applied to investigate how spatial variability occurring at the field scale affects groundwater pH and dechlorination rates. Numerical simulations were conducted to examine the influence of heterogeneous hydraulic conductivity on the distribution of the injected, fermentable substrate and on the accumulation/dilution of the acidic products of reductive dehalogenation. The influence of the geometry of the DNAPL source zone was studied, as well as the spatial distribution of soil minerals. The results of this study showed that the heterogeneous distribution of the soil properties have a potentially large effect on the remediation efficiency. For example, zones of high hydraulic conductivity can prevent the accumulation of acids and alleviate the problem of groundwater acidification. The conclusions drawn and insights gained from this modeling study will be useful to design improved in situ enhanced dehalogenation remediation schemes.

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

Training the next generation of scientists: Modeling Infectious Disease and Water Quality of Montana Streams

NASA Astrophysics Data System (ADS)

Fytilis, N.; Wyman, S.; Lamb, R.; Stevens, L.; Kerans, B.; Rizzo, D. M.

2010-12-01

The University of Vermont College of Engineering and Mathematical Sciences and the Barrett Foundation have established a scholarship program for undergraduate students. The Barrett Scholarship program, aware of the importance of developing research quantitative and writing skills for undergraduate students, provides scholarships to outstanding undergraduate students for environmental engineering research projects. The intent is to help retain student interest early in their undergraduate engineering careers when few of their first or second year classes have little engineering or real-world application. We focus on one Barrett research project, derived from a NSF Biodiversity and Infectious Disease grant, because of the multiple disciplines (engineering, ecology, biology) and education levels (spanning secondary to graduate) involved. In this research, students across three departments at two universities (University of Vermont, Montana State University) and one independent high school (Vermont Commons School) formed a cohesive collaboration with faculty members to identify different worm taxa of T. Tubifex. Whirling disease has had a severe impact on the native population of salmonids in the upper Madison River MT, USA, resulting in the death of most fish that contract the parasite. T. Tubifex is the intermediate host for Myxobolus cerebralis, the parasite that causes whirling disease in salmonids. Samples collected from eight locations along the Madison River varied in the prevalence of whirling disease. The site-specific worm community structure has been measured and identified using molecular genetic probes and a taxonomic key to link worm communities to geochemical features (e.g. site elevation, slope, pH, conductivity, temperature, dissolved oxygen and percent of organic soil matter). Using a unique clustering algorithm, we group geochemical features to discriminate over a range of water quality gradients (i.e., “clean” to “dirty”). The link between water quality and the presence of these taxa is important in determining stream health. In addition, system dynamics software (STELLA) is used to model the non-linear relationships and feedback between worm prevalence and disease dynamics. These types of collaborations between engineers, biologists, field ecologists and geneticists from secondary, post-secondary and higher institutions proved useful in linking complex geochemical data, worm community structure and molecular genetics to develop the next-generation scientists and better understand disease dynamics.

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.

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.

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.

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.

NETPATH-WIN: an interactive user version of the mass-balance model, NETPATH

USGS Publications Warehouse

El-Kadi, A. I.; Plummer, Niel; Aggarwal, P.

2011-01-01

NETPATH-WIN is an interactive user version of NETPATH, an inverse geochemical modeling code used to find mass-balance reaction models that are consistent with the observed chemical and isotopic composition of waters from aquatic systems. NETPATH-WIN was constructed to migrate NETPATH applications into the Microsoft WINDOWS® environment. The new version facilitates model utilization by eliminating difficulties in data preparation and results analysis of the DOS version of NETPATH, while preserving all of the capabilities of the original version. Through example applications, the note describes some of the features of NETPATH-WIN as applied to adjustment of radiocarbon data for geochemical reactions in groundwater systems.

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.

Software for analysis of chemical mixtures--composition, occurrence, distribution, and possible toxicity

USGS Publications Warehouse

Scott, Jonathon C.; Skach, Kenneth A.; Toccalino, Patricia L.

2013-01-01

The composition, occurrence, distribution, and possible toxicity of chemical mixtures in the environment are research concerns of the U.S. Geological Survey and others. The presence of specific chemical mixtures may serve as indicators of natural phenomena or human-caused events. Chemical mixtures may also have ecological, industrial, geochemical, or toxicological effects. Chemical-mixture occurrences vary by analyte composition and concentration. Four related computer programs have been developed by the National Water-Quality Assessment Program of the U.S. Geological Survey for research of chemical-mixture compositions, occurrences, distributions, and possible toxicities. The compositions and occurrences are identified for the user-supplied data, and therefore the resultant counts are constrained by the user’s choices for the selection of chemicals, reporting limits for the analytical methods, spatial coverage, and time span for the data supplied. The distribution of chemical mixtures may be spatial, temporal, and (or) related to some other variable, such as chemical usage. Possible toxicities optionally are estimated from user-supplied benchmark data. The software for the analysis of chemical mixtures described in this report is designed to work with chemical-analysis data files retrieved from the U.S. Geological Survey National Water Information System but can also be used with appropriately formatted data from other sources. Installation and usage of the mixture software are documented. This mixture software was designed to function with minimal changes on a variety of computer-operating systems. To obtain the software described herein and other U.S. Geological Survey software, visit http://water.usgs.gov/software/.

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.

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.

Extensive decarbonation of continuously hydrated subducting slabs

NASA Astrophysics Data System (ADS)

Arzilli, Fabio; Burton, Mike; La Spina, Giuseppe; Macpherson, Colin G.

2017-04-01

CO2 release from subducting slabs is a key element of Earth's carbon cycle, consigning slab carbon either to mantle burial or recycling to the surface through arc volcanism, however, what controls subducted carbon's fate is poorly understood. Fluids mobilized by devolatilization of subducting slabs play a fundamental role in the melting of mantle wedges and in global geochemical cycles [1]. The effect of such fluids on decarbonation in subducting lithologies has been investigated recently [2-5], but several thermodynamic models [2-3], and experimental studies [6] suggest that carbon-bearing phases are stable at sub-arc depths (80-140 km; 2.6-4.5 GPa), implying that this carbon can be carried to mantle depths of >140 km. This is inconsistent with observations of voluminous CO2 release from arc volcanoes [7-10], located above slabs that are at 2.6-4.5 GPa pressure. The aim of this study is to re-evaluate the role of metamorphic decarbonation, showing if decarbonation reactions could be feasible at sub-arc depths combined with a continuous hydration scenario. We used the PerpleX software combined with a custom-designed algorithm to simulate a pervasive fluid infiltration characterized by "continuous hydration" combined with a distillation model, in which is possible to remove CO2 when decarbonation occurs, to obtain an open-system scenario. This is performed by repeatedly flushing the sediment with pure H2O at 0.5, 1.0 or 5 wt.% until no further decarbonation occurs. Here we show that continuous hydrated of sediment veneers on subducting slabs by H2O released from oceanic crust and serpentinised mantle lithosphere [11-13], produces extensive slab decarbonation over a narrow, sub-arc pressure range, even for low temperature subduction pathways. This explains the location of CO2-rich volcanism, quantitatively links the sedimentary composition of slab material to the degree of decarbonation and greatly increases estimates for the magnitude of carbon flux through the arc in subduction zones. [1] Hilton, D.R. et al. (2002) Rev. Mineral. Geochem. 47, 319-370. [2] Gorman, P.J. et al. (2006) Geochem. Geophys. Geosyst. 7. [3] Kerrick, D.M. and Connolly, J.A.D. (2001) Nature 411, 293-296. [4] Cook-Kollars, J. et al. (2014) Chem. Geol. 386, 31-48. [5] Collins, N.C. et al. (2015) Chem. Geol. 412, 132-150. [6] Poli, S. et al. (2009) Earth Planet. Sci. Lett. 278, 350-360. [7] Sano, Y. and Williams, S.N. (1996) Geophys. Res. Lett. 23, 2749-2752. [8] Marty, B. and Tolstikhin, I.N. (1998) Chem. Geol. 145, 233-248. [9] Wallace, P.J. (2005) J. Volcanol. Geoth. Res. 140, 217-240. [10] Burton, M.R. et al. (2013) Rev. Mineral. Geochem. 75, 323-354. [11] Ulmer, P. and Trommsdorff, V. (1995) Science 268, 858-861. [12] Schmidt, M.W. and Poli, S. (1998) Earth Planet. Sci. Lett. 163, 361-379. [13] van Keken, P. E. et al. (2011) J. Geophys. Res. 116.

Index map of cross sections through parts of the Appalachian basin (Kentucky, New York, Ohio, Pennsylvania, Tennessee, Virginia, West Virginia): Chapter E.1 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Ryder, Robert T.; Trippi, Michael H.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

The appendixes in chapters E.4.1 and E.4.2 include (1) Log ASCII Standard (LAS) files, which encode gamma-ray, neutron, density, and other logs in text files that can be used by most well-logging software programs; and (2) graphic well-log traces. In the appendix to chapter E.4.1, the well-log traces are accompanied by lithologic descriptions with formation tops.

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.

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.

Zircon age-temperature-compositional spectra in plutonic rocks

DOE PAGES

Samperton, Kyle M.; Bell, Elizabeth A.; Barboni, Mélanie; ...

2017-08-23

We present that geochronology can resolve dispersed zircon dates in plutonic rocks when magma cooling time scales exceed the temporal precision of individual U-Pb analyses; such age heterogeneity may indicate protracted crystallization between the temperatures of zircon saturation (T sat) and rock solidification (T solid). Diffusive growth models predict asymmetric distributions of zircon dates and crystallization temperatures in a cooling magma, with volumetrically abundant old, hot crystallization at T sat decreasing continuously to volumetrically minor young, cold crystallization at T solid. We present integrated geochronological and geochemical data from Bergell Intrusion tonalites (Central Alps, Europe) that document zircon compositional changemore » over hundreds of thousands of years at the hand-sample scale, indicating melt compositional evolution during solidification. Ti-in-zircon thermometry, crystallization simulation using MELTS software, and U-Pb dates produce zircon mass-temperature-time distributions that are in excellent agreement with zircon growth models. In conclusion, these findings provide the first quantitative validation of longstanding expectations from zircon saturation theory by direct geochronological investigation, underscoring zircon’s capacity to quantify supersolidus cooling rates in magmas and resolve dynamic differentiation histories in the plutonic rock record.« less

Zircon age-temperature-compositional spectra in plutonic rocks

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

Samperton, Kyle M.; Bell, Elizabeth A.; Barboni, Mélanie

We present that geochronology can resolve dispersed zircon dates in plutonic rocks when magma cooling time scales exceed the temporal precision of individual U-Pb analyses; such age heterogeneity may indicate protracted crystallization between the temperatures of zircon saturation (T sat) and rock solidification (T solid). Diffusive growth models predict asymmetric distributions of zircon dates and crystallization temperatures in a cooling magma, with volumetrically abundant old, hot crystallization at T sat decreasing continuously to volumetrically minor young, cold crystallization at T solid. We present integrated geochronological and geochemical data from Bergell Intrusion tonalites (Central Alps, Europe) that document zircon compositional changemore » over hundreds of thousands of years at the hand-sample scale, indicating melt compositional evolution during solidification. Ti-in-zircon thermometry, crystallization simulation using MELTS software, and U-Pb dates produce zircon mass-temperature-time distributions that are in excellent agreement with zircon growth models. In conclusion, these findings provide the first quantitative validation of longstanding expectations from zircon saturation theory by direct geochronological investigation, underscoring zircon’s capacity to quantify supersolidus cooling rates in magmas and resolve dynamic differentiation histories in the plutonic rock record.« less

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

Integrating Geochemical Reactions with a Particle-Tracking Approach to Simulate Nitrogen Transport and Transformation in Aquifers

NASA Astrophysics Data System (ADS)

Cui, Z.; Welty, C.; Maxwell, R. M.

2011-12-01

Lagrangian, particle-tracking models are commonly used to simulate solute advection and dispersion in aquifers. They are computationally efficient and suffer from much less numerical dispersion than grid-based techniques, especially in heterogeneous and advectively-dominated systems. Although particle-tracking models are capable of simulating geochemical reactions, these reactions are often simplified to first-order decay and/or linear, first-order kinetics. Nitrogen transport and transformation in aquifers involves both biodegradation and higher-order geochemical reactions. In order to take advantage of the particle-tracking approach, we have enhanced an existing particle-tracking code SLIM-FAST, to simulate nitrogen transport and transformation in aquifers. The approach we are taking is a hybrid one: the reactive multispecies transport process is operator split into two steps: (1) the physical movement of the particles including the attachment/detachment to solid surfaces, which is modeled by a Lagrangian random-walk algorithm; and (2) multispecies reactions including biodegradation are modeled by coupling multiple Monod equations with other geochemical reactions. The coupled reaction system is solved by an ordinary differential equation solver. In order to solve the coupled system of equations, after step 1, the particles are converted to grid-based concentrations based on the mass and position of the particles, and after step 2 the newly calculated concentration values are mapped back to particles. The enhanced particle-tracking code is capable of simulating subsurface nitrogen transport and transformation in a three-dimensional domain with variably saturated conditions. Potential application of the enhanced code is to simulate subsurface nitrogen loading to the Chesapeake Bay and its tributaries. Implementation details, verification results of the enhanced code with one-dimensional analytical solutions and other existing numerical models will be presented in addition to a discussion of implementation challenges.

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

A modeling study of the long-term mineral trapping in deep saline marine sands aquifers (Invited)

NASA Astrophysics Data System (ADS)

Aagaard, P.; Pham, V.; Hellevang, H.

2009-12-01

Simulation of geochemical processes due to CO2 injection and storage are dependent on sediment petrography and the kinetics of mineral fluid reactions. Mineral trapping of CO2 in the Utsira sand and similar marine sand reservoirs have been revisited based on critical review of rate data and geochemical constraints on formation waters. Reaction paths calculations were done with the PHREEQC modeling software at relevant reservoir conditions covering a temperature range of 30-100 °C and corresponding reservoir pressures. Initial CO2 saturation was determined by the fluid fugacity corresponding with reservoir conditions. The mineral dissolution kinetics was expressed with a chemical affinity term (Aagaard & Helgeson,1982) while a critical super-saturation for mineral growth was included in the precipitation rate expression. The redox conditions and the H2S fugacity in the simulations were constrained by the acetic/propionic acid buffer trend and the magnetite-pyrite buffer (Aagaard et al. 2001) respectively. We used a revised mineralogical composition for the Utsira sand also performed a sensitivity analyses with respect to mineral content. The simulations were run over a period of 10000 years. The main simulation results included dissolution of glauconite, smectite, pyrite, muscovite and albite, with precipitation of the carbonates siderite, ankerite, and minor dawsonite, as well as kaolinite, silica (either chalcedony or quartz), and K-feldspar. The uncertainties in the simulations are specially connected with initial mineral abundances. The effect of critical super-saturation and reactive surface area for precipitation needs to be further evaluated and tested. Aagaard, P. and H.C. Helgeson (1982). Thermodynamic and Kinetic Constraints on Reaction Rates among Minerals and Aqueous Solutions. I. Theoretical Considerations. Am. J. Sci., v. 282, p. 257-285. P. Aagaard, J. Jahren & S.N. Ehrenberg (2001) H2S controling reactions in clastic hydrocarbon reservoirs from the Norwegian Shelf and Gulf Coast, in Cidu, R.(ed) Water-Rock Interaction, WRI-10, Balkema, p. 129-132.

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.

Hydrogeochemical processes and geochemical modeling in a coastal aquifer: Case study of the Marathon coastal plain, Greece

NASA Astrophysics Data System (ADS)

Papazotos, Panagiotis; Koumantakis, Ioannis; Kallioras, Andreas; Vasileiou, Eleni; Perraki, Maria

2017-04-01

Determining the hydrogeochemical processes has always been a challenge for scientists. The aim of this work is the study of the principal hydrogeochemical processes controlling groundwater quality in the Marathon coastal plain, Greece, with emphasis on the origin of the solutes. Various physicochemical parameters and major ions of twenty-five groundwater samples were analyzed. The hydrogeochemical data of groundwater were studied in order to determine the major factors controlling the chemical composition and hydrogeochemical evolution. In the Marathon coastal plain, three different zones of the alluvial granular aquifer system have been detected, considering the geochemical processes and recharge, which affect its hydrochemical characteristics. The alluvial granular aquifer system is divided eastwards into three zones: a) the natural recharge zone, b) the reverse ion exchange zone and c) the diffusion sea water zone. Cl-is the dominant anion and Na+and Ca2+ are the dominant cations, as determined by plotting the analyses on the respective Piper diagram. Near the coastline high concentrations of Na+ and Cl- were observed indicating a zone of seawater intrusion. On the other hand, westward there is increasing concentration of HCO3- with simultaneous decrease of Na+is indication of a recharge zone from karstic aquifers of the study area. Between the aforementioned zones there is an intermediate one, where reverse ion exchange takes place due to high concentrations of dissolved Na+ and Ca2+ adsorption. The saturation indices (SI) were calculated using the geochemical modeling software PHREEQC. Mineral phases of halite, sylvite, gypsum and anhydrite were estimated to be undersaturated in the water samples, suggesting these phases are minor or absent in the host rock. On the other hand, calcite, aragonite and dolomite are close to equilibrium; these minerals are present in the host rocks or in the unsaturated zone, possibly increasing the Ca2+, Mg2+ and HCO3- concentrations when carbonates are dissolved. The analyses of the bivariate scatter plots, the ionic ratios, the Indices of Base Exchange (IBE), the Gibbs diagram and the dissolution/precipitation reactions show that evaporation and water-rock interaction mechanisms such as dissolution of carbonates, followed by reverse ion exchange, have affected the groundwater chemistry in the study area. The results revealed that groundwater chemistry and therefore the origin of the solutes in the coastal alluvial granular aquifer system of the Marathon coastal plain is primarily affected by a number of factors such as groundwater and mineral equilibrium, seawater intrusion, reverse ion exchange and nitrate concentration. A possible future research could focus on the interaction among hydrogeochemistry, mineral phases and chemical thermodynamic modeling.

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.

Gold potential in the Dalradian rocks of NW Northern Ireland: GIS-based prospectivity analysis using Tellus data

NASA Astrophysics Data System (ADS)

Lusty, P. A. J.; McDonnell, P. M.; Gunn, A. G.; Chacksfield, B. C.; Cooper, M.

2009-04-01

Geographic Information Systems (GIS) are essential tools for the management and integration of the large amounts of multivariate spatial data used in mineral exploration. Prospectivity analysis combines these datasets, in the context of a mineral deposit model, to produce a map showing the distribution of potential for a particular type of mineral deposit. In this example Arc-Spatial Data Modeller software has been used to analyse the prospectivity for orogenic vein gold mineralisation in the Dalradian rocks of north-western Northern Ireland. A knowledge-driven (fuzzy logic) approach was used because of the small number of gold deposits within the area. Fuzzy logic is used in situations where information is inexact and the use of classical set theory is inappropriate. Fuzzy logic allows assignment of weightings to exploration data on a continuous scale from 1 (full membership) to 0 (full non-membership). This allows a level of uncertainty or 'fuzziness' to be incorporated into the modelling. The key stages of prospectivity analysis are: (1) analysis of the deposit model to determine key exploration indicators; (2) data processing, interpretation and analysis to extract key indicators; (3) assignment of weightings, zones and styles of influence to key indicators; and (4) calculation of prospectivity. This research is based largely on new geochemical and geophysical data resulting from the Tellus Project in Northern Ireland. The Tellus Project involved geochemical and airborne geophysical surveys over the whole of Northern Ireland carried out between 2004-6 with funding from the Government of Northern Ireland. The study area (3074 km2) is underlain mainly by Neoproterozoic rocks of the Dalradian Supergroup (ca. 590 Ma) which form part of the Caledonide orogenic belt. The Dalradian Supergroup comprises a thick succession of semi-pelites, psammites and pelites, with graphitic pelite horizons that host much of the known gold mineralisation. In the Sperrin Mountains two advanced gold projects, Curraghinalt and Golan Burn, and an operating mine, Omagh (Cavanacaw), are hosted by Dalradian rocks. The Dalradian rocks of Northern Ireland were extensively deformed during the Grampian Orogeny and at least four phases of deformation are recognised. On a deposit scale distinct differences are observed between Curraghinalt and Cavanacaw including orientation of principal structures, mineralogy and geochemistry. The mineral prospectivity analysis integrated a range of datasets including: 1:250 000 scale geological mapping, mineral occurrences, fault vectors, Tellus drainage geochemistry and geophysics. In order to assess the differences between the Curraghinalt and Cavanacaw mineral deposit models separate prospectivity models were parameterised. The key differences between the two models relate to the weightings attached to the structural vectors and the significance given to particular chemical elements in the drainage geochemical dataset. The prospectivity analysis identified several areas prospective for orogenic-style gold mineralisation. Many of the prospective zones either coincide with known occurrences of orogenic gold mineralisation or areas considered highly prospective based upon previous work, validating the process and the model results. The use of specific models for the Curraghinalt and Cavanacaw deposits has produced significant differences in the location and extent of the target areas which provide a basis for focusing exploration for orogenic gold deposits in this region. Detailed examination of all available data and ground truthing is required in order to validate these targets. In addition to delineating new areas favourable for gold mineralisation the prospectivity mapping has provided new insight into possible regional controls on the location of mineralisation and into the geology of this area. Please fill in your abstract text.

Direct measurement of neon production rates by (α,n) reactions in minerals

NASA Astrophysics Data System (ADS)

Cox, Stephen E.; Farley, Kenneth A.; Cherniak, Daniele J.

2015-01-01

The production of nucleogenic neon from alpha particle capture by 18O and 19F offers a potential chronometer sensitive to temperatures higher than the more widely used (U-Th)/He chronometer. The accuracy depends on the cross sections and the calculated stopping power for alpha particles in the mineral being studied. Published 18O(α,n)21Ne production rates are in poor agreement and were calculated from contradictory cross sections, and therefore demand experimental verification. Similarly, the stopping powers for alpha particles are calculated from SRIM (Stopping Range of Ions in Matter software) based on a limited experimental dataset. To address these issues we used a particle accelerator to implant alpha particles at precisely known energies into slabs of synthetic quartz (SiO2) and barium tungstate (BaWO4) to measure 21Ne production from capture by 18O. Within experimental uncertainties the observed 21Ne production rates compare favorably to our predictions using published cross sections and stopping powers, indicating that ages calculated using these quantities are accurate at the ∼3% level. In addition, we measured the 22Ne/21Ne ratio and (U-Th)/He and (U-Th)/Ne ages of Durango fluorapatite, which is an important model system for this work because it contains both oxygen and fluorine. Finally, we present 21Ne/4He production rate ratios for a variety of minerals of geochemical interest along with software for calculating neon production rates and (U-Th)/Ne ages.

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.