Modeling Background Radiation in our Environment Using Geochemical Data

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

Malchow, Russell L.; Marsac, Kara; Burnley, Pamela

2015-02-01

Radiation occurs naturally in bedrock and soil. Gamma rays are released from the decay of the radioactive isotopes K, U, and Th. Gamma rays observed at the surface come from the first 30 cm of rock and soil. The energy of gamma rays is specific to each isotope, allowing identification. For this research, data was collected from national databases, private companies, scientific literature, and field work. Data points were then evaluated for self-consistency. A model was created by converting concentrations of U, K, and Th for each rock and soil unit into a ground exposure rate using the following equation:more » D=1.32 K+ 0.548 U+ 0.272 Th. The first objective of this research was to compare the original Aerial Measurement System gamma ray survey to results produced by the model. The second objective was to improve the method and learn the constraints of the model. Future work will include sample data analysis from field work with a goal of improving the geochemical model.« less

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.

COTHERM: Modelling fluid-rock interactions in Icelandic geothermal systems

NASA Astrophysics Data System (ADS)

Thien, Bruno; Kosakowski, Georg; Kulik, Dmitrii

2014-05-01

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

Microbial Response in Peat Overlying Kimberlite Pipes in The Attawapiskat Area, Northern Ontario

NASA Astrophysics Data System (ADS)

Donkervoort, L. J.; Southam, G.

2009-05-01

Exploration for ore deposits occurring under thick, post-mineralized cover requires innovative methods and instrumentation [1]. Buried kimberlite pipes 'produce' geochemical conditions such as increased pH and decreased Eh in overlying peat [2] that intuitively select for bacterial populations that are best able to grow and, which in turn affect the geochemistry producing a linked signal. A microbiological study of peat was conducted over the Zulu kimberlite in the Attawapiskat area of the James Bay Lowlands to determine if the type of underlying rock influences the diversity and populations of microorganisms living in the overlying peat. Peat was sampled along an 800 m transect across the Zulu kimberlite, including samples underlain by limestone. Microbial populations and carbon source utilization patterns of peat samples were compared between the two underlying rock types. Results demonstrate an inverse relationship of increased anaerobic populations and lower biodiversity directly above the kimberlite pipe. These results support a reduced 'column' consistent with the model presented by Hamilton [3]. The combination of traditional bacterial enumeration and community- level profiling represents a cost-effective and efficient exploration technique that can serve to compliment both geophysical and geochemical surveys. [1] Goldberg (1998) J. Geochem. Explor. 61, 191-202 [2] Hattori and Hamilton (2008) Appl. Geochem. 23, 3767-3782 [3] Hamilton (1998) J. Geochem. Explor. 63, 155-172

Spatial variation of volcanic rock geochemistry in the Virunga Volcanic Province: Statistical analysis of an integrated database

NASA Astrophysics Data System (ADS)

Barette, Florian; Poppe, Sam; Smets, Benoît; Benbakkar, Mhammed; Kervyn, Matthieu

2017-10-01

We present an integrated, spatially-explicit database of existing geochemical major-element analyses available from (post-) colonial scientific reports, PhD Theses and international publications for the Virunga Volcanic Province, located in the western branch of the East African Rift System. This volcanic province is characterised by alkaline volcanism, including silica-undersaturated, alkaline and potassic lavas. The database contains a total of 908 geochemical analyses of eruptive rocks for the entire volcanic province with a localisation for most samples. A preliminary analysis of the overall consistency of the database, using statistical techniques on sets of geochemical analyses with contrasted analytical methods or dates, demonstrates that the database is consistent. We applied a principal component analysis and cluster analysis on whole-rock major element compositions included in the database to study the spatial variation of the chemical composition of eruptive products in the Virunga Volcanic Province. These statistical analyses identify spatially distributed clusters of eruptive products. The known geochemical contrasts are highlighted by the spatial analysis, such as the unique geochemical signature of Nyiragongo lavas compared to other Virunga lavas, the geochemical heterogeneity of the Bulengo area, and the trachyte flows of Karisimbi volcano. Most importantly, we identified separate clusters of eruptive products which originate from primitive magmatic sources. These lavas of primitive composition are preferentially located along NE-SW inherited rift structures, often at distance from the central Virunga volcanoes. Our results illustrate the relevance of a spatial analysis on integrated geochemical data for a volcanic province, as a complement to classical petrological investigations. This approach indeed helps to characterise geochemical variations within a complex of magmatic systems and to identify specific petrologic and geochemical investigations that should be tackled within a study area.

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.

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.

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.

