Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

In situ aerobic methane oxidation rates in a stratified lake

Abstract Microbial aerobic methane oxidation is an important sink for aquatic methane worldwide. Despite its importance to global methane fluxes, few aerobic methane oxidation rates have been obtained in freshwater or marine environments without imposing changes to the microbial community through use of ex situ methods. A novel in situ incubation method for continuous time‐series measurements was used in Jordan Lake, North Carolina, during 2020–2021, to determine reaction kinetics for aerobic methane oxidation rates across a wide range of naturally varying methane (55–1833 nM) and dissolved oxygen (DO; 28–366 μM) concentrations and temperatures (17–30°C). Methane oxidation began immediately at the start of each of 21 incubations and methane oxidation rates were 1 st order with respect to methane. The data density allowed for accurate calculation of 1 st ‐order rate constants,k, that ranged from 0.018 to 0.462 h −1 (R 2 > 0.967). Addition of ammonium (20–45 μM) to natural concentrations ranging from 0.057 to 2.4 μM did not change aerobic methane oxidation rate kinetics, suggesting that the natural population of aerobic methane oxidizers in this eutrophic lake was not nitrogen limited. Values ofkinversely correlated most strongly with initial DO concentrations (R 2 = 0.82) rather than temperature. Values forkincreased with Julian day throughout our sampling period, suggesting seasonal influences on methane oxidation via responses to geochemical changes or shifts in microbial community abundance and composition. These experiments demonstrate a high variability in the enzymatic capacity for 1 st ‐order methane oxidation rates in this eutrophic lake that is tightly and inversely coupled to oxygen concentrations. Measurements of in situ aerobic methane oxidation rate constants allow for the direct quantification and modeling of the microbial community's capacity for methane oxidation over a wide range of natural methane concentrations.

Marine & Freshwater Biology↗
DOE OSTI2025

Linking the Spin Transition of Ferric Iron in δ‐(Al,Fe)OOH to Water Storage in the Lower Mantle

As the most massive geochemical reservoir, the lower mantle affects the Earth's budget of volatile elements, including hydrogen or H 2 O. The properties of minerals in the lower mantle are further affected by changes in the electronic configurations of iron cations, that is, by spin transitions. The feedback between spin transitions and potential storage of H 2 O in solid hydrous phases in the lower mantle, however, remains unexplored. By combining high‐pressure nuclear resonant inelastic X‐ray scattering and high‐pressure high‐temperature X‐ray diffraction experiments, we constrained the thermal equation of state of δ‐(Al,Fe)OOH, a member of the phase H solid solution. Based on the derived thermal equation of state of δ‐(Al,Fe)OOH and the underlying thermodynamic model, we calculate the excess Gibbs free energy that arises from the spin transition of ferric iron in this compound and evaluate the effect on phase equilibria. The results of our analysis show that the spin transition of ferric iron in phase H may significantly reduce the thermodynamic activity and hence the concentration of H 2 O in a coexisting hydrous melt. As a consequence, nominally anhydrous minerals of the lower mantle may become dehydrated in the presence of phase H. Our analysis further suggests that, under certain conditions, the spin transition may expand the thermal stability of Fe 3+ ‐bearing phase H and create a geochemical link between the storage of H 2 O in phase H and ferric iron in the lower mantle.

