iTOUGH2: A multiphysics simulation-optimization framework for analyzing subsurface systems

NASA Astrophysics Data System (ADS)

Finsterle, S.; Commer, M.; Edmiston, J. K.; Jung, Y.; Kowalsky, M. B.; Pau, G. S. H.; Wainwright, H. M.; Zhang, Y.

2017-11-01

iTOUGH2 is a simulation-optimization framework for the TOUGH suite of nonisothermal multiphase flow models and related simulators of geophysical, geochemical, and geomechanical processes. After appropriate parameterization of subsurface structures and their properties, iTOUGH2 runs simulations for multiple parameter sets and analyzes the resulting output for parameter estimation through automatic model calibration, local and global sensitivity analyses, data-worth analyses, and uncertainty propagation analyses. Development of iTOUGH2 is driven by scientific challenges and user needs, with new capabilities continually added to both the forward simulator and the optimization framework. This review article provides a summary description of methods and features implemented in iTOUGH2, and discusses the usefulness and limitations of an integrated simulation-optimization workflow in support of the characterization and analysis of complex multiphysics subsurface systems.

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

Continuum-based DFN-consistent numerical framework for the simulation of oxygen infiltration into fractured crystalline rocks

NASA Astrophysics Data System (ADS)

Trinchero, Paolo; Puigdomenech, Ignasi; Molinero, Jorge; Ebrahimi, Hedieh; Gylling, Björn; Svensson, Urban; Bosbach, Dirk; Deissmann, Guido

2017-05-01

We present an enhanced continuum-based approach for the modelling of groundwater flow coupled with reactive transport in crystalline fractured rocks. In the proposed formulation, flow, transport and geochemical parameters are represented onto a numerical grid using Discrete Fracture Network (DFN) derived parameters. The geochemical reactions are further constrained by field observations of mineral distribution. To illustrate how the approach can be used to include physical and geochemical complexities into reactive transport calculations, we have analysed the potential ingress of oxygenated glacial-meltwater in a heterogeneous fractured rock using the Forsmark site (Sweden) as an example. The results of high-performance reactive transport calculations show that, after a quick oxygen penetration, steady state conditions are attained where abiotic reactions (i.e. the dissolution of chlorite and the homogeneous oxidation of aqueous iron(II) ions) counterbalance advective oxygen fluxes. The results show that most of the chlorite becomes depleted in the highly conductive deformation zones where higher mineral surface areas are available for reactions.

Geobiology: A Conceptual Framework for Understanding Earth's Surface

NASA Astrophysics Data System (ADS)

Sumner, D. Y.

2016-12-01

A topic of study becomes a new field when it provides a useful conceptual framework for understanding suites of important processes. Geobiology integrates microbial biology with Earth sciences in a way that allows us to ask - and answer - deeper questions about Earth and the life on it. Recent studies of the oxidation of Earth's surface exemplify the impact of Geobiology as a new field. For decades, scientists have understood that Earth's surface was oxidized by photosynthesis. Geochemical records indicate dramatic redox changes both globally, e.g. the loss of MIF sulfur signatures due to formation of an ozone layer, and locally, as preserved in sedimentary rocks. However, these records depend critically on the dynamics of both the global biosphere and local microbial ecology. For example, an increase in global redox due to photosynthetic iron oxidation has different biogeochemical implications than an increase from oxygenic photosynthesis; O2 reacts very differently with organic matter and minerals than iron oxyhydroxides do, influencing microbial ecology as well as potential geochemical signatures in sedimentary rocks. Thus, studies of modern microbial communities provide insights into the interactions among metabolisms and geochemical gradients that have shaped Earth's redox history. For example, the ability of cyanobacteria to create O2 oases in benthic mats and soils on land provides a new framework for evaluating redox-sensitive elemental fluxes to the ocean. Similarly, genomic studies of Cyanobacteria have revealed close relatives, Melainabacteria, that are mostly obligate anaerobes. The evolutionary relationships between these two groups, as preserved in their genomes, reflect important microbial processes that led to oxidation of Earth's surface. By combining insights from microbial biology and sedimentary geochemistry, geobiologists will develop significantly more accurate models of the interactions between life and Earth.

Quantifying reactive transport processes governing arsenic mobility in a Bengal Delta aquifer

NASA Astrophysics Data System (ADS)

Rawson, Joey; Neidhardt, Harald; Siade, Adam; Berg, Michael; Prommer, Henning

2017-04-01

Over the last few decades significant progress has been made to characterize the extent and severity of groundwater arsenic pollution in S/SE Asia, and to understand the underlying geochemical processes. However, comparably little effort has been made to merge the findings from this research into quantitative frameworks that allow for a process-based quantitative analysis of observed arsenic behavior and predictions of its future fate. Therefore, this study developed and tested field-scale numerical modelling approaches to represent the primary and secondary geochemical processes associated with the reductive dissolution of Fe-oxy(hydr)oxides and the concomitant release of sorbed arsenic. We employed data from an in situ field experiment in the Bengal Delta Plain, which investigated the influence of labile organic matter (sucrose) on the mobility of Fe, Mn, and As. The data collected during the field experiment were used to guide our model development and to constrain the model parameterisation. Our results show that sucrose oxidation coupled to the reductive dissolution of Fe-oxy(hydr)oxides was accompanied by multiple secondary geochemical reactions that are not easily and uniquely identifiable and quantifiable. Those secondary reactions can explain the disparity between the observed Fe and As behavior. Our modelling results suggest that a significant fraction of the released As is scavenged through (co-)precipitation with newly formed Fe-minerals, specifically magnetite, rather than through sorption to pre-existing and freshly precipitated iron minerals.

Assessing the Role of Seafloor Weathering in Global Geochemical Cycling

NASA Astrophysics Data System (ADS)

Farahat, N. X.; Abbot, D. S.; Archer, D. E.

2015-12-01

Low-temperature alteration of the basaltic upper oceanic crust, known as seafloor weathering, has been proposed as a mechanism for long-term climate regulation similar to the continental climate-weathering negative feedback. Despite this potentially far-reaching impact of seafloor weathering on habitable planet evolution, existing modeling frameworks do not include the full scope of alteration reactions or recent findings of convective flow dynamics. We present a coupled fluid dynamic and geochemical numerical model of low-temperature, off-axis hydrothermal activity. This model is designed to explore the the seafloor weathering flux of carbon to the oceanic crust and its responsiveness to climate fluctuations. The model's ability to reproduce the seafloor weathering environment is evaluated by constructing numerical simulations for comparison with two low-temperature hydrothermal systems: A transect east of the Juan de Fuca Ridge and the southern Costa Rica Rift flank. We explore the sensitivity of carbon uptake by seafloor weathering on climate and geology by varying deep ocean temperature, seawater dissolved inorganic carbon, continental weathering inputs, and basaltic host rock in a suite of numerical experiments.

Mantle ingredients for making the fingerprint of Etna alkaline magmas: implications for shallow partial melting within the complex geodynamic framework of Eastern Sicily

NASA Astrophysics Data System (ADS)

Viccaro, Marco; Zuccarello, Francesco

2017-09-01

Mantle ingredients responsible for the signature of Etnean Na- and K-alkaline magmas and their relationships with short-term geochemical changes of the erupted volcanic rocks have been constrained through a partial melting model that considers major, trace elements and water contents in the produced liquids. Characteristics of the Etnean source for alkaline magmas have been supposed similar to those of the mantle accessible at a regional scale, namely below the Hyblean Plateau. The assumption that the Etnean mantle resembles the one beneath the Hyblean Plateau is justified by the large geochemical affinities of the Etnean hawaiites/K-trachybasalts and the Hyblean hawaiites/alkali basalts for what concerns both trace elements and isotope systematics. We have modeled partial melting of a composite source constituted by two rock types, inferred by lithological and geochemical features of the Hyblean xenoliths: 1) a spinel lherzolite bearing metasomatic, hydrous phases and 2) a garnet pyroxenite in form of veins intruded into the spinel lherzolite. The partial melting modeling has been applied to each rock type and the resulting primary liquids have been then mixed in various proportions. These compositions have been compared with some Etnean alkaline magmas of the post ∼60 ka activity, which were firstly re-equilibrated to mantle conditions through mass balance calculations. Our results put into evidence that concentrations of major and trace elements along with the water obtained from the modeling are remarkably comparable with those of Etnean melts re-equilibrated at primary conditions. Different proportions of the spinel lherzolite with variable modal contents of metasomatic phases and of the garnet pyroxenite can therefore account for the signature of a large spectrum of Etnean alkaline magmas and for their geochemical variability through time, emphasizing the crucial role played by compositional small-scale heterogeneity of the source. These heterogeneities are able to produce magmas with variable compositions and volatile contents, which can then undergo distinct histories of ascent and evolution, leading to the wide range of eruptive styles observed at Mt. Etna volcano. Being partial melting confined in the spinel facies of the mantle, our model implies that the source of Mt. Etna magmas might be rather shallow (<2 GPa; i.e., lesser than ca. 60 km), excluding the presence of deep, plume-like mantle structures responsible for magma generation. Partial melting should occur consequently as a response of mantle decompression within the framework of regional tectonics affecting the Eastern Sicily, which could be triggered by extensional tectonics and/or subduction-induced mantle upwelling.

Continuum-based DFN-consistent numerical framework for the simulation of oxygen infiltration into fractured crystalline rocks.

PubMed

Trinchero, Paolo; Puigdomenech, Ignasi; Molinero, Jorge; Ebrahimi, Hedieh; Gylling, Björn; Svensson, Urban; Bosbach, Dirk; Deissmann, Guido

2017-05-01

We present an enhanced continuum-based approach for the modelling of groundwater flow coupled with reactive transport in crystalline fractured rocks. In the proposed formulation, flow, transport and geochemical parameters are represented onto a numerical grid using Discrete Fracture Network (DFN) derived parameters. The geochemical reactions are further constrained by field observations of mineral distribution. To illustrate how the approach can be used to include physical and geochemical complexities into reactive transport calculations, we have analysed the potential ingress of oxygenated glacial-meltwater in a heterogeneous fractured rock using the Forsmark site (Sweden) as an example. The results of high-performance reactive transport calculations show that, after a quick oxygen penetration, steady state conditions are attained where abiotic reactions (i.e. the dissolution of chlorite and the homogeneous oxidation of aqueous iron(II) ions) counterbalance advective oxygen fluxes. The results show that most of the chlorite becomes depleted in the highly conductive deformation zones where higher mineral surface areas are available for reactions. Copyright © 2017 Elsevier B.V. All rights reserved.

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.

Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

Professional Paper 1708 is intended primarily for geoscientists in academia, industry, and government who are interested in Appalachian basin geology and its coal and petroleum resources. Other users, however, may find the topics, papers, and digital images valuable for land-use and policy planning. Among the anticipated benefits of the report are improvements in (1) resource assessment estimates and methodology, (2) exploration strategies, (3) basin models, and (4) energy use policies.

Fluid-rock geochemical interaction for modelling calibration in geothermal exploration in Indonesia

NASA Astrophysics Data System (ADS)

Deon, Fiorenza; Barnhoorn, Auke; Lievens, Caroline; Ryannugroho, Riskiray; Imaro, Tulus; Bruhn, David; van der Meer, Freek; Hutami, Rizki; Sibarani, Besteba; Sule, Rachmat; Saptadij, Nenny; Hecker, Christoph; Appelt, Oona; Wilke, Franziska

2017-04-01

Indonesia with its large, but partially unexplored geothermal potential is one of the most interesting and suitable places in the world to conduct geothermal exploration research. This study focuses on geothermal exploration based on fluid-rock geochemistry/geomechanics and aims to compile an overview on geochemical data-rock properties from important geothermal fields in Indonesia. The research carried out in the field and in the laboratory is performed in the framework of the GEOCAP cooperation (Geothermal Capacity Building program Indonesia- the Netherlands). The application of petrology and geochemistry accounts to a better understanding of areas where operating power plants exist but also helps in the initial exploration stage of green areas. Because of their relevance and geological setting geothermal fields in Java, Sulawesi and the sedimentary basin of central Sumatra have been chosen as focus areas of this study. Operators, universities and governmental agencies will benefit from this approach as it will be applied also to new green-field terrains. By comparing the characteristic of the fluids, the alteration petrology and the rock geochemistry we also aim to contribute to compile an overview of the geochemistry of the important geothermal fields in Indonesia. At the same time the rock petrology and fluid geochemistry will be used as input data to model the reservoir fluid composition along with T-P parameters with the geochemical workbench PHREEQC. The field and laboratory data are mandatory for both the implementation and validation of the model results.

Cryptotephras: the revolution in correlation and precision dating1

PubMed Central

DAVIES, SIWAN M

2015-01-01

From its Icelandic origins in the study of visible tephra horizons, tephrochronology took a remarkable step in the late 1980 s with the discovery of a ca. 4300-year-old microscopic ash layer in a Scottish peat bog. Since then, the search for these cryptotephra deposits in distal areas has gone from strength to strength. Indeed, a recent discovery demonstrates how a few fine-grained glass shards from an Alaskan eruption have been dispersed more than 7000 km to northern Europe. Instantaneous deposition of geochemically distinct volcanic ash over such large geographical areas gives rise to a powerful correlation tool with considerable potential for addressing a range of scientific questions. A prerequisite of this work is the establishment of regional tephrochronological frameworks that include well-constrained age estimates and robust geochemical signatures for each deposit. With distal sites revealing a complex record of previously unknown volcanic events, frameworks are regularly revised, and it has become apparent that some closely timed eruptions have similar geochemical signatures. The search for unique and robust geochemical fingerprints thus hinges on rigorous analysis by electron microprobe and laser ablation-inductively coupled plasma-mass spectrometry. Historical developments and significant breakthroughs are presented to chart the revolution in correlation and precision dating over the last 50 years using tephrochronology and cryptotephrochronology. PMID:27512240

Hydrologic Responses to CO2 Injection in Basalts Based on Flow-through Experiments

NASA Astrophysics Data System (ADS)

Thomas, D.; Hingerl, F.; Garing, C.; Bird, D. K.; Benson, S. M.; Maher, K.

2015-12-01

Experimental studies of basalt-CO2 interactions have increased our ability to predict geochemical responses within a mafic reservoir during geologic CO2 sequestration. However, the lack of flow-through experiments prevents the use of coupled hydrologic-geochemical models to predict evolution of permeability and porosity, critical parameters for assessing storage feasibility. We present here results of three flow-through experiments on an intact basalt core during which we employed X-ray Computed Tomography (CT) to quantify porosity evolution and fluid flow. Using a single core of glassy basaltic tuff from the Snake River Plain (Menan Buttes complex), we performed tracer tests using a solution of NaI (~100,000 ppm) before and after injection of CO2-saturated water at reservoir conditions (90 bar, 50°C) to image porosity and flow path distribution. During the tracer tests, CT scans were taken at 2.5-minute intervals, and outlet fluid was discretely sampled at the same intervals and subsequently measured via ICP-MS, enabling interpretation of the tracer breakthrough curve through both imaging and geochemical analyses. Comparison of the porosity distribution from before and after injection of CO2 shows an overall decrease in core-averaged porosity from 34% to 31.1%. Permeability decreased exponentially from ~4.9x10-12 m2 to 1.18 x10-12 m2. The decrease in porosity and permeability suggests geochemical transformations in the mineral assemblage of the core, which we observe through petrographic analysis of an unaltered sample of the same lithology in contrast with the altered core. There is a significant increase in grain coatings, as well as reduction in the grain size, suggesting dissolution re-precipitation mechanisms. Finally, to develop a framework for the coupled geochemical and hydrologic responses observed experimentally, we have calibrated a reactive transport model at the core scale using the TOUGHREACT simulator [1]. [1] Xu et al. (2011) Comput. Geosci.

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

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

