Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional shift”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Lithium's low-temperature phase transitions: Insights into quantum lattice dynamics and superconductivity

The large lattice dynamics of lithium, driven by its low atomic mass, results in energetically similar structures and significant isotope effects under pressure, posing challenges to current theoretical models. Above 20 GPa and at low temperatures, lithium's electronic properties deviate from simple metallic behavior, with superconductivity emerging in a complex, pressure-dependent manner, alongside an unusual isotope effect. The structural phases of 7 Li reported under these conditions are inconsistent across studies, and the structures of 6 Li remain unexamined. These gaps limit our understanding of the effects of pressure on lithium's electronic properties and the role of quantum lattice effects on its structural behavior under pressure. Here, we integrate experimental and theoretical approaches to investigate the low-temperature structural phase boundaries in lithium isotopes. We map the structural phase diagram of 7 Li from 5 to 55 GPa and 15–75 K, identifying the sequence fcc → ℎ⁡R⁢1 → cI16. A pronounced isotope effect is observed, with 6 Li shifting the fcc → ℎ⁡R⁢1 phase boundary to lower pressures at 15 K. Density functional theory calculations further clarify how these structural changes affect superconducting properties, particularly emphasizing the role of the fcc → ℎ⁡R⁢1 transition in lithium's superconductivity. Furthermore, our findings offer insights into the unique behavior of lithium isotopes under pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Bayesian Exploration and Surrogate Emulation of Nonlinear Beam-Response Geometry in the LBNF Beamline

Next-generation long-baseline neutrino experiments aim to achieve multi-MW proton beam power while reducing accelerator-induced systematic uncertainties. At Fermilab, the LBNF beamline is designed for 1.2 MW operation with PIP-II and is upgradeable to 2.4 MW. DUNE will probe the three-flavor neutrino paradigm and search for CP violation, requiring precise neutrino-flux normalization and improved control of accelerator-related uncertainties. Within the LBNF beamline, the System for On-Axis Neutrino Detection (SAND) will constrain flux uncertainties using precision near-detector measurements, while the Muon Monitor System (MuMS) will provide beamline diagnostics sensitive to the proton beam, target, and horn configuration. However, the pion phase space relevant for DUNE depends simultaneously on many correlated parameters, including beam centroid, beam width, horn current and alignment, target position, optics shifts, and radiation-induced changes. Consequently, MuMS observables exhibit nonlinear and coupled responses that are difficult to characterize using traditional one-parameter scans. To address this challenge, we are developing a Bayesian Exploration framework coupled to physics-informed surrogate emulators trained on Geant4 beamline simulations. Gaussian-process emulators provide both fast predictions and uncertainty estimates, enabling adaptive selection of new simulation points in beam-parameter space. As an initial demonstration, we construct surrogate emulators for MuMS response observables using a verified simulation campaign spanning proton-beam steering conditions. The emulators reproduce the simulated dependence of MuMS centroid and gradient observables while providing predictive uncertainties, and serve as the foundation for future multidimensional exploration including beam width, horn current, and additional beamline parameters. This work establishes a framework for uncertainty-aware beam monitoring, adaptive simulation campaigns, and rapid beam-response inference for future DUNE operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗
DOE OSTI2020

Smaller than Expected Bright-spot Offsets in Spitzer Phase Curves of the Hot Jupiter Qatar-1b

We present Spitzer full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature—and intermediate surface gravity and orbital period—as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse depths of 0.21 ± 0.02% at 3.6 μm and 0.30 ± 0.02% at 4.5 μm, corresponding to dayside brightness temperatures of 1542{sub −31}{sup +32} K and 1557{sub −36}{sup +35} K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are 1117{sub −71}{sup +76} and 1167{sub −74}{sup +69} K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of 0.12{sub −0.16}{sup +0.14} and a moderate day–night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with Spitzer full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: 11° ± 7° at 3.6 μm and −4° ± 7° at 4.5 μm. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reliable bi-functional nickel-phosphate /TiO2 integration enables stable n-GaAs photoanode for water oxidation under alkaline condition

