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 649 records · Page 36

Earth survey applications division: Research leading to the effective use of space technology in applications relating to the Earth's surface and interior

Accomplishments and future plans are described for the following areas: (1) geology - geobotanical indicators and geopotential data; (2) modeling magnetic fields; (3) modeling the structure, composition, and evolution of the Earth's crust; (4) global and regional motions of the Earth's crust and earthquake occurrence; (5) modeling geopotential from satellite tracking data; (6) modeling the Earth's gravity field; (7) global Earth dynamics; (8) sea surface topography, ocean dynamics; and geophysical interpretation; (9) land cover and land use; (10) physical and remote sensing attributes important in detecting, measuring, and monitoring agricultural crops; (11) prelaunch studies using LANDSAT D; (12) the multispectral linear array; (13) the aircraft linear array pushbroom radiometer; and (14) the spaceborne laser ranging system.

Carpenter, L.↗

Seasonal thermal energy storage in aquifers: Mathematical modeling studies in 1979

A numerical model of water and heat flow in geologic media was developed, verified, and tested. The hydraulic parameters (transmittivity and storativity) and the location of a linear hydrologic barrier were simulated and compared with results from field experiments involving two injection-storage-recovery cycles. For both cycles, the initial simulated and observed temperatures agree (55c).

Tsang, C. F.↗

GRM as a follow-on to MAGSAT

MAGSAT was the first near-Earth spacecraft dedicated to measurement of the vector geomagnetic field. Its objectives were to use the vector measurements to measure and model both the main field and the crustal field of the Earth. The main field measurements were to be used for the 1980 world charts and to study the Earth's core. The crustal fields were to be utilized to aid in modeling large-scale variations in the geologic and geophysical characteristics of the crust.

Langel, R. A.↗

Rainy Lake wrench zone: An example of an Archaean subprovince boundary in northwestern Ontario

The Superior Province of the Canadian Shield comprises an alternation of subprovinces with contrasting lithological, structural and metamorphic styles. Rocks of the Rainly Lake area form a fault bounded wedge between two of these subprovinces, the Wabigoon granite-greenstone terrain to the north and the Quetico metasedimentary terrain to the south. The Quetico and Seine River-Rainy Lake Faults bound this wedge within which interpretation of the stratigraphy has been historically contentious. In the eastern part of the wedge, volcanic rocks and coeval tonalitic sills are unconformably overlain by fluviatile conglomerate and arenite of the Seine Group; in the western part of the wedge, metamorphosed wacke and mudstone of the Coutchiching Group are cut by granodioritic plutons. The Coutchiching Group has previously been correlated with the Seine Group and with the turbiditic Quetico metasediments of the Quetico Subprovince and these correlations are the cornerstone of earlier tectonic models which relate the subprovinces. The structural geology of the Rainy Lake area is characterized by attributes which compare favourably with the known characteristics of dextral wrench or 'transpressive zones based both on experimental data and natural examples. Much of this deformation involved the Seine Group, the youngest stratigraphic unit in the area, and predates the emplacement of late-to-post-tectonic granodioritic plutons for which radiometric data indicate a Late Archean age.

Poulsen, K. H.↗

Maximum likelihood classification of synthetic aperture radar imagery

Classification of synthetic aperture radar (SAR) images has important applications in geology, agriculture, and the military. A statistical model for SAR images is reviewed and a maximum likelihood classification algorithm developed for the classification of agricultural fields based on the model. It is first assumed that the target feature information is known a priori. The performance of the algorithm is then evaluated in terms of the probability of incorrect classification. A technique is also presented to extract the needed feature information from a SAR image; then both the feature extraction and the maximum likelihood classification algorithms are tested on a SEASAT-A SAR image.

Frost, V. S.↗

Volcanic geology of Tyrrhena Patera, Mars

Consideration is given to the geology of Tyrrhena Patera, a large low-relief volcano in the southern cratered highlands of Mars. The general geology of Tyrrhena Patera is outlined and models for the formation of the volcano are described. Models derived from studies of terrestrial pyroclastic flows are applied to deposits at Tyrrhena Patera, showing that the characteristics of the deposits are consistent with an origin by the emplacement of gravity-driven ash flows generated by hydromagmatic or magmatic explosive eruptions.

Greeley, Ronald↗

Terrestrial glacial eskers: Analogs for Martian sinuous ridges

