SEARCH · Engineering Papers

Results for “RICHLAND”

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 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Automated Particle Analysis of Hanford Tank Wastes 241-AN-106, 241-AN-101, and 241-AW-105

This effort is involved in developing particle size and density distribution (PSDD) for the major phases within Hanford tank waste sludge. The object of this work has been to use automated particle analysis (APA) methods on the Scanning Electron Microscope (SEM) to provide PSDDs. The PSDD can be used to quantify the proportion of gibbsite particles likely to settle quickly and enable blending of wastes that may be prone to pipeline plugging such as the uranium phase, clarkeite with other wastes. The Hanford tank waste solids were analyzed with SEM combined with x-ray Energy Dispersive Spectroscopy (EDS) using APA. This method can allow thousands of individual particles to be characterized and so may provide a more representative view of the samples and information PSDDs. We characterized several as-received sludge samples and developed PSDD with APA. Sample preparation methods were found to impact the collected results, as it was more difficult to collect representative samples of the larger particles. Much of the material was also dried and cemented together that further impacted results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Coordinated Natural Gas and Electric Planning: Case Studies of Current Approaches and Practices

This paper examines how natural gas and electric utilities across eight U.S. states and two Canadian provinces are beginning to coordinate historically separate planning processes in response to growing system interdependencies, policy mandates, aging infrastructure, and changing customer energy choices. Electricity planning has long relied on robust integrated resource planning frameworks that weigh numerous objectives, risks, and costs. Natural gas planning, on the other hand, is typically less transparent and more narrowly focused on safety and system integrity. As economic, reliability, and policy drivers place new and shared demands on both systems, jurisdictions and utilities are experimenting with approaches such as coordinated forecasting, non pipeline alternatives, and combined planning pilots. A few have issued regulatory or statutory directives for greater data sharing and methodological alignment. Case studies from British Columbia, California, Colorado, Illinois, Massachusetts, Minnesota, New York, Québec, Rhode Island, and Washington illustrate a wide range of approaches to rethinking siloed planning. The paper identifies several common themes across the jurisdictions examined. It provides observations on the methods, processes, and organizational steps that may be needed in the future to address the challenges being faced by states and utilities. Lastly, it identifies some initial steps that states and utilities can consider if they would like to pursue more integrated, cost-effective, policy-aligned energy system planning.

03 NATURAL GAS↗

DOE OSTI2026

MatLib-1.2.2: Nuclear Material Properties Library

The U.S. Nuclear Regulatory Commission (NRC) uses the computer code Fuel Analysis under Steady-state and Transients (FAST) to model steady-state and transient fuel behavior to support regulatory decisions. FAST relies on a material properties library (MatLib) that contains the thermal and mechanical properties of the nuclear materials and coolants of interest to support the U.S. commercial nuclear industry. MatLib contains properties for a variety of nuclear fuels, cladding and other structural materials, gases, and coolants. In this document, material property correlations for the materials contained within MatLib are presented and discussed. When available, comparisons are made between the material property correlations and available data. Additionally, uncertainties are quantified on the material properties, which is then used by the NRC to support uncertainty quantification for best-estimate plus uncertainty safety evaluation reviews. This document describes MatLib-1.2.2, updated from MatLib-1.2.1 to include additional properties for metallic fuel. It is one of a series of documents on FAST; the other documents detail the models used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

CMIP6-based Multi-model Hydropower Projection over the Conterminous US, Version 1.1

This dataset presents a suite of hydropower projections for the conterminous United States (CONUS), derived from multiple downscaled and bias-corrected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). The CMIP6 GCMs are downscaled using either statistical (DBCCA) or dynamical (RegCM) approaches, based on two meteorological reference datasets (Daymet and Livneh). The resulting downscaled precipitation, temperature, and wind speed data are then used to drive two calibrated hydrologic models (VIC and PRMS), enabling simulations of projected future hydrologic responses across the CONUS. Simulated total runoff is subsequently employed to drive two hydropower models (WMP, now implemented as mosartwmpy-power, and WRES) to evaluate how climate change may affect future hydropower production for both federal and non-federal hydropower fleets. This dataset was developed to support the SECURE Water Act Section 9505 Assessment for the U.S. Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, see Broman et al. (2024), Thurber et al. (2024), Kao et al. (2022), and Zhou et al. (2023).

