SEARCH · Engineering Papers

Results for “evaporation”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Advancing Continuous Manufacturing of Metal Organic Frameworks [Abstract]

InnaVenture has succeeded in developing a pilot-scale facility implementing a radically new concept for synthesis of metal-organic framework (MOF) materials based upon a patent-pending evaporative aerosolization-condensation process developed at PNNL. The successful synthesis of kilogram quantities of Ni-MOF-74 in InnaVenture’s facility represents a major step in the advancement of MOFs from being novelty materials to industrial adsorbents. Using this approach, fifteen kilograms of Ni-MOF-74 were synthesized for use in a demonstration adsorption chiller developed for the U.S. Navy. However, the evaporative condensation method for MOF synthesis has so far only been demonstrated for synthesizing MOFs under conditions where the reactants are part of a single homogeneous liquid phase solvent system. For InnaVenture to expand opportunities to support broader commercial applications of MOFs and related nanoporous materials (e.g., metal-aromatic frameworks (MAFs), covalent organic frameworks (COFs)), extension of the method to heterogeneous mixed phase conditions is essential. Examples of commercially promising MOFs that require such synthesis conditions include but are not limited to MIL-101, MIL-88B, MIL-100, MOF-801, MIL-160, MOF-841, and many others. InnaVenture has also informed PNNL that they need to address problems with the low density of MOF particles being generated in their fluidized reactor system. The purpose of this effort is to address both of these challenges.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Alternative Hanford ETF Grout Evaluation Report

The Hanford Effluent Treatment Facility (ETF) ammonia tolerant grout waste form was developed by the Vitreous State Laboratory (VSL) to stabilize ammonium and thereby prevent emission of ammonia vapor during processing of a final cementitious waste form. This ETF evaporator bottom waste is a concentrated sodium sulfate aqueous solution that also contains dissolved ammonium. The VSL waste form was produced in a two-step process in which struvite (MgNH 4 PO 4.6 H 2 O) is precipitated and then the resulting slurry is solidified in a cementitious matrix. SRNL researchers were tasked to evaluate the effect of compositional and struvite stability in simulated ETF evaporator brine. In 2022 to 2023, samples prepared by SRNL experienced significant expansion (up to 20 volume %) which was not observed by VSL.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Nanosecond Laser Drying in Wafer Foundries for Significant Energy Efficiency Improvement

Surface drying is a critical process in microelectronics wafer fabrication. In order to reduce and eliminate watermarks, many surface drying processes have been studied and used. Most of the wafer foundries currently use drying processes based on Isopropyl Alcohol (IPA), which causes serious safety (fire), health, environmental, and energy efficiency issues. The long-term goal of this project is to develop a fast, effective, chemical-free, and extremely energy efficient wafer drying technology that employs nanosecond laser-induced sub-surface evaporation and explosion. The objectives of this project are to 1) experimentally develop and verify the laser surface drying method, explore and identify how and to what extent various physical parameters influence the drying efficiency and effectiveness; and 2) understand and optimize the underlying thermodynamics in the proposed laser surface drying, including micro/nanoscale sub-surface superheating, evaporation, explosion, and photo/thermochemical bond breaking via atomistic modeling and experimental characterization. Toward these objectives, five tasks have been accomplished successfully, including 1) develop in-lab laser surface drying setup and study the effect of varied laser parameters, 2) use multiscale hybrid-modeling to study and understand the physics of water behavior during drying toward process optimization, 3) investigate the effect of drying conditions on the drying effectiveness, 4) develop a large-scale laser drying scanning system with automatic scanning, and 5) conduct large-scale simulation to study the effect of dopant level and develop strategies for system development toward industry deployment.

42 ENGINEERING↗

DOE OSTI2025

Advanced Brine Processing to Enable U.S. Lithium Independence (CRADA Report)

