Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability prediction”

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 919 records · Page 51

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

Mg-Ion Conduction in Antiperovskite Solid Electrolytes Revealed by 25 Mg Ultrahigh Field NMR and First-Principles Calculations

Magnesium-ion batteries hold the potential to outperform the energy density of lithium-ion batteries, given the divalent charge carried by each Mg 2+ cation, but remain in an early stage of development. Here, in this study, 25 Mg solid-state nuclear magnetic resonance (ssNMR) is used to gain insight into the local structure and Mg-ion dynamics of candidate Mg-ion solid electrolytes, the antiperovskites Mg 3 SbN and Mg 3 AsN. Using the highest available magnetic field (35.2 T) for high-resolution solid-state NMR, the largest 25 Mg quadrupole coupling constants (C Q ) yet measured of up to 22 MHz are reported and corroborated by first-principles calculations. Predicted C Q values are shown to correlate with the antiperovskite’s tolerance factor; thus, 25 Mg NMR linewidths can report on lattice distortions and phase stability of these antiperovskites. Variable-temperature 25 Mg NMR spectra demonstrate changes at elevated temperatures, ascribed to Mg-ion motional effects. 25 Mg T 1 relaxometry measurements at ultrahigh field reveal a lower activation energy for the more distorted Mg 3 AsN phase, matching computational predictions of a lower energy barrier for Mg 2+ ion migration and suggesting that additional scrutiny of antiperovskites as Mg-ion conductors is warranted. Given the inherent challenges of 25 Mg NMR, this work demonstrates the benefits of combining ultrahigh field NMR spectroscopy, advanced pulse sequences, modern signal processing, and first-principles calculations to facilitate NMR of quadrupolar nuclei as a tool to probe the local structure and ion dynamics in beyond-Li battery materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting oxidation damage in ultra high-temperature borides: A machine learning approach

Ultra-high temperature (UHT) borides are ceramics materials with melting points above 3000 °C for structural applications in extreme environments. However, at temperatures exceeding 1600 °C and under oxidizing conditions, the material suffers from detrimental degradation. Optimized design and performance of diboride materials under such extreme conditions requires filling the missing composition-microstructure-oxidation gap. This study proposes a computational data-driven framework to connect the processing and microstructure of Ultra-high temperature borides with the oxidation damage. Random Forest Regressor (RFR) model is adopted to forecast the oxide scale thickness developed after oxidation testing based on processing variables and microstructural features. The model trained on a dataset consisting of 107 samples of experimental data extracted from the literature aims to predict oxidation damage. With proper data manipulation and fine model tuning, the predictor could forecast the oxide scale thickness of UHT diborides with a Mean Absolute Error of 37.45 μm and an R-square of 0.83. This model could be used as a high-throughput scheme to design and test new UHT diborides materials computationally. Furthermore, a model with larger composition capabilities could also be developed in the future as more experimental data become available.

36 MATERIALS SCIENCE↗

Juvenile life history diversity is associated with lifetime individual heterogeneity in a migratory fish

Abstract Differences in the life history pathways (LHPs) of juvenile animals are often associated with differences in demographic rates in later life stages. For migratory animals, different LHPs often result in animals from the same population occupying distinct habitats subjected to different environmental drivers. Understanding how demographic rates differ among animals expressing different LHPs may reveal fitness trade‐offs that drive the expression of alternative LHPs and enable better prediction of population dynamics in a changing environment. To understand how demographic outcomes and their relationships with environmental variables differ among animals with different LHPs, we analyzed a long‐term (2006–2021) mark–recapture dataset for Chinook salmon ( Oncorhynchus tshawytscha ) from the Wenatchee River, Washington, USA. Distinct LHPs represented in this population include either remaining in the natal stream until emigrating to the ocean as a 1‐year‐old (natal‐reach rearing) or emigrating from the natal stream and rearing in downstream habitats for several months before completing the emigration to the ocean as a 1‐year‐old (downstream rearing). We found that downstream‐rearing fish emigrated to the ocean 19 days earlier on average and returned as adults from the ocean at higher rates. We detected a positive correlation between rate of return from the ocean by downstream‐rearing fish and coastal upwelling in their spring of outmigration, whereas for natal‐reach‐rearing fish we detected a positive correlation with sea surface temperature during their first marine summer. Different responses to environmental variability should lead to asynchrony in adult abundance among juvenile LHPs. A higher proportion of downstream‐rearing fish returned at younger ages compared with natal‐reach‐rearing fish, which contributed to variability in age at reproduction and greater mixing across generations. Our results demonstrate how diversity in juvenile LHPs is associated with heterogeneity in demographic rates during subsequent life stages, which can in turn affect variance in aggregate population abundance and response to environmental change. Our findings underscore the importance of considering life history diversity in demographic analyses and provide insights into the effects of life history diversity on population dynamics and trade‐offs that contribute to the maintenance of life history diversity.

