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.

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Global Characterization of Polygonally Fractured Terrain on Venus and Implications for a Climate Change Origin

On Earth, polygonal fractures with spacings of 1-100 cm form in lava flows that cool slowly under isotropic stress conditions. On Venus, polygonal fractures are observed with typical fracture spacings of 1-2 km, covering regions 10s to 100s of km across. These features have been proposed to form via cooling of lava flows or above subsurface intrusions, or due to surface cooling caused by climate change. A watershed algorithm was modified to search the entire Magellan image database autonomously for terrains with polygonal fractures. Here we characterize the 204 polygon terrains found using this method (including over 170 newly identified regions) with respect to size, morphology, stratigraphy, and geologic setting. These results are used to examine different models of origin.

Smrekar, S. E.↗

NASA NTRSDate not supplied

Chesapeake Bay Water Resources: Examining Turbidity and Sediment Dynamics in the Chesapeake Bay

An increase in total suspended sediment (TSS) concentrations and turbidity have contributed to poor water quality in the Chesapeake Bay since the 1970s. Although turbidity and TSS have been moderately improving over the past few decades, poor water quality is detrimental to the Chesapeake Bay’s ecosystems and the surrounding watersheds. The Summer 2022 Chesapeake Bay Water Resources project observed sediment dynamics and turbidity in the York River watershed using remote sensing tools, models, and Earth observations including NASA and United States Geological Survey’s (USGS) Landsat satellite series, NASA DEVELOP’s Optical Reef Coastal Area Assessment Tool 2.0 (ORCAA), and Soil and Water Assessment Tool (SWAT). The collaborators on this project were the Chesapeake Bay National Estuarine Research Reserve (CBNERR), Group on Earth Observations (GEO) AquaWatch, the Committee on Earth Observation Satellites Coastal Observations, Applications, Services and Tools (CEOS COAST), and the Virginia Department of Environmental Quality (VA DEQ). The team concluded that the upstream sections of the York River watershed increased in TSS from 2009–2019. However, the seasonal TSS patterns often correlated with higher precipitation levels, although not necessarily major storm events. This may be due to factors like wind and waves contributing to the sedimentation trends observed. Although projects have been conducted to improve water quality, additional efforts like planting riparian buffers along areas with high TSS are needed to reduce the intensity of runoff. The team’s results will allow the end user, the VA DEQ, to inform their policymaking regarding future Bay conservation efforts.

Katherine Hahn↗

NASA NTRS2022

Chesapeake Bay Water Resources: Characterization of Sediment Dynamics for Enhanced Water Quality Monitoring in the Chesapeake Bay

Katherine Hahn↗

NASA NTRS1982

On isostasy at Atlantic-type continental margins

The concept of isostasy describes the manner in which topographic features on the earth's surface are compensated at depth. The present investigation is concerned with the isostatic mechanism at Atlantic-type continental margins. Particular attention is given to the question whether the flexure model of isostasy, which has successfully been used at other geological features in oceans, is applicable at margins. Cross-spectral techniques are used to analyze the relationship between free air gravity and topography at Atlantic-type continental margins. The relatively old eastern North America is found to be associated with the highest value of the effective elastic thickness in the range 10-20 km, while the relatively young Coral Sea/Lord Howe rise is associated with the lowest value of less than 5 km. The differences in estimates of effective elastic thickness between margins can be explained by a simple model in which the flexural strength of the basement increases with age.

Karner, G. D.↗

NASA NTRS1992

Origin of giant Martian polygons

Extensive areas of the Martian northern plains in Utopia and Acidalia planitiae are characterized by 'polygonal terrane'. Polygonal terrane consists of material cut by complex troughs defining a pattern resembling mudcracks, columnar joints, or frost-wedge polygons on earth. However, the Martian polygons are orders of magnitude larger than these potential earth analogues, leading to severe mechanical difficulties for genetic models based on simple analogy arguments. Plate-bending and finite element models indicate that shrinkage of desiccating sediment or cooling volcanics accompanied by differential compaction over buried topography can account for the stresses responsible for polygon troughs as well as the large size of the polygons. Although trough widths and depths relate primarily to shrinkage, the large scale of the polygonl pattern relates to the spacing between topographic elevations on the surface buried beneath polygonal terrane material. Geological relationships favor a sedimentary origin for polygonal terrane material, but our model is not dependent on the specific genesis. Our analysis also suggests that the polygons must have formed at a geologically rapid rate.

Mcgill, George E.↗

NASA NTRS1990

Mars atmospheric D/H - Consistent with polar volatile theory?

