Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

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

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.

Boyce, Jeremy W.↗

Application of a Geochemically Informed Leak Detection (GILD) Model to CO 2 Injection Sites on the United States Gulf Coast

The Gulf Coast region possesses great potential for CO 2 enhanced oil recovery (EOR) and CO 2 storage. A geochemically informed leak detection (GILD) model has been applied to CO 2 injection sites on the Gulf Coast with considerations of measurement variability. The Jasper aquifer in Montgomery County, Texas, was chosen to demonstrate the method. Based on background data from wells in the area, combinations of mineral and fluid compositions were used to create 23 scenarios for the geochemical model. The output from the geochemical model was used to identify sensitive monitoring species, and response functions were generated for these as a function of the CO 2 leakage concentration. The sources of measurement variability for background conditions were characterized from the Jasper aquifer background data, and then normalized using the coefficient of variation of each species across the monitoring wells. Bayesian belief network (BBN) models were constructed, and measurement variability of different levels were added to compare leak detection probabilities. Increasing measurement variability decreased the power to detect a leak of a given size. For a moderately high CO 2 concentration of 0.2 mol/kg, the probability of detecting this leakage effect using pH as the monitoring variable in an aquifer with calcite decreases from 98% (no measurement variability) to 61% (medium variability) to 33% (high variability). The loss in power of the sampling protocol with increasing measurement variability is similar in magnitude when Ca 2+ or HCO 3 - is used as the monitoring parameter, but only for aquifers with calcite.

58 GEOSCIENCES↗

A hybrid model of the CO2 geochemical cycle and its application to large impact events

The effects of a large asteriod or comet impact on modern and ancient marine biospheres are analyzed. A hybrid model of the carbonate-silicate geochemical cycle, which is capable of calculating the concentrations of carbon dioxide in the atmosphere, ocean, and sedimentary rocks, is described. The differences between the Keir and Berger (1983) model and the hybrid model are discussed. Equilibrium solutions are derived for the preindustrial atmosphere/ocean system and for a system similar to that of the late Cretaceous Period. The model data reveal that globl darkening caused by a stratospheric dust veil could destroy the existing phytoplankton within a period of several weeks or months, nd the dissolution of atmospheric NO(x) compounds would lower the pH of ocean surface waters and release CO2 into the atmosphere. It is noted that the surface temperatures could be increased by several degrees and surface oceans would be uninhabitable for calcaerous organisms for approximately 20 years.

Kasting, J. F.↗

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

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.

Spudis, P. D.↗

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

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.

Plumlee, Geoffrey S.↗

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗

Petrogenetic modeling of Hawaiian tholeiitic basalts - A geochemical approach

The abundances of 29 elements in 33 samples of tholeiitic basalts from five volcanoes in Hawaii are determined by neutron-activation analysis; the results are presented in tables and graphs with the data of Murali et al. (1977) for a sixth volcano and characterized in detail; and geochemical partial-melting models are constructed to explain the origin of these basalts. The models proposed are based on olivine crystallization from three distinct source compositions (determined from the REE and Sc partitioning) and correspond to the volcano groups Mauna Kea, Kohala, and Kilauea; Mauna Loa and Lanai; and Koolau.

Budahn, J. R.↗

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗

Impact of oxidation and carbonation on the release rates of iodine, selenium, technetium, and nitrogen from a cementitious waste form

Evaluation of the long-term retention mechanisms and potential release rates for the primary constituents of potential concern (COPCs) (i.e., Tc, I, Se, and nitrate) is necessary to determine if Cast Stone, a radioactive waste form, can meet performance objectives under near-surface disposal scenarios. In this work, a mineral and parameter set accounting for the solubility of I and Se in Cast Stone was developed based on pH-dependent and monolithic diffusion leaching test results, to extend a geochemical speciation model previously developed. The impact of oxidation and carbonation as environmental aging processes on the retention properties of Cast Stone for primary COPCs was systematically estimated. Physically, the effective diffusion coefficients of 4 COPCs in Cast Stone were increased after carbonation and/or oxidation, reflecting an increase in permeability to diffusion. Chemically, i) pH & pe conditions in the original Cast Stone were favorable for the stabilization of Tc, but not for I, Se, and N; ii) oxidation (with/without carbonation) of Cast Stone changed the pe & pH conditions to be detrimental for Tc stabilization; and iii) carbonation (with/without oxidation) of Cast Stone modified the pH & pe conditions to be beneficial for the stabilization of I (in system with Ag added) and Se.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Prediction of Long-Term Geochemical Change in Bentonite Based on the Interpretative THMC Model of the FEBEX In Situ Test

Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.

Zheng, Liange↗

Impact of carbonation on leaching of constituents from a cementitious waste form for treatment of low activity waste at the DOE Hanford site

Carbonation can be a major aging process during disposal of alkaline cementitious waste forms and can impact constituent leaching by changes in material alkalinity, pore structure, and controlling mineral phases. The effect of carbonation on the leaching of major and trace constituents from Cast Stone, a cementitious waste form developed to treat high salt content low activity waste, was studied through a combination of leaching experiments and reactive transport simulations. Diffusive transport of constituents in the waste form was evaluated using reactive transport modeling of diffusion-controlled leaching test results and a geochemical speciation model derived from pH-dependent leaching. Comparisons between Cast Stone materials aged under nitrogen, air, and 2% carbon dioxide in nitrogen showed that carbonation impacts solubility, physical retention and observed diffusivity of major and trace constituents. Carbonation under 2% CO 2 decreased the diffusion-controlled leaching of chromium by two orders of magnitude. Finally, modeling results suggest that carbonation may also decrease solubility of technetium while changes to microstructure by carbonation increases effective diffusivity of constituents in Cast Stone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Little Ice Age and Human-Environmental Interactions in the Central Balkans: Insights from a New Serbian Paleorecord