Current production of LiOH, which is needed to make Li-ion battery cathode active materials, utilizes a multistep process including solar evaporation, precipitation with Na 2 CO 3 and then conversion to LiOH using Ca(OH) 2 . This process requires a large amount of land area for solar evaporation, the right weather conditions, and chemicals for the conversion process that result in NaCl and CaCO 3 waste products. The production of Ca(OH) 2 is very energy intensive and evolves significant quantities of CO 2 . An alternative process flow utilizing direct lithium extraction techniques, followed by a chemical free conversion process can have benefits in reducing the needed land requirements and chemicals for traditional brine processing. There are many potential direct lithium extraction technologies that are currently being developed. The direct lithium extraction process from typical brine sources will produce a LiCl solution with some impurities including typically high concentrations of Na. This brine then needs to be converted to LiOH for use in battery cathode production. Ideally this conversion could occur without the use of additional chemicals. Electrochemistry can do this conversion either via electrolysis or bipolar membrane electrodialysis (BPED) to produce LiOH and HCl in solution. BPED utilized bipolar membranes to split water, which has a reduced potential as compared to splitting water at electrodes into hydrogen and oxygen gas. This reduced potential required results in a significant energy savings for BPED over electrolysis methods. This CRADA project aimed to develop such an integrated process using direct lithium extraction followed by BPED to produce a LiOH solution. That solution can then be crystallized into battery grade LiOH. In particular, Albemarle utilized a direct lithium extraction process to produce a concentrated LiCl solution that could be used for the BPED process. The BPED process was first tested using various LiCl solutions with impurity ions added at bench scale to understand the effects of impurities and determine processing parameters. Then testing was performed using the direct lithium extracted brine at the bench scale before scaling the process up. After the process was scaled up a long duration test was carried out to estimate the lifetime of the membranes, which is key to the economics of the BPED process.

25 ENERGY STORAGE↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Germanium-based planar heterostructures are emerging as versatile platforms for realizing quantum devices. In particular, planar Ge/SiGe quantum wells (QWs) host hole states with exceptionally high mobility and strong, electrically tunable spin–orbit interactions, enabling full electrical control of quantum information. A key requirement for these systems is the growth of microscopic, defect-free, atomically flat Ge QW heterostructures on relaxed or reverse-graded SiGe buffer layers, as well as on commercial Ge substrates. While several physical deposition techniques have demonstrated high-quality planar Ge/SiGe QWs, a major challenge remains: minimizing defect density typically requires high growth temperatures, which are incompatible with standard CMOS process flows. Here, we present a convenient low-temperature process for realizing high-quality planar SiGe/Ge heterostructures using a combination of thermal and electron-beam evaporation. We systematically map the effects of ex-situ and in-situ substrate preparation protocols, growth temperature, and post-deposition annealing conditions, and correlate these parameters with surface roughness and defect density. We find that in-situ oxide desorption conditions and post-annealing parameters have the most pronounced impact on improving surface quality. Under optimized conditions, we achieve atomically smooth surfaces with root-mean-square roughness σrms ≤ 10 Å and negligible defect density. Interestingly, thermally evaporated Ge layers exhibit oriented triangular crystallites that elongate upon post-annealing in the presence of high Ge vapor pressure. These results demonstrate that this simple, low-temperature deposition approach is a viable and effective route for achieving high-quality Ge QW heterostructures, with strong potential for scalable quantum computing and sensing applications.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2024

Lignin molecular weights of Populus trichocarpa residues after CELF pretreatment