Current theories on the Martian water history are discussed on the basis of data for the D/H ratio in the Martian atmosphere (which was found to be enhanced by a factor of 6 + or - relative to the value on earth). In particular, the calculations of Yung et al. (1988), on the basis of their photochemical model, are found to be inconsistent with current theories for the evolution of the polar deposits on Mars. While the photochemical model predicts that about 95 percent of the Martian water has escaped over geologic time and the nonatmospheric reservoir is now only 02 m thick, the polar volatile theory suggest that polar deposits contain the equivalent of tens of meters of water, requiring a more rapid escape of H and D to space than predicted by the photochemical model to match the observed D/H ratio. The paper examines the behavior of the polar caps and the atmospheric water vapor over geologic time and the loss of H and D to space. It is concluded that, at present, it is premature to accept a specific history for water on Mars.

Jakosky, Bruce M.↗

DOE OSTI2024

Direct Geologic Constraints on the Timing of Late Holocene Ice Thickening in the Amundsen Sea Embayment, Antarctica

Abstract Constraining past West Antarctic Ice Sheet (WAIS) change helps validate numerical models simulating future ice sheet dynamics. Following rapid deglaciation during the mid‐Holocene, ice near Thwaites Glacier was ∼35 m thinner than present; however, the timing of ice regrowth to its present configuration remains unknown. To fill this knowledge gap, we present cosmogenic nuclide exposure ages of cobbles from the surface of a moraine situated between Thwaites and Pope glaciers. We infer that the moraine formed and stabilized in the Late Holocene (∼1.4 ka) when a small glacier thickened. We also present a novel reconstruction of WAIS volume constrained by sea‐level data, which demonstrates that moraine formation coincided with a large‐scale WAIS readvance. Our new geologic constraints will help inform models of the solid Earth response to surface mass loading, improving our understanding of ice sheet dynamics in a vulnerable part of WAIS.