NASA Astrophysics Data System (ADS)

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

2010-12-01

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

Modelling of the petroleum formation in the Mahakam sediments (Indonesia): Organic geochemical controls of the results

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

Brosse, E.; Burris, J.; Ouidin, J.L.

1990-06-01

Since the Miocene, the delta of the Mahakam River has accumulated thousands of meters of sediments in the eastern part of the Kutei Basin (Kalimantan, Indonesia). Source-rock candidates are the coals of the deltaic plain and several types of shales, mainly the delta front/prodelta area. Organic matter basically derives from higher plants, but each source facies presents important intrinsic variations of petroleum potential. These variations are overprinted by subsequent maturation trends. Geochemical and petrographical data are integrated on the general framework provided by a new synthetic interpretation of the sedimentary sequences, relying upon the concepts of seismic stratigraphy. From coremore » samples at a given level of maturation, the variations of several organic parameters are discussed in relation to the depositional paleoenvironment and to the possible precursors. 1D and 2D numerical routines are used to reconstruct the maturation history of source rocks. These tools are based upon a kinetic modeling of kerogen cracking. Model outputs are compared with observed maturation trends. The understanding of the initial organic facies distribution provides precise constraints in the selection of a homogenous samples set for this comparison purpose.« less

Is rhizosphere remediation sufficient for sustainable revegetation of mine tailings?

PubMed

Huang, Longbin; Baumgartl, Thomas; Mulligan, David

2012-07-01

Revegetation of mine tailings (fine-grained waste material) starts with the reconstruction of root zones, consisting of a rhizosphere horizon (mostly topsoil and/or amended tailings) and the support horizon beneath (i.e. equivalent to subsoil - mostly tailings), which must be physically and hydro-geochemically stable. This review aims to discuss key processes involved in the development of functional root zones within the context of direct revegetation of tailings and introduces a conceptual process of rehabilitating structure and function in the root zones based on a state transition model. Field studies on the revegetation of tailings (from processing base metal ore and bauxite residues) are reviewed. Particular focus is given to tailings' properties that limit remediation effectiveness. Aspects of root zone reconstruction and vegetation responses are also discussed. When reconstructing a root zone system, it is critical to restore physical structure and hydraulic functions across the whole root zone system. Only effective and holistically restored systems can control hydro-geochemical mobility of acutely and chronically toxic factors from the underlying horizon and maintain hydro-geochemical stability in the rhizosphere. Thereafter, soil biological capacity and ecological linkages (i.e. carbon and nutrient cycling) may be rehabilitated to integrate the root zones with revegetated plant communities into sustainable plant ecosystems. A conceptual framework of system transitions between the critical states of root zone development has been proposed. This will illustrate the rehabilitation process in root zone reconstruction and development for direct revegetation with sustainable plant communities. Sustainable phytostabilization of tailings requires the systematic consideration of hydro-geochemical interactions between the rhizosphere and the underlying supporting horizon. It further requires effective remediation strategies to develop hydro-geochemically stable and biologically functional root zones, which can facilitate the recovery of the microbial community and ecological linkages with revegetated plant communities.

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

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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.

In-Drift Microbial Communities

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

D. Jolley

2000-11-09

As directed by written work direction (CRWMS M and O 1999f), Performance Assessment (PA) developed a model for microbial communities in the engineered barrier system (EBS) as documented here. The purpose of this model is to assist Performance Assessment and its Engineered Barrier Performance Section in modeling the geochemical environment within a potential repository drift for TSPA-SR/LA, thus allowing PA to provide a more detailed and complete near-field geochemical model and to answer the key technical issues (KTI) raised in the NRC Issue Resolution Status Report (IRSR) for the Evolution of the Near Field Environment (NFE) Revision 2 (NRC 1999).more » This model and its predecessor (the in-drift microbial communities model as documented in Chapter 4 of the TSPA-VA Technical Basis Document, CRWMS M and O 1998a) was developed to respond to the applicable KTIs. Additionally, because of the previous development of the in-drift microbial communities model as documented in Chapter 4 of the TSPA-VA Technical Basis Document (CRWMS M and O 1998a), the M and O was effectively able to resolve a previous KTI concern regarding the effects of microbial processes on seepage and flow (NRC 1998). This document supercedes the in-drift microbial communities model as documented in Chapter 4 of the TSPA-VA Technical Basis Document (CRWMS M and O 1998a). This document provides the conceptual framework of the revised in-drift microbial communities model to be used in subsequent performance assessment (PA) analyses.« less

Geochemical heterogeneity in a small, stratigraphically complex moraine aquifer system (Ontario, Canada): Interpretation of flow and recharge using multiple geochemical parameters

USGS Publications Warehouse

Stotler, R.L.; Frape, S.K.; El Mugammar, H.T.; Johnston, C.; Judd-Henrey, I.; Harvey, F.E.; Drimmie, R.; Jones, J.P.

2011-01-01

The Waterloo Moraine is a stratigraphically complex system and is the major water supply to the cities of Kitchener and Waterloo in Ontario, Canada. Despite over 30 years of investigation, no attempt has been made to unify existing geochemical data into a single database. A composite view of the moraine geochemistry has been created using the available geochemical information, and a framework created for geochemical data synthesis of other similar flow systems. Regionally, fluid chemistry is highly heterogeneous, with large variations in both water type and total dissolved solids content. Locally, upper aquifer units are affected by nitrate and chloride from fertilizer and road salt. Typical upper-aquifer fluid chemistry is dominated by calcium, magnesium, and bicarbonate, a result of calcite and dolomite dissolution. Evidence also suggests that ion exchange and diffusion from tills and bedrock units accounts for some elevated sodium concentrations. Locally, hydraulic "windows" cross connect upper and lower aquifer units, which are typically separated by a clay till. Lower aquifer units are also affected by dedolomitization, mixing with bedrock water, and locally, upward diffusion of solutes from the bedrock aquifers. A map of areas where aquifer units are geochemically similar was constructed to highlight areas with potential hydraulic windows. ?? 2010 Springer-Verlag.

Introduction to TETHYS—an interdisciplinary GIS database for studying continental collisions

NASA Astrophysics Data System (ADS)

Khan, S. D.; Flower, M. F. J.; Sultan, M. I.; Sandvol, E.

2006-05-01

The TETHYS GIS database is being developed as a way to integrate relevant geologic, geophysical, geochemical, geochronologic, and remote sensing data bearing on Tethyan continental plate collisions. The project is predicated on a need for actualistic model 'templates' for interpreting the Earth's geologic record. Because of their time-transgressive character, Tethyan collisions offer 'actualistic' models for features such as continental 'escape', collision-induced upper mantle flow magmatism, and marginal basin opening, associated with modern convergent plate margins. Large integrated geochemical and geophysical databases allow for such models to be tested against the geologic record, leading to a better understanding of continental accretion throughout Earth history. The TETHYS database combines digital topographic and geologic information, remote sensing images, sample-based geochemical, geochronologic, and isotopic data (for pre- and post-collision igneous activity), and data for seismic tomography, shear-wave splitting, space geodesy, and information for plate tectonic reconstructions. Here, we report progress on developing such a database and the tools for manipulating and visualizing integrated 2-, 3-, and 4-d data sets with examples of research applications in progress. Based on an Oracle database system, linked with ArcIMS via ArcSDE, the TETHYS project is an evolving resource for researchers, educators, and others interested in studying the role of plate collisions in the process of continental accretion, and will be accessible as a node of the national Geosciences Cyberinfrastructure Network—GEON via the World-Wide Web and ultra-high speed internet2. Interim partial access to the data and metadata is available at: http://geoinfo.geosc.uh.edu/Tethys/ and http://www.esrs.wmich.edu/tethys.htm. We demonstrate the utility of the TETHYS database in building a framework for lithospheric interactions in continental collision and accretion.

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

NASA Astrophysics Data System (ADS)

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

2017-06-01

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

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.

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.

Modeling low-temperature geochemical processes: Chapter 2

USGS Publications Warehouse

Nordstrom, D. Kirk; Campbell, Kate M.

2014-01-01

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

Exploiting virtual sediment deposits to explore conceptual foundations

NASA Astrophysics Data System (ADS)

Dietze, Michael; Fuchs, Margret; Kreutzer, Sebastian

2017-04-01

Geomorphic concepts and hypotheses are usually formulated based on empiric data from the field or the laboratory (deduction). After translation into models they can be applied to case study scenarios (induction). However, the other way around - expressing hypotheses explicitly by models and test these by empiric data - is a rarely touched trail. There are several models tailored to investigate the boundary conditions and processes that generate, mobilise, route and eventually deposit sediment in a landscape. Thereby, the last part, sediment deposition, is usually omitted. Essentially, there is no model that explicitly focuses on mapping out the characteristics of sedimentary deposits - the material that is used by many disciplines to reconstruct landscape evolution. This contribution introduces the R-package sandbox, a model framework that allows creating and analysing virtual sediment sections for exploratory, explanatory, forecasting and inverse research questions. The R-package sandbox is a probabilistic and rule-based model framework for a wide range of possible applications. The model framework is used here to discuss a set of conceptual questions revolving around geochemical and geochronological methods, such as: How does sample size and sample volume affect age uncertainty? What determines the robustness of sediment fingerprinting results? How does the prepared grain size of the material of interest affect the analysis outcomes? Most of the concepts used in geosciences are underpinned by a set of assumptions, whose robustness and boundary conditions need to be assessed quantitatively. The R-package sandbox is a universal and flexible tool to engage with this challenge.

Is rhizosphere remediation sufficient for sustainable revegetation of mine tailings?

PubMed Central

Huang, Longbin; Baumgartl, Thomas; Mulligan, David

2012-01-01

Background Revegetation of mine tailings (fine-grained waste material) starts with the reconstruction of root zones, consisting of a rhizosphere horizon (mostly topsoil and/or amended tailings) and the support horizon beneath (i.e. equivalent to subsoil – mostly tailings), which must be physically and hydro-geochemically stable. This review aims to discuss key processes involved in the development of functional root zones within the context of direct revegetation of tailings and introduces a conceptual process of rehabilitating structure and function in the root zones based on a state transition model. Scope Field studies on the revegetation of tailings (from processing base metal ore and bauxite residues) are reviewed. Particular focus is given to tailings' properties that limit remediation effectiveness. Aspects of root zone reconstruction and vegetation responses are also discussed. Conclusions When reconstructing a root zone system, it is critical to restore physical structure and hydraulic functions across the whole root zone system. Only effective and holistically restored systems can control hydro-geochemical mobility of acutely and chronically toxic factors from the underlying horizon and maintain hydro-geochemical stability in the rhizosphere. Thereafter, soil biological capacity and ecological linkages (i.e. carbon and nutrient cycling) may be rehabilitated to integrate the root zones with revegetated plant communities into sustainable plant ecosystems. A conceptual framework of system transitions between the critical states of root zone development has been proposed. This will illustrate the rehabilitation process in root zone reconstruction and development for direct revegetation with sustainable plant communities. Sustainable phytostabilization of tailings requires the systematic consideration of hydro-geochemical interactions between the rhizosphere and the underlying supporting horizon. It further requires effective remediation strategies to develop hydro-geochemically stable and biologically functional root zones, which can facilitate the recovery of the microbial community and ecological linkages with revegetated plant communities. PMID:22648878

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.

Shale gas characterization based on geochemical and geophysical analysis: Case study of Brown shale, Pematang formation, Central Sumatra Basin

NASA Astrophysics Data System (ADS)

Haris, A.; Nastria, N.; Soebandrio, D.; Riyanto, A.

2017-07-01

Geochemical and geophysical analyses of shale gas have been carried out in Brown Shale, Middle Pematang Formation, Central Sumatra Basin. The paper is aimed at delineating the sweet spot distribution of potential shale gas reservoir, which is based on Total Organic Carbon (TOC), Maturity level data, and combined with TOC modeling that refers to Passey and Regression Multi Linear method. We used 4 well data, side wall core and 3D pre-stack seismic data. Our analysis of geochemical properties is based on well log and core data and its distribution are constrained by a framework of 3D seismic data, which is transformed into acoustic impedance. Further, the sweet spot of organic-rich shale is delineated by mapping TOC, which is extracted from inverted acoustic impedance. Our experiment analysis shows that organic materials contained in the formation of Middle Pematang Brown Shale members have TOC range from 0.15 to 2.71 wt.%, which is classified in the quality of poor to very good. In addition, the maturity level of organic material is ranging from 373°C to 432°C, which is indicated by vitrinite reflectance (Ro) of 0.58. In term of kerogen type, this Brown shale formation is categorized as kerogen type of II I III, which has the potential to generate a mixture of gasIoil on the environment.

Snowmelt induced hydrologic perturbations drive dynamic microbiological and geochemical behaviors across a shallow riparian aquifer

NASA Astrophysics Data System (ADS)

Danczak, Robert; Yabusaki, Steven; Williams, Kenneth; Fang, Yilin; Hobson, Chad; Wilkins, Michael

2016-05-01

Shallow riparian aquifers represent hotspots of biogeochemical activity in the arid western US. While these environments provide extensive ecosystem services, little is known of how natural environmental perturbations influence subsurface microbial communities and associated biogeochemical processes. Over a six-month period we tracked the annual snowmelt-driven incursion of groundwater into the vadose zone of an aquifer adjacent to the Colorado River, leading to increased dissolved oxygen (DO) concentrations in the normally suboxic saturated zone. Strong biogeochemical heterogeneity was measured across the site, with abiotic reactions between DO and sulfide minerals driving rapid DO consumption and mobilization of redox active species in reduced aquifer regions. Conversely, extensive DO increases were detected in less reduced sediments. 16S rRNA gene surveys tracked microbial community composition within the aquifer, revealing strong correlations between increases in putative oxygen-utilizing chemolithoautotrophs and heterotrophs and rising DO concentrations. The gradual return to suboxic aquifer conditions favored increasing abundances of 16S rRNA sequences matching members of the Microgenomates (OP11) and Parcubacteria (OD1) that have been strongly implicated in fermentative processes. Microbial community stability measurements indicated that deeper aquifer locations were relatively less affected by geochemical perturbations, while communities in shallower locations exhibited the greatest change. Reactive transport modeling of the geochemical and microbiological results supported field observations, suggesting that a predictive framework can be applied to develop a greater understanding of such environments.

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

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