A glacial model was introduced last year for the Argyre region, a concept which is now extended, and which was recently integrated with a Global Hydrologic Model incorporating many other aspects of Martian geology. Despite wide agreement that the Martian ridges strongly resemble glacial eskers, this hypothesis has been presented with great equivocation due to a perceived lack of other glacial landforms. Quite to the contrary, it is shown that the Martian ridges actually do occur in logical ordered sequences with many other types of characteristically glacial appearing landforms. Herein, the esker hypothesis is further supported in isolation from considerations of regional landform assemblages. It is concluded that Martian sinuous ridges are similar in every respect to terrestrial eskers: scale, morphology, planimetric pattern, and associations with other probable glaciogenic landforms. It is found that the esker hypothesis is well supported. Eskers are glaciofluvial structures, and owe their existence to large scale melting of stagnant temporate glaciers. Thus, eskers are indicators of an ameliorating climatic regime after a protracted episode of cold, humid conditions.

Kargel, Jeffrey S.↗

BOREAS Regional DEM in Raster Format and AEAC Projection

This data set is based on the GTOPO30 Digital Elevation Model (DEM) produced by the United States Geological Survey EROS Data Center (USGS EDC). The BOReal Ecosystem-Atmosphere Study (BOREAS) region (1,000 km x 1000 km) was extracted from the GTOPO30 data and reprojected by BOREAS staff into the Albers Equal-Area Conic (AEAC) projection. The pixel size of these data is 1 km. The data are stored in binary, image format files.

Knapp, David↗

Low-Altitude AVIRIS Data for Mapping Land Cover in Yellowstone National Park: Use of Isodata Clustering Techniques

Northeast Yellowstone National Park (YNP) has a diversity of forest, range, and wetland cover types. Several remote sensing studies have recently been done in this area, including the NASA Earth Observations Commercial Applications Program (EOCAP) hyperspectral project conducted by Yellowstone Ecosystems Studies (YES) on the use of hyperspectral imaging for assessing riparian and in-stream habitats. In 1999, YES and NASA's Commercial Remote Sensing Program Office began collaborative study of this area, assessing the potential of synergistic use of hyperspectral, synthetic aperture radar (SAR), and multiband thermal data for mapping forest, range, and wetland land cover. Since the beginning, a quality 'reference' land cover map has been desired as a tool for developing and validating other land cover maps produced during the project. This paper recounts an effort to produce such a reference land cover map using low-altitude Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and unsupervised classification techniques. The main objective of this study is to assess ISODATA classification for mapping land cover in Northeast YNP using select bands of low-altitude AVIRIS data. A secondary, more long-term objective is to assess the potential for improving ISODATA-based classification of land cover through use of principal components analysis and minimum noise fraction (MNF) techniques. This paper will primarily report on work regarding the primary research objective. This study focuses on an AVIRIS cube acquired on July 23, 1999, by the confluence of Soda Butte Creek with the Lamar River. Range and wetland habitats dominate the image with forested habitats being a comparatively minor component of the scene. The scene generally tracks from southwest to northeast. Most of the scene is valley bottom with some lower side slopes occurring on the western portion. Elevations within the AVIRIS scene range from approximately 1998 to 2165 m above sea level, based on US Geological Survey (USGS) 30-m digital elevation model (DEM) data. Despain and the National Park Service (NPS) provide additional description of the study area.

Spruce, Joseph P.↗

On-Orbit Cross-Calibration of AM Satellite Remote Sensing Instruments using the Moon

On April 14,2003, three Earth remote sensing spacecraft were maneuvered enabling six satellite instruments operating in the visible through shortwave infrared wavelength region to view the Moon for purposes of on-orbit cross-calibration. These instruments included the Moderate Resolution Imaging Spectroradiometer (MODIS), the Multi-angle Imaging SpectroRadiometer (MISR), the Advanced Spaceborne Thermal Emission and Reflection (ASTER) radiometer on the Earth Observing System (EOS) Terra spacecraft, the Advanced Land Imager (ALI) and Hyperion instrument on Earth Observing-1 (EO-1) spacecraft, and the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) on the SeaStar spacecraft. Observations of the Moon were compared using a spectral photometric mode for lunar irradiance developed by the Robotic Lunar Observatory (ROLO) project located at the United States Geological Survey in Flagstaff, Arizona. The ROLO model effectively accounts for variations in lunar irradiance corresponding to lunar phase and libration angles, allowing intercomparison of observations made by instruments on different spacecraft under different time and location conditions. The spacecraft maneuvers necessary to view the Moon are briefly described and results of using the lunar irradiance model in comparing the radiometric calibration scales of the six satellite instruments are presented here.

Butler, James J.↗

Using Virtual Reality for Science Missions At The Lunar South Pole

The Lunar VR toolkit, built on Mixed Reality Exploration Toolkit (MRET), combines LOLA 5m topographic data, CAD models (e.g., lander, rover), procedural textures, and geologic features (e.g., rocks, craters) for sub-5m simulation of the lunar surface. Future capabilities include the incorporation of real-time mission telemetry to provide a full mission lifecycle tool, test instrument design and operation concepts, pre-operations planning and walk-throughs of traverses, and situational awareness support during actual mission operations. We will discuss the critical features we think are needed for lunar south pole mission planning and the development of our Lunar VR toolkit.