Nichols, Keir A. [Imperial College London London U↗

DOE OSTI2022

NRAP-Open-IAM: Generic Aquifer Component Development and Testing

The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Quantum algorithms for geologic fracture networks

Abstract Solving large systems of equations is a challenge for modeling natural phenomena, such as simulating subsurface flow. To avoid systems that are intractable on current computers, it is often necessary to neglect information at small scales, an approach known as coarse-graining. For many practical applications, such as flow in porous, homogenous materials, coarse-graining offers a sufficiently-accurate approximation of the solution. Unfortunately, fractured systems cannot be accurately coarse-grained, as critical network topology exists at the smallest scales, including topology that can push the network across a percolation threshold. Therefore, new techniques are necessary to accurately model important fracture systems. Quantum algorithms for solving linear systems offer a theoretically-exponential improvement over their classical counterparts, and in this work we introduce two quantum algorithms for fractured flow. The first algorithm, designed for future quantum computers which operate without error, has enormous potential, but we demonstrate that current hardware is too noisy for adequate performance. The second algorithm, designed to be noise resilient, already performs well for problems of small to medium size (order 10–1000 nodes), which we demonstrate experimentally and explain theoretically. We expect further improvements by leveraging quantum error mitigation and preconditioning.

58 GEOSCIENCES↗

DOE OSTI2024

Comparison of Machine Learning-Based Predictive Models of the Nutrient Loads Delivered from the Mississippi/Atchafalaya River Basin to the Gulf of Mexico

Predicting nutrient loads is essential to understanding and managing one of the environmental issues faced by the northern Gulf of Mexico hypoxic zone, which poses a severe threat to the Gulf’s healthy ecosystem and economy. The development of hypoxia in the Gulf of Mexico is strongly associated with the eutrophication process initiated by excessive nutrient loads. Due to the complexities in the excessive nutrient loads to the Gulf of Mexico, it is challenging to understand and predict the underlying temporal variation of nutrient loads. The study was aimed at identifying an optimal predictive machine learning model to capture and predict nonlinear behavior of the nutrient loads delivered from the Mississippi/Atchafalaya River Basin (MARB) to the Gulf of Mexico. For this purpose, monthly nutrient loads (N and P) in tons were collected from US Geological Survey (USGS) monitoring station 07373420 from 1980 to 2020. Machine learning models—including autoregressive integrated moving average (ARIMA), gaussian process regression (GPR), single-layer multilayer perceptron (MLP), and a long short-term memory (LSTM) with the single hidden layer—were developed to predict the monthly nutrient loads, and model performances were evaluated by standard assessment metrics—Root Mean Square Error (RMSE) and Correlation Coefficient (R). The residuals of predictive models were examined by the Durbin–Watson statistic. The results showed that MLP and LSTM persistently achieved better accuracy in predicting monthly TN and TP loads compared to GPR and ARIMA. In addition, GPR models achieved slightly better test RMSE score than ARIMA models while their correlation coefficients are much lower than ARIMA models. Moreover, MLP performed slightly better than LSTM in predicting monthly TP loads while LSTM slightly outperformed for TN loads. Furthermore, it was found that the optimizer and number of inputs didn’t show effects on the LSTM performance while they exhibited impacts on MLP outcomes. This study explores the capability of machine learning models to accurately predict nonlinearly fluctuating nutrient loads delivered to the Gulf of Mexico. Further efforts focus on improving the accuracy of forecasting using hybrid models which combine several machine learning models with superior predictive performance for nutrient fluxes throughout the MARB.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

CO 2 plume evolution in a depleted natural gas reservoir: Modeling of conformance uncertainty reduction over time

Uncertainty in the long-term fate of CO 2 injected for geologic carbon sequestration (GCS) is a significant barrier to the adoption of GCS as a greenhouse-gas emission-mitigation for industry and regulatory agencies alike. Here, we present a modeling study that demonstrates that the uncertainty in forecasts of GCS site performance decreases over time as monitoring data are used to update operational models. We consider a case study of GCS in a depleted natural gas reservoir, with CO 2 injection occurring over 20 years, with a 50-year post-injection site care period. We constructed a detailed model to generate the actual model output, which is considered synthetic observation data. A series of simpler operational models based on limited data and assumptions about how an operator would model such a site are then run and compared against actual model output at specific monitoring points after one year, two years, etc. The operational model is updated and improved using the synthetic observation data from the actual model at the same time intervals. Model parameter values and model features needed to be updated over time to improve matches to the actual model. These kinds of model adjustments would be a normal part of reservoir engineering and site management at GCS sites. Uncertainty in two key measures related to site performance decreases with time: extent of the CO 2 plume up-dip migration, and radial extent of the pressure pulse. This conclusion should help allay the concerns of industry and regulators about uncertainty in long-term fate of CO 2 at GCS sites.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Enhancing 2D hydrodynamic flood models through machine learning and urban drainage integration

Two-dimensional hydrodynamic flood models are commonly employed for simulating flood extent and inundation depth. However, the influence of urban drainage network (UDN) is frequently overlooked in these models, potentially compromising their accuracy. Furthermore, the expensive computational costs and longer processing times make them challenging for large-scale hydrodynamic simulation. To address these challenges, this paper develops a machine learning (ML)-driven emulator for an open-source flood model, the Two-dimensional Runoff Inundation Toolkit for Operational Needs (TRITON). A TRITON-ML Emulator (TR-Emulator) that utilizes Convolutional Long Short-Term Memory is developed to capture the spatiotemporal features of flood events based on the outputs from TRITON. We further enhance the emulator by integrating UDN parameters (TR-UDN), such as the flow capacity of drainage pipes, pipe size, and pipe length, via an ML stacking technique to improve the water surface elevation (WSE) simulation. Hurricane Harvey 2017 in Houston, TX is used as the case study. We compare WSE results from TRITON, TR-Emulator, TR-UDN, and the United States Geological Survey (USGS) observations to evaluate the performance of these models. The results indicate that the TR-Emulator effectively replicates the WSE simulated by TRITON. Additionally, TR-UDN performs well in capturing WSE patterns and peak flows, aligning more closely with USGS observations, except in areas with milder slopes where conveyance discrepancies are observed. We further test the generalizability of our ML-based models using another smaller event. This paper shows that the TR-Emulator is effective for users and engineers to emulate a 2D hydrodynamic model, and the enhanced version of the TR-Emulator, TR-UDN, can be an efficient tool for predicting WSEs during urban flooding.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1991

Integration of diverse remote sensing data sets for geologic mapping and resource exploration

The use of high-quality multispectral images in the visible, near-infrared, shortwave infrared, thermal infrared, and microwave regions of the spectrum for producing thematic maps showing details of the surface geology is reported. The airborne data sets used in the study include the Airborne Visible/Infrared Imaging Spectrometer, the Thermal Infrared Multispectral Scanner, and the airborne SAR. Ancillary data include a digital elevation model, National High Altitude Photography, Landsat Multispectral Scanner data, Landsat Thematic Mapper data, laboratory and field spectral measurements, and traditional geologic mapping. The integrated, multispectral images are shown to provide new geologic information that can be used in mineral deposit models to provide exploration targets.

Kruse, Fred A.↗

DOE OSTI2024

Efficient Dimension Reduction of Complex Three-dimensional CO2 Saturation using Deep Learning Models

In the domain of deep learning (DL), dimension reduction is crucial for enhancing training efficiency and mitigating overfitting, particularly when managing complex data such as three-dimensional (3D) saturation data. The 3D saturation data in the context of geological carbon storage (GCS) presents unique challenges due to its inherent sparsity and the abrupt transitions at plume boundaries, known as shock fronts. To address the challenges, we proposed a novel DL framework that integrates dimension reduction with advanced 3D reconstruction techniques. Our model leveraged latent variables derived from 2D average saturation data, offering a robust and efficient solution tailored to the intricate dynamics of 3D saturation fields. The proposed framework can extract the critical features of the high-dimensional data while reducing the variable numbers, which is more tractable for DL models and enhances the model robustness and accuracy. Therefore, it provides a novel approach for modeling and analyses in complex geological scenarios, which finds great potential applications in environmental monitoring and energy storage.

Wang, Hongsheng↗

DOE OSTI2022

Cost Impacts of Risk-Based Methods for Defining AoR and PISC Duration of a CO 2 Storage Project Using NRAP Tools and FE/NETL CO 2 Saline Storage Cost Model

The purpose of this analysis is to assess the fiscal impact of using risk-based methods to define the area of review (AoR) and post-injection site care (PISC) duration for a geologic carbon dioxide (CO 2 ) storage project. The Office of Fossil Energy (FE)/National Energy Technology Laboratory (NETL) CO 2 Saline Storage Cost Model (“cost model”) was utilized to evaluate CO 2 storage economics in this study. The cost model is a widely used tool for evaluating CO 2 storage costs in geologic settings. It estimates costs for storage operations conducted under compliance to the United States Environmental Protection Agency (EPA) Underground Injection Control (UIC) Class VI regulations. This analysis evaluates and compares the costs of storage (on a first-year break-even dollar per tonne [2018$] basis) at the proposed FutureGen 2.0 site as a case study storage location by implementing different approaches to determine AoR and PISC given the prevailing geologic conditions and injection volume and duration considerations. The approaches include the use of 1) risk-based methods to define AoR and PISC duration, 2) an EPA-approved AoR and PISC documented in the FutureGen 2.0 UIC Class VI permit applications, and 3) cost model default settings that utilize uncertainty multipliers to estimate CO 2 plume and pressure front extent as part of AoR determination, as well as a 50-year PISC default. From all three sources, the AoR and PISC specified in the Class VI permit applications for FutureGen 2.0 were the largest and longest relative to the other two sources.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Representing Global Soil Erosion and Sediment Flux in Earth System Models

