Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic 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 595 records · Page 33

A kill curve for Phanerozoic marine species

A kill curve for Phanerozoic species is developed from an analysis of the stratigraphic ranges of 17,621 genera, as compiled by Sepkoski. The kill curve shows that a typical species' risk of extinction varies greatly, with most time intervals being characterized by very low risk. The mean extinction rate of 0.25/m.y. is thus a mixture of long periods of negligible extinction and occasional pulses of much higher rate. Because the kill curve is merely a description of the fossil record, it does not speak directly to the causes of extinction. The kill curve may be useful, however, to li inverted question markmit choices of extinction mechanisms.

NASA Discipline Number 52-40↗

Modelling the Phanerozoic carbon cycle and climate: constraints from the 87Sr/86Sr isotopic ratio of seawater

A numerical model describing the coupled evolution of the biogeochemical cycles of carbon, sulfur, calcium, magnesium, phosphorus, and strontium has been developed to describe the long-term changes of atmospheric carbon dioxide and climate during the Phanerozoic. The emphasis is on the effects of coupling the cycles of carbon and strontium. Various interpretations of the observed Phanerozoic history of the seawater 87Sr/86Sr ratio are investigated with the model. More specifically, the abilities of continental weathering, volcanism, and surface lithology in generating that signal are tested and compared. It is suggested that the observed fluctuations are mostly due to a changing weatherability over time. It is shown that such a conclusion is very important for the modelling of the carbon cycle. Indeed, it implies that the conventional belief that the evolution of atmospheric carbon dioxide and climate on a long time scale is governed by the balance between the volcanic input of CO2 and the rate of silicate weathering is not true. Rather carbon exchanges between the mantle and the exogenic system are likely to have played a key role too. Further, the increase of the global weathering rates with increasing surface temperature and/or atmospheric CO2 pressure usually postulated in long-term carbon cycle and climate modelling is also inconsistent with the new model. Other factors appear to have modulated the weatherability of the continents through time, such as mountain building and the existence of glaciers and ice sheets. Based on these observations, a history of atmospheric carbon dioxide and climate during Phanerozoic time, consistent with the strontium isotopic data, is reconstructed with the model and is shown to be compatible with paleoclimatic indicators, such as the timing of glaciation and the estimates of Cretaceous paleotemperatures.

NASA Discipline Number 52-30↗

Large-body impact and extinction in the Phanerozoic

The kill curve for Phanerozoic marine species is used to investigate large-body impact as a cause of species extinction. Current estimates of Phanerozoic impact rates are combined with the kill curve to produce an impact-kill curve, which predicts extinction levels from crater diameter, on the working assumption that impacts are responsible for all "pulsed" extinctions. By definition, pulsed extinction includes the approximately 60% of Phanerozoic extinctions that occurred in short-lived events having extinction rates greater than 5%. The resulting impact-kill curve is credible, thus justifying more thorough testing of the impact-extinction hypothesis. Such testing is possible but requires an exhaustive analysis of radiometric dating of Phanerozoic impact events.

Non-NASA Center↗

Applications and Results of MODIS Lunar Observations

Both Terra and Aqua MODIS have been making periodic lunar observations since their launch in December 1999 and May 2002, respectively. The primary objective of the MODIS lunar observations is to monitor on-orbit response changes for the reflective solar bands (RSB) and to track long-term radiometric calibration stability. The Moon is an extremely stable surface reflectance reference and, together with the lunar radiometric model developed recently by the US Geological Survey (USGS), the applications of lunar observations for Earth remote sensing missions have been greatly enhanced and expanded. In particular, calibration inter-comparisons among sensors on different platforms and operated in different time frames have been enabled in the VIS/NIR/SWIR spectral regions. In this paper, we discuss various applications of MODIS lunar observations, such as monitoring the RSB radiometric stability and calibration consistency between the two sensors, evaluating the calibration differences among detectors within the same spectral band, and tracking the on-orbit changes of band-to-band registration (BBR). Other applications designed for the thermal emissive bands (TEB) and special sensor characterization purposes are also illustrated. Results from multi-year Terra and Aqua MODIS lunar observations demonstrate clearly that the two sensors have been calibrated consistently to within 1% for most RSB compared to the 2% radiometric calibration requirement. In addition the spatial characterization results derived from lunar observations are in good agreement with those derived from the instrument's on-board calibrator. The applications and results of MODIS lunar observations presented in this paper provide useful references for other sensors that also make use of lunar surface observations.

Xiong, X.↗

Combining Multi-Faceted Laboratory Studies of 74001-2 and Regional Remote Sensing to Address How Pyroclastic Eruptions Record and Affect the Lunar Volatile Budget