Here we present a dataset of molecular weights of lignin from a woody energy crop (Populus trichocarpa) residues after a series of co-solvent enhanced lignocellulosic fractionation (CELF) pretreatment. The natural poplar variant GW-9947 from the Center for Bioenergy Innovation (CBI) was used. The poplar was knife milled and passed through a 1 mm sieve and CELF pretreatment was performed in a Parr autoclave reactor with 7.5 wt % solids loading, 0.5 wt% H2SO4 as catalyst at 150°C with various time. Tetrahydrofuran was added in a 1:1 mass ratio with water as the pretreatment solvent. Lignin was isolated from the pretreated samples after ball-milling in a porcelain jar with ceramic balls via Retsch PM 200 at 580 rpm for 2.5 h followed by enzymatic hydrolysis in acetate buffer (pH 4.8, 50 °C) for 48 h. The solid residue was isolated by centrifugation and hydrolyzed again with freshly added buffer and enzymes for another 48 h. After filtration, the solid residue was extracted twice with 96% (v/v) 1,4-dioxane/water mixture at room temperature overnight. The extracts were combined, rotary evaporated, and freeze-dried to recover lignin. The lignin samples were then derivatized in an acetic anhydride/pyridine (1:1, v/v) mixture and stirred at room temperature for 24 h. Ethanol was added to the reaction mixture, left for 30 min and then removed with a rotary evaporator. The addition and removal of ethanol was repeated at least 3 times until all traces of acetic acid were removed. Acetylated lignin samples were then dissolved in tetrahydrofuran (THF) at a concentration of 1.0 mg/mL. The molecular weight of acetylated lignin was measured by a gel permeation chromatography (GPC) on a PSS-Polymer Standards Service (Warwick, RI, USA) GPC SECurity 1200 system featuring Agilent HPLC 1200 components equipped with four Waters Styragel columns (HR1, HR2, HR4 and HR6) and an UV detector (270 nm). Tetrahydrofuran was used as the mobile phase and flow rate was 0.3 mL/min. The Polymer Standards Service WinGPC Unity software (Build 6807) was used for data processing for all the samples. The data provides information about the effects of CELF pretreatment time at 150 ºC on lignin molecular weights.

09 BIOMASS FUELS↗

DOE OSTI2022

Heterogeneous Ice Growth in Micron-Sized Water Droplets Due to Spontaneous Freezing

Understanding how ice nucleates and grows into larger crystals is of crucial importance for many research fields. The purpose of this study was to shed light on the phase and structure of ice once a nucleus is formed inside a metastable water droplet. Wide-angle X-ray scattering (WAXS) was performed on micron-sized droplets evaporatively cooled to temperatures where homogeneous nucleation occurs. We found that for our weak hits ice grows more cubic compared to the strong hits that are completely hexagonal. Due to efficient heat removal caused by evaporation, we propose that the cubicity of ice at the vicinity of the droplet’s surface is higher than for ice formed within the bulk of the droplet. Moreover, the Bragg peaks were classified based on their geometrical shapes and positions in reciprocal space, which showed that ice grows heterogeneously with a significant population of peaks indicative of truncation rods and crystal defects. Frequent occurrences of the (100) reflection with extended in-planar structure suggested that large planar ice crystals form at the droplet surface, then fracture into smaller domains to accommodate to the curvature of the droplets. Planar faulting due to misaligned domains would explain the increased cubicity close to the droplet surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Comparative Energetics of Various Membrane Distillation Configurations and Guidelines for Design and Operation

This paper presents a comparative performance study of single-stage desalination processes with major configurations of membrane distillation (MD) modules. MD modules covered in this study are (a) direct contact MD (DCMD), (b) vacuum MD (VMD), (c) sweeping gas MD (SGMD), and (d) air gap MD (AGMD). MD-based desalination processes are simulated with rigorous theoretical MD models supported by molecular thermodynamic property models for the accurate calculation of performance metrics. The performance metrics considered in MD systems are permeate flux and energy efficiency, i.e., gained output ratio (GOR). A general criterion is established to determine the critical length of these four MDs (at fixed width) for the feasible operation of desalination in a wide range of feed salinities. The length of DCMD and VMD is restricted by the feed salinity and permeate flux, respectively, while relatively large AGMD and SGMD are allowed. The sensitivity of GOR flux with respect to permeate conditions is investigated for different MD configurations. AGMD outperforms other configurations in terms of energy efficiency, while VMD reveals the highest permeate production. With larger MD modules, utilization of thermal energy supplied by the hot feed for evaporation is in the order of VMD > AGMD > SGMD > DCMD. Simulation results highlight that energy efficiency of the overall desalination process relies on the efficient recovery of spent for evaporation, suggesting potential improvement in energy efficiency for VMD-based desalination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Effects of Dark Matter Self-interactions on Sagittarius and its Stream

