Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global sensitivity 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The Determinants of Offshore Wind's Role in a Future U.S. Energy System: A Preliminary Modeling Sensitivity Analysis [Slides]

Offshore wind is an emerging industry in the with rapidly changing technology advancements and an accelerating global deployment. Policy demands and interest in the United States has grown, but the long-term future for offshore wind remains uncertain. This study is designed to evaluate offshore wind's potential role in the future U.S. energy system. The analysis uses long-term power system models to assess a wide-range of power system possibilities. Specific objectives of the analysis are to: (1) identify the determinants of offshore wind deployment under different scenarios and (2) understand the impacts of substantial levels of offshore wind deployment on the power system. We also highlight future research needs for more robust understanding the opportunities and challenges to U.S. offshore wind development.

17 WIND ENERGY↗
DOE OSTI2022

Global dipole symmetry, compact Lifshitz theory, tensor gauge theory, and fractons

We study field theories with global dipole symmetries and gauge dipole symmetries. The famous Lifshitz theory is an example of a theory with a global dipole symmetry. We study in detail its 1+1d version with a compact field. When this global symmetry is promoted to a $U(1)$ dipole gauge symmetry, the corresponding gauge field is a tensor gauge field. This theory is known to lead to fractons. In order to resolve various subtleties in the precise meaning of these global or gauge symmetries, we place these 1+1d theories on a lattice and then take the continuum limit. Interestingly, the continuum limit is not unique. Different limits lead to different continuum theories, whose operators, defects, global symmetries, etc. are different. We also consider a lattice gauge theory with a $\mathbb Z_N$ dipole gauge group. Surprisingly, several physical observables, such as the ground state degeneracy and the mobility of defects depend sensitively on the number of sites in the lattice. Our analysis forces us to think carefully about global symmetries that do not act on the standard Hilbert space of the theory, but only on the Hilbert space in the presence of defects. We refer to them as time-like global symmetries and discuss them in detail. Furthermore, these time-like global symmetries allow us to phrase the mobility restrictions of defects (including those of fractons) as a consequence of a global symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Understanding PDF uncertainty in W boson mass measurements*

Abstract We study the dependence of the transverse mass distribution of charged leptons and the missing energy on parton distributions (PDFs) adapted to W boson mass measurements at the CDF and ATLAS experiments. We compare the shape variations of the distribution induced by different PDFs and find that the spread of predictions from different PDF sets can be significantly larger than the PDF uncertainty predicted by a specific PDF set. We suggest analyzing the experimental data using up-to-date PDFs to gain a better understanding of the PDF uncertainties in W boson mass measurements. We also perform a series of Lagrange multiplier scans to identify the constraints on the transverse mass distribution imposed by individual data sets in the CT18 global analysis. In the case of the CDF measurement, the distribution is mostly sensitive to d -quark PDFs in the intermediate x region, which are largely constrained by DIS and Drell-Yan data on deuteron targets and Tevatron lepton charge asymmetry data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gluon helicity from global analysis of experimental data and lattice QCD Ioffe time distributions

We perform a new global analysis of spin-dependent parton distribution functions with the inclusion of Ioffe time pseudodistributions computed in lattice QCD (LQCD), which are directly sensitive to the gluon helicity distribution, Δ g . These lattice data have an analogous relationship to parton distributions as do experimental cross sections, and can be readily included in global analyses. We focus in particular on the constraining capability of current LQCD data on the sign of Δ g at intermediate parton momentum fractions x , which was recently brought into question by analysis of data in the absence of parton positivity constraints. We find that present LQCD data cannot discriminate between positive and negative Δ g solutions, although significant changes in the solutions for both the gluon and quark sectors are observed. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Techno-Economic and life cycle assessment of standalone Single-Stream material recovery facilities in the United states

