Complete suite of geochemical values computed using wireline logs

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

Lancaster, J.R.; Atkinson, A.

1996-12-31

Geochemical values of {open_quotes}black shale{close_quotes} source rocks can be computed from a complete suite of wireline log data. The computed values are: Total Organic Carbon (Wt%). S1, S2, S3, Hydrogen Index, Oxygen Index, Atomic H/C and O/C ratios, Genetic Potential (S1+S2), S2/S3, and Transfomation Ratio (S1/(S1+S2)). The results are most reliable when calibrated to laboratory analyses of samples in the study area. However, in the absence of samples, reasonable estimates can be made using calibration data from analogous depositional and thermal environments and/or professional judgement and experience. The evaluations provide answers to critical geochemical questions relative to: (1) Organic Mattermore » Quantity; T.O.C. (Wt%), S1, and S2. (2) Kerogen Types; I, II, and III, based on T.O.C. vs S2 cross plot and the van Krevelen diagram of Atomic O/C vs Atomic H/C ratios. (3) Thermal Maturation levels; Transfomation Ratio can be converted to Level of Organic Metamorphism (LOM), pyrolysis Tmax (degC), Vitrinite Reflectance (Ro), Time Temperature Index (TTI) and others. Various analog plots and cross plots can be prepared for interpretation. Case history examples are shown and discussed. Lowstand fan deposits on Barbados were studied in outcrop to construct a conceptual reservoir model for prediction of facies assemblages.« less

Complete suite of geochemical values computed using wireline logs

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

Lancaster, J.R.; Atkinson, A.

1996-01-01

Geochemical values of [open quotes]black shale[close quotes] source rocks can be computed from a complete suite of wireline log data. The computed values are: Total Organic Carbon (Wt%). S1, S2, S3, Hydrogen Index, Oxygen Index, Atomic H/C and O/C ratios, Genetic Potential (S1+S2), S2/S3, and Transfomation Ratio (S1/(S1+S2)). The results are most reliable when calibrated to laboratory analyses of samples in the study area. However, in the absence of samples, reasonable estimates can be made using calibration data from analogous depositional and thermal environments and/or professional judgement and experience. The evaluations provide answers to critical geochemical questions relative to: (1)more » Organic Matter Quantity; T.O.C. (Wt%), S1, and S2. (2) Kerogen Types; I, II, and III, based on T.O.C. vs S2 cross plot and the van Krevelen diagram of Atomic O/C vs Atomic H/C ratios. (3) Thermal Maturation levels; Transfomation Ratio can be converted to Level of Organic Metamorphism (LOM), pyrolysis Tmax (degC), Vitrinite Reflectance (Ro), Time Temperature Index (TTI) and others. Various analog plots and cross plots can be prepared for interpretation. Case history examples are shown and discussed. Lowstand fan deposits on Barbados were studied in outcrop to construct a conceptual reservoir model for prediction of facies assemblages.« less

An alternative to fully coupled reactive transport simulations for long-term prediction of chemical reactions in complex geological systems

NASA Astrophysics Data System (ADS)

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

2014-05-01

Fully-coupled reactive transport simulations involving multiphase hydrodynamics and chemical reactions in heterogeneous settings are extremely challenging from a computational point of view. This often leads to oversimplification of the investigated system: coarse spatial discretization, to keep the number of elements in the order of few thousands; simplified chemistry, disregarding many potentially important reactions. A novel approach for coupling non-reactive hydrodynamic simulations with the outcome of single batch geochemical simulations was therefore introduced to assess the potential long-term mineral trapping at the Ketzin pilot site for underground CO2 storage in Germany [1],[2]. The advantage of the coupling is the ability to use multi-million grid non-reactive hydrodynamics simulations on one side and few batch 0D geochemical simulations on the other, so that the complexity of both systems does not need to be reduced. This contribution shows the approach which was taken to validate this simplified coupling scheme. The procedure involved batch simulations of the reference geochemical model, then performing both non-reactive and fully coupled 1D and 3D reactive transport simulations and finally applying the simplified coupling scheme based on the non-reactive and geochemical batch model. The TOUGHREACT/ECO2N [3] simulator was adopted for the validation. The degree of refinement of the spatial grid and the complexity and velocity of the mineral reactions, along with a cut-off value for the minimum concentration of dissolved CO2 allowed to originate precipitates in the simplified approach were found out to be the governing parameters for the convergence of the two schemes. Systematic discrepancies between the approaches are not reducible, simply because there is no feedback between chemistry and hydrodynamics, and can reach 20 % - 30 % in unfavourable cases. However, even such discrepancy is completely acceptable, in our opinion, given the amount of uncertainty underlying the geochemical models. References [1] Klein, E., De Lucia, M., Kempka, T. Kühn, M. 2013. Evaluation of longterm mineral trapping at the Ketzin pilot site for CO2 storage: an integrative approach using geochemical modelling and reservoir simulation. International Journal of Greenhouse Gas Control 19: 720-730, doi:10.1016/j.ijggc.2013.05.014 [2] Kempka, T., Klein, E., De Lucia, M., Tillner, E. Kühn, M. 2013. Assessment of Long-term CO2 Trapping Mechanisms at the Ketzin Pilot Site (Germany) by Coupled Numerical Modelling. Energy Procedia 37: 5419-5426, doi:10.1016/j.egypro.2013.06.460 [3] Xu, T., Spycher, N., Sonnenthal, E., Zhang, G., Zheng, L., Pruess, K. 2010. TOUGHREACT Version 2.0: A simulator for subsurface reactive transport under non-isothermal multiphase flow conditions, Computers & Geosciences 37(6), doi:10.1016/j.cageo.2010.10.007

SeSBench - An initiative to benchmark reactive transport models for environmental subsurface processes

NASA Astrophysics Data System (ADS)

Jacques, Diederik

2017-04-01

As soil functions are governed by a multitude of interacting hydrological, geochemical and biological processes, simulation tools coupling mathematical models for interacting processes are needed. Coupled reactive transport models are a typical example of such coupled tools mainly focusing on hydrological and geochemical coupling (see e.g. Steefel et al., 2015). Mathematical and numerical complexity for both the tool itself or of the specific conceptual model can increase rapidly. Therefore, numerical verification of such type of models is a prerequisite for guaranteeing reliability and confidence and qualifying simulation tools and approaches for any further model application. In 2011, a first SeSBench -Subsurface Environmental Simulation Benchmarking- workshop was held in Berkeley (USA) followed by four other ones. The objective is to benchmark subsurface environmental simulation models and methods with a current focus on reactive transport processes. The final outcome was a special issue in Computational Geosciences (2015, issue 3 - Reactive transport benchmarks for subsurface environmental simulation) with a collection of 11 benchmarks. Benchmarks, proposed by the participants of the workshops, should be relevant for environmental or geo-engineering applications; the latter were mostly related to radioactive waste disposal issues - excluding benchmarks defined for pure mathematical reasons. Another important feature is the tiered approach within a benchmark with the definition of a single principle problem and different sub problems. The latter typically benchmarked individual or simplified processes (e.g. inert solute transport, simplified geochemical conceptual model) or geometries (e.g. batch or one-dimensional, homogeneous). Finally, three codes should be involved into a benchmark. The SeSBench initiative contributes to confidence building for applying reactive transport codes. Furthermore, it illustrates the use of those type of models for different environmental and geo-engineering applications. SeSBench will organize new workshops to add new benchmarks in a new special issue. Steefel, C. I., et al. (2015). "Reactive transport codes for subsurface environmental simulation." Computational Geosciences 19: 445-478.

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.

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

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

USGS Publications Warehouse

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

1994-01-01

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

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.

Hydrological and geochemical processes constraining groundwater salinity in wetland areas related to evaporitic (karst) systems. A case study from Southern Spain

NASA Astrophysics Data System (ADS)

Gil-Márquez, J. M.; Barberá, J. A.; Andreo, B.; Mudarra, M.

2017-01-01

Chemical and isotopic evolution of groundwater in an evaporite karst plateau (including wetland areas and saline to hyper-saline springs) located at S Spain was studied. Physicochemical parameters, major ions and stable isotopes were analyzed in rain, brine spring, wetland and leakage water samples, from which the most common mineral saturation indexes were computed and geochemical and isotopic modelling were performed. Results show an apparent relationship between the elevation of brine springs and their water mineralization, indicating that drainage at higher altitude may be associated to gravity-driven flows, since brackish groundwater is isotopically fractionated due to evaporation. On the other hand, the lower altitude springs could drain deeper flows with longer residence time, resulting in highly mineralized and warmer (briny) groundwater. The dissolution of halite and gypsum has proved to be the main geochemical processes, which are favored by the great ionic strength of groundwater. Calcite precipitation occurs in brackish waters draining wetlands, being boosted by common ion effect (when CaSO4 waters are present) and solute concentration caused by evaporation. Modelling results strongly support the hypothesis that most of the selected springs geochemically evolve in a common (S-N) flowpath. The methods used in this research contribute to a better understanding of the hydrogeological processes occurring in the studied evaporitic system, but also in equivalent hydrological environments worldwide.

Geomechanical/Geochemical Modeling Studies Conducted within theInternational DECOVALEX Project

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

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

2005-10-19

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

Geomechanical/ Geochemical Modeling Studies onducted Within the International DECOVALEX Project

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

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

2006-02-01

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

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

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.

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

PubMed

Bartzas, Georgios; Komnitsas, Kostas

2010-11-15

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

Modelling radionuclide transport in fractured media with a dynamic update of K d values

DOE PAGES

Trinchero, Paolo; Painter, Scott L.; Ebrahimi, Hedieh; ...

2015-10-13

Radionuclide transport in fractured crystalline rocks is a process of interest in evaluating long term safety of potential disposal systems for radioactive wastes. Given their numerical efficiency and the absence of numerical dispersion, Lagrangian methods (e.g. particle tracking algorithms) are appealing approaches that are often used in safety assessment (SA) analyses. In these approaches, many complex geochemical retention processes are typically lumped into a single parameter: the distribution coefficient (Kd). Usually, the distribution coefficient is assumed to be constant over the time frame of interest. However, this assumption could be critical under long-term geochemical changes as it is demonstrated thatmore » the distribution coefficient depends on the background chemical conditions (e.g. pH, Eh, and major chemistry). In this study, we provide a computational framework that combines the efficiency of Lagrangian methods with a sound and explicit description of the geochemical changes of the site and their influence on the radionuclide retention properties.« less

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

A flow-simulation model of the tidal Potomac River

USGS Publications Warehouse

Schaffranek, Raymond W.

1987-01-01

A one-dimensional model capable of simulating flow in a network of interconnected channels has been applied to the tidal Potomac River including its major tributaries and embayments between Washington, D.C., and Indian Head, Md. The model can be used to compute water-surface elevations and flow discharges at any of 66 predetermined locations or at any alternative river cross sections definable within the network of channels. In addition, the model can be used to provide tidal-interchange flow volumes and to evaluate tidal excursions and the flushing properties of the riverine system. Comparisons of model-computed results with measured watersurface elevations and discharges demonstrate the validity and accuracy of the model. Tidal-cycle flow volumes computed by the calibrated model have been verified to be within an accuracy of ? 10 percent. Quantitative characteristics of the hydrodynamics of the tidal river are identified and discussed. The comprehensive flow data provided by the model can be used to better understand the geochemical, biological, and other processes affecting the river's water quality.

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

NASA Astrophysics Data System (ADS)

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

2014-08-01

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

PHAST Version 2-A Program for Simulating Groundwater Flow, Solute Transport, and Multicomponent Geochemical Reactions

USGS Publications Warehouse

Parkhurst, David L.; Kipp, Kenneth L.; Charlton, Scott R.

2010-01-01

The computer program PHAST (PHREEQC And HST3D) simulates multicomponent, reactive solute transport in three-dimensional saturated groundwater flow systems. PHAST is a versatile groundwater flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. Major enhancements in PHAST Version 2 allow spatial data to be defined in a combination of map and grid coordinate systems, independent of a specific model grid (without node-by-node input). At run time, aquifer properties are interpolated from the spatial data to the model grid; regridding requires only redefinition of the grid without modification of the spatial data. PHAST is applicable to the study of natural and contaminated groundwater systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock/water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, or density-dependent flow. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux (specified-flux), and leaky (head-dependent) conditions, as well as the special cases of rivers, drains, and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association or Pitzer specific interaction thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, ion exchange sites, surface complexation sites, solid solutions, and gases; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, exchangers, surfaces, gases, kinetic reactants, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and postprocessing programs; and in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST.

Data collection handbook to support modeling the impacts of radioactive material in soil

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

Yu, C.; Cheng, J.J.; Jones, L.G.

1993-04-01

A pathway analysis computer code called RESRAD has been developed for implementing US Department of Energy Residual Radioactive Material Guidelines. Hydrogeological, meteorological, geochemical, geometrical (size, area, depth), and material-related (soil, concrete) parameters are used in the RESRAD code. This handbook discusses parameter definitions, typical ranges, variations, measurement methodologies, and input screen locations. Although this handbook was developed primarily to support the application of RESRAD, the discussions and values are valid for other model applications.

Combining particle-tracking and geochemical data to assess public supply well vulnerability to arsenic and uranium

USGS Publications Warehouse

Hinkle, S.R.; Kauffman, L.J.; Thomas, M.A.; Brown, C.J.; McCarthy, K.A.; Eberts, S.M.; Rosen, Michael R.; Katz, B.G.

2009-01-01

Flow-model particle-tracking results and geochemical data from seven study areas across the United States were analyzed using three statistical methods to test the hypothesis that these variables can successfully be used to assess public supply well vulnerability to arsenic and uranium. Principal components analysis indicated that arsenic and uranium concentrations were associated with particle-tracking variables that simulate time of travel and water fluxes through aquifer systems and also through specific redox and pH zones within aquifers. Time-of-travel variables are important because many geochemical reactions are kinetically limited, and geochemical zonation can account for different modes of mobilization and fate. Spearman correlation analysis established statistical significance for correlations of arsenic and uranium concentrations with variables derived using the particle-tracking routines. Correlations between uranium concentrations and particle-tracking variables were generally strongest for variables computed for distinct redox zones. Classification tree analysis on arsenic concentrations yielded a quantitative categorical model using time-of-travel variables and solid-phase-arsenic concentrations. The classification tree model accuracy on the learning data subset was 70%, and on the testing data subset, 79%, demonstrating one application in which particle-tracking variables can be used predictively in a quantitative screening-level assessment of public supply well vulnerability. Ground-water management actions that are based on avoidance of young ground water, reflecting the premise that young ground water is more vulnerable to anthropogenic contaminants than is old ground water, may inadvertently lead to increased vulnerability to natural contaminants due to the tendency for concentrations of many natural contaminants to increase with increasing ground-water residence time.

Exploration computer applications to primary dispersion halos: Kougarok tin prospect, Seward Peninsula, Alaska, USA

USGS Publications Warehouse

Reid, Jeffrey C.

1989-01-01

Computer processing and high resolution graphics display of geochemical data were used to quickly, accurately, and efficiently obtain important decision-making information for tin (cassiterite) exploration, Seward Peninsula, Alaska (USA). Primary geochemical dispersion patterns were determined for tin-bearing intrusive granite phases of Late Cretaceous age with exploration bedrock lithogeochemistry at the Kougarok tin prospect. Expensive diamond drilling footage was required to reach exploration objectives. Recognition of element distribution and dispersion patterns was useful in subsurface interpretation and correlation, and to aid location of other holes.

A geochemical module for "AMDTreat" to compute caustic quantity, effluent quantity, and sludge volume

USGS Publications Warehouse

Cravotta, Charles A.; Parkhurst, David L.; Means, Brent P; McKenzie, Bob; Morris, Harry; Arthur, Bill

2010-01-01

