Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Predictive Control”

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,297 records · Page 72

GrainNN: A neighbor-aware long short-term memory network for predicting microstructure evolution during polycrystalline grain formation

High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.

36 MATERIALS SCIENCE↗

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗

Predictive Risk Mapping of Schistosomiasis in Madagascar Using Ecological Niche Modeling and Precision Mapping

Schistosomiasis is a neglected tropical disease (NTD) found throughout tropical and subtropical Africa. In Madagascar, the condition is widespread and endemic in 74% of all administrative districts in the country. Despite the significant burden of the disease, high-resolution risk maps have yet to be produced to guide national control programs. This study used an ecological niche modeling (ENM) and precision mapping approach to estimate environmental suitability and disease transmission risk. The results show that suitability for schistosomiasis is widespread and covers 264,781 km 2 (102,232 sq miles). Covariates of significance to the model were the accessibility to cities, distance to water, enhanced vegetation index (EVI), annual mean temperature, land surface temperature (LST), clay content, and annual precipitation. Disease transmission risk is greatest in the central highlands, tropical east coast, arid-southwest, and northwest. An estimated 14.9 million people could be at risk of schistosomiasis; 11.4 million reside in rural areas, while 3.5 million are in urban areas. This study provides valuable insight into the geography of schistosomiasis in Madagascar and its potential risk to human populations. Because of the focal nature of the disease, these maps can inform national surveillance programs while improving understanding of areas in need of medical interventions.

60 APPLIED LIFE SCIENCES↗

A Roadmap for the Future of Systems Biology in Cancer Research

Cancer systems biology seeks to understand how cancer arises as a system of interconnected molecules, cells, and tissues, with the goal of understanding, predicting, and controlling the disease. In the last decade, the field has rapidly grown as advances in experimental, computational, and analytic technologies have improved our ability to capture and recapitulate the complexities of cancer at multiple scales. However, the field’s promise to understand how specific molecular changes give rise to altered cancer outcomes remains incompletely fulfilled. Fortunately, an opportunity exists to accelerate progress by better coordinating modeling and data-gathering efforts across the cancer systems biology community. This will create the foundation for building accurate, multiscale cancer models that can better predict and identify improved therapeutic interventions. Here, in this study, we outline some of the current challenges in cancer systems biology research, how they can be addressed, and actions that the community can take to accelerate progress in the field.

Modeling & Simulation↗

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗

Harnessing dimethyl ether with ultra-low-grade heat for scaling-resistant brine concentration and fractional crystallization

Solvent-driven separations may enable scalable concentration of hypersaline brines, supporting a circular resource economy from the extraction of lithium and rare earth elements from spent battery and magnet leachates. This work analyses a novel solvent-driven water extraction (SDWE) system employing dimethyl ether (DME) and ultra-low-grade heat for brine concentration and fractional crystallization. SDWE exploits DME’s unique properties: (1) a low dielectric constant that promotes water solubility over charged solutes by a factor of 10 3 , and (2) a high volatility that facilitate efficient DME reconcentration with ultra-low-grade heat. The techno-economic viability of SDWE is assessed with a computational framework that encompasses a liquid–liquid separator and a solvent concentrator. We integrate the extended universal quasichemical model with the virial equation of state to predict the compositions of the complex three-phase DME-water mixture at vapor–liquid and liquid–liquid equilibrium. Subsequently, we optimize the thermodynamic and economic performance of SDWE, by controlling the interstage flash pressure, heat source temperature, and the number of concentrating stages. DME-based SDWE concentrates an input saline feed to 5.5 M and regenerates over 99 % of the DME using ultra-low-grade heat below 50 °C, with a DME/water selectivity ratio of 125. Here our calculations reveal that optimal performance is achieved at interstage flash pressures of 0.4 – 0.5 bar for heat source temperatures between 323–373 K, with improved exergetic efficiencies at lower temperatures. At a heat source temperature of 323 K and an interstage pressure of 0.489 bar, DME-driven SDWE achieves an optimal thermodynamic efficiency of 20.5 % and a projected specific cost of US$ 1.93 m -3 . These specific costs suggest that SDWE is competitive with commercialized thermal distillation technologies, while mitigating the traditional risks associated with scaling in heat and mass exchangers with hypersaline brines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Air compressor load forecasting using artificial neural network

