Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic model”

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 883 records · Page 49

Factors influencing organic carbon preservation in marine sediments

The organic matter that escapes decomposition is buried and preserved in marine sediments, with much debate as to whether the amount depends on bottom-water O2 concentration. One group argues that decomposition is more efficient with O2, and hence, organic carbon will be preferentially oxidized in its presence, and preserved in its absence. Another group argues that the kinetics of organic matter decomposition are similar in the presence and absence of O2, and there should be no influence of O2 on preservation. A compilation of carbon preservation shows that both groups are right, depending on the circumstances of deposition. At high rates of deposition, such as near continental margins, little difference in preservation is found with varying bottom-water O2. It is important that most carbon in these sediments decomposes by anaerobic pathways regardless of bottom-water O2. Hence, little influence of bottom-water O2 on preservation would, in fact, be expected. As sedimentation rate drops, sediments deposited under oxygenated bottom water become progressively more aerobic, while euxinic sediments remain anaerobic. Under these circumstances, the relative efficiencies of aerobic and anaerobic decomposition could affect preservation. Indeed, enhanced preservation is observed in low-O2 and euxinic environments. To explore in detail the factors contributing to this enhanced carbon preservation, aspects of the biochemistries of the aerobic and anaerobic process are reviewed. Other potential influences on preservation are also explored. Finally, a new model for organic carbon decomposition, the "pseudo-G" model, is developed. This model couples the degradation of refractory organic matter to the overall metabolic activity of the sediment, and has consequences for carbon preservation due to the mixing together of labile and refractory organic matter by bioturbation.

Non-NASA Center↗

Identifying geologic characteristics and operational decisions to meet global carbon sequestration goals

Geologic carbon sequestration is the process of injecting and storing CO 2 in subsurface reservoirs and is an essential technology for global environmental security (e.g., climate change mitigation) and economic security (e.g., CO 2 tax credits). To meet energy, economic, and environmental goals, society will have to identify vast volumes of high-capacity, low-cost, and viable storage reservoirs for sequestering CO 2 . In turn, this requires understanding how major geologic characteristics (such as reservoir depth, thickness, permeability, porosity, and temperature) and design and operational decisions (such as injection well spacing) impact CO 2 injection rates, storage capacity, and economics. Although many numerical simulation tools exist, they cannot repeat the required thousands or millions of simulations to identify ideal reservoir properties and the sensitivity and interaction between geologic parameters and operational decisions. Here, we use SCO 2 T—a fast-running, reduced-order modeling framework—to explore the sensitivity of major geologic parameters and operational decisions to engineering (CO 2 injection rates, plume dimensions, and storage capacities and effectiveness) and costs. Our results show, for the first time, benefits and impacts such as allowing CO 2 plumes to overlap, how different well spacing patterns affect CO 2 sequestration, the effects on costs of including brine treatment and disposal, and the effect of restricting injection rates to 1 MtCO 2 per y based on well limitations. We reveal multiple novel and unintuitive findings including: (i) deeper reservoirs have reduced carbon sequestration costs until injection rates reach 1 MtCO 2 per y, at which point deeper reservoirs become more expensive, (ii) thicker formations allow for increased injection rates and storage capacity, but thickness barely impacts plume areas, (iii) higher geothermal gradients result in reduced sequestration costs, unless brine treatment/disposal costs are included, at which point reservoirs having lower geothermal gradients are more economical because they produce less brine for each unit of injected CO 2 , and (iv) allowing plumes to overlap has a significantly positive impact of increasing storage capacities but has only a small influence on reducing sequestration costs. Altogether, our results illustrate new scientific conclusions to help identify suitable sites to inject and store CO 2 , to help understand the complex interaction between geology and resulting costs, and to help support the pursuit of meeting global sequestration targets.

58 GEOSCIENCES↗

Geology Before Pluto: Pre-Encounter Considerations

Pluto, its large satellite Charon, and its four known satellites represent the first trans-Neptunian Kuiper Belt objects populating the outer-most solar system beyond the gas giant planets to be studied in detail from a spacecraft (New Horizons). A complete picture of the solar nebula, and solar system formation cannot be confidently formulated until representatives of this group of bodies at the edge of solar space have been examined. The Pluto system is composed of unique lunar- and intermediate-sized objects that can tell us much about how objects with volatile icy compositions evolve. Modeling of the interior suggests that geologic activity may have been to some degree, and observations of frost on the surface could imply the need for a geologic reservoir for the replenishment of these phases. However, the putative indicators of Pluto's geologic history are inconclusive and unspecific. Detailed examination of Pluto's geologic record is the only plausible means of bridging the gap between theory and observations. In this talk I will examine the potential importance of these tentative indications of geologic activity and how specific spacecraft observations have been designed and used to constrain the Pluto system's geologic history. The cameras of New Horizons will provide robust data sets that should be immanently amenable to geological analysis of the Pluto System's landscapes. In this talk, we begin with a brief discussion of the planned observations by New Horizons' cameras that will bear most directly on geological interpretability. Then I will broadly review major geological processes that could potentially operate of the surfaces of Pluto and its moons. I will first survey exogenic processes (i.e., those for which energy for surface modification is supplied externally to the planetary surface): impact cratering, sedimentary processes (including volatile migration) and the work of wind. I will conclude with an assessment of prospects for endogenic activity in the form of tectonics and cryo-volcanism.