Treatment with caustic chemicals typically is used to increase pH and decrease concentrations of dissolved aluminum, iron, and/or manganese in largevolume, metal-laden discharges from active coal mines. Generally, aluminum and iron can be removed effectively at near-neutral pH (6 to 8), whereas active manganese removal requires treatment to alkaline pH (~10). The treatment cost depends on the specific chemical used (NaOH, CaO, Ca(OH)2, Na2CO3, or NH3) and increases with the quantities of chemical added and sludge produced. The pH and metals concentrations do not change linearly with the amount of chemical added. Consequently, the amount of caustic chemical needed to achieve a target pH and the corresponding effluent composition and sludge volume can not be accurately determined without empirical titration data or the application of geochemical models to simulate the titration of the discharge water with caustic chemical(s). The AMDTreat computer program (http://amd.osmre.gov/ ) is widely used to compute costs for treatment of coal-mine drainage. Although AMDTreat can use results of empirical titration with industrial grade caustic chemicals to compute chemical costs for treatment of net-acidic or net-alkaline mine drainage, such data are rarely available. To improve the capability of AMDTreat to estimate (1) the quantity and cost of caustic chemicals to attain a target pH, (2) the concentrations of dissolved metals in treated effluent, and (3) the volume of sludge produced by the treatment, a titration simulation is being developed using the geochemical program PHREEQC (wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/) that will be coupled as a module to AMDTreat. The simulated titration results can be compared with or used in place of empirical titration data to estimate chemical quantities and costs. This paper describes the development, evaluation, and potential utilization of the PHREEQC titration module for AMDTreat.

Biogeochemical metabolic modeling of methanogenesis by Methanosarcina barkeri

NASA Astrophysics Data System (ADS)

Jensvold, Z. D.; Jin, Q.

2015-12-01

Methanogenesis, the biological process of methane production, is the final step of natural organic matter degradation. In studying natural methanogenesis, important questions include how fast methanogenesis proceeds and how methanogens adapt to the environment. To address these questions, we propose a new approach - biogeochemical reaction modeling - by simulating the metabolic networks of methanogens. Biogeochemical reaction modeling combines geochemical reaction modeling and genome-scale metabolic modeling. Geochemical reaction modeling focuses on the speciation of electron donors and acceptors in the environment, and therefore the energy available to methanogens. Genome-scale metabolic modeling predicts microbial rates and metabolic strategies. Specifically, this approach describes methanogenesis using an enzyme network model, and computes enzyme rates by accounting for both the kinetics and thermodynamics. The network model is simulated numerically to predict enzyme abundances and rates of methanogen metabolism. We applied this new approach to Methanosarcina barkeri strain fusaro, a model methanogen that makes methane by reducing carbon dioxide and oxidizing dihydrogen. The simulation results match well with the results of previous laboratory experiments, including the magnitude of proton motive force and the kinetic parameters of Methanosarcina barkeri. The results also predict that in natural environments, the configuration of methanogenesis network, including the concentrations of enzymes and metabolites, differs significantly from that under laboratory settings.

PHAST--a program for simulating ground-water flow, solute transport, and multicomponent geochemical reactions

USGS Publications Warehouse

Parkhurst, David L.; Kipp, Kenneth L.; Engesgaard, Peter; Charlton, Scott R.

2004-01-01

The computer program PHAST simulates multi-component, reactive solute transport in three-dimensional saturated ground-water flow systems. PHAST is a versatile ground-water flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. PHAST is applicable to the study of natural and contaminated ground-water systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock-water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, density-dependent flow, or waters with high ionic strengths. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux, and leaky conditions, as well as the special cases of rivers and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, gases, surface complexation sites, ion exchange sites, and solid solutions; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, gases, exchangers, surfaces, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and post-processing programs; or in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST. Operator splitting of the flow, transport, and geochemical equations is used to separate the three processes into three sequential calculations. No iterations between transport and reaction calculations are implemented. A three-dimensional Cartesian coordinate system and finite-difference techniques are used for the spatial and temporal discretization of the flow and transport equations. The non-linear chemical equilibrium equations are solved by a Newton-Raphson method, and the kinetic reaction equations are solved by a Runge-Kutta or an implicit method for integrating ordinary differential equations. The PHAST simulator may require large amounts of memory and long Central Processing Unit (CPU) times. To reduce the long CPU times, a parallel version of PHAST has been developed that runs on a multiprocessor computer or on a collection of computers that are networked. The parallel version requires Message Passing Interface, which is currently (2004) freely available. The parallel version is effective in reducing simulation times. This report documents the use of the PHAST simulator, including running the simulator, preparing the input files, selecting the output files, and visualizing the results. It also presents four examples that verify the numerical method and demonstrate the capabilities of the simulator. PHAST requires three input files. Only the flow and transport file is described in detail in this report. The other two files, the chemistry data file and the database file, are identical to PHREEQC files and the detailed description of these files is found in the PHREEQC documentation.

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

USGS Publications Warehouse

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

1988-01-01

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

A Comparison of Methods for Modeling Geochemical Variability in the Earth's Mantle

NASA Astrophysics Data System (ADS)

Kellogg, J. B.; Tackley, P. J.

2004-12-01

Numerial models of isotopic and chemical heterogeneity of the Earth's mantle fall into three categories, in decreasing order of computational demand. First, several authors have used chemical tracers within a full thermo-chemical convection calculation (e.g., Christensen and Hofmann, 1994, van Keken and Ballentine, 1999; Xie and Tackley, 2004). Second, Kellogg et al. (2002) proposed an extension of the traditional geochemical box model calculations in which numerous subreservoirs were tracked within the bulk depleted mantle reservoir. Third, Allègre and Lewin (1995) described a framework in which the variance in chemical and isotopic ratios were treated as quantities intrinsic to the bulk reservoirs, complete with sources and sinks. Results from these three methods vary, particularly with respect to conclusions drawn about the meaning of the Pb-Pb pseudo-isochron. We revisit these methods in an attempt to arrive at a common understanding. By considering all three we better identify the strengths and weaknesses of each approach and allow each to inform the other. Finally, we present results from a new hybrid model that combines the complexity and regional-scale variability of the thermochemical convection models with the short length-scale sensitivity of the Kellogg et al. approach.

Application of a new model for groundwater age distributions: Modeling and isotopic analysis of artificial recharge in the Rialto-Colton basin, California

USGS Publications Warehouse

Ginn, T.R.; Woolfenden, L.

2002-01-01

A project for modeling and isotopic analysis of artificial recharge in the Rialto-Colton basin aquifer in California, is discussed. The Rialto-Colton aquifer has been divided into four primary and significant flowpaths following the general direction of groundwater flow from NW to SE. The introductory investigation include sophisticated chemical reaction modeling, with highly simplified flow path simulation. A comprehensive reactive transport model with the established set of geochemical reactions over the whole aquifer will also be developed for treating both reactions and transport realistically. This will be completed by making use of HBGC123D implemented with isotopic calculation step to compute Carbon-14 (C14) and stable Carbon-13 (C13) contents of the water. Computed carbon contents will also be calibrated with the measured carbon contents for assessment of the amount of imported recharge into the Linden pond.

User's guide to PHREEQC (Version 2) : a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations

USGS Publications Warehouse

Parkhurst, David L.; Appelo, C.A.J.

1999-01-01

PHREEQC version 2 is a computer program written in the C programming language that is designed to perform a wide variety of low-temperature aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations; (2) batch-reaction and one-dimensional (1D) transport calculations involving reversible reactions, which include aqueous, mineral, gas, solid-solution, surface-complexation, and ion-exchange equilibria, and irreversible reactions, which include specified mole transfers of reactants, kinetically controlled reactions, mixing of solutions, and temperature changes; and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for differences in composition between waters, within specified compositional uncertainty limits.New features in PHREEQC version 2 relative to version 1 include capabilities to simulate dispersion (or diffusion) and stagnant zones in 1D-transport calculations, to model kinetic reactions with user-defined rate expressions, to model the formation or dissolution of ideal, multicomponent or nonideal, binary solid solutions, to model fixed-volume gas phases in addition to fixed-pressure gas phases, to allow the number of surface or exchange sites to vary with the dissolution or precipitation of minerals or kinetic reactants, to include isotope mole balances in inverse modeling calculations, to automatically use multiple sets of convergence parameters, to print user-defined quantities to the primary output file and (or) to a file suitable for importation into a spreadsheet, and to define solution compositions in a format more compatible with spreadsheet programs. This report presents the equations that are the basis for chemical equilibrium, kinetic, transport, and inverse-modeling calculations in PHREEQC; describes the input for the program; and presents examples that demonstrate most of the program's capabilities.

Hydrologic and geochemical approaches for determining ground-water flow components

USGS Publications Warehouse

Hjalmarson, H.W.; Robertson, F.N.

1991-01-01

Lyman Lake is an irrigation-storage reservoir on the Little Colorado River near St. Johns, Arizona. The main sources of water for the lake are streamflow in the Little Colorado River and ground-water inflow from the underlying Coconino aquifer. Two approaches, a hydrologic analysis and a geochemical analysis, were used to compute the quantity of ground-water flow to and from Lyman Lake. Hydrologic data used to calculate a water budget were precipitation on the lake, evaporation from the lake, transpiration from dense vegetation, seepage through the dam, streamflow in and out of the lake, and changes in lake storage. Geochemical data used to calculate the ground-water flow components were major ions, trace elements, and the stable isotopes of hydrogen and oxygen. During the study, the potentiometric level of the Coconino aquifer was above the lake level at the upstream end of the lake and below the lake level at the downstream end. Hydrologic and geochemical data indicate that about 10 percent and 8 percent, respectively, of the water in the lake is ground-water inflow and that about 35 percent of the water in the Little Colorado River 6 miles downgradient from the lake near Salado Springs is ground water. These independent estimates of ground-water flow derived from each approach are in agreement and support a conceptual model of the water budget.

Application of artificial neural networks to chemostratigraphy

NASA Astrophysics Data System (ADS)

Malmgren, BjöRn A.; Nordlund, Ulf

1996-08-01

Artificial neural networks, a branch of artificial intelligence, are computer systems formed by a number of simple, highly interconnected processing units that have the ability to learn a set of target vectors from a set of associated input signals. Neural networks learn by self-adjusting a set of parameters, using some pertinent algorithm to minimize the error between the desired output and network output. We explore the potential of this approach in solving a problem involving classification of geochemical data. The data, taken from the literature, are derived from four late Quaternary zones of volcanic ash of basaltic and rhyolithic origin from the Norwegian Sea. These ash layers span the oxygen isotope zones 1, 5, 7, and 11, respectively (last 420,000 years). The data consist of nine geochemical variables (oxides) determined in each of 183 samples. We employed a three-layer back propagation neural network to assess its efficiency to optimally differentiate samples from the four ash zones on the basis of their geochemical composition. For comparison, three statistical pattern recognition techniques, linear discriminant analysis, the k-nearest neighbor (k-NN) technique, and SIMCA (soft independent modeling of class analogy), were applied to the same data. All of these showed considerably higher error rates than the artificial neural network, indicating that the back propagation network was indeed more powerful in correctly classifying the ash particles to the appropriate zone on the basis of their geochemical composition.

Mineral and Geochemical Classification From Spectroscopy/Diffraction Through Neural Networks

NASA Astrophysics Data System (ADS)

Ferralis, N.; Grossman, J.; Summons, R. E.

2017-12-01

Spectroscopy and diffraction techniques are essential for understanding structural, chemical and functional properties of geological materials for Earth and Planetary Sciences. Beyond data collection, quantitative insight relies on experimentally assembled, or computationally derived spectra. Inference on the geochemical or geophysical properties (such as crystallographic order, chemical functionality, elemental composition, etc.) of a particular geological material (mineral, organic matter, etc.) is based on fitting unknown spectra and comparing the fit with consolidated databases. The complexity of fitting highly convoluted spectra, often limits the ability to infer geochemical characteristics, and limits the throughput for extensive datasets. With the emergence of heuristic approaches to pattern recognitions though machine learning, in this work we investigate the possibility and potential of using supervised neural networks trained on available public spectroscopic database to directly infer geochemical parameters from unknown spectra. Using Raman, infrared spectroscopy and powder x-ray diffraction from the publicly available RRUFF database, we train neural network models to classify mineral and organic compounds (pure or mixtures) based on crystallographic structure from diffraction, chemical functionality, elemental composition and bonding from spectroscopy. As expected, the accuracy of the inference is strongly dependent on the quality and extent of the training data. We will identify a series of requirements and guidelines for the training dataset needed to achieve consistent high accuracy inference, along with methods to compensate for limited of data.

Pattern recognition analysis and classification modeling of selenium-producing areas

USGS Publications Warehouse

Naftz, D.L.

1996-01-01

Established chemometric and geochemical techniques were applied to water quality data from 23 National Irrigation Water Quality Program (NIWQP) study areas in the Western United States. These techniques were applied to the NIWQP data set to identify common geochemical processes responsible for mobilization of selenium and to develop a classification model that uses major-ion concentrations to identify areas that contain elevated selenium concentrations in water that could pose a hazard to water fowl. Pattern recognition modeling of the simple-salt data computed with the SNORM geochemical program indicate three principal components that explain 95% of the total variance. A three-dimensional plot of PC 1, 2 and 3 scores shows three distinct clusters that correspond to distinct hydrochemical facies denoted as facies 1, 2 and 3. Facies 1 samples are distinguished by water samples without the CaCO3 simple salt and elevated concentrations of NaCl, CaSO4, MgSO4 and Na2SO4 simple salts relative to water samples in facies 2 and 3. Water samples in facies 2 are distinguished from facies 1 by the absence of the MgSO4 simple salt and the presence of the CaCO3 simple salt. Water samples in facies 3 are similar to samples in facies 2, with the absence of both MgSO4 and CaSO4 simple salts. Water samples in facies 1 have the largest selenium concentration (10 ??gl-1), compared to a median concentration of 2.0 ??gl-1 and less than 1.0 ??gl-1 for samples in facies 2 and 3. A classification model using the soft independent modeling by class analogy (SIMCA) algorithm was constructed with data from the NIWQP study areas. The classification model was successful in identifying water samples with a selenium concentration that is hazardous to some species of water-fowl from a test data set comprised of 2,060 water samples from throughout Utah and Wyoming. Application of chemometric and geochemical techniques during data synthesis analysis of multivariate environmental databases from other national-scale environmental programs such as the NIWQP could also provide useful insights for addressing 'real world' environmental problems.

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.

Numerical investigation of coupled density-driven flow and hydrogeochemical processes below playas

NASA Astrophysics Data System (ADS)

Hamann, Enrico; Post, Vincent; Kohfahl, Claus; Prommer, Henning; Simmons, Craig T.

2015-11-01

Numerical modeling approaches with varying complexity were explored to investigate coupled groundwater flow and geochemical processes in saline basins. Long-term model simulations of a playa system gain insights into the complex feedback mechanisms between density-driven flow and the spatiotemporal patterns of precipitating evaporites and evolving brines. Using a reactive multicomponent transport model approach, the simulations reproduced, for the first time in a numerical study, the evaporite precipitation sequences frequently observed in saline basins ("bull's eyes"). Playa-specific flow, evapoconcentration, and chemical divides were found to be the primary controls for the location of evaporites formed, and the resulting brine chemistry. Comparative simulations with the computationally far less demanding surrogate single-species transport models showed that these were still able to replicate the major flow patterns obtained by the more complex reactive transport simulations. However, the simulated degree of salinization was clearly lower than in reactive multicomponent transport simulations. For example, in the late stages of the simulations, when the brine becomes halite-saturated, the nonreactive simulation overestimated the solute mass by almost 20%. The simulations highlight the importance of the consideration of reactive transport processes for understanding and quantifying geochemical patterns, concentrations of individual dissolved solutes, and evaporite evolution.

Quantification of CO2 generation in sedimentary basins through carbonate/clays reactions with uncertain thermodynamic parameters

NASA Astrophysics Data System (ADS)

Ceriotti, G.; Porta, G. M.; Geloni, C.; Dalla Rosa, M.; Guadagnini, A.

2017-09-01

We develop a methodological framework and mathematical formulation which yields estimates of the uncertainty associated with the amounts of CO2 generated by Carbonate-Clays Reactions (CCR) in large-scale subsurface systems to assist characterization of the main features of this geochemical process. Our approach couples a one-dimensional compaction model, providing the dynamics of the evolution of porosity, temperature and pressure along the vertical direction, with a chemical model able to quantify the partial pressure of CO2 resulting from minerals and pore water interaction. The modeling framework we propose allows (i) estimating the depth at which the source of gases is located and (ii) quantifying the amount of CO2 generated, based on the mineralogy of the sediments involved in the basin formation process. A distinctive objective of the study is the quantification of the way the uncertainty affecting chemical equilibrium constants propagates to model outputs, i.e., the flux of CO2. These parameters are considered as key sources of uncertainty in our modeling approach because temperature and pressure distributions associated with deep burial depths typically fall outside the range of validity of commonly employed geochemical databases and typically used geochemical software. We also analyze the impact of the relative abundancy of primary phases in the sediments on the activation of CCR processes. As a test bed, we consider a computational study where pressure and temperature conditions are representative of those observed in real sedimentary formation. Our results are conducive to the probabilistic assessment of (i) the characteristic pressure and temperature at which CCR leads to generation of CO2 in sedimentary systems, (ii) the order of magnitude of the CO2 generation rate that can be associated with CCR processes.

A Metastable Equilibrium Model for the Relative Abundances of Microbial Phyla in a Hot Spring

PubMed Central

Dick, Jeffrey M.; Shock, Everett L.

2013-01-01

Many studies link the compositions of microbial communities to their environments, but the energetics of organism-specific biomass synthesis as a function of geochemical variables have rarely been assessed. We describe a thermodynamic model that integrates geochemical and metagenomic data for biofilms sampled at five sites along a thermal and chemical gradient in the outflow channel of the hot spring known as “Bison Pool” in Yellowstone National Park. The relative abundances of major phyla in individual communities sampled along the outflow channel are modeled by computing metastable equilibrium among model proteins with amino acid compositions derived from metagenomic sequences. Geochemical conditions are represented by temperature and activities of basis species, including pH and oxidation-reduction potential quantified as the activity of dissolved hydrogen. By adjusting the activity of hydrogen, the model can be tuned to closely approximate the relative abundances of the phyla observed in the community profiles generated from BLAST assignments. The findings reveal an inverse relationship between the energy demand to form the proteins at equal thermodynamic activities and the abundance of phyla in the community. The distance from metastable equilibrium of the communities, assessed using an equation derived from energetic considerations that is also consistent with the information-theoretic entropy change, decreases along the outflow channel. Specific divergences from metastable equilibrium, such as an underprediction of the relative abundances of phototrophic organisms at lower temperatures, can be explained by considering additional sources of energy and/or differences in growth efficiency. Although the metabolisms used by many members of these communities are driven by chemical disequilibria, the results support the possibility that higher-level patterns of chemotrophic microbial ecosystems are shaped by metastable equilibrium states that depend on both the composition of biomass and the environmental conditions. PMID:24023738

A rapid local singularity analysis algorithm with applications

NASA Astrophysics Data System (ADS)

Chen, Zhijun; Cheng, Qiuming; Agterberg, Frits

2015-04-01

The local singularity model developed by Cheng is fast gaining popularity in characterizing mineralization and detecting anomalies of geochemical, geophysical and remote sensing data. However in one of the conventional algorithms involving the moving average values with different scales is time-consuming especially while analyzing a large dataset. Summed area table (SAT), also called as integral image, is a fast algorithm used within the Viola-Jones object detection framework in computer vision area. Historically, the principle of SAT is well-known in the study of multi-dimensional probability distribution functions, namely in computing 2D (or ND) probabilities (area under the probability distribution) from the respective cumulative distribution functions. We introduce SAT and it's variation Rotated Summed Area Table in the isotropic, anisotropic or directional local singularity mapping in this study. Once computed using SAT, any one of the rectangular sum can be computed at any scale or location in constant time. The area for any rectangular region in the image can be computed by using only 4 array accesses in constant time independently of the size of the region; effectively reducing the time complexity from O(n) to O(1). New programs using Python, Julia, matlab and C++ are implemented respectively to satisfy different applications, especially to the big data analysis. Several large geochemical and remote sensing datasets are tested. A wide variety of scale changes (linear spacing or log spacing) for non-iterative or iterative approach are adopted to calculate the singularity index values and compare the results. The results indicate that the local singularity analysis with SAT is more robust and superior to traditional approach in identifying anomalies.

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.

A computer program for geochemical analysis of acid-rain and other low-ionic-strength, acidic waters

USGS Publications Warehouse

Johnsson, P.A.; Lord, D.G.

1987-01-01

ARCHEM, a computer program written in FORTRAN 77, is designed primarily for use in the routine geochemical interpretation of low-ionic-strength, acidic waters. On the basis of chemical analyses of the water, and either laboratory or field determinations of pH, temperature, and dissolved oxygen, the program calculates the equilibrium distribution of major inorganic aqueous species and of inorganic aluminum complexes. The concentration of the organic anion is estimated from the dissolved organic concentration. Ionic ferrous iron is calculated from the dissolved oxygen concentration. Ionic balances and comparisons of computed with measured specific conductances are performed as checks on the analytical accuracy of chemical analyses. ARCHEM may be tailored easily to fit different sampling protocols, and may be run on multiple sample analyses. (Author 's abstract)

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

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

Chen, Jinsong

2013-05-01

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

Experimental design applications for modeling and assessing carbon dioxide sequestration in saline aquifers

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

Rogers, John

2014-11-29

This project was a computer modeling effort to couple reservoir simulation and ED/RSM using Sensitivity Analysis, Uncertainty Analysis, and Optimization Methods, to assess geologic, geochemical, geomechanical, and rock-fluid effects and factors on CO 2 injectivity, capacity, and plume migration. The project objective was to develop proxy models to simplify the highly complex coupled geochemical and geomechanical models in the utilization and storage of CO 2 in the subsurface. The goals were to investigate and prove the feasibility of the ED/RSM processes and engineering development, and bridge the gaps regarding the uncertainty and unknowns of the many geochemical and geomechanical interactingmore » parameters in the development and operation of anthropogenic CO 2 sequestration and storage sites. The bottleneck in this workflow is the high computational effort of reactive transport simulation models and large number of input variables to optimize with ED/RSM techniques. The project was not to develop the reactive transport, geomechanical, or ED/RSM software, but was to use what was commercially and/or publically available as a proof of concept to generate proxy or surrogate models. A detailed geologic and petrographic mineral assemblage and geologic structure of the doubly plunging anticline was defined using the USDOE RMOTC formations of interest data (e.g., Lower Sundance, Crow Mountain, Alcova Limestone, and Red Peak). The assemblage of 23 minerals was primarily developed from literature data and petrophysical (well log) analysis. The assemblage and structure was input into a commercial reactive transport simulator to predict the effects of CO 2 injection and complex reactions with the reservoir rock. Significant impediments were encountered during the execution phase of the project. The only known commercial reactive transport simulator was incapable of simulating complex geochemistry modeled in this project. Significant effort and project funding was expended to determine the limitations of both the commercial simulator and the Lawrence Berkeley National Laboratory (LBNL) R&D simulator, TOUGHREACT available to the project. A simplified layer cake model approximating the volume of the RMOTC targeted reservoirs was defined with 1-3 minerals eventually modeled with limited success. Modeling reactive transport in porous media requires significant computational power. In this project, up to 24 processors were used to model a limited mineral set of 1-3 minerals. In addition, geomechanical aspects of injecting CO 2 into closed, semi-open, and open systems in various well completion methods was simulated. Enhanced Oil Recovery (EOR) as a storage method was not modeled. A robust and stable simulation dataset or base case was developed and used to create a master dataset with embedded instructions for input to the ED/RSM software. Little success was achieved toward the objective of the project using the commercial simulator or the LBNL simulator versions available during the time of this project. Several hundred realizations were run with the commercial simulator and ED/RSM software, most having convergence problems and terminating prematurely. A proxy model for full field CO 2 injection sequestration utilization and storage was not capable of being developed with software available for this project. Though the chemistry is reasonably known and understood, based on the amount of effort and huge computational time required, predicting CO 2 sequestration storage capacity in geologic formations to within the program goals of ±30% proved unsuccessful.« less

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

NASA Astrophysics Data System (ADS)

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

2012-12-01

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

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.

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

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

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

Cook, J.R.

1994-04-15

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

Geochemical signatures of possible deep-seated ore deposits in Tertiary volcanic centers, Arizona and New Mexico, U.S.A.

USGS Publications Warehouse

Watts, K.C.; Hassemer, J.R.

1989-01-01

A reconnaissance geochemical survey of stream drainages within 21,000 km2 of southeastern Arizona and southwestern New Mexico shows broad zones of low-level to moderate contrast anomalies, many associated with mid-Tertiary eruptive centers and Tertiary fault zones. Of these eruptive centers, few are known to contain metallic deposits, and most of those known are minor. This, however, may be more a function of shallow erosion level than an indication of the absence of mineralization, since hydrothermal alteration and Fe-Mn-oxide staining are widespread, and geochemical anomalies are pervasive over a larger part of the region than outcrop observations would predict. Accordingly, interpretations of the geochemical data use considerations of relative erosion levels, and inferred element zonalities, to focus on possible undiscovered deposits in the subsurface of base-, precious-, and rare-metal deposits of plutonic-volcanic association. In order to enhance the identification of specific deep targets, we use the empirically determined ratio: Ag+Mn+Pb+Zn+Ba Au+Mo+Cu+Bi+W This ratio is based on reported metal contents of nonmagnetic heavy-mineral samples from the drainage sediment, determined by emission spectrographic analysis. Before the ratio was computed for each sample site, the data were normalized to a previously estimated regional threshold value. A regional isopleth map was then prepared, using a cell-averaging computer routine, with contours drawn at the 25th, 50th, 75th, 80th, 90th, 95th and 99th percentiles of the computed data. ?? 1989.

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.

The Merensky Reef in the Chineisky Pluton (Siberia)? A myth or a reality?

NASA Astrophysics Data System (ADS)

Zhitova, L.; Sharapov, V.; Zhukova, I.

2006-12-01

It is a dream of each geologist to find a `Merensky Reef' in each layered basic intrusion. Scientists have been trying many various techniques to come this dream to reality. The most perspective way to do so is probably a combination of physicochemical and computer modeling of layered basic intrusion crystallization together with fluid and melt inclusions studies in situ. This combination allows us to do the following: 1. To study boundary conditions for separation of low density gas phase and salt melt from the crystallizing primary basic melt in large magma chamber. 2. To determine correct quantitative parameters for formation of residual fluid-bearing brines extracting high metal concentrations. 3. To compute critical levels for substance differentiation at phase, geochemical and other `barriers' in those basic mantle-crust ore magmatic systems. 4. To model metal extraction, transportation and deposition at these `barriers' for systems of various `silicate melt - residual salt brines' ratios under the conditions of continental lithosphere. Comparison of real and modeled data allows us to conclude if a formation of a narrow zone of high metal concentration is possible at those critical levels (phase and geochemical `barriers'). The above-mentioned algorithm has been used for the Chineisky Pluton (the Transbaikal region, Siberia). Fortunately we have found our own `Merensky Reef', which happened to be a PGE enrichment marginal zone of the Chineisky Pluton due to specific fluid regime of crystallization! This work was supported by the Ministry for Russian Science and Education, Grant #DSP.2.1.1.702.

Sensitivity analysis of 1-D dynamical model for basin analysis

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

Cao, S.

1987-01-01

Geological processes related to petroleum generation, migration and accumulation are very complicated in terms of time and variables involved, and it is very difficult to simulate these processes by laboratory experiments. For this reasons, many mathematic/computer models have been developed to simulate these geological processes based on geological, geophysical and geochemical principles. The sensitivity analysis in this study is a comprehensive examination on how geological, geophysical and geochemical parameters influence the reconstructions of geohistory, thermal history and hydrocarbon generation history using the 1-D fluid flow/compaction model developed in the Basin Modeling Group at the University of South Carolina. This studymore » shows the effects of some commonly used parameter such as depth, age, lithology, porosity, permeability, unconformity (eroded thickness and erosion time), temperature at sediment surface, bottom hole temperature, present day heat flow, thermal gradient, thermal conductivity and kerogen type and content on the evolutions of formation thickness, porosity, permeability, pressure with time and depth, heat flow with time, temperature with time and depth, vitrinite reflectance (Ro) and TTI with time and depth, and oil window in terms of time and depth, amount of hydrocarbons generated with time and depth. Lithology, present day heat flow and thermal conductivity are the most sensitive parameters in the reconstruction of temperature history.« less

Analytical techniques for identification and study of organic matter in returned lunar samples

NASA Technical Reports Server (NTRS)

Burlingame, A. L.

1974-01-01

The results of geochemical research are reviewed. Emphasis is placed on the contribution of mass spectrometric data to the solution of specific structural problems. Information on the mass spectrometric behavior of compounds of geochemical interest is reviewed and currently available techniques of particular importance to geochemistry, such as gas chromatograph-mass spectrometer coupling, modern sample introduction methods, and computer application in high resolution mass spectrometry, receive particular attention.

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.

COTHERM: Modelling fluid-rock interactions in Icelandic geothermal systems

NASA Astrophysics Data System (ADS)

Thien, Bruno; Kosakowski, Georg; Kulik, Dmitrii

2014-05-01

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

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

Numerical Simulation of Permeability Change in Wellbore Cement Fractures after Geomechanical Stress and Geochemical Reactions Using X-ray Computed Tomography Imaging.

PubMed

Kabilan, Senthil; Jung, Hun Bok; Kuprat, Andrew P; Beck, Anthon N; Varga, Tamas; Fernandez, Carlos A; Um, Wooyong

2016-06-21

X-ray microtomography (XMT) imaging combined with three-dimensional (3D) computational fluid dynamics (CFD) modeling technique was used to study the effect of geochemical and geomechanical processes on fracture permeability in composite Portland cement-basalt caprock core samples. The effect of fluid density and viscosity and two different pressure gradient conditions on fracture permeability was numerically studied by using fluids with varying density and viscosity and simulating two different pressure gradient conditions. After the application of geomechanical stress but before CO2-reaction, CFD revealed fluid flow increase, which resulted in increased fracture permeability. After CO2-reaction, XMT images displayed preferential precipitation of calcium carbonate within the fractures in the cement matrix and less precipitation in fractures located at the cement-basalt interface. CFD estimated changes in flow profile and differences in absolute values of flow velocity due to different pressure gradients. CFD was able to highlight the profound effect of fluid viscosity on velocity profile and fracture permeability. This study demonstrates the applicability of XMT imaging and CFD as powerful tools for characterizing the hydraulic properties of fractures in a number of applications like geologic carbon sequestration and storage, hydraulic fracturing for shale gas production, and enhanced geothermal systems.

Numerical Simulation of Permeability Change in Wellbore Cement Fractures after Geomechanical Stress and Geochemical Reactions Using X-ray Computed Tomography Imaging

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

Kabilan, Senthil; Jung, Hun Bok; Kuprat, Andrew P.

X-ray microtomography (XMT) imaging combined with a three-dimensional (3D) computational fluid dynamics (CFD) modeling technique was used to study the effect of geochemical and geomechanical processes on fracture properties in composite Portland cement–basalt caprock core samples. The effect of fluid properties and flow conditions on fracture permeability was numerically studied by using fluids with varying physical properties and simulating different pressure conditions. CFD revealed that the application of geomechanical stress led to increased fluid flow, which resulted in increased fracture permeability. After CO2-reaction, XMT images displayed preferential precipitation of calcium carbonate within the fractures in the cement matrix and lessmore » precipitation in fractures located at the cement–basalt interface. CFD predicted changes in flow characteristics and differences in absolute values of flow properties due to different pressure gradients. CFD was able to highlight the profound effect of fluid properties on flow characteristics and hydraulic properties of fractures. This study demonstrates the applicability of XMT imaging and CFD as powerful tools for characterizing the hydraulic properties of fractures in a number of applications like geologic carbon sequestration and storage, hydraulic fracturing for shale gas production, and enhanced geothermal systems.« less

Origin and thermal evolution of Mars

NASA Technical Reports Server (NTRS)

Schubert, G.; Solomon, Sean C.; Turcotte, D. L.; Drake, M. J.; Sleep, N. H.

1993-01-01

The thermal evolution of Mars is governed by subsolidus mantle convection beneath a thick lithosphere. Models of the interior evolution are developed by parameterizing mantle convective heat transport in terms of mantle viscosity, the superadiabatic temperature rise across the mantle and mantle heat production. Geological, geophysical, and geochemical observations of the composition and structure of the interior and of the timing of major events in Martian evolution, such as global differentiation, atmospheric outgassing and the formation of the hemispherical dichotomy and Tharsis, are used to constrain the model computations. Isotope systematics of SNC meteorites suggest core formation essentially contemporaneously with the completion of accretion. Other aspects of this investigation are discussed.

Parameters on reconstructions of geohistory, thermal history, and hydrocarbon generation history in a sedimentary basin

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

Cao, S.; Lerche, I.

1988-01-01

Geological processes related to petroleum generation, migration, and accumulation are very complicated in terms of time and variables involved, and are very difficult to simulate by laboratory experiments. For this reason, many mathematic/computer models have been developed to simulate these geological processes based on geological, geophysical, and geochemical principles. Unfortunately, none of these models can exactly simulate these processes because of the assumptions and simplifications made in these models and the errors in the input for the models. The sensitivity analysis is a comprehensive examination on how geological, geophysical, and geochemical parameters affect the reconstructions of geohistory, thermal history, andmore » hydrocarbon generation history. In this study, a one-dimensional fluid flow/compaction model has been used to run the sensitivity analysis. The authors will show the effects of some commonly used parameters such as depth, age, lithology, porosity, permeability, unconformity (time and eroded thickness), temperature at sediment surface, bottom hole temperature, present day heat flow, thermal gradient, thermal conductivity and kerogen type, and content on the evolutions of formation thickness, porosity, permeability, pressure with time and depth, heat flow with time, temperature with time and depth, vitrinite reflectance (R/sub 0/) and TTI with time and depth, oil window in terms of time and depth, and amount of hydrocarbon generated with time and depth.« less

GEOCHEMICAL MODELING OF ARSENIC SPECIATION AND MOBILIZATION: IMPLICATIONS FOR BIOREMEDIATION

EPA Science Inventory

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

Multifractal and Singularity Maps of soil surface moisture distribution derived from 2D image analysis.

NASA Astrophysics Data System (ADS)

Cumbrera, Ramiro; Millán, Humberto; Martín-Sotoca, Juan Jose; Pérez Soto, Luis; Sanchez, Maria Elena; Tarquis, Ana Maria

2016-04-01

Soil moisture distribution usually presents extreme variation at multiple spatial scales. Image analysis could be a useful tool for investigating these spatial patterns of apparent soil moisture at multiple resolutions. The objectives of the present work were (i) to describe the local scaling of apparent soil moisture distribution and (ii) to define apparent soil moisture patterns from vertical planes of Vertisol pit images. Two soil pits (0.70 m long × 0.60 m width × 0.30 m depth) were excavated on a bare Mazic Pellic Vertisol. One was excavated in April/2011 and the other pit was established in May/2011 after 3 days of a moderate rainfall event. Digital photographs were taken from each Vertisol pit using a Kodak™ digital camera. The mean image size was 1600 × 945 pixels with one physical pixel ≈373 μm of the photographed soil pit. For more details see Cumbrera et al. (2012). Geochemical exploration have found with increasingly interests and benefits of using fractal (power-law) models to characterize geochemical distribution, using the concentration-area (C-A) model (Cheng et al., 1994; Cheng, 2012). This method is 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. We have applied it to each soil image. The results show that, in spite of some computational and practical limitations, image analysis of apparent soil moisture patterns could be used to study the dynamical change of soil moisture sampling in agreement with previous results (Millán et al., 2016). 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. Cumbrera, R., Ana M. Tarquis, Gabriel Gascó, Humberto Millán (2012) Fractal scaling of apparent soil moisture estimated from vertical planes of Vertisol pit images. Journal of Hydrology (452-453), 205-212. Martin Sotoca; J.J. Antonio Saa-Requejo, Juan Grau and Ana M. Tarquis (2016). Segmentation of singularity maps in the context of soil porosity. Geophysical Research Abstracts, 18, EGU2016-11402. Millán, H., Cumbrera, R. and Ana M. Tarquis (2016) Multifractal and Levy-stable statistics of soil surface moisture distribution derived from 2D image analysis. Applied Mathematical Modelling, 40(3), 2384-2395.

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.

East-China Geochemistry Database (ECGD):A New Networking Database for North China Craton

NASA Astrophysics Data System (ADS)

Wang, X.; Ma, W.

2010-12-01

North China Craton is one of the best natural laboratories that research some Earth Dynamic questions[1]. Scientists made much progress in research on this area, and got vast geochemistry data, which are essential for answering many fundamental questions about the age, composition, structure, and evolution of the East China area. But the geochemical data have long been accessible only through the scientific literature and theses where they have been widely dispersed, making it difficult for the broad Geosciences community to find, access and efficiently use the full range of available data[2]. How to effectively store, manage, share and reuse the existing geochemical data in the North China Craton area? East-China Geochemistry Database(ECGD) is a networking geochemical scientific database system that has been designed based on WebGIS and relational database for the structured storage and retrieval of geochemical data and geological map information. It is integrated the functions of data retrieval, spatial visualization and online analysis. ECGD focus on three areas: 1.Storage and retrieval of geochemical data and geological map information. Research on the characters of geochemical data, including its composing and connecting of each other, we designed a relational database, which based on geochemical relational data model, to store a variety of geological sample information such as sampling locality, age, sample characteristics, reference, major elements, rare earth elements, trace elements and isotope system et al. And a web-based user-friendly interface is provided for constructing queries. 2.Data view. ECGD is committed to online data visualization by different ways, especially to view data in digital map with dynamic way. Because ECGD was integrated WebGIS technology, the query results can be mapped on digital map, which can be zoomed, translation and dot selection. Besides of view and output query results data by html, txt or xls formats, researchers also can generate classification thematic maps using query results, according different parameters. 3.Data analysis on-line. Here we designed lots of geochemical online analysis tools, including geochemical diagrams, CIPW computing, and so on, which allows researchers to analyze query data without download query results. Operation of all these analysis tools is very easy; users just do it by click mouse one or two time. In summary, ECGD provide a geochemical platform for researchers, whom to know where various data are, to view various data in a synthetic and dynamic way, and analyze interested data online. REFERENCES [1] S. Gao, R.L. Rudnick, and W.L. Xu, “Recycling deep cratonic lithosphere and generation of intraplate magmatism in the North China Craton,” Earth and Planetary Science Letters,270,41-53,2008. [2] K.A. Lehnert, U. Harms, and E. Ito, “Promises, Achievements, and Challenges of Networking Global Geoinformatics Resources - Experiences of GeosciNET and EarthChem,” Geophysical Research Abstracts, Vol.10, EGU2008-A-05242,2008.

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.

Effect of CO2 Solubility on Dissolution Rates of Minerals in Porous Media Imbibed with Brine: Actual Efficiency of CO2 Sequestration

NASA Astrophysics Data System (ADS)

Alizadeh Nomeli, M.; Riaz, A.

2016-12-01

A new model is developed for geochemical reactions to access dissolution rate of minerals in saline aquifers with respect to saturated concentration of dissolved CO2 as a function of parameters that are dynamically available during computer program execution such as pressure, temperature, and salinity. A general Arrhenius-type equation, with an explicit dependence on the pH of brine, is employed to determine the rates of mineral dissolution. The amount of dissolved CO2 is determined with the help of an accurate PVTx model for the temperature range of 50-100C and pressures up to 600 bar relevant to the geologic sequestration of CO2. We show how activity coefficients for a given salinity condition alters solubility, pH, and reaction rates. We further evaluate the significance of the pre-exponential factor and the reaction order associated with the modified Arrhenius equation to determine the sensitivity of the reaction rates as a function to the pH of the system. It is found that the model can reasonably reproduce experimental data with new parameters that we obtain from sensitivity studies. Using the new rate equation, we investigate geochemically induced alterations of fracture geometry due to mineral dissolution. Finally, we use our model to evaluate the effects of temperature, pressure, and salinity on the actual efficiency of CO2 storage.

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

An Integrated 3D Hydrogeological, Geophysical, and Microbiological Investigation of Geochemical Gradients in a Pristine Aquifer Located in Laurentian Hills, ON, Canada

NASA Astrophysics Data System (ADS)

Shirokova, V.; Graves, L.; Stojanovic, S.; Enright, A. M.; Bank, C.; Ferris, F. G.

2013-12-01

A pristine glaciofluvial aquifer displaying naturally occurring geochemical gradients was investigated using hydrogeological, geophysical, and microbiological methods. A network of 25 piezometers was used to collect samples for groundwater chemical analysis, including parameters such as total iron (Fe), ferrous iron (Fe2+), sulphate (SO42-), sulfur (S2-), ammonium (NH4+), nitrate (NO3-), nitrite (NO2-), silica (SiO2), phosphate (PO43-), pH, and oxidation reduction potential (ORP). Ion concentration values between piezometers were interpolated using kriging and inverse distance weighting. Yearly analysis of the network shows spatially and temporally persistent plumes of iron and sulfur. A 3D model of the aquifer was compiled to aid in the understanding of the nature and origin of the geochemical gradients. The resulting maps showed zones with high concentrations of dissolved total iron (predominantly soluble ferric iron and complexed iron compounds), followed immediately downgradient by a high concentration of ferrous iron. Similarly, zones of high sulfide concentration were followed by areas of high sulfate concentration. There was some overlap between the iron and sulfur plumes, and ion concentrations were higher in years with a lower water table elevation. Metagenomic analysis revealed a diverse microbial community in the sediment, capable of the biogeochemical cycling of iron, sulfur, and nitrogen. The aquifer basin, as bounded by a till aquitard, was delineated using ground penetrating radar tomography from 45 lines. The plumes corresponded to an area where there is large, channel-like depression in the till boundary. Flow vectors from hydrogeological modelling indicated increased velocity followed by a slowing and convergence of groundwater in this location. Resistivity values from 20 lines varied in general from high values (2000-6000 Ohm.m) above 1-2 m to lower values (less than 1000 Ohm.m) below 2 to a 5m depth. The resistivity surveys consistently showed low resistivity values in areas of ionic enrichment, the location of the geochemical plumes, and high resistivity values at the top of the vadose zone including below dry sand outcrops. Fluorescent microscopy suggests the plumes are associated with attached subsurface bacteria dominated by species such as Gallionella and Leptothrix. These bacteria are likely responsible for conductive anomalies (<200 Ohm.m), observed in the resistivity models, that were at the centre of areas with high ionic concentrations. The above aquifer chemical network is currently being computationally simulated, and attempts are being made to determine the extents to which biotic and abiotic processes contribute to the formation of the geochemical gradients.

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.

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.

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.

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

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.

Coastal California's Fog Aerobiology and Ecology: Designing and Testing an Optimal Passive Impactor Collection Unit

NASA Astrophysics Data System (ADS)

Gentry, D.; Whinnery, J. T.; Ly, V. T.; Travers, S. V.; Sagaga, J.; Dahlgren, R. P.

2017-12-01

Microorganisms play a major role in our biosphere due to their ability to alter water, carbon and other geochemical cycles. Fog and low-level cloud water can play a major role in dispersing and supporting such microbial diversity. An ideal region to gather these microorganisms for characterization is the central coast of California, where dense fog is common. Fog captured from an unmanned aerial vehicle (UAV) at different altitudes will be analyzed to better understand the nature of microorganisms in the lower atmosphere and their potential geochemical impacts. The capture design consists of a square-meter hydrophobic mesh that hangs from a carbon fiber rod attached to a UAV. The DJI M600, a hexacopter, will be utilized as the transport for the payload, the passive impactor collection unit (PICU). The M600 will hover in a fog bank at altitudes between 10 and 100 m collecting water samples via the PICU. A computational flow dynamics (CFD) model will optimize the PICU's size, shape and placement for maximum capture efficiency and to avoid contamination from the UAV downwash. On board, there will also be an altitude, temperature and barometric pressure sensor whose output is logged to an SD card. A scale model of the PICU has been tested with several different types of hydrophobic meshes in a fog chamber at 90-95% humidity; polypropylene was found to capture the fog droplets most efficiently at a rate of .0042 g/cm2/hour. If the amount collected is proportional to the area of mesh, the estimated amount of water collected under optimal fog and flight conditions by the impactor is 21.3 g. If successful, this work will help identify the organisms living in the lower atmosphere as well as their potential geochemical impacts.

The U.S. geological survey rass-statpac system for management and statistical reduction of geochemical data

USGS Publications Warehouse

VanTrump, G.; Miesch, A.T.

1977-01-01

RASS is an acronym for Rock Analysis Storage System and STATPAC, for Statistical Package. The RASS and STATPAC computer programs are integrated into the RASS-STATPAC system for the management and statistical reduction of geochemical data. The system, in its present form, has been in use for more than 9 yr by scores of U.S. Geological Survey geologists, geochemists, and other scientists engaged in a broad range of geologic and geochemical investigations. The principal advantage of the system is the flexibility afforded the user both in data searches and retrievals and in the manner of statistical treatment of data. The statistical programs provide for most types of statistical reduction normally used in geochemistry and petrology, but also contain bridges to other program systems for statistical processing and automatic plotting. ?? 1977.

Geochemical maps of the Cornplanter Roadless Area, Warren County, Pennsylvania

USGS Publications Warehouse

Lesure, Frank G.; Day, Gordon W.

1984-01-01

The U.S. Geological Survey (USGS) made a reconnaissance geochemical survey of the Cornplanter Roadless Area (fig. 1) to test for indistinct or unexposed mineral deposits that might be recognized by their geochemical halos or patterns formed by the distribution of trace elements. Lesure, assisted by Andrew E. Grosz, collected 22 stream-sediment, 63 soil, and 23 rock samples from within and dear the study area during October 1980. All samples were analyzed for 31 elements using semi-quantitative spectrographic methods by Day in USGS laboratories, Denver, Colo. (table 1). In addition, the samples were also analyzed for zinc by means of an atomic absorption method by B.F. Arbogast and W.C. Martin, USGS laboratories, Denver Colo. J.T. Hanley and P.G. Schruben formatted the analytical data by computer methods for table 1. 

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

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.

A kinetics database and scripts for PHREEQC

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

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

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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

NASA Astrophysics Data System (ADS)

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

2016-07-01

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

Geoneutrinos and Heat Production in the Earth: Constraints and Implications

ScienceCinema

McDonough, Bill

2017-12-29

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

Effects of Injected CO2 on Geomechanical Properties Due to Mineralogical Changes

NASA Astrophysics Data System (ADS)

Nguyen, B. N.; Hou, Z.; Bacon, D. H.; Murray, C. J.; White, J. A.

2013-12-01

Long-term injection and storage of CO2 in deep underground reservoirs may significantly modify the geomechanical behavior of rocks since CO2 can react with the constituent phases of reservoir rocks and modify their composition. This can lead to modifications of their geomechanical properties (i.e., elastic moduli, Biot's coefficients, and permeability). Modifications of rock geomechanical properties have important consequences as these directly control stress and strain distributions, affect conditions for fracture initiation and development and/or fault healing. This paper attempts to elucidate the geochemical effects of CO2 on geomechanical properties of typical reservoir rocks by means of numerical analyses using the STOMP-ABAQUS sequentially coupled simulator that includes the capability to handle geomechanics and the reactive transport of CO2 together with a module (EMTA) to compute the homogenized rock poroelastic properties as a function of composition changes. EMTA, a software module developed at PNNL, implements the standard and advanced Eshelby-Mori-Tanaka approaches to compute the thermoelastic properties of composite materials. In this work, EMTA will be implemented in the coupled STOMP-ABAQUS simulator as a user subroutine of ABAQUS and used to compute local elastic stiffness based on rock composition. Under the STOMP-ABAQUS approach, STOMP models are built to simulate aqueous and CO2 multiphase fluid flows, and relevant chemical reactions of pore fluids with minerals in the reservoirs. The ABAQUS models then read STOMP output data for cell center coordinates, gas pressures, aqueous pressures, temperatures, saturations, constituent volume fractions, as well as permeability and porosity that are affected by chemical reactions. These data are imported into ABAQUS meshes using a mapping procedure developed for the exchange of data between STOMP and ABAQUS. Constitutive models implemented in ABAQUS via user subroutines then compute stiffness, stresses, strains, pore pressure, permeability, porosity, and capillary pressure, and return updated permeability, porosity, and capillary pressure to STOMP at selected times. In preliminary work, the enhanced STOMP-ABAQUS sequentially coupled approach is validated and illustrated in an example analysis of a cylindrical rock specimen subjected to axial loading, confining pressure, and CO2 fluid injection. The geomechanical analysis accounting for CO2 reactions with rock constituents is compared to that without chemical reactions to elucidate the geochemical effects of injected CO2 on the response of the reservoir rock to stress.

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.

Chemical trends in ocean islands explained by plume–slab interaction

NASA Astrophysics Data System (ADS)

Dannberg, Juliane; Gassmöller, Rene

2018-04-01

Earth's surface shows many features, of which the genesis can be understood only through their connection with processes in Earth's deep interior. Recent studies indicate that spatial geochemical patterns at oceanic islands correspond to structures in the lowermost mantle inferred from seismic tomographic models. This suggests that hot, buoyant upwellings can carry chemical heterogeneities from the deep lower mantle toward the surface, providing a window to the composition of the lowermost mantle. The exact nature of this link between surface and deep Earth remains debated and poorly understood. Using computational models, we show that subducted slabs interacting with dense thermochemical piles can trigger the ascent of hot plumes that inherit chemical gradients present in the lowermost mantle. We identify two key factors controlling this process: (i) If slabs induce strong lower-mantle flow toward the edges of these piles where plumes rise, the pile-facing side of the plume preferentially samples material originating from the pile, and bilaterally asymmetric chemical zoning develops. (ii) The composition of the melt produced reflects this bilateral zoning if the overlying plate moves roughly perpendicular to the chemical gradient in the plume conduit. Our results explain some of the observed geochemical trends of oceanic islands and provide insights into how these trends may originate.

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

NASA Astrophysics Data System (ADS)

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

2017-04-01

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

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

USGS Publications Warehouse

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

1978-01-01

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

Modeling microbial communities: current, developing, and future technologies for predicting microbial community interaction.

PubMed

Larsen, Peter; Hamada, Yuki; Gilbert, Jack

2012-07-31

Never has there been a greater opportunity for investigating microbial communities. Not only are the profound effects of microbial ecology on every aspect of Earth's geochemical cycles beginning to be understood, but also the analytical and computational tools for investigating microbial Earth are undergoing a rapid revolution. This environmental microbial interactome, the system of interactions between the microbiome and the environment, has shaped the planet's past and will undoubtedly continue to do so in the future. We review recent approaches for modeling microbial community structures and the interactions of microbial populations with their environments. Different modeling approaches consider the environmental microbial interactome from different aspects, and each provides insights to different facets of microbial ecology. We discuss the challenges and opportunities for the future of microbial modeling and describe recent advances in microbial community modeling that are extending current descriptive technologies into a predictive science. Copyright © 2012 Elsevier B.V. All rights reserved.

Reduced order models for prediction of groundwater quality impacts from CO₂ and brine leakage

DOE PAGES

Zheng, Liange; Carroll, Susan; Bianchi, Marco; ...

2014-12-31

A careful assessment of the risk associated with geologic CO₂ storage is critical to the deployment of large-scale storage projects. A potential risk is the deterioration of groundwater quality caused by the leakage of CO₂ and brine leakage from deep subsurface reservoirs. In probabilistic risk assessment studies, numerical modeling is the primary tool employed to assess risk. However, the application of traditional numerical models to fully evaluate the impact of CO₂ leakage on groundwater can be computationally complex, demanding large processing times and resources, and involving large uncertainties. As an alternative, reduced order models (ROMs) can be used as highlymore » efficient surrogates for the complex process-based numerical models. In this study, we represent the complex hydrogeological and geochemical conditions in a heterogeneous aquifer and subsequent risk by developing and using two separate ROMs. The first ROM is derived from a model that accounts for the heterogeneous flow and transport conditions in the presence of complex leakage functions for CO₂ and brine. The second ROM is obtained from models that feature similar, but simplified flow and transport conditions, and allow for a more complex representation of all relevant geochemical reactions. To quantify possible impacts to groundwater aquifers, the basic risk metric is taken as the aquifer volume in which the water quality of the aquifer may be affected by an underlying CO₂ storage project. The integration of the two ROMs provides an estimate of the impacted aquifer volume taking into account uncertainties in flow, transport and chemical conditions. These two ROMs can be linked in a comprehensive system level model for quantitative risk assessment of the deep storage reservoir, wellbore leakage, and shallow aquifer impacts to assess the collective risk of CO₂ storage projects.« less

How the Geothermal Community Upped the Game for Computer Codes

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

None

The Geothermal Technologies Office Code Comparison Study brought 11 research institutions together to collaborate on coupled thermal, hydrologic, geomechanical, and geochemical numerical simulators. These codes have the potential to help facilitate widespread geothermal energy development.

Automatic measurements and computations for radiochemical analyses

USGS Publications Warehouse

Rosholt, J.N.; Dooley, J.R.

1960-01-01

In natural radioactive sources the most important radioactive daughter products useful for geochemical studies are protactinium-231, the alpha-emitting thorium isotopes, and the radium isotopes. To resolve the abundances of these thorium and radium isotopes by their characteristic decay and growth patterns, a large number of repeated alpha activity measurements on the two chemically separated elements were made over extended periods of time. Alpha scintillation counting with automatic measurements and sample changing is used to obtain the basic count data. Generation of the required theoretical decay and growth functions, varying with time, and the least squares solution of the overdetermined simultaneous count rate equations are done with a digital computer. Examples of the complex count rate equations which may be solved and results of a natural sample containing four ??-emitting isotopes of thorium are illustrated. These methods facilitate the determination of the radioactive sources on the large scale required for many geochemical investigations.

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.

Inverse modeling of geochemical and mechanical compaction in sedimentary basins

NASA Astrophysics Data System (ADS)

Colombo, Ivo; Porta, Giovanni Michele; Guadagnini, Alberto

2015-04-01

We study key phenomena driving the feedback between sediment compaction processes and fluid flow in stratified sedimentary basins formed through lithification of sand and clay sediments after deposition. Processes we consider are mechanic compaction of the host rock and the geochemical compaction due to quartz cementation in sandstones. Key objectives of our study include (i) the quantification of the influence of the uncertainty of the model input parameters on the model output and (ii) the application of an inverse modeling technique to field scale data. Proper accounting of the feedback between sediment compaction processes and fluid flow in the subsurface is key to quantify a wide set of environmentally and industrially relevant phenomena. These include, e.g., compaction-driven brine and/or saltwater flow at deep locations and its influence on (a) tracer concentrations observed in shallow sediments, (b) build up of fluid overpressure, (c) hydrocarbon generation and migration, (d) subsidence due to groundwater and/or hydrocarbons withdrawal, and (e) formation of ore deposits. Main processes driving the diagenesis of sediments after deposition are mechanical compaction due to overburden and precipitation/dissolution associated with reactive transport. The natural evolution of sedimentary basins is characterized by geological time scales, thus preventing direct and exhaustive measurement of the system dynamical changes. The outputs of compaction models are plagued by uncertainty because of the incomplete knowledge of the models and parameters governing diagenesis. Development of robust methodologies for inverse modeling and parameter estimation under uncertainty is therefore crucial to the quantification of natural compaction phenomena. We employ a numerical methodology based on three building blocks: (i) space-time discretization of the compaction process; (ii) representation of target output variables through a Polynomial Chaos Expansion (PCE); and (iii) model inversion (parameter estimation) within a maximum likelihood framework. In this context, the PCE-based surrogate model enables one to (i) minimize the computational cost associated with the (forward and inverse) modeling procedures leading to uncertainty quantification and parameter estimation, and (ii) compute the full set of Sobol indices quantifying the contribution of each uncertain parameter to the variability of target state variables. Results are illustrated through the simulation of one-dimensional test cases. The analyses focuses on the calibration of model parameters through literature field cases. The quality of parameter estimates is then analyzed as a function of number, type and location of data.

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.

A Physically Based Coupled Chemical and Physical Weathering Model for Simulating Soilscape Evolution

NASA Astrophysics Data System (ADS)

Willgoose, G. R.; Welivitiya, D.; Hancock, G. R.

2015-12-01

A critical missing link in existing landscape evolution models is a dynamic soil evolution models where soils co-evolve with the landform. Work by the authors over the last decade has demonstrated a computationally manageable model for soil profile evolution (soilscape evolution) based on physical weathering. For chemical weathering it is clear that full geochemistry models such as CrunchFlow and PHREEQC are too computationally intensive to be couplable to existing soilscape and landscape evolution models. This paper presents a simplification of CrunchFlow chemistry and physics that makes the task feasible, and generalises it for hillslope geomorphology applications. Results from this simplified model will be compared with field data for soil pedogenesis. Other researchers have previously proposed a number of very simple weathering functions (e.g. exponential, humped, reverse exponential) as conceptual models of the in-profile weathering process. The paper will show that all of these functions are possible for specific combinations of in-soil environmental, geochemical and geologic conditions, and the presentation will outline the key variables controlling which of these conceptual models can be realistic models of in-profile processes and under what conditions. The presentation will finish by discussing the coupling of this model with a physical weathering model, and will show sample results from our SSSPAM soilscape evolution model to illustrate the implications of including chemical weathering in the soilscape evolution model.

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.

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

NASA Astrophysics Data System (ADS)

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

2003-12-01

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

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

Origin and thermal evolution of Mars

NASA Technical Reports Server (NTRS)

Schubert, Gerald; Soloman, S. C.; Turcotte, D. L.; Drake, M. J.; Sleep, N. H.

1990-01-01

The thermal evolution of Mars is governed by subsolidus mantle convection beneath a thick lithosphere. Models of the interior evolution are developed by parameterizing mantle convective heat transport in terms of mantle viscosity, the superadiabatic temperature rise across the mantle, and mantle heat production. Geological, geophysical, and geochemical observations of the compositon and structure of the interior and of the timing of major events in Martian evolution are used to constrain the model computations. Such evolutionary events include global differentiation, atmospheric outgassing, and the formation of the hemispherical dichotomy and Tharsis. Numerical calculations of fully three-dimensional, spherical convection in a shell the size of the Martian mantle are performed to explore plausible patterns of Martian mantel convection and to relate convective features, such as plumes, to surface features, such as Tharsis. The results from the model calculations are presented.

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.

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.

Geologic field notes and geochemical analyses of outcrop and drill core from Mesoproterozoic rocks and iron-oxide deposits and prospects of southeast Missouri

USGS Publications Warehouse

Day, Warren C.; Granitto, Matthew

2014-01-01

The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources/Missouri Geological Survey, undertook a study from 1988 to 1994 on the iron-oxide deposits and their host Mesoproterozoic igneous rocks in southeastern Missouri. The project resulted in an improvement of our understanding of the geologic setting, mode of formation, and the composition of many of the known deposits and prospects and the associated rocks of the St. Francois terrane in Missouri. The goal for this earlier work was to allow the comparison of Missouri iron-oxide deposits in context with other iron oxide-copper ± uranium (IOCG) types of mineral deposits observed globally. The raw geochemical analyses were released originally through the USGS National Geochemical Database (NGDB, http://mrdata.usgs.gov). The data presented herein offers all of the field notes, locations, rock descriptions, and geochemical analyses in a coherent package to facilitate new research efforts in IOCG deposit types. The data are provided in both Microsoft Excel (Version Office 2010) spreadsheet format (*.xlsx) and MS-DOS text formats (*.txt) for ease of use by numerous computer programs.

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.

Characteristic analysis-1981: Final program and a possible discovery

USGS Publications Warehouse

McCammon, R.B.; Botbol, J.M.; Sinding-Larsen, R.; Bowen, R.W.

1983-01-01

The latest ornewest version of thecharacteristicanalysis (NCHARAN)computer program offers the exploration geologist a wide variety of options for integrating regionalized multivariate data. The options include the selection of regional cells for characterizing deposit models, the selection of variables that constitute the models, and the choice of logical combinations of variables that best represent these models. Moreover, the program provides for the display of results which, in turn, makes possible review, reselection, and refinement of a model. Most important, the performance of the above-mentioned steps in an interactive computing mode can result in a timely and meaningful interpretation of the data available to the exploration geologist. The most recent application of characteristic analysis has resulted in the possible discovery of economic sulfide mineralization in the Grong area in central Norway. Exploration data for 27 geophysical, geological, and geochemical variables were used to construct a mineralized and a lithogeochemical model for an area that contained a known massive sulfide deposit. The models were applied to exploration data collected from the Gjersvik area in the Grong mining district and resulted in the identification of two localities of possible mineralization. Detailed field examination revealed the presence of a sulfide vein system and a partially inverted stratigraphic sequence indicating the possible presence of a massive sulfide deposit at depth. ?? 1983 Plenum Publishing Corporation.

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.

Origin of geochemical mantle components: Role of spreading ridges and thermal evolution of mantle

NASA Astrophysics Data System (ADS)

Kimura, Jun-Ichi; Gill, James B.; van Keken, Peter E.; Kawabata, Hiroshi; Skora, Susanne

2017-02-01

We explore the element redistribution at mid-ocean ridges (MOR) using a numerical model to evaluate the role of decompression melting of the mantle in Earth's geochemical cycle, with focus on the formation of the depleted mantle component. Our model uses a trace element mass balance based on an internally consistent thermodynamic-petrologic computation to explain the composition of MOR basalt (MORB) and residual peridotite. Model results for MORB-like basalts from 3.5 to 0 Ga indicate a high mantle potential temperature (Tp) of 1650-1500°C during 3.5-1.5 Ga before decreasing gradually to ˜1300°C today. The source mantle composition changed from primitive (PM) to depleted as Tp decreased, but this source mantle is variable with an early depleted reservoir (EDR) mantle periodically present. We examine a two-stage Sr-Nd-Hf-Pb isotopic evolution of mantle residues from melting of PM or EDR at MORs. At high-Tp (3.5-1.5 Ga), the MOR process formed extremely depleted DMM. This coincided with formation of the majority of the continental crust, the subcontinental lithospheric mantle, and the enriched mantle components formed at subduction zones and now found in OIB. During cooler mantle conditions (1.5-0 Ga), the MOR process formed most of the modern ocean basin DMM. Changes in the mode of mantle convection from vigorous deep mantle recharge before ˜1.5 Ga to less vigorous afterward is suggested to explain the thermochemical mantle evolution.

NetpathXL - An Excel Interface to the Program NETPATH

USGS Publications Warehouse

Parkhurst, David L.; Charlton, Scott R.

2008-01-01

NetpathXL is a revised version of NETPATH that runs under Windows? operating systems. NETPATH is a computer program that uses inverse geochemical modeling techniques to calculate net geochemical reactions that can account for changes in water composition between initial and final evolutionary waters in hydrologic systems. The inverse models also can account for the isotopic composition of waters and can be used to estimate radiocarbon ages of dissolved carbon in ground water. NETPATH relies on an auxiliary, database program, DB, to enter the chemical analyses and to perform speciation calculations that define total concentrations of elements, charge balance, and redox state of aqueous solutions that are then used in inverse modeling. Instead of DB, NetpathXL relies on Microsoft Excel? to enter the chemical analyses. The speciation calculation formerly included in DB is implemented within the program NetpathXL. A program DBXL can be used to translate files from the old DB format (.lon files) to NetpathXL spreadsheets, or to create new NetpathXL spreadsheets. Once users have a NetpathXL spreadsheet with the proper format, new spreadsheets can be generated by copying or saving NetpathXL spreadsheets. In addition, DBXL can convert NetpathXL spreadsheets to PHREEQC input files. New capabilities in PHREEQC (version 2.15) allow solution compositions to be written to a .lon file, and inverse models developed in PHREEQC to be written as NetpathXL .pat and model files. NetpathXL can open NetpathXL spreadsheets, NETPATH-format path files (.pat files), and NetpathXL-format path files (.pat files). Once the speciation calculations have been performed on a spreadsheet file or a .pat file has been opened, the NetpathXL calculation engine is identical to the original NETPATH. Development of models and viewing results in NetpathXL rely on keyboard entry as in NETPATH.

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

DOE Data Explorer

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

2018-06-13

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

Strategies to predict metal mobility in surficial mining environments

USGS Publications Warehouse

Smith, Kathleen S.

2007-01-01

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

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

NASA Astrophysics Data System (ADS)

Chen, Guoxiong; Cheng, Qiuming

2016-02-01

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

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.

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.

Yucca Mountain Area Saturated Zone Dissolved Organic Carbon Isotopic Data

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

Thomas, James; Decker, David; Patterson, Gary

2007-06-25

Groundwater samples in the Yucca Mountain area were collected for chemical and isotopic analyses and measurements of water temperature, pH, specific conductivity, and alkalinity were obtained at the well or spring at the time of sampling. For this project, groundwater samples were analyzed for major-ion chemistry, deuterium, oxygen-18, and carbon isotopes of dissolved inorganic carbon (DIC) and dissolved organic carbon (DOC). The U.S. Geological Survey (USGS) performed all the fieldwork on this project including measurement of water chemistry field parameters and sample collection. The major ions dissolved in the groundwater, deuterium, oxygen-18, and carbon isotopes of dissolved inorganic carbon (DIC)more » were analyzed by the USGS. All preparation and processing of samples for DOC carbon isotopic analyses and geochemical modeling were performed by the Desert Research Institute (DRI). Analysis of the DOC carbon dioxide gas produced at DRI to obtain carbon-13 and carbon-14 values was conducted at the University of Arizona Accelerator Facility (a NSHE Yucca Mountain project QA qualified contract facility). The major-ion chemistry, deuterium, oxygen-18, and carbon isotopes of DIC were used in geochemical modeling (NETPATH) to determine groundwater sources, flow paths, mixing, and ages. The carbon isotopes of DOC were used to calculate groundwater ages that are independent of DIC model corrected carbon-14 ages. The DIC model corrected carbon-14 calculated ages were used to evaluate groundwater travel times for mixtures of water including water beneath Yucca Mountain. When possible, groundwater travel times were calculated for groundwater flow from beneath Yucca Mountain to down gradient sample sites. DOC carbon-14 groundwater ages were also calculated for groundwaters in the Yucca Mountain area. When possible, groundwater travel times were estimated for groundwater flow from beneath Yucca Mountain to down gradient groundwater sample sites using the DOC calculated groundwater ages. The DIC calculated groundwater ages were compared with DOC calculated groundwater ages and both of these ages were compared to travel times developed in ground-water flow and transport models. If nuclear waste is stored in Yucca Mountain, the saturated zone is the final barrier against the release of radionuclides to the environment. The most recent rendition of the TSPA takes little credit for the presence of the saturated zone and is a testament to the inadequate understanding of this important barrier. If radionuclides reach the saturated zone beneath Yucca Mountain, then there is a travel time before they would leave the Yucca Mountain area and flow down gradient to the Amargosa Valley area. Knowing how long it takes groundwater in the saturated zone to flow from beneath Yucca Mountain to down gradient areas is critical information for potential radionuclide transport. Radionuclide transport in groundwater may be the quickest pathway for radionuclides in the proposed Yucca Mountain repository to reach land surface by way of groundwater pumped in Amargosa Valley. An alternative approach to ground-water flow and transport models to determine the travel time of radionuclides from beneath Yucca Mountain to down gradient areas in the saturated zone is by carbon-14 dating of both inorganic and organic carbon dissolved in the groundwater. A standard method of determining ground-water ages is to measure the carbon-13 and carbon-14 of DIC in the groundwater and then correct the measured carbon-14 along a flow path for geochemical reactions that involve carbon containing phases. These geochemical reactions are constrained by carbon-13 and isotopic fractionations. Without correcting for geochemical reactions, the ground-water ages calculated from only the differences in carbon-14 measured along a flow path (assuming the decrease in carbon-14 is due strictly to radioactive decay) could be tens of thousands of years too old. The computer program NETPATH, developed by the USGS, is the best geochemical program for correcting carbon-14 activities for geochemical reactions. The DIC carbon-14 corrected ages can be further constrained by measuring the carbon isotopes of DOC. Because the only source of organic carbon in aquifers is almost always greater than 40,000 years old, any organic carbon that may be added to the groundwater would contain no carbon-14. Thus, ground-water ages determined by carbon isotopes of DOC should be maximum ages that can be used to constrain DIC corrected ages.« less

Geochemical orientation for mineral exploration in the Hashemite Kingdom of Jordan

USGS Publications Warehouse

Overstreet, W.C.; Grimes, D.J.; Seitz, J.F.

1982-01-01

This report is a supplement to previous accounts of geochemical exploration conducted in the Hashemite Kingdom of Jordan by the Natural Resources Authority of the Royal Government of Jordan and the U.S. Geological Survey. The field work on which this report is based was sponsored by the U.S. Agency for International Development, U.S. Department of State. Procedures used in collecting various kinds of rocks, ores, slags, eluvial and alluvial sediments, heavy-mineral concentrates, and organic materials for use as geochemical sample media are summarized, as are the laboratory procedures followed for the analysis of these sample materials by semiquantitative spectrographic, atomic absorption, fluorometric, and X-ray diffraction methods. Geochemical evaluations of the possibilities for economic mineral deposits in certain areas are presented. The results of these preliminary investigations open concepts for further use in geochemical exploration in the search for metallic mineral deposits in Jordan. Perhaps the most desirable new activity would be hydrogeochemical exploration for uranium and base metals, accompanied by interpretation of such remote-sensing data as results of airborne radiometric surveys and computer-enhanced LANDSAT imagery. For more conventional approaches to geochemical exploration, however, several fundamental problems regarding proper choice of geochemical sample media for different geologic and geographic parts of the Country must be solved before effective surveys can be made. The present results also show that such common geochemical exploration techniques as the determination of the trace-element contents of soils, plant ash, and slags have direct application also toward the resolution of several archaeological problems in Jordan. These include the relation of trace-elements chemistry of local soils to the composition of botanic remains, the trace-elements composition of slags to the technological development of the extractive metallurgy of copper and iron in the region, and the use of charcoal from slags for the C-14 dating of periods of archaeometallurgical activity. Less directly, interpretations based on the distribution in time and space of the archaeometallurgical activities of the region might add to the knowledge of early climatic conditions and vegetative cover of the area.

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

Sewage contamination in the upper Mississippi River as measured by the fecal sterol, coprostanol

USGS Publications Warehouse

Writer, J.H.; Leenheer, J.A.; Barber, L.B.; Amy, G.L.; Chapra, S.C.

1995-01-01

The molecular sewage indicator, coprostanol, was measured in bed sediments of the Mississippi River for the purpose of determining sewage contamination. Coprostanol is a non-ionic, non-polar, organic molecule that associates with sediments in surface waters, and concentrations of coprostanol in bed sediments provide an indication of long-term sewage loads. Because coprostanol concentrations are dependent on particle size and percent organic carbon, a ratio between coprostanol (sewage sources) and cholestanol + cholesterol (sewage and non-sewage sources) was used to remove the biases related to particle size and percent organic carbon. The dynamics of contaminant transport in the Upper Mississippi River are influenced by both hydrologic and geochemical parameters. A mass balance model incorporating environmental parameters such as river and tributary discharge, suspended sediment concentration, fraction of organic carbon, sedimentation rates, municipal discharges and coprostanol decay rates was developed that describes coprostanol concentrations and therefore, expected patterns of municipal sewage effects on the Upper Mississippi River. Comparison of the computed and the measured coprostanol concentrations provides insight into the complex hydrologic and geochemical processes of contaminant transport and the ability to link measured chemical concentrations with hydrologic characteristics of the Mississippi River.

New Numerical Approaches for Modeling Thermochemical Convection in a Compositionally Stratified Fluid

NASA Astrophysics Data System (ADS)

Puckett, E. G.; Turcotte, D. L.; He, Y.; Lokavarapu, H. V.; Robey, J.; Kellogg, L. H.

2017-12-01

Geochemical observations of mantle-derived rocks favor a nearly homogeneous upper mantle, the source of mid-ocean ridge basalts (MORB), and heterogeneous lower mantle regions.Plumes that generate ocean island basalts are thought to sample the lower mantle regions and exhibit more heterogeneity than MORB.These regions have been associated with lower mantle structures known as large low shear velocity provinces below Africa and the South Pacific.The isolation of these regions is attributed to compositional differences and density stratification that, consequently, have been the subject of computational and laboratory modeling designed to determine the parameter regime in which layering is stable and understanding how layering evolves.Mathematical models of persistent compositional interfaces in the Earth's mantle may be inherently unstable, at least in some regions of the parameter space relevant to the mantle.Computing approximations to solutions of such problems presents severe challenges, even to state-of-the-art numerical methods.Some numerical algorithms for modeling the interface between distinct compositions smear the interface at the boundary between compositions, such as methods that add numerical diffusion or `artificial viscosity' in order to stabilize the algorithm. We present two new algorithms for maintaining high-resolution and sharp computational boundaries in computations of these types of problems: a discontinuous Galerkin method with a bound preserving limiter and a Volume-of-Fluid interface tracking algorithm.We compare these new methods with two approaches widely used for modeling the advection of two distinct thermally driven compositional fields in mantle convection computations: a high-order accurate finite element advection algorithm with entropy viscosity and a particle method.We compare the performance of these four algorithms on three problems, including computing an approximation to the solution of an initially compositionally stratified fluid at Ra = 105 with buoyancy numbers {B} that vary from no stratification at B = 0 to stratified flow at large B.

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.

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

NASA Astrophysics Data System (ADS)

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

2017-04-01

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

Computer finds ore

NASA Astrophysics Data System (ADS)

Bell, Peter M.

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

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

Novel proxies for reconstructing paleohydrology from ombrotrophic peatlands: biomarker and compound-specific H and C stable isotope ratios

NASA Astrophysics Data System (ADS)

Wang, J.; Nichols, J. E.; Huang, Y.

2008-12-01

Ombrotrophic peatlands are excellent archives for paleohydrologic information because they are hydrologically isolated from their surroundings. However, quantitative proxies for deciphering peatland archives are lacking. Here, we present development and application of novel organic geochemical methods for quantitative reconstruction of paleohydrology from the ombrotrophic sediments, and comparison of organic geochemical data with conventional paleoecological proxies. Application of these methods to the sediments of several North American and European peatlands has revealed significant changes in the hydroclimate throughout the Holocene. The plant assemblage living at the surface of the peatland is tightly controlled by surface moisture. Under wet conditions, Sphagnum mosses, with no active mechanism for drawing water from below the surface of the peatland, are dominant. During dry conditions, vascular plants are more productive relative to Sphagnum. A ratio of the abundance of two biomarkers representing Sphagnum and vascular plants sensitively records changes in hydrologic balance (Nichols et al., 2006, Org. Geochem. 37, 1505-1513). We have further developed stable isotope models to compute climate parameters from compound-specific H and C isotope ratios of biomarkers to create a more comprehensive climate reconstruction. Vascular plant leaf waxes carry the D/H ratio signature of precipitation that is little affected by evaporation, whereas the Sphagnum biomarker records isotopic ratios of the water at the peatland surface, which is strongly enriched by evaporation. Evaporation amount can be calculated using the differences between D/H ratios of the two types of biomarkers. C isotope ratios of Sphagnum biomarkers can also be used to quantify surface wetness. Methanotrophic bacteria live symbiotically with Sphagnum, providing isotopically light carbon for photosynthesis. These bacteria are more active when the Sphagnum is wet, thus providing more 13C-depleted CO2. Using a mass balance model we can use the carbon isotope ratios of Sphagnum biomarkers to assess the contribution of methane-derived CO2, and hence, the wetness of the peatland surface.

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

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

NASA Astrophysics Data System (ADS)

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

2017-06-01

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

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

USGS Publications Warehouse

Charlton, Scott R.; Parkhurst, David L.

2011-01-01

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

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

USGS Publications Warehouse

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

2011-01-01

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

Updated Reference Model for Heat Generation in the Lithosphere

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

Long-term climate change and the geochemical cycle of carbon

NASA Technical Reports Server (NTRS)

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

1988-01-01

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

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

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

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

2005-11-15

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

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

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

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.

Geochemical Exploration Techniques Applicable in the Search for Copper Deposits

USGS Publications Warehouse

Chaffee, Maurice A.

1975-01-01

Geochemical exploration is an important part of copper-resource evaluation. A large number of geochemical exploration techniques, both proved and untried, are available to the geochemist to use in the search for new copper deposits. Analyses of whole-rock samples have been used in both regional and local geochemical exploration surveys in the search for copper. Analyses of mineral separates, such as biotite, magnetite, and sulfides, have also been used. Analyses of soil samples are widely used in geochemical exploration, especially for localized surveys. It is important to distinguish between residual and transported soil types. Orientation studies should always be conducted prior to a geochemical investigation in a given area in order to determine the best soil horizon and the best size of soil material for sampling in that area. Silty frost boils, caliche, and desert varnish are specialized types of soil samples that might be useful sampling media. Soil gas is a new and potentially valuable geochemical sampling medium, especially in exploring for buried mineral deposits in arid regions. Gaseous products in samples of soil may be related to base-metal deposits and include mercury vapor, sulfur dioxide, hydrogen sulfide, carbon oxysulfide, carbon dioxide, hydrogen, oxygen, nitrogen, the noble gases, the halogens, and many hydrocarbon compounds. Transported materials that have been used in geochemical sampling programs include glacial float boulders, glacial till, esker gravels, stream sediments, stream-sediment concentrates, and lake sediments. Stream-sediment sampling is probably the most widely used and most successful geochemical exploration technique. Hydrogeochemical exploration programs have utilized hot- and cold-spring waters and their precipitates as well as waters from lakes, streams, and wells. Organic gel found in lakes and at stream mouths is an unproved sampling medium. Suspended material and dissolved gases in any type of water may also be useful media. Samples of ice and snow have been used for limited geochemical surveys. Both geobotanical and biogeochemical surveys have been successful in locating copper deposits in many parts of the world. Micro-organisms, including bacteria and algae, are other unproved media that should be studied. Animals can be used in geochemical-prospecting programs. Dogs have been used quite successfully to sniff out hidden and exposed sulfide minerals. Tennite mounds are commonly composed of subsurface material, but have not as yet proved to be useful in locating buried mineral deposits. Animal tissue and waste products are essentially unproved but potentially valuable sampling media. Knowledge of the location of areas where trace-element-associated diseases in animals and man are endemic as well as a better understanding of these diseases, may aid in identifying regions that are enriched in or depleted of various elements, including copper. Results of analyses of gases in the atmosphere are proving valuable in mineral-exploration surveys. Studies involving metallic compounds exhaled by plants into the atmosphere, and of particulate matter suspended in the atmosphere are reviewed these methods may become important in the future. Remote-sensing techniques are useful for making indirect measurements of geochemical responses. Two techniques applicable to geochemical exploration are neutron-activation analysis and gamma-ray spectrometry. Aerial photography is especially useful in vegetation surveys. Radar imagery is an unproved but potentially valuable method for use in studies of vegetation in perpetually clouded regions. With the advent of modern computers, many new techniques, such as correlation analysis, regression analysis, discriminant analysis, factor analysis, cluster analysis, trend-surface analysis, and moving-average analysis can be applied to geochemical data sets. Selective use of these techniques can provide new insights into the interpretatio

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.

New numerical approaches for modeling thermochemical convection in a compositionally stratified fluid

NASA Astrophysics Data System (ADS)

Puckett, Elbridge Gerry; Turcotte, Donald L.; He, Ying; Lokavarapu, Harsha; Robey, Jonathan M.; Kellogg, Louise H.

2018-03-01

Geochemical observations of mantle-derived rocks favor a nearly homogeneous upper mantle, the source of mid-ocean ridge basalts (MORB), and heterogeneous lower mantle regions. Plumes that generate ocean island basalts are thought to sample the lower mantle regions and exhibit more heterogeneity than MORB. These regions have been associated with lower mantle structures known as large low shear velocity provinces (LLSVPS) below Africa and the South Pacific. The isolation of these regions is attributed to compositional differences and density stratification that, consequently, have been the subject of computational and laboratory modeling designed to determine the parameter regime in which layering is stable and understanding how layering evolves. Mathematical models of persistent compositional interfaces in the Earth's mantle may be inherently unstable, at least in some regions of the parameter space relevant to the mantle. Computing approximations to solutions of such problems presents severe challenges, even to state-of-the-art numerical methods. Some numerical algorithms for modeling the interface between distinct compositions smear the interface at the boundary between compositions, such as methods that add numerical diffusion or 'artificial viscosity' in order to stabilize the algorithm. We present two new algorithms for maintaining high-resolution and sharp computational boundaries in computations of these types of problems: a discontinuous Galerkin method with a bound preserving limiter and a Volume-of-Fluid interface tracking algorithm. We compare these new methods with two approaches widely used for modeling the advection of two distinct thermally driven compositional fields in mantle convection computations: a high-order accurate finite element advection algorithm with entropy viscosity and a particle method that carries a scalar quantity representing the location of each compositional field. All four algorithms are implemented in the open source finite element code ASPECT, which we use to compute the velocity, pressure, and temperature associated with the underlying flow field. We compare the performance of these four algorithms on three problems, including computing an approximation to the solution of an initially compositionally stratified fluid at Ra =105 with buoyancy numbers B that vary from no stratification at B = 0 to stratified flow at large B .

A reappraisal of instrumental magnetic measurements made in Western Europe before AD 1750: confronting historical geomagnetism and archeomagnetism

NASA Astrophysics Data System (ADS)

Le Goff, Maxime; Gallet, Yves

2017-02-01

We present a new compilation and analysis of historical geomagnetic measurements made in Western Europe before AD 1750. The dataset in its ensemble provides a coherent evolution of magnetic field directions. Several data points excluded from previous analyses actually appear very consistent with most of the present compilation. A new average historical curve is computed for Paris, which is in very good agreement with the archeomagnetic data obtained in France, while significantly differing from the directional curve expected for Paris before AD 1675 based on the gufm1 model (Jackson et al. in Philos Trans R Soc Lond A 358:957-990, 2000). This finding suggests that the older segment of the gufm1 model lacks reliability and should be improved. Similarly, the historical part of the regional geomagnetic field model built for Europe by Pavón-Carrasco et al. (Geochem Geophys Geosyst 10:Q03013, 2009) should be revised because it erroneously incorporates directions derived from the gufm1 model.

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.

Coupling R and PHREEQC: an interactive and extensible environment for efficient programming of geochemical models

NASA Astrophysics Data System (ADS)

De Lucia, Marco; Kühn, Michael

2013-04-01

PHREEQC [1] is a widely used non-interactive open source software for speciation, batch-reactions, one-dimensional transport, and inverse geochemical caclulations. It represents the tool of choice for many researchers and practicioners for a broad set of geochemical problems, underground CO2 storage among others. Its open source nature, the flexibility to program arbitrary kinetic laws for the chemical reactions, as well as a thorough implementation of the Pitzer formalism explain its success and longevity. However, its non-interactive architecture make it cumbersome to couple PHREEQC to transport programs to achieve reactive transport simulations [2], but also to overcome the limitations of PHREEQC itself regarding the setup of large numbers of simulations - for example exploring wide ranges of conditions - and the graphical evaluation of the results. This has been the main motivation leading to the development of an interface with the high level language and environment for statistical computing and graphics GNU R [3]. The interface consists of minor modifications in PHREEQC's C source code, only affecting data I/O, plus on the R side a bunch of helper functions used to setup the simulations - basically automated manipulation of PHREEQC's input files, which are text files - and to collect and visualize the results. The most relevant subset of PHREEQC's capabilities and features are fully usable through the interface. Illustratory examples for the utility of this programmable interface were given in the framework of the research project this developement originated from: CLEAN [4], a project investigating the feasibility of enhanced gas recovery combined with CO2 storage. This interface allowed us to successfully and easily manipulate, compare and refit chemical databases, perform sensitivity analysis by combinatory variations of parameters, and all that in an environment which is both scriptable and interactive, with all results directly available for further manipulations and visualization in a powerful high level language, and benefiting from an enormous amount of third-party open source R extensions. The possibility to rapidly prototype complex algorithms involving geochemical modelling is in our opinion a huge advantage. A demonstration is given by the successful evaluation of a strategy to reduce the CPU-time needed to perform reactive transport simulations in a sequential coupling scheme. The idea is the "reduction" of the number of actual chemical simulations to perform at every time step, by searching for "duplicates" of each chemical simulations in the grid: such comparison involves typically a huge number of elements (one chemical simulation for grid element for time step) and a quite large number of variables (concentrations and mineral abundances). However, through the straightforward implementation of the prototype algorithm through the R/PHREEQC interface, we found out that the scan is extremely cost-effective in terms of CPU-time and typically allows a relevant speedup for simulations starting from a homogeneous or zone-homogeneous state. This speedup can even greatily exceed that of parallelization in some favorable but not unfrequent case. This feature should therefore be implemented in reactive transport simulators. References [1] Parkhurst D, Appelo C (1999) Users guide to PHREEQC (version 2). Tech. rep, U.S. Geological Survey. [2] Beyer C, Li D, De Lucia M, Kühn M, Bauer S (2012): Modelling CO2-induced fluid-rock interactions in the Altensalzwedel gas reservoir. Part II: coupled reactive transport simulation. Environ. Earth Sci., 67, 2, 573-588. [3] R Core Team (2012) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/. [4] Kühn M, Münch U (2012) CLEAN: CO2 Large-Scale Enhanced Gas Recovery. GEOTECHNOLOGIEN Science Report No. 19. Series: Advanced. Technologies in Earth Sciences, 199 p, ISBN 978-3-642-31676-0.

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.

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

PubMed

Laabidi, Ezzeddine; Bouhlila, Rachida

2015-07-01

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

Reconstructing mantle volatile contents through the veil of degassing

NASA Astrophysics Data System (ADS)

Tucker, J.; Mukhopadhyay, S.; Gonnermann, H. M.

2014-12-01

The abundance of volatile elements in the mantle reveals critical information about the Earth's origin and evolution such as the chemical constituents that built the Earth and material exchange between the mantle and exosphere. However, due to magmatic degassing, volatile element abundances measured in basalts usually do not represent those in undegassed magmas and hence in the mantle source of the basalts. While estimates of average mantle concentrations of some volatile species can be obtained, such as from the 3He flux into the oceans, volatile element variability within the mantle remains poorly constrained. Here, we use CO2-He-Ne-Ar-Xe measurements in basalts and a new degassing model to reconstruct the initial volatile contents of 8 MORBs from the Mid-Atlantic Ridge and Southwest Indian Ridge that span a wide geochemical range from depleted to enriched MORBs. We first show that equilibrium degassing (e.g. Rayleigh degassing), cannot simultaneously fit the measured CO2-He-Ne-Ar-Xe compositions in MORBs and argue that kinetic fractionation between bubbles and melt lowers the dissolved ratios of light to heavy noble gas species in the melt from that expected at equilibrium. We present a degassing model (after Gonnermann and Mukhopadhyay, 2007) that explicitly accounts for diffusive fractionation between melt and bubbles. The model computes the degassed composition based on an initial volatile composition and a diffusive timescale. To reconstruct the undegassed volatile content of a sample, we find the initial composition and degassing timescale which minimize the misfit between predicted and measured degassed compositions. Initial 3He contents calculated for the 8 MORB samples vary by a factor of ~7. We observe a correlation between initial 3He and CO2 contents, indicating relatively constant CO2/3He ratios despite the geochemical diversity and variable gas content in the basalts. Importantly, the gas-rich popping rock from the North Atlantic, as well as the average mantle ratio computed from the ridge 3He flux and independently estimated CO2 content fall along the same correlation. This observation suggests that undegassed CO2 and noble gas concentrations can be reconstructed in individual samples through measurement of noble gases and CO2 in erupted basalts.

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.

Geochemical Atlas of the San Jose and Golfito quadrangles, Costa Rica. Atlas Geoquimico de los cuadrangulos de San Jose y Golfito, Costa Rica (in English and Spanish)

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

Not Available

The Geochemical Atlas of the San Jose and Golfito 1:200,000-scale quadrangles, Costa Rica, was produced to help stimulate the growth of the Costa Rican mining industry and, thus, to benefit the economy of the country. As a result of the geochemical data presented in the Atlas, future exploration for metallic minerals in Costa Rica can be focused on specific areas that have the highest potential for mineralization. Stream-sediment samples were collected from drainage basins within the two quadrangles. These samples were analyzed for 50 elements and the results were displayed as computer-generated color maps. Each map shows the variation inmore » abundance of a single element within the quadrangle. Basic statistics, geological and cultural data are included as insets in each map to assist in interpretation. In the Golfito quadrangle, the geochemical data do not clearly indicate undiscovered gold mineralization. The areas known to contain placer (alluvial) gold are heavily affected by mining activity. Statistical treatment of the geochemical data is necessary before it will be possible to determine the gold potential of this quadrangle. In San Jose quadrangle, gold and the pathfinder elements, arsenic and antimony, are indicators of the gold mineralization characteristic of the Costa Rican gold district located in the Tilaran-Montes del Aguacate Range. This work shows that high concentrations of these elements occur in samples collected downstream from active gold mines. More importantly, the high concentrations of gold, arsenic, and antimony in sediment samples from an area southeast of the known gold district suggest a previously unknown extension of the district. This postulated extension underlain by Tertiary volcanic rocks which host the gold deposits within the gold district. The geochemical data, displayed herein, also indicate that drainage basins north of Ciudad Quesada on the flanks of Volcan Platanar have high gold potential.« less

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.

Regionally and seasonally differentiated primary production in the North Atlantic

NASA Astrophysics Data System (ADS)

Sathyendranath, Shubha; Longhurst, Alan; Caverhill, Carla M.; Platt, Trevor

1995-10-01

A bio-geochemical classification of the N. Atlantic Basin is presented according to which the basin is first divided into four primary algal domains: Polar, West-Wind, Trades and Coastal. These are in turn sub-divided into smaller provinces. The classification is based on differences in the physical environment which are likely to influence regional algal dynamics. The seasonally-differentiated parameters of the photosynthesis-light curve ( P-I curve) and parameters that define the vertical structure in chlorophyll profile are then established for each province, based on an analysis of an archive of over 6000 chlorophyll profiles, and over 1800 P-I curves. These are then combined with satellite-derived chlorophyll data for the N. Atlantic, and information on cloud cover, to compute primary production at the annual scale. using a model that computes spectral transmission of light underwater, and spectral, photosynthetic response of phytoplankton to available light. The results are compared with earlier, satellite-derived, estimates of basin-scale primary production.

Scientific Communication of Geochemical Data and the Use of Computer Databases.

ERIC Educational Resources Information Center

Le Bas, M. J.; Durham, J.

1989-01-01

Describes a scheme in the United Kingdom that coordinates geochemistry publications with a computerized geochemistry database. The database comprises not only data published in the journals but also the remainder of the pertinent data set. The discussion covers the database design; collection, storage and retrieval of data; and plans for future…

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.

Application of the Geo-Anomaly Unit Concept in Quantitative Delineation and Assessment of Gold Ore Targets in Western Shandong Uplift Terrain, Eastern China

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

Chen Yongqing, E-mail: ydonglai@mail.cgs.gov.cn; Zhao Pengda; Chen Jianguo

2001-03-15

A number of large and giant ore deposits have been discovered within the relatively small areas of lithospheric structure anomalies, including various boundary zones of tectonic plates. The regions have become the well-known intercontinental ore-forming belts, such as the circum-Pacific gold-copper, copper-molybdenum, and tungsten-tin metallogenic belts. These belts are typical geological anomalous areas. An investigation into the hydrothermal ore deposits in different regions in the former Soviet Union illustrated that the geologic structures of ore fields of almost all major commercial deposits have distinct features compared with the neighboring areas. These areas with distinct features are defined as geo-anomalies. Amore » geo-anomaly refers to such a geologic body or a combination of bodies that their composition, texture-structure, and genesis are significantly different from those of their surroundings. A geo-anomaly unit (GU) is an area containing distinct features that can be delineated with integrated ore-forming information using computer techniques on the basis of the geo-anomaly concept. Herein, the GU concept is illustrated by a case study of delineating the gold ore targets in the western Shandong uplift terrain, eastern China. It includes: (1) analyses of gold ore-forming factors; (2) compilation of normalized regional geochemical map and extraction of geochemical anomalies; (3) compilation of gravitational and aeromagnetic tectonic skeleton map and extraction of gravitational and aeromagnetic anomalies; (4) extraction of circular and linear anomalies from remote-sensing Landsat TM images; (5) establishment of a geo-anomaly conceptual model associated with known gold mineralization; (6) establishment of gold ore-forming favorability by computing techniques; and (7) delineation and assessment of ore-forming units. The units with high favorability are suggested as ore targets.« less

A geochemical atlas of North Carolina, USA

USGS Publications Warehouse

Reid, J.C.

1993-01-01

A geochemical atlas of North Carolina, U.S.A., was prepared using National Uranium Resource Evaluation (NURE) stream-sediment data. Before termination of the NURE program, sampling of nearly the entire state (48,666 square miles of land area) was completed and geochemical analyses were obtained. The NURE data are applicable to mineral exploration, agriculture, waste disposal siting issues, health, and environmental studies. Applications in state government include resource surveys to assist mineral exploration by identifying geochemical anomalies and areas of mineralization. Agriculture seeks to identify areas with favorable (or unfavorable) conditions for plant growth, disease, and crop productivity. Trace elements such as cobalt, copper, chromium, iron, manganese, zinc, and molybdenum must be present within narrow ranges in soils for optimum growth and productivity. Trace elements as a contributing factor to disease are of concern to health professionals. Industry can use pH and conductivity data for water samples to site facilities which require specific water quality. The North Carolina NURE database consists of stream-sediment samples, groundwater samples, and stream-water analyses. The statewide database consists of 6,744 stream-sediment sites, 5,778 groundwater sample sites, and 295 stream-water sites. Neutron activation analyses were provided for U, Br, Cl, F, Mn, Na, Al, V, Dy in groundwater and stream water, and for U, Th, Hf, Ce, Fe, Mn, Na, Sc, Ti, V, Al, Dy, Eu, La, Sm, Yb, and Lu in stream sediments. Supplemental analyses by other techniques were reported on U (extractable), Ag, As, Ba, Be, Ca, Co, Cr, Cu, K, Li, Mg, Mo, Nb, Ni, P, Pb, Se, Sn, Sr, W, Y, and Zn for 4,619 stream-sediment samples. A small subset of 334 stream samples was analyzed for gold. The goal of the atlas was to make available the statewide NURE data with minimal interpretation to enable prospective users to modify and manipulate the data for their end use. The atlas provides only very general indication of geochemical distribution patterns and should not be used for site specific studies. The atlas maps for each element were computer-generated at the state's geographic information system (Center for Geographic Information and Analysis [CGIA]). The Division of Statistics and Information Services provided input files. The maps in the atlas are point maps. Each sample is represented by a symbol generally corresponding to a quartile class. Other reports will transmit sample and analytical data for state regions. Data are tentatively planned to be available on disks in spreadsheet format for personal computers. During the second phase of this project, stream-sediment samples are being assigned to state geologic map unit names using a GIS system to determine background and anomaly values. Subsequent publications will make this geochemical data and accompanying interpretations available to a wide spectrum of interdisciplinary users. ?? 1993.

Computational methods for reactive transport modeling: A Gibbs energy minimization approach for multiphase equilibrium calculations

NASA Astrophysics Data System (ADS)

Leal, Allan M. M.; Kulik, Dmitrii A.; Kosakowski, Georg

2016-02-01

We present a numerical method for multiphase chemical equilibrium calculations based on a Gibbs energy minimization approach. The method can accurately and efficiently determine the stable phase assemblage at equilibrium independently of the type of phases and species that constitute the chemical system. We have successfully applied our chemical equilibrium algorithm in reactive transport simulations to demonstrate its effective use in computationally intensive applications. We used FEniCS to solve the governing partial differential equations of mass transport in porous media using finite element methods in unstructured meshes. Our equilibrium calculations were benchmarked with GEMS3K, the numerical kernel of the geochemical package GEMS. This allowed us to compare our results with a well-established Gibbs energy minimization algorithm, as well as their performance on every mesh node, at every time step of the transport simulation. The benchmark shows that our novel chemical equilibrium algorithm is accurate, robust, and efficient for reactive transport applications, and it is an improvement over the Gibbs energy minimization algorithm used in GEMS3K. The proposed chemical equilibrium method has been implemented in Reaktoro, a unified framework for modeling chemically reactive systems, which is now used as an alternative numerical kernel of GEMS.

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.

PEaCH4 v.2.0: A modelling platform to predict early diagenetic processes in marine sediments with a focus on biogenic methane - Case study: Offshore Namibia

NASA Astrophysics Data System (ADS)

Arning, Esther T.; Häußler, Steffen; van Berk, Wolfgang; Schulz, Hans-Martin

2016-07-01

The modelling of early diagenetic processes in marine sediments is of interest in marine science, and in the oil and gas industry, here, especially with respect to methane occurrence and gas hydrate formation as resources. Early diagenesis in marine sediments evolves from a complex web of intertwining (bio)geochemical reactions. It comprises microbially catalysed reactions and inorganic mineral-water-gas interactions. A model that will describe and consider all of these reactions has to be complex. However, it should be user-friendly, as well as to be applicable for a broad community and not only for experts in the field of marine chemistry. The presented modelling platform PeaCH4 v.2.0 combines both aspects, and is Microsoft Excel©-based. The modelling tool is PHREEQC (version 2), a computer programme for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. The conceptual PEaCH4 model is based on the conversion of sediment-bound degradable organic matter. PEaCH4 v.2.0 was developed to quantify and predict early diagenetic processes in marine sediments with the focus on biogenic methane formation and its phase behaviour, and allows carbon mass balancing. In regard to the irreversible degradation of organic matter, it comprises a "reaction model" and a "kinetic model" to predict methane formation. Both approaches differ in their calculations and outputs as the "kinetic model" considers the modelling time to integrate temperature dependent biogenic methane formation in its calculations, whereas the "reaction model" simply relies on default organic matter degradation. With regard to the inorganic mineral-water-gas interactions, which are triggered by irreversible degradation of organic matter, PEaCH4 v.2.0 is based on chemical equilibrium thermodynamics, appropriate mass-action laws, and their temperature dependent equilibrium constants. The programme is exemplarily presented with the example of upwelling sediments off Namibia, ODP Leg 175, Site 1082. The application demonstrates that the modelling platform PEaCH4 v.2.0 provides a user-friendly, but complex scientific tool that delivers retraceable information about early diagenetic processes and products in marine sediments.

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

NASA Astrophysics Data System (ADS)

Asavin, A. M.

2001-12-01

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

User's guide to PHREEQC, a computer program for speciation, reaction-path, advective-transport, and inverse geochemical calculations

USGS Publications Warehouse

Parkhurst, D.L.

1995-01-01

PHREEQC is a computer program written in the C pwgranuning language that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations, (2) reaction-path and advective-transport calculations involving specified irreversible reactions, mixing of solutions, mineral and gas equilibria surface-complex-ation reactions, and ion-exchange reactions, and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for composition differences between waters, within specified compositional uncertainties. PHREEQC is derived from the Fortran program PHREEQE, but it has been completely rewritten in C with the addition many new capabilities. New features include the capabilities to use redox couples to distribute redox elements among their valence states in speciation calculations; to model ion-exchange and surface-compiexation reactions; to model reactions with a fixed-pressure, multicomponent gas phase (that is, a gas bubble); to calculate the mass of water in the aqueous phase during reaction and transport calculations; to keep track of the moles of minerals present in the solid phases and determine antomaticaHy the thermodynamically stable phase assemblage; to simulate advective transport in combination with PHREEQC's reaction-modeling capability; and to make inverse modeling calculations that allow for uncertainties in the analytical data. The user interface is improved through the use of a simplified approach to redox reactions, which includes explicit mole-balance equations for hydrogen and oxygen; the use of a revised input that is modular and completely free format; and the use of mineral names and standard chemical symbolism rather than index numbers. The use of (2 eliminates nearly all limitations on army sizes, including numbers of elements, aqueous species, solutions, phases, and lengths of character strings. A new equation solver that optimizes a set of equalities subject to both equality and inequality constraints is used to determine the thermodynamically stable set of phases in equilibrium with a solution. A more complete Newton-Raphson formulation, master-species switching, and scaling of the algebraic equations reduce the number of failures of the nunmrical method in PHREEQC relative to PHREEQE. This report presents the equations that are the basis for chemical equilibrium and inverse-modeling calculations in PHREEQC, describes the input for the program, and presents twelve examples that demonstrate most of the program's capabilities.

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

PubMed

Zheng, Liange; Samper, Javier; Montenegro, Luis

2011-09-25

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

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

NASA Technical Reports Server (NTRS)

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

1974-01-01

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

The hills are alive: Earth surface dynamics in the University of Arizona Landscape Evolution Observatory

NASA Astrophysics Data System (ADS)

DeLong, S.; Troch, P. A.; Barron-Gafford, G. A.; Huxman, T. E.; Pelletier, J. D.; Dontsova, K.; Niu, G.; Chorover, J.; Zeng, X.

2012-12-01

To meet the challenge of predicting landscape-scale changes in Earth system behavior, the University of Arizona has designed and constructed a new large-scale and community-oriented scientific facility - the Landscape Evolution Observatory (LEO). The primary scientific objectives are to quantify interactions among hydrologic partitioning, geochemical weathering, ecology, microbiology, atmospheric processes, and geomorphic change associated with incipient hillslope development. LEO consists of three identical, sloping, 333 m2 convergent landscapes inside a 5,000 m2 environmentally controlled facility. These engineered landscapes contain 1 meter of basaltic tephra ground to homogenous loamy sand and contains a spatially dense sensor and sampler network capable of resolving meter-scale lateral heterogeneity and sub-meter scale vertical heterogeneity in moisture, energy and carbon states and fluxes. Each ~1000 metric ton landscape has load cells embedded into the structure to measure changes in total system mass with 0.05% full-scale repeatability (equivalent to less than 1 cm of precipitation), to facilitate better quantification of evapotraspiration. Each landscape has an engineered rain system that allows application of precipitation at rates between3 and 45 mm/hr. These landscapes are being studied in replicate as "bare soil" for an initial period of several years. After this initial phase, heat- and drought-tolerant vascular plant communities will be introduced. Introduction of vascular plants is expected to change how water, carbon, and energy cycle through the landscapes, with potentially dramatic effects on co-evolution of the physical and biological systems. LEO also provides a physical comparison to computer models that are designed to predict interactions among hydrological, geochemical, atmospheric, ecological and geomorphic processes in changing climates. These computer models will be improved by comparing their predictions to physical measurements made in LEO. The main focus of our iterative modeling and measurement discovery cycle is to use rapid data assimilation to facilitate validation of newly coupled open-source Earth systems models. LEO will be a community resource for Earth system science research, education, and outreach. The LEO project operational philosophy includes 1) open and real-time availability of sensor network data, 2) a framework for community collaboration and facility access that includes integration of new or comparative measurement capabilities into existing facility cyberinfrastructure, 3) community-guided science planning and 4) development of novel education and outreach programs.Artistic rendering of the University of Arizona Landscape Evolution Observatory

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