Air compressor systems are responsible for approximately 10% of the electricity consumed in United States and European Union industry. As many researches have proven the effectiveness of using Artificial Neural Network inair compressor performance prediction, there is still a need to forecast the air compressor electrical load profile. The objective of this study is to predict compressed air systems’ electrical load profile, which is valuable to industry practitioners as well as software providers in developing better practice and tools for load managementand look-ahead scheduling programs. Two artificial neural networks, Two-Layer Feed-Forward Neural Networkand Long Short-Term Memory were used to predict an air compressors electrical load. Compressors with three different control mechanisms are evaluated with a total number of 11,874 observations. Here, the forecasts were validated using out-of-sample datasets with 5-fold cross-validation. Models produced average coefficient ofdetermination values from 0.24 to 0.94, average root-mean-square errors from 0.05 kW - 5.83 kW, and meanabsolute scaled errors from 0.20 to 1.33. The results indicate that both artificial neural networks yield goodresults for compressors using variable speed drive (average R 2 = 0.8 and no naïve forecasting), only the longshort-term memory model gives acceptable results for compressors using on/off control (average R 2 = 0.82 and no naïve forecasting), and no satisfactory results are obtained for load/unload type air compressors (models constituting active forecasting).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arctic shrub expansion generates regional variation in litter decomposition by altering litter quality and the decomposition environment

Abstract The expansion of deciduous shrubs into the graminoid‐dominated arctic tundra is expected to alter litter decomposition by changing litter quality and local abiotic and biotic conditions. However, it remains unclear how shrub expansion will affect litter decomposition at regional scales, where macroclimate is expected to be the dominant regulator of decomposition. To determine the relative influence of macroclimate and local controls on regional patterns of litter mass loss and nitrogen release, we conducted two hierarchical litter decomposition experiments across spatial scales. We decomposed leaf and root litter from a prominent graminoid ( Eriophorum vaginatum ) and three genera of deciduous shrubs ( Betula , Alnus and Salix ) for 1 year within replicated plots at five sites spanning a 160 km latitudinal gradient in northern Alaska. Using Eriophorum litter as a substrate, we found that macroclimate was the primary regulator of mass loss but had opposing effects on leaf and root litter. As summer temperature increased along the latitudinal gradient (11.9 to 13.9°C), leaf litter mass loss increased by 20% whereas root litter mass loss decreased by 33%. Leaf nitrogen release also increased with summer temperature. Conversely, root nitrogen release was controlled by the vegetation type of the decomposition environment. Using different shrub litters as substrates, we found that litter quality and its interaction with soil microclimate and macroclimate controlled decomposition. Overall, shrub root litters decomposed faster than Eriophorum root litter, losing 53% more mass and 190% more nitrogen across all sites and decomposition environments. For leaf litter, however, patterns varied by litter genus, with Salix losing more mass and Betula and Alnus losing less mass than Eriophorum . Our findings demonstrate that shrub expansion in the Arctic can regulate leaf and root litter decomposition at the regional scale through its effects on local controls, primarily litter quality. Ongoing increases in shrub cover are likely to accelerate the turnover of root litter carbon and nitrogen pools in tundra ecosystems. Therefore, including shrub‐related processes in Earth system models will improve our ability to predict regional‐scale litter decomposition and its effects on carbon and nutrient cycling in a warming Arctic. Read the free Plain Language Summary for this article on the Journal blog.