Abstract Hydrogen is one of the most widely used essential chemicals worldwide, and it is also employed in the production of many other chemicals, especially carbon-free energy fuels produced via photoelectrochemical (PEC) water splitting. At present, gallium arsenide represents the most efficient photoanode material for PEC water oxidation, but it is known to either be anodically photocorroded or photopassivated by native metal oxides in the competitive reaction, limiting efficiency and stability. Here, we report chemically etched GaAs that is decorated with thin titanium dioxide (~30 nm-thick, crystalline) surface passivation layer along with nickel-phosphate (Ni-Pi) cocatalyst as a surface hole-sink layer. The integration of Ni-Pi bifunctional co-catalyst results in a highly efficient GaAs electrode with a ~ 100 mV cathodic shift of the onset potential. In this work, the electrode also has enhanced photostability under 110 h testing for PEC water oxidation at a steady current density J ph > 25 mA·cm −2 . The Et-GaAs/TiO 2 /Ni-Pi║Ni-Pi tandem configuration results in the best unassisted bias-free water splitting device with the highest J ph (~7.6 mA·cm −2 ) and a stable solar-to-hydrogen conversion efficiency of 9.5%.

08 HYDROGEN↗
DOE OSTI2024

Processing and characterization of the homologous Zr x Ta 2 O 2x+5 series

In this study, Zr x Ta 2 O 2x+5 (ZTOx) was systematically studied to determine ideal solid-state synthesis and pressureless sintering conditions. For Zr 6 Ta 2 O 17 (ZTO6) and other compositions with values of x in ZTOx the optimal synthesis temperature of 1100 °C with 1 h dwell time was optimal. XRD phase analysis of the homologous series indicated that for ZTOx, x = 5, 6, 7, and 8 were within the phase stability field while x = 4 and x = 9 were outside the phase stability field. Unit cell lattice parameters changed anisotropically with Zr:Ta ratio with an overall <0.1% change in unit cell volume across the solid solution. Some Raman modes red shifted while others blue shifted, indicating competing tensile and compressive strains. Sintering of ZTOx required a progressively higher temperature as Zr-content increased despite the same starting powder grain sizes and synthesis temperature. Higher Vickers hardness for Zr-rich (12.7–12.8 GPa) than Ta-rich Zr x Ta 2 O 2x+5 (11.6–11.7 GPa) disappeared when accounting for different sintering temperatures. Indentation toughness was similar for all compositions (1.2–1.4 MPa m 1/2 ). Properties of Zr x Ta 2 O 2x+5 series indicate tunability of Zr:Ta composition ratio while maintaining the structure and room temperature mechanical properties comparable to Hf 6 Ta 2 O 17 and 8 mol% yttria stabilized zirconia ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Using taxa-based approaches to delineate stream macroinvertebrate assemblage responses to stressor gradients in modified alluvial agroecosystems

Alluvial plain landscapes are some of the most agriculturally productive lands in the world but often have modified stream ecosystems due to cultivation history. This context requires consideration when establishing water quality management goals. We analyzed state water quality databases to demonstrate that Mississippi Alluvial Plain (MAP) ecoregion streams have elevated specific conductivity (SC) and nutrients and lower macroinvertebrate local and regional taxa pools compared to streams in other ecoregions, potentially reducing the efficacy of traditional biomonitoring approaches within the region. To overcome these challenges, we used threshold indicator taxa analysis (TITAN) to compare macroinvertebrate assemblage responses to water quality gradients among ecoregions in Mississippi. We identified individual taxa and assemblage-level responses to increasing water quality degradation in MAP streams. Observed responses occurred at higher concentrations for SC, total organic carbon (TOC) and total phosphorus (TP), but not total nitrogen (TN) relative to other ecoregions. These responses appeared to be driven by a large proportion of indicator taxa considered tolerant or unresponsive in other ecoregions, responding negatively to increasing water quality stressors in MAP streams. Our observed assemblage-level stressor responses to WQ gradients in MAP streams demonstrate shifting tolerance in highly altered ecosystems may require adjustments to recovery expectations but also provide useful measures for monitoring improvements in regional water quality. For example, our observed macroinvertebrate assemblage response to increasing TP identified a management goal similar to guidance based on distributional analysis of water quality data within the MAP ecoregion (0.11 vs 0.128 mg L –1 ) and thus provide some biological basis for previously identified nutrient goals for the region. Our approach can guide and monitor success of nutrient reduction efforts in MAP watersheds and other alluvial plain agroecosystems where reference conditions do not exist, and local and regional taxa pools are less diverse and may not support full recovery of ecological assemblages. While our results are promising, they should also be compared with more sensitive and less habitat-limited biological assemblages (e.g., algae or bacteria) to better understand complex ecological responses to best management practices designed to increase sustainability of high production agricultural regions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Observing differential spin currents by resonant inelastic X-ray scattering