Buchen, Johannes [Univ. of Bayreuth (Germany); Cal↗
DOE OSTI2020

Geochemical narrowing of cement fracture aperture during multiphase flow of supercritical CO2 and brine

For carbon capture and storage operations, wells are used as a necessary conduit for injecting CO2 at depth, but they can also act as leakage conduits if the cement seals are compromised. The specific objective of this research was to investigate the nature of self-healing of a fracture within cement under multiphase flow of CO2 and brine, and to compare this to a newly developed model. This was accomplished through a laboratory investigation, flowing reacting supercritical CO2 and brine through a cement fracture with multiphase flow. The influent end of the fracture was degraded as evidenced by the formation of cracks across the surface. At the effluent end of the fracture, the initial aperture of 137 µm was reduced to 50 µm. This reduction by 87 µm was similar to results from a recently developed model tested in this study, which predicted an aperture reduction of 80 µm. Our results demonstrated that the model accurately predicted the reactions which occurred in this multiphase flow system.

wellbore, cement, fracture, multiphase-flow, super↗
DOE OSTI2021

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗
DOE OSTI2025

Elasticity of davemaoite as a primary contributor to lower-mantle heterogeneities

Geophysical detection of subducted mid–ocean ridge basalt (MORB) in the lower mantle is hindered by uncertainties in the elasticity of Fe,Al,Mg,Ti–bearing davemaoite, a key MORB component. Using Brillouin spectroscopy and x-ray diffraction, we determined the elasticity of a Ca 0.906(1) Fe 2+ 0.027(1) Fe 3+ 0.042(1) Mg 0.033(1) Al 0.072(1) Ti 0.020(1) Si 0.912(1) O 3 davemaoite up to 113 gigapascals and 2294 K. Here, we found that it exhibited a shear wave velocity 10 to 20% slower than end-member davemaoite, making it the slowest phase among major lower-mantle minerals. Our models show that MORB, containing 20 to 25 volume percent davemaoite, potentially contributes to large low-shear-velocity provinces (LLSVPs), whereas a cumulate layer enriched in davemaoite crystallized from basal magma ocean may comprise ultralow-velocity zones (ULVZs). Davemaoite’s ability to host incompatible and heat-producing elements possibly links LLSVPs and ULVZs to mantle plume initiation and geochemical signatures of ocean island basalts.

Zhou, Wen-Yi [Texas A & M Univ., College Station, ↗
DOE OSTI2025

Assessing the suitability of sites near Pine Island Glacier for subglacial bedrock drilling aimed at detecting Holocene retreat–readvance

Abstract. Unambiguous identification of past episodes of ice sheet thinning below the modern surface and grounding line retreat inboard of present requires recovery and exposure dating of subglacial bedrock. Such efforts are needed to understand the significance and potential future reversibility of ongoing and projected change in Antarctica. Here we evaluate the suitability for subglacial bedrock drilling of sites in the Hudson Mountains, which are located in the Amundsen Sea sector of West Antarctica. We use an ice sheet model and field data – geological observations, glaciological observations and bedrock samples from nunataks, and ground-penetrating radar from subglacial ridges – to rate each site against four key criteria: (i) presence of ridges extending below the ice sheet, (ii) likelihood of increased exposure of those ridges if the grounding line was inboard of present, (iii) suitability of bedrock for drilling and geochemical analysis, and (iv) accessibility for aircraft and drilling operations. Our results demonstrate that although no site in the Hudson Mountains is perfect for this study when assessed against all criteria, the accessibility, N–S orientation and basaltic bedrock lithology of Winkie Nunatak's southernmost ridge (74.86° S, 99.77° W) make it a feasible site both for drilling and subsequent cosmogenic nuclide analysis. Furthermore, the ridge is strewn with glacial erratics at all elevations, providing valuable constraints on its early Holocene deglacial history. Based on our experiences during this study, we conclude with a series of recommendations for assessing site suitability for future bedrock drilling campaigns. We emphasise the importance of consulting a range of expertise prior to drilling and ensuring that sufficient field reconnaissance is undertaken (including obtaining detailed grids of radar survey data and bedrock samples).

58 GEOSCIENCES↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Can Glacial Sea‐Level Drop‐Induced Gas Hydrate Dissociation Cause Submarine Landslides?

Abstract We conducted two‐dimensional numerical simulations to investigate the mechanisms underlying the strong spatiotemporal correlation observed between submarine landslides and gas hydrate dissociation due to glacial sea‐level drops. Our results suggest that potential plastic deformation or slip could occur at localized and small scales in the shallow‐water portion of the gas hydrate stability zone (GHSZ). This shallow‐water portion of the GHSZ typically lies within the area enclosed by three points: the BGHSZ–seafloor intersection, the seafloor at ∼600 m below sea level (mbsl), and the base of the GHSZ (BGHSZ) at ∼1,050 mbsl in low‐latitude regions. The deep BGHSZ (>1,050 mbsl) could not slip; therefore, the entire BGHSZ was not a complete slip surface. Glacial hydrate dissociation alone is unlikely to cause large‐scale submarine landslides. Observed deep‐water (much greater than 600 mbsl) turbidites containing geochemical evidence of glacial hydrate dissociation potentially formed from erosion or detachment in the GHSZ pinch‐out zone.

58 GEOSCIENCES↗
DOE OSTI2025

High-Temperature Aquifer Thermal Energy Storage (HT-ATES) Projects in Germany and the Netherlands—Review and Lessons Learned

Aquifer thermal energy storage (ATES) is a concept that can help to address heating and cooling needs through the use of the subsurface as a seasonal thermal energy storage (STES) system. Over 2800 ATES systems have been deployed with storage temperatures typically below 25 °C and only a few with higher temperatures (>40 °C), which would increase the energy density and utility of the stored thermal fluids. Until now, only a few high-temperature aquifer thermal energy storage (HT-ATES) projects have been initiated and are still in operation. These HT-ATES projects have encountered a range of technical and non-technical challenges. This study reviews ten such projects: four in Germany and six in the Netherlands. The non-technical issues include public acceptance, a lack of regulatory framework for these systems, managing overlapping uses of the subsurface, managing changes with the providers and off-takers of thermal energy, and obtaining financing to implement these projects. Common technical issues include geological factors such as incomplete characterization of the subsurface and reservoir heterogeneity; geochemical issues such as mineral scaling, corrosion, and biofouling; lower than expected thermal recovery; and issues with system design and reliability. This review highlights benefits and challenges faced by HT-ATES projects with the goal to use the lessons learned to improve the siting, design, development, and operation of such systems. Recommendations include improved initial subsurface site characterization, use of coupled process models to optimize system design and predict system performance, cascaded uses of stored thermal energy to better utilize the stored heat, monitoring networks to provide feedback on system performance, and expanded system scale to allow for continued operation even when maintenance of some system components is required. Techno-economic modeling and risk analysis could be used to optimize such HT-ATES project design and identify key factors that will affect sustained economic viability. In addition, design flexibility is important for these systems to allow for changing conditions regarding the supply and demand of thermal energy. Adopting these findings should improve the performance and reduce the risks for future HT-ATES projects worldwide.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2021

The Colorado East River Community Observatory Data Collection

Abstract The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow‐dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long‐term and spatially‐extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross‐disciplinary process‐based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value‐added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Influence of clay mineral weathering on green rust formation at iron-reducing conditions

Green rusts (GR) are important drivers for trace metal and nutrient cycling in suboxic environments. We investigated whether green rusts would incorporate aluminum (Al) or other elements from naturally-formed clay minerals containing easily-weatherable clay minerals (e.g. mica, interlayered clays). We isolated the clay minerals from a Matapeake silt loam soil by removal of silt and sand, organic matter, and reducible oxides to study mechanisms of interaction between Fe(II) and soil-sourced clay minerals. We conducted batch Fe(II) sorption experiments at multiple near-neutral pHs (6.5–7.5) and reaction times (2h–365 days). Mineral transformations were characterized by selective extractions, X-ray diffraction (XRD), and Fe X-ray absorption spectroscopy (XAS) analyzed by shell-fitting and linear combination fitting (LCF) with natural and synthetic standards. Clay mineral fraction contained a mixture of quartz, kaolinite, interlayered vermiculite, mica, and chlorite with significant structural Fe (2.6% wt). Uptake of Fe(II) increased with pH and kinetics were rapid until 5 days, followed by slow continuous Fe(II) uptake. Citrate-bicarbonate desorption kinetics from Fe(II) sorbed clay released more Al and silicon (Si) compared with unreacted soil clay fraction whereas magnesium (Mg) and potassium (K) were unaffected. Citrate-bicarbonate extracted Fe contained more Fe(II) than an ideal GR with an Fe(II)/Fe(III) molar ratio of 5.50. Analysis of the Fe EXAFS by both LCF and shell fitting was best modeled as a combination of Fe(III)-clay reduction to Fe(II) and precipitation of GR and Fe(II)-Al LDH. After 7 days of Fe(II) sorption, LCF identified 55.2% total Fe in clay, 33.4% GR(Cl) and 11.4% Fe(II)-Al LDH. Further, these results provide novel evidence of Fe(II)-Al LDHs precipitating on naturally-formed soil clay minerals as a minor phase to GR. The geochemical implications are that GRs formed in soils and sediments should be considered to have Al and Si as well as Mg substitutions affecting their structure and reactivity.

58 GEOSCIENCES↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

North‐East Peri‐Tethyan Water Column Deoxygenation and Euxinia at the Paleocene Eocene Thermal Maximum

Abstract The Paleocene–Eocene Thermal Maximum (PETM) is associated with climatic change and biological turnover. It shares features with the Oceanic Anoxic Events (OAEs) of the Mesozoic, such as transient global warming and biogeochemical perturbations. However, the PETM experienced a more muted expansion of marine anoxia compared to the Mesozoic OAEs (especially OAE2), with geographically limited evidence for photic zone euxinia (PZE). We explore the extent and drivers of marine deoxygenation during the PETM using biomarkers for water column euxinia and anoxia as well as an intermediate complexity Earth system model (cGEnIE). These reveal that the water column in the North‐East Peri‐Tethys became anoxic, with euxinic conditions reaching the photic zone (PZE) during the PETM. Our model shows that euxinia developed due to a global increase in the ocean nutrient inventory with concomitant oxygen consumption, similar to findings for OAE2. The particularly strong regional response in the NE Peri‐Tethys appears to arise from a combination of global CO 2 ‐weathering forcing, regionally restricted circulation and upwelling of sulphidic thermocline waters. Unlike OAE2, anoxia and PZE do not become widespread in our PETM simulations, consistent with new and existing geochemical and biological proxy data. This globally muted response could result from reduced biogeochemical feedbacks to climate forcing relative to the mid‐Cretaceous climate. Our observations suggest that similar mechanisms operated in response to disparate Cenozoic (PETM) and Mesozoic (OAEs) transient global warming events, while also highlighting that background conditions are crucial in modulating the sensitivity of Earth's system to them.

Behrooz, L.↗
DOE OSTI2021

Modeling and Analysis of the Transport and Disposal of Beryllium Moderator Blocks and Greater than Class C (GTCC) Waste

Radioactive waste in the United states is categorized based on its source, radioactivity, and security risk. The categorization method employed by the Nuclear Regulatory Committee (NRC) for low-level radioactive waste is not currently applied to any waste generated by the Department of Energy (DOE) or disposed of in DOE facilities. As a result, there is a large volume of DOE-generated waste with characteristics similar to the NRC’s “Greater than Class C” (GTCC) waste category which presently do not have a path for long-term disposal. Much of this waste cannot undergo the standard commercial disposition process due to it being categorized as Transuranic (TRU) waste under DOE guidelines, due to elevated concentrations of key fission products. Possible solutions to this dilemma are explored in the Environmental Impact Statement for the “Disposal of Greater-Than-Class-C Low Level Radioactive Wave” (EIS-0375). This paper analyzes several EIS possible paths for disposal for this “orphaned” waste. One example that highlights this categorization issue is the spent beryllium cladding that has been extracted from the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) throughout its operational lifetime. These beryllium blocks surround the reactor’s main chambers and components and serve as neutron moderators. This paper will use this particular waste form to analyze the economic and technological viability of transporting and storing this material under the storage and transportation criteria of the Waste isolation Pilot Plant (WIPP), which was deemed the most feasible disposition path according to the EIS. The example of this waste form used in this paper that highlights this categorization issue is the spent beryllium cladding extracted from Idaho National Laboratory’s (INL) Advanced Test Reactor (ATR). These beryllium blocks surround the reactor’s components and serve as neutron moderators, and thus have accumulated high concentration of high-activity transuranic isotopes. This paper analyzed the viability of transporting and storing this material under the storage and transportation criteria of the five primary disposition paths covered in EIS-0375. An array of calculations was carried out to assess the viability of the various disposition paths for the beryllium shipments as well as other waste shipments that fall within the GTCC category. Modeling with MCNP 6.2 was conducted to determine whether the beryllium blocks could be safely stored and transported within a 72-B cask; the standard shipping container for remote-handled, transuranic waste, while also meeting regulatory limits at various disposition paths. A cost analysis of the transportation and long-term disposal paths mentioned in the EIS was also carried out using available data and information from similar waste shipments. Lastly, a geochemical analysis of the various geological repository discussed in the EIS was also carried out using the Geochemists Workbench Release 14. Long-term disposal of the beryllium blocks at the Waste Isolation Pilot Plant (WIPP) proved to be the most cost-effective long-term disposition option of the ones considered in the EIS. A dose rate calculation at both contact and remote-handling distances indicate that the analyzed beryllium shipments should not exceed the exposure limits at any of the considered locations. Greater-than-class-C waste has been left in a regulatory state of limbo for years, resulting in backlogs of inventory across several research sites in the United States. Due to its high activity and presence of transuranic isotopes, it is imperative to ensure that it remains inaccessible and sequestered both in its short-term interim as well as a long-term geologic time scale. The information and assessments done in the paper could potentially serve as a reference for any future shipments of this waste form at the WIPP facility as well as other disposition paths that may be considered in the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Tracking natural CO 2 migration through a sandstone aquifer using Sr, U and C isotopes: Chimayó, New Mexico, USA

The geochemical and isotopic characteristics of groundwaters in Chimayó, New Mexico, reflect processes that affect water quality in the Tesuque Aquifer, which overlies a leaking natural CO 2 source in a structurally complex region. In this study, select isotopes (δ 13 C, 87 Sr/ 86 Sr, 234 U/ 238 U) are applied to groundwaters to better understand CO 2 transport mechanisms and related water-rock interactions that impact Chimayó groundwater. Carbon stable isotope ratios of dissolved inorganic carbon (DIC; δ 13 CDIC = -15.10‰ to 4.50‰) identify a distinct source of upward-migrating CO 2 that interacts with the groundwater. Additionally, groundwater 87 Sr/ 86 Sr compositions (0.7098 to 0.7154) reflect intrusion of varying amounts of saline water associated with the high CO 2 source, while 234 U/ 238 U ratios corroborate presence of deep groundwater and suggest impacts from distinctive natural uranium sources are affecting groundwater. Previous work proposed two CO 2 transport mechanisms at the site: (1) dissolved in deep brine that underlies the aquifer and (2) in the gas phase; this study uses isotope mixing models to identify wells that are affected by these CO 2 transport mechanisms and demonstrates that both transport mechanisms are associated with impaired groundwaters. Overall, this study demonstrates that applying multiple isotope systems (δ 13 CDIC, 87 Sr/ 86 Sr, 234 U/ 238 U) is a dynamic tool for identifying and measuring the impact of CO 2 leakage from a sequestration site.

58 GEOSCIENCES↗