SEARCH · Engineering Papers

Results for “Reanalysis”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Fill and dump measurement of the neutron lifetime using an asymmetric magneto-gravitational trap

The past two decades have yielded several new measurements and reanalysis of older measurements of the neutron lifetime. These have led to a 4.4 standard deviation discrepancy between the most precise measurements of the neutron decay rate producing protons in cold neutron beams and the most precise lifetime measured in neutron storage experiments. Measurements using different techniques are important for investigating whether there are unidentified systematic effects in any of the measurements. In this paper we report a new measurement using the Los Alamos asymmetric magneto-gravitational trap where the surviving neutrons are counted external to the trap using the fill and dump method. In this work, the new measurement gives a free neutron lifetime of $τ_n$ = 876.3(2.4) stat (0.8) syst . Although this measurement is not as precise, it is in statistical agreement with previous results using in situ counting in the same apparatus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Polarization transfer observables in elastic electron-proton scattering at Q 2 = 2.5 , 5.2, 6.8, and 8.5 GeV 2

In this paper, interest in the behavior of nucleon electromagnetic form factors at large momentum transfers has steadily increased since the discovery, using polarization observables, of the rapid decrease of the ratio G p E /G p M of the proton's electric and magnetic form factors for momentum transfers Q 2 ≳ 1 GeV 2 , in strong disagreement with previous extractions of this ratio using the traditional Rosenbluth separation technique. The GEp-III and GEp-2$\gamma$ experiments were carried out in Jefferson Lab's (JLab's) Hall C from 2007-2008, to extend the knowledge of $G_E^p/G_M^p$ to the highest practically achievable $Q^2$ and to search for effects beyond the Born approximation in polarization transfer observables of elastic $\vec{e}p$ scattering. This article reports an expanded description of the common experimental apparatus and data analysis procedure, and the results of a final reanalysis of the data from both experiments, including the previously unpublished results of the full-acceptance data of the GEp-2$\gamma$ experiment. The Hall C High Momentum Spectrometer detected and measured the polarization of protons recoiling elastically from collisions of JLab's polarized electron beam with a liquid hydrogen target. A large-acceptance electromagnetic calorimeter detected the elastically scattered electrons in coincidence to suppress inelastic backgrounds. The final GEp-III data are largely unchanged relative to the originally published results. The statistical uncertainties of the final GEp-2$\gamma$ data are significantly reduced at $\epsilon = 0.632$ and $0.783$ relative to the original publication. The decrease with $Q^2$ of $G_E^p/G_M^p$ continues to $Q^2 = 8.5$ GeV$^2$, but at a slowing rate relative to the approximately linear decrease observed in earlier Hall A measurements. At $Q^2 = 2.5$ GeV$^2$, the proton form factor ratio $G_E^p/G_M^p$ shows no statistically significant $\epsilon$-dependence, as expected in the Born approximation. The ratio $P_\ell/P_\ell^{Born}$ of the longitudinal polarization transfer component to its Born value shows an enhancement of roughly 1.7\% at $\epsilon = 0.783$ relative to $\epsilon = 0.149$, with $\approx 2.2\sigma$ significance based on the total uncertainty, implying a similar effect in the transverse component $P_t$ that cancels in the ratio $R$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Searching for states analogous to the C 12 Hoyle state in heavier nuclei using the thick target inverse kinematics technique

Identification of α-cluster states analogous to the C 12 Hoyle state in heavier α-conjugate nuclei can provide tests of the existence of α condensates in nuclei. Such states are predicted for O 16 , Ne 20 , Mg 24 , Si 28 , etc., at excitation energies slightly above the multi- α-particle decay threshold but have not yet been experimentally identified. The thick target inverse kinematics (TTIK) technique can be used to study the breakup of excited self-conjugate nuclei into many α particles. The reaction Ne 20 + α was studied using a Ne 20 beam at 12 MeV/nucleon from the K150 cyclotron at Texas A&M University. The TTIK method was used to study both single α-particle emission and multiple α-particle decays. Events with α multiplicity up to four were analyzed. The analysis of the three α-particle emission data allowed the identification of the Hoyle state and other C 12 excited states decaying into three α particles. The results are shown and compared with other data available in the literature. Although the statistics for events with α-multiplicity four is low, the data show a structure at about 15.2 MeV that could indicate the existence in O 16 of a state analogous to the C 12 Hoyle state. Moreover, the reconstructed excitation energy of Mg 24 for these events peaks at around 34 MeV, very close to the predicted excitation energy for an excited state analogous to the C 12 Hoyle state in Mg 24 . The structure is further confirmed by the reanalysis of α-multiplicity-four events from a previous experiment performed at 9.7 MeV/nucleon with a similar, but lower granularity, experimental setup.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

QFT in curved spacetime from quantum gravity: Proper WKB decomposition of the gravitational component

Starting from a reanalysis of previous work, we construct the proper low-energy quantum field theory (QFT) limit of a full quantum gravity theory in the Born-Oppenheimer approach. We separate the gravitational sector into a classical background, given by a vacuum diagonal Bianchi I cosmology, and its quantum perturbations represented by the two graviton degrees of freedom; we further include quantum matter in the form of a test scalar field. We then implement a Born-Oppenheimer separation, where the gravitons and matter play the roles of “slow” and “fast” quantum components, respectively, and perform a WKB expansion in a Planckian parameter. The functional Schrödinger evolution for matter is recovered after averaging over quantum-gravitational effects, provided that a condition is imposed on the gravitons’ wave functional. Such a condition fixes the graviton dynamics and is equivalent to the purely gravitational Wheeler-DeWitt constraint imposed in previous approaches. Here, the main accomplishment of the present work is to clarify that QFT in curved spacetime can be recovered in the low-energy limit of quantum gravity only after averaging over the graviton degrees of freedom, in the spirit of effective field theory. Furthermore, it justifies a posteriori the implementation of the gravitational Wheeler-DeWitt equation on the “slow” gravitons’ wave functional rather than assuming its validity a priori.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Precise cosmological constraints from BOSS galaxy clustering with a simulation-based emulator of the wavelet scattering transform

For this study, we perform a reanalysis of the BOSS CMASS DR12 galaxy dataset using a simulation-based emulator for the wavelet scattering transform (WST) coefficients. Moving beyond our previous works, which laid the foundation for the first galaxy clustering application of this estimator, we construct a neural net-based emulator for the cosmological dependence of the WST coefficients and the 2-point correlation function multipoles, trained from the state-of-the-art suite of abacussummit simulations combined with a flexible halo occupation distribution (HOD) galaxy model. In order to confirm the accuracy of our pipeline, we subject it to a series of thorough internal and external mock parameter recovery tests, before applying it to reanalyze the CMASS observations in the redshift range 0.46 < z < 0.57. We find that a joint WST+2-point correlation function likelihood analysis allows us to obtain marginalized 1⁢σ errors on the Λ⁢ CDM parameters that are tighter by a factor of 2.5–6, compared to the 2-point correlation function, and by a factor of 1.4–2.5 compared to the WST-only results. This corresponds to a competitive 0.9%, 2.3% and 1% level of determination for parameters ω c , ⁢σ 8 &n s , respectively, and also to a 0.7% and 2.5% constraint on derived parameters h and ƒ⁡(z)⁢⁢σ 8 ⁡(z), in agreement with the Planck 2018 results. Our results reaffirm the constraining power of the WST and highlight the exciting prospect of employing higher-order statistics in order to fully exploit the power of upcoming stage-IV spectroscopic observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

New constraints on ultraheavy dark matter from the LZ experiment

Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV / c 2 to a few TeV / c 2 . Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal freeze-out, but they have generally been less explored experimentally. In this work, we present a reanalysis of the first science run of the LZ experiment, with an exposure of 0.9 tonne × yr , to search for ultraheavy particle dark matter. The signal topology consists of multiple energy deposits in the active region of the detector forming a straight line, from which the velocity of the incoming particle can be reconstructed on an event-by-event basis. Zero events with this topology were observed after applying the data selection calibrated on a simulated sample of signal-like events. New experimental constraints are derived, which rule out previously unexplored regions of the dark matter parameter space of spin-independent interactions beyond a mass of 10 17 GeV / c 2 . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Noble dark matter: Surprising elusiveness of dark baryons

Dark matter could be a baryonic composite of stronglyc upled constituents transforming under SU⁡(2) 𝐿 . We classify the SU⁡(2) 𝐿 representations of baryons in a class of simple confining dark sectors and find that the lightest state can be a pure singlet or a singlet that mixes with other neutral components of SU⁡(2) 𝐿 representations, which strongly suppresses the dark matter candidate’s interactions with the Standard Model. We focus on models with a confining SU⁡(𝑁 𝑐 ) and heavy dark quarks constituting vectorlike 𝑁 𝑓 -plet of SU⁢(2) 𝐿 . For benchmark 𝑁 𝑐 and 𝑁 𝑓 , we calculate baryon mass spectra, incorporating electroweak gauge boson exchange in the nonrelativistic quark model, and demonstrate that above TeV mass scales, dark matter is dominantly a singlet state. The combination of this singlet nature with the recently discovered ℋ-parity results in an inert state analogous to noble gases, hence we coin the term noble dark matter. Our results can be understood in the nonrelativistic effective theory that treats the dark baryons as elementary states, where we find singlets accompanying triplets, 5-plets, or more exotic representations. This generalization of WIMP-like theories is more difficult to find or rule out than dark matter models that include only a single SU⁡(2) 𝐿 multiplet (such as a wino), motivating new searches in colliders and a reanalysis of direct and indirect detection prospects in astrophysical observations.

dark matter direct detection↗

DOE OSTI2025

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Reducing Uncertainties in Large-Scale Solar Resource Data: The Impact of Aerosols

This study examines the impact of various aerosol databases on the accuracy of irradiance predictions pertaining to the National Solar Radiation Database (NSRDB) over the Americas. Compared to the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) reanalysis and moderate resolution imaging spectroradiometer (MODIS) space-borne observations, the various improvements that have been made to develop the SOLSUN aerosol database are described, using a monthly temporal resolution and a 4-km spatial resolution. Relative to MERRA-2, these improvements result in lower uncertainty in the solar irradiance components. The most recent version of the NSRDB now uses hourly MERRA-2 aerosol data with bias and elevation correction. This significantly improves the frequency distribution of hourly direct irradiance while having only a small impact on the annual solar resource, as suggested by the validation undertaken here at 12 control sites.