Pluto↗

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. We developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. We tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.

58 GEOSCIENCES↗

Coldspots and hotspots - Global tectonics and mantle dynamics of Venus

Based on geologic observations provided by Magellan's first cycle of data collection and recent models of mantle convection in spherical shells and crustal deformation, the major topographic and geologic features of Venus are incorporated into a model of global mantle dynamics. Consideration is given to volcanic rises, such as Beta Regio and Atla Regio, plateau-shaped highlands dominated by complex ridged terrain (e.g., Ovda Regio and Alpha Regio), and circular lowland regions, such as Atalanta Planitia. Each of these features is related to either mantle plumes (hotspots) or mantle downwellings (coldspots).

Bindschadler, Duane L.↗

The Role of Chemistry in Fracture Pattern Development and Opportunities to Advance Interpretations of Geological Materials

Abstract Fracture pattern development has been a challenging area of research in the Earth sciences for more than 100 years. Much has been learned about the spatial and temporal complexity inherent to these systems, but severe challenges remain. Future advances will require new approaches. Chemical processes play a larger role in opening‐mode fracture pattern development than has hitherto been appreciated. This review examines relationships between mechanical and geochemical processes that influence the fracture patterns recorded in natural settings. For fractures formed in diagenetic settings (~50 to 200 °C), we review evidence of chemical reactions in fractures and show how a chemical perspective helps solve problems in fracture analysis. We also outline impediments to subsurface pattern measurement and interpretation, assess implications of discoveries in fracture history reconstruction for process‐based models, review models of fracture cementation and chemically assisted fracture growth, and discuss promising paths for future work. To accurately predict the mechanical and fluid flow properties of fracture systems, a processes‐based approach is needed. Progress is possible using observational, experimental, and modeling approaches that view fracture patterns and properties as the result of coupled mechanical and chemical processes. A critical area is reconstructing patterns through time. Such data sets are essential for developing and testing predictive models. Other topics that need work include models of crystal growth and dissolution rates under geological conditions, cement mechanical effects, and subcritical crack propagation. Advances in machine learning and 3‐D imaging present opportunities for a mechanistic understanding of fracture formation and development, enabling prediction of spatial and temporal complexity over geologic timescales. Geophysical research with a chemical perspective is needed to correctly identify and interpret fractures from geophysical measurements during site characterization and monitoring of subsurface engineering activities.

58 GEOSCIENCES↗

Preliminary Geothermal Resource Assessment of the St. Helens Seismic Zone Using the Results from the Geothermal Play-Fairway Analysis of Washington State Prospects

The northern St. Helens Seismic Zone was one of the areas of interest during Phase 2 of the Washington State Play-Fairway analysis. The area is assessed for geothermal resource potential using the volumetric reserve estimation method developed by the US Geological Survey. Multiple hypotheses for reservoir volume are proposed and analyzed based on interpretation of data gathered during Phase 2 within a geothermal lease. Here, the data utilized in the analysis includes new geologic mapping, ground-based gravity survey, magnetotelluric survey, and passive seismic survey. The hypotheses will be tested, and the resource model will be updated after temperature gradient holes are drilled in the area during Phase 3 at sites selected from Phase 2. The analysis also makes use of the prior geologic knowledge of the area, including the increased seismic and geodetic monitoring following the eruption of Mount St. Helens on May 18th, 1980, existing geologic mapping and models developed from seismic and geodetic data, and past geothermal prospecting efforts by the Washington Division of Geology and Earth Resources. The analysis uses reservoir temperature from geochemistry and a single temperature gradient hole for a low value, and a temperature required for steam flash as a high value. Subsurface temperatures, and heat flows are expected to be better constrained after Phase 3 drilling. The preliminary results suggest power generation potential ranging from 65 MWe to 1.5 GWe within the lease area.

geomechanics↗

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

The origin of fluvial valleys and early geologic history, Aeolis Quadrangle, Mars

