SEARCH · Engineering Papers

Results for “mass balance”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

In this work, we report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for magneto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six vertices and the center of a hexagon covering approximately one-tenth of the surface area of a 9-ft-diameter spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, and Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging, leading to shock formation rather than potentially more desirable shockless jet merging; and (3) the liner Mach number can remain high (≳10), as required for plasma liners to be an effective compression driver. Secondly, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet mass balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. Ultimately, for both experiments, we present comparisons between experimental and synthetic data from three-dimensional hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗

DOE OSTI2022

N-dependent dynamics of root growth and nitrate and ammonium uptake are altered by the bacterium Herbaspirillum seropedicae in the cereal model Brachypodium distachyon

Nitrogen (N) fixation in cereals by root-associated bacteria is a promising solution for reducing use of chemical N fertilizers in agriculture. However, plant and bacterial responses are unpredictable across environments. We hypothesized that cereal responses to N-fixing bacteria are dynamic, depending on N supply and time. To quantify the dynamics, a gnotobiotic, fabricated ecosystem (EcoFAB) was adapted to analyse N mass balance, to image shoot and root growth, and to measure gene expression of Brachypodium distachyon inoculated with the N-fixing bacterium Herbaspirillum seropedicae. Phenotyping throughput of EcoFAB-N was 25–30 plants h -1 with open software and imaging systems. Herbaspirillum seropedicae inoculation of B. distachyon shifted root and shoot growth, nitrate versus ammonium uptake, and gene expression with time; directions and magnitude depended on N availability. Primary roots were longer and root hairs shorter regardless of N, with stronger changes at low N. At higher N, H. seropedicae provided 11% of the total plant N that came from sources other than the seed or the nutrient solution. The time-resolved phenotypic and molecular data point to distinct modes of action: at 5 mM NH 4 NO 3 the benefit appears through N fixation, while at 0.5 mM NH 4 NO 3 the mechanism appears to be plant physiological, with H. seropedicae promoting uptake of N from the root medium.Future work could fine-tune plant and root-associated microorganisms to growth and nutrient dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Clinical pharmacology information in regulatory submissions and labeling: A comparative analysis of orphan and non-orphan drugs approved by the FDA

Clinical pharmacology is an integral discipline supporting the development, regulatory evaluation, and clinical use of drugs for the treatment of both common and rare diseases. Here, we evaluated the recommendations and information available from select clinical pharmacology studies in the therapeutic product labeling of new molecular entities (NMEs) approved from 2017 to 2019 for both common and rare diseases. A total of 151 NMEs, including 72 orphan and 79 non-orphan drugs, were analyzed for recommendations and information available related to food–drug interaction, drug–drug interaction, renal impairment, hepatic impairment, QT assessment, and human radiolabeled mass balance studies using data collected from the original labeling and other regulatory documents. The analysis showed no statistically significant difference in the recommendations between orphan and non-orphan drugs except for renal impairment related recommendations in section 8 of the labeling. Although not significant, fewer hepatic impairment labeling recommendations were available for orphan drugs when compared with non-orphan drugs. At the time of initial approval, 79 postmarketing requirements (PMRs) and postmarketing commitments (PMCs) for 33 orphan drugs and 39 PMRs and PMCs for 19 non-orphan drugs were established; with most difference observed for drug–drug interaction, hepatic impairment, and QT assessment. Overall, although there was a trend for more labeling recommendations and fewer postmarketing studies and clinical trials for non-orphan drugs, there appeared to be no substantial differences in how these select clinical pharmacology studies are leveraged during the development and approval of orphan and non-orphan drugs.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Reduction of a Heme Cofactor Initiates N -Nitroglycine Degradation by NnlA

