Search NASA⌕ Search

SEARCH · Search NASA

Results for “Space Science Data Systems”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Radiation-induced alteration of apatite on the surface of Mars: first in situ observations with SuperCam Raman onboard Perseverance

Abstract Planetary exploration relies considerably on mineral characterization to advance our understanding of the solar system, the planets and their evolution. Thus, we must understand past and present processes that can alter materials exposed on the surface, affecting space mission data. Here, we analyze the first dataset monitoring the evolution of a known mineral target in situ on the Martian surface, brought there as a SuperCam calibration target onboard the Perseverance rover. We used Raman spectroscopy to monitor the crystalline state of a synthetic apatite sample over the first 950 Martian days (sols) of the Mars2020 mission. We note significant variations in the Raman spectra acquired on this target, specifically a decrease in the relative contribution of the Raman signal to the total signal. These observations are consistent with the results of a UV-irradiation test performed in the laboratory under conditions mimicking ambient Martian conditions. We conclude that the observed evolution reflects an alteration of the material, specifically the creation of electronic defects, due to its exposure to the Martian environment and, in particular, UV irradiation. This ongoing process of alteration of the Martian surface needs to be taken into account for mineralogical space mission data analysis.

Science & Technology - Other Topics↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reduced Erosion Augments Soil Carbon Storage Under Cover Crops

ABSTRACT Cover crops, a promising strategy to increase soil organic carbon (SOC) storage in croplands and mitigate climate change, have typically been shown to benefit soil carbon (C) storage from increased plant C inputs. However, input‐driven C benefits may be augmented by the reduction of C outputs induced by cover crops, a process that has been tested by individual studies but has not yet been synthesized. Here we quantified the impact of cover crops on organic C loss via soil erosion (SOC erosion) and revealed the geographical variability at the global scale. We analyzed the field data from 152 paired control and cover crop treatments from 57 published studies worldwide using meta‐analysis and machine learning. The meta‐analysis results showed that cover crops widely reduced SOC erosion by an average of 68% on an annual basis, while they increased SOC stock by 14% (0–15 cm). The absolute SOC erosion reduction ranged from 0 to 18.0 Mg C −1 ha −1 year −1 and showed no correlation with the SOC stock change that varied from −8.07 to 22.6 Mg C −1 ha −1 year −1 at 0–15 cm depth, indicating the latter more likely related to plant C inputs. The magnitude of SOC erosion reduction was dominantly determined by topographic slope. The global map generated by machine learning showed the relative effectiveness of SOC erosion reduction mainly occurred in temperate regions, including central Europe, central‐east China, and Southern South America. Our results highlight that cover crop‐induced erosion reduction can augment SOC stock to provide additive C benefits, especially in sloping and temperate croplands, for mitigating climate change.