This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.

Milky Way dark matter halo↗

DOE OSTI2021

Measurement report: Cloud condensation nuclei activity and its variation with organic oxidation level and volatility observed during an aerosol life cycle intensive operational period (ALC-IOP)

Abstract. Cloud condensation nuclei (CCN) spectrum and the CCN activated fraction of size-resolved aerosols (SR-CCN) were measured at a rural site on Long Island during the Department of Energy (DOE) aerosol life cycle intensive operational period (ALC-IOP) from 15 July to 15 August 2011. During the last week of the ALC-IOP, the dependence of the activated fraction on aerosol volatility was characterized by sampling downstream of a thermodenuder (TD) operated at temperatures up to 100 ∘C. Here we present aerosol properties, including aerosol total number concentration, CCN spectrum, and the CCN hygroscopicity, for air masses of representative origins during the ALC-IOP. The hygroscopicity of organic species in the aerosol is derived from CCN hygroscopicity and chemical composition. The dependence of organic hygroscopicity on the organic oxidation level (e.g., atomic O:C ratio) agrees well with theoretical predictions and results from previous laboratory and field studies. The derived κorg and O:C ratio first increase as TD temperature increases from 20 ∘C (i.e., ambient temperature) to 50 or 75 ∘C and then decrease as TD temperature further increases to 100 ∘C. The initial increases of O:C and κorg with TD temperature below 50 ∘C are likely due to evaporation of more volatile organics with relatively lower O:C and hygroscopicity such as primary organic aerosol. At the high TD temperatures, the decreases of O:C and κorg indicate that evaporated organics were more oxygenated and had lower molecular weights. These trends are different from previous laboratory experiments and field observations, which reported that organic O:C increased monotonically with increasing TD temperature, whereas κorg decreased with the TD temperature. One possible reason is that previous studies were either focused on laboratory-generated secondary organic aerosol (SOA) or based on field observations at locations more dominated by SOA.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Dilution impacts on smoke aging: evidence in Biomass Burning Observation Project (BBOP) data

Abstract. Biomass burning emits vapors and aerosols into the atmosphere that can rapidly evolve as smoke plumes travel downwind and dilute, affecting climate- and health-relevant properties of the smoke. To date, theory has been unable to explain observed variability in smoke evolution. Here, we use observational data from the Biomass Burning Observation Project (BBOP) field campaign and show that initial smoke organic aerosol mass concentrations can help predict changes in smoke aerosol aging markers, number concentration, and number mean diameter between 40–262 nm. Because initial field measurements of plumes are generally >10 min downwind, smaller plumes will have already undergone substantial dilution relative to larger plumes and have lower concentrations of smoke species at these observations closest to the fire. The extent to which dilution has occurred prior to the first observation is not a directly measurable quantity. We show that initial observed plume concentrations can serve as a rough indicator of the extent of dilution prior to the first measurement, which impacts photochemistry, aerosol evaporation, and coagulation. Cores of plumes have higher concentrations than edges. By segregating the observed plumes into cores and edges, we find evidence that particle aging, evaporation, and coagulation occurred before the first measurement. We further find that on the plume edges, the organic aerosol is more oxygenated, while a marker for primary biomass burning aerosol emissions has decreased in relative abundance compared to the plume cores. Finally, we attempt to decouple the roles of the initial concentrations and physical age since emission by performing multivariate linear regression of various aerosol properties (composition, size) on these two factors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

An attribution of the low single-scattering albedo of biomass burning aerosol over the southeastern Atlantic