14 SOLAR ENERGY↗

DOE OSTI2023

TROPHY: A Topologically Robust Physics-Informed Tracking Framework for Tropical Cyclones

Tropical cyclones (TCs) are among the most destructive weather systems. Realistically and efficiently detecting and tracking TCs are critical for assessing their impacts and risks. In particular, the eye is a signature feature of a mature TC. Therefore, knowing the eyes’ locations and movements is crucial for both operational weather forecasts and climate risk assessments. Recently, a multilevel robustness framework has been introduced to study the critical points of time-varying vector fields. The framework quantifies the robustness (i.e., structural stability) of critical points across varying neighborhoods. By relating the multilevel robustness with critical point tracking, the framework has demonstrated its potential in cyclone tracking. An advantage is that it identifies cyclonic features using only 2D wind vector fields, which is encouraging as most tracking algorithms require multiple dynamic and thermodynamic variables at different altitudes. A disadvantage is that the framework does not scale well computationally for datasets containing a large number of cyclones. Herein this paper introduces a topologically robust physics-informed tracking framework (TROPHY) for TC tracking. The main idea is to integrate physical knowledge of TC to drastically improve the computational efficiency of multilevel robustness framework for large-scale climate datasets. First, during preprocessing, we propose a physics-informed feature selection strategy to filter 90% of critical points that are short-lived and have low stability, thus preserving good candidates for TC tracking. Second, during in-processing, we impose constraints during the multilevel robustness computation to focus only on physics-informed neighborhoods of TCs. We apply TROPHY to 30 years of 2D wind fields from reanalysis data in ERA5 and generate a number of TC tracks. In comparison with the observed tracks, we demonstrate that TROPHY can capture TC characteristics (e.g., frequency, intensity, duration, latitudes with maximum intensity, and genesis) that are comparable to and sometimes even better than a well-validated TC tracking algorithm that requires multiple dynamic and thermodynamic scalar fields.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Climate regime shift and forest loss amplify fire in Amazonian forests

