Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOURCE TERMS”

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.

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

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Immunophenotypic characterization and therapeutics effects of human bone marrow- and umbilical cord-derived mesenchymal stromal cells in an experimental model of sepsis

Highlights: • Mesenchymal stromal cells (MSCs) may be a therapeutic option for sepsis. • Bone marrow (BM) vs. umbilical cord (UC) -MSCs induce different M1/M2 gene profiles. • In a CLP model of sepsis, only BM-MSCs improved survival and bacterial clearance. Sepsis is a complicated multi-system disorder characterized by a dysregulated host response to infection. Despite substantial progress in the understanding of mechanisms of sepsis, translation of these advances into clinically effective therapies remains challenging. Mesenchymal Stromal Cells (MSCs) possess immunomodulatory properties that have shown therapeutic promise in preclinical models of sepsis. The therapeutic effects of MSCs may vary depending on the source and type of these cells. In this comparative study, the gene expression pattern and surface markers of bone marrow-derived MSCs (BM-MSCs) and umbilical cord-derived MSCs (UC-MSCs) as well as their therapeutic effects in a clinically relevant mouse model of polymicrobial sepsis, cecal ligation and puncture (CLP), were investigated. The results showed remarkable differences in gene expression profile, surface markers and therapeutic potency in terms of enhancing survival and pro/anti-inflammatory responses between the two MSC types. BM-MSCs improved survival concomitant with an enhanced systemic bacterial clearance and improved inflammatory profile post CLP surgery. Despite some improvement in the inflammatory profile of the septic animals, treatment with UC-MSCs did not enhance survival or bacterial clearance. Overall, the beneficial therapeutic effects of BM-MSCs over UC-MSCs may likely be attributed to their pro-inflammatory function, and to some extent anti-inflammatory features, reflected in their gene expression pattern enhancing macrophage polarization to M1/M2 phenotypes resulting in a balanced pro- and anti-inflammatory response against polymicrobial sepsis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

QCOR; A Language Extension Specification for the Heterogeneous Quantum-Classical Model of Computation

Quantum computing (QC) is an emerging computational paradigm that leverages the laws of quantum mechanics to perform elementary logic operations. Existing programming models for QC were designed with fault-tolerant hardware in mind, envisioning stand-alone applications. However, the susceptibility of near-term quantum computers to noise limits their stand-alone utility. To better leverage limited computational strengths of noisy quantum devices, hybrid algorithms have been suggested whereby quantum computers are used in tandem with their classical counterparts in a heterogeneous fashion. This modus operandi calls out for a programming model and a high-level programming language that natively and seamlessly supports heterogeneous quantum-classical hardware architectures in a single-source-code paradigm. Motivated by the lack of such a model, we introduce a language extension specification, called QCOR, which enables single-source quantum-classical programming. Programs written using the QCOR library–based language extensions can be compiled to produce functional hybrid binary executables. After defining QCOR’s programming model, memory model, and execution model, we discuss how QCOR enables variational, iterative, and feed-forward QC. Additionally, QCOR approaches quantum-classical computation in a hardware-agnostic heterogeneous fashion and strives to build on best practices of high-performance computing. The high level of abstraction in the language extension is intended to accelerate the adoption of QC by researchers familiar with classical high-performance computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Net-Zero Carbon Microgrids

