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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Oxygen Depletion on the Surface of Mercury: Evidence of Silicon Smelting?

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft collected data that provided important insights into the structure, chemical makeup, and compositional diversity of Mercury. The X-Ray Spectrometer (XRS) and Gamma-Ray Spectrometer (GRS) onboard MESSENGER provided the first detailed chemical analyses of Mercury's surface. Among the many discoveries included several surprising characteristics about the surface of Mercury, including elevated S abundances (up to 4 percent by weight), low Fe abundances (less than 4 percent by weight), and relatively low O abundances (O/Si ratio of 1.40 plus or minus 0.03). The surface chemistry as determined by MESSENGER has been used to identify up to nine distinct geochemical terranes on Mercury. Numerous modeling and experimental efforts have been undertaken to infer the mineralogy and petrology of mercurian lavas and surface materials. However, all of these efforts have presumed valence states for each of the elements according to the following: Si4+, Ti4+, Al3+, Cr2+, Fe2+, Mn2+, Mg2+, Ca2+, Na+, K+, S2-, Cl-. Based on these valence assignments, cations are charged balanced with the anions O2-, S2-, and Cl- and the compositions are recast in terms of oxides, sulfides, and chlorides. Based on these assumptions, the geochemical terranes that have been identified on Mercury yield O/Si wt. ratios ranging from 1.61 to 1.84, which is substantially higher than the preliminary O/Si ratio of 1.40 plus or minus 0.03 determined by the MESSENGER GRS]. We have re-evaluated the O/Si ratio using the entire MESSENGER dataset to reassess its implications for the geochemistry of Mercury.

McCubbin, F. M.↗
DOE OSTI2024

Hydrologic windows into the crystalline basement and their controls on groundwater flow patterns across the Paradox Basin, western USA

Conceptual models of sedimentary basin groundwater flow systems typically assume that the crystalline basement acts as an impermeable boundary and can be neglected. In this study, we use hydrologic models constrained by isotopic and geochemical datasets to argue that the La Sal Mountains, Utah, USA, act as a hydrologic window into the Paradox Basin’s lower aquifer system and underlying crystalline basement. We conducted a sensitivity study in which we varied crystalline basement/laccolith permeability as well as fault zone connectivity along a cross-sectional transect from the La Sal Mountains to Lisbon Valley. When the crystalline basement/laccolith units are set at relatively permeable levels (10–14 m 2 ), simulated tracers that include total dissolved solids, oxygen isotopic composition of pore fluids (δ 18 O), and groundwater residence times are in closest agreement with field measurements. Model results indicate that pore fluids in the basal aquifer system underlying the Paradox Formation confining unit are a mixture of relatively young meteoric fluids and older Paradox Formation brines. The presence of faults did not significantly modify fluid exchange between the upper and lower aquifer systems. This was due, in part, to underpressuring within the Paradox Formation. Our study concludes that the Paradox Basin represents a regional recharge area for the Colorado Plateau, with groundwater discharge occurring along the Colorado River within the Grand Canyon some 375 km away to the southwest. This is only possible with a permeable crystalline basement. Our findings help explain the genesis of Mississippi Valley-type ore deposits of the US Midcontinent, where the presence of a permeable basement may be useful in addressing issues related to solute mass and energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗
DOE OSTI2023

Tonian Carbonates Record Phosphate-Rich Shallow Seas

The early-middle Neoproterozoic is thought to have witnessed significant perturbations to marine P cycling, in turn facilitating the rise of eukaryote-dominated primary production. However, with few robust constraints on aqueous P concentrations, current understanding of Neoproterozoic P cycling is generally model-dependent. To provide new geochemical constraints, we combined microanalytical data sets with solid-state Nuclear Magnetic Resonance, synchrotron-based X-ray Absorption Near Edge Structure spectroscopy, and micro-X-ray Fluorescence imaging to characterize the speciation and distribution of P in Tonian shallow-water carbonate rocks. These data reflect shallow water phosphate concentrations 10–100× higher than modern systems, supporting the hypothesis that tectonically-driven influxes in P periodically initiated kinetically-controlled CaCO 3 deposition, in turn destabilizing marine carbonate chemistry, climate, and nutrient inventories. Alongside these observations, a new compilation and statistical analysis of mudstone geochemistry data indicates that, in parallel, C org and P burial increased across later Tonian continental margins until becoming decoupled at the close of the Tonian, implicating widespread N-limitation triggered by increasing atmospheric O 2 .