Huang, Wenjuan [Department of Ecology, Evolution, ↗
DOE OSTI2025

Modeling of the high column density systems in the Lyman-Alpha forest

The Lyman-α forests observed in the spectra of high-redshift quasars can be used as a tracer of the cosmological matter density to study baryon acoustic oscillations (BAO) and the Alcock-Paczynski effect. Extraction of cosmological information from these studies requires modeling of the forest correlations. While the models depend most importantly on the bias parameters of the intergalactic medium (IGM), they also depend on the numbers and characteristics of high-column-density systems (HCDs) ranging from Lyman-limit systems with column densities log N HI /1cm -2 > 17 to damped Lyman-α systems (DLAs) with log N HI /1cm -2 > 20.2. These HCDs introduce broad damped absorption characteristic of a Voigt profile. Consequently they imprint a component on the power spectrum whose modes in the radial direction are suppressed, leading to a scale-dependent bias. Using mock data sets of known HCD content, we test a model that describes this effect in terms of the distribution of column densities of HCDs, the Fourier transforms of their Voigt profiles and the bias of the halos containing the HCDs. Our results show that this physically well-motivated model describes the effects of HCDs with an accuracy comparable to that of the ad-hoc models used in published forest analyses. We also discuss the problems of applying the model to real data, where the HCD content and their bias is uncertain.

Lyman alpha forest↗
DOE OSTI2025

Anisotropic electron damping and energy gap in Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝛿

The many-body electron-electron interaction in cuprates causes broadening of the electronic bands in 𝒌 space, leading to a deviation from the standard Fermi liquid. While a 𝒌-dependent anisotropic electronic scattering (𝒌-DAES) has been assessed by photoemission, its fingerprint in 𝑸 space has been scarcely considered. Here, we explore the 𝑸-dependent electron dynamics in optimally doped Bi 2⁢ Sr 2 ⁢CaCu 2 ⁢O 8+𝛿 through the evolution of low-energy charge excitations as measured by resonant inelastic x-ray scattering (RIXS). In the normal state, the RIXS spectra display a continuum of excitations down to 0 meV, while the superconducting state features a spectral weight suppression below 80 meV without any enhancement at higher energies. To interpret the energy and 𝑸 evolution of our data, we introduce a phenomenological expression of the charge susceptibility by including the 𝒌-DAES. We show that only the charge susceptibility with 𝒌-DAES captures the RIXS data, highlighting the importance of 𝒌-DAES when describing the 𝑸 dependence of charge excitations from 0 to a few eV scale. Furthermore, we also find that the inclusion of 𝒌-DAES is essential when quantitative parameters such as the electronic energy gap are extracted from RIXS data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

State of Innovation 2025: Progress in Accelerating Next-Generation Cement and Concrete Technologies

The cement and concrete sectors are entering a decisive period as next-generation technologies advance from laboratory research to demonstration, early deployment, and first-of-a-kind commercial plants. Building on the 2024 State of Innovation report, the 2025 outlook highlights both the rapid acceleration of innovation and the urgent need for coordinated action across the value chain. Venture capital activity into the cement and concrete space stabilized following the record surge of 2022-2023, yet landmark financings, such as Sublime Systems' $200 million round and Terra CO2's $124 million Series B, signal continued investor confidence in companies approaching commercialization. Corporate procurement has become a powerful new catalyst, with Microsoft, Amazon, and CRH (Cement Roadstone Holdings) Ventures providing long-term commitments that underpin the first wave of next-generation cementitious products. The sector is shifting from early-stage experimentation toward the scaling of well-capitalized leaders capable of bridging the critical "capitalization gap." Early innovators continue to expand the toolkit through novel binders, electrochemical cements, biogenic limestone, and carbonate mineralization pathways. Going into 2026, cost competitiveness, durability validation, and scalability enabled by resilient supply chains remain the decisive factors for market adoption. At the 2025 Next Generation Cement and Concrete Critical Technologies Meeting, attendees emphasized dual-track funding strategies that integrate federal grants with private capital as key to enabling market breakthrough. State programs and corporate demand are sustaining momentum, while successful companies increasingly demonstrate both economic value and reduced dependence on imported materials. The National Concrete Pavement Technology Center and others underscored that broad integration of next-generation materials will hinge on standards compatibility, verified field performance, and workforce readiness. Colorado continues to serve as a proving ground through pilot programs that combine supplier training, phased implementation, and real-world data to de-risk innovation and provide replicable models for other regions. The 2025 Cement and Concrete Critical Technologies Workshop reinforced that scaling next-generation materials will require alignment among technology innovation, performance validation, and market demand. Stakeholders must move beyond siloed efforts toward collaborative frameworks that coordinate standards, funding, and infrastructure deployment. As a neutral convener and technical validator, the National Laboratory of the Rockies (NLR) plays a pivotal role in bridging innovation and market adoption through collaborative research, technology validation, and entrepreneurship programs. By uniting innovators, incumbents, policymakers, contractors, and investors, NLR and its partners are helping chart a credible pathway toward widespread commercialization in the decade ahead.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The ambiguous AT2022rze: changing-look AGN mimicking a supernova in a merging galaxy system

AT2022rze is a luminous, ambiguous transient located south-east of the geometric centre of its host galaxy at redshift $z = 0.08$. The host appears to be formed by a merging galaxy system. The observed characteristics of AT2022rze are reminiscent of active galactic nuclei (AGNs), tidal disruption events, and superluminous supernovae. The transient reached a peak absolute magnitude of $-$20.2 $\pm$ 0.2 mag, showing a sharp rise (t$_{\mathrm{rise,1/e}} = 27.5 \pm 0.6$ d) followed by a slow decline (t$_{\mathrm{dec,1/e}} = 382.9 \pm 0.6$). Its bumpy light curve and narrow Balmer lines indicate the presence of gas (and dust). Its light curve shows rather red colours, indicating that the transient could be affected by significant host extinction. The spectra reveal coronal lines, indicative of high-energy (X-ray/UV) emission. Archival data reveal no prior activity at this location, disfavouring a steady-state AGN, although an optical spectrum obtained prior to the transient is consistent with an AGN classification of the host. Based on this, we conclude that the transient most likely represents a changing-look AGN at the centre of the smallest component of the merging system.

galaxies: active↗
DOE OSTI2025

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A real-time search for Type Ia Supernovae with late-time interactions with circumstellar material in ZTF data

While it is generally accepted that Type Ia supernovae (SNe Ia) are the terminal explosions of white dwarfs (WDs), the nature of their progenitor systems and the mechanisms that lead up to these explosions remain widely debated. In rare cases, the SN ejecta interact with circumstellar material (CSM) that had previously been ejected from the progenitor system. The longer the delay between the creation of the CSM and the SN explosion, the greater the distance between the SN explosion site and the CSM and the later the onset of the interaction. The unknown distance between the CSM and SN explosion site makes it impossible to predict when the interaction will start. If the time between the SN explosion and the onset of the CSM interaction is of the order of several months to years, the SN has generally faded and it is no longer actively followed up on. This makes it even more difficult to detect the interaction while it is happening. In this work, we report on a real-time monitoring programme running between 13 November 2023 and 9 July 2024. It monitored 6914 SNe Ia for signs of late-time rebrightening using the Zwicky Transient Facility (ZTF). Flagged candidates were rapidly followed up on with photometry and spectroscopy to confirm the late-time excess and its position. We report the discovery of a ∼50 day rebrightening event in SN 2020qxz around 1200 rest-frame days after the peak of its light curve. SN 2020qxz exhibited signs of an early CSM interaction, but had faded from view over two years before its reappearance. Initial follow-up spectroscopy revealed the presence of four emission lines, while later follow-up spectroscopy showed that these had faded shortly after the end of the ZTF-detected rebrightening event. Our best match for these emission lines are H β (blueshifted by ∼5900 km s −1 ) and Ca II λ8542 , N I λ8567 , and K I λλ8763, 8767 (all blueshifted by 5100 km s −1 ; although we note that the line identifications are uncertain). This shows that catching and following up on late-time interactions as they occur can offer new clues on the nature of the progenitor systems that produce these SNe by putting constraints on the possible type of donor star. The only way to do this systematically is to use large sky surveys such as ZTF and the upcoming Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) to monitor a large sample of objects for the rare events that reappear long after the object has faded from view.