A valley genesis model of the southern Aeolis Quadrangle, eastern Mars, is proposed. In response to widespread effusive volcanism in interstratified, ice-reach terrains, local subsidence occurred along fractures/faults, and produced scarps that intersected local aquifers. Heated spring discharges issuing from these scarps may have carved the valleys through a combination of downvalley water, ice, and debris flows, and headward erosion along stratal discontinuities and individual conduit faults and fractures. Fifty-six asymmetric scarps or ridges that are probable thrust faults are mapped. The faults exhibit an orientation vector mean of N 63 deg W +/- 11 deg and they transect the lava plains and the older plateau sequence units. The vector mean for the 264 valleys mapped is N 48 deg W +/- 12 deg, with a larger dispersion about the mean. The similar orientations displayed by thrust faults and valley axes suggest that valley locations are partly controlled by preexisting thrust faults and fracture systems.

Brakenridge, G. Robert↗

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗

Utilizing Earth Observations to Model Probable Coastal Wetland Extent, Sea-Level Rise Inundation Risk, and Assess Impacts on Historic Hawaiian Lands

Climate induced sea-level rise poses a risk to coastal areas on the Island of Hawai’i, and many of the island’s historic cultural lands are in danger of becoming overtaken by wetlands or inundation. In partnership with the County of Hawai’i, State of Hawai’i Department of Land and Natural Resources, and Arizona State University, NASA DEVELOP mapped wetland extent and short-term sea-level rise inundation risk. We utilized Earth observations over a 10-year span (2013 – 2022) that included the NASA MEaSUREs Gridded Sea Surface Height Anomalies and MEaSUREs Group for High Resolution Sea Surface Temperature datasets, United States Geological Survey (USGS) Hawaii Digital Elevation Models (DEM), and in situ tidal gauge data. Flood risk index values were acquired for 5 known Hawai’i flood events between 2019 – 2021 from the Global Flood Mapper tool on Google Earth Engine. We used a random forest model to predict short-term sea-level rise inundation risk along the entire coast of Hawai’i. Current wetland extents and probabilistic locations of new wetlands were modeled with the most recently available data from PlanetScope Surface Reflectance optical imagery (2022), USGS 3D Elevation Program (3DEP) 10m DEM (2020), temperature and precipitation data from the Hawai’i Climate Atlas, and soils data from the Hawai’i Soil Atlas (2014) using the Wetland Intrinsic Potential tool. Results indicated locations that had the highest probability of wetland creation. The end products aimed to help the partners prioritize efforts to meeting regulation requirements for wetlands protection, evaluate the inundation risk to historical features, and support decision-making for their Shoreline Setback and Climate Adaption plans.

Lisa Tanh↗

The deep structure of Venusian plateau highlands

Magellan gravity data confirm that several of the large, tectonically deformed, plateau-like highlands on venus are shallowly compensated, most likely by crustal thickness variations. Apparent depths of isostatic compensation, computed in the spatial domain, range from 30 to 50 km for Alpha, Tellus, Ovda, and Thetis Regiones. Using a two-layer model for isostatic compensation, Alpha, Tellus, and Ovda are best represented as nearly completely compensated in crust that is regionally 20-40 km thick around these highlands, with little contribution from deeper mantle sources. In contrast to these three areas, a stronger regional gravity high associated with Thetis requires a significant upper mantle component to compensation. This is evident in the spectral admittance as a pronounced deep, long-wavelength anomaly. In the two-layer isostatic model, a broad, deeply compensated upland underlies a shallowly compensated central block of Thetis. If this deep component is interpreted as a thermal anomaly, the loci of maximum upwelling agree well with sites of recent extension. The plateau highlands are thus physiographically and isostatically equivalent to terrestrial continents, though probably not compositionally. They also share the record of a long tectonic history. The large regional gravity anomaly of Thetis indicates that active mantle proceses continue even beneath some areas (tessera) thought to be a relic of a former geological regime. The excellent agreement of modeled crustal thicknesses around Alpha, Tellus, and Ovda Regiones suggests that 20-40 km is a representative global value for the plains. Such a crust is thicker than previously estimated and about twice as thick as the expected thickness of crust produced at venusian spreading centers

Grimm, Robert E.↗

Discovering Blind Geothermal Systems in the Great Basin Region: An Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways: All Phases