NASA Astrophysics Data System (ADS)

Shamberger, Patrick J.; Garcia, Michael O.

2007-02-01

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

Predicting Biological Information Flow in a Model Oxygen Minimum Zone

NASA Astrophysics Data System (ADS)

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

2016-02-01

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

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

PubMed

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

2016-12-01

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

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

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

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

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

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

DOE PAGES

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

2016-09-16

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

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

NASA Astrophysics Data System (ADS)

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

2009-12-01

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

MetPetDB: A database for metamorphic geochemistry

NASA Astrophysics Data System (ADS)

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

2009-12-01

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

New insight into unstable hillslopes hydrology from hydrogeochemical modelling.

NASA Astrophysics Data System (ADS)

Bertrand, C.; Marc, V.; Malet, J.-P.

2010-05-01

In the black marl outcrops of the French South Alps, sub surface flow conditions are considered as the main triggering factor for initiation and reactivation of landslides. The problem is traditionally addressed in term of hydrological processes (how does percolation to the water table occur?) but in some cases the origin of water is also in question. Direct rainfall is generally assumed as the only water source for groundwater recharge in shallow hillslope aquifers. The bedrock is also supposed impervious and continuous. Yet the geological environment of the study area is very complex owing to the geological history of this alpine sector. The autochthonous callovo-oxfordian black marl bedrock is highly tectonized (Maquaire et al., 2003) and may be affected by large, possibly draining discontinuities. A deep water inflow at the slip surface may at least locally result in increase the pore pressure and decrease the effective shearing resistance of the landslide material. In the active slow-moving landslide of Super-Sauze (Malet and Maquaire, 2003), this question has been addressed using hydrochemical investigations. The groundwater was sampled during five field campaigns uniformly spread out over the year from a network of boreholes. Water chemistry data were completed by geochemical and mineralogical analyses of the marl material. The major hydro-geochemical processes over area proved (1) mixing processes, (2) pyrite alteration, (3) dissolution/precipitation of carbonates and (4) cations exchange (de Montety et al., 2007). A geochemical modelling was carried out using the model Phreeqc (Parkhurst and Appelo, version 2.15, 2008) to check how suitable was observed water chemistry with the reservoir characteristics. The modelling exercise was based on a kinetics approach of soil-water interactions. The model simulates the rock alteration by the dissolution of the primary minerals and the precipitation of new phases. Initial parameters were obtained from geochemical and mineralogical analyses or from the literature (kinetics constants). The simulations showed that pH, sulphate and calcium concentrations in groundwater could be reproduced from reasonable assumptions. However, the observed high concentrations in magnesium and sodium were not correctly simulated by the model. Furthermore, a particular anomaly in the Na+ concentration was observed in the most active part of the landslide. Lastly, isotopic investigation showed that groundwater 3H content in this sector was significantly lower than groundwater content in the other parts of the landslide and lower than the mean rainwater content. This result showed that the mean groundwater age in the active part was probably higher than elsewhere in the landslide. All these arguments led us to conclude that groundwater was locally recharged with saline waters from areas outside the watershed, coming up through the bedrock using major discontinuities. This assumption is in agreement with the geological context. de Montety, V., V. Marc, C. Emblanch, J.-P. Malet, C. Bertrand, O. Maquaire, and T. A. Bogaard, 2007, Identifying the origin of groundwater and flow processes in complex landslides affecting black marls: insights from a hydrochemical survey.: Earth Surface Processes and Landforms, v. 32, p. 32-48. Malet, J.-P. and Maquaire, O., 2003. Black marl earthflows mobility and long-term seasonal dynamic in southeastern France. In: Picarelli, L. (Ed). Proceedings of the International Conference on Fast Slope Movements: Prediction and Prevention for Risk Mitigation. Patron Editore, Bologna: 333-340. Maquaire, O., Malet, J.-P., Remaître, A., Locat, J., Klotz, S. and Guillon, J., 2003. Instability conditions of marly hillslopes: towards landsliding or gullying? The case of the Barcelonnette Bassin, South East France. Engineering Geology, 70(1-2): 109-130. Parkhurst, D.L. and Appelo, C.A.J., 1999, User's guide to PHREEQC (version 2)--A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations: U.S. Geological Survey Water-Resources Investigations Report 99-4259, 312 p.

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.

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