Sorel, Mark H.↗

Development of a variable tensioning system to reduce separation force in large scale stereolithography

Projection micro stereolithography (PµSL) is an additive manufacturing tool that offers multiple advantages, including unparalleled resolution and throughput, but the ability to print high viscosity resin for large-scale parts is limited. One of the key challenges in PµSL is to separate a newly polymerized layer from the vat floor without damaging the part. Since the separation force scales with the printing area, the risk of damaging the part increases significantly with larger-scale systems and must be addressed. In this paper, a novel roll-to-roll, variable tensioning system is proposed to reduce the separation force during printing. A mathematical model is proposed to predict the separation force for different 2D geometries, and a set of experiments is conducted on an experimental prototype to validate the model. Finally, the effects of different separation parameters including peel rate and peel angle are discussed in detail. The results show that the proposed tensioning system reduces the separation force significantly.

36 MATERIALS SCIENCE↗

Charge Carrier Transport in Iron Pyrite Thin Films: Disorder-Induced Variable-Range Hopping

Here, the origin of p-type conductivity and the mechanism responsible for low carrier mobility were investigated in pyrite (FeS 2 ) thin films. Temperature-dependent resistivity measurements (10-400 K) were performed on polycrystalline and nanostructured thin films prepared by three different methods: (1) spray pyrolysis, (2) hot-injection synthesized and spin-coated nanocubes, and (3) pulsed laser deposition. The films have a high hole density (10 18 -10 19 ) cm -3 and low mobility (0.1-4 cm 2 V -1 s -1 ) regardless of the method used for their preparation. The charge transport mechanism is determined to be thermally activated conduction (TAC) at near room temperature, with Mott-type variable-range hopping (VRH) of holes via localized states occurring at lower temperatures. The density functional theory (DFT) predicts that sulfur vacancy induces localized defect states within the band gap and the charge remains localized around the defect. The data indicates that the electronic properties including hopping transport in pyrite thin films can be correlated to sulfur vacancy-related defects. The results provide insights into the electronic properties of pyrite thin films and their implications for charge transport.

36 MATERIALS SCIENCE↗

Discovering type I cis-AT polyketides through computational mass spectrometry and genome mining with Seq2PKS

Type 1 polyketides are a major class of natural products used as antiviral, antibiotic, antifungal, antiparasitic, immunosuppressive, and antitumor drugs. Analysis of public microbial genomes leads to the discovery of over sixty thousand type 1 polyketide gene clusters. However, the molecular products of only about a hundred of these clusters are characterized, leaving most metabolites unknown. Characterizing polyketides relies on bioactivity-guided purification, which is expensive and time-consuming. To address this, we present Seq2PKS, a machine learning algorithm that predicts chemical structures derived from Type 1 polyketide synthases. Seq2PKS predicts numerous putative structures for each gene cluster to enhance accuracy. The correct structure is identified using a variable mass spectral database search. Benchmarks show that Seq2PKS outperforms existing methods. Applying Seq2PKS to Actinobacteria datasets, we discover biosynthetic gene clusters for monazomycin, oasomycin A, and 2-aminobenzamide-actiphenol.

60 APPLIED LIFE SCIENCES↗

Spatial variability in Arctic–boreal fire regimes influenced by environmental and human factors

