Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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 1,693 records · Page 94

Grain Size Variations in the Murray Formation: Stratigraphic Evidence for Changing Depositional Environments in Gale Crater, Mars

The lowermost exposure of the Murray formation in Gale crater, Mars, was interpreted as sediment deposited in an ancient lake based on data collected by the Curiosity rover. Constraining the stratigraphic extent and duration of this environment has important implications for the paleohydrology of Gale. Insights into early Martian environments and paleofluid flow velocity can be obtained from grain size in rocks. Visual inspection of grain size is not always available for rocks investigated at field sites on Mars due to limited image coverage. But grain sizes can also be estimated from the Gini Index Mean Score, a grain-size proxy that uses point-to-point chemical variations in ChemCam Laser Induced Breakdown Spectroscopy data. Here, new Gini Index Mean Score results indicate that the Murray formation is dominated by mudstones with grains smaller than the spatial resolution of all rover cameras. Intervals of fine to coarse sandstone also are present, some of which are verified using observations of grain size and sedimentary structures in associated images. Overall, results demonstrate that most of the Murray consists of mudstone, suggesting settling of grains from suspension in low energy depositional environments such as lakes. Some of the mudstones contain desiccation cracks indicating periods of drying with a lowering of lake water level. However, beds and lenses of cross-bedded sandstones are common at specific intervals, suggesting episodes of fluvial and possibly eolian deposition. The persistence of lacustrine deposits interspersed with fluvial deposits suggests that liquid water was sustained on the Martian surface for tens of thousands to millions of years.

58 GEOSCIENCES↗

Spectral, Compositional, and Physical Properties of the Upper Murray Formation and Vera Rubin Ridge, Gale Crater, Mars

During 2018 and 2019, the Mars Science Laboratory Curiosity rover investigated the chemistry, morphology, and stratigraphy of Vera Rubin ridge (VRR). Using orbital data from the Compact Reconnaissance Imaging Spectrometer for Mars, scientists attributed the strong 860 nm signal associated with VRR to the presence of red crystalline hematite. However, Mastcam multispectral data and CheMin X-ray diffraction (XRD) measurements show that the depth of the 860 nm absorption is negatively correlated with the abundance of red crystalline hematite, suggesting that other mineralogical or physical parameters are also controlling the 860 nm absorption. Here, we examine Mastcam and ChemCam passive reflectance spectra from VRR and other locations to link the depth, position, and presence or absence of iron-related mineralogic absorption features to the XRD-derived rock mineralogy. Correlating CheMin mineralogy to spectral parameters showed that the ~860 nm absorption has a strong positive correlation with the abundance of ferric phyllosilicates. New laboratory reflectance measurements of powdered mineral mixtures can reproduce trends found in Gale crater. We hypothesize that variations in the 860 nm absorption feature in Mastcam and ChemCam observations of VRR materials are a result of three factors: (1) variations in ferric phyllosilicate abundance due to its ~800–1,000 nm absorption; (2) variations in clinopyroxene abundance because of its band maximum at ~860 nm; and (3) the presence of red crystalline hematite because of its absorption centered at 860 nm. We also show that relatively small changes in Ca-sulfate abundance is one potential cause of the erosional resistance and geomorphic expression of VRR.

58 GEOSCIENCES↗

Synergistic ground and orbital observations of iron oxides on Mt. Sharp and Vera Rubin ridge

Visible/short-wave infrared spectral data from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) show absorptions attributed to hematite at Vera Rubin ridge (VRR), a topographic feature on northwest Mt. Sharp. The goals of this study are to determine why absorptions caused by ferric iron are strongly visible from orbit at VRR, and to improve interpretation of CRISM data throughout lower Mt. Sharp. These goals are achieved by analyzing coordinated CRISM and in situ spectral data along the Curiosity Mars rover’s traverse. VRR bedrock within areas that have the deepest ferric absorptions in CRISM data also have the deepest ferric absorptions measured in situ. This suggests strong ferric absorptions are visible from orbit at VRR because of the unique spectral properties of VRR bedrock. Dust and mixing with basaltic sand additionally inhibit the ability to measure ferric absorptions in bedrock stratigraphically below VRR from orbit. There are two implications of these findings: (1) Ferric absorptions in CRISM data initially dismissed as noise could be real, and ferric phases are more widespread in lower Mt. Sharp than previously reported, (2) Patches with the deepest ferric absorptions in CRISM data are, like VRR, reflective of deeper absorptions in the bedrock. One model to explain this spectral variability is late-stage diagenetic fluids that changed the grain size of ferric phases, deepening absorptions. Curiosity’s experience highlights the strengths of using CRISM data for spectral absorptions and associated mineral detections, and the caveats in using these data for geologic interpretations and strategic path planning tools.