Frequent Amazonian fires over the last decade have raised the alarm about the fate of the Earth's most biodiverse forest. The increased fire frequency has been attributed to altered hydrological cycles. However, observations over the past few decades have demonstrated hydrological changes that may have opposing impacts on fire, including higher basin-wide precipitation and increased drought frequency and severity. Here, we use multiple satellite observations and climate reanalysis datasets to demonstrate compelling evidence of increased fire susceptibility in response to climate regime shifts across Amazonia. In this work, we show that accumulated forest loss since 2000 warmed and dried the lower atmosphere, which reduced moisture recycling and resulted in increased drought extent and severity, and subsequent fire. Extremely dry and wet events accompanied with hot days have been more frequent in Amazonia due to climate shift and forest loss. Simultaneously, intensified water vapor transport from the tropical Pacific and Atlantic increased high-altitude atmospheric humidity and heavy rainfall events, but those events did not alleviate severe and long-lasting droughts. Amazonia fire risk is most significant in the southeastern region where tropical savannas undergo long seasonally dry periods. We also find that fires have been expanding through the wet-dry transition season and northward to savanna-forest transition and tropical seasonal forest regions in response to increased forest loss at the "Arc of Deforestation." Tropical forests, which have adapted to historically moist conditions, are less resilient and easily tip into an alternative state. Our results imply forest conservation and fire protection options to reduce the stress from positive feedback between forest loss, climate change, and fire.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Recent increases in annual, seasonal, and extreme methane fluxes driven by changes in climate and vegetation in boreal and temperate wetland ecosystems