Basaltic magmatism is an efficient process for bringing volatiles from a planetary interior to its surface, with the possibility of generation of a transient lunar atmosphere as abundant volcanic materials de-gassed. However, pyroclastic deposits are locations where trapped gases may be studied [e.g., 2,3]. Volatile-rich pyroclastic deposits occur over a wide surface area of the Moon, indicating that the transport of volatiles and associated pyroclastic materials from the Moon’s mantle to the surface was a wide-spread phenomenon. Numerous studies analyzing remotely sensed data and using empirical modeling have demonstrated that various stages of pyroclastic eruptions significantly influence gas release patterns, morphology, and mineralogy of the deposit. Many observations based on mare basalts and pyroclastic deposits have identified potential histories of gas release [e.g., 2-10] and their influence on volatiles and their stable isotopes [e.g., 11-13]. The best representation of these pyroclastic deposits in the sample collection is core sample 74001-74002 that was collected during the Apollo 17 mission to the Taurus-Littrow Valley (TLV). The double drive tube penetrated a part of a regional-scale pyroclastic deposit and sampled approximately 68.1 cm of that deposit in the TLV. Remnants of this and other pyroclastic depos-its are represented throughout and beyond the TLV [e.g., 14-16]. The stratigraphy of this core has been investigated and defined by numerous studies. The CASA Moon SSERVI research team is conducting a multi-faceted analytical study of this deposit. Data generated from revisiting the stratigraphy of 74001-74002 will be used to place stable isotopes (H, B, Cl, S, Zn, Cu, Rb, Ga, Pb), Ar-Ar and U-Pb chronology, geochemistry, nanometer-scale observations of mineral surfaces, orbital observations, and experiments and modeling within a stratigraphic, eruptive, and geologic context. It is important to place these data into such a context. For example, recent S isotope measurements reported by Dottin et al. show differences within this core that may be related to either changes in source or eruptive process over the course of the eruption (vs. multiple eruptions). A fuller understanding of the stratigraphy is fundamental to resolving this interpretation. This comprehensive approach can only be achieved within the context of a program such as SSERVI. In addition, imaging produced in this project will be incorporated into a citizen scientist program to further identify many of the textural features of this double drive tube.

SSERVI↗

Investigation and Analysis of Dolomite Dissolution in Variable Strength Systems Relevant to the WIPP

The Waste Isolation Pilot Plant (WIPP) was opened in 1999 as a solution to the long-term disposal of the nation's nuclear defense-generated transuranic (TRU) waste. The salt beds were recommended as the ideal location because they are free from flowing fresh water, easily mined, impermeable, and geologically stable. However, for performance assessment modeling, release scenarios must be considered. The most likely release scenario was determined to be human intrusion leading to direct and/or long-term brine release (US DOE 1995, US DOE 1996, Perkins et al., 1999). The focus of this research is the Rustler formation because it is located above the WIPP and is the most transmissive layer (Perkins et al., 1999). Objective: To determine the effect of chelating ligands ( e.g. EDTA) on dissolution of dolomite in high ionic strength systems. Do brines impact dolomite dissolution and how is it affected by EDTA binding with cations? By varying the ionic strengths, different brine conditions can be simulated similar to those in and around the Culebra member and evaluate trends of dissolved Ca and Mg in the aqueous phase from the quantified data. Conclusions: <2% of dolomite dissolved over one week at pH 8.5. Greater dissolution occurred at higher ionic strength and in the presence of NaCl with approximately 0.4 vs. 0.9% of Mg removed from dolomite in 0.1 and 1.0 M Na. Further, 0.3 vs. 0.9% of Mg was removed from dolomite in 1.0 M IS CaCl{sub 2} versus NaCl. EDTA increased dissolution of dolomite at low ionic strength due to its strong complexation affinity for Ca/Mg with slightly greater removal of Ca than Mg. Significant differences were not observed at high ionic strength due to the presence of large concentrations of competing cations (Na, Ca, Mg). Slightly higher concentrations of Ca were removed from dolomite than Mg after one week, suggesting either incongruent dissolution or release through ion exchange of Ca from trace mineral impurities. Further, when EDTA was present, it may have increased Ca removal due to its greater complexation for Ca than Mg. Environmental Implications: High ionic strength brines and low ionic strength solutions with organic ligands increase dissolution of dolomite. Due to the long half-lives of radionuclides in the TRU waste, it may be important to consider the long term impact of dolomite dissolution for potential release scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluating Probability of Containment Effectiveness at a GCS Sites using integrated assessment modeling approach with Bayesian decision Networks