Abstract Wildfire activity in Arctic and boreal regions is rapidly increasing, with severe consequences for climate and human health. Regional long-term variations in fire frequency and intensity characterize fire regimes. The spatial variability in Arctic–boreal fire regimes and their environmental and anthropogenic drivers, however, remain poorly understood. Here we present a fire tracking system to map the sub-daily evolution of all circumpolar Arctic–boreal fires between 2012 and 2023 using 375 m Visible Infrared Imaging Radiometer Suite active fire detections and the resulting dataset of the ignition time, location, size, duration, spread and intensity of individual fires. We use this dataset to classify the Arctic–boreal biomes into seven distinct ‘pyroregions’ with unique climatic and geographic environments. We find that these pyroregions exhibit varying responses to environmental drivers, with boreal North America, eastern Siberia and northern tundra regions showing the highest sensitivity to climate and lightning density. In addition, anthropogenic factors play an important role in influencing fire number and size, interacting with other factors. Understanding the spatial variability of fire regimes and its interconnected drivers in the Arctic–boreal domain is important for improving future predictions of fire activity and identifying areas at risk for extreme events.

Geology↗

A Multi-Instrument Cloud Condensation Nuclei Spectrum Product (Final Technical Report)

A wealth of observational data exists on the characteristics of atmospheric particulate matter, over multiple years, at the DOE ARM Southern Great Plains (SGP) Central Facility site. This site is located in a region of the country that frequently experiences weather extremes, and that is removed from many local sources of pollution but is affected by transported smoke, dust, and urban emissions. The relationships between particulate matter, cloud formation and evolution, and precipitation are therefore of strong interest, and are being explored via modeling on a variety of scales. These models require as input detailed information on the characteristics of particles capable of serving as the nuclei for cloud formation. Sufficient data exist to be able to put together a picture of the nature of the total aerosol and the cloud condensation nuclei (CCN) subset, and their variability, through merged data products. This study was aimed at exploiting the multiple measurement types at SGP to develop the first such multi-year estimates. Further, the resulting data were analyzed to understand temporal patterns ranging from hourly to seasonal, thereby gaining insights into the particle sources affecting the atmosphere in this region. DOE-funded datasets that were analyzed in this study include total particle number concentrations, submicron aerosol scattering coefficients, dry aerosol size distributions, and more recently, time-resolved submicron aerosol chemical composition. Data are also available for the number concentrations of particles that are activated in a cloud condensation nucleus instrument at a series of setpoint supersaturations, providing direct observations of the number concentrations of “CCN”. This variable is the quantity that is generally desired for inclusion in numerical models that seek to represent and predict the impacts of varying aerosol characteristics on the formation and microphysical properties of clouds. One limitation of the use of direct CCN observations is that they are not available for supersaturations larger than about 1%, which is insufficient for deep convection and may be insufficient even for shallow convection, depending on the nature of the available CCN and the dynamics of the cloud. We developed a data-based approach to representing the full aerosol size spectrum with size-dependent hygroscopicity, and used this to extrapolate CCN spectra beyond the limited measurements. Five years of SGP aerosol data (2009 -2013) were analyzed. As a side product of our work, we identified and communicated several previously-unflagged data quality issues. The resulting merged aerosol distributions, along with fits for seasonal averages, were published and submitted to the ARM archive as a special value-added product (VAP; submitted as a PI product). CCN spectra were computed for the same data period and will similarly be published and submitted to the archive for use by the community. We also note that our methodologies and findings have been discussed at several Joint ARM User Facility/Atmospheric System Research (ASR) Principal Investigators Meetings and that recent ARM/ASR aerosol data reporting strategies have included similar ideas for data merging, indicating that this work has had a lasting impact on ARM aerosol data acquisition and reporting. The proposed work advances the science of the interactions of aerosols, clouds and precipitation, with direct application to improve representation of such interactions for clouds in regional and global climate models. The archived data will continue to serve research studies in the future.

54 ENVIRONMENTAL SCIENCES↗

Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]

Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

Quantifying the drivers and predictability of seasonal changes in African fire