Climate warming is expected to increase global methane (CH 4 ) emissions from wetland ecosystems. Although in situ eddy covariance (EC) measurements at ecosystem scales can potentially detect CH 4 flux changes, most EC systems have only a few years of data collected, so temporal trends in CH 4 remain uncertain. Here, we use established drivers to hindcast changes in CH 4 fluxes (FCH 4 ) since the early 1980s. We trained a machine learning (ML) model on CH 4 flux measurements from 22 [methane-producing sites] in wetland, upland, and lake sites of the FLUXNET-CH 4 database with at least two full years of measurements across temperate and boreal biomes. The gradient boosting decision tree ML model then hindcasted daily FCH 4 over 1981-2018 using meteorological reanalysis data. We found that, mainly driven by rising temperature, half of the sites (n = 11) showed significant increases in annual, seasonal, and extreme FCH 4 , with increases in FCH 4 of ca. 10% or higher found in the fall from 1981–1989 to 2010–2018. The annual trends were driven by increases during summer and fall, particularly at high-CH 4 -emitting fen sites dominated by aerenchymatous plants. We also found that the distribution of days of extremely high FCH 4 (defined according to the 95th percentile of the daily FCH 4 values over a reference period) have become more frequent during the last four decades and currently account for 10–40% of the total seasonal fluxes. The share of extreme FCH 4 days in the total seasonal fluxes was greatest in winter for boreal/taiga sites and in spring for temperate sites, which highlights the increasing importance of the non-growing seasons in annual budgets. Our results shed light on the effects of climate warming on wetlands, which appears to be extending the CH 4 emission seasons and boosting extreme emissions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modeling Plant Nutrient Acquisition Strategies Alters Projections of Carbon and Nitrogen Dynamics in Bioenergy Agroecosystems

