Search NASA⌕ Search

SEARCH · Search NASA

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 343 records · Page 19

Seasonal mixing from intermittent flow drives concentration-discharge behaviour in a stream affected by coal mine drainage

Abandoned mining operations continue to severely degrade many ecosystems worldwide by releasing acidic water and/or heavy metals into surface and groundwater. Contaminant concentrations in affected streams vary with discharge in patterns that reflect both geochemical reactions and variable mixing of contaminated and non-contaminated waters. However, controls on concentration-discharge (C-Q) patterns remain unclear, particularly for constituents that experience changing solubility across redox and pH gradients. Understanding the C-Q behaviour of contaminants aids in predicting both downstream transport and effects on aquatic life under variable flow. Here, we examined the C-Q behaviours of non-reactive (Na, K, Ca, Mg, Cl - ) and reactive (Fe, Mn, Al, H + , SO 4 2- ) solutes in a stream contaminated with acid mine drainage in northeastern Ohio, USA. Concentration-discharge patterns at the watershed outlet primarily reflected mixing of contaminated baseflow with intermittent inputs of high pH water draining from a passive limestone treatment system into the stream. The treatment system acted as an ephemeral tributary that mitigated contamination in the stream by diluting solutes, raising pH, and driving metal precipitation, but only when flow was present during wet seasons. Consequently, AMD-derived reactive solutes (H + , Fe, Mn, Al) decreased with increasing stream discharge while relatively conservative solutes (e.g., Ca, Mg, K, Na) decreased only slightly or were chemostatic. This study highlights both the unique C-Q patterns of reactive solutes when compared to those of non-reactive solutes and the potential for intermittent streams to control C-Q behaviour in headwater catchments.

58 GEOSCIENCES↗

Vegetation water sources in California's Sierra Nevada (USA) are young and change over time, a multi-isotope ( δ 18 O, δ 2 H, 3 H) tracer approach.

Sierra Nevada forests transpire a significant amount of California's water resources, sparking interest in applying forest management to improve California's water supply. Determining the source water of evapotranspiration enables forest managers to make informed decisions. To this end, a significant interest in critical zone science is to develop new methods to work across time scales to predict subsurface water storage and use. Here, forest vegetation accessed young water and switched sources depending on availability, suggesting that forest drought vulnerability may depend on the range of water sources available (rain, snowmelt and deeply stored water). This finding also suggests that changes in transpiration rates may have immediate effects on water sources in close proximity to vegetation, and delayed effects on storage and runoff. New δ 18 O, δ 2 H and 3 H data were used to track precipitation, runoff, evapotranspiration and storage through the critical zone seasonally, including seasons where evapotranspiration and snowmelt were in phase (winter snowmelt) and out of phase (seasonally dry summer). The main source of this headwater catchment's runoff is derived from its meadow saturated zone water, which was dominated by snowmelt. Water that originated as snowmelt contributed to transpiration, unless other sources, such as recent rain, became available. In cases where xylem δ 18 O and δ 2 H signatures matched those of deeper saturated zone water, 3 H data showed that xylem water was distinctly younger than the deep saturated zone water. During 2016, which experienced relatively normal snowpack in winter and seasonally dry summer conditions, mean summer saturated zone water and vegetation water were similar in δ 18 O, –12.4 ± 0.04 ‰ and – 12.5 ± 0.3 ‰, respectively, but were distinctly different in 3 H, 5.5 ± 0.2 pCi/L and 13.7 ± 1.1 pCi/L, respectively. While δ 18 O shows that vegetation and meadow saturated zone water have similar origins, 3 H shows they have dissimilar ages.

58 GEOSCIENCES↗

Chemostatic concentration–discharge behaviour observed in a headwater catchment underlain with discontinuous permafrost