Most geothermal resources in the Great Basin region of the western USA are blind, and thus the discovery of new commercial-grade systems requires synthesis of favorable characteristics for geothermal activity. The geothermal play fairway concept involves integration of multiple parameters indicative of geothermal activity to identify promising areas for new development. This project integrated multiple datasets to apply the play fairway concept and assess geothermal potential in a large region of the Great Basin in Nevada. It is therefore referred to as the Nevada play fairway project. This project was a strong collaborative effort between several organizations, led by the Nevada Bureau of Mines and Geology at the University of Nevada, Reno, but with key support from the U.S. Geological Survey, ATLAS Geosciences, Inc,, Hi-Q Geophysical, Inc., Lawrence Berkeley National Laboratory, Utah Geological Survey, and Innovative Geothermal Ltd. In Budget Period 1 of this project, available data for nine geologic, geochemical, and geophysical parameters were initially synthesized to produce a new detailed geothermal potential map of 96,000 km2 from west-central to eastern Nevada. These parameters were grouped into subsets and individually weighted to delineate rankings for local permeability, intermediate permeability, regional permeability, and thermal potential, which collectively defined geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). This initial work was aimed at reducing the risks in regional exploration and therefore facilitating discovery of new commercial-grade systems in blind settings, as well as in areas with surface expressions of geothermal activity. Budget Period 2 of the project involved detailed analysis of some of the most promising areas identified in Phase 1. Twenty-four highly prospective areas, including both known undeveloped systems and previously undiscovered potential blind systems, were identified for further analysis. After reconnaissance of these areas, five of the most promising sites were selected for detailed studies. Multiple techniques were employed in the detailed studies, including geologic mapping, shallow temperature surveys, gravity surveys, Lidar, geochemical studies, seismic reflection analysis, and 3D modeling. The goal of the detailed studies was to identify specific areas with the highest likelihood for high permeability and thermal fluids, such that drill sites could be targeted. Three main sets of predictive maps were generated for each detailed study area: 1) play fairway maps, 2) play fairway error maps, and 3) direct evidence maps. Local- and intermediate-scale permeability models were revised to reflect results of the detailed geologic, geophysical, and geochemical analyses. Budget Period 3 of the project involved more detailed geophysical analyses and temperature-gradient (TG) drilling in southeastern Gabbs Valley and northern Granite Springs Valley, deemed the two most promising sites, with the goal of providing preliminary validation of the play fairway methodology. In southeastern Gabbs Valley, the collocation of a favorable structural setting (displacement transfer zone and fault intersections), Quaternary faults, intersecting and terminating gravity gradients, magnetic low, shallow (2 m) temperature anomaly, low resistivity anomaly, and promising geothermometry from nearby water wells provided evidence for a blind system. Drilling of six TG holes defines an apparent geothermal system at this locality with temperatures as high as 124°C at 152 m. This system is blind, with no surface hot springs, fumaroles, or paleo-geothermal deposits. For northern Granite Springs Valley, a favorable structural setting (termination of a major Quaternary normal fault), terminating gravity gradient, magnetic gradient, newly discovered sinter deposits, nearby warm water wells, previously drilled TG holes in the vicinity, and promising geothermometry suggest a hidden system. Drilling of six new TG holes yields temperatures of ~96°C at ~250 m, suggesting the presence of a geothermal system. Major lessons learned in the course of this project include: 1) initially identified sites commonly include multiple favorable structural settings at a finer scale; 2) promising sites in Cenozoic basins cannot be recognized without detailed geophysical surveys; and 3) play fairway analysis should be refined as the exploration program vectors into the most promising sites and finer-scale data are acquired. In addition to producing copious amounts of data, this project resulted in 16 published papers, 10 abstracts, more than 40 presentations across the U.S. and abroad (including several keynote addresses), 2 Masters theses, and 7 media reports.

15 GEOTHERMAL ENERGY↗

Advances in GDSA Framework Development and Process Model Integration

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and highlevel nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2020 advances of the Geologic Disposal Safety Assessment (GDSA) Framework and PFLOTRAN development groups of the SFWST Campaign. The common mission of these groups is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of disposal options and generic sites. The capability is a framework called GDSA Framework that employs high-performance computing (HPC) capable codes PFLOTRAN and Dakota.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Constraints on Thermal Evolution of Mars from Relaxation Models of Crustal and Topographic Dichotomy

The early thermal evolution of Mars is largely unconstrained. Models such as degree one convection [1,2,3], plate tectonics [4], and a transition to stagnant lid [5] have been proposed to explain formation of the dichotomy, the Tharsis rise, crustal production, and dynamo evolution. Here we model both the early deformation of the dichotomy and the long-term preservation as a means of examining the plausibility of a range of early thermal evolution models. Constraints include the preservation of crustal thickness and topographic differences between the northern and southern hemispheres and the geologic history of the dichotomy [6]). Our previous modeling indicates that the lower crust must have been weak enough to allow for relaxation early on, but the Martian interior had to cool fast enough to preserve the crustal difference and the associated topographic difference (~5 km) over approx. 3-3.5 Gyr [7].

Guest, A.↗