Continuum-based DFN-consistent simulations of oxygen ingress in fractured crystalline rocks

NASA Astrophysics Data System (ADS)

Trinchero, P.; Puigdomenech, I.; Molinero, J.; Ebrahimi, H.; Gylling, B.; Svensson, U.; Bosbach, D.; Deissmann, G.

2016-12-01

The potential transient infiltration of oxygenated glacial meltwater into initially anoxic and reducing fractured crystalline rocks during glaciation events is an issue of concern for some of the prospected deep geological repositories for spent nuclear fuel. Here, this problem is assessed using reactive transport calculations. First, a novel parameterisation procedure is presented, where flow, transport and geochemical parameters (i.e. hydraulic conductivity, effective/kinetic porosity, and mineral specific surface and abundance) are defined on a finite volume numerical grid based on the (spatially varying) properties of an underlying Discrete Fracture Network (DFN). Second, using this approach, a realistic reactive transport model of Forsmark, i.e. the selected site for the proposed Swedish spent nuclear fuel repository, is implemented. The model consists of more than 70 million geochemical transport degrees of freedom and simulates the ingress of oxygen-rich water from the recharge area of the domain and its depletion due to reactions with the Fe(II) mineral chlorite. Third, the calculations are solved in the supercomputer JUQUEEN of the Jülich Supercomputing Centre. The results of the simulations show that oxygen infiltrates relatively quickly along fractures and deformation zones until a steady state profile is reached, where geochemical reactions counterbalance advective transport processes. Interestingly, most of the iron-bearing minerals are consumed in the highly conductive zones, where larger mineral surfaces are available for reactions. An analysis based on mineral mass balance shows that the considered rock medium has enough capacity to buffer oxygen infiltration for a long period of time (i.e. some thousand years).

Parameter Sensitivity and Laboratory Benchmarking of a Biogeochemical Process Model for Enhanced Anaerobic Dechlorination

NASA Astrophysics Data System (ADS)

Kouznetsova, I.; Gerhard, J. I.; Mao, X.; Barry, D. A.; Robinson, C.; Brovelli, A.; Harkness, M.; Fisher, A.; Mack, E. E.; Payne, J. A.; Dworatzek, S.; Roberts, J.

2008-12-01

A detailed model to simulate trichloroethene (TCE) dechlorination in anaerobic groundwater systems has been developed and implemented through PHAST, a robust and flexible geochemical modeling platform. The approach is comprehensive but retains flexibility such that models of varying complexity can be used to simulate TCE biodegradation in the vicinity of nonaqueous phase liquid (NAPL) source zones. The complete model considers a full suite of biological (e.g., dechlorination, fermentation, sulfate and iron reduction, electron donor competition, toxic inhibition, pH inhibition), physical (e.g., flow and mass transfer) and geochemical processes (e.g., pH modulation, gas formation, mineral interactions). Example simulations with the model demonstrated that the feedback between biological, physical, and geochemical processes is critical. Successful simulation of a thirty-two-month column experiment with site soil, complex groundwater chemistry, and exhibiting both anaerobic dechlorination and endogenous respiration, provided confidence in the modeling approach. A comprehensive suite of batch simulations was then conducted to estimate the sensitivity of predicted TCE degradation to the 36 model input parameters. A local sensitivity analysis was first employed to rank the importance of parameters, revealing that 5 parameters consistently dominated model predictions across a range of performance metrics. A global sensitivity analysis was then performed to evaluate the influence of a variety of full parameter data sets available in the literature. The modeling study was performed as part of the SABRE (Source Area BioREmediation) project, a public/private consortium whose charter is to determine if enhanced anaerobic bioremediation can result in effective and quantifiable treatment of chlorinated solvent DNAPL source areas. The modelling conducted has provided valuable insight into the complex interactions between processes in the evolving biogeochemical systems, particularly at the laboratory scale.

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

Negus-De Wys, J.; Dixon, J. M.; Evans, M. A.