Controlling spin currents, that is, the flow of spin angular momentum, in small magnetic devices, is the principal objective of spin electronics, a main contender for future energy-efficient information technologies. A pure spin current has never been measured directly because the associated electric stray fields and/or shifts in the non-equilibrium spin-dependent distribution functions are too small for conventional experimental detection methods optimized for charge transport. Here we report that resonant inelastic X-ray scattering (RIXS) can bridge this gap by measuring the spin current carried by magnons—the quanta of the spin wave excitations of the magnetic order—in the presence of temperature gradients across a magnetic insulator. This is possible due to the sensitivity of the momentum- and energy-resolved RIXS intensity to minute changes in the magnon distribution under non-equilibrium conditions. Furthermore, we use the Boltzmann equation in the relaxation time approximation to extract transport parameters, such as the magnon lifetime at finite momentum, essential for the realization of magnon spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Setting the H-mode pedestal structure: variations of particle source location using gas puff and pellet fueling

Experiments in DIII-D show that a particle source location inside the top of the H-mode pedestal (pellets) maintains a higher pedestal pressure than an edge source (gas fueling) through a widening of the electron temperature pedestal with reduction of the temperature gradient. The effect of these two fueling schemes on the H-mode pedestal structure was examined in DIII-D by comparing controlled pellet-fueled and gas-fueled discharges across a fueling scan up to 40 torr 1 s –1 . High resolution electron profiles reveal that gas fueling lowers the pedestal pressure as the density profile shifts radially outwards and the separatrix density increases, while pellet fueling maintains a constant pedestal pressure. The neutral source locations from pellets and gas are determined with the PELLET and UEDGE codes, respectively, and quantify the particle source localization. Pellets provide significant ionization inside the pedestal top while gas puffing localizes ionization in the scrape-off-layer and pedestal foot, broadly consistent with the density profile structure influenced by the source. ELMs are observed to increase in frequency and reduce impurity content as fueling is increased. Stability analysis with ELITE shows that both conditions are near the type-I ELM corner of the peeling–ballooning stability diagram, which is altered significantly by the introduction of pellets. Since transport mechanisms are not observed to change substantially with particle source location, wider pedestals allow the pellet-fueled discharges to retain higher pedestal temperatures at similar pedestal densities. EPED1 is tested to capture the pedestal pressure, under-predicting the height with pellets and over-predicting the height with gas by ~15%. Overall, these results have important implications for future reactors where pellet fueling will be the primary particle source due to an opaque scrape-off-layer by showing that the neutral source location plays a role in setting the structure of the H-mode pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimizing Magnetic Force Microscopy Resolution and Sensitivity to Visualize Nanoscale Magnetic Domains

