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 901 records · Page 50

Physical characteristics and evolutionary trends of continental rifts

Rifts may be defined as zones beneath which the entire lithosphere has ruptured in extension. They are widespread and occur in a variety of tectonic settings, and range up to 2,600 m.y. in age. The object of this review is to highlight characteristic features of modern and ancient rifts, to emphasize differences and similarities in order to help characterize evolutionary trends, to identify physical conditions favorable for initiation as well as termination of rifting, and to provide constraints for future modeling studies of rifting. Rifts are characterized on the basis of their structural, geomorphic, magmatic and geophysical features and the diverse character of these features and their evolutionary trends through time are discussed. Mechanisms of rifting are critically examined in terms of the physical characteristics and evolutionary trends of rifts, and it is concluded that while simple models can give valuable insight into specific processes of rifting, individual rifts can rarely, if ever, be characterized by well defined trends predicted by these models. More data are required to clearly define evolutionary trends, and the models require development to incorporate the effects of lithospheric heterogeneities and complex geologic histories.

Ramberg, I. B.↗

A simple tectonic model for crustal accretion in the Slave Province: A 2.7-2.5 Ga granite greenstone terrane

A prograding (direction unspecified) trench-arc system is favored as a simple yet comprehensive model for crustal generation in a 250,000 sq km granite-greenstone terrain. The model accounts for the evolutionary sequence of volcanism, sedimentation, deformation, metamorphism and plutonism, observed througout the Slave province. Both unconformable (trench inner slope) and subconformable (trench outer slope) relations between the volcanics and overlying turbidities; and the existence of relatively minor amounts of pre-greenstone basement (microcontinents) and syn-greenstone plutons (accreted arc roots) are explained. Predictions include: a varaiable gap between greenstone volcanism and trench turbidite sedimentation (accompanied by minor volcanism) and systematic regional variations in age span of volcanism and plutonism. Implications of the model will be illustrated with reference to a 1:1 million scale geological map of the Slave Province (and its bounding 1.0 Ga orogens).

Hoffman, P. F.↗

Lithospheric loading and tectonics of the lunar irregular maria

Many of the tectonic features associated with the circular maria are the result of stresses generated by loading of the lunar lithosphere by mare basalt units. The hypothesis is tested that the irregular lunar maria share with the circular maria the processes leading to the formation of associated graben and mare ridge systems. A formulation of the lithospheric flexure problem that accounts for the variable distribution of basalt loads for the irregular maria is applied. On the basis of the tectonic structures and geologic history of each of the major irregular maria, as well as models for the distribution of mare basalt for each region, the predicted stress fields are compared with the distribution of tectonic features for a range of assumed values of the thickness of the elastic lithosphere. The best fitting values for lithospheric thickness beneath the regular maria are then compared with those previously inferred for the circular maria.

Hall, J. Lynn↗

Global Stratigraphy of Venus: Analysis of a Random Sample of Thirty-Six Test Areas