This document consists of the following papers: inorganic geochemistry studies of the Eastern Kentucky Gas Field; lithology studies of upper Devonian well cuttings in the Eastern Kentucky Gas Field; possible effects of plate tectonics on the Appalachian Devonian black shale production in eastern Kentucky; preliminary depositional model for upper Devonian Huron age organic black shale in the Eastern Kentucky Gas Field; the anatomy of a large Devonian black shale gas field; the Cottageville (Mount Alto) Gas Field, Jackson County, West Virginia: a case study of Devonian shale gas production; the Eastern Kentucky Gas Field: a geological study of the relationshipsmore » of Ohio Shale gas occurrences to structure, stratigraphy, lithology, and inorganic geochemical parameters; and a statistical analysis of geochemical data for the Eastern Kentucky Gas Field.« less

Late-Quaternary recharge determined from chloride in shallow groundwater in the central Great Plains

USGS Publications Warehouse

Macfarlane, P.A.; Clark, J.F.; Davisson, M.L.; Hudson, G.B.; Whittemore, Donald O.

2000-01-01

An extensive suite of isotopic and geochemical tracers in groundwater has been used to provide hydrologic assessments of the hierarchy of flow systems in aquifers underlying the central Great Plains (southeastern Colorado and western Kansas) of the United States and to determine the late Pleistocene and Holocene paleotemperature and paleorecharge record. Hydrogeologic and geochemical tracer data permit classification of the samples into late Holocene, late Pleistocene-early Holocene, and much older Pleistocene groups. Paleorecharge rates calculated from the Cl concentration in the samples show that recharge rates were at least twice the late Holocene rate during late Pleistocene-early Holocene time, which is consistent with their relative depletion in 16O and D. Noble gas (Ne, Ar, Kr, Xe) temperature calculations confirm that these older samples represent a recharge environment approximately 5??C cooler than late Holocene values. These results are consistent with the global climate models that show a trend toward a warmer, more arid climate during the Holocene. (C) 2000 University of Washington.

Mantle convection and plate tectonics: toward an integrated physical and chemical theory

PubMed

Tackley

2000-06-16

Plate tectonics and convection of the solid, rocky mantle are responsible for transporting heat out of Earth. However, the physics of plate tectonics is poorly understood; other planets do not exhibit it. Recent seismic evidence for convection and mixing throughout the mantle seems at odds with the chemical composition of erupted magmas requiring the presence of several chemically distinct reservoirs within the mantle. There has been rapid progress on these two problems, with the emergence of the first self-consistent models of plate tectonics and mantle convection, along with new geochemical models that may be consistent with seismic and dynamical constraints on mantle structure.

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

USGS Publications Warehouse

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

2000-01-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 O2 and SO4/2- in the aerobic zone, and to the reduction of SO4/2- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55 x 10-4 to 48.6 x 10-4 mmol 1-1 yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO4/2- (1.31 x 10-4 to 15 x 10-4 mmol 1-1 yr-1). The overall increase in SO4/2- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO4/2- loss through microbial reduction is exceeded by SO4/2- gain through diffusion from sediments and through the oxidation of FeS2. Geochemichal and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO4/2- -rich zones have localized SO4/2- -reducing zones in which the formation of iron disulfides been depleted by microbial reduction and resulted in decreases dissolved iron concentrations. These localized zones of SO4/2- 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. (C) 2000 Elsevier Science B.V.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 O2 and SO42- in the aerobic zone, and to the reduction of SO42- in the anaerobic zone; estimated rates of CO2 production through reduction of Fe(III) were relatively minor by comparison. The rates of CO2 production calculated from dissolved inorganic carbon mass transfer (2.55??10-4 to 48.6??10-4mmol l-1yr-1) generally were comparable to the calculated rates of CO2 production by the combined reduction of O2, Fe(III) and SO42- (1.31??10-4 to 15??10-4mmol l-1yr-1). The overall increase in SO42- concentrations along the flow path, together with the results of mass-balance calculations, and variations in ??34S values along the flow path indicate that SO42- loss through microbial reduction is exceeded by SO42- gain through diffusion from sediments and through the oxidation of FeS2. Geochemical and microbial data on cores indicate that Fe(III) oxyhydroxide coatings on sediment grains in local, organic carbon- and SO42--rich zones have been depleted by microbial reduction and resulted in localized SO42--reducing zones in which the formation of iron disulfides decreases dissolved iron concentrations. These localized zones of SO42- 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.