Chen, Jinsong

2013-05-01

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

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

GoPhast: a graphical user interface for PHAST

USGS Publications Warehouse

Winston, Richard B.

2006-01-01

GoPhast is a graphical user interface (GUI) for the USGS model PHAST. PHAST simulates multicomponent, reactive solute transport in three-dimensional, saturated, ground-water flow systems. PHAST can model both equilibrium and kinetic geochemical reactions. PHAST is derived from HST3D (flow and transport) and PHREEQC (geochemical calculations). The flow and transport calculations are restricted to constant fluid density and constant temperature. The complexity of the input required by PHAST makes manual construction of its input files tedious and error-prone. GoPhast streamlines the creation of the input file and helps reduce errors. GoPhast allows the user to define the spatial input for the PHAST flow and transport data file by drawing points, lines, or polygons on top, front, and side views of the model domain. These objects can have up to two associated formulas that define their extent perpendicular to the view plane, allowing the objects to be three-dimensional. Formulas are also used to specify the values of spatial data (data sets) both globally and for individual objects. Objects can be used to specify the values of data sets independent of the spatial and temporal discretization of the model. Thus, the grid and simulation periods for the model can be changed without respecifying spatial data pertaining to the hydrogeologic framework and boundary conditions. This report describes the operation of GoPhast and demonstrates its use with examples. GoPhast runs on Windows 2000, Windows XP, and Linux operating systems.

Contaminant source identification using semi-supervised machine learning

NASA Astrophysics Data System (ADS)

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

2018-05-01

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

Contaminant source identification using semi-supervised machine learning

DOE PAGES

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

2017-11-08

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

Contaminant source identification using semi-supervised machine learning

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

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

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

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

USGS Publications Warehouse

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

2014-01-01

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

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

NASA Astrophysics Data System (ADS)

Ghezelbash, Reza; Maghsoudi, Abbas

2018-05-01

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

The influence of precipitation kinetics on trace element partitioning between solid and liquid solutions: A coupled fluid dynamics/thermodynamics framework to predict distribution coefficients

NASA Astrophysics Data System (ADS)

Kavner, A.

2017-12-01

In a multicomponent multiphase geochemical system undergoing a chemical reaction such as precipitation and/or dissolution, the partitioning of species between phases is determined by a combination of thermodynamic properties and transport processes. The interpretation of the observed distribution of trace elements requires models integrating coupled chemistry and mechanical transport. Here, a framework is presented that predicts the kinetic effects on the distribution of species between two reacting phases. Based on a perturbation theory combining Navier-Stokes fluid flow and chemical reactivity, the framework predicts rate-dependent partition coefficients in a variety of different systems. We present the theoretical framework, with applications to two systems: 1. species- and isotope-dependent Soret diffusion of species in a multicomponent silicate melt subjected to a temperature gradient, and 2. Elemental partitioning and isotope fractionation during precipitation of a multicomponent solid from a multicomponent liquid phase. Predictions will be compared with results from experimental studies. The approach has applications for understanding chemical exchange in at boundary layers such as the Earth's surface magmatic systems and at the core/mantle boundary.

A compendium of geochemical information from the Saanich Inlet water column

NASA Astrophysics Data System (ADS)

Torres-Beltrán, Mónica; Hawley, Alyse K.; Capelle, David; Zaikova, Elena; Walsh, David A.; Mueller, Andreas; Scofield, Melanie; Payne, Chris; Pakhomova, Larysa; Kheirandish, Sam; Finke, Jan; Bhatia, Maya; Shevchuk, Olena; Gies, Esther A.; Fairley, Diane; Michiels, Céline; Suttle, Curtis A.; Whitney, Frank; Crowe, Sean A.; Tortell, Philippe D.; Hallam, Steven J.

2017-10-01

Extensive and expanding oxygen minimum zones (OMZs) exist at variable depths in coastal and open ocean waters. As oxygen levels decline, nutrients and energy are increasingly diverted away from higher trophic levels into microbial community metabolism, resulting in fixed nitrogen loss and production of climate active trace gases including nitrous oxide and methane. While ocean deoxygenation has been reported on a global scale, our understanding of OMZ biology and geochemistry is limited by a lack of time-resolved data sets. Here, we present a historical dataset of oxygen concentrations spanning fifty years and nine years of monthly geochemical time series observations in Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia, Canada that undergoes recurring changes in water column oxygenation status. This compendium provides a unique geochemical framework for evaluating long-term trends in biogeochemical cycling in OMZ waters.

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.

Water age and stream solute dynamics at the Hubbard Brook Experimental Forest (US)

NASA Astrophysics Data System (ADS)

Botter, Gianluca; Benettin, Paolo; McGuire, Kevin; Rinaldo, Andrea

2016-04-01

The contribution discusses experimental and modeling results from a headwater catchment at the Hubbard Brook Experimental Forest (New Hampshire, USA) to explore the link between stream solute dynamics and water age. A theoretical framework based on water age dynamics, which represents a general basis for characterizing solute transport at the catchment scale, is used to model both conservative and weathering-derived solutes. Based on the available information about the hydrology of the site, an integrated transport model was developed and used to estimate the relevant hydrochemical fluxes. The model was designed to reproduce the deuterium content of streamflow and allowed for the estimate of catchment water storage and dynamic travel time distributions (TTDs). Within this framework, dissolved silicon and sodium concentration in streamflow were simulated by implementing first-order chemical kinetics based explicitly on dynamic TTD, thus upscaling local geochemical processes to catchment scale. Our results highlight the key role of water stored within the subsoil glacial material in both the short-term and long-term solute circulation at Hubbard Brook. The analysis of the results provided by the calibrated model allowed a robust estimate of the emerging concentration-discharge relationship, streamflow age distributions (including the fraction of event water) and storage size, and their evolution in time due to hydrologic variability.

Snowmelt Induced Hydrologic Perturbations Drive Dynamic Microbiological and Geochemical Behaviors across a Shallow Riparian Aquifer

DOE PAGES

Danczak, Robert E.; Yabusaki, Steven B.; Williams, Kenneth H.; ...

2016-05-11

Shallow riparian aquifers represent hotspots of biogeochemical activity in the arid western US. While these environments provide extensive ecosystem services, little is known of how natural environmental perturbations influence subsurface microbial communities and associated biogeochemical processes. Over a 6-month period we tracked the annual snowmelt-driven incursion of groundwater into the vadose zone of an aquifer adjacent to the Colorado River, leading to increased dissolved oxygen (DO) concentrations in the normally suboxic saturated zone. Strong biogeochemical heterogeneity was measured across the site, with abiotic reactions between DO and sulfide minerals driving rapid DO consumption and mobilization of redox active species inmore » reduced aquifer regions. Conversely, extensive DO increases were detected in less reduced sediments. 16S rRNA gene surveys tracked microbial community composition within the aquifer, revealing strong correlations between increases in putative oxygen-utilizing chemolithoautotrophs and heterotrophs and rising DO concentrations. The gradual return to suboxic aquifer conditions favored increasing abundances of 16S rRNA sequences matching members of the Microgenomates (OP11) and Parcubacteria (OD1) that have been strongly implicated in fermentative processes. Microbial community stability measurements indicated that deeper aquifer locations were relatively less affected by geochemical perturbations, while communities in shallower locations exhibited the greatest change. Thus, reactive transport modeling of the geochemical and microbiological results supported field observations, suggesting that a predictive framework can be applied to develop a greater understanding of such environments.« less

Paleoclimatological study using stalagmites from Java Island, Indonesia

NASA Astrophysics Data System (ADS)

Watanabe, Y.; Matsuoka, H.; Ohsawa, S.; Yamada, M.; Kitaoka, K.; Kiguchi, M.; Ueda, J.; Yoshimura, K.; Kurisaki, K.; Nakai, S.; Brahmantyo, B.; Maryunani, K. A.; Tagami, T.; Takemura, K.; Yoden, S.

2006-12-01

In the last decade, decoding geochemical records in stalagmites has been widely recognized as a powerful tool for the elucidation of paleoclimate/environment of the terrestrial areas. The previous data are mainly reported from areas that are located in middle latitude. However, this study aims at reconstructing past climate variations in the Asian equatorial regions by using oxygen isotopes and other geochemical proxies recorded in Indonesian stalagmites.. Especially, we focus on the detection of the precipitation anomaly that reflects the El Niño Southern Oscillation (ENSO). We performed geological surveys in Buniayu limestone caves, Sukabumi, West Java, and Karangbolong, Central Java, Indonesia and collected a series of stalagmites/stalactites and drip water samples. Detailed textures of stalagmite samples were observed using thin sections to identify "annual" bandings. Moreover, we also measured both (1) annual luminescent banding that can be viewed by ultraviolet-light stimulation and (2) uranium series disequilibrium ages using the MC-ICP-MS for each stalagmite to construct the age model. We also carried out 3H-3He dating and stable isotope measurements of drip water samples to understand hydrogeology in study areas. Based on these frameworks, oxygen isotopes and other geochemical proxies will be analyzed for annual or sub-annual time scales. The proxy data will then be compared with meteorological data set, such as local precipitation, in the past 50 years. Finally, we will reconstruct for longer timescales the past climate, particularly the precipitation anomaly, in the region to detect ancient ENSO.

Snowmelt Induced Hydrologic Perturbations Drive Dynamic Microbiological and Geochemical Behaviors across a Shallow Riparian Aquifer

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

Danczak, Robert E.; Yabusaki, Steven B.; Williams, Kenneth H.

Shallow riparian aquifers represent hotspots of biogeochemical activity in the arid western US. While these environments provide extensive ecosystem services, little is known of how natural environmental perturbations influence subsurface microbial communities and associated biogeochemical processes. Over a 6-month period we tracked the annual snowmelt-driven incursion of groundwater into the vadose zone of an aquifer adjacent to the Colorado River, leading to increased dissolved oxygen (DO) concentrations in the normally suboxic saturated zone. Strong biogeochemical heterogeneity was measured across the site, with abiotic reactions between DO and sulfide minerals driving rapid DO consumption and mobilization of redox active species inmore » reduced aquifer regions. Conversely, extensive DO increases were detected in less reduced sediments. 16S rRNA gene surveys tracked microbial community composition within the aquifer, revealing strong correlations between increases in putative oxygen-utilizing chemolithoautotrophs and heterotrophs and rising DO concentrations. The gradual return to suboxic aquifer conditions favored increasing abundances of 16S rRNA sequences matching members of the Microgenomates (OP11) and Parcubacteria (OD1) that have been strongly implicated in fermentative processes. Microbial community stability measurements indicated that deeper aquifer locations were relatively less affected by geochemical perturbations, while communities in shallower locations exhibited the greatest change. Thus, reactive transport modeling of the geochemical and microbiological results supported field observations, suggesting that a predictive framework can be applied to develop a greater understanding of such environments.« less

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

EPA Science Inventory

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

Evolution of spreading rate and H2 production by serpentinization at mid-ocean ridges from 200 Ma to Present

NASA Astrophysics Data System (ADS)

Andreani, M.; García del Real, P.; Daniel, I.; Wright, N.; Coltice, N.

2017-12-01

Mid-oceanic ridge (MOR) spreading rate spatially varies today from 20 to 200 mm/yr and geological records attest of important temporal variations, at least during the past 200 My. The spreading rate has a direct impact on the mechanisms accomodating extension (magmatic vs tectonic), hence on the nature of the rocks forming the oceanic lithosphere. The latter is composed of variable amount of magmatic and mantle rocks, that dominate at fast and (ultra-) slow spreading ridges, respectively. Serpentinization of mantle rocks contributes to global fluxes and notably to those of hydrogen and carbon by providing a pathways for dihydrogen (H2) production, carbone storage by mineralization, and carbon reduction to CH4 and possibly complex organic compounds. Quantification of the global chemical impact of serpentinization through geological time requires a coupling of geochemical parameters with plate-tectonic reconstructions. Here we quantify serpentinization extent and concurrent H2 production at MOR from the Jurassic (200 Ma) to present day (0 Ma). We coupled mean values of relevant petro-chemical parameters such as the proportion of mantle rocks, initial iron in olivine, iron redox state in serpentinites, % of serpentinization to high-resolution models of plate motion within the GPlates infrastructure to estimate the lengths in 1 Myr intervals for the global MOR plate boundary (spreading and transform components), and spreading ridges as a function of their rate. The model sensitivity to selected parameters has been tested. The results show that fragmentation of Pangea resulted in elevated H2 rates (>1012 to 1013 mol/yr) starting at 160 Ma compared to Late Mesozoic (<160 Ma) rates (<1011-1012 mol/yr). From 160 Ma to present, the coupled opening of the Atlantic and Indian oceans as well as the variation in spreading rates maintained H2 generation in the 1012 mol/yr level, but with significant excursions mainly related to the length of ultra-slow spreading segments. For the first time, this model offers a framework toward flux modeling at MOR through time. The model can be further implemented by adding supplementary geochemical parameters or serve other geochemical issues.

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

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

NASA Astrophysics Data System (ADS)

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

2016-04-01

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

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.

HYDROGIOLOGIC FRAMEWORK, GROUND-WATER GEOCHEMISTRY, AND ASSESSMENT OF NITROGEN YIELD FROM BASE FLOW IN TWO AGRICULTURAL WATERSHEDS, KENT COUNTY, MARYLAND

EPA Science Inventory

Hydrostratigraphic and geochemical data collected in two adjacent watersheds on the Delmarva Peninsula, in Kent County, Maryland, indicate that shallow subsurface stratigraphy is an important factor that affects the concentrations of nitrogen in ground water discharging as stream...

Biological reduction of chlorinated solvents: Batch-scale geochemical modeling

NASA Astrophysics Data System (ADS)

Kouznetsova, Irina; Mao, Xiaomin; Robinson, Clare; Barry, D. A.; Gerhard, Jason I.; McCarty, Perry L.

2010-09-01

Simulation of biodegradation of chlorinated solvents in dense non-aqueous phase liquid (DNAPL) source zones requires a model that accounts for the complexity of processes involved and that is consistent with available laboratory studies. This paper describes such a comprehensive modeling framework that includes microbially mediated degradation processes, microbial population growth and decay, geochemical reactions, as well as interphase mass transfer processes such as DNAPL dissolution, gas formation and mineral precipitation/dissolution. All these processes can be in equilibrium or kinetically controlled. A batch modeling example was presented where the degradation of trichloroethene (TCE) and its byproducts and concomitant reactions (e.g., electron donor fermentation, sulfate reduction, pH buffering by calcite dissolution) were simulated. Local and global sensitivity analysis techniques were applied to delineate the dominant model parameters and processes. Sensitivity analysis indicated that accurate values for parameters related to dichloroethene (DCE) and vinyl chloride (VC) degradation (i.e., DCE and VC maximum utilization rates, yield due to DCE utilization, decay rate for DCE/VC dechlorinators) are important for prediction of the overall dechlorination time. These parameters influence the maximum growth rate of the DCE and VC dechlorinating microorganisms and, thus, the time required for a small initial population to reach a sufficient concentration to significantly affect the overall rate of dechlorination. Self-inhibition of chlorinated ethenes at high concentrations and natural buffering provided by the sediment were also shown to significantly influence the dechlorination time. Furthermore, the analysis indicated that the rates of the competing, nonchlorinated electron-accepting processes relative to the dechlorination kinetics also affect the overall dechlorination time. Results demonstrated that the model developed is a flexible research tool that is able to provide valuable insight into the fundamental processes and their complex interactions during bioremediation of chlorinated ethenes in DNAPL source zones.

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.