Concentration–discharge dynamics were evaluated in a small (~ 2.25 km2) headwater catchment underlain with discontinuous permafrost on the Seward Peninsula of western Alaska. A large storm, during which 48 mm of rain fell over a 24-h period, enabled the evaluation of solute concentration–discharge response to a sizeable hydrological event, while water stable isotopes enabled an appraisal of the contributions of event water. Under normal catchment conditions, chemostatic behaviour was observed for solutes typically derived from mineral weathering (e.g. calcium, magnesium, sodium and silica). The chemostatic behaviour observed for most solutes under normal catchment conditions indicated that catchment storage and residence times are sufficiently long for many solute generating reactions to approach equilibrium. Following the storm however, most solutes exhibited dilutive and highly variable behaviour. This likely indicated the exceedance of a discharge threshold where chemostatic behaviour could no longer be maintained for most solutes. Dissolved organic carbon and silica were the only solutes monitored to exhibit chemostatic behaviour during all time periods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying subsurface parameter and transport uncertainty using surrogate modelling and environmental tracers

Here, we combine physics-based groundwater reactive transport modelling with machine-learning techniques to quantify hydrogeological model and solute transport predictive uncertainties. We train an artificial neural network (ANN) on a dataset of groundwater hydraulic heads and 3 H concentrations generated using a high-fidelity groundwater reactive transport model. Using the trained ANN as a surrogate model to reproduce the input–output response of the high-fidelity reactive transport model, we quantify the posterior distributions of hydrogeological parameters and hydraulic forcing conditions using Markov chain Monte Carlo calibration against field observations of groundwater hydraulic heads and 3 H concentrations. We demonstrate the methodology with a model application that predicts Chlorofluorocarbon-12 (CFC-12) solute transport at a contaminated field site in Wyoming, United States. Our results show that including 3 H observations in the calibration dataset reduced the uncertainty in the estimated permeability field and infiltration rates, compared to calibration against hydraulic heads alone. However, predictive uncertainty quantification shows that CFC-12 transport predictions conditioned to the parameter posterior distributions cannot reproduce the field measurements. We found that calibrating the model to hydraulic head and 3 H observations results in groundwater mean ages that are too large to explain the observed CFC-12 concentrations. The coupling of the physics-based reactive transport model with the machine-learning surrogate model allows us to efficiently quantify model parameter and predictive uncertainties, which is typically computationally intractable using reactive transport models alone.

58 GEOSCIENCES↗

Seasonal and tidal variations in hydrologic inputs drive salt marsh porewater nitrate dynamics

Salt marshes remove terrestrially derived nutrients en route to coasts. While these systems play a critical role in improving water quality, we still have a limited understanding of the spatiotemporal variability of biogeochemically reactive solutes and processes within salt marshes. We implemented a high-frequency sampling system to monitor sub-hourly nitrate (NO$^{–}_{3}$) concentrations in salt marsh porewater at Elkhorn Slough in central California, USA. We instrumented three marsh positions along an elevation gradient subjected to different amounts of tidal inundation, which we predicted would lead to varied biogeochemical characteristics and hydrological interactions. At each marsh position, we continuously monitored porewater NO$^{–}_{3}$ concentrations at depths of 10, 30, and 50 cm and porewater levels measured at 70 cm depth over seven deployments of ~10 days each that spanned seasonal wet/dry periods common to Mediterranean climates. We quantified tidal event hysteresis between NO$^{–}_{3}$ and water level to understand how NO$^{–}_{3}$ concentrations and sources fluctuate across tidal cycles. In dry periods, the NO$^{–}_{3}$-porewater level relationship indicated that the NO$^{–}_{3}$ source was likely estuarine surface water that flooded the transect during high tides and the salt marsh was a NO$^{–}_{3}$ sink. In wet periods, the NO$^{–}_{3}$-porewater level relationship suggested the salt marsh was a source of NO$^{–}_{3}$. Finally, these findings suggest that tidal and seasonal hydrologic fluxes together control NO$^{–}_{3}$ porewater dynamics and export and influence ecological processes in coastal environments.

54 ENVIRONMENTAL SCIENCES↗

How Rains and Floods Become Groundwater: Understanding Recharge Pathways With Stable and Cosmogenic Isotopes