NASA Astrophysics Data System (ADS)

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

2017-06-01

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

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.

Mineralogic and petrologic implications of Viking geochemical results from Mars - Interim report

NASA Technical Reports Server (NTRS)

Baird, A. K.; Toulmin, P., III; Rose, H. J., Jr.; Christian, R. P.; Clark, B. C.; Keil, K.; Gooding, J. L.

1976-01-01

Chemical results from four samples of Martian fines delivered to Viking landers 1 and 2 are remarkably similar in that they all have high iron; moderate magnesium, calcium, and sulfur; low aluminum; and apparently very low alkalies and trace elements. This composition is best interpreted as representing the weathering products of mafic igneous rocks. A mineralogic model, derived from computer mixing studies and laboratory analog preparations, suggests that Mars fines could be an intimate mixture of about 80% iron-rich clay, about 10% magnesium sulfate (kieserite), about 5% carbonate (calcite), and about 5% iron oxides (hematite, magnetite, maghemite, goethite). The mafic nature of the present fines (distributed globally) and their probable source rocks seems to preclude large-scale planetary differentiation of a terrestrial nature.

Mineralogic and petrologic implications of viking geochemical results from Mars: interim report.

PubMed

Baird, A K; Toulmin, P; Clark, B C; Rose, H J; Keil, K; Christian, R P; Gooding, J L