Orbital forcing of climate 1.4 billion years ago.

PubMed

Zhang, Shuichang; Wang, Xiaomei; Hammarlund, Emma U; Wang, Huajian; Costa, M Mafalda; Bjerrum, Christian J; Connelly, James N; Zhang, Baomin; Bian, Lizeng; Canfield, Donald E

2015-03-24

Fluctuating climate is a hallmark of Earth. As one transcends deep into Earth time, however, both the evidence for and the causes of climate change become difficult to establish. We report geochemical and sedimentological evidence for repeated, short-term climate fluctuations from the exceptionally well-preserved ∼1.4-billion-year-old Xiamaling Formation of the North China Craton. We observe two patterns of climate fluctuations: On long time scales, over what amounts to tens of millions of years, sediments of the Xiamaling Formation record changes in geochemistry consistent with long-term changes in the location of the Xiamaling relative to the position of the Intertropical Convergence Zone. On shorter time scales, and within a precisely calibrated stratigraphic framework, cyclicity in sediment geochemical dynamics is consistent with orbital control. In particular, sediment geochemical fluctuations reflect what appear to be orbitally forced changes in wind patterns and ocean circulation as they influenced rates of organic carbon flux, trace metal accumulation, and the source of detrital particles to the sediment.

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

EPA Science Inventory

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

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

NASA Astrophysics Data System (ADS)

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

2007-12-01

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

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

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

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

2003-09-01

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

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

NASA Astrophysics Data System (ADS)

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

2018-05-01

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

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.

Stratigraphic framework of Cambrian and Ordovician rocks in the central Appalachian basin from Campbell County, Kentucky, to Tazewell County, Virginia: Chapter E.2.4 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Ryder, Robert T.; Repetski, John E.; Harris, Anita G.; Lentz, Erika E.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

This chapter is a re-release of U.S. Geological Survey Miscellaneous Investigations Series Map I-2530, of the same title, by Ryder and others (1997; online version 2.0 revised and digitized by Erika E. Lentz, 2004). Version 2.0 is a digital version of the original and also includes the gamma-ray well log traces.

Geochemical Reaction Mechanism Discovery from Molecular Simulation

DOE PAGES

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

2014-11-10

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

Modeling Low-temperature Geochemical Processes

NASA Astrophysics Data System (ADS)

Nordstrom, D. K.

2003-12-01

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

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

USGS Publications Warehouse

Han, L. F; Plummer, Niel

2016-01-01

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

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.

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

DTIC Science & Technology

2006-09-01

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

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

EPA Science Inventory

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

A kinetics database and scripts for PHREEQC

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

Quantifying Risks and Uncertainties Associated with Induced Seismicity due to CO2 Injection into Geologic Formations with Faults

NASA Astrophysics Data System (ADS)

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

2016-12-01

A multiphase flow and reactive transport simulator named STOMP-CO2-R has been developed and coupled to the ABAQUS® finite element package for geomechanical analysis enabling comprehensive thermo-hydro-geochemical-mechanical (THMC) analyses. The coupled THMC simulator has been applied to analyze faulted CO2 reservoir responses (e.g., stress and strain distributions, pressure buildup, slip tendency factor, pressure margin to fracture) with various complexities in fault and reservoir structures and mineralogy. Depending on the geological and reaction network settings, long-term injection of CO2 can have a significant effect on the elastic stiffness and permeability of formation rocks. In parallel, an uncertainty quantification framework (UQ-CO2), which consists of entropy-based prior uncertainty representation, efficient sampling, geostatistical reservoir modeling, and effective response surface analysis, has been developed for quantifying risks and uncertainties associated with CO2 sequestration. It has been demonstrated for evaluating risks in CO2 leakage through natural pathways and wellbores, and for developing predictive reduced order models. Recently, a parallel STOMP-CO2-R has been developed and the updated STOMP/ABAQUS model has been proven to have a great scalability, which makes it possible to integrate the model with the UQ framework to effectively and efficiently explore multidimensional parameter space (e.g., permeability, elastic modulus, crack orientation, fault friction coefficient) for a more systematic analysis of induced seismicity risks.

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

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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

NASA Astrophysics Data System (ADS)

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

2016-07-01

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

Geoneutrinos and Heat Production in the Earth: Constraints and Implications

ScienceCinema

McDonough, Bill

2017-12-29

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

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

USGS Publications Warehouse

Ellefsen, Karl J.; Smith, David

2016-01-01

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

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.

Archean crust-mantle geochemical differentiation

NASA Astrophysics Data System (ADS)

Tilton, G. R.

Isotope measurements on carbonatite complexes and komatiites can provide information on the geochemical character and geochemical evolution of the mantle, including the sub-continental mantle. Measurements on young samples establish the validity of the method. These are based on Sr, Nd and Pb data from the Tertiary-Mesozoic Gorgona komatiite and Sr and Pb data from the Cretaceous Oka carbonatite complex. In both cases the data describe a LIL element-depleted source similar to that observed presently in MORB. Carbonatite data have been used to study the mantle beneath the Superior Province of the Canadian Shield one billion years (1 AE) ago. The framework for this investigation was established by Bell et al., who showed that large areas of the province appear to be underlain by LIL element-depleted mantle (Sr-85/Sr-86=0.7028) at 1 AE ago. Additionally Bell et al. found four complexes to have higher initial Sr ratios (Sr-87/Sr-86=0.7038), which they correlated with less depleted (bulk earth?) mantle sources, or possibly crustal contamination. Pb isotope relationships in four of the complexes have been studied by Bell et al.

Archean crust-mantle geochemical differentiation

NASA Technical Reports Server (NTRS)

Tilton, G. R.

1983-01-01

Isotope measurements on carbonatite complexes and komatiites can provide information on the geochemical character and geochemical evolution of the mantle, including the sub-continental mantle. Measurements on young samples establish the validity of the method. These are based on Sr, Nd and Pb data from the Tertiary-Mesozoic Gorgona komatiite and Sr and Pb data from the Cretaceous Oka carbonatite complex. In both cases the data describe a LIL element-depleted source similar to that observed presently in MORB. Carbonatite data have been used to study the mantle beneath the Superior Province of the Canadian Shield one billion years (1 AE) ago. The framework for this investigation was established by Bell et al., who showed that large areas of the province appear to be underlain by LIL element-depleted mantle (Sr-85/Sr-86=0.7028) at 1 AE ago. Additionally Bell et al. found four complexes to have higher initial Sr ratios (Sr-87/Sr-86=0.7038), which they correlated with less depleted (bulk earth?) mantle sources, or possibly crustal contamination. Pb isotope relationships in four of the complexes have been studied by Bell et al.

Alchemy in the underworld - recent progress and future potential of organic geochemistry applied to speleothems.

NASA Astrophysics Data System (ADS)

Blyth, Alison

2016-04-01

Speleothems are well used archives for chemical records of terrestrial environmental change, and the integration of records from a range of isotopic, inorganic, and organic geochemical techniques offers significant power in reconstructing both changes in past climates and identifying the resultant response in the overlying terrestrial ecosystems. The use of organic geochemistry in this field offers the opportunity to recover new records of vegetation change (via biomarkers and compound specific isotopes), temperature change (via analysis of glycerol dialkyl glycerol tetraethers, a compound group derived from microbes and varying in structure in response to temperature and pH), and changes in soil microbial behaviour (via combined carbon isotope analysis). However, to date the use of organic geochemical techniques has been relatively limited, due to issues relating to sample size, concerns about contamination, and unanswered questions about the origins of the preserved organic matter and rates of transport. Here I will briefly review recent progress in the field, and present a framework for the future research needed to establish organic geochemical analysis in speleothems as a robust palaeo-proxy approach.

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.

Snowmelt Induced Hydrologic Perturbations Drive Dynamic Microbiological and Geochemical Behaviors across a Shallow Riparian Aquifer

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

Danczak, Robert E.; Yabusaki, Steven B.; Williams, Kenneth H.

Shallow riparian aquifers represent hotspots of biogeochemical activity in the arid western US. While these environments provide extensive ecosystem services, little is known of how natural environmental perturbations influence subsurface microbial communities and associated biogeochemical processes. Over a six-month period we tracked the annual snowmelt-driven incursion of groundwater into the vadose zone of an aquifer adjacent to the Colorado River, leading to increased dissolved oxygen (DO) concentrations in the normally suboxic saturated zone. Strong biogeochemical heterogeneity was measured across the site, with abiotic reactions between DO and sulfide minerals driving rapid DO consumption and mobilization of redox active species inmore » reduced aquifer regions. Conversely, extensive DO increases were detected in less reduced sediments. 16S rRNA gene surveys tracked microbial community composition within the aquifer, revealing strong correlations between increases in putative oxygen-utilizing chemolithoautotrophs and heterotrophs and rising DO concentrations. The gradual return to suboxic aquifer conditions favored increasing abundances of 16S rRNA sequences matching members of the Microgenomates (OP11) and Parcubacteria (OD1) that have been strongly implicated in fermentative processes. Microbial community stability measurements indicated that deeper aquifer locations were relatively less affected by geochemical perturbations, while communities in shallower locations exhibited the greatest change. Reactive transport modeling of the geochemical and microbiological results supported field observations, suggesting that a predictive framework can be applied to develop a greater understanding of such environments. Frontiers in Earth Science Journal Impact & Description - ResearchGate - Impact Rankings ( 2015 and 2016 ). Available from: https://www.researchgate.net/journal/2296-6463_Frontiers_in_Earth_Science [accessed Jul 25, 2016].« less

Potential accumulation of contaminated sediments in a reservoir of a high-Andean watershed: Morphodynamic connections with geochemical processes

NASA Astrophysics Data System (ADS)

Contreras, María. Teresa; Müllendorff, Daniel; Pastén, Pablo; Pizarro, Gonzalo E.; Paola, Chris; Escauriaza, Cristián.

2015-05-01

Rapid changes due to anthropic interventions in high-altitude environments, such as the Altiplano region in South America, require new approaches to understand the connections between physical and geochemical processes. Alterations of the water quality linked to the river morphology can affect the ecosystems and human development in the long term. The future construction of a reservoir in the Lluta River, located in northern Chile, will change the spatial distribution of arsenic-rich sediments, which can have significant effects on the lower parts of the watershed. In this investigation, we develop a coupled numerical model to predict and evaluate the interactions between morphodynamic changes in the Lluta reservoir, and conditions that can potentially desorb arsenic from the sediments. Assuming that contaminants are mobilized under anaerobic conditions, we calculate the oxygen concentration within the sediments to study the interactions of the delta progradation with the potential arsenic release. This work provides a framework for future studies aimed to analyze the complex connections between morphodynamics and water quality, when contaminant-rich sediments accumulate in a reservoir. The tool can also help to design effective risk management and remediation strategies in these extreme environments. This article was corrected on 15 JUNE 2015. See the end of the full text for details.

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.

Geothermal systems: Principles and case histories

NASA Astrophysics Data System (ADS)

Rybach, L.; Muffler, L. J. P.

The classification of geothermal systems is considered along with the geophysical and geochemical signatures of geothermal systems, aspects of conductive heat transfer and regional heat flow, and geothermal anomalies and their plate tectonic framework. An investigation of convective heat and mass transfer in hydrothermal systems is conducted, taking into account the mathematical modelling of hydrothermal systems, aspects of idealized convective heat and mass transport, plausible models of geothermal reservoirs, and preproduction models of hydrothermal systems. Attention is given to the prospecting for geothermal resources, the application of water geochemistry to geothermal exploration and reservoir engineering, heat extraction from geothermal reservoirs, questions of geothermal resource assessment, and environmental aspects of geothermal energy development. A description is presented of a number of case histories, taking into account the low enthalpy geothermal resource of the Pannonian Basin in Hungary, the Krafla geothermal field in Northeast Iceland, the geothermal system of the Jemez Mountains in New Mexico, and extraction-reinjection at the Ahuachapan geothermal field in El Salvador.

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.

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

NASA Astrophysics Data System (ADS)

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

2017-04-01

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

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

PubMed

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

2014-01-01

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

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

NASA Technical Reports Server (NTRS)

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

2017-01-01

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

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

NASA Astrophysics Data System (ADS)

Gonçalves, Mario; Mateus, Antonio

2016-04-01

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

Boron Isotopic Composition Correlates with Ultra-Structure in a - Sea Coral Lophelia Pertusa: Implications for Biomineralization and - PH

NASA Astrophysics Data System (ADS)

Blamart, D.; Rollion-Bard, C.; Meibom, A.; Cuif, J.; Juillet-Leclerc, A.; Dauphin, Y.; Douarin, M.

2007-12-01

The geochemistry (stable isotopes and trace elements) of biogenic carbonates has been widely used for more than fifty years to reconstruct past climatic variability. During this time, the studies were mainly based on bulk sampling limiting sometimes the interpretations of the geochemical data as paleoclimatic proxies. Recently, high spatial resolution sampling techniques, such as micro-mill and SIMS, have been employed in the study of C, O and B isotopic compositions and trace elements (Mg, Sr) in the skeletons of a variety of (deep-sea) coral species. These studies have documented dramatic 'vital effects' and uncovered a systematic relationship between skeletal ultra-structure and stable isotopic composition. The formation of skeleton corals follows a universal two-step growth process. At the tips of the skeletal structures, the mineralizing cell layer produces centers of calcification (COC) or, equivalently, Early Mineralization Zone (EMZ). These EMZ are subsequently overgrown by fibrous aragonite(FA) consisting of cyclically added layers. The EMZ are characterized by systematically lighter C and O isotopic compositions compared with the adjacent FA. A number of geochemical models have been proposed, in which this systematic stable isotopic difference between EMZ and FA is ascribed to a biologically induced variation in the pH of a proposed Extra-cytoplasmic Calcifying Fluid (ECF) reservoir. In these models, relatively high pH conditions during the formation of EMZ result in relatively light C and O isotopic compositions compared with FA, which form under generally lower pH conditions. A direct test of such models would be possible if the Boron isotopic composition, which is pH sensitive, of EMZ and FA could be measured. We performed ion microprobe d11B measurements for EMZ and FA in Lophelia pertusa, a deep-sea coral common in the North-East Atlantic Ocean. We observe a systematic difference in B isotopic composition between the EMZ and FA skeleton. In EMZ, the measured δ11B values are consistently low. Fibrous aragonite is characterized by systematically higher d11B values, but also display B isotopic heterogeneity associated with specific growth bands in the calyx wall. The magnitude of the observed B isotopic variations cannot be explained by changes in environmental conditions and are likely caused by biological processes involved in the biomineralization of new skeleton; i.e. 'vital' effects. The observed B isotopic variations are opposite to the predictions of geochemical models for vital effects. Our data indicate that pH variations are not responsible for the observed stable isotopic fractionations. Geochemical models therefore do not provide an adequate framework within which to understand coral skeletal formation. Without a better understanding of these processes, which require experiments, the use of B isotopic composition to reconstruct paleo-pH variations in the oceans must be considered problematic - at least as far as Lophelia pertusa is concerned.