Anthropogenic climate change leads to increased precipitation intensity and exacerbated droughts in California, challenging the reliability and drought resiliency of water supply. Storing floodwater underground via managed aquifer recharge can mitigate these effects through direct infiltration or streambed infiltration. Seasonally dry streams (arroyos) already play an important part in managing groundwater recharge to the Livermore basin (CA). Understanding how, when and where stormwater and arroyo water infiltrate is critical to effectively utilise this strategy. To track water from recent storms (water year 2022–2023, WY23) into the Livermore Valley Groundwater Basin, we analysed stable water isotopes (δ 18 O and δ 2 H) in combination with naturally occurring radioactive isotopic tracers, sulphur-35 ( 35 S, t ½ = 87 days) and tritium ( 3 H, t ½ = 12.3 years). By comparing measurements of δ 18 O, 35 S and 3 H in arroyos to precipitation and groundwater, we classified the relative age and identified source of recharge to 16 wells near two arroyos. Two wells contained water with recent recharge (from WY23) from local precipitation. One well had recent recharge from variable (precipitation and imported water) sources. One well contained imported water recharge. Three wells contained water from mixed recent and older (pre-WY23) waters, from local precipitation sources. Two wells contained recent recharge from local mine settling ponds. Seven wells had older recharge from local precipitation sources. This combination of isotopes allows us to delineate where local and imported water recharges in this highly managed basin and identify locations where managed aquifer recharge is contributing to rapid groundwater infiltration. Our combined interpretation of isotopic water ages and sources in the context of land use shows that local infiltration of precipitation in open spaces is an important recharge mechanism, in addition to the managed arroyo recharge. Finally, a broader familiarity with 35 S will enable more extensive research on the infiltration of urban floodwaters.

58 GEOSCIENCES↗

Novel color center platforms enabling fundamental scientific discovery

Abstract Color centers are versatile systems that generate quantum light, sense magnetic fields and produce spin‐photon entanglement. We review how these properties have pushed the limits of fundamental knowledge in a variety of scientific disciplines, from rejecting local‐realistic theories to sensing superconducting phase transitions. In the light of recent progress in material processing and device fabrication, we identify new opportunities for interdisciplinary fundamental discoveries in physics and geochemistry. image

24 POWER TRANSMISSION AND DISTRIBUTION↗

Climatic controls on soil and saprock nitrogen distribution and persistence in the Sierra Nevada

Nitrogen (N) is an essential nutrient in soil that regulates plant growth, terrestrial sequestration of atmospheric carbon dioxide, and persistence of organic compounds. However, major knowledge gaps remain about how climate change may impact N accumulation and persistence, especially in deep soil and saprock (friable weakly weathered bedrock). Our objective was to understand how climate impacts the accumulation of N in soil and saprock and how climate impacts soil N distribution, persistence, and mechanisms of N persistence. We investigated N concentration in bulk soil and density fractions. We estimated N persistence along a bio-climatic sequence—sites range from a low-elevation oak savannah, mid-elevation pine-oak/mixed-conifer forest, to a high-elevation subalpine forest—in the southern Sierra Nevada in California. A combination of radiocarbon and elemental composition measurements along with a first-order kinetic model was used. The N concentration in the bulk soil and density fractions declined with depth, and there was a relatively greater mineral-associated heavy fraction (HF) N in deeper samples. The cooler/wetter mixed conifer site held 37% of profile N in saprock, which was greater than that of the entire soil profile at the drier/hotter oak savannah. The majority of N in soil, which was in the HF, was not influenced by climate proxies tested. However, both unprotected and occluded fractions of N were strongly influenced by climate. Soil N mean residence time (MRT) showed that drier/hotter climates have a shorter MRT, compared to mid-elevation sites with cooler/wetter climates. The effect of climate on deep saprock N storage might be indirect, primarily through climatic influence on the thickness of saprock. Overall, our findings suggest the mineral-associated HF N pool will not be vulnerable to changes in climate and will continue to contribute to the persistent soil N pool. The amount of topsoil and subsoil unprotected and occluded N can be explained by gross primary productivity and mean annual precipitation indicating that changes in climate can influence N partitioning. N stored in deep soil and saprock may be less vulnerable to climate than N stored in drier/hotter climates with less deeply stored N. It is critical to dig deeper to understand terrestrial ecosystems’ response to climate.

58 GEOSCIENCES↗

A stable isotope record of late Quaternary hydrologic change in the northwestern Brooks Range, Alaska (eastern Beringia)