Linear nitramines are potentially carcinogenic environmental contaminants. The NnlA enzyme from Variovorax sp. strain JS1663 degrades the nitramine N-nitroglycine (NNG)—a natural product produced by some bacteria—to glyoxylate and nitrite (NO 2 –). Ammonium (NH 4 +) was predicted as the third product of this reaction. A source of nonheme Fe II was shown to be required for initiation of NnlA activity. However, the role of this Fe II for NnlA activity was unclear. This work reveals that NnlA contains a b-type heme cofactor. Reduction of this heme—either by a nonheme iron source or dithionite—is required to initiate NnlA activity. Therefore, Fe II is not an essential substrate for holoenzyme activity. Our data show that reduced NnlA (Fe II -NnlA) catalyzes at least 100 turnovers and does not require O 2 . Finally, NH 4 + was verified as the third product, accounting for the complete nitrogen mass balance. Size exclusion chromatography showed that NnlA is a dimer in solution. Additionally, Fe II -NnlA is oxidized by O 2 and NO 2 – and stably binds carbon monoxide (CO) and nitric oxide (NO). These are characteristics shared with heme-binding PAS domains. Furthermore, a structural homology model of NnlA was generated using the PAS domain from Pseudomonas aeruginosa Aer2 as a template. The structural homology model suggested His 73 is the axial ligand of the NnlA heme. Site-directed mutagenesis of His 73 to alanine decreased the heme occupancy of NnlA and eliminated NNG activity, validating the homology model. We conclude that NnlA forms a homodimeric heme-binding PAS domain protein that requires reduction for initiation of the activity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗

DOE OSTI2021

WaterTAP3 (The Water Technoeconomic Assessment Pipe-Parity Platform)

The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience.

Miara, Ariel↗

DOE OSTI2021

Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3)

Miara, Ariel↗

DOE OSTI2024

The Extraordinary March 2022 East Antarctica “Heat” Wave. Part II: Impacts on the Antarctic Ice Sheet

Abstract Between 15 and 19 March 2022, East Antarctica experienced an exceptional heat wave with widespread 30°–40°C temperature anomalies across the ice sheet. In Part I, we assessed the meteorological drivers that generated an intense atmospheric river (AR) that caused these record-shattering temperature anomalies. Here, we continue our large collaborative study by analyzing the widespread and diverse impacts driven by the AR landfall. These impacts included widespread rain and surface melt that was recorded along coastal areas, but this was outweighed by widespread high snowfall accumulations resulting in a largely positive surface mass balance contribution to the East Antarctic region. An analysis of the surface energy budget indicated that widespread downward longwave radiation anomalies caused by large cloud-liquid water contents along with some scattered solar radiation produced intense surface warming. Isotope measurements of the moisture were highly elevated, likely imprinting a strong signal for past climate reconstructions. The AR event attenuated cosmic ray measurements at Concordia, something previously never observed. Last, an extratropical cyclone west of the AR landfall likely triggered the final collapse of the critically unstable Conger Ice Shelf while further reducing an already record low sea ice extent. Significance Statement Using our diverse collective expertise, we explored the impacts from the March 2022 heat wave and atmospheric river across East Antarctica. One key takeaway is that the Antarctic cryosphere is highly sensitive to meteorological extremes originating from the midlatitudes and subtropics. Despite the large positive temperature anomalies driven from strong downward longwave radiation, this event led to huge amounts of snowfall across the Antarctic interior desert. The isotopes in this snow of warm airmass origin will likely be detectable in future ice cores and potentially distort past climate reconstructions. Even measurements of space activity were affected. Also, the swells generated from this storm helped to trigger the final collapse of an already critically unstable Conger Ice Shelf while further degrading sea ice coverage.

Wille, Jonathan D.↗

DOE OSTI2024

Data for Harnessing the Potential of Oilcane Waste Mud for Recovering Biobased Waxes

Oilcane is an engineered sugarcane with the ability to hyper-accumulate vegetative lipids. It is processed to obtain juice and bagasse as a potential substrate for the production of biofuels and biochemicals. The juice comprises solid particles that are separated as waste mud before the fermentation of the juice. In this study, the oilcane waste mud (OWM) generated from 1000 liters of oilcane juice was quantified and evaluated as a potential resource for recovering biobased waxes. Hexane and ethyl acetate were evaluated as two different solvents for extracting waxes from OWM followed by its purification using acetone. The extracted biobased wax samples were characterized for their chemical and thermal profiles which were then compared with commercial natural waxes. Detailed mass balance shows that 53.6 ± 2.6 kg (dry basis) of solid OWM gets generated upon processing 1000 L (~1068 kg) of oilcane juice. Hexane and ethyl acetate led to a crude wax yield of 25.6 ± 0.2% and 16.6 ± 0.4% (wt/wt, dry basis) respectively from OWM at the end of 8 h. The relative purification of the wax samples was reported in the range of 58%–65% (wt/wt). The purified OWM wax has a melting point of 74.7°C. The waste mud was valorized as a source of biobased waxes with characteristic chemical and thermal profiles comparable to commercial natural waxes (carnauba and beeswax). Considering the decline in the supply of petroleum wax in the future coupled with the switch to “greener” alternative products by consumers, OWM could be a valuable source of natural wax in the industrial sector reducing the dependence on petroleum waxes. Eventually, recovering biobased wax as a co-product from OWM would bring in an additional stream of revenue leading to the development of a zero-waste biorefinery based on bioenergy crops.