Voisin, Nathalie [Pacific Northwest National Labor↗

DOE OSTI2025

Data and Code for Understanding Generative AI Content with Embedding Models

This repository contains code for the experiments in the paper "Understanding Generative AI Content with Embedding Models". Constructing high-quality features is critical to any quantitative data analysis. While feature engineering was historically addressed by carefully hand-crafting data representations based on domain expertise, deep neural networks (DNNs) now offer a radically different approach. DNNs implicitly engineer features by transforming their input data into hidden feature vectors called embeddings. For embedding vectors produced by foundation models -- which are trained to be useful across many contexts -- we demonstrate that simple and well-studied dimensionality-reduction techniques such as Principal Component Analysis uncover inherent heterogeneity in input data concordant with human-understandable explanations. Of the many applications for this framework, we find empirical evidence that there is intrinsic separability between real samples and those generated by artificial intelligence (AI).

Vargas, Max [Pacific Northwest National Laboratory↗

DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗

DOE OSTI2025

Multiscale ACI Satellite Database

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

DOE OSTI2026

Sampling Microbial Dynamics in the Salish Sea Estuary: Evaluating Methods to Capture Cyanobacteria and Cyanophage

Introduction: Picocyanobacteria from the genera Prochlorococcus and Synechococcus thrive across the globe in aquatic environments, have relatively small genomes, and have growth dynamics regulated by both viral interactions and abiotic conditions, making them excellent model organisms for exploring host-pathogencoevolution. Methods: We developed and refined methods to sample and sequence cyanobacteria, cyanophages, and measured features of their abiotic environment. Results: The protocol described herein can successfully discriminate large-cell eukaryotic organisms, but size fractionation of picocyanobacteria appears to be affected by the presence of free DNA, multicellular structures, and abundant tycheposons. Our preferred final protocol from this exploratory effort included a combination of in-line and single vacuum flask filtrations, which reduced filtration processing time by over threefold in some cases compared to other tested methods, such as a fully in-line sequence or in-site filtrations. We successfully extracted an average of approximately 400–1200 ng for all filter fractions, with some variations between kits. Discussion: The protocol described herein can successfully discriminate large-cell eukaryotic organisms, but size fractionation of picocyanobacteria appears to be affected by the presence of free DNA, multicellular structures, and abundant tycheposons.

Salish Sea↗

DOE OSTI2025

Gene network centrality analysis identifies key regulators coordinating day-night metabolic transitions in Synechococcus elongatus PCC 7942 despite limited accuracy in predicting direct regulator-gene interactions

Synechococcus elongatus PCC 7942 is a model organism for studying circadian regulation and bioproduction, where precise temporal control of metabolism significantly impacts photosynthetic efficiency and CO 2 -to-bioproduct conversion. Despite extensive research on core clock components, our understanding of the broader regulatory network orchestrating genome-wide metabolic transitions remains incomplete. We address this gap by applying machine learning tools and network analysis to investigate the transcriptional architecture governing circadian-controlled gene expression. While our approach showed moderate accuracy in predicting individual transcription factor-gene interactions - a common challenge with real expression data - network-level topological analysis successfully revealed the organizational principles of circadian regulation. Our analysis identified distinct regulatory modules coordinating day-night metabolic transitions, with photosynthesis and carbon/nitrogen metabolism controlled by day-phase regulators, while nighttime modules orchestrate glycogen mobilization and redox metabolism. Through network centrality analysis, we identified potentially significant but previously understudied transcriptional regulators: HimA as a putative DNA architecture regulator, and TetR and SrrB as potential coordinators of nighttime metabolism, working alongside established global regulators RpaA and RpaB. This work demonstrates how network-level analysis can extract biologically meaningful insights despite limitations in predicting direct regulatory interactions. The regulatory principles uncovered here advance our understanding of how cyanobacteria coordinate complex metabolic transitions and may inform metabolic engineering strategies for enhanced photosynthetic bioproduction from CO 2 .

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Tracing priming effects in palsa peat carbon dynamics using a stable isotope-assisted metabolomics approach

Introduction: Peatlands store up to a third of global soil carbon, and in high latitudes their litter inputs are increasing and changing in composition under climate change. Although litter significantly influences peatland carbon and nutrient dynamics by changing the overall lability of peatland organic matter, the physicochemical mechanisms of this impact—and thus its full scope—remain poorly understood. Methods: We applied multimodal metabolomics (UPLC-HRMS, 1 H NMR) paired with 13 C Stable Isotope-Assisted Metabolomics (SIAM) to track litter carbon and its potential priming effects on both existing soil organic matter and carbon gas emissions. Through this approach, we achieved molecule-specific tracking of carbon transformations at unprecedented detail. Results: Our analysis revealed several key findings about carbon dynamics in palsa peat. Microbes responded rapidly to litter addition, producing a short-term increase in CO 2 emissions, fueled nearly exclusively by transformations of litter carbon. Litter inputs significantly contributed to the organic nitrogen pool through amino acids and peptide derivatives, which served as readily accessible nutrient sources for microbial communities. We traced the fate of plant-derived polyphenols including flavonoids like rutin, finding evidence of their degradation through heterocyclic C-ring fission, while accumulation of some polyphenols suggested their role in limiting overall decomposition. The SIAM approach detected subtle molecular changes indicating minimal and transient priming activity that was undetectable through conventional gas measurements alone. This transient response was characterized by brief microbial stimulation followed by rapid return to baseline metabolism. Pre-existing peat organic matter remained relatively stable; significant priming of its consumption was not observed, nor was its structural alteration. Discussion: This suggests that while litter inputs temporarily increase CO 2 emissions, they don’t sustain long-term acceleration of stored carbon decomposition or substantially decrease peat’s carbon store capacity. Our findings demonstrate how technological advancements in analytical tools can provide a more detailed view of carbon cycling processes in complex soil systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Pultrusion and Vitrimer Composites: Emerging Pathways for Sustainable Structural Materials

Pultrusion is a manufacturing process used to produce fiber-reinforced polymer composites with excellent mechanical, thermal, and chemical properties. The resulting materials are lightweight, durable, and corrosion-resistant, making them valuable in aerospace, automotive, construction, and energy sectors. However, conventional thermoset composites remain difficult to recycle due to their infusible and insoluble cross-linked structure. This review explores integrating vitrimer technology a novel class of recyclable thermosets with dynamic covalent adaptive networks into the pultrusion process. As only limited studies have directly reported vitrimer pultrusion to date, this review provides a forward-looking perspective, highlighting fundamental principles, challenges, and opportunities that can guide future development of recyclable high-performance composites. Vitrimers combine the mechanical strength (tensile strength and modulus) of thermosets with the reprocessability and reshaping of thermoplastics through dynamic bond exchange mechanisms. These polymers offer high-temperature reprocessability, self-healing, and closed-loop recyclability, where recycling efficiency can be evaluated by the recovery yield retention of mechanical properties and reuse cycles meeting the demand for sustainable manufacturing. Key aspects discussed include resin formulation, fiber impregnation, curing cycles, and die design for vitrimer systems. The temperature-dependent bond exchange reactions present challenges in achieving optimal curing and strong fiber–matrix adhesion. Recent studies indicate that vitrimer-based composites can maintain structural integrity while enabling recycling and repair, with mechanical performance such as flexural and tensile strength comparable to conventional composites. Incorporating vitrimer materials into pultrusion could enable high-performance, lightweight products for a circular economy. The remaining challenges include optimizing curing kinetics, improving interfacial adhesion, and scaling production for widespread industrial adoption.

Fiber composites↗

DOE OSTI2025

Stacked networks improve physics-informed training: Applications to neural networks and deep operator networks

Physics-informed neural networks and operator networks have shown promise for effectively solving equations modeling physical systems. However, these networks can happen to be difficult or impossible to train accurately. Here, we present a novel multifidelity framework for stacking physics-informed neural networks and operator networks that facilitates training. We successively build a chain of networks, where the output at one step can act as a low-fidelity input for training a longer chain, gradually increasing the expressivity of the learnt model. The equations imposed at each step of the iterative process can be the same or different (akin to simulated annealing). The iterative (stacking) nature of the proposed method allows us to learn progressively features of a solution which could have been hard to learn directly. Through benchmark problems including a nonlinear pendulum, the wave equation, and the viscous Burgers equation, we show how stacking can be used to improve the accuracy and reduce the required size of physics-informed neural networks and operator networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Computational Fluid Dynamics Modeling to Facilitate Qualification of Stack Sampling Probe Location

Computational fluid dynamics (CFD) modeling was used to help evaluate modifications to a radiological effluent stack and assist with establishing a stack sampling location that met the mixing criteria for qualification. Requirements for stack sampling location are listed in the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-2021 standard. Modeling was performed to help develop a suitable design for increasing building ventilation for the radiological effluent stack. The ANSI/HPS N13.1-2021 criteria for the air monitoring probe location are that the coefficient of variation of velocity uniformity, gaseous tracer uniformity, and particulate tracer uniformity must be less than or equal to 20%. Furthermore, no point in the sampling location may have a gaseous tracer concentration that varies from the mean concentration by more than 30%. Additionally, the flow angle at the sampling location must not be more than 20 degrees. The ANSI/HPS N13.1-2021 standard allows for models (physical or computational) to be employed to perform the full suite of qualification tests, followed by a more limited set of verification tests on the actual stack to qualify the stack sampling location. Here, a series of computational model simulations were employed to evaluate the stack qualification criteria. Significant time and re-source savings are achieved using CFD modeling. CFD modeling demonstrated that the stack meets the criteria at the sample probe location. Verification tests were performed on the modified stack to measure the velocity uniformity and flow angle at the stack sampling location, and results demonstrated that the CFD model results may be used to support the qualification of the stack sampling location.

Air Monitoring↗

DOE OSTI2025

African dust transported to Barbados in the wintertime lacks indicators of chemical aging

The chemical processing (“aging”) of mineral dust is thought to increase dust light-scattering efficiency, cloud droplet activation, and nutrient solubility. However, the extent of African dust aging during long-range transport to the western Atlantic is poorly understood. Here, we explore African dust aging in wintertime samples collected from Barbados when dust is transported at lower altitudes. Ion chromatography (IC) analysis indicates that the mass concentrations of bulk nitrate, sulfate, and oxalate increase, relative to background conditions, when African dust reaches Barbados after transatlantic transport, indicating dust aging. However, aerosol mixing-state analysis from computer-controlled scanning electron microscopy with energy dispersive X-ray spectroscopy (CCSEM/EDX) indicates that approximately 67 % of dust particles are internally mixed with sea salt, while only about 26 % of dust particles contain no internally mixed components. SEM/EDX elemental mapping and time-of-flight secondary ion mass spectrometry (TOF-SIMS) reveals that within internally mixed dust and sea salt particles, only sea salt components contain signs of aging as indicated by the loss of chloride and acquisition of nitrate and/or sulfate. Our results suggest that chemical aging may only modestly increase the solubility of nutrients in African dust during long-range transport. Because most dust that we measured was internally mixed with sea salt, chemical aging is not necessarily required to increase the hygroscopicity of dust, at least in the lower boundary layer. Further, our findings have implications for understanding the release of halogens from sea salts, which may be enhanced in internally mixed dust and sea salt particles.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Machine learning reveals strong grid-scale dependence in the satellite N d –LWP relationship

The relationship between cloud droplet number concentration ( N d ) and liquid water path (LWP) is highly uncertain yet crucial for determining the impact of aerosol-cloud interactions (ACI) on Earth's radiation budget. The N d -LWP relationship is examined using a machine learning (ML) random forest model applied to five years of satellite data at grid resolutions ranging from 10° to 0.05° in 12 distinct regions. In the subtropics, the shape of the N d -LWP relationship switches from an inverted-V at 1° grid-resolution to an “M” shape at 0.1° resolution with decreased $\frac{\textrm{dln⁡LWP}}{\textrm{dln⁡}N_d}$ sensitivity. Tropical and midlatitude regions generally show a more positive sensitivity. Cloud sampling and filtering also influence this slope, wherein the exclusion of thin clouds, as commonly performed to reduce retrieval uncertainty, leads to strongly negative sensitivity across all regions. Precipitation is primarily responsible for driving the strength of the sensitivity, with strong positive slopes in raining clouds and negative and/or neutral responses found in non-raining clouds. A new method to compute radiative forcing from the ML model shows a robust Twomey radiative forcing across all regions and grid resolutions. However, LWP and cloud fraction adjustments to the radiative forcing, which are ∼50 % or smaller than the Twomey effect, decrease to negligible values with higher spatial resolution data. As Earth system models move toward higher spatial resolutions in the future, evaluating the LWP and CF adjustment contributions to the radiative forcing budget at these finer resolutions will be essential for evaluation and model development.

Aerosol-Cloud Interactions↗

DOE OSTI2025

Transformation rate maps of dissolved organic carbon in the contiguous US

Riverine dissolved organic carbon (DOC) plays a vital role in regional and global carbon cycles. However, the processes of DOC conversion from soil organic carbon (SOC) and leaching into rivers are insufficiently understood, inconsistently represented, and poorly parameterized, particularly in land surface and Earth system models. As a first attempt to fill this gap, we propose a generic formula that directly connects SOC concentration with DOC concentration in headwater streams, where a single parameter, the transformation rate from SOC in the soil to DOC leaching flux (P r ), accounts for the overall processes governing SOC conversion to DOC and leaching from soils (along with runoff) into headwater streams. We then derive high-resolution P r maps over the contiguous US (CONUS) using SOC data from two different sources: the Harmonized World Soil Database v1.2 (HWSD) and SoilGrids 2.0. Both maps are developed following the same five major steps: (1) selecting independent catchments where observed riverine DOC data are available with reasonable quality; (2) estimating catchment-average SOC for the independent catchments; (3) estimating the P r values for these catchments based on the generic formula and catchment-average SOC; (4) developing a predictive model of P r with machine learning (ML) techniques and catchment-scale climate, hydrology, geology, and other attributes; and (5) deriving a national map of P r based on the ML model. For evaluation, we compare the DOC concentration derived using the P r map and the observed DOC concentration values at evaluation catchments. The resulting mean absolute scaled error and coefficient of determination are 0.73 and 0.47 for the HWSD-based model and 0.58 and 0.72 for the SoilGrids-based model, respectively, suggesting the effectiveness of the overall methodology. Efforts to constrain uncertainty and evaluate sensitivity of P r to different factors are discussed. To illustrate the use of such maps, we derive a riverine DOC concentration reanalysis dataset over CONUS. The two P r maps, robustly derived and empirically validated, lay a critical cornerstone for better simulating the terrestrial carbon cycle in land surface and Earth system models. Our findings not only set a foundation for improving our predictive understanding of the terrestrial carbon cycle at the regional and global scales, but also hold promises for informing policy decisions related to decarbonization and climate change mitigation. The data presented in this study are publicly available at https://doi.org/10.5281/zenodo.14563816 (Li et al., 2024).

54 ENVIRONMENTAL SCIENCES↗