The microgrid concept has been effective in creating aggregations of distributed energy resources—generation, storage and loads—for resiliency, in the form of energy security. The success of microgrids in bringing energy security to a wide range of customers—from individual residences to commercial and industrial installations to military bases—has been exemplified during power disruptions and extended outages due to extreme weather events, cybersecurity attacks, and equipment failures. Now microgrids have an opportunity to meet the challenges of climate change and contribute to a carbon-free power delivery system. The transition to net-zero starts within microgrids themselves. In fact, today’s microgrids are largely dominated by generators using fossil fuels, natural gas and diesel, with high greenhouse gas emissions. In short, the transition to net-zero means replacing fossil fueled generators with renewable generation in microgrids. This transition is extended by including new dispatchable generation technologies that are 100% carbon-free and that offer additional advantage of a more-dependable and sustainable source of energy and power. Basically, the decarbonization of microgrids requires three elements: 1) maximizing generation from renewable energy resources, 2) management of storage and flexible loads to balance the variability and intermittency of renewable energy resources, and 3) introducing new clean power sources, including hydrogen-based generation and small modular reactors. This report affirms a need for specific focus by governmental agencies at national, regional, and local levels to establish technology, policy, and investment in this area. The intention of the Net-Zero Microgrid (NZM) Program is to inform these constituencies with cross-cutting research and tools for the reduction of GHG in microgrids – to net-zero in the near term eventually to zero in the longer term. The NZM Program is committed to achieving decarbonization for resiliency and for providing clean energy at the local or distribution level, from remote communities to underserved communities, and large industrial and military facilities. The NZM Planning and Design Platform is a core tool to be developed as an early deliverable of the NZM Program because only a fully integrated microgrid-design approach will ensure maximum carbon reduction in energy production.

13 HYDRO ENERGY↗
DOE OSTI2022

Sustainable Biofuels for Low-Carbon Maritime Transportation

The marine shipping sector heavily depends on fossil fuels and is one of the largest petroleum fuel consumers [1,2]. The annual global marine fuel consumption was estimated to be around 400 million metric tons in 2019 (2.5 billion barrels). Moreover, ocean shipping is one of the most significant contributors to sulfur oxides, nitrogen oxides, and particulate matter emissions. Global shipping contributes 13% of human-caused sulfur emissions and 2.6% of anthropogenic carbon dioxide emissions. As a major source of pollutant emissions, the marine industry faces several challenges related to emission regulations. The International Maritime Organization (IMO) has established a framework for reducing the carbon intensity of shipping: 40% reduction relative to 2008 levels by 2030 and 70% reduction by 2050. As the aviation sector, the maritime shipping sector is difficult to decarbonize through electrification. Biofuels offer the best opportunities for decarbonizing marine shipping in the near and medium-term. Advanced biofuels such as pyrolysis bio-oil offer the low-cost potential for meeting carbon reduction goals. For instance, the pyrolysis bio-oil exhibited promising marginal CO2 abatement costs at less than $100/tonne CO2-equivalent at a heavy fuel oil price greater than $1.10/gal [1]. As a potential biofuel option for low-carbon maritime shipping, this presentation focuses on a comparative techno-economic analysis (TEA) of bio-oils produced via a fast pyrolysis-based conversion pathway. The pathway converts a 50/50 blend of forest residues and clean pine to bio-oil via three process options: fast pyrolysis without vapor upgrading, and fast pyrolysis with vapor phase upgrading over ZSM-5 zeolite catalyst and Pt/TiO2 catalyst. The process configuration and operation variation led to different capital and operating costs, as well as the resulting raw bio-oil’s yield and quality, e.g., the water content, total acid number, and carboxylic acid number. The study also determined the minimum upgrading of bio-oils required to enable blending with very low sulfur fuel oil (VLSFO), with the associated costs reflected in TEA. This study shows that bio-oil could be a cost-effective fuel option for decarbonizing maritime shipping. Further research is required with respect to biofuel blending properties, such as compatibility with existing fuel system infrastructure and suitable engine performance.

BIOMASS FUELS↗
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Phase field modeling of dislocations and obstacles in InSb