Material Recovery Facilities (MRFs) are crucial players in achieving a circular economy. MRFs receive complex waste streams and separate valuable recyclables from these mixtures. This study conducts techno-economic analysis (TEA) to estimate the net present value (NPV) and life cycle assessment (LCA) to estimate different environmental impacts of a commercial scale standalone, single-stream MRF to assess the economic feasibility and environmental impacts of recovering valuable recyclables from an MRF processing 120,000 tonnes per year (t/y). The TEA employs a discounted cash flow rate of return (DCFROR) analysis over a 20-year facility lifetime, along with a sensitivity analysis on the impact of different operating and economic parameters. Results show that the total fixed cost of building the MRF facility is $\$ $23 MM, and the operating cost is $\$ $45.48/tonne. Further the NPV of the MRF can vary from $\$ $3.57 MM to $\$ $60 MM, while 100-year global warming potential can range from 5.98 to 8.53 kg carbon dioxide equivalents (CO2-eq) per tonne of MSW. We have also found that MSW composition (arising from regional effects) significantly impacts costs, 100-year global warming potential, and other impact categories such as acidification potential, eutrophication potential, ecotoxicity, ozone depletion, photochemical oxidation, carcinogenic effects, and non-carcinogenic effects. Sensitivity and uncertainty analysis indicate that waste composition and market prices significantly impact the profitability of the MRF, and the waste composition mostly impacts global warming potential. Our analysis also indicates that facility capacity, fixed capital cost, and waste tipping fees are vital parameters that affect the economic viability of MRF operations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vegetation-induced asymmetric diurnal land surface temperatures changes across global climate zones

Unprecedented global vegetation greening during past decades is well known to affect annual and seasonal land surface temperatures (LST). However, the impact of observed vegetation cover change on diurnal LST across global climatic zones is not well understood. In this study, using global climatic time-series datasets, we investigated the long-term growing season daytime and nighttime LST changes globally and explored associated dominant contributors including vegetation and climate factors including air temperature, precipitation, and solar radiation. Results revealed asymmetric growing season mean daytime and nighttime LST warming (0.16 °C/10a and 0.30 °C/10a, respectively) globally from 2003 to 2020, as a result, the diurnal LST range (DLSTR) declined at 0.14 °C/10a. The sensitivity analysis indicated the LST response to changes in LAI, precipitation, and SSRD mainly concentrated during daytime instead of nighttime, however, which showed comparable sensitivities for air temperature. Combining the sensitivities results and the observed LAI and climate trends, we found rising air temperature contributes to 0.24 ± 0.11 °C/10a global daytime LST warming and 0.16 ± 0.07 °C/10a nighttime LST warming, turns to be the dominant contributor to the LST changes. Increased LAI cooled global daytime LST (–0.068 ± 0.096 °C/10a) while warmed nighttime LST (0.064 ± 0.046 °C/10a); hence LAI dominates declines in DLSTR trends (–0.12 ± 0.08 °C/10a), despite some daynight process variations across climate zones. In Boreal regions, reduced DLSTR was due to nighttime warming from LAI increases. In other climatic zones, daytime cooling, and DLSTR decline, was induced by increased LAI. Biophysically, the pathway from air temperature heats the surface through sensible heat and increased downward longwave radiation during day and night, while the pathway from LAI cools the surface by enhancing energy redistribution into latent heat rather than sensible heat during the daytime. These empirical findings of diverse asymmetric responses could help calibrate and improve biophysical models of diurnal surface temperature feedback in response to vegetation cover changes in different climate zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nano-Engineered Catalyst Supported on Ceramic Hollow Fibers for the Utilization of CO 2 in Dry Reforming to Produce Syngas

