Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Extreme scale”

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 271 records · Page 15

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗

Structure of the sunspot penumbra

An exceptionally highly resolved sunspot photograph is presented which reveals the fine-scale structure of the photospheric penumbra down to 0.2 arcsec. The photograph was taken through a glass filter in a 170-A FWHM band centered on 4660 A on the 65-cm vacuum Gregorian telescope at Big Bear Solar Observatory during practically perfect seeing. The photograph reveals the penumbra to be full of extremely fine-scale fibril structures, with apparent widths near the resolution limit. The three-dimensional structure of the photospheric penumbra is seen to be similar to that of the chromospheric penumbra and superpenumbra, with dark fibrils arching above a lower brighter structure and many of the bright fibrils representing portions of wider bright structures shining through the narrow spaces between the elevated dark fibrils. Thus, in contrast to previous pictures, much of the dark component of the photospheric penumbra is not physically analogous to the intergranular dark lanes in the normal photosphere.

Moore, R. L.↗

Seasonal and Geographical Variations in Fundamental Weather Patterns during Extreme Precipitation as Identified from Omega Equation Forcing

Abstract We investigate the large-scale weather patterns during extreme precipitation (PEx) events over the conterminous United States (CONUS) by applying a version of the quasigeostrophic (QG) omega equation. This work aims to develop a climatology of the weather patterns most related to PEx events during current climate. Extreme events are examined for each of seven regions defined by consistent annual cycles of precipitation and spanning the CONUS. For the CONUS we train several self-organizing maps (SOM) on a pressure–time series of vertical velocity from each of the advective forcing terms in the QG omega equation for each extreme event. The unsupervised learning of the SOM allows us to identify the most descriptive set of nine patterns in vertical velocity associated with precipitation extremes. This method finds multiple frontal- and cyclone-driven patterns while grouping primarily convective events into one pattern. Frontal events include a synoptic pattern consistent with West Coast atmospheric river events as well as pattern groups linked to developing and to mature (“occluded”) frontal cyclones. The primary patterns found during PEx events vary seasonally and geographically. Frontal cyclone patterns are most common during PEx events during summer in the part of the Great Plains and during winter for the Northeast, Southeast, Pacific Northwest, and Southwest. Convection is the most common pattern during summer in all regions. Except in the Southeast, the annual cycles of monthly number of PEx events and average precipitation match well, partially validating our choice of regions to aggregate PEx events.

Meteorology & Atmospheric Sciences↗

Emulating climate extreme indices

We use simple pattern scaling and time-shift to emulate changes in a set of climate extreme indices under future scenarios, and evaluate the emulators' accuracy. We propose a metric for the error in emulation in the context of initial condition ensembles, to specifically characterize the role of internal variability in the emulation performance. Our metric separates systematic emulation errors from unavoidable discrepancies between emulated and target values due to internal variability. We compute the metricis at grid-point scale, and we show geographically resolved results, or aggregate them at global scale. We demonstrate the use of our error metric in the emulation of a suite of temperature and precipitation extreme indices. We test and compare simple pattern scaling and time-shift using a range of trajectories spanning targets inspired by the Paris agreement -- warming to 1.5C and 2.0C from the pre-industrial baseline -- and two of the longer-established trajectories, RCP4.5 and RCP8.5. With this suite of scenarios we can test the effects on the performance of the size of the temperature gap between emulation origin and target. We find that for most indices emulation the dominant source of discrepancy is internal variability. For at least one index, however, counting exceedances of a high temperature threshold, significant portions of the globally aggregated discrepancy and its regional pattern originate from the systematic emulation error. This error exceeds internal variability of both the target and the emulated quantities in large coherent regions at low latitudes, and the explanation can be found in the differential behavior of temperature distributions across latitudes. The metric also highlights a fundamental difference in the two methods related to the simulation of internal variability, which is dampened significantly by simple pattern scaling. This aspect is of consequence when using these methods for specific applications, where preserving variability for uncertainty quantification is deemed important. With this study we offer our metric as a diagnostic tool, facilitating the formulation of scientific hypotheses on the reasons for the error. In the meantime, we show that for many impact relevant indices by now traditional emulation techniques can be accurate within the variations unavoidably introduced by internal variability, establishing the fundamental condition for using their emulation in impact modeling.

54 ENVIRONMENTAL SCIENCES↗

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. Here, in the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radar characteristics of Viking 1 landing sites