A model for the oceanic mass balance of rhenium and implications for the extent of Proterozoic ocean anoxia

NASA Astrophysics Data System (ADS)

Sheen, Alex I.; Kendall, Brian; Reinhard, Christopher T.; Creaser, Robert A.; Lyons, Timothy W.; Bekker, Andrey; Poulton, Simon W.; Anbar, Ariel D.

2018-04-01

Emerging geochemical evidence suggests that the atmosphere-ocean system underwent a significant decrease in O2 content following the Great Oxidation Event (GOE), leading to a mid-Proterozoic ocean (ca. 2.0-0.8 Ga) with oxygenated surface waters and predominantly anoxic deep waters. The extent of mid-Proterozoic seafloor anoxia has been recently estimated using mass-balance models based on molybdenum (Mo), uranium (U), and chromium (Cr) enrichments in organic-rich mudrocks (ORM). Here, we use a temporal compilation of concentrations for the redox-sensitive trace metal rhenium (Re) in ORM to provide an independent constraint on the global extent of mid-Proterozoic ocean anoxia and as a tool for more generally exploring how the marine geochemical cycle of Re has changed through time. The compilation reveals that mid-Proterozoic ORM are dominated by low Re concentrations that overall are only mildly higher than those of Archean ORM and significantly lower than many ORM deposited during the ca. 2.22-2.06 Ga Lomagundi Event and during the Phanerozoic Eon. These temporal trends are consistent with a decrease in the oceanic Re inventory in response to an expansion of anoxia after an interval of increased oxygenation during the Lomagundi Event. Mass-balance modeling of the marine Re geochemical cycle indicates that the mid-Proterozoic ORM with low Re enrichments are consistent with extensive seafloor anoxia. Beyond this agreement, these new data bring added value because Re, like the other metals, responds generally to low-oxygen conditions but has its own distinct sensitivity to the varying environmental controls. Thus, we can broaden our capacity to infer nuanced spatiotemporal patterns in ancient redox landscapes. For example, despite the still small number of data, some mid-Proterozoic ORM units have higher Re enrichments that may reflect a larger oceanic Re inventory during transient episodes of ocean oxygenation. An improved understanding of the modern oceanic Re cycle and a higher temporal resolution for the Re compilation will enable further tests of these hypotheses regarding changes in the surficial Re geochemical cycle in response to variations in atmosphere-ocean oxygenation. Nevertheless, the existing Re compilation and model results are in agreement with previous Cr, Mo, and U evidence for pervasively anoxic and ferruginous conditions in mid-Proterozoic oceans.

Biokinetics of yttrium and comparison with its geochemical twin holmium

DOE PAGES

Leggett, Rich

2017-06-01

The transition metal yttrium (Y, atomic number 39) is chemically similar to elements in the lanthanide family (atomic numbers 57-71, lanthanum through lutetium) and is always present with the lanthanides in rare earth ores. Yttrium and the lanthanide holmium are particularly close chemical and physical analogues and are referred to as geochemical twins because they typically show little fractionation in geological material. Extensive measurements on rocks, soils, and meteorites indicate that the Y/Ho mass concentration ratio rarely falls far from the “chondritic” or “solar system” ratio of ~26. Our paper presents a new biokinetic model for yttrium in adult humansmore » and examines whether yttrium and holmium may be biological as well as geochemical twins. Collected data on yttrium and holmium in plants and human tissues do not allow precise derivations of Y/Ho concentration ratios but with occasional exceptions yield ratios that are reasonably consistent with chondritic values. Predictions of the time-dependent behavior of yttrium in adult humans based on the yttrium model presented here closely approximate predictions of the behavior of holmium based on a previously developed model for holmium. We know that yttrium and holmium are close biological analogues, but the available comparative data are too limited and imprecise to reveal whether there are any significant differences in their biological behavior.« less

Geologic, geophysical, and geochemical aspects of site-specific studies of the geopressured-geothermal energy resource of southern Louisiana. Final report

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

Pilger, R.H. Jr.

1985-01-01

The report consists of four sections dealing with progress in evaluating geologic, geochemical, and geophysical aspects of geopressured-geothermal energy resources in Louisiana. Separate abstracts have been prepared for the individual sections. (ACR)

Self-gravity, self-consistency, and self-organization in geodynamics and geochemistry

NASA Astrophysics Data System (ADS)

Anderson, Don L.