We present a phase-field dislocation dynamics (PFDD) model informed by first-principle calculations to elucidate the competitive dislocation nucleation and propagation between the glide and shuffle sets in InSb diamond cubic crystal. The calculations are directly informed with generalized stacking fault energy curves on the (111) slip plane for both the “glide set,” with the smaller interplanar spacing, and the “shuffle set,” with the larger interplanar spacing. The formulation also includes elastic anisotropy and the gradient term associated with the dislocation core. The PFDD calculations show that under no stress the equilibrium structure of screw glide set dislocations dissociates into Shockley partials, while those of the shuffle set dislocations do not dissociate, remaining compact. The calculated dislocation core widths of these InSb dislocations agree well with the measured values for other semiconductor materials, such as Si and GaN. We find that a shuffle set dislocation emits from a dislocation source at an applied stress about three times smaller than that needed to emit leading and trailing partials successively on the glide set plane. Once the partial dislocations in the glide set are emitted, they propagate faster than the shuffle set perfect dislocation at the same stress level.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model-Free Probabilistic Forecasting of Nodal Voltages in Distribution Systems

As the penetration of distributed energy resources (DERs) into distribution systems increases, so does the interest in forecasting relevant system variables to help mitigate the associated challenges. One such challenge is the more frequent occurrence of excessive voltages in distribution systems with higher shares of DERs. Accurate and reliable estimates together with forecasts of system states (i.e., nodal voltages) will therefore play a key role in improving the utilization of these variable and uncertain sources while mitigating potential operational risks. Whilst recent literature has explored machine learning (ML) methods for voltage estimation and their extrapolation for a short-time period into the future, few have taken uncertainty quantification into account, and these methods have not yet been translated into operations. This paper discusses the advantages offered by probabilistic voltage forecasts and proposes a non-parametric Bayesian method suitable for forecasting nodal voltages at short-term time horizons while accounting for uncertainties in load and distributed photovoltaic (PV) generation. We demonstrate the value of the proposed Gaussian process (GP) model for a case study using historical forecasts and observation data.

distribution system↗
DOE OSTI2020

Indoor and Outdoor Test Results for 'DUSST', a Low-Cost, Low-Maintenance PV Soiling Sensor

Soiling can cause significant losses to photovoltaic systems, and therefore it is often measured for the purposes of predicting long-term energy forecasts or for monitoring real time performance and triggering maintenance events as needed. Currently, the most common soiling monitoring technologies are soiling stations that use the electrical outputs of a regularly cleaned PV device and of a naturally soiled PV device to quantify soiling. As part of a new class of low-cost and low-maintenance soiling stations NREL has previously presented 'DUSST'. DUSST projects a collimated monochromatic light source through a glass surface (exposed similarly to the PV modules that need to be monitored) and on to a light detector to measure the intensity of the transmitted light. As the glass surface naturally soils, the losses are quantified by comparing this soiled reading with a calibrated reading under baseline clean condition. This work presents the ongoing improvement of DUSST and the ongoing indoor and outdoor validation of this new soiling sensor.

14 SOLAR ENERGY↗
DOE OSTI2025

Reflections on the role of MOOSE in Unleashing a Golden Era of American Energy Innovation

As the use of artificial intelligence spreads across most professions, projections for the expected energy consumption of these virtual assistants/workers dictate a large increase in energy production. Recent published estimates refer to up to a necessary tripling of nuclear energy production~\cite{iaea}. % TODO add citations This is supported by the new administration, calling on February 5th to "unleash Commercial Nuclear Power in the United States" and to launch the "long-awaited American nuclear renaissance"~\cite{doe-order}. The Department of Energy (DOE) is tasked with "enabl[ing] the rapid deployment and export of next generation nuclear technology". In terms of research of development, we must seek to "unleash American Energy Innovation", by "focus[ing] our time and resources on technologies that will advance basic science, grow America’s scientific leadership, reduce costs for American families [...]". To that end, the Department will "prioritize affordable, reliable, and secure energy technologies, including fossil fuels, advanced nuclear, geothermal, and hydropower." Following this call to action, we elect to examine the role of the open-source Multiphysics Object Oriented Simulation Environment (MOOSE)~\cite{moose3} in meeting the stated goals.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Hydrodynamical Simulations of Black Hole Binary Formation in AGN Disks

