Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux analysis”

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,045 records · Page 58

Atmospheric 81 Kr as an Integrator of Cosmic-Ray Flux on the Hundred-Thousand-Year Time Scale

The atmospheric abundance of 81 Kr is a global integrator of cosmic rays. It is insensitive to climate shifts, geographical variations, and short-term solar cycle activity, making it an ideal standard to test models of cosmic-ray flux on the time scale of 10 5 years. In this work, we present the first calculation of absolute 81 Kr production rates in the atmosphere, and a measurement of the atmospheric 81 Kr/Kr abundance via the Atom Trap Trace Analysis method. The measurement result significantly deviates from previously reported values. The agreement between measurement and model prediction supports the current understanding of the production mechanisms. Additionally, the calculated 81 Kr atmospheric inventory over the past 1.5 Myr provides a more accurate input function for radiokrypton dating.

58 GEOSCIENCES↗

Overview of the Nephele Perturbed Parameter Ensemble for Aerosol‐Cloud Interactions in E3SMv3

Aerosol-cloud interactions (aci) are the leading source of uncertainty in inferring climate sensitivity from the historical record. Earth system models (ESMs) struggle to represent aci because the processes responsible for these phenomena occur at much finer time and space scales than can be resolved by any ESM. Observational constraints provide key benchmarks to test ESMs, but cannot be used alone to fully understand aci processes except in very specific cases where causality is controlled; some degree of modeling is required to infer aci and estimate radiative forcing. Here, we generate and characterize a perturbed parameter ensemble (PPE) in version 3 of the Energy Exascale ESM (E3SMv3). We perturb 25 parameters that govern aci processes over 250 members and integrate the model over present-day and preindustrial aerosol emissions. We find that the process representation in E3SMv3 is flexible and can generate global-mean effective radiative forcings due to aci (ERFaci) ranging from −3.0 to +0.9 W m −2 . The positive ERFaci values simulated by a portion of the PPE are implausible and result from parameter combinations that produce unrealistic top-of-atmosphere energy fluxes. While global-mean cloud droplet number concentration always increases in response to anthropogenic aerosol, cloud liquid water path can both increase and decrease, suggesting that precipitation suppression is not the only aerosol-cloud adjustment represented by E3SMv3. Analysis of which processes control liquid cloud adjustment in the PPE points toward stratiform precipitation processes and aerosol activation, which is consistent with many previous ESMs, as well as the new two-moment convective cloud microphysics in E3SMv3.