The objective of this project was to develop a novel catalytic reactor containing nano-engineered catalysts for the utilization of CO 2 (captured from coal-fired power plants and other CO 2 emitting sources) in dry methane reforming (DMR) (CO 2 + CH 4 → 2 H 2 + 2 CO) to produce synthesis gas (syngas). The technology aims to reduce CO2 emissions by developing beneficial uses for CO 2 from coal-fired power plants. It also offers an alternative to mitigate CO 2 emissions in areas where geologic storage may not be an optimal solution and/or utilization could significantly offset the costs of carbon capture and sequestration. The nano-engineered Ni-based catalyst was prepared by atomic layer deposition (ALD). The Ni particles were as small as ~2-4 nm. The nano-engineered catalyst showed CH 4 conversion >95%, H 2 /CO ratio in the range of 0.7-1.0, and CH 4 reforming rate as high as 2,500 L/h/gNi at 850 ºC and pressure of 15-25 psia. The Ni-based ALD catalyst also showed good stability in DMR reaction during a 200-h continuous operation at 850 °C. This is due to strong bonding between the nanoparticles and substrates since the Ni nanoparticles were chemically bonded to the substrate during the ALD process. The high thermal stability maintains the high dispersion of Ni nanoparticles, which can inhibit coke formation because their step edges are small enough to limit carbon nucleation and growth. Technoeconomic analysis (TEA) indicates the levelized cost of syngas (LCOS) is $172/ton with our technology, which is lower than the equivalent (molar) cost of hydrogen produced by steam methane reforming (SMR) or autothermal reforming (ATR). The major operating cost is natural gas feed and fuel, and the levelized cost is highly sensitive to the price of natural gas and relatively insensitive to the CAPEX. Revenues from syngas could have a significant impact on the net cost of electricity (COE), depending on the cost of natural gas and the selling price of syngas, estimated at $36 per MWh if the syngas were sold at $195 per ton. Following DOE NETL’s guidance, a lifecycle analysis (LCA) was conducted to compare with SMR. The functional unit for the basis of comparison was defined as 1kg carbon monoxide in the product stream. The global warming potential (GWP) of our process was found to be 40% lower than the state-of-the-art SMR process. The sensitivity analysis confirms the emissions are most sensitive to the natural gas fuel requirements to deliver heat to the process.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Current and Future Global Lake Methane Emissions: A Process‐Based Modeling Analysis

Abstract Freshwater ecosystem contributions to the global methane budget remains the most uncertain among natural sources. With warming and accompanying carbon release from thawed permafrost and thermokarst lake expansion, the increase of methane emissions could be large. However, the impact and relative importance of various factors related to warming remain uncertain. Based on diverse lake characteristics incorporated in modeling and observational data, we calibrate and verify a lake biogeochemistry model. The model is then applied to estimate global lake methane emissions and examine the impacts of temperature increase for the first and the last decades of the 21st century under different climate scenarios. We find that current emissions are 24.0 ± 8.4 Tg CH 4 yr −1 from lakes larger than 0.1 km 2 , accounting for 11% of the global total natural source as estimated based on atmospheric inversion. Future projections under the RCP8.5 scenario suggest a 58%–86% growth in emissions from lakes. Our model sensitivity analysis indicates that additional carbon substrates from thawing permafrost may enhance methane production under warming in the Arctic. Warming enhanced methane oxidation in lake water can be an effective sink to reduce the net release from global lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ensemble modeling of global lake evaporation under climate change

Global projections of lake evaporation are typically based on simulations using single mechanistic models. However, because of its complex interactions with various lake physical properties, environmental and anthropogenic drivers, lake evaporation is highly variable and sensitive to the choice of model used. In this study, we present a multi-model analysis to investigate differences across global simulations of lake evaporation during the warm-season using three different lake models driven by outputs from four general circulation models (GCM) (i.e. 12 model combinations in total) for historic and future scenarios. Our results suggest substantial differences among lake-climate model simulations of lake evaporation. These differences varied throughout the 20th and 21st century, with model driver data explaining 74% of the variance in future projections of warm- season lake evaporation. Our projections indicate that, by the end of the 21st century (2070–2099), global annual lake evaporation rates will increase by 10–27% under Representative Concentration Pathways (RCPs) 2.6–8.5. Finally, we highlight the importance of using a multi-model approach for the prediction of future global lake evaporation responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Type II Hamiltonian Variational Principle and Adjoint Systems for Lie Groups