Abstract We study close encounters between two single black holes (BHs) embedded in an AGN disk using a series of global 2D hydrodynamics simulations. We find that when the disk density is sufficiently high, bound BH binaries can be formed by the collision of their circum-single disks. Our analysis demonstrates that, after a BH pair passes the pericenter of their relative trajectory, post-collision gas drag may slow down the BHs, possibly forcing the two BHs to stay tightly bound. A binary formed by a close encounter can have a compact semimajor axis, large eccentricity, and retrograde orbital angular momentum. We provide a fitting formula that can accurately predict whether a close encounter can form a binary based on the gas mass and the incoming energy of the encounter. This fitting formula can be easily implemented in other long-term simulations that study the dynamical evolution of BHs in active galactic nucleus disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

FAST discovery of an extremely radio-faint millisecond pulsar from the Fermi -LAT unassociated source 3FGL J0318.1+0252

High sensitivity radio searches of unassociated γ-ray sources have proven to be an effective way of finding new pulsars. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) during its commissioning phase, we have carried out a number of targeted deep searches of Fermi Large Area Telescope (LAT) γ-ray sources. On February 27, 2018 we discovered an isolated millisecond pulsar (MSP), PSR J0318+0253, coincident with the unassociated γ-ray source 3FGL J0318.1+0252. PSR J0318+0253 has a spin period of 5.19 ms, a dispersion measure (DM) of 26 pc cm –3 corresponding to a DM distance of about 1.3 kpc, and a period-averaged flux density of (~11±2) µJy at L-band (1.05–1.45 GHz). Among all high energy MSPs, PSR J0318+0253 is the faintest ever detected in radio bands, by a factor of at least ~4 in terms of L-band fluxes. With the aid of the radio ephemeris, an analysis of 9.6 years of Fermi-LAT data revealed that PSR J0318+0253 also displays strong γ-ray pulsations. Follow-up observations carried out by both Arecibo and FAST suggest a likely spectral turn-over around 350 MHz. This is the first result from the collaboration between FAST and the Fermi-LAT teams as well as the first confirmed new MSP discovery by FAST, raising hopes for the detection of many more MSPs. Finally, such discoveries will make a significant contribution to our understanding of the neutron star zoo while potentially contributing to the future detection of gravitational waves, via pulsar timing array (PTA) experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Differential Trends in Iron Concentrations of Boreal Streams Linked to Catchment Characteristics

Increasing iron (Fe) concentrations have been reported for freshwaters across northern Europe over the last decades. This increase, together with elevated concentrations of dissolved organic carbon (DOC), leads to browning of freshwaters, which affects aquatic organisms, ecosystem functioning, biogeochemical cycles, and brings challenges to drinking water production. However, how such increasing trends in stream Fe concentrations reflect the contribution of different catchment sources remains poorly resolved. Here, we explored how catchment characteristics, that is, mires and coniferous soils, regulate spatial and temporal patterns of Fe in a boreal stream network. For this, we determined Fe speciation in riparian and mire soils, and studied temporal Fe dynamics in soil-water and stream-water over a span of 18 years. Positive Fe trends were found in the solution of the riparian soil, while no long-term trend was observed in the mire. These differences were reflected in stream-water, where three headwater streams dominated by coniferous cover also displayed positive Fe trends, whereas the mire dominated stream showed no trend. Surprisingly, the majority of higher order streams showed declining Fe trends, despite long-term increases in DOC. In addition, we found that an extreme drought event led to a prolonged release of Fe and DOC from the riparian soils, that could have long-term effects on stream Fe concentrations. Our results show that riparian forest soils can be major contributors to ongoing increases in freshwater Fe concentrations and that drought can further promote the release of Fe from organic soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

AMPP Newsletter (Jun 2021)