circumstellar matter↗
DOE OSTI2025

DESI Mg II Absorbers: Extinction Characteristics and Quasar Redshift Accuracy

In this paper, we study how absorption-line systems affect the spectra and redshifts of quasi-stellar objects (QSOs), using catalogs of Mg II absorbers from the early data release and first data release of the Dark Energy Spectroscopic Instrument. We determine the reddening effect of an absorption system by fitting an unreddened template spectrum to a sample of 50,674 QSO spectra that contain Mg II absorbers. We find that reddening caused by intervening absorbers (v off > 3500 km s −1 ) has an average color excess of $\overline{E(B-V)}$ = 0.04 mag. We find that the E(B − V) tends to be greater for absorbers at low redshifts, or those having Mg II absorption lines with higher equivalent widths, but shows no clear trend with v off for intervening systems. However, the $\overline{E(B-V)}$ of associated absorbers, those at v off < 3500 km s −1 , shows a strong trend with v off , increasing rapidly with decreasing v off and peaking (∼0.15 mag) around v off = 0 km s −1 . We demonstrate that Mg II absorbers impact redshift estimation for QSOs by investigating the distributions of v off for associated absorbers. We find that at z > 1.5, these distributions broaden and bifurcate in a nonphysical manner. In an effort to mitigate this effect, we mask pixels associated with the Mg II absorption lines and recalculate the QSO redshifts. We find that we can recover voff populations in better agreement with those for z < 1.5 absorbers and in doing so typically shift background QSO redshifts by Δz ≈ ± 0.005.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

JWST’s PEARLS: A z ≃ 6 quasar in a train-wreck galaxy merger system