We present a novel Type II variational principle on the cotangent bundle of a Lie group which enforces Type II boundary conditions, i.e., fixed initial position and final momentum. In general, such Type II variational principles are only globally defined on vector spaces or locally defined on general manifolds; however, by left translation, we are able to define this variational principle globally on cotangent bundles of Lie groups. Type II boundary conditions are particularly important for adjoint sensitivity analysis, which is our motivating application. As such, we additionally discuss adjoint systems on Lie groups, their properties, and how they can be used to solve optimization problems subject to dynamics on Lie groups.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sensitivity analysis for characterizing the impact of HNGD model on the prediction of hydrogen redistribution in Zircaloy cladding using BISON code

Hydrogen in zirconium cladding is able to precipitate into zirconium hydrides which impacts cladding integrity. The Hydride Nucleation-Growth-Dissolution (HNGD) model in the BISON code accounts for the precipitation and dissolution kinetics of hydride in Zircaloy material. This paper presents global sensitivity analyses of the HNGD model aiming to enhance our understanding of the hydride precipitation phenomena by quantifying the variance that key parameters have on the prediction of hydrogen behavior under various environmental conditions. Here, model predictions are compared to experimental data obtained under two different conditions: 1) with uniformly precharged specimens subjected to a linear thermal gradient, and 2) specimens precharged with a cathodically applied hydride rim at one end of the sample and subjected to an asymmetric thermal gradient. The Sobol sensitivity analysis identifies the key parameters in the HNGD model for both types of specimens. For linear temperature cases, the heat of transport dominates the accuracy of predictions when no precipitation occurs at the cold end, while Terminal Solid Solubility for Dissolution (TSSD) is the most important parameter when precipitation occurs. A large variation in the predicted hydrogen concentration profiles is found in the range of high TSSD due to the occurrence of precipitation. For asymmetric temperature cases, the solubility coefficient gives the largest impact on the predicted hydrogen distribution, as it determines the amount of solute hydrogen dissolved from the initially applied hydride rim. A large discrepancy in hydrogen distribution between simulations and experiments exists with the asymmetric specimens because BISON simulations fail to predict the precipitation of hydride at the cooler end. Comparative studies using former and updated models verifies the significant impact of the hydride growth mechanism on predicted hydrogen concentration profiles. In particular, when hydride initially exists, changes in TSSD generate a large variation in the predicted amount of precipitation by hydride growth, giving large uncertainty in predicting the hydrogen distribution over the sample length. The outputs characterize the significant impact of the hydride growth mechanism in the HNGD model on predicting hydrogen behavior, and improve the understanding of the precipitation of hydride in Zircaloy cladding within a range of expected environmental conditions. The analyses indicate work is still needed to improve the hydride solvus models in the BISON code to accurately predict experimentally observed hydride concentrations and distributions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗
DOE OSTI2026

Jet charge with global event shapes: Probing quark flavor dynamics