Executive summary: Chapter A.1 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Ruppert, Leslie F.; Ryder, Robert T.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

This publication supplements and updates older USGS regional studies of Appalachian basin coal and petroleum resources such as those by Arndt and others (1968) and the numerous contributors to USGS Miscellaneous Map Series I−917 (for example, Harris and others, 1978), respectively. USGS Professional Paper 1708 is intended primarily for geoscientists in academia, industry, and government who are interested in Appalachian basin geology and its coal and petroleum resources. Other users, however, may find the wide variety of topics, papers, and digital images of value for landuse and policy planning issues. Among the anticipated benefits of the report are improvements in (1) resource assessment estimates and methodology, (2) exploration strategies, (3) basin models, and (4) energy use policies.

Hydraulic and geochemical framework of the Idaho National Engineering and Environmental Laboratory vadose zone

USGS Publications Warehouse

Nimmo, John R.; Rousseau, Joseph P.; Perkins, Kim S.; Stollenwerk, Kenneth G.; Glynn, Pierre D.; Bartholomay, Roy C.; Knobel, LeRoy L.

2004-01-01

Questions of major importance for subsurface contaminant transport at the Idaho National Engineering and Environmental Laboratory (INEEL) include (i) travel times to the aquifer, both average or typical values and the range of values to be expected, and (ii) modes of contaminant transport, especially sorption processes. The hydraulic and geochemical framework within which these questions are addressed is dominated by extreme heterogeneity in a vadose zone and aquifer consisting of interbedded basalts and sediments. Hydraulically, major issues include diverse possible types of flow pathways, extreme anisotropy, preferential flow, combined vertical and horizontal flow, and temporary saturation or perching. Geochemically, major issues include contaminant mobility as influenced by redox conditions, the concentration of organic and inorganic complexing solutes and other local variables, the interaction with infiltrating waters and with the contaminant source environment, and the aqueous speciation of contaminants such as actinides. Another major issue is the possibility of colloid transport, which inverts some of the traditional concepts of mobility, as sorbed contaminants on mobile colloids may be transported with ease compared with contaminants that are not sorbed. With respect to the goal of minimizing aquifer concentrations of contaminants, some characteristics of the vadose zone are essentially completely favorable. Examples include the great thickness (200 m) of the vadose zone, and the presence of substantial quantities of fine sediments that can retard contaminant transport both hydraulically and chemically. Most characteristics, however, have both favorable and unfavorable aspects. For example, preferential flow, as promoted by several notable features of the vadose zone at the INEEL, can provide fast, minimally sorbing pathways for contaminants to reach the aquifer easily, but it also leads to a wide dispersal of contaminants in a large volume of subsurface material, thus increasing the opportunity for dilution and sorption.

Segmentation of singularity maps in the context of soil porosity

NASA Astrophysics Data System (ADS)

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

2016-04-01

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

A geochemical examination of humidity cell tests

USGS Publications Warehouse

Maest, Ann; Nordstrom, D. Kirk

2017-01-01

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

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

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

USGS Publications Warehouse

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

2009-01-01

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

Can treatment and disposal costs be reduced through metal recovery?

USGS Publications Warehouse

Smith, Kathleen S.; Figueroa, Linda; Plumlee, Geoffrey S.

2015-01-01

This paper describes a framework to conduct a “metal-recovery feasibility assessment” for mining influenced water (MIW) and associated treatment sludge. There are multiple considerations in such a determination, including the geologic/geochemical feasibility, market feasibility, technical feasibility, economic feasibility, and administrative feasibility. Each of these considerations needs to be evaluated to determine the practicality of metal recovery from a particular MIW.

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.

Coal and coalbed-methane resources in the Appalachian and Black Warrior basins: maps showing the distribution of coal fields, coal beds, and coalbed-methane fields: Chapter D.1 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Trippi, Michael H.; Ruppert, Leslie F.; Milici, Robert C.; Kinney, Scott A.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

The study area for most reports in this volume is the Appalachian basin. The term “Appalachian basin study area” (shortened from “Appalachian basin geologic framework study area”) includes all of the Appalachian Basin Province (Province 67) and part of the neighboring Black Warrior Basin Province (Province 65) of Dolton and others (1995). The boundaries for these two provinces and the study area are shown on figure 1.

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.

Composition of natural gas and crude oil produced from 10 wells in the Lower Silurian "Clinton" Sandstone, Trumbull County, Ohio: Chapter G.7 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Burruss, Robert A.; Ryder, Robert T.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

Natural gases and associated crude oils in the “Clinton” sandstone, Medina Group sandstones, and equivalent Tuscarora Sandstone in the northern Appalachian basin are part of a regional, continuous-type or basin-centered accumulation. The origin of the hydrocarbon charge to regional continuoustype accumulations is poorly understood. We have analyzed the molecular and stable isotopic composition of gases and oils produced from 10 wells in the “Clinton” sandstone in Trumbull County, Ohio, in an initial attempt to identify the characteristics of the accumulated fluids. The analyses show that the fluids have remarkably uniform compositions that are similar to previously published analyses of oils (Cole and others, 1987) and gases (Laughrey and Baldasarre, 1998) in Early Silurian reservoirs elsewhere in Ohio; however, geochemical parameters in the oils and gases suggest that the fluids have experienced higher levels of thermal stress than the present-day burial conditions of the reservoir rocks. The crude oils have an unusual geochemical characteristic: they do not contain detectable levels of sterane and triterpane biomarkers. The origin of these absences is unknown.

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

NASA Astrophysics Data System (ADS)

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

2018-05-01

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

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

NASA Technical Reports Server (NTRS)

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

1982-01-01

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

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

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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

NASA Technical Reports Server (NTRS)

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

2015-01-01

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

Calculation of individual isotope equilibrium constants for implementation in geochemical models

USGS Publications Warehouse

Thorstenson, Donald C.; Parkhurst, David L.

2002-01-01

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

Tracing Marine Cryptotephras in the North Atlantic during the Last Glacial Period

NASA Astrophysics Data System (ADS)

Abbott, Peter; Davies, Siwan; Griggs, Adam; Bourne, Anna

2017-04-01

Tephrochronology is a powerful technique that can be utilised for the independent correlation and synchronisation of disparate palaeoclimatic records from different depositional environments. There is a high potential to utilise this technique to integrate ice, marine and terrestrial records to study climatic phasing within the North Atlantic region due to the high eruptive frequency of Icelandic volcanic systems. However, until now North Atlantic marine records have been relatively understudied. Here we report on investigations to define a tephra framework integrating new studies of cryptotephra horizons within a wide network of North Atlantic marine cores with horizons identified in prior work. This framework has the potential to underpin the correlation of the marine records to the Greenland ice-core records and European terrestrial sequences. Tephrochronological investigations were conducted on 13 marine sequences from a range of locations and depositional settings using cryptotephra extraction techniques, including density and magnetic separation, to gain high resolution glass shard concentration profiles and rigorous single-shard major element geochemical analysis to characterise identified deposits. Cryptotephras with an Icelandic source were identified in many records and displayed diversity in shard concentration profiles and the geochemical homo/heterogeneity of shards within the deposits. These differences reflect spatial and temporal variability in the operation of a range of transport processes, e.g. airfall, sea-ice and iceberg rafting, and post-depositional processes, e.g. bioturbation and secondary redeposition. The operation of these processes within the marine environment can potentially impart a temporal delay on tephra deposition and hamper the placement of the isochron, therefore, it is crucial to assess their influence. To aid this assessment a range of deposit types with common transport and depositional histories have been defined. Spatial patterns in the occurrence of these deposit types have been detected, the dominant controls at different sites explored and key regions of the North Atlantic with a greater likelihood for preserving isochronous deposits identified. Overall, these investigations have allowed a framework of isochronous marine cryptotephras to be defined for the last glacial period. The most widespread deposit is the rhyolitic phase of North Atlantic Ash Zone II, identified in 9 of the marine sequences and providing a direct tie-line to the Greenland ice-cores records. The framework is dominated by horizons with a basaltic composition, predominantly sourced from the Icelandic Grímsvötn volcanic system but horizons with Katla, Hekla, Kverkfjöll, Veidivötn and Vestmannaeyjar like compositions have also been isolated. Correlations to horizons in the Greenland ice-core tephra framework are being explored, however, this is a challenging process due to the large number of horizons with similar geochemical signatures in the records and the difference in temporal resolution and stratigraphic control between the ice and marine sequences.

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.

Studying ocean acidification in the Arctic Ocean

USGS Publications Warehouse

Robbins, Lisa

2012-01-01

The U.S. Geological Survey (USGS) partnership with the U.S. Coast Guard Ice Breaker Healey and its United Nations Convention Law of the Sea (UNCLOS) cruises has produced new synoptic data from samples collected in the Arctic Ocean and insights into the patterns and extent of ocean acidification. This framework of foundational geochemical information will help inform our understanding of potential risks to Arctic resources due to ocean acidification.

Final Report for the ZERT Project: Basic Science of Retention Issues, Risk Assessment & Measurement, Monitoring and Verification for Geologic Sequestration

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

Spangler, Lee; Cunningham, Alfred; Lageson, David

2011-03-31

ZERT has made major contributions to five main areas of sequestration science: improvement of computational tools; measurement and monitoring techniques to verify storage and track migration of CO{sub 2}; development of a comprehensive performance and risk assessment framework; fundamental geophysical, geochemical and hydrological investigations of CO{sub 2} storage; and investigate innovative, bio-based mitigation strategies.

Dolomitization of Quaternary reef limestones, Aitutaki, Cook Islands

USGS Publications Warehouse

Hein, J.R.; Gray, S.C.; Richmond, B.M.; White, L.D.

1992-01-01

The primary reef framework is considered to have been deposited during several highstands of sea level. Following partial to local recrystallization of the limestone, a signle episode of dolomitization occurred. Both tidal and thermal pumping drove large quantities of seawater through the porous rocks and perhaps maintained a wide mixing zone. However, the isotopic, geochemical and petrographic data do not clearly indicate the extent of seawater mixing. -from Authors

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

NASA Astrophysics Data System (ADS)

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

1981-03-01

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

Modeling biogechemical reactive transport in a fracture zone

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

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

2005-01-14

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

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

PubMed Central

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

2014-01-01

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

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

EPA Science Inventory

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

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

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

Mattson, Earl; Smith, Robert; Fujita, Yoshiko

2015-03-01

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

Biogeochemical Transformations in the History of the Ocean.

PubMed

Lenton, Timothy M; Daines, Stuart J

2017-01-03

The ocean has undergone several profound biogeochemical transformations in its 4-billion-year history, and these were an integral part of the coevolution of life and the planet. This review focuses on changes in ocean redox state as controlled by changes in biological activity, nutrient concentrations, and atmospheric O 2 . Motivated by disparate interpretations of available geochemical data, we aim to show how quantitative modeling-spanning microbial mats, shelf seas, and the open ocean-can help constrain past ocean biogeochemical redox states and show what caused transformations between them. We outline key controls on ocean redox structure and review pertinent proxies and their interpretation. We then apply this quantitative framework to three key questions: How did the origin of oxygenic photosynthesis transform ocean biogeochemistry? How did the Great Oxidation transform ocean biogeochemistry? And how was ocean biogeochemistry transformed in the Neoproterozoic-Paleozoic?

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.

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

PubMed Central

2011-01-01

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

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

PubMed Central

2011-01-01

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

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

NASA Astrophysics Data System (ADS)

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

2015-12-01

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

High-temperature life without photosynthesis as a model for Mars

NASA Technical Reports Server (NTRS)

Shock, E. L.

1997-01-01

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

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

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

NONE

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

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

NASA Astrophysics Data System (ADS)

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

2014-12-01

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

Characteristics and habitat of deep vs. shallow slow slip events

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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.

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

USGS Publications Warehouse

Day-Lewis, F. D.; Singha, K.; Binley, A.M.

2005-01-01

Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core-scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel-scale estimates. We refer to this loss of information as "correlation loss." Our approach upscales the core-scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical-resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near-surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications. Copyright 2005 by the American Geophysical Union.

New Comprehensive System to Construct Speleothem Fabrics Time Series

NASA Astrophysics Data System (ADS)

Frisia, S.; Borsato, A.

2014-12-01

Speleothem fabrics record processes that influence the way geochemical proxy data are encoded in speleothems, yet, there has been little advance in the use of fabrics as a complement to palaeo-proxy datasets since the fabric classification proposed by us in 2010. The systematic use of fabrics documentation in speleothem science has been limited by the absence of a comprehensive, numerical system that would allow constructing fabric time series comparable with the widely used geochemical time series. Documentation of speleothem fabrics is fundamental for a robust interpretation of speleothem time series where stable isotopes and trace elements are used as proxy, because fabrics highlight depositional as well as post-depositional processes whose understanding complements reconstructions based on geochemistry. Here we propose a logic system allowing transformation of microscope observations into numbers tied to acronyms that specify each fabric type and subtype. The rationale for ascribing progressive numbers to fabrics is based on the most up-to-date growth models. In this conceptual framework, the progression reflects hydrological conditions, bio-mediation and diagenesis. The lowest numbers are given to calcite fabrics formed at relatively constant drip rates: the columnar types (compact and open). Higher numbers are ascribed to columnar fabrics characterized by presence of impurities that cause elongation or lattice distortion (Elongated, Fascicular Optic and Radiaxial calcites). The sequence progresses with the dendritic fabrics, followed by micrite (M), which has been observed in association with microbial films. Microsparite (Ms) and mosaic calcite (Mc) have the highest numbers, being considered as diagenetic. Acronyms and subfixes are intended to become universally acknowledged. Thus, fabrics can be plotted vs. age to yield time series, where numbers are replaced by the acronyms. This will result in a visual representation of climate- or environmental-related parameters underpinning speleothem crystal growth. The Fabric log thus becomes a useful tool providing robustness to the geochemical data or test the overall utility of the speleothem record.

Design Support of an Above Cap-rock Early Detection Monitoring System using Simulated Leakage Scenarios at the FutureGen2.0 Site

DOE PAGES

Williams, Mark D.; USA, Richland Washington; Vermuel, Vince R.; ...

2014-12-31