1976-12-11

Chemical results from four samples of martian fines delivered to Viking landers 1 and 2 are remarkably similar in that they all have high iron; moderate magnesium, calcium, and sulfur; low aluminum; and apparently very low alkalies and trace elements. This composition is best interpreted as representing the weathering products of mafic igneous rocks. A mineralogic model, derived from computer mixing studies and laboratory analog preparations, suggests that Mars fines could be an intimate mixture of about 80 percent iron-rich clay, about 10 percent magnesium sulfate (kieserite?), about 5 percent carbonate (calcite), and about 5 percent iron oxides (hematite, magnetite, maghemite, goethite?). The mafic nature of the present fines (distributed globally) and their probable source rocks seems to preclude large-scale planetary differentiation of a terrestrial nature.

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

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

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

1997-05-01

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

Databases, data integration, and expert systems: new directions in mineral resource assessment and mineral exploration

USGS Publications Warehouse

McCammon, Richard B.; Ramani, Raja V.; Mozumdar, Bijoy K.; Samaddar, Arun B.

1994-01-01

Overcoming future difficulties in searching for ore deposits deeper in the earth's crust will require closer attention to the collection and analysis of more diverse types of data and to more efficient use of current computer technologies. Computer technologies of greatest interest include methods of storage and retrieval of resource information, methods for integrating geologic, geochemical, and geophysical data, and the introduction of advanced computer technologies such as expert systems, multivariate techniques, and neural networks. Much experience has been gained in the past few years in applying these technologies. More experience is needed if they are to be implemented for everyday use in future assessments and exploration.