A submillennial-resolution record of lake water oxygen isotope composition (δ 18 O) from chironomid head capsules is presented from Burial Lake, northwest Alaska. The record spans the Last Glacial Maximum (LGM; ~20–16k cal a bp) to the present and shows a series of large lake δ 18 O shifts (~5‰). Relatively low δ 18 O values occurred during a period covering the LGM, when the lake was a shallow, closed-basin pond. Higher values characterize deglaciation (~16–11.5k cal a Bp) when the lake was still closed but lake levels were higher. A rapid decline between ~11 and 10.5k cal a Bp indicates that lake levels rose to overflowing. Lake δ 18 O values are interpreted to reflect the combined effects of changes in lake hydrology, growing season temperature and meteoric source water as well as large-scale environmental changes impacting this site, including opening of the Bering Strait and shifts in atmospheric circulation patterns related to ice-sheet dynamics. The results indicate significant shifts in precipitation minus evaporation across the late Pleistocene to early Holocene transition, which are consistent with temporal patterns of vegetation change and paludification. This study provides new perspectives on the paleohydrology of eastern Beringia concomitant with human migration and major turnover in megafaunal assemblages.

54 ENVIRONMENTAL SCIENCES↗

Early glacier advance in New Zealand during the Antarctic Cold Reversal

ABSTRACT Glacial landscapes preserve records of past climate change. Investigating the glacier–climate system over the Late Quaternary provides information about past climate change and context for present‐day glacier response to climate warming. Using 28 beryllium‐10 ( 10 Be) surface exposure dates and snowline reconstructions, we present glacier fluctuations and climate changes for the Antarctic Cold Reversal in the Ahuriri River catchment, Southern Alps of New Zealand (44°7′50″S, 169°38′29″E). Prominent terminal and lateral moraine features from the upper right tributary of the Ahuriri River valley have exposure ages of 14.5 ± 0.3, 13.6 ± 0.3 and 12.6 ± 0.2 ka, suggesting retreat of the glacier during the Antarctic Cold Reversal. Maximum elevation of lateral moraines (MELM) and accumulation area ratio (AAR) suggest snowline elevations at these ages were ≤700, ≤630 and ~360 m lower than today, respectively. This equates to air temperatures ≤3.9, ≤3.5 and 2.3 ± 0.7 °C lower than today (1981–2010), assuming no changes in past precipitation. Ice‐sculpted bedrock surfaces bound by a lateral moraine at nearby Canyon Creek have an age of 13.1 ± 0.3 ka, indicating the moraine correlates with those in the Ahuriri upper right tributary. MELM and AAR reconstructions from the Canyon Creek suggest that snowline elevations at 14.5–13.6 ka were ≤500 or ~380 m lower than today, corresponding to air temperatures ≤2.8 or 2.4 ± 0.7 °C lower than the present‐day (1981–2010). Our results provide insight into the structure of the Antarctic Cold Reversal in the Southern Alps, showing that the largest glacier advance occurred at the start of this interval at c . 14.5 ± 0.3 ka and was followed by gradual retreat. We hypothesize that the early cooling and glacier readvance in New Zealand at the onset of the Antarctic Cold Reversal were triggered by a latitudinal shift of the Southern Hemisphere westerly wind belt.

58 GEOSCIENCES↗

Effect of Particle Size on Raman Signal Strength of Silicate Minerals

Understanding the effects of particle size is necessary for quantifying minerals in mixtures using Raman spectroscopy. Raman signal intensity is evaluated using six common silicate minerals (two olivines, two pyroxenes, and two feldspars) at 10 particle size ranges. For olivines and feldspars, the highest peak intensities are observed in samples with 38–63 and 63–106 μm particle sizes. There is no such consistent trend for the pyroxene samples, although the overall low signal strength complicates those measurements. In conclusion, Raman spectra of samples with varying particle sizes appear to be influenced by two competing effects: scattering from particle boundaries and effective sampling volume.

58 GEOSCIENCES↗

Establishing fluvial silicon regimes and their stability across the Northern Hemisphere