The FutureGen 2.0 Project will design and build a first-of-its-kind, near-zero emissions coal-fueled power plant with carbon capture and storage (CCS). To assess storage site performance and meet the regulatory requirements of the Class VI Underground Injection Control (UIC) Program for CO 2 Geologic Sequestration, the FutureGen 2.0 project will implement a suite of monitoring technologies designed to evaluate CO 2 mass balance and detect any unforeseen loss in CO 2 containment. The monitoring program will include direct monitoring of the reservoir, and early-leak-detection monitoring directly above the primary confining zone. This preliminary modeling study described here focuses on hypotheticalmore » leakage scenarios into the first permeable unit above the primary confining zone (Ironton Sandstone) and is used to support assessment of early-leak detection capabilities. Future updates of the model will be used to assess potential impacts on the lowermost underground source of drinking water (Saint Peter Sandstone) for a range of theoretical leakage scenarios. This preliminary modeling evaluation considers both pressure response and geochemical signals in the overlying Ironton Sandstone. This model is independent of the FutureGen 2.0 reservoir model in that it does not simulate caprock discontinuities, faults, or failure scenarios. Instead this modeling effort is based on theoretical, volumetric-rate based leakage scenarios. The scenarios include leakage of 1% of the total injected CO 2 mass, but spread out over different time periods (20, 100, and 500 years) with each case yielding a different mass flux (i.e., smaller mass fluxes for longer duration leakage cases]. A brine leakage scenario using a volumetric leakage similar to the 20 year 1% CO 2 case was also considered. A framework for the comparison of the various cases was developed based on the exceedance of selected pressure and geochemical thresholds at different distances from the point of leakage and at different vertical positions within the Ironton Sandstone. These preliminary results, and results from an updated models that incorporate additional site-specific characterization data, support development/refinement of the monitoring system design.« less

The impact of pre-restoration land-use and disturbance on sediment structure, hydrology and the sediment geochemical environment in restored saltmarshes.

PubMed

Spencer, Kate L; Carr, Simon J; Diggens, Lucy M; Tempest, James A; Morris, Michelle A; Harvey, Gemma L

2017-06-01

Saltmarshes are being lost or degraded as a result of human activity resulting in loss of critical ecosystem services including the provision of wild species diversity, water quality regulation and flood regulation. To compensate, saltmarshes are being restored or re-created, usually driven by legislative requirements for increased habitat diversity, flood regulation and sustainable coastal defense. Yet, there is increasing evidence that restoration may not deliver anticipated ecosystem services; this is frequently attributed to poor drainage and sediment anoxia. However, physical sediment characteristics, hydrology and the sediment geochemical environment are rarely examined in restoration schemes, despite such factors being critical for plant succession. This study presents the novel integration of 3D-computed X-ray microtomography to quantify sediment structure and porosity, with water level and geochemical data to understand the impact of pre-restoration land use and disturbance on the structure and functioning of restored saltmarshes. The study combines a broad-scale investigation of physical sediment characteristics in nine de-embanked saltmarshes across SE England, with an intensive study at one site examining water levels, sediment structure and the sediment geochemical environment. De-embankment does not restore the hydrological regime, or the physical/chemical framework in the saltmarshes and evidence of disturbance includes a reduction in microporosity, pore connectivity and water storage capacity, a lack of connectivity between the sub-surface environment and overlying floodwaters, and impeded sub-surface water flow and drainage. This has significant consequences for the sediment geochemical environment. This disturbance is evident for at least two decades following restoration and is likely to be irreversible. It has important implications for plant establishment in particular, ecosystem services including flood regulation, nutrient cycling and wild species diversity and for future restoration design. Copyright © 2016 Elsevier B.V. All rights reserved.

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.

Reconstruction of paleostorm history using geochemical proxies in sediment cores from Eastern Lake, Florida

NASA Astrophysics Data System (ADS)

Das, O.; Wang, Y.; Donoghue, J. F.; Coor, J. L.; Kish, S.; Elsner, J.; Hu, X. B.; Niedoroda, A. W.; Ye, M.; Xu, Y.

2009-12-01

Analysis of geochemical proxies of coastal lake sediments provides a useful tool for reconstructing paleostorm history. Such paleostorm records can help constrain models that are used to predict future storm events. In this study, we collected two sediment cores (60 and 103 cm long, respectively) from the center of Eastern Lake located on the Gulf coast of NW Florida. These cores, which are mainly composed of organic-rich mud and organic-poor sand, were sub-sampled at 2-3mm intervals for analyses of their organic carbon and nitrogen concentrations as well as δ13C and δ15N isotopic signatures. Selected samples were submitted for radiocarbon dating in order to establish a chronological framework for the interpretation of the geochemical data. There are significant variations in δ13C, δ15N, C%, N% and C/N with depth. The δ13C and δ15N values vary from -21.8‰ to -26.7‰ and 2.6‰ to 5‰, respectively. The stable isotopic signatures of carbon and nitrogen indicate that the sources of organic matter in sediments include terrestrial C3 type vegetation, marine input from Gulf of Mexico and biological productivity within the lake, such as phytoplankton and zooplankton growing in the lacustrine environment. The δ13C and δ15N values exhibit significant negative excursions by 2‰ in a 30 cm thick sand layer, bounded by a rapid return to the base value. A positive shift in the δ15N record observed in the upper part of the cores likely reflects increased anthropogenic input of N such as sewage or septic tank effluents associated with recent development of areas around the lake for human habitation. Similarly, organic C% and N% range from 5.8 to 0.4 and 0.4 to 0.1, respectively. A prominent negative shift by 2σ relative to the baseline in C% and N% has been observed at approx. 55 to 58 cm depth, consisting of an organic-poor sand layer. This shift in C% and N% can be correlated with the negative shift in the δ13C and δ15N values, indicating a major storm event. Radiocarbon dates will allow establishment of the timing of the event and the interpretation of geochemical data in terms of changes in the lake environment.

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.

Interactive client side data visualization with d3.js

NASA Astrophysics Data System (ADS)

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

2015-12-01

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

Off-axis volcano-tectonic activity during continental rifting: Insights from the transversal Goba-Bonga lineament, Main Ethiopian Rift (East Africa)

NASA Astrophysics Data System (ADS)

Corti, Giacomo; Sani, Federico; Agostini, Samuele; Philippon, Melody; Sokoutis, Dimitrios; Willingshofer, Ernst

2018-03-01

The Main Ethiopian Rift, East Africa, is characterized by the presence of major, enigmatic structures which strike approximately orthogonal to the trend of the rift valley. These structures are marked by important deformation and magmatic activity in an off-axis position in the plateaus surrounding the rift. In this study, we present new structural data based on a remote and field analysis, complemented with analogue modelling experiments, and new geochemical analysis of volcanic rocks sampled in different portions of one of these transversal structures: the Goba-Bonga volcano-tectonic lineament (GBVL). This integrated analysis shows that the GBVL is associated with roughly E-W-trending prominent volcano-tectonic activity affecting the western plateau. Within the rift floor, the approximately E-W alignment of Awasa and Corbetti calderas likely represent expressions of the GBVL. Conversely, no tectonic or volcanic features of similar (E-W) orientation have been recognized on the eastern plateau. Analogue modelling suggests that the volcano-tectonic features of the GBVL have probably been controlled by the presence of a roughly E-W striking pre-existing discontinuity beneath the western plateau, which did not extend beneath the eastern plateau. Geochemical analysis supports this interpretation and indicates that, although magmas have the same sub-lithospheric mantle source, limited differences in magma evolution displayed by products found along the GBVL may be ascribed to the different tectonic framework to the west, to the east, and in the axial zone of the rift. These results support the importance of the heterogeneous nature of the lithosphere and the spatial variations of its structure in controlling the architecture of continental rifts and the distribution of the related volcano-tectonic activity.

Inverse Geochemical Reaction Path Modelling and the Impact of Climate Change on Hydrologic Structure in Snowmelt-Dominated Catchments in the Southwestern USA

NASA Astrophysics Data System (ADS)

Driscoll, J. M.; Meixner, T.; Molotch, N. P.; Sickman, J. O.; Williams, M. W.; McIntosh, J. C.; Brooks, P. D.

2011-12-01

Snowmelt from alpine catchments provides 70-80% of the American Southwest's water resources. Climate change threatens to alter the timing and duration of snowmelt in high elevation catchments, which may also impact the quantity and the quality of these water resources. Modelling of these systems provides a robust theoretical framework to process the information extracted from the sparse physical measurement available in these sites due to their remote locations. Mass-balance inverse geochemical models (via PHREEQC, developed by the USGS) were applied to two snowmelt-dominated catchments; Green Lake 4 (GL4) in the Rockies and Emerald Lake (EMD) in the Sierra Nevada. Both catchments primarily consist of granite and granodiorite with a similar bulk geochemistry. The inputs for the models were the initial (snowpack) and final (catchment output) hydrochemistry and a catchment-specific suite of mineral weathering reactions. Models were run for wet and dry snow years, for early and late time periods (defined hydrologically as 1/2 of the total volume for the year). Multiple model solutions were reduced to a representative suite of reactions by choosing the model solution with the fewest phases and least overall phase change. The dominant weathering reactions (those which contributed the most solutes) were plagioclase for GL4 and albite for EMD. Results for GL4 show overall more plagioclase weathering during the dry year (214.2g) than wet year (89.9g). Both wet and dry years show more weathering in the early time periods (63% and 56%, respectively). These results show that the snowpack and outlet are chemically more similar during wet years than dry years. A possible hypothesis to explain this difference is a change in contribution from subsurface storage; during the wet year the saturated catchment reduces contact with surface materials that would result in mineral weathering reactions by some combination of reduced infiltration and decreased subsurface transit time. By contrast, during the dry year infiltration and subsequent displacement of stored water that has had longer contact time with minerals and therefore has become more geochemically evolved to produce a greater difference between snowmelt and catchment outlet hydrochemistry. The results for EMD show little distinction between albite weathering for wet and dry years (55.9g and 66.0g, relatively). A hypothesis for this lack of difference in mineral phase changes may be due to less subsurface storage capacity in EMD relative to GL4. The spatial distribution of snowmelt has also been shown to influence the integrated watershed response, and future work includes using the Alpine Hydrochemical Model (AHM) to further investigate catchment response to these spatial data. The AHM will also provide further insight of surface-groundwater interactions through a more integrated model which includes hydrochemical, biological and physical processes to elucidate catchment response to changes in snowmelt dynamics.

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

NASA Astrophysics Data System (ADS)

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

2009-12-01

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

Comment on "Geochemistry of the Early Miocene volcanic succession of Northland, New Zealand, and implications for the evolution of subduction in the Southwest Pacific" by M.A. Booden, I.E.M. Smith, P.M. Black and J.L. Mauk

NASA Astrophysics Data System (ADS)

Schellart, W. P.

2012-01-01

In a recent paper Booden et al. (2011) present new geochemical and petrological data of Early Miocene volcanics from the Northland region (Northland volcanic belt) in New Zealand, and interpret these data to support a particular regional tectonic model. This tectonic model involves Early Miocene westward subduction of Cretaceous Pacific oceanic lithosphere below the Northland volcanic belt and the authors interpret the volcanic belt as a continental magmatic arc. Although the new data are not in disagreement with such a tectonic model, they provide more support for an alternative interpretation that involves a northeast-dipping subduction zone. Furthermore, geometric and plate kinematic data show that the west-dipping subduction model is unviable, geological and geophysical data contradict the model, while geodynamic arguments indicate that the model is implausible. Here it will be shown that a subduction model, involving a northeast-dipping southwestward retreating slab (made of the Late Cretaceous-Paleocene South Loyalty backarc basin lithosphere) that subsequently detaches, is in agreement with the local geology, geophysics and geochemistry, is geometrically, kinematically and geodynamically viable, and fits within the regional Southwest Pacific tectonic framework.

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.

Geochronology and geochemistry of tuff beds from the Shicaohe Formation of Shennongjia Group and tectonic evolution in the northern Yangtze Block, South China

NASA Astrophysics Data System (ADS)

Du, Qiuding; Wang, Zhengjiang; Wang, Jian; Deng, Qi; Yang, Fei

2016-03-01

Meso- to Neoproterozoic magmatic events are widespread in the Yangtze Block. The geochronology and tectonic significance of the Shennongjia Group in the Yangtze Block are still highly controversial. An integrated geochronology and geochemistry approach provides new insights into the geochronological framework, tectonic setting, magmatic events, and basin evolution of the northern Yangtze Block. Our new precise sensitive high-resolution ion microprobe U-Pb data indicate a deposition age of 1180 ± 15 Ma for the Shicaohe Formation subalkaline basaltic tuff that is geochemically similar to modern intracontinental rift volcanic rocks. The integration of available geochemical data together with our new U-Pb ages indicates the Shicaohe Formation subalkaline basaltic tuff formed ca. 1180 in a continental rift-related setting on a passive continental margin. The Shennongjia Group is topped by the Zhengjiaya Formation volcanic sequence, indicating arc-related igneous events at 1103 Ma. The transition of the late Mesoproterozoic tectonic regime from intracontinental extension to convergence occurred between ca. 1180 and 1103 Ma in the northern Yangtze Block. Tectonic evolution in the Neoproterozoic led to accretion along the northern margin of the Yangtze Block. These results provide geochronological evidence, which is of utmost importance for reconfiguration of the chronostratigraphic framework and for promoting research on Mesoproterozoic strata in China, thereby increasing understanding of magmatic events and basin evolutionary history in the northern Yangtze Block.

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

USGS Publications Warehouse

Washburn, Kathryn E.; Birdwell, Justin E.

2013-01-01

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

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

USGS Publications Warehouse

Parkhurst, David L.; Wissmeier, Laurin

2015-01-01

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

Combined time-series resistivity and geochemical tracer techniques to examine submarine groundwater discharge at Dor Beach, Israel

USGS Publications Warehouse

Swarzenski, P.W.; Burnett, W.C.; Greenwood, W.J.; Herut, B.; Peterson, R.; Dimova, N.; Shalem, Y.; Yechieli, Y.; Weinstein, Y.

2006-01-01

A high-resolution, stationary geophysical and geochemical survey was conducted at Dor Beach, Israel, to examine the shallow coastal hydrogeology and its control on the exchange of submarine groundwater with the shallow Mediterranean Sea. Time-series resistivity profiles using a new 56 electrode (112-m long) marine cable produced detailed profiles of the fresh water/salt water interface and the subtle response of this interface to tidal excursions and other forcing factors. Such information, when ground-truthed with representative pore water salinities and formation resistivity values, can provide unique information of the extent and rates of submarine groundwater discharge (SGD). Time-series 222Rn measurements of the adjacent coastal water column complemented these geophysical techniques and were modeled to yield integrated advective flow rates across the sediment/water interface, which at Dor Beach ranged from about 0 to 30 cm day-1 (mean = 7.1 cm d-1), depending on the tidal range. Such results suggest that the underlying hydrogeologic framework at Dor is favorable for substantial SGD. Extrapolating these SGD estimates across a 100-m wide coastal zone implies that the Rn-derived SGD rate would equal ???7.1 m3 d-1 per m of shoreline, and that the source of this discharging groundwater is a complex mixture of fresh groundwater derived from the upland Kurkar deposits, as well as locally recycled seawater. Copyright 2006 by the American Geophysical Union.