Resident research associateships, postdoctoral research awards 1989: opportunities for research at the U.S. Geological Survey, U.S. Department of the Interior

USGS Publications Warehouse

,; ,

1989-01-01

The scientists of the U.S. Geological Survey are engaged in a wide range of geologic, geophysical, geochemical, hydrologic, and cartographic programs, including the application of computer science to them. These programs offer exciting possibilities for scientific achievement and professional growth to young scientists through participation as Research Associates.

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.

Geochemical processes governing the compositional features of the crater fumarolic field at Mt. Etna

NASA Astrophysics Data System (ADS)

Liotta, Marcello; Paonita, Antonio; Caracausi, Antonio; Martelli, Mauro; Rizzo, Andrea; Favara, Rocco

2010-05-01

Mt Etna is one of the most-active volcanoes in the world. It is characterized by major eruptions, frequent Strombolian activity, and ash emissions. The volcano summit consists of the central crater of Voragine surrounded by the three active cones of the North-East Crater, Bocca Nuova, and the South-East Crater. They are characterized by very fractured and unstable edges. Under these conditions most of the fractures represent preferential degassing pathways for volcanic fluids, so that the main fumarolic fields develop in such fractured areas. The geochemistry of the fumaroles at the summit area of Mt. Etna was investigated. Fumarolic samples were collected between June 2008 and August 2009. Gas samples were usually collected as 'dry gas' and analyzed for the concentrations of He, H2, O2, N2, CO, CH4, and CO2. Fumarolic gases were also sampled a few times using the classical Giggenbach bottles and Giggenbach-type bottles filled with ammonia and silver nitrate in order to determine the SO2/H2S ratio. In addition a novel method was employed in order to sample fumaroles characterized by high content of atmospheric gases. Two types of fumaroles were identified: low-temperature fumaroles, which are dominated by CO2 with minor amounts of SO2 and H2S, and negligible halogen contents, and high-temperature fumaroles, which are strongly air-contaminated and characterized by appreciable amounts of volcanogenic carbon, sulfur, and chlorine. Our data clearly indicate that secondary processes modify the composition of the fluids once they leave the magma body. A model based on thermodynamic data is proposed to explore such postmagmatic processes. We computed the equilibrium composition of magmatic gases that cool starting from magmatic temperatures under several pressure conditions. The model, which uses Etnean plume geochemistry as starting composition of fluids exsolved from magma, shows that SO2 and H2S control the redox conditions of the gas mixture during the cooling, until the reactions involving CO/CO2 and H2/H2O ratios are fully quenched. The scrubbing processes occurring subsequent to condensation and gas-liquid water interaction allow total removal of HCl and partial removal of sulfur species. During the ascent toward the surface, the concentration of CH4 increases in all fumaroles due to a modest contribution from hydrothermal fluid. A geochemical model for the interaction of pristine magmatic fluids with shallower systems is proposed. The model explains geochemical changes at the crater fumaroles in terms of variable hydrothermal and magmatic contributions, and modified thermodynamic conditions.

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.