The results of seismology and geochemistry for mantle structure are widely believed to be discordant, the former favoring whole-mantle convection and the latter favoring layered convection with a boundary near 650 km. However, a different view arises from recognizing effects usually ignored in the construction of these models, including physical plausibility and dimensionality. Self-compression and expansion affect material properties that are important in all aspects of mantle geochemistry and dynamics, including the interpretation of tomographic images. Pressure compresses a solid and changes physical properties that depend on volume and does so in a highly nonlinear way. Intrinsic, anelastic, compositional, and crystal structure effects control seismic velocities; temperature is not the only parameter, even though tomographic images are often treated as temperature maps. Shear velocity is not a good proxy for density, temperature, and composition or for other elastic constants. Scaling concepts are important in mantle dynamics, equations of state, and wherever it is necessary to extend laboratory experiments to the parameter range of the Earth's mantle. Simple volume-scaling relations that permit extrapolation of laboratory experiments, in a thermodynamically self-consistent way, to deep mantle conditions include the quasiharmonic approximation but not the Boussinesq formalisms. Whereas slabs, plates, and the upper thermal boundary layer of the mantle have characteristic thicknesses of hundreds of kilometers and lifetimes on the order of 100 million years, volume-scaling predicts values an order of magnitude higher for deep-mantle thermal boundary layers. This implies that deep-mantle features are sluggish and ancient. Irreversible chemical stratification is consistent with these results; plausible temperature variations in the deep mantle cause density variations that are smaller than the probable density contrasts across chemical interfaces created by accretional differentiation and magmatic processes. Deep-mantle features may be convectively isolated from upper-mantle processes. Plate tectonics and surface geochemical cycles appear to be entirely restricted to the upper ˜1,000 km. The 650-km discontinuity is mainly an isochemical phase change but major-element chemical boundaries may occur at other depths. Recycling laminates the upper mantle and also makes it statistically heterogeneous, in agreement with high-frequency scattering studies. In contrast to standard geochemical models and recent modifications, the deeper layers need not be accessible to surface volcanoes. There is no conflict between geophysical and geochemical data, but a physical basis for standard geochemical and geodynamic mantle models, including the two-layer and whole-mantle versions, and qualitative tomographic interpretations has been lacking.

Petrographic and geochemical characteristic of volcanic rocks from Tasik Kenyir and Kampung Awah, East Malaya block, Peninsular Malaysia

NASA Astrophysics Data System (ADS)

Roselee, Muhammad Hatta; Umor, Mohd Rozi; Ghani, Azman Abdul; Badruldin, Muhamad Hafifi; Quek, Long Xiang

2018-04-01

Kampung Awah and Tasik Kenyir are geologically located in East Malaya Blocks. These block is also known as western margin of Indochina terrane. Apart from sedimentary formations, East Malaya Blocks is also dominated by plutonic and volcanic rocks of mafic to rhyolitic compositions. Petrography and geochemical data suggest that Kampung Awah and Tasik Kenyir are one of locations which consists of volcanic rocks of generally basaltic to basaltic andesite compositions. Volcanic rocks from both area consists of plagioclcase, clinopyroxene, orthpyroxene as main mineral constituents with minor occurrences of hornblende. Geochemical data also indicate that volcanic rocks from both area were formed during subduction of the Paleo-tethys oceanic underneath the East Malaya Block or Indochina terrane. Most of the samples are metaluminous which indicate the volcanics are derived from igneous origin. This paper will contribute new geochemical data of mafic volcanics from Kampung Awah and Tasik Kenyir with the support of petrographic and field evidence to deduce the magma evolution and the tectonic setting.

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.

Hydrochemical characteristics and groundwater evolution modeling in sedimentary rocks of the Tono mine, Japan

NASA Astrophysics Data System (ADS)

Sasamoto, Hiroshi; Yui, Mikazu; Arthur, Randolph C.