Nugent, Jacqueline M. [Univ. of Wyoming, Laramie, ↗

Temporal changes in global soil respiration since 1987

As the second-largest terrestrial carbon (C) flux, soil respiration (R S ) has been stimulated by climate warming. However, the magnitude and dynamics of such stimulations of soil respiration are highly uncertain at the global scale, undermining our confidence in future climate projections. Here, we present an analysis of global R S observations from 1987–2016. R S increased (P < 0.001) at a rate of 27.66 g C m -2 yr -2 (equivalent to 0.161 Pg C yr -2 ) in 1987–1999 globally but became unchanged in 2000–2016, which were related to complex temporal variations of temperature anomalies and soil C stocks. However, global heterotrophic respiration (R h ) derived from microbial decomposition of soil C increased in 1987–2016 (P < 0.001), suggesting accumulated soil C losses. Given the warmest years on records after 2015, our modeling analysis shows a possible resuscitation of global R S rise. This study of naturally occurring shifts in R S over recent decades has provided invaluable insights for designing more effective policies addressing future climate challenges.

54 ENVIRONMENTAL SCIENCES↗

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗

Constraints on neutrino emission from nearby galaxies using the 2MASS redshift survey and IceCube

The distribution of galaxies within the local universe is characterized by anisotropic features. Observatories searching for the production sites of astrophysical neutrinos can take advantage of these features to establish directional correlations between a neutrino dataset and overdensities in the galaxy distribution in the sky. The results of two correlation searches between a seven-year time-integrated neutrino dataset from the IceCube Neutrino Observatory, and the 2MASS Redshift Survey (2MRS) catalog are presented here. The first analysis searches for neutrinos produced via interactions between diffuse intergalactic Ultra-High Energy Cosmic Rays (UHECRs) and the matter contained within galaxies. The second analysis searches for low-luminosity sources within the local universe, which would produce subthreshold multiplets in the IceCube dataset that directionally correlate with galaxy distribution. No significant correlations were observed in either analyses. Herein, constraints are presented on the flux of neutrinos originating within the local universe through diffuse intergalactic UHECR interactions, as well as on the density of standard candle sources of neutrinos at low luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗

Fractionalized excitations probed by ultrasound

In this work, we study magnetoelastic interactions by means of ultrasound experiments in 𝛼−RuCl 3 —a prototypical material for the Kitaev spin model on the honeycomb lattice, with a possible spin-liquid state featuring Majorana fermions and ℤ 2 -flux excitations. We present results of the temperature and in-plane magnetic-field dependence of the sound velocity and sound attenuation for several longitudinal and transverse phonon modes propagating along high-symmetry crystallographic directions. A comprehensive data analysis above the magnetically ordered state provides strong evidence of phonon scattering by Majorana fermions. This scattering depends sensitively on the value of the phonon velocities relative to the characteristic velocity of the low-energy fermionic excitations describing the spin dynamics of the underlying Kitaev magnet. Moreover, our data displays a distinct reduction of anisotropy of the sound attenuation, consistent with the presence of thermally excited ℤ 2 visons. Here, we demonstrate the potential of phonon dynamics as a promising probe for uncovering fractionalized excitations in 𝛼−RuCl 3 and provide new insights into the 𝐻−𝑇 phase diagram of this material.

Kitaev magnet↗

A Field Enhancement Integration Design Featuring Misalignment Tolerance for Wireless EV Charging Using LCL Topology

This article proposes a magnetic integration design for EV wireless power transfer (WPT) systems, where the compensation and transmitting coils overlap one on top of the other to share the ferrite layer without any decoupling consideration. The magnetic field generated by both the compensation and transmitting coils are exploited to transfer power. To this end, a compensation method is proposed to enable magnetic field enhancement without any reactive power flow between the compensation and the transmitting coils, and to achieve input zero phase angle and constant current output. In addition, an efficient finite element analysis-based coil optimization algorithm is proposed to improve the coil misalignment tolerance on the horizontal plane, in which a reversely connected inner coil is used to stabilize the system output under misalignment conditions. Analytical and simulation results confirm transmission flux density enhancement and reduced leakage field characteristics of the proposed coil design. Finally, a scaled-down WPT prototype is built and tested to verify the performance and effectiveness of the design. The proposed design achieves 91.17% efficiency and misalignment tolerance up to 200 mm in any XY-direction while maintaining the power transfer and its efficiency. The results of this work provide insights into magnetic integration design of high-order compensation topologies featuring higher compactness, less ferrite usage, magnetic field enhancement, and misalignment tolerance for WPT systems.

42 ENGINEERING↗

Linking soil phosphorus with forest litterfall resistance and resilience to cyclone disturbance: A pantropical meta-analysis

While tropical cyclone regimes are shifting with climate change, the mechanisms underpinning the resistance (ability to withstand disturbance-induced change) and resilience (capacity to return to pre-disturbance reference) of tropical forest litterfall to cyclones remain largely unexplored pantropically. Single-site studies in Australia and Hawaii suggest that litterfall on low-phosphorus (P) soils is more resistant and less resilient to cyclones. Here, we conducted a meta-analysis to investigate the pantropical importance of total soil P in mediating forest litterfall resistance and resilience to 22 tropical cyclones. We evaluated cyclone-induced and post-cyclone litterfall mass (g/m 2 /day), and P and nitrogen (N) fluxes (mg/m 2 /day) and concentrations (mg/g), all indicators of ecosystem function and essential for nutrient cycling. Across 73 case studies in Australia, Guadeloupe, Hawaii, Mexico, Puerto Rico, and Taiwan, total litterfall mass flux increased from ~2.5 ± 0.3 to 22.5 ± 3 g/m 2 /day due to cyclones, with large variation among studies. Litterfall P and N fluxes post-cyclone represented ~5% and 10% of the average annual fluxes, respectively. Post-cyclone leaf litterfall N and P concentrations were 21.6 ± 1.2% and 58.6 ± 2.3% higher than pre-cyclone means. Mixed-effects models determined that soil P negatively moderated the pantropical litterfall resistance to cyclones, with a 100 mg P/kg increase in soil P corresponding to a 32% to 38% decrease in resistance. Based on 33% of the resistance case studies, total litterfall mass flux reached pre-disturbance levels within one-year post-disturbance. A GAMM indicated that soil P, gale wind duration and time post-cyclone jointly moderate the short-term resilience of total litterfall, with the nature of the relationship between resilience and soil P contingent on time and wind duration. Across pantropical forests observed to date, our results indicate that litterfall resistance and resilience in the face of intensifying cyclones will be partially determined by total soil P.

59 BASIC BIOLOGICAL SCIENCES↗

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗

Effect of Impurities on the Redox Properties of Goethite

Iron oxide minerals regulate the flux of electrons in the environment and are important hosts for trace and minor, yet critical, elements. Here, we present the first evidence of a direct link between the local coordination environments of Ni and Zn and the redox properties of their host phase goethite (α-FeOOH), the most abundant Fe(III) (oxyhydr)oxide at Earth’s surface. Here, we used aqueous redox measurements to show that the redox potential E H 0 , and hence the mineral’s stability, follows the order: pure goethite ≥ Zn-goethite > Ni-goethite. Parallel X-ray absorption and scattering measurements demonstrate, using quantum-informed analysis, that the local coordination environment of the smaller impurity, Ni, causes more bulk strain energy than Zn, which nearly accounts for the difference in E H 0 between Ni- and Zn-goethite. Our theory-informed, experimental study reveals how two common impurities affect the stability of goethite with implications for the biogeochemical reactivity of Fe(III) (oxyhydr)oxide in mediating elemental and electron fluxes in the environment.

42 ENGINEERING↗

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗

Material transfers detection with seismic observations

We are exploring the use of data from a seismic network around a research nuclear reactor and isotope production facility at Oak Ridge National Laboratory to study activity patterns related to the transfer of nuclear material. A sensor network with eight seismometers was installed around the High Flux Isotope Reactor and the Radiochemical Engineering Development Center and started operating in October 2019. These data are used to detect and characterize operational events around facility. We are using those data to extract signals related to the movement of vehicles involved in the transport of nuclear materials (e.g., transportation of reactor fuel, targets, and product isotopes). Large vehicles produce mechanical energy that can be observed as seismic signals, either as the result of sound emanating from the vehicle engine or by the generation of surface waves as a response of the ground to the load of the vehicles. In particular, the station located near the entrance gate of the facility displays a clear seismic signature as vehicles cross a metal platform on the ground. This signal is characterized by sharp energy bursts that correspond to the number of axles on the vehicle. We developed and are testing an algorithm to identify sequences of energy bursts to count and characterize vehicles of different sizes. We are also exploring the use of seismic polarization analysis to measure the degree of polarization of the observed signals. These two techniques for detection and characterization of vehicles used for nuclear material transfers will be validated with seismic data from a targeted collection using different vehicles and routes similar to the ones used in real-life scenarios at the High Flux Isotope Reactor.

Marcillo, Omar↗

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗

Influences of Environmental Relative Humidity and Horizontal Scale of Subcloud Ascent on Deep Convective Initiation

This study examines two factors impacting initiation of moist deep convection: free-tropospheric environmental relative humidity (φ E ) and horizontal scale of subcloud ascent (R sub ), the latter exerting a dominant control on cumulus cloud width. A simple theoretical model is used to formulate a “scale selection” hypothesis: that a minimum R sub is required for moist convection to go deep, and that this minimum scale decreases with increasing φ E . Specifically, the ratio of R 2 sub to saturation deficit (1 - φ E ) must exceed a certain threshold value that depends on cloud-layer environmental lapse rate. Idealized, large-eddy simulations of moist convection forced by horizontally varying surface fluxes show strong sensitivity of maximum cumulus height to both φ E and R sub consistent with the hypothesis. Increasing R sub by only 300–400 m can lead to a large increase (>5 km) in cloud height. A passive tracer analysis shows that the bulk fractional entrainment rate decreases rapidly with R sub but depends little on φ E . However, buoyancy dilution increases as either R sub or φ E decreases; buoyancy above the level of free convection is rapidly depleted in dry environments when R sub is small. While deep convective initiation occurs with an increase in relative humidity of the near environment from moistening by earlier convection, the importance of this moisture preconditioning is inconclusive as it is accompanied by an increase in R sub . Overall, it is concluded that small changes to R sub driven by external forcing or by convection itself could be a dominant regulator of deep convective initiation.

54 ENVIRONMENTAL SCIENCES↗

Modeling and Manipulating Phenylpropanoid Pathway Flux for Bioenergy

We have made important strides toward our project goals over the past five years. First, we developed and published a novel isotopic feeding system to investigate metabolic fluxes of the lignin biosynthesis. Next, we have developed and published a kinetic model of the phenylpropanoid pathway. As an extension of that work, we developed a multilayer kinetic model of lignin biosynthesis. Next we conducted a systematic transcriptome and metabolome analysis for phenylpropanoid metabolism. Finally, we have identified three strategies that enhance the accumulation of 2-phenylethanol in plants.

09 BIOMASS FUELS↗

Second Target Station High-Fidelity Target Activation Comparison

The development of the Second Target Station (STS) target system at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is well underway. The target system at STS consists of a rotating target disk that contains 21 segments of tungsten clad in tantalum clad in steel. A key aspect of the design of the target system is to account for the delayed heating and material damage caused by the delayed dose from decaying radionuclides. These radionuclides are a product of either spallation reactions or transmutation of the nuclei in the target system. These radionuclides build up in the target system components over the lifetime of the facility, and the radiation that is emitted can deposit energy in the components causing significant component heating and material damage. Monte Carlo N-Particle (MCNP) Version 6.2 transports the various particle species and calculates the spallation products and neutron fluxes throughout the target system. These spallation products and neutron fluxes along with the material definition of each component are relayed to the CINDER2008 transmutation code to calculate the radionuclide inventories and the corresponding decay gamma emission spectra. MCNP6.2 coupled with CINDER2008 is the computational method-of-choice for the analysis discussed in the following sections of this report. The analysis focuses on validating major assumptions in calculating the radionuclide inventory in the STS target system: all of the target segments are fresh, unirradiated material when the protons are incident on the segment, the average of the 21 segments of the target is sufficient to represent a single segment, and that averaging the proton pulse structure over time does not significantly affect the radionuclide inventory. The position-averaged, high-fidelity, and single-tally computational methods are used to validate the assumptions and provide a point of comparison to evaluate how the assumptions impact the radionuclide inventories. A more detailed explanation of the three computational methods is provided in Section 2. The position-averaged and single-tally methods are less computationally expensive when compared with the high-fidelity method where 54,000 individual calculations are needed to calculate 1 hr of STS operation. Section 3 details the comparison of the three methods to show that the assumptions made in the position-averaged method do not significantly impact the radionuclide inventory after 1 hr of operation. The discussions and results in this report are for 1 hr of operation. Due to the computational cost associated with calculating the transmutation and activation using the high-fidelity method, only 1 hr of operation has been calculated. The discrepancies observed after 1 hr of operation are not extrapolated out to longer operational times, and this report does not address how the discrepancies between the computational methods may manifest for longer operational periods.

43 PARTICLE ACCELERATORS↗

King City Asbestos Corporation (KCAC) Mine Carbon Mineralization Field Test (Final Combined Report)

The overall objective of the KCAC carbon mineralization project is to determine the most successful approaches to CO 2 mineralization of asbestos-bearing serpentinite rocks and mine waste at the closed KCAC asbestos mine. The objectives of Budget Period 1 (BP1) are to characterize the site including, mineralogical and geochemical variability in site materials, baseline conditions including CO 2 flux and local climate, and to establish the designs for the second, experimental phase of the project. As described in Section 3 of this report, all tasks outlined in the Statement of Project Objective (SOPO) have been completed or are on schedule to be completed by the end of BP1. As described in Section 4, all success criteria and milestones have been met or will be met by the end of BP1. We conclude that the findings of BP1 strongly support moving to the second phase of the project to investigate methods to increase the rate of CO 2 mineralization. CO 2 flux monitoring during BP1 established an estimated baseline uptake of CO 2 over the 142,193 m 2 KCAC mine site of 147 and 179 metric tons CO 2 /yr using dynamic closed chamber and eddy covariance monitoring methods, respectively. Mineralogical analysis shows that the loose surficial waste material contains pyroaurite and minor hydromagnesite, indicating past carbon mineralization. Sediment samples at depths below 60 cm contained ~3 wt% brucite and localized high brucite contents of ~24 wt% were measured at depth, indicating significant unreacted material at relatively shallow depths. Together, these results demonstrate that despite prolonged exposure to the atmosphere, reactive material is still present and accessible and is naturally taking up atmospheric CO 2 . However, these natural rates are still well below rates that can be achieved by applying practices that can accelerate the CO 2 mineralization reaction. Designs for the Field-Testing Phase (BP2) of the project have been established, incorporating the findings of BP1. Based on results to date we have selected two methods to accelerate the rate of CO 2 mineralization, which we refer to as the Tilling Method and the Greenhouse Method. We hypothesize that these methods can increase the CO 2 mineralization rate by up to 5 times greater than the natural background rate measured in BP1. An Initial Lifecyle Assessment (LCA) and techno-economic analysis (TEA) were conducted for these two methods. The LCA concluded that the baseline version of the Greenhouse and Tilling methods at commercial scale are expected to achieve a CDR efficiency of 97% and 71%, respectively (cumulative, over 30 years of operation). Incorporating the net CO 2 removal rates calculated in the LCA, the TEA found the levelized, 30-year costs of the baseline Greenhouse and Tilling methods to be $\$$95/t-CO 2 and $\$$119/t-CO 2 , respectively, and that costs could be brought to well below $\$$100/t-CO 2 for either method by scaling up the volume of rock utilized. An initial EH&S Risk Assessment, evaluating the potential EH&S risks of the proposed carbon mineralization methods at commercial scale, was completed by LLNL EH&S professionals for these two methods.

58 GEOSCIENCES↗