Biomass Analytics↗

DOE OSTI2025

Utah FORGE: EGS Reservoir Produced Fluids Geochemistry 2022-2024

This dataset contains geochemical analyses of produced fluids from the Utah FORGE site, specifically from wells 16A(78)-32, 16B(78)-32, and 58-32, collected during various stimulation, flowback, and circulation tests conducted between 2022 and 2024. The data contains element concentrations, pH, and dissolved gas compositions. Geochemical analyses for 2022 and 2023 were performed at the Brigham Young University geochemistry laboratory, while the 2024 results were obtained from Thermochem. Dry gas samples were collected using a mini-separator attached to the single-phase production line between the wellhead and separator, where fluid was flashed to atmospheric pressure. Gas concentrations were recalculated to a single phase reservoir liquid based on heat and mass balance expressions. Additional contextual information, including interpretations of geochemical trends and reservoir behavior, is available in the included report from Simmons et al. (2025), which was presented at the Stanford Geothermal Workshop in February 2025.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Overview of the MOSAiC expedition: Snow and sea ice

Year-round observations of the physical snow and ice properties and processes that govern the ice pack evolution and its interaction with the atmosphere and the ocean were conducted during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to September 2020. This work was embedded into the interdisciplinary design of the 5 MOSAiC teams, studying the atmosphere, the sea ice, the ocean, the ecosystem, and biogeochemical processes. The overall aim of the snow and sea ice observations during MOSAiC was to characterize the physical properties of the snow and ice cover comprehensively in the central Arctic over an entire annual cycle. This objective was achieved by detailed observations of physical properties and of energy and mass balance of snow and ice. By studying snow and sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the variability across scales can be considered. On-ice observations of in situ and remote sensing properties of the different surface types over all seasons will help to improve numerical process and climate models and to establish and validate novel satellite remote sensing methods; the linkages to accompanying airborne measurements, satellite observations, and results of numerical models are discussed. We found large spatial variabilities of snow metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly variable snow cover needs to be considered in more detail (in observations, remote sensing, and models) to better understand snow-related feedback processes. The ice pack revealed rapid transformations and motions along the drift in all seasons. The number of coupled ice–ocean interface processes observed in detail are expected to guide upcoming research with respect to the changing Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Spatiotemporal evolution of melt ponds on Arctic sea ice

Melt ponds on sea ice play an important role in the Arctic climate system. Their presence alters the partitioning of solar radiation: decreasing reflection, increasing absorption and transmission to the ice and ocean, and enhancing melt. The spatiotemporal properties of melt ponds thus modify ice albedo feedbacks and the mass balance of Arctic sea ice. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition presented a valuable opportunity to investigate the seasonal evolution of melt ponds through a rich array of atmosphere-ice-ocean measurements across spatial and temporal scales. In this study, we characterize the seasonal behavior and variability in the snow, surface scattering layer, and melt ponds from spring melt to autumn freeze-up using in situ surveys and auxiliary observations. We compare the results to satellite retrievals and output from two models: the Community Earth System Model (CESM2) and the Marginal Ice Zone Modeling and Assimilation System (MIZMAS). During the melt season, the maximum pond coverage and depth were 21% and 22 ± 13 cm, respectively, with distribution and depth corresponding to surface roughness and ice thickness. Compared to observations, both models overestimate melt pond coverage in summer, with maximum values of approximately 41% (MIZMAS) and 51% (CESM2). This overestimation has important implications for accurately simulating albedo feedbacks. During the observed freeze-up, weather events, including rain on snow, caused high-frequency variability in snow depth, while pond coverage and depth remained relatively constant until continuous freezing ensued. Both models accurately simulate the abrupt cessation of melt ponds during freeze-up, but the dates of freeze-up differ. MIZMAS accurately simulates the observed date of freeze-up, while CESM2 simulates freeze-up one-to-two weeks earlier. This work demonstrates areas that warrant future observation-model synthesis for improving the representation of sea-ice processes and properties, which can aid accurate simulations of albedo feedbacks in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