58 GEOSCIENCES↗

Boron and Lithium in Calcium Sulfate Veins: Tracking Precipitation of Diagenetic Materials in Vera Rubin Ridge, Gale Crater

The NASA Curiosity rover's ChemCam instrument suite has detected boron in calcium-sulfate-filled fractures throughout the sedimentary strata of Gale crater including Vera Rubin ridge. The presence of elevated B concentration provides insights into Martian subsurface aqueous processes. Here, we extend the data set of B in Ca-sulfate veins across Gale crater, comparing the detection frequency and relative abundances with Li. We report 33 new detections of B within veins analyzed between Sols 1548 and 2311 where detections increase in Pettegrove Point and Jura members, which form Vera Rubin ridge. The presence of B and Li in the Ca-sulfate veins is possibly due to dissolution of preexisting B in clays of the bedrock by acids or neutral water and redistribution of the elements into the veins. Elevated frequency of B detection in veins of Gale crater correlates with presence of dehydration features such as desiccation cracks, altered clay minerals and detections of evaporites such as Mg-sulfates and chloride salts in the host rocks. The increased observations of B also coincide with decreased Li concentration in the veins (average Li concentration of veins drops by ~15 ppm). Boron and Li have varying solubilities, and Li does not form salts as readily upon dehydration as B, causing it to remain in the solution. So the weak negative correlation between B and Li may reflect the crystallization sequence during dehydration on Vera Rubin ridge.