Improved scientific and engineering understanding of the behavior of geologic CO2 storage together with established regulatory framework and incentive structures raise the prospects for accelerated, large-scale deployment of this greenhouse gas emissions reduction approach. Incentive structures call for the establishment of appropriate verification and accounting approaches to support claims of the integrity of a geologic storage complex and to justify taking credit for long-term storage. In this study, we present a framework for assessing the probability of containment effectiveness over the lifetime of a geologic carbon storage site (e.g., after 70 years of injection and post-injection site performance) using forward stochastic model realizations based on site characterization data and using a monitoring-informed Bayesian network based on hypothetical detectability from surface seismic surveys over the site injection and post-injection phases. The National Risk Assessment Partnership’s open-source Integrated Assessment Model (NRAP-Open-IAM) was utilized to develop an ensemble of 10,000 a priori stochastic forecasts of CO2 containment. Those simulations were used to train the Bayesian network model to estimate the prior probabilities of the CO2 leakage mass into overlying, monitorable aquifers considering the uncertainties in the reservoir properties, permeability of potentially leaky wells and the overlying aquifers. The conditional probabilities in the Bayesian network were either learned from the NRAP-Open-IAM simulations or derived from the predefined detection thresholds for the monitoring method. Observations obtained from monitoring, over time during the site operation phases were then used to generate updated posterior probabilities of containment (and any loss from containment) in the Bayesian network by propagating the prior probabilities through the conditional probabilities. We demonstrate how to construct and use the Bayesian network for verifying the long-term storage complex effectiveness informed by monitoring based on the NRAP-Open-IAM simulations previously developed for the FutureGen 2.0 site. This approach may have relevance for stake holders to demonstrate secure geologic storage, provide a defensible, probabilistic approach to claim credit for geologic storage, and to estimate the likelihood that any fraction of the claimed credit may need to be refunded to the creditor based on available monitoring information.

Bayesian network, Risk assessment, Monitoring, car↗

The geology of Europa

The map units and lineations of Europa are detailed, and the geologic processes, and history, and thick and thin ice models of the satellite are discussed. It is concluded that Europa lacks evidence of a horizontally stratified crust; the geology appears characterized by disruption of the crust and intrusions into an icy shell. The surface consists of plains and mottled terrain, the former being older. Numerous straight and curved lineations, streaks, stripes, and bands cross EuropA's surface on a global and surface scale. Most lineations appear related to fractures in the crust. Five fresh craters in the 10 to 30 km diameter range are visible. The dark spots, stripes, and bands that appear to have replaced sections of the crust suggest that material was transported to the surface from the subjacent silicate lithosphere. The apparent low density of craters superposed on Europa's surface suggests that the surface is about 100 million years old.

Lucchitta, B. K.↗

Impact of GRM: New evidence from the Soviet Union

Gravity information released by the Soviet Union allows the quantitative assessment of how the geopotential research mission (GRM) mission might effect the ability to use global gravity data for continental tectonic interpretation. The information is of an isostatic response spectra for eight individual tectonic units in the USSR. The regions examined include the Caroathians, Caucasus, Urals, Pamirs, Tien-Shan, Altal, Chersky Ridge, and East Siberian Platform. The 1 deg x 1 deg gravity data are used to calculate the admittances are used in two different sorts of tectonic studies of mountain belts in the USSR: (1) interpretation of isostatic responses in terms of plate models of compensation for mountainous terrain. Using geologic information concerning time of the orogeny, lithospheric plates involved, and polarity of subduction in collision zones, they convert the best-fitting flexural rigidity to an elastic plate thickness for the lithospheric plate inferred to underlie the mountains; the isostatic admittance functions is an attempt to directly model gravity and topography data for a few select regions in the Soviet Union. By knowing the value of the expected correlation between topography and gravity from the admittances, the Artemjev's map in mountainous areas can be calibrated, and the maps are converted back to Bouguer gravity. This procedure is applied to the Caucasus and southern Urals.

Mcnutt, M.↗

Optimal carbon storage reservoir management through deep reinforcement learning

Model-based optimization plays a central role in energy system design and management. The complexity and high-dimensionality of many process-level models, especially those used for geosystem energy exploration and utilization, often lead to formidable computational costs when the dimension of decision space is also large. This work adopts elements of recently advanced deep learning techniques to solve a sequential decision-making problem in applied geosystem management. Specifically, a deep reinforcement learning framework was formed for optimal multiperiod planning, in which a deep Q-learning network (DQN) agent was trained to maximize rewards by learning from high-dimensional inputs and from exploitation of its past experiences. To expedite computation, deep multitask learning was used to approximate high-dimensional, multistate transition functions. Both DQN and deep multitask learning are pattern based. As a demonstration, the framework was applied to optimal carbon sequestration reservoir planning using two different types of management strategies: monitoring only and brine extraction. Both strategies are designed to mitigate potential risks due to pressure buildup. Results show that the DQN agent can identify the optimal policies to maximize the reward for given risk and cost constraints. Finally, experiments also show that knowledge the agent gained from interacting with one environment is largely preserved when deploying the same agent in other similar environments.