Abstract Soil erosion produces enormous amounts of sediment, carbon, and nutrient fluxes from land to rivers, thus playing crucial roles in global biogeochemical cycles and food security. To predict soil erosion in the context of climate and land use changes, we explicitly parameterize cropland management actions (i.e., irrigation, conserved agriculture, and crop residue management) and geological factors (i.e., lithology and glacier) in the Energy Exascale Earth System Model (E3SM) soil erosion module. The erosion model is calibrated using a global‐scale regionalized parameter calibration method. The spatial variabilities of the modeled and the Revised Universal Soil Loss Equation (RUSLE) based benchmark soil erosion are consistent across vegetation, climate and soil properties. Compared with independent data, our model shows a bias reduction in 59% of the observations relative to the RUSLE‐based soil erosion, with 53% of the bias reduction exceeding 50%. This improvement is mainly due to a better representation of the topographic effect on soil erosion. Our results indicate that conserved agriculture practices have effectively reduced soil erosion in cropland by over 25% in the United States and Argentina. In contrast, irrigation has increased soil erosion in many Asian countries. For upland sediment flux, our model is consistent with the WBMsed benchmark data in inter‐basin variability but could be more skillful in simulating intra‐basin variability because it couples soil erosion and sediment flux explicitly. The developed model provides useful skills for more realistic predictions of soil erosion and river sediment dynamics under environmental changes.

58 GEOSCIENCES↗

DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗

DOE OSTI2021

Machine Learning Analysis of Hydrologic Exchange Flows and Transit Time Distributions in a Large Regulated River

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. High-resolution numerical models were often used to resolve the spatial and temporal variations of exchange flows, which are computationally expensive. In this study, we adopt Random Forest (RF) and Extreme Gradient Boosting (XGB) approaches for deriving reduced order models of hydrologic exchange flows and associated transit time distributions, with integrated field observations (e.g., bathymetry) and hydrodynamic simulation data (e.g., river velocity, depth). The setup allows an improved understanding of the influences of various physical, spatial, and temporal factors on the hydrologic exchange flows and transit times. The predictors also contain those derived using hybrid clustering, leveraging our previous work on river corridor system hydromorphic classification. The machine learning-based predictive models are developed and validated along the Columbia River Corridor, and the results show that the top parameters are the thickness of the top geological formation layer, the flow regime, river velocity, and river depth; the RF and XGB models can achieve 70% to 80% accuracy and therefore are effective alternatives to the computational demanding numerical models of exchange flows and transit time distributions. Each machine learning model with its favorable configuration and setup have been evaluated. The transferability of the models to other river reaches and larger scales, which mostly depends on data availability, is also discussed.

97 MATHEMATICS AND COMPUTING↗