Vozzo, Justin T. [Department of Natural Resources ↗

The effect of blood viscosity on shear‐induced hemolysis using a magnetically levitated shearing device

Abstract Introduction Blood contacting medical devices, including rotary blood pumps, can cause shear‐induced blood damage that may lead to adverse effects in patients. Due in part to an inadequate understanding of how cell‐scale fluid mechanics impact red blood cell membrane deformation and damage, there is currently not a uniformly accepted engineering model for predicting blood damage caused by complex flow fields within ventricular assist devices (VADs). Methods We empirically investigated hemolysis in a magnetically levitated axial Couette flow device typical of a rotary VAD. The device is able to accurately control the shear rate and exposure time experienced by blood and to minimize the effects of other uncharacterized stresses. Using this device, we explored the effects of both hematocrit and plasma viscosity on shear‐induced hemolysis to characterize blood damage based on the viscosity‐independent shear rate, rather than on shear stress. Results Over a shear rate range of 20 000 – 80 000 1/s, the Index of Hemolysis (IH) was found to be dependent upon and well‐predicted by the shear rate alone. IH was independent of hematocrit, bulk viscosity, or the suspension media viscosity and less correlated to shear stress (MSE = 0.46–0.75) than to shear rate (MSE = 0.06–0.09). Conclusion This study recommends that future investigations of shear‐induced blood damage report findings with respect to the viscosity‐neutral term of shear rate, in addition to the bulk whole blood viscosity measured at an appropriate shear rate relevant to the flow conditions of the device.

Krisher, James A.↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transition Core Modeling for Extended Enrichment & Accident-Tolerant Fuels Using Polaris/PARCS

Commercial light water reactor (LWR) operators and fuel vendors are currently interested in increasing the low-enriched uranium (LEU) fuel enrichments from the current limit of 5.0 $^w/_o$ $^{235}U$ up to 10 $^ w/_o$ $^{235}U$ (referred to as "LEU+") in their current fleets; they are also interested in using accident-tolerant fuel (ATF) with both LEU and LEU+ fuel. This report aims to identify modeling challenges and accuracy concerns in transition core analysis using the SCALE Polaris lattice physics code and U.S. Nuclear Regulatory Commission core simulator PARCS. At the time this study was started, no publicly available LEU+ core designs existed for boiling water reactor (BWR) or pressurized water reactor (PWR) systems. Therefore, fuel lattices were shuffled within a multi-assembly model to mimic neutronically challenging lattice combinations seen in transition cores, such as a fresh LEU+ lattice next to depleted LEU lattices. In addition to multi-assembly models, whole-core BWR transition core calculations were performed for ATF and LEU+ fuel using an existing Hatch-1 Cycle 3 core model. A whole-core BWR model was chosen due to the more heterogeneous core designs compared to those for a PWR core. Since the original core is an old checkerboard core design and no core or fuel design optimization was performed for the modeled fuel types, these core calculations were intended only to provide: (1) Comparisons of core characteristics of interest, such as the pin power distributions and peaking factors, Doppler temperature coefficients (DTCs), and control blade worths (CBWs) under challenging core designs, (2) Identification of reactor physics challenges in modeling LEU+ and ATF cores, and (3) A stress test for the Polaris/PARCS two-step modeling approach, including characterization of the relative accuracy for predicting characteristics of interest such as pin power distributions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Continental-Scale Controls on Hyporheic Respiration Revealed by Knowledge-Guided Machine Learning

Hyporheic zone sediments regulate organic matter turnover and in-stream respiration, yet controls on sediment respiration remain poorly constrained across heterogeneous river networks, limiting prediction of stream metabolism and carbon processing at continental scales. Here, we integrate observations from ~90 river corridors across the United States in the WHONDRS consortium with a knowledge-guided machine learning (KGML) framework that couples thermodynamic rate theory with machine learning to identify dominant controls on hyporheic respiration. Diagnostic analyses show that organic matter concentration and thermodynamic favorability define an upper bound on respiration potential, whereas biological catalytic capacity and physical accessibility jointly govern realized respiration rates through interaction effects. To represent unmeasurable accessibility constraints, we use the mechanistic model as a scaffold for KGML, allowing machine learning to target residual structure not explained by process theory. This hybrid framework improves predictive skill relative to both the mechanistic model alone and fully data-driven models while preserving interpretability. These results indicate that variability in hyporheic respiration is largely mechanistically structured and demonstrate how integrating process theory with explainable AI enhances predictive performance while enabling scalable synthesis of river corridor observations.

Zheng, Jianqiu↗

Predicting the tensile and compressive modulus of electrospun fiber mat‐reinforced hydrogels using the Halpin–Tsai equations

Abstract The reinforcement of mechanically‐weak hydrogels to yield composites with increased stiffness, strength, or toughness is a well‐established approach. In particular, introducing electrospun nanofibers into hydrogels is a common strategy for biomedical applications, as the resulting hierarchical structure mimics biology and allows for control over fiber diameter and alignment and tuning of mechanical properties. However, further study of the link between the constituent materials and the mechanical properties of the composite is uncommon. One potential model to understand the mechanical properties of fiber‐reinforced hydrogels involves the Halpin–Tsai equations, which relate the modulus values of the fibers and hydrogel matrix and the fiber volume fraction, to the modulus of the composite. To assess the application of this model to fiber‐reinforced hydrogels, predicted values were compared with experimental values from mechanical testing of a poly(ethylene glycol) (PEG) matrix reinforced with an electrospun polycaprolactone (PCL) fiber mat. Although the equations described these systems well in tension, providing a facile approach to identify a fiber volume fraction that will achieve a desired modulus, the Halpin–Tsai approach was less successful under compression. This study motivates additional investigation of the role of structural features of hydrogel composites in determining mechanical properties to enable design of materials for specific applications.

Polymer Science↗

Controls on Soil Organic Matter Degradation and Subsequent Greenhouse Gas Emissions Across a Permafrost Thaw Gradient in Northern Sweden

Warming-induced permafrost thaw could enhance microbial decomposition of previously stored soil organic matter (SOM) to carbon dioxide (CO 2 ) and methane (CH 4 ), one of the most significant potential feedbacks from terrestrial ecosystems to the atmosphere in a changing climate. The environmental parameters regulating microbe-organic matter interactions and greenhouse gas (GHG) emissions in northern permafrost peatlands are however still largely unknown. The objective of this work is to understand controls on SOM degradation and its impact on porewater GHG concentrations across the Stordalen Mire, a thawing peat plateau in Northern Sweden. Here, we applied high-resolution mass spectrometry to characterize SOM molecular composition in peat soil samples from the active layers of a Sphagnum-dominated bog and rich fen sites in the Mire. Microbe-organic matter interactions and porewater GHG concentrations across the thaw gradient were controlled by aboveground vegetation and soil pH. An increasingly high abundance of reduced organic compounds experiencing greater humification rates due to enhanced microbial activity were observed with increasing thaw, in parallel with higher CH 4 and CO 2 porewater concentrations. Bog SOM however contained more Sphagnum-derived phenolics, simple carbohydrates, and organic- acids. The low degradation of bog SOM by microbial communities, the enhanced SOM transformation by potentially abiotic mechanisms, and the accumulation of simple carbohydrates in the bog sites could be attributed in part to the low pH conditions of the system associated with Sphagnum mosses. We show that Gibbs free energy of C half reactions based on C oxidation state for OM can be used as a quantifiable measure for OM decomposability and quality to enhance current biogeochemical models to predict C decomposition rates. We found a direct association between OM chemical diversity and δ 13 C-CH 4 in peat porewater; where higher substrate diversity was positively correlated with enriched δ 13 C-CH 4 in fen sites. Oxidized sulfur-containing compounds, produced by Sphagnum, were further hypothesized to control GHG emissions by acting as electron acceptors for a sulfate-reducing electron transport chain, inhibiting methanogenesis in peat bogs. These results suggest that warming-induced permafrost thaw might increase organic matter lability, in subset of sites that become wetlands, and shift biogeochemical processes toward faster decomposition with an increasing proportion of carbon released as CH 4 .

54 ENVIRONMENTAL SCIENCES↗