15 GEOTHERMAL ENERGY↗

Geoanalytical Evaluation of Saline Storage (GEESS) Geodatabase v2.0

The Geoanalytical Economic Evaluation of Saline Storage (GEESS) geodatabase was developed to support the United States Department of Energy (DOE) and National Energy Technology Laboratory (NETL) in their geologic carbon storage efforts by characterizing saline geologic formations present in the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM) [1]. Using publicly available literature and data, the GEESS geodatabase characterizes 57 geologic formations across the lower-48 U.S. states in what are called Fully Integrated Geodatabases (FIGs). The FIG is a vector polygon feature containing thousands or tens of thousands of individual polygons, which each contain discrete geologic parameter values. A list of the critical geologic parameters that are characterized in the GEESS geodatabase are described in the “Processing Steps and Workflow” part of the ReadMe file, as well as the Data Catalog accompanying the GEESS geodatabase. The FIG is the basis of the GEESS geodatabase and is the direct representation of the collected geologic data. In addition to the FIG, the GEESS system contains grid files. Due to the complexity of the FIGs, grids are used to sample the geologic data so they can be exercised within CO2_S_COM. The grid files contain the geologic data sampled from the FIG, as well as estimates of “Plume Uncertainty Diameter” and “First-year Break-even Price of CO2” derived from CO2_S_COM based on the GEESS grid data.

carbon↗

Venus climate stability and volcanic resurfacing rates

The climate of Venus is to a large degree controlled by the radiative properties of its massive atmosphere. In addition, outgassing due to volcanic activity, exospheric escape processes, and surface/atmosphere interactions may all be important in moderating the abundances of atmospheric CO2 and other volatiles. We have developed an evolutionary climate model for Venus using a systems approach that emphasizes feedbacks between elements in the climate system. Modules for atmospheric radiative transfer, surface/atmosphere interactions, tropospheric chemistry, and exospheric escape processes have so far been developed. Climate feedback loops result from interconnections between modules, in the form of the environmental parameters pressure, temperature, and atmospheric mixing ratios. The radiative transfer module has been implemented by using Rosseland mean opacities in a one dimensional grey radiative-convective model. The model has been solved for the static (time independent) case to determine climate equilibrium points. The dynamics of the model have also been explored by employing reaction/diffusion kinetics for possible surface atmosphere heterogeneous reactions over geologic timescales. It was found that under current conditions, the model predicts that the climate of Venus is at or near an unstable equilibrium point. The effects of constant rate volcanism and corresponding exsolution of volatiles on the stability of the climate model were also explored.

Bullock, M. A.↗

The Formation of Life-sustaining Planets in Extrasolar Systems

The spatial exploration is providing us a large quantity of information about the composition of the planets and satellites crusts. However, most of the experiences that are proposed in the guides of activities in Planetary Geology are based exclusively on the images utilization: photographs, maps, models or artistic reconstructions [1,2]. That things help us to recognize shapes and to deduce geological processes, but they says us little about the materials that they are implicated. In order to avoid this dicotomy between shapes and materials, we have designed an experience in the one which, employing of rocks and landscapes of our geological environment more next, the pupils be able to do an exercise of compared planetology analyzing shapes, processes and material of several planetary bodies of the Solar System.

Chambers, J. E.↗

Data for The Value of Reversible Carbon Storage in a Zero-Emissions World

Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.

Carbon↗

Dislocation model for aseismic fault slip in the transverse ranges of Southern California

Geodetic data at a plate boundary can reveal the pattern of subsurface displacements that accompany plate motion. These displacements are modelled as the sum of rigid block motion and the elastic effects of frictional interaction between blocks. The frictional interactions are represented by uniform dislocation on each of several rectangular fault patches. The block velocities and fault parameters are then estimated from geodetic data. Bayesian inversion procedure employs prior estimates based on geological and seismological data. The method is applied to the Transverse Ranges, using prior geological and seismological data and geodetic data from the USGS trilateration networks. Geodetic data imply a displacement rate of about 20 mm/yr across the San Andreas Fault, while the geologic estimates exceed 30 mm/yr. The prior model and the final estimates both imply about 10 mm/yr crustal shortening normal to the trend of the San Andreas Fault. Aseismic fault motion is a major contributor to plate motion. The geodetic data can help to identify faults that are suffering rapid stress accumulation; in the Transverse Ranges those faults are the San Andreas and the Santa Susana.

Cheng, A.↗

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