Radar observations of Mars at centimeter wavelengths in May, June, and July 1976 provided estimates of surface roughness and reflectivity in three potential landing areas for Viking 1. Surface roughness is characterized by the distribution of surface landing slopes or tilts on lateral scales of the order of 1 to 10 meters; measurements of surface reflectivity are indicators of bulk surface density in the uppermost few centimeters. By these measures, the landing site at 47.5 deg W, 22.4 deg N is rougher than the Martian average, although it may be near the average for elevations accessible to Viking, and is estimated to be near the Mars average in reflectivity. The site at the center of Chryse Planitia, 43.5 deg W, 23.4 deg N, may be an area of anomalous radar characteristics, indicative of extreme, small-scale roughness, very low surface density, or a combination of these two characteristics. Observations of the original Chryse site at 34 deg W, 19.5 deg N indicate that that area is at least twice as rough as the Mars average.

Tyler, G. L.↗

Sizing Up the Situation

Hailstorm damage to the Space Shuttle's External Tank inspired a NASA innovation with extensive photography applications. In order to measure the defects caused by the storm, Kennedy Space Center used telephoto lenses to zoom in on the tank and view the damage clearly. However, since there was no reference object in the image, the engineers could not determine the scale of the damage. In photographic situations similar to this, an object, such as a ruler, is placed within the field of view. This allows a person to look at a photograph and have a visual indication of the scale of the objects in it. In the External Tank situation, however, this procedure was not possible. As a solution, Kennedy developed the Scaling and Measurement Device for Photographic Images, which provides a non-intrusive means of adding a scale to a photograph. In addition to meeting Kennedy's needs, scaling images is extremely important in crime and accident scene investigations, oil and chemical tank monitoring, and aerial photography. The innovation consists of a tool that attaches directly to a camera or charge coupled device using a standard screw. Two lasers fitted to the device provide parallel beams that are set 1 inch apart. These lasers enable the device to project a pattern into the field of view. When a photograph is taken, the image of the laser pattern appears, along with the image of the object under investigation, allowing the viewer quantifiable information as to the size of the object. The laser beams are accurate to approximately 200 feet. Windows-based software was developed to work with the scaling device tool. The software provides further techniques to measure objects in photographs and digital images. By using the software, any object in the image can be measured diagonally, vertically, and horizontally. The device and its software enable the user to determine two-dimensional measurements within a photograph.

Source record↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Early Season Large-Area Winter Crop Mapping Using MODIS NDVI Data, Growing Degree Days Information and a Gaussian Mixture Model

Knowledge on geographical location and distribution of crops at global, national and regional scales is an extremely valuable source of information applications. Traditional approaches to crop mapping using remote sensing data rely heavily on reference or ground truth data in order to train/calibrate classification models. As a rule, such models are only applicable to a single vegetation season and should be recalibrated to be applicable for other seasons. This paper addresses the problem of early season large-area winter crop mapping using Moderate Resolution Imaging Spectroradiometer (MODIS) derived Normalized Difference Vegetation Index (NDVI) time-series and growing degree days (GDD) information derived from the Modern-Era Retrospective analysis for Research and Applications (MERRA-2) product. The model is based on the assumption that winter crops have developed biomass during early spring while other crops (spring and summer) have no biomass. As winter crop development is temporally and spatially non-uniform due to the presence of different agro-climatic zones, we use GDD to account for such discrepancies. A Gaussian mixture model (GMM) is applied to discriminate winter crops from other crops (spring and summer). The proposed method has the following advantages: low input data requirements, robustness, applicability to global scale application and can provide winter crop maps 1.5-2 months before harvest. The model is applied to two study regions, the State of Kansas in the US and Ukraine, and for multiple seasons (2001-2014). Validation using the US Department of Agriculture (USDA) Crop Data Layer (CDL) for Kansas and ground measurements for Ukraine shows that accuracies of greater than 90% can be achieved in mapping winter crops 1.5-2 months before harvest. Results also show good correspondence to official statistics with average coefficients of determination R(exp. 2) greater than 0.85.

mixture model↗

High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) towards CMIP7

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models.

54 ENVIRONMENTAL SCIENCES↗

SIF-based GPP modeling for evergreen forests considering the seasonal variation in maximum photochemical efficiency