Based on geochemical data collected by Japan Nuclear Cycle Development Institute (JNC) in the Tono uranium mine, a conceptual groundwater evolution model developed by JNC is tested to evaluate whether equilibrium-based concepts of water-rock interaction are consistent with observed variations in the mineralogy and hydrochemistry of the Tono mine area. The chemical evolution of the groundwaters is modeled assuming local equilibrium for selected mineral-fluid reactions, taking into account the rainwater origin of these solutions. Results suggest that it is possible to interpret approximately the actual groundwater chemistry (i.e., pH, Eh, total dissolved concentrations of Si, Na, Ca, K, Al, carbonate and sulfate) if the following assumptions are adopted (a) CO 2 concentration in the gas phase contacting pore solutions in the overlying soil zone=10 -1 atm, and (b) minerals in the rock zone that control the solubility of respective elements in the groundwater include: chalcedony (Si), albite (Na), kaolinite (Al), calcite (Ca and carbonate), muscovite (K) and pyrite (Eh and sulfate). This result helps to build confidence in the use of simplified geochemical modeling techniques to develop an understanding of dominant geochemical reactions controlling groundwater chemistry in rocks similar to those that could be used for the geological disposal of radioactive wastes. It is noted, however, that the available field data may not be sufficient to adequately constrain parameters in the groundwater evolution model. In particular, more detailed information characterizing certain site properties are needed to improve the model. For this reason, a model that accounts for ion-exchange reactions among clay minerals, and which is based on the results of laboratory experiments, has also been evaluated in the present study. Further improvement of model considering ion-exchange reactions are needed in future, however.

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.

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.

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

USGS Publications Warehouse

Siegel, M.D.; Anderholm, S.

1994-01-01

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

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.

Geochemical landscapes of the conterminous United States; new map presentations for 22 elements

USGS Publications Warehouse

Gustavsson, N.; Bolviken, B.; Smith, D.B.; Severson, R.C.

2001-01-01

Geochemical maps of the conterminous United States have been prepared for seven major elements (Al, Ca, Fe, K, Mg, Na, and Ti) and 15 trace elements (As, Ba, Cr, Cu, Hg, Li, Mn, Ni, Pb, Se, Sr, V, Y, Zn, and Zr). The maps are based on an ultra low-density geochemical survey consisting of 1,323 samples of soils and other surficial materials collected from approximately 1960-1975. The data were published by Boerngen and Shacklette (1981) and black-and-white point-symbol geochemical maps were published by Shacklette and Boerngen (1984). The data have been reprocessed using weighted-median and Bootstrap procedures for interpolation and smoothing.

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.

Multivariate analysis of the heterogeneous geochemical processes controlling arsenic enrichment in a shallow groundwater system.

PubMed

Huang, Shuangbing; Liu, Changrong; Wang, Yanxin; Zhan, Hongbin

2014-01-01

The effects of various geochemical processes on arsenic enrichment in a high-arsenic aquifer at Jianghan Plain in Central China were investigated using multivariate models developed from combined adaptive neuro-fuzzy inference system (ANFIS) and multiple linear regression (MLR). The results indicated that the optimum variable group for the AFNIS model consisted of bicarbonate, ammonium, phosphorus, iron, manganese, fluorescence index, pH, and siderite saturation. These data suggest that reductive dissolution of iron/manganese oxides, phosphate-competitive adsorption, pH-dependent desorption, and siderite precipitation could integrally affect arsenic concentration. Analysis of the MLR models indicated that reductive dissolution of iron(III) was primarily responsible for arsenic mobilization in groundwaters with low arsenic concentration. By contrast, for groundwaters with high arsenic concentration (i.e., > 170 μg/L), reductive dissolution of iron oxides approached a dynamic equilibrium. The desorption effects from phosphate-competitive adsorption and the increase in pH exhibited arsenic enrichment superior to that caused by iron(III) reductive dissolution as the groundwater chemistry evolved. The inhibition effect of siderite precipitation on arsenic mobilization was expected to exist in groundwater that was highly saturated with siderite. The results suggest an evolutionary dominance of specific geochemical process over other factors controlling arsenic concentration, which presented a heterogeneous distribution in aquifers. Supplemental materials are available for this article. Go to the publisher's online edition of the Journal of Environmental Science and Health, Part A, to view the supplemental file.

Generation of mantle heterogeneity by oceanic crust recycling: how well can we match geochemical and geophysical observations? (Invited)

NASA Astrophysics Data System (ADS)

van Keken, P. E.; Brandenburg, J. P.; Hauri, E. H.; Ballentine, C. J.

2009-12-01