Current State of an Intelligent System to Aid in Tephra Layer Correlation

NASA Astrophysics Data System (ADS)

Hanson-Hedgecock, S.; Bursik, M.; Rogova, G.

2007-12-01

We are developing a computer based intelligent system to correlate tephra layers by using the lithologic, mineralogic, and geochemical characteristics of field samples, to aid geologists in interpreting eruption patterns of volcanic chains and fields. The intelligent system is used to define groups of tephra source vents by utilizing geochemical data, and to correlate tephra layers based on lithostratigraphic characteristics. Understanding the eruption history of a volcano from stratigraphic studies is important for forecasting future eruptive behavior and hazards. In volcanic chains and fields with a complex eruptive history and no central vent, determining the spatio- temporal eruption patterns is difficult. Sedimentologic and chemical variability, and sparse sampling often result in relatively large variances and imprecision in the dataset. Lithostratigraphic and geochemical interpretation also depends on ones' level of expertise and can be subjective. The processing of lithostratigraphic features is conducted by a hybrid classifier, composed of supervised artificial neural networks (ANNs) combined within the framework of the Dempster-Shafer theory of evidence. Since lithostratigraphic features vary with distance from source, hypothetical vent locations are determined by using expert domain knowledge and geostatistical methods. Geochemical data are processed by a suit of fuzzy k- means classifiers. Each fuzzy k-means classifier assigns observations to multiple clusters with various degrees, called membership coefficients. The assignment minimizes a function of the total distance between the centers of clusters and the individual geochemical data patterns weighed by the membership coefficients. Improved clustering results of geochemical data are achieved by the fusion of individual clustering results with an evidential combination method. Lithostratigraphic data from individual tephra beds of the North Mono eruption sequence are used to test the effectiveness of the intelligent system for tephra layer correlation. Geochemical data from tephra bedsets of the Mono and Inyo Craters, CA, are used to test the effectiveness of the intelligent system for eruption sequence correlation. The intelligent system aids correlation by showing matches and disparities between data patterns from different outcrops that may have been overlooked in initial interpretations. Initial results show that the lithostratigraphic classifier is able to accurately differentiate known layers 76% of the time. Output from the lithostratigraphic classifier can furthermore be plotted directly as isopleth maps that can aid in rapid recognition of tephra layers as well as determination of eruption characteristics, e.g. eruption volume, plume height, etc. The intelligent system produces a useful recognition result, while dealing with the uncertainty from sparse data and the imprecise description of layer characteristics.

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.

Using an intelligent system to aid in tephra layer correlation of the tephra beds of the Mono-Inyo Craters, California

NASA Astrophysics Data System (ADS)

Hanson-Hedgecock, S.; Bursik, M.; Rogova, G.

2008-12-01

We are developing an intelligent system to correlate tephra layers by using the lithologic and geochemical characteristics of field samples, to aid geologists in interpreting eruption patterns in volcanic fields. Understanding the eruption history of a volcanic field from stratigraphic studies is important for forecasting future eruptive behavior and hazards. The intelligent system is used to define groups of tephra source vents and to correlate tephra layers based on a combination of geochemical data and lithostratigraphic characteristics. The tephra beds of the Mono-Inyo Craters, California, are used to test the ability of the intelligent system for tephra layer correlation. The data processing is performed by a suite of both unsupervised and supervised classifiers, built and combined within the framework of the Dempster-Shafer theory of evidence. We have developed algorithms to calculate isopleth maps of thickness, lithic and pumice size that are used in the processing of the lithostratigraphic data. This spatial information is important in the determination of eruption patterns and is used by an evidential nearest neighbor classifier to correlate tephra layers. Integrating a better isopleth approximation function and expert knowledge about stratigraphic order of the tephra layers into the classifier improves the lithostratigraphic correlation from 56% to 87% of layers correctly identified. Geochemical data for defining groups of tephra sources are processed by a suit of fuzzy k-means classifiers. Improved clustering results of geochemical data are achieved by the fusion of individual clustering results with an evidential combination method. The intelligent system aids correlation by showing matches and disparities between data patterns from different outcrops that may have been overlooked. The intelligent system produces a useful recognition result, while dealing with the uncertainty from sparse data and the imprecise description of layer characteristics.

Baseline Geochemistry of Natural Occurring Methane and Saline Groundwater in an Area of Unconventional Shale Gas Development Through Time

NASA Astrophysics Data System (ADS)

Harkness, J.; Darrah, T.; Warner, N. R.; Whyte, C. J.; Moore, M. T.; Millot, R.; Kloppmann, W.; Jackson, R. B.; Vengosh, A.

2017-12-01

Naturally occurring methane is nearly ubiquitous in most sedimentary basins and delineating the effects of anthropogenic contamination sources from geogenic sources is a major challenge for evaluating the impact of unconventional shale gas development on water quality. This study employs a broadly integrated study of various geochemical techniques to investigate the geochemical variations of groundwater and surface water before, during, and after hydraulic fracturing.This approache combines inorganic geochemistry (major cations and anions), stable isotopes of select inorganic constituents including strontium (87Sr/86Sr), boron (δ11B), lithium (δ7Li), and carbon (δ13C-DIC), select hydrocarbon molecular (methane, ethane, propane, butane, and pentane) and isotopic tracers (δ13C-CH4, δ13C-C2H6), tritium (3H), and noble gas elemental and isotopic composition (He, Ne, Ar) to apportion natural and anthropogenic sources of natural gas and salt contaminants both before and after drilling. Methane above 1 ccSTP/L in groundwater samples awas strongly associated with elevated salinity (chloride >50 mg/L).The geochemical and isotopic analysis indicate saline groundwater originated via naturally occurring processes, presumably from the migration of deeper methane-rich brines that have interacted extensively with coal lithologies. The chemistry and gas compostion of both saline and fresh groundwater wells did not change following the installation of nearby shale-gas wells.The results of this study emphasize the value of baseline characterization of water quality in areas of fossil fuel exploration. Overall this study presents a comprehensive geochemical framework that can be used as a template for assessing the sources of elevated hydrocarbons and salts to water resources in areas potentially impacted by oil and gas development.

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

USGS Publications Warehouse

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

2012-01-01

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

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

NASA Astrophysics Data System (ADS)

Asavin, A. M.

2001-12-01

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

Using 3D geological modelling and geochemical mixing models to characterise alluvial aquifer recharge sources in the upper Condamine River catchment, Queensland, Australia.

PubMed

Martinez, Jorge L; Raiber, Matthias; Cendón, Dioni I

2017-01-01

The influence of mountain front recharge on the water balance of alluvial valley aquifers located in upland catchments of the Condamine River basin in Queensland, Australia, is investigated through the development of an integrated hydrogeological framework. A combination of three-dimensional (3D) geological modelling, hydraulic gradient maps, multivariate statistical analyses and hydrochemical mixing calculations is proposed for the identification of hydrochemical end-members and quantification of the relative contributions of each end-member to alluvial aquifer recharge. The recognised end-members correspond to diffuse recharge and lateral groundwater inflows from three hydrostratigraphic units directly connected to the alluvial aquifer. This approach allows mapping zones of potential inter-aquifer connectivity and areas of groundwater mixing between underlying units and the alluvium. Mixing calculations using samples collected under baseflow conditions reveal that lateral contribution from a regional volcanic aquifer system represents the majority (41%) of inflows to the alluvial aquifer. Diffuse recharge contribution (35%) and inflow from two sedimentary bedrock hydrostratigraphic units (collectively 24%) comprise the remainder of major recharge sources. A detailed geochemical assessment of alluvial groundwater evolution along a selected flowpath of a representative subcatchment of the Condamine River basin confirms mixing as a key process responsible for observed spatial variations in hydrochemistry. Dissolution of basalt-related minerals and dolomite, CO 2 uptake, ion-exchange, precipitation of clay minerals, and evapotranspiration further contribute to the hydrochemical evolution of groundwater in the upland alluvial aquifer. This study highlights the benefits of undertaking an integrated approach that combines multiple independent lines of evidence. The proposed methods can be applied to investigate processes associated with inter-aquifer mixing, including groundwater contamination resulting from depressurisation of underlying geological units hydraulically connected to the shallower water reservoirs. Copyright © 2016 Elsevier B.V. All rights reserved.

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.

Experimental investigations and geochemical modelling of site-specific fluid-fluid and fluid-rock interactions in underground storage of CO2/H2/CH4 mixtures: the H2STORE project

NASA Astrophysics Data System (ADS)

De Lucia, Marco; Pilz, Peter

2015-04-01

Underground gas storage is increasingly regarded as a technically viable option for meeting the energy demand and environmental targets of many industrialized countries. Besides the long-term CO2 sequestration, energy can be chemically stored in form of CO2/CH4/H2 mixtures, for example resulting from excess wind energy. A precise estimation of the impact of such gas mixtures on the mineralogical, geochemical and petrophysical properties of specific reservoirs and caprocks is crucial for site selection and optimization of storage depth. Underground gas storage is increasingly regarded as a technically viable option for meeting environmental targets and the energy demand through storage in form of H2 or CH4, i.e. resulting from excess wind energy. Gas storage in salt caverns is nowadays a mature technology; in regions where favorable geologic structures such as salt diapires are not available, however, gas storage can only be implemented in porous media such as depleted gas and oil reservoirs or suitable saline aquifers. In such settings, a significant amount of in-situ gas components such as CO2, CH4 (and N2) will always be present, making the CO2/CH4/H2 system of particular interest. A precise estimation of the impact of their gas mixtures on the mineralogical, geochemical and petrophysical properties of specific reservoirs and caprocks is therefore crucial for site selection and optimization of storage depth. In the framework of the collaborative research project H2STORE, the feasibility of industrial-scale gas storage in porous media in several potential siliciclastic depleted gas and oil reservoirs or suitable saline aquifers is being investigated by means of experiments and modelling on actual core materials from the evaluated sites. Among them are the Altmark depleted gas reservoir in Saxony-Anhalt and the Ketzin pilot site for CO2 storage in Brandenburg (Germany). Further sites are located in the Molasse basin in South Germany and Austria. In particular, two work packages hosted at the German Research Centre for Geosciences (GFZ) focus on the fluid-fluid and fluid-rock interactions triggered by CO2, H2 and their mixtures. Laboratory experiments expose core samples to hydrogen and CO2/hydrogen mixtures under site-specific conditions (temperatures up to 200 °C and pressure up to 300 bar). The resulting qualitative and, whereas possible, quantitative data are expected to ameliorate the precision of predictive geochemical and reactive transport modelling, which is also performed within the project. The combination of experiments, chemical and mineralogical analyses and models is needed to improve the knowledge about: (1) solubility model and mixing rule for multicomponent gas mixtures in high saline formation fluids: no data are namely available in literature for H2-charged gas mixtures in the conditions expected in the potential sites; (2) chemical reactivity of different mineral assemblages and formation fluids in a broad spectrum of P-T conditions and composition of the stored gas mixtures; (3) thermodynamics and kinetics of relevant reactions involving mineral dissolution or precipitation. The resulting amelioration of site characterization and the overall enhancement in understanding the potential processes will benefit the operational reliability, the ecological tolerance, and the economic efficiency of future energy storing plants, crucial aspects for public acceptance and for industrial investors.

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.

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

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

NASA Astrophysics Data System (ADS)

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

2015-04-01

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

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

NASA Astrophysics Data System (ADS)

Shamberger, Patrick J.; Garcia, Michael O.

2007-02-01

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

Predicting Biological Information Flow in a Model Oxygen Minimum Zone

NASA Astrophysics Data System (ADS)

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

2016-02-01

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

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

PubMed

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

2016-12-01

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

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

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

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

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

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

NASA Astrophysics Data System (ADS)

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

2011-12-01

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

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

DOE PAGES

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

2016-09-16

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

Gaining the necessary geologic, hydrologic, and geochemical understanding for additional brackish groundwater development, coastal San Diego, California, USA

USGS Publications Warehouse

Danskin, Wesley R.

2012-01-01

Local water agencies and the United States Geological Survey are using a combination of techniques to better understand the scant freshwater resources and the much more abundant brackish resources in coastal San Diego, California, USA. Techniques include installation of multiple-depth monitoring well sites; geologic and paleontological analysis of drill cuttings; geophysical logging to identify formations and possible seawater intrusion; sampling of pore-water obtained from cores; analysis of chemical constituents including trace elements and isotopes; and use of scoping models including a three-dimensional geologic framework model, rainfall-runoff model, regional groundwater flow model, and coastal density-dependent groundwater flow model. Results show that most fresh groundwater was recharged during the last glacial period and that the coastal aquifer has had recurring intrusions of fresh and saline water. These intrusions disguise the source, flowpaths, and history of ground water near the coast. The flow system includes a freshwater lens resting on brackish water; a 100-meter-thick flowtube of freshwater discharging under brackish estuarine water and above highly saline water; and broad areas of fine-grained coastal sediment filled with fairly uniform brackish water. Stable isotopes of hydrogen and oxygen indicate the recharged water flows through many kilometers of fractured crystalline rock before entering the narrow coastal aquifer.

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

USGS Publications Warehouse

Amos, Richard T.; Mayer, K. Ulrich

2006-01-01

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

Origin of the cataclysmic Late Heavy Bombardment period of the terrestrial planets.

PubMed

Gomes, R; Levison, H F; Tsiganis, K; Morbidelli, A

2005-05-26

The petrology record on the Moon suggests that a cataclysmic spike in the cratering rate occurred approximately 700 million years after the planets formed; this event is known as the Late Heavy Bombardment (LHB). Planetary formation theories cannot naturally account for an intense period of planetesimal bombardment so late in Solar System history. Several models have been proposed to explain a late impact spike, but none of them has been set within a self-consistent framework of Solar System evolution. Here we propose that the LHB was triggered by the rapid migration of the giant planets, which occurred after a long quiescent period. During this burst of migration, the planetesimal disk outside the orbits of the planets was destabilized, causing a sudden massive delivery of planetesimals to the inner Solar System. The asteroid belt was also strongly perturbed, with these objects supplying a significant fraction of the LHB impactors in accordance with recent geochemical evidence. Our model not only naturally explains the LHB, but also reproduces the observational constraints of the outer Solar System.

Empirical kinetics and their role in elucidating the utility of transition-state theory to mineral–water reactions. A comment upon, ''Evidence and Potential Implications of Exponential Tails to Concentration Versus Time Plots for the Batch Dissolution of Calcite'' by V. W. Truesdale

DOE PAGES

Icenhower, Jonathan P.

2015-06-23