58 GEOSCIENCES↗
NASA NTRS1974

The evolution of the moon and the terrestrial planets

The thermal evolutions of the Moon, Mars, Venus and Mercury are calculated theoretically starting from cosmochemical condensation models. An assortment of geological, geochemical and geophysical data are used to constrain both the present day temperatures and the thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history. The moon, smallest in size, is characterized as a differentiated body with a crust, a thick solid mantle and an interior region which may be partially molten. Mars, intermediate in size, is assumed to have differentiated an Fe-FeS core. Venus is characterized as a planet not unlike the earth in many respects. Core formation has occurred probably during the first billion years after the formation. Mercury, which probably has a large core, may have a 500 km thick solid lithosphere and a partially molten core if it is assumed that some heat sources exist in the core.

Toksoez, M. N.↗
NASA NTRS1977

The evolution of the moon and the terrestrial planets

The thermal evolutions of the Moon, Mars, Venus, and Mercury were calculated theoretically starting from cosmochemical condensation models. An assortment of geological, geochemical, and geophysical data were used to constrain both the present day temperature and the thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history.

Toksoez, M. N.↗
NASA NTRS1987

Digital cartography of the planets - New methods, its status, and its future

A system has been developed that establishes a standardized cartographic database for each of the 19 planets and major satellites that have been explored to date. Compilation of the databases involves both traditional and newly developed digital image processing and mosaicking techniques, including radiometric and geometric corrections of the images. Each database, or digital image model (DIM), is a digital mosaic of spacecraft images that have been radiometrically and geometrically corrected and photometrically modeled. During compilation, ancillary data files such as radiometric calibrations and refined photometric values for all camera lens and filter combinations and refined camera-orientation matrices for all images used in the mapping are produced. The system is designed to aid in the following: access to and utilization of spatial data returned by planetary spacecraft; coregistration of cartographic image and topographic elevation models with geological, geophysical, and geochemical databases; and the design and preparation of both standard and custom map products. A major goal of the project is to publish the DIMs on compact computer-compatible media such as digital-optical (CD-ROM) disks.

Batson, R. M.↗
NASA NTRS1992

Data management for JGOFS: Theory and design

The Joint Global Ocean Flux Study (JGOFS), currently being organized under the auspices of the Scientific Committee for Ocean Research (SCOR), is intended to be a decade long internationally coordinated program. The main goal of JGOFS is to determine and understand on a global scale the processes controlling the time-varying fluxes of carbon and associated biogenic elements in the ocean and to evaluate the related exchanges with the atmosphere, sea floor and continental boundaries. 'A long-term goal of JGOFS will be to establish strategies for observing, on long time scales, changes in ocean biogeochemical cycles in relation to climate change'. Participation from a large number of U.S. and foreign institutions is expected. JGOFS investigators have begun a set of time-series measurements and global surveys of a wide variety of biological, chemical and physical quantities, detailed process-oriented studies, satellite observations of ocean color and wind stress and modeling of the bio-geochemical processes. These experiments will generate data in amounts unprecedented in the biological and chemical communities; rapid and effortless exchange of these data will be important to the success of JGOFS.

Flierl, Glenn R.↗
NASA NTRS1991

Simulations of the HDO and H2O-18 atmospheric cycles using the NASA GISS general circulation model - Sensitivity experiments for present-day conditions