Abstract Fluvial silicon (Si) plays a critical role in controlling primary production, water quality, and carbon sequestration through supporting freshwater and marine diatom communities. Geological, biogeochemical, and hydrological processes, as well as climate and land use, dictate the amount of Si exported by streams. Understanding Si regimes—the seasonal patterns of Si concentrations—can help identify processes driving Si export. We analyzed Si concentrations from over 200 stream sites across the Northern Hemisphere to establish distinct Si regimes and evaluated how often sites moved among regimes over their period of record. We observed five distinct regimes across diverse stream sites, with nearly 60% of sites exhibiting multiple regime types over time. Our results indicate greater spatial and interannual variability in Si seasonality than previously recognized and highlight the need to characterize the watershed and climate variables that affect Si cycling across diverse ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of preservation protocols for oxygen-sensitive minerals within laminated aquatic sediments

Laminated sediments can record seasonal changes in sedimentation of material from anoxic waters, including minerals of the redox-sensitive elements Fe, Mn, and S that form under varying oxygen levels, mineral saturation conditions, and from microbial metabolism. However, preserving the oxygen-sensitive minerals for identification is challenging when preservation of the spatial arrangement of laminae is also required. In this study, we compare methods for embedding sedimentary materials from anoxic waters and sediments from Brownie Lake, Minnesota, USA for analysis of the speciation for Fe, Mn, and S using synchrotron-based X-ray absorption near edge spectroscopy (XANES). We found that acetone dehydration and resin replacement in a 100% N 2 glovebox successfully preserved the speciation of Fe and Mn minerals within laminated sediments. However, acetone removed some sulfur species from sediments, and epoxies contained sulfur species, which challenged identification of native sulfur species. Results from this study will aid researchers who are interested in spatial analysis of oxygen sensitive sediments, soils, or microbial mats in choosing a preservation method.

58 GEOSCIENCES↗

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗

An extended numerical manifold method for unsaturated soil-water interaction analysis at micro-scale

To investigate unsaturated soil-water interaction at micro-scale, this work extends the numerical manifold method (NMM) by incorporating a soil-water coupling model considering specific capillary water distribution and capillary force calculation. The soil skeleton is constructed by a soil skeleton generation algorithm with random polygons. To more realistically capture the interaction between soil grains and capillary water, a capillary mechanics-based geometric algorithm is proposed to iteratively calculate the capillary water distribution. The capillary forces corresponding to the capillary water distribution are calculated based on the Young-Laplace equation. The proposed capillary water solving framework is first verified by reproducing the soil-water characteristic curve and the capillary water distribution of an ideal contact-disk model against analytical solutions. To further validate the ability of the capillary water solving framework to predict hydraulic behavior of the real soil, a laboratory test on the Toyoura sand is reproduced numerically. Then an ideal direct shear test is performed to further validate the two-way soil-water coupling procedure, in which a comparison between the numerical and analytical results regarding the shear strength and matric suction is presented. Finally, microscopic hydraulic and compression tests are conducted on two soil specimens with the same porosity and mean grain diameter but different uniformity coefficients. The results elucidate that the extended method is a potential tool to explore unsaturated soil behaviors at micro-scale.

58 GEOSCIENCES↗

Spectral deferred correction methods for high-order accuracy in poroelastic problems

In this work, we investigate high-order accuracy in time integration by examining two operator splitting methods for poroelastic problems: the two-pass and the spectral deferred correction (SDC) methods. To enhance the order of accuracy, the two-pass method partitions a coupled operator symmetrically, whereas the SDC method corrects truncation errors by establishing an error equation. These high-order methods are applied to underlying solution strategies, i.e., monolithic, fixed-stress sequential, and undrained sequential methods. We observe that semi-discretized systems from spatial discretization have forms similar to those of index-1 differential algebraic equations (DAEs), causing order reduction against the two-pass method when it is used in conjunction with either the monolithic or sequential method. On the other hand, the SDC in conjunction with the monolithic method exhibits the desired second-order accuracy in poroelastic problems while increasing the order of accuracy for index-1 DAEs. However, the SDC in conjunction with either of the two sequential methods does not achieve the desired order of accuracy, and maintains first order because the flow equation for poroelasticity has an additional approximation associated with the volumetric strain rate term, which does not yield exactly the same forms as those of conventional DAEs. Thus, the monolithic SDC method can achieve higher-order accuracy, but may require higher computational costs because it involves solving matrix systems larger than those for the sequential methods.

02 PETROLEUM↗