Abstract. Aerosol over the remote southeastern Atlantic is some of the most sunlight-absorbing aerosol on the planet: the in situ free-tropospheric single-scattering albedo at the 530 nm wavelength (SSA530 nm) ranges from 0.83 to 0.89 within ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) aircraft flights from late August–September. Here we seek to explain the low SSA. The SSA depends strongly on the black carbon (BC) number fraction, which ranges from 0.15 to 0.4. Low organic aerosol (OA)-to-BC mass ratios of 8–14 and modified combustion efficiency values >0.975 point indirectly to the dry, flame-efficient combustion of primarily grass fuels, with back trajectories ending in the miombo woodlands of Angola. The youngest aerosol, aged 4–5 d since emission, occupied the top half of a 5 km thick plume sampled directly west of Angola with a vertically consistent BC:ΔCO (carbon monoxide) ratio, indicating a homogenization of the source emissions. The younger aerosol, transported more quickly off of the continent by stronger winds, overlaid older, slower-moving aerosol with a larger mean particle size and fraction of BC-containing particles. This is consistent with ongoing gas condensation and the coagulation of smaller non-BC particles upon the BC-containing particles. The particle volumes and OA:BC mass ratios of the older aerosol were smaller, attributed primarily to evaporation following fragmentation, instead of dilution or thermodynamics. The CLARIFY (CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017) aircraft campaign sampled aerosols that had traveled further to reach the more remote Ascension Island. CLARIFY reported higher BC number fractions, lower OA:BC mass ratios, and lower SSA yet larger mass absorption coefficients compared to this study's. Values from one ORACLES 2017 flight, held midway to Ascension Island, are intermediate, confirming the long-range changes. Overall the data are most consistent with continuing oxidation through fragmentation releasing aerosols that subsequently enter the gas phase, reducing the OA mass, rather than evaporation through dilution or thermodynamics. The data support the following best fit: SSA530nm=0.801+0055⋅(OA:BC) (r=0.84). The fires of southern Africa emit approximately one-third of the world's carbon; the emitted aerosols are distinct from other regional smoke emissions, and their composition needs to be represented appropriately to realistically depict regional aerosol radiative effects.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Photoaging of phenolic secondary organic aerosol in the aqueous phase: evolution of chemical and optical properties and effects of oxidants

Abstract. While gas-phase reactions are well established to have significant impacts on the mass concentration, chemical composition, and optical properties of secondary organic aerosol (SOA), the aqueous-phase aging of SOA remains poorly understood. In this study, we performed a series of long-duration photochemical aging experiments to investigate the evolution of the composition and light absorption of the aqueous SOA (aqSOA) from guaiacyl acetone (GA), a semivolatile phenolic carbonyl that is common in biomass burning smoke. The aqSOA was produced from reactions of GA with hydroxyl radical (•OH-aqSOA) or a triplet excited state of organic carbon (3C∗-aqSOA) and was then photoaged in water under conditions that simulate sunlight exposure in northern California for up to 48 h. The effects of increasing aqueous-phase •OH or 3C∗ concentration on the photoaging of the aqSOA were also studied. High-resolution aerosol mass spectrometry (HR-AMS) and UV–Vis spectroscopy were utilized to characterize the composition and the light absorptivity of the aqSOA and to track their changes during aging. Compared to •OH-aqSOA, the 3C∗-aqSOA is produced more rapidly and shows less oxidation, a greater abundance of oligomers, and higher light absorption. Prolonged photoaging promotes fragmentation and the formation of more volatile and less light-absorbing products. More than half of the initial aqSOA mass is lost, and substantial photobleaching occurs after 10.5 h of prolonged aging under simulated sunlight illumination for 3C∗-aqSOA and 48 h for •OH-aqSOA. By performing positive matrix factorization (PMF) analysis of the combined HR-AMS and UV–Vis spectral data, we resolved three generations of aqSOA with distinctly different chemical and optical properties. The first-generation aqSOA shows significant oligomer formation and enhanced light absorption at 340–400 nm. The second-generation aqSOA is enriched in functionalized GA species and has the highest mass absorption coefficients in 300–500 nm, while the third-generation aqSOA contains more fragmented products and is the least light absorbing. These results suggest that intermediately aged phenolic aqSOA is more light absorbing than other generations, and that the light absorptivity of phenolic aqSOA results from a competition between brown carbon (BrC) formation and photobleaching, which is dependent on aging time. Although photoaging generally increases the oxidation of aqSOA, a slightly decreased O/C of the •OH-aqSOA is observed after 48 h of prolonged photoaging with additional •OH exposure. This is likely due to greater fragmentation and evaporation of highly oxidized compounds. Increased oxidant concentration accelerates the transformation of aqSOA and promotes the decay of BrC chromophores, leading to faster mass reduction and photobleaching. In addition, compared with •OH, photoaging by 3C∗ produces more low-volatility functionalized products, which counterbalances part of the aqSOA mass loss due to fragmentation and evaporation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