Characterisation and modelling of mixing processes in groundwaters of a potential geological repository for nuclear wastes in crystalline rocks of Sweden.

PubMed

Gómez, Javier B; Gimeno, María J; Auqué, Luis F; Acero, Patricia

2014-01-15

This paper presents the mixing modelling results for the hydrogeochemical characterisation of groundwaters in the Laxemar area (Sweden). This area is one of the two sites that have been investigated, under the financial patronage of the Swedish Nuclear Waste and Management Co. (SKB), as possible candidates for hosting the proposed repository for the long-term storage of spent nuclear fuel. The classical geochemical modelling, interpreted in the light of the palaeohydrogeological history of the system, has shown that the driving process in the geochemical evolution of this groundwater system is the mixing between four end-member waters: a deep and old saline water, a glacial meltwater, an old marine water, and a meteoric water. In this paper we put the focus on mixing and its effects on the final chemical composition of the groundwaters using a comprehensive methodology that combines principal component analysis with mass balance calculations. This methodology allows us to test several combinations of end member waters and several combinations of compositional variables in order to find optimal solutions in terms of mixing proportions. We have applied this methodology to a dataset of 287 groundwater samples from the Laxemar area collected and analysed by SKB. The best model found uses four conservative elements (Cl, Br, oxygen-18 and deuterium), and computes mixing proportions with respect to three end member waters (saline, glacial and meteoric). Once the first order effect of mixing has been taken into account, water-rock interaction can be used to explain the remaining variability. In this way, the chemistry of each water sample can be obtained by using the mixing proportions for the conservative elements, only affected by mixing, or combining the mixing proportions and the chemical reactions for the non-conservative elements in the system, establishing the basis for predictive calculations. © 2013 Elsevier B.V. All rights reserved.

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.