The age relations between 36 impact craters with dark paraboloids and other geologic units and structures at these localities have been studied through photogeologic analysis of Magellan SAR images of the surface of Venus. Geologic settings in all 36 sites, about 1000 x 1000 km each, could be characterized using only 10 different terrain units and six types of structures. These units and structures form a major stratigraphic and geologic sequence (from oldest to youngest): (1) tessera terrain; (2) densely fractured terrains associated with coronae and in the form of remnants among plains; (3) fractured and ridged plains and ridge belts; (4) plains with wrinkle ridges; (5) ridges associated with coronae annulae and ridges of arachnoid annulae which are contemporary with wrinkle ridges of the ridged plains; (6) smooth and lobate plains; (7) fractures of coronae annulae, and fractures not related to coronae annulae, which disrupt ridged and smooth plains; (8) rift-associated fractures; and (9) craters with associated dark paraboloids, which represent the youngest 1O% of the Venus impact crater population (Campbell et al.), and are on top of all volcanic and tectonic units except the youngest episodes of rift-associated fracturing and volcanism; surficial streaks and patches are approximately contemporary with dark-paraboloid craters. Mapping of such units and structures in 36 randomly distributed large regions (each approximately 10(exp 6) sq km) shows evidence for a distinctive regional and global stratigraphic and geologic sequence. On the basis of this sequence we have developed a model that illustrates several major themes in the history of Venus. Most of the history of Venus (that of its first 80% or so) is not preserved in the surface geomorphological record. The major deformation associated with tessera formation in the period sometime between 0.5-1.0 b.y. ago (Ivanov and Basilevsky) is the earliest event detected. In the terminal stages of tessera fon-nation, extensive parallel linear graben swarms representing a change in the style of deformation from shortening to extension were formed on the tessera and on some volcanic plains that were emplaced just after, and perhaps also during the latter stages of the major compressional phase of tessera emplacement. Our stratigraphic analyses suggest that following tessera formation, extensive volcanic flooding resurfaced at least 85% of the planet in the form of the presently-ridged and fractured plains. Several lines of evidence favor a high flux in the post-tessera period but we have no independent evidence for the absolute duration of ridged plains emplacement. During this time, the net state of stress in the lithosphere apparently changed from extensional to compressional, first in the form of extensive ridge belt development, followed by the formation of extensive wrinkle ridges on the flow units. Subsequently, there occurred local emplacement of smooth and lobate plains units which are presently essentially undefortned. The major events in the latest 10% of the presently preserved history of Venus (less than 50 m.y. ago) are continued rifting and some associated volcanism, and the redistribution of eolian material largely derived from impact crater deposits.

Basilevsky, Alexander T.↗

Seafloor hydrothermal activity and spreading rates: the Eocene carbon dioxide greenhouse revisted

A suggestion has been made that enhanced rates of hydrothermal activity during the Eocene could have caused a global warming by adding calcium to the ocean and pumping CO2 into the atmosphere (Owen and Rea, 1984). This phenomenon was purported to be consistent with the predictions of the CO2 geochemical cycle model of Berner, Lasaga and Garrels (1983) (henceforth BLAG). In fact, however, the BLAG model predicts only a weak connection between hydrothermal activity and atmospheric CO2 levels. By contrast, it predicts a strong correlation between seafloor spreading rates and pCO2, since the release rate of CO2 from carbonate metamorphism is assumed to be proportional to the mean spreading rate. The Ecocene warming can be conveniently explained if the BLAG model is extended by assuming that the rate of carbonate metamorphism is also proportional to the total length of the midocean ridges from which the spreading originates.

Seawater/chemistry↗

Spectroscopy of Mars from 2.04 to 2.44 micrometers during the 1993 opposition: absolute calibration and atmospheric vs mineralogic origin of narrow absorption features

We present moderate-resolution (lambda/delta lambda = 300 to 370) reflectance spectra of Mars from 2.04 to 2.44 micrometers that were obtained at UKIRT during the 1993 opposition. Seven narrow absorption features were detected and found to have a Mars origin. By comparison with solar and Mars atmospheric spectra, five of these features were attributed all or in part to Mars atmospheric CO2 or CO(2.052 +/- 0.003, 2.114 +/- 0.002, 2.150 +/- 0.003, 2.331 +/- 0.001, and 2.357 +/- 0.002 micrometers). Two of the bands (2.331 +/- 0.001 and 2.357 +/- 0.002 micrometers) appear to have widths and depths that are consistent with additional, nonatmospheric absorptions, although a solar contribution cannot be entirely ruled out. Two other weak bands centered at 2.278 +/- 0.002 and 2.296 +/- 0.002 micrometers may be at least partially mineralogic in origin. The data provide no conclusive identification of the mineralogy responsible for these absorption features. However, examination of terrestrial spectral libraries and previous moderate spectral resolution mineral studies indicates that the most likely origin of these features is either (bi)carbonate or (bi)sulfate anions in framework silicates or (Fe, Mg)-OH bonds in sheet silicates. If the bands are caused by phyllosilicate minerals, then an explanation must be found for the extremely narrow widths of the cation-OH features in the Mars spectra as compared to terrestrial minerals.

Mars↗

Constraints on Early Mars Evolution and Dichotomy Origin from Relaxation Modeling of Dichotomy Boundary in the Ismenius Region