We present JWST NIRSpec integral field spectroscopy observations of the z = 5.89 quasar NDWFS J1425+3254 from 0.6–5.3 μm, covering the rest-frame ultraviolet and optical at a spectral resolution of R ∼ 100. The quasar has a black hole mass of M BH = (1.4 +3.1 −1.0 ) × 10 9 M ⊙ and an Eddington ratio of L Bol /L Edd = 0.3 +0.6 −0.2 , as implied from the broad Balmer Hα and Hβ lines. The quasar host has significant ongoing obscured star formation, as well as a quasar-driven outflow with velocity 6050 +460 −630 km s −1 and ionised outflow rate of 1650 +130 −1230 M ⊙ yr −1 . This is possibly one of the most extreme outflows in the early Universe. The data also reveal that two companion galaxies are merging with the quasar host. The north-eastern companion galaxy is relatively old and very massive, with a luminosity-weighted stellar age of 65 +9 −4 Myr, stellar mass of (3.6 +0.6 −0.3 #x00D7; 10 11 M ⊙ , and star-formation rate (SFR) of ∼15–30 M ⊙ yr −1 . A bridge of gas connects this companion galaxy and the host, confirming their ongoing interaction. A second merger is occurring between the quasar host and a much younger companion galaxy to the south, with a stellar age of 6.7 ± 1.8 Myr, stellar mass of (1.9 ± 0.4)×10 10 M ⊙ , and SFR of ∼40–65 M ⊙ yr −1 . There is also another galaxy in the field, likely in the foreground at z = 1.135, which could be gravitationally lensing the quasar with a magnification of 1 < μ < 2 and, thus, < 0.75 mag. Overall, the system is a ‘train-wreck’ merger of three galaxies, with star formation and extreme quasar activity that were likely triggered by these ongoing interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Estimating the Electric Fields Driving Lightning Dart Leader Development With BIMAP‐3D Observations

Abstract In this paper, a numerical dart leader model has been implemented to understand the leader's development and the corresponding electric field changes observed by the 3D broadband mapping and polarization (BIMAP‐3D) system. The model assumes the extending leader channel is equipotential and has a linear charge distribution induced by an ambient electric field. The charge distribution induced by the ambient field can be used to model the electric field change at the ground. We then find the ambient electric field which best fits the field change measurements at the two BIMAP stations. The estimated ambient electric field decreases in the direction of dart leader propagation. Our observations and modeling results are consistent with our earlier hypothesis that dart leader speed is proportional to the electric field at the leader tip. The model also supports our earlier analysis that leader speed variations near branch junctions were due to previous charge deposits near the junctions. The modeled tip electric field is generally lower than the breakdown field unless the pre‐dart‐leader channel has a significant temperature of 3,000 K. This is consistent with the fact that dart leaders typically do not form new branches into the virgin air. Furthermore, the tip field is generally close to the negative streamer stability field at ambient temperatures, explaining the nature of the narrow and well‐defined channel structure. In addition to the charge distribution and the ambient and tip electric field, the development of the channel potential and current distribution are also presented.

3D lightning mapping↗
DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗
DOE OSTI2026

Measurements of Quasar Proximity Zones with the Lyα Forest of DESI Y1 Quasars

Abstract The intergalactic medium (IGM) around a quasar is shaped by its dense environment and by its excess ionizing radiation, which form a “quasar proximity zone” whose size and anisotropy depend on the quasar’s halo mass, luminosity, age, and radiation geometry. Using over 10,000 quasar pairs from the Dark Energy Spectroscopic Instrument (DESI) Year 1 data, with projected comoving separations r ⊥ < 2 h −1 Mpc, we investigate how the proximity zone of foreground quasars at z ∼ 2–3.5 affects Lyα absorption in their background quasars. The large DESI sample enables unprecedented precision in measuring this “transverse proximity” effect, allowing a detailed investigation of the signal’s dependence on the projected separation of quasar pairs and the luminosity of the foreground quasar. We find that enhanced gas clustering near quasars dominates over their ionizing effect, leading to stronger absorption on neighboring sightlines. Under the assumption that quasar ionizing luminosity is isotropic and steady, we infer the IGM overdensity profile in the vicinity of quasars, finding overdensities as high as Δ ∼ 10 at comoving distance ∼1 h −1 Mpc from the most luminous systems. Surprisingly, however, we find no significant dependence of the proximity profile on the luminosity of the foreground quasar. This lack of luminosity dependence could reflect a cancellation between higher ionizing flux and higher gas overdensity, or it could indicate that quasar emission is highly time-variable or anisotropic, so that the observed luminosity does not trace the ionizing flux on nearby sightlines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Multi‐Decadal Dynamics of Wetland Methane Emissions Revealed by Knowledge‐Guided Machine Learning

Measurement of methane fluxes (FCH 4 ) from natural systems, such as wetlands, has lagged far behind carbon dioxide fluxes. Short and fragmented wetland FCH 4 data limit our ability to assess its long-term dynamics and potential climate feedbacks. Extrapolating short-term FCH 4 records to recent decades remains challenging for both process-based models and data-driven machine learning (ML) approaches. Here, we develop a knowledge-guided ML framework that integrates eddy covariance (EC) FCH 4 observations, field warming experiments, and biogeochemical knowledge to reconstruct the long-term FCH 4 budgets and trends. Focusing on the 11 longest EC monitoring sites in the AmeriFlux network, we found considerable variability in multi-decadal trends of wetland FCH 4 , with increases up to 14% per decade from 2000 to 2024. We also found that the strength of these increasing trends declines from high to low latitudes, highlighting the vulnerability of northern wetlands. This work presents novel and robust reconstructions of long-term wetland FCH 4 , offering critical benchmark datasets for bottom-up ecosystem models and advancing fundamental understanding of wetland biogeochemistry.

AmeriFlux site↗
DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗
DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