Magnetic force microscopy (MFM) enables mapping local magnetic fields across a sample surface with nanoscale resolution. To perform MFM, an atomic force microscopy (AFM) probe whose tip has been magnetized vertically (i.e., perpendicular to the probe cantilever) is oscillated at a fixed height above the sample surface. The resultant shifts in the oscillation phase or frequency, which are proportional to the magnitude and sign of the vertical magnetic force gradient at each pixel location, are then tracked and mapped. Although the spatial resolution and sensitivity of the technique increases with decreasing lift height above the surface, this seemingly straightforward path to improved MFM images is complicated by considerations such as minimizing topographical artifacts due to shorter range van der Waals forces, increasing the oscillation amplitude to further improve sensitivity, and the presence of surface contaminants (in particular water due to humidity under ambient conditions). In addition, due to the orientation of the probe's magnetic dipole moment, MFM is intrinsically more sensitive to samples with an out-of-plane magnetization vector. Here, high-resolution topographical and magnetic phase images of single and bicomponent nanomagnet artificial spin-ice (ASI) arrays obtained in an inert (argon) atmosphere glovebox with <0.1 ppm O 2 and H 2 O are reported. Further, optimization of lift height and drive amplitude for high resolution and sensitivity while simultaneously avoiding the introduction of topographical artifacts is discussed, and detection of the stray magnetic fields emanating from either end of the nanoscale bar magnets (~250 nm long and <100 nm wide) aligned in the plane of the ASI sample surface is shown. Likewise, using the example of a Ni-Mn-Ga magnetic shape memory alloy (MSMA), MFM is demonstrated in an inert atmosphere with magnetic phase sensitivity capable of resolving a series of adjacent magnetic domains each ~200 nm wide.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Impacts of directed evolution and soil management legacy on the maize rhizobiome

Domestication and agricultural intensification dramatically altered maize and its cultivation environment. Changes in maize genetics (G) and environmental (E) conditions increased productivity under high-synthetic-input conditions. However, novel selective pressures on the rhizobiome may have incurred undesirable tradeoffs in organic agroecosystems, where plants obtain nutrients via microbially mediated processes including mineralization of organic matter. Using twelve maize genotypes representing an evolutionary transect (teosintes, landraces, inbred parents of modern elite germplasm, and modern hybrids) and two agricultural soils with contrasting long-term management, here we integrated analyses of rhizobiome community structure, potential microbe-microbe interactions, and N-cycling functional genes to better understand the impacts of maize evolution and soil management legacy on rhizobiome recruitment. We show complex shifts in rhizobiome communities during directed evolution of maize (defined as the transition from teosinte to modern hybrids), with a larger effect of domestication (teosinte to landraces) than modern breeding (inbreds to hybrids) on rhizobiome structure and greater impacts of modern breeding on potential microbe-microbe interactions. Rhizobiome structure was significantly correlated with plant nutrient composition. Furthermore, plant biomass and nutrient content were affected by G x E interactions in which teosinte and landrace genotypes had better relative performance in the organic legacy soil than inbred and modern genotypes. The abundance of six N-cycling genes of relevance for plant nutrition and N loss pathways did not significantly differ between teosinte and modern rhizospheres in either soil management legacy. These results provide insight into the potential for improving maize adaptation to organic systems and contribute to interdisciplinary efforts toward developing resource-efficient, biologically based agroecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Systematic Characterization on Optical and Phonon Properties of Cs[Pb x Sn 1– x ]I 3 Alloyed Perovskites via First-Principles Modeling

Here, understanding the optical properties and phonon thermodynamic features of materials enables the determination of their viability for optoelectronic applications including optical-based gas sensors. Metal halide perovskites ABX 3 as well as their mixed alloys can be attractive in applications of this kind. In particular, the formation of alloyed perovskites A[B x B’ 1–x ]X 3 modulates the functional characteristics that can overtake the single-halide ABX 3 . To this end, we use the first-principles simulations to model the effect of B-site alloying on the optical and phonon properties of cubic Cs[Pb x Sn 1–x ]I 3 (x = 0.75, 0.50, and 0.25). We identified the systematic shifts of optical properties as a function of alloying concentration, such as dielectric constants, refraction indices, and absorption coefficients. It was also found that a linear correlation exists between the optical constants and the band gaps of these materials, which satisfies the Penn model. Modeling of phonons in Cs[Pb x Sn 1–x ]I 3 showed that dissipating the dynamically unstable phonons can be feasible at a carefully chosen alloying concentration, which are originally present in their single-halide Pb- and Sn-ends. Our findings pave a way to optimize the alloying conditions of the alloyed perovskite systems of interest for tuning their optical and phonon stability characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Altering the localization and toxicity of arsenic in rice grain