Transition-state theory (TST) is a successful theory for understanding many different types of reactions, but its application to mineral-water systems has not been successful, especially as the system approaches saturation with respect to a rate-limiting phase. A number of investigators have proposed alternate frameworks for using the kinetic rate data to construct models of dissolution, including Truesdale (Aquat Geochem, 2015; this issue). This alternate approach has been resisted, in spite of self-evident discrepancies between TST expectations and the data. The failure of TST under certain circumstances is a result of the presence of metastable intermediaries or reaction layers that formmore » on the surface of reacting solids, and these phenomena are not anticipated by the current theory. Furthermore, alternate approaches, such as the shrinking object model advocated by Truesdale, represent a potentially important avenue for advancing the science of dissolution kinetics.« less

MetPetDB: A database for metamorphic geochemistry

NASA Astrophysics Data System (ADS)

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

2009-12-01

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

Community-Based Development of Standards for Geochemical and Geochronological Data

NASA Astrophysics Data System (ADS)

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

2007-12-01

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

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

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

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

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

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.

The developing framework of marine ecotoxicology: Pollutants as a variable in marine ecosystems?

USGS Publications Warehouse

Luoma, Samuel N.

1996-01-01

Marine ecosystems include a subset in which at least some interrelated geochemical, biochemical, physiological, population and community characteristics are changed by pollutants. Moderate contamination is relatively widespread in coastal and estuarine ecosystems, so the subset of ecosystems with at least some processes affected could be relatively large. Pollutant influences have changed and will probably continue to change on time scales of decades. Biological exposures and dose in such ecosystems are species-specific and determined by how the species is exposed to different environmental media and the geochemistry of individual pollutants within those media. Bioaccumulation models offer significant promise for interpreting such exposures. Biological responses to pollutants need to be more directly linked to exposure and dose. At the level of the individual this might be improved by better understanding relationships between tissue concentrations of pollutants and responses to pollutants. Multi-discipline field and laboratory studies combined with advanced understanding of some basic processes have reduced the ambiguities in interpreting a few physiological/organismic responses to pollutants in nature. Recognition of pollutant-induced patterns in population responses could lead to similar advances. A rational framework for ecotoxicology is developing, but its further advance is dependent upon better integration of ecotoxicology with basic marine ecology and biology.

Stratigraphic framework of Cambrian and Ordovician rocks in the central Appalachian basin from Medina County, Ohio, through southwestern and south-central Pennsylvania to Hampshire County, West Virginia: Chapter E.2.2 in Coal and petroleum resources in the Appalachian basin: distribution, geologic framework, and geochemical character

USGS Publications Warehouse

Ryder, Robert T.; Harris, Anita G.; Repetski, John E.; Crangle, Robert D.; Ruppert, Leslie F.; Ryder, Robert T.

2014-01-01

This chapter is a re-release of U.S. Geological Survey Bulletin 1839-K, of the same title, by Ryder and others (1992; online version 2.0 revised and digitized by Robert D. Crangle, Jr., 2003). It consists of one file of the report text as it appeared in USGS Bulletin 1839-K and a second file containing the cross section, figures 1 and 2, and tables 1 and 2 on one oversized sheet; the second file was digitized in 2003 as version 2.0 and also includes the gamma-ray well log traces.

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.

Modelling coupled microbial processes in the subsurface: Model development, verification, evaluation and application

NASA Astrophysics Data System (ADS)

Masum, Shakil A.; Thomas, Hywel R.

2018-06-01

To study subsurface microbial processes, a coupled model which has been developed within a Thermal-Hydraulic-Chemical-Mechanical (THCM) framework is presented. The work presented here, focuses on microbial transport, growth and decay mechanisms under the influence of multiphase flow and bio-geochemical reactions. In this paper, theoretical formulations and numerical implementations of the microbial model are presented. The model has been verified and also evaluated against relevant experimental results. Simulated results show that the microbial processes have been accurately implemented and their impacts on porous media properties can be predicted either qualitatively or quantitatively or both. The model has been applied to investigate biofilm growth in a sandstone core that is subjected to a two-phase flow and variable pH conditions. The results indicate that biofilm growth (if not limited by substrates) in a multiphase system largely depends on the hydraulic properties of the medium. When the change in porewater pH which occurred due to dissolution of carbon dioxide gas is considered, growth processes are affected. For the given parameter regime, it has been shown that the net biofilm growth is favoured by higher pH; whilst the processes are considerably retarded at lower pH values. The capabilities of the model to predict microbial respiration in a fully coupled multiphase flow condition and microbial fermentation leading to production of a gas phase are also demonstrated.

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.

Do weirs affect the physical and geochemical mobility of toxic metals in mining-impacted floodplain sediments?

NASA Astrophysics Data System (ADS)

Bulcock, Amelia; Coleman, Alexandra; Whitfield, Elizabeth; Andres Lopez-Tarazon, Jose; Byrne, Patrick; Whitfield, Greg

2015-04-01

Weirs are common river structures designed to modify river channel hydraulics and hydrology for purposes of navigation, flood defence, irrigation and hydrometry. By design, weirs constrain natural flow processes and affect sediment flux and river channel forms leading to homogenous river habitats and reduced biodiversity. The recent movement towards catchment-wide river restoration, driven by the EU Water Framework Directive, has recognised weirs as a barrier to good ecological status. However, the removal of weirs to achieve more 'natural' river channels and flow processes is inevitably followed by a period of adjustment to the new flow regime and sediment flux. This period of adjustment can have knock-on effects that may increase flood risk, sedimentation and erosion until the river reaches a state of geomorphological equilibrium. Many catchments in the UK contain a legacy of toxic metals in floodplain sediments due to historic metal mining activities. The consequences of weir removal in these catchments may be to introduce 'stored' mine wastes into the river system with severe implications for water quality and biodiversity. The purpose of this study is to investigate the potential impact of a weir on the physical and geochemical mobilisation of mine wastes in the formerly mined River Twymyn catchment, Wales. Our initial investigations have shown floodplain and riverbed sediments to be grossly contaminated (up to 15,500 mg/kg Pb) compared to soil from a pre-mining Holocene terrace (180 mg/kg Pb). Geomorphological investigations also suggest that weir removal will re-establish more dynamic river channel processes resulting in lateral migration of the channel and erosion of contaminated floodplain sediments. These data will be used as a baseline for more detailed investigations of the potential impact of weirs on the physical and geochemical mobilisation of contaminated sediments. We have two specific objectives. (1) Geomorphological assessments will use unmanned aerial vehicle (UAV) photographic surveys, historical aerial photographs, ground-based topographic surveys, surface and subsurface particle size determination, bed stability and sediment entrainment assessment, together with discharge and sediment (both suspended and bedload) monitoring to establish the effect of the weir on patterns of sediment flux and the physical transport of metal contaminants. 2D and 1D models (IBER, HEC-RAS) of the weir-affected reach will investigate sediment and metal flux following weir removal. (2) The physicochemical speciation and geochemical stability of contaminated floodplain sediments will be characterised using bulk chemistry, mineralogical (XRD, SEM) and speciation methods (sequential extractions, electron microprobe analysis).

Hydrogeothermal Convective Circulation Model for the Formation of the Chicxulub Ring of Cenotes in the Yucatan Peninsula, Mexico.

NASA Astrophysics Data System (ADS)

Monroy-Rios, E.; Beddows, P. A.

2015-12-01

Despite being deeply buried, the topography and geophysical characteristics of the multi-ring Chicxulub impact structure are reflected on the now subaerial Yucatan Peninsula with aligned arcs of cenotes (sinkholes), forming the "Ring of Cenotes". A pending question is the determination of the geological, geochemical, structural features and associated processes that have led to void development, and the upwards propagation of the voids, cross cutting over 1000 m of super-deposited carbonate sequences. Drawing from the published literature on drill core and geophysical surveys undertaken by Pemex, UNAM, and IODP/ICDP, numerical modeling, and general carbonate platform hydrothermal reactive transport models, we provide a conceptual model for the genesis of the Ring of Cenotes. In horizontally bedded carbonate platforms, geothermal gradients will drive convective flow, with strong vertical components specifically in the platform center. In the Yucatan Platform, a high occurrence of anhydrite and dolomite at depth evokes early burial dolomitization and anhydritization, sourcing Mg from seawater. The Chicxulub impact near the center of the platform produced a low permeability and high thermal conductivity melt rock that arguably extends to the basement rock at 3.5 km below surface. Heat of impact enforced the pre-existing geothermal circulation pattern, and even with depletion of the heat of impact, the high thermal conductivity of the crystalline melt would lead to enhanced geothermal gradients in the center of the platform. The cenotes overlying the crater are deep (150+ m) vertical shafts with most (but not all) breaching the surface. The pit geomorphology suggests a bottom-up formation. Excess Si in the shallow groundwater points to a convective circulation with strong vertical components geochemically linking the granodioritic basement rock to the surface. Water temperature and conductivity profiles support ongoing vertical flux in some deep pit cenotes. Within this framework, we argue for the formation of the Ring of Cenotes by hydrogeothermal convective circulation in the post-impact carbonate sequences, leading to spatially focused dissolution at depth, with voids initiated along the crater edge effectively propagating upwards, often breaching the surface.

The origin and nature of thermal evolution during Granite emplacement and differentiation and its influence on upper crustal dynamics.

NASA Astrophysics Data System (ADS)

Buchwaldt, R.; Toulkeridis, T.; Todt, W.

2014-12-01

Structural geological, geochemical and geochronological data were compiled with the purpose to exercise models for the construction of upper crustal batholith. Models for pulsed intrusion of small magma batches over long timescales versus transfer of larger magma bodies on a shorter time scales are able to predict a different thermal, metamorphic, and rheological state of the crust. For this purpose we have applied the chronostratigraphic framework for magma differentiation on three granite complexes namely the St. Francois Mountain granite pluton (Precambrian), the Galway granite (Cambrian), and the Sithonia Plutonic Complex (Eocene). These plutons have similar sizes and range in composition from quartz diorites through granodiorites and granites to alkali granites, indicating multiple intrusive episodes. Thermobarometric calculations imply an upper crustal emplacement. Geochemical, isotopic and petrological data indicate a variety of pulses from each pluton allowing to be related through their liquid line of decent, which is supported by fractional crystallization of predominantly plagioclase, K-feldspar, biotite, hornblende and some minor accessory mineral phases, magma mingling and mixing as well as crustal contamination. To obtain the temporal relationship we carried out high-precision CA-TIMS zircon geochronology on selected samples along the liquid line of decent. The obtained data indicate a wide range of rates: such as different pulses evolved on timescales of about only 10-30ka, although, the construction time of the different complexes ranges from millions of years with prolonged tectonically inactive phases to relatively short lived time ranges of about ~300 ka. For a better understanding how these new data were used and evaluated in order to reconstruct constraints on the dynamics of the magmatic plumbing system, we integrated the short-lived, elevated heat production, due to latent heat of crystallization, into a 2D numerical model of the thermal evolution of segments of continental crust. Our model indicates that during the stage of enhanced fractional crystallization, the crustal viscosity decreases by several orders of magnitude, playing hereby a fundamental role in the thermal, magmatic, and tectonic evolution of the studied areas and most probably in similar regions too.

Quantifying Reactive Transport Processes Governing Arsenic Mobility after Injection of Reactive Organic Carbon into a Bengal Delta Aquifer.

PubMed

Rawson, Joey; Siade, Adam; Sun, Jing; Neidhardt, Harald; Berg, Michael; Prommer, Henning

2017-08-01

Over the last few decades, significant progress has been made to characterize the extent, severity, and underlying geochemical processes of groundwater arsenic (As) pollution in S/SE Asia. However, comparably little effort has been made to merge the findings into frameworks that allow for a process-based quantitative analysis of observed As behavior and for predictions of its long-term fate. This study developed field-scale numerical modeling approaches to represent the hydrochemical processes associated with an in situ field injection of reactive organic carbon, including the reductive dissolution and transformation of ferric iron (Fe) oxides and the concomitant release of sorbed As. We employed data from a sucrose injection experiment in the Bengal Delta Plain to guide our model development and to constrain the model parametrization. Our modeling results illustrate that the temporary pH decrease associated with the sucrose transformation and mineralization caused pronounced, temporary shifts in the As partitioning between aqueous and sorbed phases. The results also suggest that while the reductive dissolution of Fe(III) oxides reduced the number of sorption sites, a significant fraction of the released As was rapidly scavenged through coprecipitation with neo-formed magnetite. These secondary reactions can explain the disparity between the observed Fe and As behavior.

Synthesis of benthic flux components in the Patos Lagoon coastal zone, Rio Grande do Sul, Brazil

NASA Astrophysics Data System (ADS)

King, J. N.

2012-12-01

The primary objective of this work is to synthesize components of benthic flux in the Patos Lagoon coastal zone, Rio Grande do Sul, Brazil. Specifically, the component of benthic discharge flux forced by the terrestrial hydraulic gradient is 0.8 m3 d-1; components of benthic discharge and recharge flux associated with the groundwater tidal prism are both 2.1 m3 d-1; components of benthic discharge and recharge flux forced by surface-gravity wave setup are both 6.3 m3 d-1; the component of benthic discharge flux that transports radium-228 is 350 m3 d-1; and components of benthic discharge and recharge flux forced by surface-gravity waves propagating over a porous medium are both 1400 m3 d-1. (All models are normalized per meter shoreline.) Benthic flux is a function of components forced by individual mechanisms and nonlinear interactions that exist between components. Constructive and destructive interference may enhance or diminish the contribution of benthic flux components. It may not be possible to model benthic flux by summing component magnitudes. Geochemical tracer techniques may not accurately model benthic discharge flux or submarine groundwater discharge (SGD). A conceptual model provides a framework on which to quantitatively characterize benthic discharge flux and SGD with a multifaceted approach.

Synthesis of benthic flux components in the Patos Lagooncoastal zone, Rio Grande do Sul, Brazil

USGS Publications Warehouse

King, Jeffrey N.

2012-01-01

The primary objective of this work is to synthesize components of benthic flux in the Patos Lagoon coastal zone, Rio Grande do Sul, Brazil. Specifically, the component of benthic discharge flux forced by the terrestrial hydraulic gradient is 0.8 m3 d-1; components of benthic discharge and recharge flux associated with the groundwater tidal prism are both 2.1 m3 d-1; components of benthic discharge and recharge flux forced by surface-gravity wave setup are both 6.3 m3 d-1; the component of benthic discharge flux that transports radium-228 is 350 m3 d-1; and components of benthic discharge and recharge flux forced by surface-gravity waves propagating over a porous medium are both 1400 m3 d-1. (All models are normalized per meter shoreline.) Benthic flux is a function of components forced by individual mechanisms and nonlinear interactions that exist between components. Constructive and destructive interference may enhance or diminish the contribution of benthic flux components. It may not be possible to model benthic flux by summing component magnitudes. Geochemical tracer techniques may not accurately model benthic discharge flux or submarine groundwater discharge (SGD). A conceptual model provides a framework on which to quantitatively characterize benthic discharge flux and SGD with a multifaceted approach.

An array processing system for lunar geochemical and geophysical data

NASA Technical Reports Server (NTRS)

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

1977-01-01

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

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

USGS Publications Warehouse

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

2003-01-01

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

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

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.