This paper presents a 600-year well-dated, high-resolution Central Balkan paleo-record including the Little Ice Age (LIA; 1450-850 CE). Utilizing pollen-based REVEALS modelling estimates, geochemical indicators, rarefaction analyses and the AMS 14C-based Bacon age model, this first-hand record from the Sava Basin reveals the transformation of the Central Balkan landscape involving linkages between changing climatic and socio-political regimes. The pre-LIA interval (1370-1418 CE) in the Sava Region reveals a wooded steppe and increased cultivation under warmer/stable climatic and socio-political conditions. In contrast, the LIA interval in the region is expressed through continuous transitions between forest and grassland, extensive land erosion and stressed agriculture, potentially as a collective artifact of the climatic variability and human impact associated with socio- political stressors of the time. The post-LIA/Industrial Era interval (1850-2012 CE) in the Sava Region shows an overall increase in woodland as well as agriculture following the exit of the Ottomans. However, increased population pressures and the subsequent onset of the Industrial Revolution as well as increased trade led to intense deforestation throughout the 20th Century, which continued during the Socialist period.

Societal impacts of climate change↗

Comments on the BLAG model - The carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years

A self-consistent method of determining initial conditions for the model presented by Berner, Lasaga, and Garrels (1983) (henceforth, the BLAG model) is derived, based on the assumption that the CO2 geochemical cycle was in steady state at t = -100 m.y. (million years). This initialization procedure leads to a dissolved magnesium concentration higher than that calculated by Berner, Lasaga, and Garrels and to a low ratio of dissolved calcium to bicarbonate prior to 60 m.y. ago. The latter prediction conflicts with the geologic record of evaporite deposits, which requires that this ratio remain greater than 0.5. The contradiction is probably caused by oversimplifications in the BLAG model, such as the neglect of the cycles of organic carbon and sulfur.

Kasting, J. F.↗

Comments on the BLAG model: the carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years

A self-consistent method of determining initial conditions for the model presented by Berner, Lasaga, and Garrels (1983) (henceforth, the BLAG model) is derived, based on the assumption that the CO2 geochemical cycle was in steady state at t = -100 my (million years). This initialization procedure leads to a dissolved magnesium concentration higher than that calculated by Berner, Lasaga, and Garrels and to a low ratio of dissolved calcium to bicarbonate prior to 60 my ago. The latter prediction conflicts with the geologic record of evaporite deposits, which requires that this ratio remain greater than 0.5. The contradiction is probably caused by oversimplifications in the BLAG model, such as the neglect of the cycles of organic carbon and sulfur.

Evolution, Planetary↗

The Martian ocean: First acid, then alkaline

In Mars' distant past, carbon dioxide and water may have been plentiful. Values of total outgassed CO2 from several to about 10 bar are consistent with present knowledge, and this amount of CO2 implies an amount of water outgassed equal to an equivalent depth of 500 to 1000 m. It is quite reasonable, therefore, to envision an early Mars in which there was a body or bodies of liquid water, perhaps in the northern plains, and a dense carbon dioxide atmosphere. Under such conditions, the pH of the water will be low, due to the dissolution of carbon dioxide in the water to form carbonic acid. This acidic water is capable of weathering the available rock quite intensely, particularly because this rock is likely to be heavily fractured (from meteorite bombardment) or even consist of fine particles (such as pyroclastic deposits). As time goes on, however, the carbon dioxide atmosphere will rapidly pass through the ocean to form carbonate deposits. As the density of the atmosphere decreases, so will the flux of carbonic acid into the ocean. Without this input of carbonic acid, the effect of the dissolved weathering products will be to increase the pH of the water. The ocean will then become alkaline. To study this process, I have developed a geochemical cycle model for the atmosphere-hydrosphere-regolith system of Mars. The treatment of geochemical cycles as complex kinetic chemical reactions has been undertaken for terrestrial systems in recent years with much success. This method is capable of elegantly handling the interactions between the simultaneous chemical reactions needed to understand such a system.

Schaefer, M. W.↗

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗

The relation between the age of the subconducting slab and the recycling of sediments into the mantle

The recycling of sediments into the mantle has become an important issue because recent papers have suggested that the geochemical inverse models of the evolution of radiogenic isotope abundances over the history of the Earth have nonunique solutions. Both the recycling of continent-derived sediments into the mantle and mixing in the mantle could produce similar geochemical effects in the mean isotopic ratios of new igneous material emplaced in continents. Recent models of Archaean heat flow and of plate tectonics during early Earth history have demonstrated that higher internal heat production of the early Earth was mainly dissipated through a higher creation rate of oceanic lithosphere. If the seafloor creation rate was higher on the early Earth, then the residence time of any one piece of oceanic lithosphere on the surface would have been shorter. It is possible that a higher rate of recycling of oceanic lithosphere into the mantle could have resulted in some transport of sediment into the mantle.

Abbott, D.↗