Previous work has shown that inorganic As localizes in rice bran whereas DMA localizes in the endosperm, but less is known about co-localization of As and S species and how they are affected by growing conditions. We used high-resolution synchrotron X-ray fluorescence imaging to image As and S species in rice grain from plants grown to maturity in soil (field and pot) and hydroponically (DMA or arsenite dosed) at field-relevant As concentrations. In hydroponics, arsenite was localized in the ovular vascular trace (OVT) and the bran while DMA permeated the endosperm and was absent from the OVT in all grains analyzed, and As species had no affect on S species. In pot studies, soil amended with Si-rich rice husk with higher DMA shifted grain As into the endosperm for both japonica and indica ecotypes. In field-grown rice from low-As soil, As localized in the OVT as arsenite glutathione, arsenite, and DMA. Results support a circumferential model of grain filling for arsenite and DMA and show Si-rich soil amendments alter grain As localization, potentially lessening risk to rice consumers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Seawater exposure causes hydraulic damage in dying Sitka-spruce trees

Sea-level rise is one of the most critical challenges facing coastal ecosystems under climate change. Observations of elevated tree mortality in global coastal forests are increasing, but the mechanism of salinity-stress induced tree mortality is currently unknown. We monitored progressive mortality and associated gas exchange and hydraulic shifts in Sitka-spruce (Picea sitchensis) trees located within a salinity gradient under an ecosystem-scale change of seawater exposure in Washington State, USA. Percentage of live foliated crown (PLFC) decreased and tree mortality increased with increasing soil salinity during the study period. A strong reduction in gas exchange and xylem hydraulic conductivity (Ks) occurred during tree death, with an increase of the percentage loss of xylem conductivity (PLC) and turgor loss point (ptlp). Hydraulic and osmotic shifts reflected that hydraulic function declined from seawater exposure and dying trees were unable to support osmotic adjustment. Constrained gas exchange was strongly related to hydraulic damage at both stem and leaf levels. Significant correlations between foliar sodium (Na+) concentration and gas exchange and key hydraulic parameters (Ks, PLC and ptlp ) suggest that cellular injury related to the toxic effects of ion accumulation impacted the physiology of these dying trees. This study provides evidence of toxic effects on cellular function that become manifest in all aspects of plant functioning, leading to unfavourable osmotic and hydraulic conditions.

Zhang, Hongxia↗
DOE OSTI2024

Analysis and Control of a Three-Port Converter with High-Frequency AC Link for Hybrid Unmanned/Autonomous Aerial Vehicles

In this study, a three-port resonant converter with a high-frequency isolation link for hybrid unmanned/autonomous aerial vehicle (UAV/AAV) applications is analyzed. The rectified engine-generator set output supplies the first input ports while the other is connected to a battery energy storage system. Input ports consist of half-bridge high-frequency inverters, whereas output port utilizes a semi-bridgeless active rectifier. This rectifier is formed by replacing rectifier low-side diodes with active switching devices. A phase-shift duty-cycle control is implemented for switching signals to regulate the power share and flow between ports. The proposed control system implemented achieves zero voltage switching (ZVS), reduces switching losses, and minimum current flow optimization ensures regulation in a wide range of output power. Compared to the existing bidirectional pulse width modulation (PWM) converters, the proposed system decreases the number of power conversion stages and provides centralized control without needing communications between the ports. The proposed control technique operating in closed-loop control achieves the power flow between ports under various load conditions. Finally, theoretical analysis demonstrates the proposed converter’s viability and laboratory experimental validations at 500 W full power.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Calibration of CDM-Based Creep Constitutive Model Using Accelerated Creep Test (ACT) Data