Welcome to our June 2021 issue of the AMPP newsletter. I hope the information you find within this newsletter allows you, as either a new addition to AMPP or as a long-term AMPP family member, to gain a better understanding of our organization and how we are working to meet not just our mission deliverables of today, but preparing for our future. In this newsletter issue the articles are dedicated to the principle of how we do work is as important as what we do. AMPP is known for excellence in our NNSA and DOE Missions execution. In this issue we are highlighting our Heat Source Final Assembly (HSFA) manufacturing efforts. Although what we produce for our HSFA customer is important their faith in us as an organization is just as much based on how we perform work to ensure that we safely and securely produce a quality product. We have added a new section to the newsletter called the “Group Leader Corner” where in this issues Chasity Kolar, the Group Leader for AMPP-3, discusses what it means to have a questioning attitude.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bringing Fault Detection and Diagnostics (FDD) Tools into the Mainstream: Retro Commissioning and Continuous Commissioning of HVAC and Refrigeration Systems

HVAC systems in the commercial building sector consume about 3.4 quads of energy annually. Of all the HVAC systems, packaged rooftop air–conditioning units (RTUs) provide cooling and heating for over 60 percent of the commercial building floorspace (about 90 billion ft2) in the U.S. and they are a significant source of energy consumption and peak demand. All HVAC systems suffer from faults that impact thermal comfort and increase energy consumption. There are several commercially available Automated Fault Detection and Diagnostic (AFDD) tools on the market that can detect and diagnose faults, and if those faults are corrected, can save significant national energy. However, there are multiple market barriers for these tools including lack of independent verification of their performance in the field in terms of their technical capabilities, ease of use and installation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Griffin Capability Improvements in Support of Ex-core Deep-Penetration Problems

Advanced reactor designs, especially portable reactors that are designed to be located closer to humans and operate autonomously, require the ability to accurately compute the ex-core neutron and gamma flux solutions in terms of shielding design optimization to reduce dose rates at the vessel boundary and detector signal prediction to drive the reactor control system. The Nuclear Energy Advanced Modeling and Simulation program has prioritized improvements to the Griffin discrete ordinates (SN) solver for deep-penetration problems in fiscal year 2025. Significant advancements have been made to the Griffin methodologies for solving ex-core deep-penetration problems for steady-state, fixed-source and transient calculations. This work presents the methodology improvements as well as a comprehensive demonstration with a Transient Test Reactor model and measurements.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Machine learning deciphers CO 2 sequestration and subsurface flowpaths from stream chemistry

Endmember mixing analysis (EMMA) is often used by hydrogeochemists to interpret the sources of stream solutes, but variations in stream concentrations and discharges remain difficult to explain. We discovered that machine learning can be used to highlight patterns in stream chemistry that reveal information about sources of solutes and subsurface groundwater flowpaths. The investigation has implications, in turn, for the balance of CO 2 in the atmosphere. For example, CO 2 -driven weathering of silicate minerals removes carbon from the atmosphere over ~106-year timescales. Weathering of another common mineral, pyrite, releases sulfuric acid that in turn causes dissolution of carbonates. In that process, however, CO 2 is released instead of sequestered from the atmosphere. Thus, understanding long-term global CO 2 sequestration by weathering requires quantification of CO 2 - versus H 2 SO 4 -driven reactions. Most researchers estimate such weathering fluxes from stream chemistry, but interpreting the reactant minerals and acids dissolved in streams has been fraught with difficulty. We apply a machine-learning technique to EMMA in three watersheds to determine the extent of mineral dissolution by each acid, without pre-defining the endmembers. The results show that the watersheds continuously or intermittently sequester CO 2 , but the extent of CO 2 drawdown is diminished in areas heavily affected by acid rain. Prior to applying the new algorithm, CO 2 drawdown was overestimated. The new technique, which elucidates the importance of different subsurface flowpaths and long-timescale changes in the watersheds, should have utility as a new EMMA for investigating water resources worldwide.

54 ENVIRONMENTAL SCIENCES↗