SEARCH · Search NASA

Results for “Background Estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

CalderaCast User Manual Version 1.0

CalderaCast is a user-friendly web-based tool for electrical-load forecast, providing stakeholders with a fully customizable decision-support framework that estimates the likely power draw from a possible future electric-vehicle (EV) charging station at a given location on a given day along an alternative fuel corridor (AFC). These EV charging profiles are accurately modeled in CalderaCast using the Caldera software framework developed by Idaho National Laboratories (INL), reflecting the realistic charging levels observed in actual charge events. This tool was developed as part of the National Electric Vehicle Infrastructure (NEVI) program, which is quickly generating substantial interest from would-be charging station operators (CSO), large and small electric utilities, and state transportation planners, some of whom had not seriously considered EV charging previously. All these entities—with or without background in EV infrastructure—must estimate the electricity load that a proposed charging station will generate. This load forecast is critically important for a utility to properly assess the capacity of their distribution network to support the proposed station or properly size grid upgrades for potential load growth due to future EV adoption, vehicle technology improvements, or station growth. This document describes each aspect of the CalderaCast tool and provides guidance to users who are interested in utilizing the tool for their work.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Deriving Stable Peak Models to Fit Complex XPS Data From Cu Contaminated Pt Electrocatalysts

X-ray Photoelectron Spectroscopy spectra peak models, designed to partition photoemission signals emanating from different elements or chemical states within an atom, are fitted to data limited to an energy interval over which inelastically scattered photoemission signal can be estimated. While the choice of background approximation and line shapes of components to the peak model requires careful consideration, the energy interval used to define the data to which the peak model is optimized has a significant impact on the final peak model. The relationship between the background intensity and data intensity at the start and end of the energy interval dictates the line shapes used in the peak model. In this work, we devise a method to peak fit a complex overlapping Cu 3p and Pt 4f XPS peak structure to perform the elemental quantification. We first use an Al 2s peak to illustrate how background curves approach data at the limits of the energy interval over which the background is defined, influencing the analysis of XPS spectra. Next, we demonstrate the nature of interactions between specific line shapes (Voigt and pseudo-Voigt profiles) suitable for photoemission peaks and a specific background curve (Shirley) and a peak model is presented that includes components to the peak model that accommodates background intensity during fitting of the peak model to data. The peak model allowed for quantification of the contributions of Pt 4f peaks emanating from the substrate that exhibits strong asymmetry in the presence of the inhomogeneously distributed Cu species, mostly of Lorentzian character.

XPS↗

DOE OSTI2025

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗

DOE OSTI2025

Evaluation of a generalized least squares algorithm for infrasound beamforming with coherent background noise

Infrasonic signals of interest can occur during periods with persistent, coherent, background noise, which may be natural or anthropogenic. For high signal-to-noise (SNR) ratio transient signals, an ‘overprinting’ of the coherent background may occur, and the signal may still be detected. However, this approach fails for low SNR signals of interest, which may be obscured by coherent noise. An infrasound beamforming method based on generalized least squares (GLS) is investigated for detecting transient signals of interest in the presence of coherent and incoherent background noise. This approach relies on an estimate of the noise covariance, captured in a covariance matrix, to effectively null contributions to the array response from noisy directions of arrival. Synthetic array data is used to investigate the performance of the GLS beamformer compared to the Bartlett beamformer when coherent and incoherent backgrounds are present. Additionally, the effects of array element number and relative strength of the interfering signal on the GLS estimates is investigated. GLS empirical area under the curve estimates suggest that the beamformer can recover coherent power for a signal of interest lower in amplitude than the coherent background, but this effectiveness degrades more quickly with SNR for a four element array compared to a six or eight element infrasound array. Finally, infrasound from the Forensic Surface Experiment, a bolide signal observed at IMS array I37NO, and a volcanic signal recorded at the Alaska Volcano Observatory array ADKI are used to evaluate GLS performance on recorded data. A ten minute window was used to capture the background noise, and the coherent background signal was nulled in all three examples.

58 GEOSCIENCES↗

DOE OSTI2024

First Constraints on the Epoch of Reionization Using the Non-Gaussianity of the Kinematic Sunyaev-Zel’dovich Effect from the South Pole Telescope and Herschel-SPIRE Observations

We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel’dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and Herschel-SPIRE experiments over a 100 deg 2 field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 μ K - arcmin in bands centered at 95, 150, and 220 GHz. For SPIRE, we include data from the 600 and 857 GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at 10.3 σ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the agora simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this nondetection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of Δ z re , 50 < 4.5 (95% confidence level). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, Planck’s optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Atmospheric pion, kaon, and muon fluxes for sub-orbital experiments