Plant strategies to acquire nutrients from limited environments help shape ecosystem carbon (C) and nitrogen (N) cycling and response to environmental change. The effects of plant strategies on ecosystem dynamics are largely uncharacterized in bioenergy agroecosystems, where the impacts could determine bioenergy's ability to meet its sustainability goals of storing C and reducing N loss. We used FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms and Protection), a plant–microbe interaction model of coupled plant nutrient uptake and soil organic matter decomposition, to simulate the effects of nutrient acquisition strategies on soil microbial activity and ecosystem nutrient cycling in bioenergy feedstocks miscanthus (Miscanthus × giganteus) and sorghum (Sorghum bicolor (L.) Moench). We examined the model's ability to reproduce the relative effects of belowground nutrient uptake on microbial activity using a reanalysis of empirical data showing that miscanthus root exudation provoked a larger soil microbial response than sorghum. From baseline model simulations, we found that the ability of miscanthus to retranslocate N resulted in higher N uptake at a lower C cost than the sorghum/soybean rotation and that soil C and N pools increased under perennial (miscanthus) and decreased under annual (sorghum/soybean) cultivation. The model also predicted that greater root exudation increased soil C accumulation, highlighting the role of roots in forming stable soil C. Overall, the baseline model was unable to reproduce field observations of miscanthus root exudation stimulating microbial activity more than sorghum. To improve the model, we updated the soil microbial parameters in miscanthus to have faster decomposition, a higher C/N ratio, and greater carbon use efficiency. These changes improved the simulated soil microbial response to miscanthus root exudation, supporting the hypothesis that miscanthus soils foster a microbial community that is more responsive to root exudation than that of sorghum.

09 BIOMASS FUELS↗

DOE OSTI2024

Estimating UV-B, UV-Erithemic, and UV-A Irradiances From Global Horizontal Irradiance and MERRA-2 Ozone Column Information

The ground ultraviolet (UV) solar radiation is relevant due to its impacts on plastics degradation (mainly UVA) and on human health (UVB and erithemic UV (UVE)). UV ground measurements are not as ubiquitous as the relatively common global horizontal irradiance (GHI) measurements. Three simple models that estimate the UVA, UVB, and UVE components of solar irradiance from GHI and ozone column information are locally adjusted and validated. Five one-minute datasets from three sites in southeastern South America and two in the United States are used for simultaneous solar irradiance and UV data. All sites correspond to temperate mid-latitude regions. Simultaneous atmospheric total ozone column information is obtained from the reanalysis modern-era retrospective analysis for research and applications (MERRA-2) database for each site. Aside from locally adjusted models, average models with a single set of coefficients are also evaluated. For instance, the best average model is able to estimate UVE with a typical uncertainty below 12% and mean biases between +-3%, relative to the average of the measurements. Similar results are reported for the UVB and UVA components. These results, which can be useful in regions with similar climate and geography, provide a simple way to estimate UV irradiance under all-sky conditions with known uncertainty. This is an alternative to global satellite-based UV estimates, which can have high uncertainties at specific locations. Because MERRA-2 information has a global coverage, when coupled with good satellite-based estimates for GHI, UV irradiances can be estimated by this method over a large territory.

environmental UV radiation↗

DOE OSTI2025

Convective potential and fuel availability complement near-surface weather in regulating global wildfire activity