Africa contains some of the most vulnerable ecosystems to fires. Successful seasonal prediction of fire activity over these fire-prone regions remains a challenge and relies heavily on in-depth understanding of various driving mechanisms underlying fire evolution. Here, we assess the seasonal environmental drivers and predictability of African fire using the analytical framework of Stepwise Generalized Equilibrium Feedback Assessment (SGEFA) and machine learning techniques (MLTs). The impacts of sea-surface temperature, soil moisture, and leaf area index are quantified and found to dominate the fire seasonal variability by regulating regional burning condition and fuel supply. Compared with previously-identified atmospheric and socioeconomic predictors, these slowly evolving oceanic and terrestrial predictors are further identified to determine the seasonal predictability of fire activity in Africa. Our combined SGEFA-MLT approach achieves skillful prediction of African fire one month in advance and can be generalized to provide seasonal estimates of regional and global fire risk.

54 ENVIRONMENTAL SCIENCES↗

Advancing Sea Ice Predictability in E3SM with Machine Learning

Focal area(s): To improve predictions of sea ice in E3SM we propose to develop a hierarchy of data-driven models using observational and simulation data to investigate the most important Earth system drivers of sea ice variability and loss, develop surrogates that build on the reduced parameter space of important drivers, and, where appropriate, couple machine learning models with standard PDE models to capture important physical behavior at different scales. This work falls under Focal Area 2. Predictive modeling through the use of AI techniques.

54 ENVIRONMENTAL SCIENCES↗

Knowledge-guided machine learning can improve carbon cycle quantification in agroecosystems

Abstract Accurate and cost-effective quantification of the carbon cycle for agroecosystems at decision-relevant scales is critical to mitigating climate change and ensuring sustainable food production. However, conventional process-based or data-driven modeling approaches alone have large prediction uncertainties due to the complex biogeochemical processes to model and the lack of observations to constrain many key state and flux variables. Here we propose a Knowledge-Guided Machine Learning (KGML) framework that addresses the above challenges by integrating knowledge embedded in a process-based model, high-resolution remote sensing observations, and machine learning (ML) techniques. Using the U.S. Corn Belt as a testbed, we demonstrate that KGML can outperform conventional process-based and black-box ML models in quantifying carbon cycle dynamics. Our high-resolution approach quantitatively reveals 86% more spatial detail of soil organic carbon changes than conventional coarse-resolution approaches. Moreover, we outline a protocol for improving KGML via various paths, which can be generalized to develop hybrid models to better predict complex earth system dynamics.

54 ENVIRONMENTAL SCIENCES↗

The Q 10 of in situ microbial soil respiration varies with mean annual temperature, precipitation, pH, and plant cover: a meta-analysis and spatial prediction of Q 10

The temperature sensitivity of soil microbial respiration, commonly quantified using the Q 10 coefficient, is a key parameter in carbon cycle models. Uncovering how environmental factors affect in situ Q 10 values can therefore provide critical insight into potential shifts in global carbon stocks under climate change. We collected data from previously published field experiments that measured soil microbial respiration across a range of temperatures. We hypothesized that the Q 10 coefficient of in situ soil microbial respiration would vary based on environmental factors including mean annual temperature (MAT), mean annual precipitation (MAP), plant cover type, pH, soil C:N, and latitude. Linear regression revealed that Q 10 correlates negatively with MAT and MAP and positively with pH and absolute latitude. Additionally, average Q 10 varied significantly across different plant cover types; it was highest in mountain grasslands and lowest in tropical moist forests. Variation in microbial Q 10 across environmental factors may arise from underlying mechanisms such as enzyme kinetics, substrate availability and complexity, and microbial adaptation. To capture patterns in Q 10 more comprehensively, we developed a multiple linear regression model of Q 10 based on the most individually significant environmental drivers and applied it to public datasets to generate a global map of predicted Q 10 . Q 10 was higher in high-latitude and high-altitude regions, where large permafrost carbon stores are vulnerable to thawing and decomposition. We also compared fits between the Q 10 equation and a model produced from macromolecular rate theory (MMRT). We found that the MMRT model had the superior fit and may be better suited to model temperature sensitivity of complex biological reactions. Overall, our results emphasize that relationships between microbial Q 10 and environmental variables should be accounted for in climate models. Incorporating these variations in the Q 10 parameter, rather than using a fixed value, will help predict whether CO 2 emissions will be buffered or exacerbated by soil microbial respiration under climate change.

54 ENVIRONMENTAL SCIENCES↗