Incorporating the full geochemical cycles of stable water isotopes (HDO and H2O-18) into an atmospheric general circulation model (GCM) allows an improved understanding of global delta-D and delta-O-18 distributions and might even allow an analysis of the GCM's hydrological cycle. A detailed sensitivity analysis using the NASA/Goddard Institute for Space Studies (GISS) model II GCM is presented that examines the nature of isotope modeling. The tests indicate that delta-D and delta-O-18 values in nonpolar regions are not strongly sensitive to details in the model precipitation parameterizations. This result, while implying that isotope modeling has limited potential use in the calibration of GCM convection schemes, also suggests that certain necessarily arbitrary aspects of these schemes are adequate for many isotope studies. Deuterium excess, a second-order variable, does show some sensitivity to precipitation parameterization and thus may be more useful for GCM calibration.

Jouzel, Jean↗
DOE OSTI2022

QA/QC of the East River, Colorado, discharge and geochemical time series datasets (Almont, BCC, and Pump House) to be used for modeling of hydrogeochemical balance

The following datasets were QA/QC-ed (Quality Assurance/Quality Control): 1. Brush Creek Confluence (BCC) discharge data (from Helen Malenda, USGS, Colorado School of Mines), which were calculated using the pressure transducer data and rating curves. The original 15 min time series data were presented as mean daily discharge. 2. Almont discharge data from United States Geological Survey (USGS). The original data were in 15 min time intervals, and were averaged to mean daily discharge time series. 3. Pump House discharge data as mean daily discharge (downloaded from the SFA portal). 4. BCC and Pump House chemistry data from SFA data portal and/or original spreadsheets provided by Roelof Versteeg. The following challenging QA/QC problems of the datasets were resolved: Missing data with the duration of gaps up to >1 month; Duplicated dates; Anomalies and outliers of discharge and concentrations; Time stamps of measurements of the discharge and concentrations are not aligned (hydrogeochemical balance calculations require the timestamps to be aligned). All QA/QC-ed datasets are given as csv files. The csv files were prepared using the xts files with multiple worksheets, which are also included in the data packages. Figures of the QA/QC-ed datasets are given in the jpeg and pdf formats. The QA/QC-ed datasets have been used to quantify discharge and chemical concentrations in river water in order to understand riverine exports of water and dissolved constituents in the East River watershed. These datasets served as a basis in the presentation given by P. Fox et al. at the 2021 Goldschmidt Conference.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Sensitivity of climate and atmospheric CO2 to deep-ocean and shallow-ocean carbonate burial

A model of the carbonate-silicate geochemical cycle is presented that distinguishes carbonate masses produced by shallow-ocean and deep-ocean carbonate burial and shows that reasonable increases in deep-ocean burial could produce substantial warmings over a few hundred million years. The model includes exchanges between crust and mantle; transients from burial shifts are found to be sensitive to the fraction of nondegassed carbonates subducted into the mantle. Without the habitation of the open ocean by plankton such as foraminifera and coccolithophores, today's climate would be substantially colder.

Volk, Tyler↗
NASA NTRS2002

SNC Meteorites and Martian Reservoirs

Jones first suggested that the inverse covariation of initial epsilon (Nd-143) and Sr-87/Sr-86 of the shergottites could be explained by interaction between mantle-derived magmas with another isotopic reservoir(s) (i.e., assimilation or contamination). In that model, magmas were generated in a source region that was isotopically similar to the Nakhla source and the second reservoir(s) was presumed to be crust. The text also permitted the second reservoir to be another type of mantle, but I can confirm that a second mantle reservoir was never seriously considered by that author. Other features of this model were that (i) it occurred at a particular time, 180 m.y. ago, and (ii) the interacting reservoirs had been separated at approximately 4.5 b.y. In a later paper Jones explored this mixing model more quantitatively and concluded that magmas from a Nakhla-like source region at 180 m.y. would fall on or near an isotopic Nd-Sr-Pb hyperplane defined by the shergottites. This criterion was a necessary prerequisite for the parent magma(s) of the shergottites to have initially been Nakhla-like isotopically. At this juncture, it is perhaps worthwhile to note that this mixing model was not presented to explain geochemical variations but as a justification for a 180 m.y. crystallization age for the shergottites and a 1.3 b.y. crystallization age for the nakhlites. In the mid-1980's crystallization ages estimated for Nakhla ranged from approximately 1.3 b.y to 4.5 b.y. Similarly, preferred crystallization ages for the shergottites ranged from 360 m.y., to 1.3 b.y., to 4.5 b.y. In all these models, the 180 m.y. event seen in the shergottites was deemed to be metamorphic. The fit between the Nakhla-like source region and the shergottite hyperplane was a validation both of the 1.3 b.y. igneous age of Nakhla and the 180 m.y. igneous age of the shergottites.