Thomas G Grubb↗

Abundant phosphorus expected for possible life in Enceladus’s ocean

Saturn’s moon Enceladus has a potentially habitable subsurface water ocean that contains canonical building blocks of life (organic and inorganic carbon, ammonia, possibly hydro-gen sulfide) and chemical energy (disequilibria for methanogenesis). However, its habit-ability could be strongly affected by the unknown availability of phosphorus (P). Here, we perform thermodynamic and kinetic modeling that simulates P geochemistry based on recent insights into the geochemistry of the ocean–seafloor system on Enceladus. We find that aqueous P should predominantly exist as orthophosphate (e.g., HPO422), and total dissolved inorganic P could reach 1027to 1022mol/kg H2O, generally increasing with lower pH and higher dissolved CO2, but also depending upon dissolved ammonia and silica. Levels are much higher than<10210mol/kg H2Ofrompreviousestimatesand close to or higher than∼1026mol/kg H2O in modern Earth seawater. The high P concentration is primarily ascribed to a high (bi)carbonate concentration, which decreases the concentrations of multivalent cations via carbonate mineral formation, allowing phosphate to accumulate. Kinetic modeling of phosphate mineral dissolution suggests that geologically rapid release of P from seafloor weathering of a chondritic rocky core could supply millimoles of total dissolved P per kilogram of H2Owithin105y, much less than the likely age of Enceladus’s ocean (108to 109y). These results provide further evidence of habitable ocean conditions and show that any oceanic life would not be inhibited by low P availability.

Enceladus ocean↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

Subsurface fluvial sediments beneath InSight on Mars from geophysical constraints

Subsurface structure investigation on Mars is crucial for understanding its geological evolution and past hydrological conditions. Elysium Planitia (EP), located near the hypothesized ancient ocean shorelines, could contain clues for past water activity and paleoclimate. Here we present better-constrained subsurface models beneath InSight extending to ~800 m depth, obtained from joint inversion of seismic and seismoacoustic coupling data, and use the well-resolved subsurface structure to explore the lithological profile through rock physics models. The derived subsurface lithology agrees well with local geological context and exhibits a shallow 60-m-thick low-rigidity layer consistent with hydrated sedimentary materials. Despite possible contributions of aeolian and volcanic deposits, we favor the interpretation that the low-rigidity layer originated from fluvial activity in EP during the Hesperian or Hesperian-to-Amazonian epoch, as supported by adjacent paleo-shoreline morphology observations. These results hint at a period of warmer paleoclimate at low latitudes, possibly during high-obliquity phases of Mars’ rotational axis.

58 GEOSCIENCES↗

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗

Ice-Covered Lakes in Gale Crater Mars: The Cold and Wet Hypothesis

Recent geological discoveries from the Mars Science Laboratory provide evidence that Gale crater may have intermittently hosted a fluvio-lacustine environment during the Hesperian, with individual lakes lasting for a period of tens to hundreds of thousands of years. (Grotzinger et al., Science, 350 (6257), 2015). Estimates of the CO2 content of the atmosphere at the time the Gale sediments formed are far less than needed by any climate model to warm early Mars (Bristow et al., Geology, submitted), given the low solar energy input available at Mars 3.5 Gya. We have therefore explored the possibility that the lakes in Gale during the Hesperian were perennially covered with ice using the Antarctic Lakes as an analog. Using our best estimate for the annual mean surface temperature at Gale at this time (approx. 230K) we computed the thickness of an ice-covered lake. These thickness range from 10-30 meters depending on the ablation rate and ice transparency and would likely inhibit sediments from entering the lake. Thus, a first conclusion is that the ice must not be too cold. Raising the mean temperature to 245K is challenging, but not quite as hard as reaching 273K. We found that a mean annual temperature of 245K ice thicknesses range from 3-10 meters. These values are comparable to the range of those for the Antarctic lakes (3-6 m), and are not implausible. And they are not so thick that sediments cannot penetrate the ice. For the ice-covered lake hypothesis to work, however, a melt water source is needed. This could come from subaqueous melting of a glacial dam in contact with the lakes (as is the case for Lake Untersee) or from seasonal melt water from nearby glaciers (as is the case for the Dry Valley lakes). More work is needed to better assess these possibilities. However, the main advantage of the ice-covered lake model (and the main reason we pursued it) is that it relaxes the requirement for a long-lived active hydrological cycle involving rainfall and runoff, which no climate model is able to produce given known constraints on the early Mars environment.

Antarctic lakes↗