Cosmic rays interacting with the Earth's atmosphere generate extensive air showers, which produce Cherenkov, fluorescence and radio emissions. These emissions are key signatures for detection by ground-based, sub-orbital, and satellite-based telescopes aiming to study high energy cosmic ray and neutrino events. However, detectors operating at ground and balloon altitudes are also exposed to a background of atmospheric charged particles, primarily pions, kaons, and muons, that can mimic or obscure the signals from astrophysical sources. In this work, we use coupled cascade equations to calculate the atmospheric pion, kaon and muon fluxes reaching detectors at various altitudes. Our analysis focuses on energies above 10 GeV, where the influence of the Earth's magnetic field on particle trajectories is minimal. We provide angular and energy-resolved flux estimates and discuss their relevance as background for extensive air shower detection. Furthermore, our results are potentially relevant for interpreting data from current and future balloon-borne experiments such as EUSO-SPB2 and for refining trigger and veto strategies in Cherenkov and fluorescence telescopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

DOE OSTI2025

Monopole Fluctuations in Galaxy Surveys

Galaxy clustering provides a powerful way to probe cosmology. This requires understanding of the background mean density of galaxy samples, which is estimated from the survey itself by averaging the observed galaxy number density over the angular position. The angle average includes not only the background mean density but also the monopole fluctuation at each redshift. Here for the first time we compute the monopole fluctuations in galaxy surveys and investigate their impact on galaxy clustering. The monopole fluctuations vary as a function of redshift, and it is correlated with other fluctuations, affecting the two-point correlation function measurements. In an idealized all-sky survey, the rms fluctuation at z = 0.5 can be as large as 7% of the two-point correlation function in amplitude at the baryonic acoustic oscillation scale, and it becomes smaller than 1% at z > 2. The monopole fluctuations are unavoidable, but they can be modeled. We discuss its relation to the integral constraint and the implications for the galaxy clustering analysis.

Baryon acoustic oscillations↗

DOE OSTI2024

Robust anomaly detection for particle physics using multi-background representation learning

Abstract Anomaly, or out-of-distribution, detection is a promising tool for aiding discoveries of new particles or processes in particle physics. In this work, we identify and address two overlooked opportunities to improve anomaly detection (AD) for high-energy physics. First, rather than train a generative model on the single most dominant background process, we build detection algorithms using representation learning from multiple background types, thus taking advantage of more information to improve estimation of what is relevant for detection. Second, we generalize decorrelation to the multi-background setting, thus directly enforcing a more complete definition of robustness for AD. We demonstrate the benefit of the proposed robust multi-background AD algorithms on a high-dimensional dataset of particle decays at the Large Hadron Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Dielectron production in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02TeV

The first measurement of the 𝑒 + ⁢𝑒 − pair production at midrapidity and low invariant mass in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider is presented. The yield of 𝑒 + ⁢𝑒 − pairs is compared with a cocktail of expected hadronic decay contributions in the invariant mass (𝑚𝑒⁢𝑒) and pair transverse momentum (𝑝 T,𝑒⁢𝑒 ) ranges 𝑚 𝑒⁢𝑒 < 3.5 GeV/𝑐 2 and 𝑝 T,𝑒⁢𝑒 < 8 GeV/𝑐. For 0.18 < 𝑚 𝑒⁢𝑒 < 0.5 GeV/𝑐 2 the ratio of data to the cocktail of hadronic contributions amounts to 1.40 ± 0.11⁢(stat.) ± 0.23 ⁢(syst.) ± 0.16⁢ (cocktail) and 1.42 ± 0.11⁢ (stat.) ± 0.23 (syst.)$^{+0.24}_{−0.29}$ ⁢(cocktail), including or not including medium effects in the estimation of the heavy-flavor background, respectively. It is consistent with predictions from two different models for an additional contribution of thermal 𝑒+⁢𝑒− pairs from the hadronic and partonic phases. In the intermediate-mass range (1.2 < 𝑚 𝑒⁢𝑒 < 2.6 GeV/𝑐 2 ), the pair transverse impact parameter of the 𝑒 + ⁢𝑒 − pairs (DCA𝑒⁢𝑒, where “DCA” denotes “distance of closest approach”) is used for the first time in Pb-Pb collisions to separate displaced dielectrons from heavy-flavor hadron decays from a possible (thermal) contribution produced at the interaction point. The data are consistent with a suppression of 𝑒 + ⁢𝑒 − pairs from $𝑐\bar{𝑐}$ and an additional prompt component. Finally, the first direct-photon measurement in the 10% most central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV is reported via the study of virtual direct photons in the transverse momentum range 1 < 𝑝 T < 5 GeV/𝑐. A model including prompt photons, as well as photons from the preequilibrium and fluid-dynamic phases, can reproduce the result, while being at the upper edge of the data uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