Wildfires are favored by hot, dry, windy, rainless conditions—this knowledge about fire weather informs both short-term forecast and long-term prediction of wildfire activity. Yet, wildfires rely on the availability of ignition and fuel, which are underrepresented in fire forecast and prediction practices. By analyzing satellite measurements and atmospheric reanalysis, here we show that near-surface weather only partially captures wildfire occurrence and intensity across the daily to seasonal timescales. Beyond near-surface weather, convection and fuel abundance play a complementary role in regulating burning processes. Specifically, enhanced atmospheric convection is identified for over 40% of the low-human-impact regions and 61% of global burnable areas during wildfire ignition and spreading periods. Meanwhile, 56% of shrublands and 54% of grasslands see higher fuel load with actual occurrence of fire. Our results highlight the role of convection and fuel in wildfire forecast, prompting a revisit of wildfire prediction under intertwined atmospheric and terrestrial changes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Snow-eater heat waves of the western United States

Abrupt snowmelt, triggered by rain-on-snow events or "snow-eater heat waves," can cause flooding, initiate or accelerate snow drought, and affect water availability. However, the characteristics (e.g., area, duration, and frequency), impacts, and trends of snow-eater heat waves have received little attention. To address this gap, we developed a method to identify snow-eater heat waves and estimate their melt potential using 20th Century Reanalysis version 3 air temperature data, the TempestExtremes algorithm, and an operational snowmelt model (SNOW-17) across 1850-2015. Melt season snow-eater heat waves typically last 3 to 5 days, with three to five events, doubling snowmelt rates. Seven of 11 spring superfloods are shown to coincide with snow-eater heat waves. Since the 1850s, snow-eater heat waves have increased in area and frequency, decreased in duration, and shifted earlier in the melt season. Incorporating snow-eater heat-wave impacts into SNOW-17 enhances extreme melt estimates, improving water management support tools.

Rhoades, Alan M↗

DOE OSTI2019

Comment on “Insulator-metal transition in dense fluid deuterium”

Celliers et al. (Reports, 17 August 2018, p. 677), in an attempt to reconcile differences in inferred metallization pressures, provide an alternative temperature analysis of the Knudson et al. experiments (Reports, 26 June 2015, p. 1455). We show here that this reanalysis implies an anomalously low specific heat for the metallic fluid that is clearly inconsistent with first-principles calculations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Addressing the dynamic nature of reference data: a new nucleotide database for robust metagenomic classification

Accurate metagenomic classification relies on comprehensive, up-to-date, and validated reference databases. While the NCBI BLAST Nucleotide (nt) database, encompassing a vast collection of sequences from all domains of life, represents an invaluable resource, its massive size—currently exceeding 10 12 nucleotides—and exponential growth pose significant challenges for researchers seeking to maintain current nt-based indices for metagenomic classification. Recognizing that no current nt-based indices exist for the widely used Centrifuge classifier, and the last public version currently available was released in 2018, we addressed this critical gap by leveraging advanced high-performance computing resources. We present new Centrifuge-compatible nt databases, meticulously constructed using a novel pipeline incorporating different quality control measures, including reference decontamination and filtering. These measures demonstrably reduce spurious classifications, as shown through our reanalysis of published metagenomic data where Plasmodium annotations were dramatically reduced using our decontaminated database, highlighting how database quality can significantly impact research conclusions. Through temporal comparisons, we also reveal how our approach minimizes inconsistencies in taxonomic assignments stemming from asynchronous updates between public sequence and taxonomy databases. These discrepancies are particularly evident in taxa such as Listeria monocytogenes and Naegleria fowleri, where classification accuracy varied significantly across database versions. These new databases, made available as pre-built Centrifuge indexes, respond to the need for an open, robust, nt-based pipeline for taxonomic classification in metagenomics. Applications such as environmental metagenomics, forensics, and clinical metagenomics, which require comprehensive taxonomic coverage, will benefit from this resource. Our work highlights the importance of treating reference databases as dynamic entities, subject to ongoing quality control and validation akin to software development best practices. This approach is crucial for ensuring accuracy and reliability of metagenomic analysis, especially as databases continue to expand in size and complexity.

59 BASIC BIOLOGICAL SCIENCES↗