The Martian dichotomy is a global feature separating the northern and southern hemispheres. The 3.5 - 4 Gyr old feature is manifested by a topographic difference of 2-6 km and crustal thickness difference of approx. 15 - 30 km between the two hemispheres. In the Ismenius region, sections of the boundary are characterized by a single scarp with a slope of approx. 20 deg. - 23 deg. and are believed to be among the most well preserved parts of the dichotomy boundary. The origin of the dichotomy is unknown. Endogenic hypotheses do not predict the steep slopes (scarps) of the dichotomy boundary. Exogenic models for forming the northern lowlands by impact cratering, associate the scarps along the dichotomy boundary with craters' rims, but are not globally consistent with the topography and gravity. In order to better understand the origin of the Martian dichotomy, it is necessary to know if the steep scarps along the boundary represent the original shape of the dichotomy. Smrekar et al. presented evidence showing that the boundary scarp in Ismenius is a fault along which the highland crust was down faulted. We test whether the relaxation process could produce faulting along the dichotomy boundary and examine the crustal and mantle conditions that would allow for faulting to occur within 1 Gyr and preserve the long wavelength topography over another 3 Gyr. We approach the problem by a combination of numerical and semi-analytical modeling. We test different viscosity profiles and crustal thicknesses by comparing our modeled magnitude, location and timing of plastic strain and displacements to detailed geologic observations in the Ismenius region.

Guest, A.↗

Aqueous Alteration of Basalts: Earth, Moon, and Mars

The geologic processes responsible for aqueous alteration of basaltic materials on Mars are modeled beginning with our knowledge of analog processes on Earth, i.e., characterization of elemental and mineralogical compositions of terrestrial environments where the alteration and weathering pathways related to aqueous activity are better understood. A key ingredient to successful modeling of aqueous processes on Mars is identification of phases that have formed by those processes. The purpose of this paper is to describe what is known about the elemental and mineralogical composition of aqueous alteration products of basaltic materials on Mars and their implications for specific aqueous environments based upon our knowledge of terrestrial systems. Although aqueous alteration has not occurred on the Moon, it is crucial to understand the behaviors of basaltic materials exposed to aqueous environments in support of human exploration to the Moon over the next two decades. Several methods or indices have been used to evaluate the extent of basalt alteration/weathering based upon measurements made at Mars by the Mars Exploration Rover (MER) Moessbauer and Alpha Particle X-Ray Spectrometers. The Mineralogical Alteration Index (MAI) is based upon the percentage of total Fe (Fe(sub T)) present as Fe(3+) in alteration products (Morris et al., 2006). A second method is the evaluation of compositional trends to determine the extent to which elements have been removed from the host rock and the likely formation of secondary phases (Nesbitt and Young, 1992; Ming et al., 2007). Most of the basalts that have been altered by aqueous processes at the two MER landing sites in Gusev crater and on Meridiani Planum have not undergone extensive leaching in an open hydrolytic system with the exception of an outcrop in the Columbia Hills. The extent of aqueous alteration however ranges from relatively unaltered to pervasively altered materials. Several experimental studies have focused upon the aqueous alteration of lunar materials and simulants (e.g., Keller and Huang, 1971; Eick et al., 1996). Lunar basalts are void of water and highly reduced, hence, these materials are initially very reactive when exposed to water under oxidizing conditions.

Ming, Douglas W.↗

Collapse of the Late Proterozoic ecosystem

Evidence accumulated over the past two decades is now sufficient to permit an initial quantitative assessment of the patterns of biotic diversity and extinction that occurred during Proterozoic time. Because of limitations in both the quality and quantity of data currently available, however, generalizations thus derived must be regarded as tentative. Nevertheless, read literally, available palaeontological data appear to indicate that the global ecosystem experienced a gradual but massive collapse between 1 000 Ma and the beginning of the Phanerozoic, a supposition consistent with other lines of geological and geochemical evidence. A possible forcing agent for such a collapse appears to have been a decrease in ambient levels of carbon dioxide and a resultant decrease in average global temperature, photosynthetic efficiency, and primary productivity.

NASA Discipline Exobiology↗

Hydrothermal systems on Mars: an assessment of present evidence