Solar-induced chlorophyll fluorescence (SIF) has shown great potential in estimating gross primary production (GPP). However, their quantitative relationship is not invariant, which undermines the reliability of empirical SIF-based GPP estimation at fine spatiotemporal scales, especially under extreme conditions. In this study, we developed a parsimonious mechanistic model for SIF-based GPP estimation in evergreen needle forests (ENF) by employing the Mechanistic Light Response framework and Eco-Evolutionary theory to describe the light and dark reactions during photosynthesis, respectively. Specifically, we found that considering the seasonal variation in a key parameter of the MLR framework, the maximum photochemical efficiency of photosystem II (Φ PSIImax ), can avoid the GPP overestimation in winter and early spring due to the relatively low environmental sensitivity of SIF. Compared to the estimates from other benchmark models, our GPP estimates were closer to the 1: 1 line and had higher accuracy (average R 2 = 0.86, RMSE=1.99 μmol m -2 s -1 ) across sites. Furthermore, the changes in the relationship between SIF and J (refers to the electron transport rate) contribute a lot to the dynamic SIF–GPP relationship in this study, while the J–GPP relationship is less variant when the temperature drops. Further, the seasonal variation in the SIF–J relationship, especially the reduction in its slope at low temperatures, is found largely explained by the Φ PSIImax . These results indicate the importance of the uncertainty caused by the variation in the SIF–J relationship for SIF-based GPP estimation, and the consideration of changes in Φ PSIImax under extreme conditions (such as severe winter in this study) is crucial for the improvement of GPP estimation via SIF.

59 BASIC BIOLOGICAL SCIENCES↗

Overview of the Alaskan Layered Pollution and Chemical Analysis (ALPACA) Field Experiment

The Alaskan Layered Pollution And Chemical Analysis (ALPACA) field experiment was a collaborative study designed to improve understanding of pollution sources and chemical processes during winter (cold climate and low-photochemical activity), to investigate indoor pollution, and to study dispersion of pollution as affected by frequent temperature inversions. A number of the research goals were motivated by questions raised by residents of Fairbanks, Alaska, where the study was held. This paper describes the measurement strategies and the conditions encountered during the January and February 2022 field experiment, and reports early examples of how the measurements addressed research goals, particularly those of interest to the residents. Outdoor air measurements showed high concentrations of particulate matter and pollutant gases including volatile organic carbon species. During pollution events, low winds and extremely stable atmospheric conditions trapped pollution below 73 m, an extremely shallow vertical scale. Tethered-balloon-based measurements intercepted plumes aloft, which were associated with power plant point sources through transport modeling. Because cold climate residents spend much of their time indoors, the study included an indoor air quality component, where measurements were made inside and outside a house to study infiltration and indoor sources. In the absence of indoor activities such as cooking and/or heating with a pellet stove, indoor particulate matter concentrations were lower than outdoors; however, cooking and pellet stove burns often caused higher indoor particulate matter concentrations than outdoors. The mass-normalized particulate matter oxidative potential, a health-relevant property measured here by the reactivity with dithiothreiol, of indoor particles varied by source, with cooking particles having less oxidative potential per mass than pellet stove particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How the Great Plains Dust Bowl drought spread heat extremes around the Northern Hemisphere

Abstract Extraordinary heat extremes occurred in the 1930s in areas of the Northern Hemisphere far from the record setting heat over the US associated with the Great Plains Dust Bowl drought. A climate model sensitivity experiment is used to identify a new mechanism involving a warm season circumglobal atmospheric teleconnection pattern that spread heat extremes over far-flung areas of the Northern Hemisphere arising from the intense heating over the desiccated Great Plains themselves. It has only been in the twenty-first century that human populations in these regions of the Northern Hemisphere have experienced heat extremes comparable to the 1930s. This demonstrates that humans influenced Northern Hemisphere temperature and heat extremes through disastrous and unprecedented regional land use practices over the Great Plains, and points to the possibility that future intense regional droughts could affect heat extremes on hemispheric scales.

54 ENVIRONMENTAL SCIENCES↗

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗

The Design of Relativistic Ultrafast Electron Diffraction and Imaging (RUEDI) Facility for Materials in Extremes

Materials in extreme conditions in terms of temperature and pressure are of great interest where novel phenomena are expected for new science and materials design. A wide range of materials are of interest, including warm dense matter, particulate materials, liquid materials, fusion materials and geophysical materials. Their environments also range from liquid phase, gas phase and cryogenic temperatures. Here, the conditions often correspond to the level of the interiors of planets or stars. Such conditions are, however, typically not available in the laboratory. In addition, novel phenomena exist as transition structures that occur at an extremely short time scale in the harsh environments, which are often difficult to be modelled.