In conventional creep testing (CCT) a specimen is subject to constant load and temperature for a long-duration until creep rupture occurs. Conventional testing can be costly when considering the number of experiments needed to characterize the creep response of a material over a range of stress and temperature. To predict long-term creep-rupture properties, the time-temperature-stress superposition principle (TTSSP) approach has been employed where stress and/or temperature is applied at an elevated level; the result of which are extrapolated down to low stress and/or temperature conditions. These methods have been successful in predicting minimum-creep-strain-rate (MCSR) and stress-rupture (SR) but suffer from an inability to predict the creep deformation curve or account for changes in deformation mechanisms or aging that occurs at long-duration. An accelerated technique, termed the Stepped isostress method (SSM) allows the accelerated testing of materials to determine their creep deformation response. Unlike TTSSP tests, the SSM test employs a single specimen where the stress is periodically step increased until rupture. The SSM creep deformation curve is processed (time and strain shifted) to produce an accelerated creep deformation curve that represent the creep deformation curve at the initial stress level in SSM. A processing procedure for metals has yet to be developed. The research objective of this study is to develop a processing procedure for SSM test data using a creep-damage constitutive model. Triplicate SSM tests were conducted on Ni-based superalloy Inconel 718 at 650°C with stress being periodically increased until rupture. Triplicate CCT tests were conducted at the initial stress level of the SSM tests. The Sine-hyperbolic (Sinh) creep-damage model was employed in this study. The Sinh creep-damage constitutive model is based on coupled creep strain rate and damage evolution equations; where both rates are dependent on the current state of damage. Calibration is two-step: analytical and numerical optimization. Each stepped creep deformation curve is tackled quasi-analytically to determine MCSR and SR related material constants and accumulated damage. The damage accumulated at the end of each step was then passed onto subsequent steps to calibrate the MCSR, rupture prediction, and damage evolution. Numerical optimization was applied to optimize model constants involved in the creep strain constitutive equations in order to generate best-fitted Sinh creep deformation curves. The Sinh model predictions were compared to the SSM and CCT data. The Sinh model satisfactorily predicts the SSM data and thus the calibrated material constants provides a good estimate of rupture found in the CCT data. Calibration using SSM data reduces the number of tests needed to calibrate a model; significantly reducing costs. A single SSM test replaces numerous creep tests at different stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tool Chain for Simulations of Bi-Filar Coil Winding for Fast Quench Protection

The advancement of high-field magnets utilizing high-temperature superconductors (HTS) brings about complex challenges, especially in quench detection and protection. Traditional methods often fall short due to the inherently slow quench propagation in HTS materials. One promising approach to overcome this involves using a bifilar winding configuration, where two conductors are placed side by side. Under normal operation, they function in series, but during a quench event, they switch to an anti-parallel mode. This shift reduces the differential inductance of the coil to near zero, enabling rapid current oscillations through a capacitor discharge. The resulting high-frequency current flow leads to swift, uniform heating, triggering a full-coil quench within microseconds. Moreover, the strong mutual coupling between the two windings significantly reduces electrical noise in voltage measurements. In this work, we explore the viability of this concept by designing, constructing, and testing a REBCO bifilar racetrack coil in liquid nitrogen. We also present a validated simulation model that closely mirrors the coil's dynamic behavior under these conditions, aligning well with experimental observations.

Jayathilaka, Rehan [Northern Illinois U.]↗
DOE OSTI2021

A mechanistic understanding of oxygen isotopic changes in the Western United States at the Last Glacial Maximum

At the Last Glacial Maximum (LGM), records suggest drier conditions in the northwest United States and wetter conditions in the southwest United States relative to present-day as well as widespread changes in the isotopic composition of water. However, the mechanisms responsible for these changes remain ambiguous. Here, we explore differences in western United States hydroclimate between the LGM and preindustrial with a water isotope tracer enabled Earth System Model. We then use proxy forward models to compare simulated and recorded δ 18 O in speleothems. We find that the pattern of hydroclimate response in the western United States at the LGM relates to a combination of 1) increased frequency and southward shifted wintertime extratropical cyclones in the North Pacific, 2) greater rainout of moisture as it moves over the continent, and 3) reduced evaporation in the cooler LGM climate. The simulated lower δ 18 O of precipitation at the LGM relates predominantly to an increase in cool season moisture removal efficiency, with a secondary contribution from relatively more cool season precipitation. Both surface temperatures and North American ice sheets contribute to these hydroclimate changes at the LGM. Comparisons between δ 18 O from proxy forward models and speleothem records in the western United States show general agreement at the LGM, with increasing depletion moving towards the continental interior. Finally, this study highlights the similarities and differences between hydrologic and δ 18 O changes at the LGM and emphasizes the utility of model-proxy comparison for interpretation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