Hydrothermal processes have been suggested to explain a number of observations for Mars, including D/H ratios of water extracted from Martian meteorites, as a means for removing CO2 from the Martian atmosphere and sequestering it in the crust as carbonates, and as a possible origin for iron oxide-rich spectral units on the floors of some rifted basins (chasmata). There are numerous examples of Martian channels formed by discharges of subsurface water near potential magmatic heat sources, and hydrothermal processes have also been proposed as a mechanism for aquifer recharge needed to sustain long term erosion of sapping channels. The following geological settings have been identified as targets for ancient hydrothermal systems on Mars: channels located along the margins of impact crater melt sheets and on the slopes of ancient volcanoes; chaotic and fretted terranes where shallow subsurface heat sources are thought to have interacted with ground ice; and the floors of calderas and rifted basins (e.g. chasmata). On Earth, such geological environments are often a locus for hydrothermal mineralization. But we presently lack the mineralogical information needed for a definitive evaluation of hypotheses. A preferred tool for identifying minerals by remote sensing methods on Earth is high spatial resolution, hyperspectral, near-infrared spectroscopy, a technique that has been extensively developed by mineral explorationists. Future efforts to explore Mars for ancient hydrothermal systems would benefit from the application of methods developed by the mining industry to look for similar deposits on Earth. But Earth-based exploration models must be adapted to account for the large differences in the climatic and geological history of Mars. For example, it is likely that the early surface environment of Mars was cool, perhaps consistently below freezing, with the shallow portions of hydrothermal systems being dominated by magma-cryosphere interactions. Given the smaller gravitational field, declining atmospheric pressure, and widespread, permeable megaregolith on Mars, volatile outgassing and magmatic cooling would have been more effective than on Earth. Thus, hydrothermal systems are likely to have had much lower average surface temperatures than comparable geological settings on Earth. The likely predominance of basaltic crust on Mars suggests that hydrothermal fluids and associated deposits should be enriched in Fe, Mg, Si and Ca, with surficial deposits being dominated by lower temperature, mixed iron oxide and carbonate mineralogies.

Review↗

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗

Slope-Aspect Induced Climate Differences Influence How Water Is Exchanged Between the Land and Atmosphere

Cross-slope climate differences in the midlatitudes are ecologically important, and impact vegetation-mediated water balance between the earth surface and the atmosphere. We made high-resolution in situ observations of air temperature, relative humidity, soil moisture, insolation, and sap velocity observations on 14 Pacific madrone trees (Arbutus menziesii) spanning adjacent north and south slopes at the University of California's Angelo Coast Range Reserve. To understand the cross-slope response of sap velocity, a proxy for transpiration, to microclimate, we modeled the sap velocity on each slope using a transpiration model driven by ambient environment and parameterized with a Markov Chain Monte Carlo parameter estimation process. The results show that trees on opposing slopes do not follow a shared pattern of physiological response to transpiration drivers. This means that the observed sap velocity differences are not due entirely to observed microclimate differences, but also due to population-level physiological differences, which indicates acclimation to inhabited microclimate. While our present data set and analytical tools do not identify mechanisms of acclimation, we speculate that differing proportions of sun-adapted and shade-adapted leaves, differences in stomatal regulation, and cross-slope root zone moisture differences could explain some of the observed and modeled differences.

54 ENVIRONMENTAL SCIENCES↗

Landform Erosion and Volatile Redistribution on Ganymede and Callisto

We have been modeling landscape evolution on the Galilean satellites driven by volatile transport. Our work directly addresses some of the most fundamental issues pertinent to deciphering icy Galilean satellite geologic histories by employing techniques currently at the forefront of terrestrial, martian, and icy satellite landscape evolution studies [e.g., 1-6], including modeling of surface and subsurface energy and volatile exchanges, and computer simulation of long-term landform evolution by a variety of processes. A quantitative understanding of the expression and rates of landform erosion, and of volatile redistribution on landforms, is especially essential in interpreting endogenic landforms that have, in many cases, been significantly modified by erosion [e.g., 7-9].

Moore, Jeffrey Morgan↗

An Eight-County Appraisal of the San Andres Residual Oil Zone (ROZ) 'Fairway' of the Permian Basin