The heterogeneity of the Earth's mantle is expressed in complementary geochemical and geophysical signatures, where the geochemistry provides a time-integrated signal and the geophysics tends to see a recent snapshot of the Earth's interior. While the geophysical evidence tends to support a form of whole mantle convection that is moderated by rheological and phase changes below the transition zone, the geochemical observations have been generally used to support the presence of long-lived and isolated reservoirs. Recent dynamical modeling (Brandenburg et al., EPSL, 2008) employed high resolution finite modeling of mantle convection using an energetically consistent simulation of tectonic plates. A suite of models was developed with a dynamic vigor similar to that of the present day earth. The recycling of oceanic crust combined with a two-stage formation of the continental crust leads to a satisfactory match to the observed spread between HIMU-DMM-EM1 in multiple isotope systems without invoking recycling of continental crust. Due to the rheological contrast between upper and lower mantle there is a natural occurrence of a well-mixed upper mantle overlaying a chemically more heterogeneous lower mantle. The pooling of dense oceanic crust provides the formation of dense piles at the base of the mantle. Together with the occurrence of slabs that thicken and/or stagnate at the 670 discontinuity we find reasonable correspondance with the present day tomographic signatures. At present the models fail to explain noble gas systematics, even when taking the suggested high compatibility of helium into account.

Application of artificial neural networks to chemostratigraphy

NASA Astrophysics Data System (ADS)

Malmgren, BjöRn A.; Nordlund, Ulf

1996-08-01

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

Mineral and Geochemical Classification From Spectroscopy/Diffraction Through Neural Networks

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

Synthesis of soil geochemical characteristics and organic carbon degradation in Arctic polygon tundra, Barrow, Alaska

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

Zheng, Jianqiu; RoyChowdhury, Taniya; Herndon, Elizabeth

This is a synthesis data product that reports (1) the results of soil geochemical characterization and (2) organic carbon degradation in low temperature soil incubations on cores collected on the NGEE Arctic Study Area near Utqiaġvik (Barrow), Alaska. The study area consists of thaw lakes, drained thaw lake basins and interstitial tundra with a polygonal landscape of microtopographic features created by ice wedges. Integrated geochemical and organic carbon degradation data from 9 individual soil cores are included in the synthesis product.

Under an Orange Sky: The Many Implications of Organic Haze for Earthlike Planets

NASA Astrophysics Data System (ADS)

Arney, Giada; Domagal-Goldman, Shawn D.; Meadows, Victoria S.; Wolf, Eric; Schwieterman, Edward W.; Charnay, Benjamin; Claire, Mark; Hebrard, Eric

2015-11-01

Geochemical evidence suggests Archean Earth was intermittently enshrouded in an organic haze resulting from methane photolysis. Hazy exoplanets may be common, and hazes can significantly impact the environment of habitable planets. Earth is frequently studied as an analog for habitable exoplanets, and Archean Earth is the most alien planet we have geochemical data for. We have used 1D photochemical-climate and radiative transfer simulations to examine the climate, surface radiation environment, and spectra of Archean Earth with fractal hydrocarbon haze. We find that haze would have strongly impacted Earth’s climate, lowering the planetary surface temperature by 20-30 K. However, this cooling can be countered by concentrations of greenhouses gases consistent with geochemical constraints. For example, an atmosphere with 2% CO2, 0.37% CH4 and a self-consistent hydrocarbon haze has a globally averaged surface temperature of 274 K, which GCM models have shown is consistent with a large open ocean fraction (Charnay et al 2013). The cooling from haze means that there exists a “hazy habitable zone” closer to the star than the traditional habitable zone boundaries. Our results suggest that the hazy habitable zone can extend to the distance of Venus. An organic haze produces strong, remotely detectable spectral features, especially at wavelengths < 0.5 μm, reddening the planet’s color. The strong absorption of UV radiation by this haze means it could have provided a UV shield for the Archean Earth prior to the rise of oxygen when there was no ozone layer: we show that an organic haze can block 97% of the surface-incident UVC (λ < 0.28 μm) radiation compared to a haze-free planet. UVC radiation directly dissociates DNA, and it is blocked by ozone in the modern atmosphere. Organic hazes may therefore benefit surface biospheres on Earth and similar exoplanets. Finally, assuming geochemical constraints on the Archean atmospheric composition, we show that abiotic levels of methane flux to the atmosphere are insufficient to form an organic haze. For Earthlike exoplanets, organic haze may therefore be a novel type of spectral biosignature.

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.

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.