We propose measuring the jet electric charge of jet regions, defined within the framework of global event shapes, as a probe of quark flavor dynamics within the nucleon and the hadronization process. In particular, we consider a measurement of the jet region charge while simultaneously keeping track of the energy flow throughout the event, as characterized by the global event shape. As a concrete example, we focus on the measurement of the (Q), the jet charge of the jet region (𝐽) defined within the framework of the 1-Jettiness global event shape (𝜏 1 ) for the Deep Inelastic Scattering (DIS) process, 𝑒 − +𝑝→𝑒 − +𝐽+𝑋, with unpolarized or longitudinally polarized protons. The 1-Jettiness distribution, binned according to jet charge, allows for enhanced quark flavor separation of the initial state unpolarized or polarized PDFs. On the other hand, the jet charge distribution binned by 1-Jettiness can serve as a probe of quark flavor dynamics in the final state hadronization process. We derive a factorization theorem for simultaneous measurements of 𝜏 1 and Q in the resummation region, 𝜏 1 ≪𝑃 𝐽 𝑇 , where 𝑃 𝐽 𝑇 denotes the transverse momentum of the jet region. The correlation between the jet region charge and the charge of the struck quark is expected to be the strongest in the resummation region. The factorization theorem contains a new universal , generalizing the standard jet function to include a jet charge measurement. Therefore these universal functions can be extracted from a global analysis of N-jettiness and thrust at 𝑒 + ⁢𝑒 − colliders. We provide simulation studies to demonstrate the sensitivity of the 1-Jettiness jet charge observable to quark flavor dynamics in nucleon structure and explore the possibility of probing the final state hadronization process. This observable is well-suited for applications with existing HERA data and the future Electron-Ion Collider (EIC).

Mantry, Sonny [University of North Georgia, Dahlon↗
DOE OSTI2021

Engineering well-expressed, V2-immunofocusing HIV-1 envelope glycoprotein membrane trimers for use in heterologous prime-boost vaccine regimens

HIV-1 vaccine immunofocusing strategies may be able to induce broadly-reactive neutralizing antibodies (NAbs). Here, we engineered a panel of diverse, membrane-resident native HIV-1 trimers vulnerable to two broad targets—the V2 apex and fusion peptide (FP). Selection criteria included i) high expression and ii) infectious function, so that trimer neutralization sensitivity can be profiled in pseudovirus (PV) assays. Initially, we boosted the expression of 17 candidate trimers by truncating gp41 and introducing a gp120-gp41 SOS disulfide to prevent gp120 shedding. "Repairs" were made to fill glycan holes and eliminate other strain-specific aberrations. A new neutralization assay allowed PV infection when our standard assay was insufficient. Trimers with exposed V3 loops, a target of non-NAbs, were discarded. To try to increase V2-sensitivity, we removed clashing glycans and modified the C-strand. Notably, a D167N mutation improved V2-sensitivity in several cases. Glycopeptide analysis of JR-FL trimers revealed near complete sequon occupation and that filling the N197 glycan hole was well-tolerated. In contrast, sequon optimization and inserting/removing glycans at other positions frequently had global "ripple" effects on glycan maturation and sequon occupation throughout the gp120 outer domain and gp41. V2 MAb CH01 selectively bound to trimers with small high mannose glycans near the base of the V1 loop, thereby avoiding clashes. Knocking in a rare N49 glycan was found to perturb gp41 glycans, increasing FP NAb sensitivity—and sometimes improving expression. Finally, a biophysical analysis of VLPs revealed that i) ~25% of particles bear Env spikes, ii) spontaneous particle budding is high and only increases 4-fold upon Gag transfection, and iii) Env+ particles express ~30–40 spikes. Taken together, we identified 7 diverse trimers with a range of sensitivities to two targets to allow rigorous testing of immunofocusing vaccine concepts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Weather Sensitive High Spatio-Temporal Resolution Transportation Electric Load Profiles For Multiple Decarbonization Pathways

Electrification of transport compounded with climate change will transform hourly load profiles and their response to weather. We present a novel approach to generating hourly electric load profiles that considers charging strategies and evolving sensitivity to temperature. The approach consists of downscaling annual state-scale sectoral load projections from the multisectoral Global Change Analysis Model (GCAM) into hourly electric load profiles leveraging high resolution climate and population datasets. Profiles are developed and evaluated at the Balancing Authority scale, with a 5-year increment until 2050 over the Western U.S. Interconnect for multiple decarbonization pathways and climate scenarios. The datasets are readily available for production cost model analysis. Our open source approach is transferable to other regions.

Decarbonization, tranportation, Electric vehicle c↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.

58 GEOSCIENCES↗