An iterative CMB lensing estimator minimizing instrumental noise bias

Noise maps from cosmic microwave background (CMB) experiments are generally statistically anisotropic, due to scanning strategies, atmospheric conditions, or instrumental effects. Any mismodeling of this complex noise can bias the reconstruction of the lensing potential and the measurement of the lensing power spectrum from the observed CMB maps. We introduce a new CMB lensing estimator based on the maximum (MAP) reconstruction that is minimally sensitive to these instrumental noise biases. By modifying the likelihood to rely exclusively on correlations between CMB map splits with independent noise realizations, we minimize autocorrelations that contribute to biases. In the regime of many independent splits, this maximum closely approximates the optimal MAP reconstruction of the lensing potential. In simulations, we demonstrate that this method is able to determine lensing observables that are immune to any noise mismodeling with a negligible cost in signal-to-noise ratio. Our estimator enables unbiased and nearly optimal lensing reconstruction for next-generation CMB surveys.

Legrand, Louis [University of Cambridge (United Ki↗

DOE OSTI2024

Upper limit on the chiral magnetic effect in isobar collisions at the Relativistic Heavy-Ion Collider

The chiral magnetic effect (CME) is a phenomenon that arises from the QCD anomaly in the presence of an external magnetic field. The experimental search for its evidence has been one of the key goals of the physics program of the Relativistic Heavy-Ion Collider. The STAR Collaboration has previously presented the results of a blind analysis of isobar collisions ( Ru 44 96 + Ru 44 96 , Zr 40 96 + Zr 40 96 ) in the search for the CME. The isobar ratio ( Y ) of CME-sensitive observable, charge separation scaled by elliptic anisotropy, is close to but systematically larger than the inverse multiplicity ratio, the naive background baseline. This indicates the potential existence of a CME signal and the presence of remaining nonflow background due to two- and three-particle correlations, which are different between the isobars. In this postblind analysis, we estimate the contributions from those nonflow correlations as a background baseline to Y , utilizing the isobar data as well as Heavy Ion Jet Interaction Generator simulations. This baseline is found consistent with the isobar ratio measurement, and an upper limit of 10% at 95% confidence level is extracted for the CME fraction in the charge separation measurement in isobar collisions at s NN = 200 GeV. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Searching for heavy neutral leptons at a future muon collider

As the planning stages for a high energy muon collider enter a more concrete era, an important question arises as to what new physics could be uncovered. A TeV-scale muon collider is also a vector boson fusion (VBF) factory with a very clean background, and as such it is a promising environment to look for new physics that couples to the electroweak (EW) sector. In this paper, we explore the ability of a future TeV-scale muon collider to search for Majorana and Dirac heavy neutral leptons (HNLs) produced via EW bosons. Employing a model-independent, conservative approach, we present an estimation of the production and decay rate of HNLs over a mass range between 200 GeV and 9.5 TeV in two benchmark collider proposals with s = 3 , 10 TeV, as well as an estimation of the dominant Standard Model (SM) background. We find that exclusion limits for the mixing between the HNLs and SM neutrinos can be as low as O ( 10 − 6 ) . Additionally, we demonstrate that a TeV-scale muon collider allows for the ability to discriminate between Majorana and Dirac type HNLs for a large range of mixing values. Published by the American Physical Society 2024

Kwok, Tsz Hong↗

DOE OSTI2025

Bias-hardened estimators of patchy screening profiles

Detecting anisotropic screening of the cosmic microwave background (CMB) holds the promise of revealing the distribution of gas in the Universe, characterizing the complex processes of galaxy formation and feedback, and studying the epoch of reionization. Estimators for inhomogeneous screening, including some recently proposed small-scale (stacked) estimators, are quadratic or higher order in the CMB temperature or polarization fields and are therefore subject to contamination from CMB lensing. We review the origin of this lensing bias and show that, when stacking on unWISE galaxies, the expected lensing bias dominates the signal if left unmitigated. Hardening techniques that null the lensing bias have been proposed for standard quadratic estimators, whereas only approximate methods have been proposed for stacked estimators. In conclusion, we review these techniques and apply the former to stacked estimators, presenting several strategies (including the optimal strategy) to null lensing contamination when stacking on any large-scale structure tracer.

Cosmic microwave background↗

DOE OSTI2025

Light in the dark forest. Part I. An efficient optimal estimator for 3D Lyman-alpha forest power spectrum

The highly anisotropic nature of the Lyman-alpha (Lyα) forest data introduces a complex survey window function that complicates the measurement of the three-dimensional power spectrum ( P 3D ). In this paper, we present the first fully optimal estimator for P 3D , which exactly deconvolves the survey window function and marginalizes contaminated modes that distort the power spectrum. Our approach adapts optimal estimator techniques developed for the 2D cosmic microwave background data to the 3D case. To achieve computational feasibility, we employ the conjugate gradient method and implement the P 3 M formalism to handle large-scale and small-scale operations separately and efficiently. We validate our estimator using Monte Carlo mocks and Gaussian simulations, demonstrating its accuracy and computational efficiency. We confirm that mode marginalization eliminates distortions arising from quasar continuum errors and delivers robust power spectrum estimation, though it also inflates errors at large scales. This first implementation works in the flat-sky case; we discuss the remaining steps needed to generalize it to the curved-sky case. This formalism offers a foundation for the Lyα forest P 3D measurements and a new path toward cosmological constraints from the Lyα forest data.

Lyman alpha forest↗

DOE OSTI2024

Search for light long-lived neutral particles from Higgs boson decays via vector-boson-fusion production from pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search is reported for long-lived dark photons with masses between 0.1 GeV and 15 GeV, from exotic decays of Higgs bosons produced via vector-boson-fusion. Events that contain displaced collimated Standard Model fermions reconstructed in the calorimeter or muon spectrometer are probed. This search uses the full LHC Run 2 (2015–2018) data sample collected in proton–proton collisions at $\sqrt{s}=13$ TeV, corresponding to an integrated luminosity of 139 fb –1 . Dominant backgrounds from Standard Model processes and non-collision sources are estimated using data-driven techniques. The observed event yields in the signal regions are consistent with the expected background. Upper limits on the Higgs boson to dark photon branching fraction are reported as a function of the dark photon mean proper decay length or of the dark photon mass and the coupling between the Standard Model and the potential dark sector. This search is combined with previous ATLAS searches obtained in the gluon–gluon fusion and WH production modes. A branching fraction above 10% is excluded at 95% CL for a 125 GeV Higgs boson decaying into two dark photons for dark photon mean proper decay lengths between 173 and 1296 mm and mass of 10 GeV.

Aad, G. (ORCID:0000000266654934)↗

DOE OSTI2025

Search for single-production of vector-like quarks decaying into Wb in the fully hadronic final state in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for T and Y vector-like quarks produced in proton-proton collisions at a centre-of-mass energy of 13 TeV and decaying into Wb in the fully hadronic final state is presented. The search uses 139 fb −1 of data collected by the ATLAS detector at the LHC from 2015 to 2018. The final state is characterised by a hadronically decaying W boson with large Lorentz boost and a b-tagged jet, which are used to reconstruct the invariant mass of the vector-like quark candidate. The main background is QCD multijet production, which is estimated using a data-driven method. Upon finding no significant excess in data, mass limits at 95% confidence level are obtained as a function of the global coupling parameter, κ. The observed lower limits on the masses of Y quarks with κ = 0.5 and κ = 0.7 are 2.0 TeV and 2.4 TeV, respectively. For T quarks, the observed mass limits are 1.4 TeV for κ = 0.5 and 1.9 TeV for κ = 0.7.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Which Plant Traits Increase Soil Carbon Sequestration? Empirical Evidence From a Long‐Term Poplar Genetic Diversity Trial

Plants play a key role in mediating soil response to global change, and breeding or engineering crops to increase soil organic carbon (SOC) storage is a potential route to land-based carbon dioxide removal in agricultural systems. However, due to limited observational datasets plus shifting paradigms of SOC stabilization, it is unclear which plant traits are most important for enhancing different types of soil organic matter. Existing long-term common gardens of genetically diverse plant populations may provide an opportunity to evaluate biological controls on SOC, separate from environmental or management variability. Here we report on soil and root chemical data collected for 24 genotypes within a 13-year-old common garden in northwestern Oregon planted with a large natural variant population of Populus trichocarpa. Fractionating surface soil (0–15 cm) revealed substantial variation in stocks of mineral-associated organic matter (MAOM; 18–67 t C/ha) and particulate organic matter (POM; 2–22 t C/ha). Tree genotype explained 24% and 26% of the MAOM and POM stock variability, respectively, after controlling for background variability. We found minimal association between SOC concentration and either aboveground tree productivity or root biomass recalcitrance (C/N ratios and lignin content). In contrast, root elemental content appeared influential for MAOM-C concentration, which showed a strong positive association with root aluminum (Al) and a strong negative association with root boron (B) and magnesium (Mg). Furthermore, root concentrations of these elements were highly heritable (57%–78%) and not simply a reflection of background variation in soil elemental concentrations. We estimate that surface SOC stocks under these 24 genotypes have diverged at rates of up to 1.2–4.3 t C/ha/year. These results suggest that long-term genetic diversity trials have value for elucidating biological controls on soil organic matter dynamics, and that traits associated with root elemental content may be a useful target for enhancing biosequestration.

biomass recalcitrance↗