Jones, J. H.↗
DOE OSTI2021

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

The Role of KREEP in the Production of Mg-Suite Magmas and Its Influence on the Extent of Mg-Suite Magmatism in the Lunar Crust

The lunar magnesian-suite, or Mg-suite, is a series of ancient plutonic rocks from the lunar crust. They have received a considerable amount of attention from lunar scientists since their discovery for three primary reasons: 1) their ages and geochemistry indicate they represent pristine magmatic samples that crystallized very soon after the formation of the Moon; 2) their ages often overlap with ages of the ferroan anorthosite (FAN) crust; and 3) planetary-scale processes are needed in formation models to account for their unique geochemical features. Taken as a whole, the Mg-suite samples, as magmatic cumulate rocks, approximate a fractional crystallization sequence in the low-pressure forsterite-anorthite-silica system, and thus these samples are generally thought to be derived from layered mafic intrusions which crystallized very slowly from magmas that intruded the anorthositic crust. However, no direct linkages have been established between different Mg-suite samples based either on field relationships or geochemistry.The model for the origin of the Mg-suite, which best fits the limited available data, is one where Mg-suite magmas form from melting of a hybrid cumulate package consisting of deep mantle dunite, crustal anorthosite, and KREEP (potassium-rare earth elements-phosphorus) at the base of the crust under the Procellarum KREEP Terrane (PKT). In this model, these three LMO (Lunar Magma Ocean) cumulate components are brought into close proximity by the cumulate overturn process. Deep mantle dunitic cumulates with an Mg number of approximately 90 rise to the base of the anorthositic crust due to their buoyancy relative to colder, more dense Fe- and Ti-rich cumulates. This hybridized source rock melts to form Mg-suite magmas, saturated in Mg-rich olivine and anorthitic plagioclase, that have a substantial KREEP component.

Elardo, S. M.↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2023

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Rock-Sourced Nitrogen in Semi-Arid, Shale-Derived California Soils

Models suggest that rock-derived nitrogen (N) inputs are of global importance to ecosystem N budgets; however, field studies demonstrating the significance of rock N inputs are rare. We examined rock-derived N fluxes in soils derived from sedimentary rocks along a catena formed under a semi-arid climate. Our measurements demonstrate that there are distinct and traceable pools of N in the soil and bedrock and that the fraction of rock-derived N declines downslope along the catena. We used geochemical mass balance weathering flux measurements to estimate a rock-derived N flux of 0.145 to 0.896 kg ha -1 yr -1 at the ridgecrest. We also developed independent N flux estimates using a 15 N-based isotope mixing model. While geochemical mass-balance-based estimates fell within the 95% confidence range derived from the isotope mixing model (-1.1 to 44.3 kg ha -1 yr -1 ), this range was large due to uncertainty in values for atmospheric 15 N deposition. Along the catena, N isotopes suggest a diminishing effect of rock-derived N downslope. Overall, we found that despite relatively large N pools within the saprolite and bedrock, slow chemical weathering and landscape denudation limit the influence of rock-derived N, letting atmospheric N deposition (7.1 kg ha -1 yr -1 ) and N fixation (0.9-3.1 kg ha -1 yr -1 ) dominate N inputs to this grassland ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Utah FORGE 5-2565: Hydrothermal Evolution of Fracture Properties - Workshop Presentation

This is a presentation on the Evolution of Permeability and Strength Recovery of Shear Fractures Under Hydrothermal Conditions project by the U.S. Geological Survey, presented by Dr. David Lockner. The project's objective was to determine how thermal, hydraulic, mechanical, and chemical processes affect the sustainability of fracture networks in geothermal reservoirs and provide strategies for improved EGS techniques that maximize thermal coupling and increase reservoir longevity. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