36 MATERIALS SCIENCE↗

EXAGRAPH: Graph and combinatorial methods for enabling exascale applications

Combinatorial algorithms in general and graph algorithms in particular play a critical enabling role in numerous scientific applications. However, the irregular memory access nature of these algorithms makes them one of the hardest algorithmic kernels to implement on parallel systems. With tens of billions of hardware threads and deep memory hierarchies, the exascale computing systems in particular pose extreme challenges in scaling graph algorithms. The codesign center on combinatorial algorithms, ExaGraph, was established to design and develop methods and techniques for efficient implementation of key combinatorial (graph) algorithms chosen from a diverse set of exascale applications. Algebraic and combinatorial methods have a complementary role in the advancement of computational science and engineering, including playing an enabling role on each other. In this paper, we survey the algorithmic and software development activities performed under the auspices of ExaGraph from both a combinatorial and an algebraic perspective. In particular, we detail our recent efforts in porting the algorithms to manycore accelerator (GPU) architectures. We also provide a brief survey of the applications that have benefited from the scalable implementations of different combinatorial algorithms to enable scientific discovery at scale. We believe that several applications will benefit from the algorithmic and software tools developed by the ExaGraph team.

97 MATHEMATICS AND COMPUTING↗

Root Influences on Mobilization and Export of Mineral-bound Soil Organic Matter

Biogeochemical cycles within mountainous watersheds are key regulators of ecosystem carbon storage and downstream nutrient loadings, and they have shown to be particularly vulnerable to climate change impacts. Increasing temperature and persistent droughts have already dramatically changed vegetation cover across the mountainous western US, with unknown consequences for soil carbon and nutrient cycles belowground. What remains elusive is to what extent associated changes in root-soil interactions may mobilize the vast pool of organic matter (OM) that has been stabilized by associations with minerals for centuries or millennia. Although plant root-driven OM mobilization from minerals may be a central control on carbon loss and nutrient export, such mechanisms are currently missing from conceptual and numerical models. The overall objective of this Exploratory Project is to identify the biogeochemical mechanisms by which roots destabilize mineral-associated organic matter (MAOM) and the cumulative impact on carbon and nutrient fate. To accomplish this goal, we integrated well-controlled laboratory experiments with in-field measurements and a scalable modeling approach. First, we conducted model system experiments to assess the vulnerability of MAOM to exudate-mediated mechanisms. Our results show that common root exudates effectively destabilize MAOM not only through direct, ligand-driven mobilization mechanisms, but also indirect, microbially-mediated mechanisms relying on secondary metabolites and enzymes. We further found that OM bound to poorly crystalline Fe and Al (hydr)oxides is more vulnerable to exudate-induced destabilization than OM bound to more crystalline phases, particularly in response to direct, ligand-promoted mechanisms. These findings demonstrate that the stability of MAOM is not just a function of their inherent properties, but also will depend in large parts on the ability of plant roots and microbes to produce exudates capable of triggering suitable mobilization mechanisms. We further employed a well-controlled rhizobox approach, combining advanced microsensor and mass spectrometry techniques, to resolve spatiotemporal variations in the composition and availability of exudates along single growing roots of grasses. Our results show that the composition of functionally relevant exudate compounds varies at extremely short time scales, seemingly shifting from ligands such as aromatic acids around root tips to less reactive metabolites such as amino acids around mature root segments. These results suggest a prevalence of direct MAOM mobilization mechanisms around the root tip, while indirect MAOM mobilization strategies may dominate around more mature root segments. Additionally, we utilized microsensor measurements in rhizoboxes to parameterize a rhizosphere (hydro)biogeochemistry reactive transport model (eSTOMP-ROOTS). The resulting model was used to assess how root exudations affects the stability of MAOM. Model simulations show that diel pulses of root exudation are strong enough to cause oscillations in biogeochemical conditions—particularly in pH, oxygen concentrations, and microbial activity— that repeatedly disrupt MAOM. Finally, we assessed field-scale impacts of root-driven MAOM destabilization over short and long time scales. A combination of in-field incubations and field-based characterizations showed that MAOM may not only be vulnerable to seasonal variations in root activity, but also to root-driven weathering over pedogenic time scales. In sum, results from our Exploratory Award highlight the strong control plant roots exert on MAOM mobilization and, thus, on the potential for carbon and nutrient export from watersheds.

54 ENVIRONMENTAL SCIENCES↗