An ensemble of 48 physically perturbed model estimates of the 1/8° terrestrial water budget over the conterminous United States, 1980–2015

Terrestrial water budget (TWB) data over large domains are of high interest for various hydrological applications. Spatiotemporally continuous and physically consistent estimations of TWB rely on land surface models (LSMs). As an augmentation of the operational North American Land Data Assimilation System Phase 2 (NLDAS-2) four-LSM ensemble, this paper describes a dataset simulated from an ensemble of 48 physics configurations of the Noah LSM with multi-physics options (Noah-MP). The 48 Noah-MP physics configurations are selected to give a representative cross-section of commonly used LSMs for parameterizing runoff, atmospheric surface layer turbulence, soil moisture limitation on photosynthesis, and stomatal conductance. The dataset spans from 1980 to 2015 over the conterminous United States (CONUS) at a monthly temporal resolution and a 1/8° spatial resolution. The dataset variables include total evapotranspiration and its constituents (canopy evaporation, soil evaporation, and transpiration), runoff (the surface and subsurface components), as well as terrestrial water storage (snow water equivalent, four-layer soil water content from the surface down to 2 m, and the groundwater storage anomaly). The dataset is available at https://doi.org/10.5281/zenodo.7109816. Evaluations carried out in this study and previous investigations show that the ensemble per forms well in reproducing the observed terrestrial water storage, snow water equivalent, soil moisture, and runoff. Noah-MP complements the NLDAS models well, and adding Noah-MP consistently improves the NLDAS es timations of the above variables in most areas of CONUS. Besides, the perturbed-physics ensemble facilitates the identification of model deficiencies. The parameterizations of shallow snow, spatially varying groundwater dynamics, and near-surface atmospheric turbulence should be improved in future model versions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Soil-related developments of the Biome-BGCMuSo v6.2 terrestrial ecosystem model

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from -0.047 to -0.007 m 3 m -3 for SWC and from -0.68 to -0.2 mm d -1 for soil evaporation). Simulations related to nitrogen balance and soil CO 2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N 2 O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N 2 O efflux simulation bias was between -0.13 and -0.1 mgNm -2 d -1 , and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO 2 efflux simulations bias was 0.04–0.17 gCm -2 d -1 , while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

58 GEOSCIENCES↗

DOE OSTI2025

Causal relationships of vegetation productivity with root zone water availability and atmospheric dryness at the catchment scale