58 GEOSCIENCES↗

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year -1 (1Tg = 10 12 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH 4 year -1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of -34 to -46 Tg CH 4 year -1 , we estimated that the global net land methane emissions are 149–176 Tg CH 4 year -1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

58 GEOSCIENCES↗

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.

58 GEOSCIENCES↗

Influences of Hillslope Biogeochemistry on Anaerobic Soil Organic Matter Decomposition in a Tundra Watershed

In this work, we investigated rates and controls on greenhouse gas (CO 2 and CH 4 ) production in two contrasting water-saturated tundra soils within a permafrost-affected watershed near Nome, Alaska, United States. Three years of field sample analysis have shown that soil from a fen-like area in the toeslope of the watershed had higher pH and higher porewater ion concentrations than soil collected from a bog-like peat plateau at the top of the hillslope. The influence of these contrasting geochemical and topographic environments on CO 2 and CH 4 production was tested in soil microcosms by incubating both the organic- and mineral-layer soils anaerobically for 55 days. Additionally, nitrogen (as NH 4 Cl) was added to half of the microcosms to test potential effects of N limitation on microbial greenhouse gas production. We found that the organic toeslope soils produced more CO 2 and CH 4 , fueled by higher pH and higher concentrations of water-extractable organic C (WEOC). Our results also indicate N limitation on CO 2 production in the peat plateau soils but not the toeslope soils. Together these results suggest that the weathering and leaching of ions and nutrients from tundra hillslopes can increase the rate of anaerobic soil organic matter decomposition in downslope soils by (1) increasing the pH of soil porewater; (2) providing bioavailable WEOC and fermentation products such as acetate; and (3) relieving microbial N limitation through nutrient runoff. We conclude that the soil geochemistry as mediated by landscape position is an important factor influencing the rate and magnitude of greenhouse gas production in tundra soils.

58 GEOSCIENCES↗

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗

Numerically Relevant Timescales in the MG2 Microphysics Model

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.

58 GEOSCIENCES↗

Four Theories of the Madden‐Julian Oscillation

Abstract Studies of the Madden‐Julian Oscillation (MJO) have progressed considerably during the past decades in observations, numerical modeling, and theoretical understanding. Many theoretical attempts have been made to identify the most essential processes responsible for the existence of the MJO. Criteria are proposed to separate a hypothesis from a theory (based on the first principles with quantitative and testable assumptions, able to predict quantitatively the fundamental scales and eastward propagation of the MJO). Four MJO theories are selected to be summarized and compared in this article: the skeleton theory, moisture‐mode theory, gravity‐wave theory, and trio‐interaction theory of the MJO. These four MJO theories are distinct from each other in their key assumptions, parameterized processes, and, particularly, selection mechanisms for the zonal spatial scale, time scale, and eastward propagation of the MJO. The comparison of the four theories and more recent development in MJO dynamical approaches lead to a realization that theoretical thinking of the MJO is diverse and understanding of MJO dynamics needs to be further advanced.

58 GEOSCIENCES↗

Inertial Effects During the Process of Supercritical CO 2 Displacing Brine in a Sandstone: Lattice Boltzmann Simulations Based on the Continuum-Surface-Force and Geometrical Wetting Models

Inertial effects during the process of supercritical CO 2 displacing brine in porous media may not be negligible according to recent studies. Capturing the inertial effects of the physical CO 2 -brine system imposes a requirement on the grid resolution and viscosity to surface tension ratio for pore-scale simulations, which some commonly used simulators may not be able to meet. To fulfill the parameter requirement, we combine the continuum-surface-force based color-gradient lattice Boltzmann (LB) multiphase model and the geometrical wetting model and extend the model to 3D under the multiple-relaxation-time framework. We validate the model via simple benchmarks which show significant improvement over the traditional models. Furthermore, we perform 3D drainage simulations in a heterogeneous micromodel where the simulation result agrees well with experimental data, while our previous work fails to reproduce certain displacement patterns in the experiment due to the use of a traditional LB model that cannot fulfill the parameter requirement. Finally, we perform high-fidelity 3D drainage simulations to study the inertial effects in a Bentheimer sandstone sample. These findings show that stronger inertial effects generally help develop more CO 2 flow pathways for the same capillary number which results in higher CO 2 saturation, consistent with the micromodel results. The phenomena can be found in both low and high capillary number cases, indicating that the inertial effects are not dependent on the mean velocity. In addition, the change of the invasion patterns is not proportional to the change of inertial effects, thus exhibiting threshold behavior.

58 GEOSCIENCES↗

Impacts of Mixed-Wettability on Brine Drainage and Supercritical CO 2 Storage Efficiency in a 2.5-D Heterogeneous Micromodel

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO 2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO 2 (scCO 2 ) is used in enhanced oil recovery. Here, we performed scCO 2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from –8.1 to –4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO 2 saturation increase at both capillary fingering and crossover flow regimes (–8.1 ≤ logCa ≤ –6.1), the one heterogeneous wetting to scCO 2 results in ~0.09 saturation increase only at the crossover flow regime (–7.1 ≤ logCa ≤ –6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO 2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > –6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO 2 saturation and distribution suggest the importance of wettability on CO 2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

58 GEOSCIENCES↗

Simplified Green-Ampt Model, Imbibition-Based Estimates of Permeability, and Implications for Leak-off in Hydraulic Fracturing

Predicting water imbibition into porous materials is important in a wide variety of fields, and hydraulic fracturing of low permeability hydrocarbon reservoirs has emerged as an application that is imposing a large water footprint. Furthermore, reliable predictions of imbibition are needed to better manage water use, yet are challenging because of uncertainties in both the permeability and capillary pressure driving force. Here, this uncertainty is reduced through evaluating correlations between the permeability and the effective capillary pressure associated with the wetting front, P c,f . These correlations allow elimination of P c,f from the Green and Ampt equation and concentrate all uncertainties in fluxes on the effective permeability k . Over a wide range of k and porosities n , imbibition scales approximately with k 1/3 . Although Leverett k 1/4 scaling for predicting P c,f is shown to be inferior when tested with data spanning a wide range of n , it nevertheless predicted imbibition fairly well. From simple imbibition measurements, both the empirical and Leverett scaling approaches allow estimates of k that have root-mean-square deviations of about one order of magnitude relative to measurements that ranged over 10 orders of magnitude in k .

58 GEOSCIENCES↗