This document assesses the size and distribution of the oil resource in an eight-county portion of the San Andres residual oil zone (ROZ) fairway in the Permian Basin. Using well log analysis, the study divided the ROZ into partitions, developed a database of geologic properties for each partition and estimated the oil in place in each partition. The study then used the FE/NETL CO 2 Prophet Model to estimate the oil that could be produced and CO 2 stored from applying CO 2 EOR to each partition. With the FE/NETL Onshore CO 2 EOR Cost Model, the commercially viable partitions were identified as those that could produce oil and store CO 2 economically at an oil price of $75/STB. The remaining geologically viable partitions could use CO 2 EOR to store CO 2 with production of oil to offset some of the costs. The database used for this study, and a report describing the database is available on NETL's website under the Collection Name Eight-County San Andres ROZ Appraisal.

02 PETROLEUM↗

A Four-County Appraisal of the San Andres Residual Oil Zone (ROZ) 'Fairway' of the Permian Basin

This document assesses the size and distribution of the oil resource in a four-county portion of the San Andres residual oil zone (ROZ) fairway in the Permian Basin. Using well log analysis, the study divided the ROZ into partitions, developed a database of geologic properties for each partition and estimated the oil in place in each partition. The study then used the FE/NETL CO 2 Prophet Model to estimate the oil that could be produced and CO 2 stored from applying CO 2 EOR to each partition. With the FE/NETL Onshore CO 2 EOR Cost Model, the commercially viable partitions were identified as those that could produce oil and store CO 2 economically at an oil price of $75/STB. The remaining geologically viable partitions could use CO 2 EOR to store CO 2 with production of oil to offset some of the costs. The database used for this study, and a report describing the database is available on NETL's website under the Collection Name: Four-County San Andres ROZ Appraisal.

02 PETROLEUM↗

Advancing wildlife connectivity in land use planning: a case study with four-toed salamanders

Stable habitat connections that wildlife can safely traverse are essential to biodiversity conservation and healthy ecosystems. We developed high-resolution landscape connectivity models to predict resistance to movement by a threatened wetland-obligate amphibian, the four-toed salamander ( Hemidactylium scutatum ), and identified priority management areas on the 13,000-ha Department of Energy Oak Ridge Reservation (ORR) from 2019 to 2022. We developed a resistance surface based on aerial light detection and ranging data (LiDAR), >30 years of field-based mapping of forest, hydrologic, and geologic features, and contemporary population surveys, alongside derived predictors at <1-m resolution. We then modeled predicted movement corridors using a circuit theory-based modeling approach. We worked closely with land management and natural resources personnel to integrate ecological modeling with broader land use priorities, monetary costs, and feasibility. We identified important terrestrial and aquatic areas on ORR and simulated management scenarios to promote stable connections for four-toed salamanders. This approach allowed us to narrow down a list of 438 potential habitat manipulation sites to 10 sites where open-bottomed culverts and buffers could be implemented. This smaller-scale restoration approach produced a similar increase in landscape connectivity while costing <20% of a larger-scale approach based on barrier removal. We successfully identified feasible, cost-effective management strategies that integrated knowledge from a variety of sources. In conclusion, we offer a strategy that permitted integration of wildlife management goals into infrastructure upgrades wherein wildlife was not an initial consideration.

60 APPLIED LIFE SCIENCES↗

Modeling Flow and Transport in Fractured Rock using Machine Learning [Slides]

Fractured systems in the subsurface play a role in many natural and engineered applications such as geologic carbon sequestration, hydraulic fracturing and underground nuclear test detection. Structural information (fracture size, orientation, etc.) plays a key role in governing the dominant physics for these systems but can only be known statistically. Traditional modeling approaches either ignore or idealize structural information at these larger scales because we lack a computational framework that utilizes it in its entirety. The work presented here integrates computational physics, machine learning and graph theory to make a paradigm shift from computationally intensive high-fidelity models to coarse-scale graphs without loss of critical structural information. We exploit the underlying discrete structure of fracture networks in systems considering flow through fractures and fracture propagation. We demonstrate that compact graph representations require significantly fewer degrees of freedom (dof) to capture micro-fracture information and further accelerate these models with Machine Learning.

58 GEOSCIENCES↗

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