Abstract. This study explores the causal relationships between catchment water availability, vapor pressure deficit, and gross primary productivity (GPP) across 341 catchments in the contiguous US. Seasonal climatic, hydrological, and vegetation characteristics were represented using the Horton index, ecological aridity index, evaporative fraction index, and carbon uptake efficiency. Statistical methods, including circularity statistics, correlation analysis, and causality tests, were employed to determine the complex interactions between catchment wetness, atmospheric dryness, and vegetation carbon uptake. The results revealed a maximum lag of 2 months in the intra-annual variability of catchment water supply–productivity and atmospheric water demand–productivity relationships, with hysteresis patterns varying with the catchment's hydrological characteristics. In catchments not permanently under water-limited or energy-limited conditions, vegetation experiences hydrological stress during the peak growing period, coinciding with the highest gross primary productivity and carbon uptake efficiency being out of phase with the Horton index and in phase with the evaporative fraction index. Causality analysis highlights strong temporal continuity in GPP seasonal characteristics, with a cause–effect relationship between catchment water supply, atmospheric demand, and vegetation productivity spanning a maximum of 2 months. These findings underscore the need for a comprehensive functional framework that integrates catchment water supply, atmospheric demand, and vegetation productivity to enhance our understanding and predictive capabilities with regard to ecosystem responses to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

IM3 Projected US Data Center Locations

IM3 Projected US Data Center Locations This dataset contains model projections of new data center facilities in the contiguous United States (CONUS) through 2035 using the CERF – Data Centers model. Data center locations are modeled across four data center electricity demand growth scenarios (low, moderate, high, higher) and five market gravity scenarios (0%, 25%, 50%, 75%, 100%). Projected locations are intended to be regional representations of feasible siting locations in the future to assess potential grid and water stress impacts. The data center load growth scenarios correspond with the rates outlined in EPRI (2024) and include 3.71%, 5%, 10%, and 15% annual growth of electricity demand for data centers from 2023 values in 37 states across the CONUS. Market gravity scenarios correspond to the relative importance of proximity to data center markets or high population areas compared to locational cost in the siting algorithm. 0% market gravity means that siting decisions were entirely determined by the locational cost in each feasible location. 100% market gravity means that only market proximity was considered when siting. Other scenarios have weight placed on both components where total weight always equals 100%. Locational cost is dependent on facility cooling type and corresponding electricity cost, taxes, and other factors. Facility cooling type is spatially determined where high water stress and/or areas with high summer wet bulb temperatures are assumed to operate with mechanical cooling for a higher fraction of the year rather than evaporative cooling. Feasible data center siting areas are based on geospatial suitability raster data developed with open-source information. The following areas are excluded from siting: Areas within 300 m of a federal airport runway Waterbodies Areas with slope >16% Areas susceptible to sinkholes High coastal or inland flood risk areas Local, state, and federal parks, leisure areas, and cemeteries Areas >2 km away from electric substations Areas >5 km away from a municipal water supplier service area Areas >2 km away from high-speed fiber provider service territory Protected Areas Database of the United States (PAD-US) areas Railroads, major roadways, and minor roadways Military areas and training grounds NLCD developed lands Areas >0.8 km (0.5 miles) from NLCD developed lands Because we use open-source information, proprietary information that can influence siting decisions such as individual tax agreements with cities, detailed fiber line connectivity, electric grid power capacity agreements, and others, are not currently accounted for in the modeling process. Using specific building locations and footprints in the dataset for local planning purposes is not advised. Technical Information Geospatial data is provided in geojson format using the Albers Equal Area Conic (ESRI:102003) coordinate reference system. The datasets contain the following parameters: id - unique identification number within given scenario file growth_scenario – data center demand growth scenario market_gravity_weight – market gravity weight scenario (%) region – name of region (i.e., US State) total_cost_million_usd – locational siting cost ($million) campus_size_square_ft – total land acquired for data center facility (square ft) data_center_it_power_mw – IT power of data center facility (MW) mechanical_cooling_frac – fraction of year when data center uses mechanical cooling system water_cooling_frac– fraction of year when data center uses evaporative cooling system cooling_energy_demand_mwh – total annual facility energy demand for cooling (MWh) cooling_water_demand_mgy – total annual facility water demand for cooling (MG) cooling_water_consumption_mgy – total annual facility water consumed (MG) normalized_locational_cost – normalized total locational cost score for location normalized_gravity_score – normalized market gravity score for location weighted_siting_score – total weighted siting score of locational cost and gravity score geometry – polygon geometry of facility Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4.0 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall (ORCID:0000000328077088)↗