The two radiative states of the Arctic atmosphere and their impacts on the surface energy budget of sea ice

The surface energy budget (SEB) is a central regulator of Arctic climate and sea ice evolution, yet its processes remain poorly constrained due to sparse observations and complex, coupled surface-atmosphere interactions. This study leverages year-long measurements from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) to provide the most comprehensive assessment to date of the central Arctic SEB and its modulation by atmospheric variability. Ship- and ice-based observations from October 2019 to September 2020 were used to directly measure or tightly constrain each term of the SEB, leading to exceptional energetic closure with the seasonal snow and ice mass balance. The analysis reveals strong seasonal transitions in atmosphere-surface energy transfer that are modulated by the atmospheric state and constrained by the ability of the surface temperature to respond. Classification of the atmosphere into its two dominant radiative states—the semi-transparent (ST) and opaque (OP)—highlights the central role of synoptic-scale variability in clouds. The ST atmospheric state dominated the long winter ice growth season, with limited cloudiness supporting persistent surface radiative cooling and ice growth. The OP state, associated with liquid-containing or thick ice clouds, became dominant in spring, with the combination of increased solar heating and cloud surface longwave warming driving ice and snow melt. Eddy covariance versus bulk approaches for deriving surface turbulent heat fluxes provide vastly different perspectives on the role of turbulence in modulating the SEB. These results establish a high-quality benchmark dataset for Arctic SEB studies and demonstrate how the balance of atmospheric radiative states exerts a first-order control on the annual evolution of the sea ice. The findings have broad implications for advancing observing technologies, understanding Arctic amplification, improving climate models, and predicting future sea ice change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Raw Data published in PeerJ, 2021 for Crested Butte decomposition field study 2017-2019.

This data package contains text files that describe geochemical measurements collected from 2017-2019 during isolated conifer needle decomposition field studies in Crested Butte, Colorado. The geochemical measurements were collected across three elevations (2,800–3,500 m) ranging from montane to subalpine ecoregions. The data sets within include total carbon and nitrogen content and fourier-transform infrared spectroscopy (FTIR) results for the initial needles collected in 2016 and after decomposition in 2019. Data collection results from August 2020 are also included from when the experimental plots were removed to understand final concentrations of soil extractable carbon and nitrogen content as well as mass balances from litter bag deployments. Soil porewater results are also provided from 2017-2019 for DOC, TN, UV254, and specific UV absorbance (SUVA) analyses at 15 cm soil depth. Gas flux raw data provides CO2, CH4, N2O, and NH3 measurements above needle decomposition over the three study years. Finally, soil samples for microbial DNA extractions were collected from the upper soil depth. This raw data is available in the NCBI SRA database under SRA accession numbers PRJNA605259 and PRJNA715914. These data sets were generated to investigate the isolated decomposition of spruce and lodgepole conifer needles. The goal of this work was to determine the roles of elevation, soil type, seasonal changes in soil moisture, and snowmelt timing on litter decomposition processes. Results from this work are detailed in the reference paper "Effect of elevation, season and accelerated snowmelt on biogeochemical processes during isolated conifer needle litter decomposition. DOI: 10.7717/peerj.11926."

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Model files for estimating snow dynamics and stable water isotopes across the East River, CO.

A coupled hydrologic and snowpack stable water isotope model is used to assesses controls on isotopic inputs across the East River, Colorado, a large, mountainous basin. The hydrologic model uses the semi-empirical, spatially distributed and publicly available U.S. Geological Survey numerical code Precipitation-Modelling Runoff System (PRMS). Water and energy are tracked daily through the atmosphere, canopy and subsurface at a 100-m grid resolution. The isotope mass balance model follows previous work by Ala-aho et al. (2017) to track stable isotopes entering the soil system as snowmelt or rain. Water stores and fluxes needed for the isotope model use hydrologic model output for each timestep and model modeled grid location. This data package provides all files related to the hydrologic model (executable, input and output files), the isotopic model calibration and the historic isotopic model (water years 2015-2020). The isotope model source code and executable are provided. The file "readme.txt" includes information on the file structure, steps to run the model, and description of included folders.

54 ENVIRONMENTAL SCIENCES↗