Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial transferability”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evidence of Alfvenic Poynting Flux as the Primary Driver of Auroral Motion During a Geomagnetic Substorm

Geomagnetic substorms are major energy transfer events where energy stored in the Earth's magnetotail is released into the ionosphere. Substorm phenomena, including auroral activities, earthward Poynting flux, magnetic field dipolarization, etc, have been extensively studied. However, the complex interplay among them is not fully understood. In a fortuitous event on June 7, 2013, the twin Van Allen Probes (separated by 0.4 h in local time) observed bursts of earthward Alfvenic Poynting flux in the vicinity of the plasma sheet boundary layer (PSBL). The Poynting flux bursts correlate with enhancements of auroral brightness around the footpoints of both spacecraft. This indicates a temporal and spatial correlation between the auroral brightening and Poynting flux bursts, and that the auroral motion is directly linked to the perpendicular expansion of the Alfven wave. These observations suggest that the Alfvenic Poynting flux is a primary driver for the auroral electron acceleration. Furthermore, around the time of auroral brightening, a dipolarization was seen to propagate more than 4 h in local time during a 20 min period. The azimuthal phase speed of this dipolarization (2 deg/min) is too small to explain the azimuthal motion of the aurora (13.6 deg/min), but the dipolarization could be related to the generation of the Alfvenic Poynting flux through phase mixing at strong density gradients like those in the PSBL.

58 GEOSCIENCES↗
DOE OSTI2026

Multiscale Cryo Electron Microscopy Reveals Interfacial Degradation and Stabilization in Battery Electrodes

Electrochemical interfaces are dynamic systems, evolving based on their local environment and reactant surface structures. The electrode-electrolyte interface in Li-ion batteries can be protective, limiting parasitic reactions with the electrolyte to passivate the surface [1]. Additionally, this interphase has an impact on the Li-ion transport through that layer based on its composition, bonding environment, and thickness. These parameters are challenging to collect and may vary depending on the electrode surface site investigated relative to its spatial position in a coin cell. This study will detail a multiscale cryogenic electron microscopy approach where millimeter-scale cross-sections through the coin cell batteries were made using a cryogenic stage within a fs-laser plasma focused ion beam (laser PFIB) with complementary energy dispersive X-ray spectroscopy able to detect variations in the composition at electrode interfaces [2]. Microscale cross-sectioning and lamella sample preparation of battery electrodes was conducted at the Center for Integrated Nanotechnologies using a Ga-ion focused ion beam (FIB) with air-free and cryo-transfer [3], followed by nanoscale mapping of composition and bonding within the CEI through cryo-scanning transmission electron microscopy (cryo-STEM) electron energy loss spectroscopy [4]. This multiscale approach enabled identification of millimeter-scale features of a battery stack with visualization of degradation in electrodes such as cracks in cathode particles, gas evolution, and SEI evolution; microscale interfacial characteristics, such as heterogeneity in the SEI or barrier layer and identification of electrolyte networks to the electrode surfaces; and nanoscale measurement of the CEI thickness, mapping of transition metal bonding within the cathode particles to identify loss of active materials, and identification of beneficial electrolyte additives incorporated into the CEI structure. This multiscale approach allows for a statistical understanding of the primary mechanisms and parasitic degradation pathways that impact performance by limiting the ion transport pathways within Li+ batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anodic Dissolution Rates Accelerate with Decreasing MoS 2 Nanoflake Thickness

Electrochemical gating of 2D transition metal dichalcogenide (TMD) electrodes is an emerging frontier in the field of semiconductor electrochemistry. In this approach, an applied bias modifies the charge carrier concentration of the 2D TMD, causing band edge shifts and drastic changes in charge transfer rates. However, leveraging this effect for (photo)electrocatalysis is practically limited by the stability of the TMD material under gating conditions. Gerischer showed anodic dissolution of bulk TMD electrodes can occur in the dark and hypothesized that the reaction proceeds via an electron tunneling mechanism from surface states to the TMD conduction band [H. Gerischer, D. Ross, and M. Lubke, Z. Physickalische Chem., 139, 1 (1984)]. Here we investigate this possibility in single MoS 2 nanoflakes using in situ optical microscopy and explore whether Gerischer’s electron tunneling mechanism can explain anodic dissolution rates of thin 2D semiconductors. Here, spatially resolved measurements show anodic dissolution initiates at perimeter edge sites and accelerates exponentially with decreasing layer thickness, consistent with Gerischer’s tunneling mechanism. Interestingly, single layer MoS 2 is impervious to anodic dissolution at applied potentials >200 mV more positive than those required to drive dissolution in bulk and multilayer-thick nanoflakes.

2D semiconductor↗
DOE OSTI2021

Assessment of Contaminant Bioaccumulation in Aquatic Biota on and Adjacent to the Oak Ridge Reservation—2020

This report provides information on contaminant concentrations in multiple wildlife prey species inhabiting or associated with water bodies on and downstream from the Oak Ridge Reservation (ORR), including regional reference sites. This information can be used to understand the nature and extent of contaminant exposure and transfer through the food web and to provide a baseline for evaluating temporal trends in contaminant exposure and accumulation. This information was gathered as part of the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA) Five-Year Review (FYR) process, and it also addresses regulator comments on the annual Remediation Effectiveness Report (RER). This report summarizes the results of biological sampling completed in the spring and summer of 2020, prior to the 2021 FYR. Sampled locations included key sites upstream from the point sources of contaminants, along with watershed integration points on and off the ORR. Integration points are spatially the most important sampling sites because they capture biological responses to upstream contaminant inputs from both aqueous and groundwater sources, as well as responses to various remedial actions. These sites are most often located at a watershed administration unit boundary and may have regulatory significance for performance measures. In many cases, there are long-term biological and other monitoring data associated with these sites; these data generally fall within the areas of exposure and effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing electrogenetic activation via a network model of biological signal propagation

Introduction: Molecular communication is the transfer of information encoded by molecular structure and activity. We examine molecular communication within bacterial consortia as cells with diverse biosynthetic capabilities can be assembled for enhanced function. Their coordination, both in terms of engineered genetic circuits within individual cells as well as their population-scale functions, is needed to ensure robust performance. We have suggested that “electrogenetics,” the use of electronics to activate specific genetic circuits, is a means by which electronic devices can mediate molecular communication, ultimately enabling programmable control. Methods: Here, we have developed a graphical network model for dynamically assessing electronic and molecular signal propagation schemes wherein nodes represent individual cells, and their edges represent communication channels by which signaling molecules are transferred. We utilize graph properties such as edge dynamics and graph topology to interrogate the signaling dynamics of specific engineered bacterial consortia. Results: We were able to recapitulate previous experimental systems with our model. In addition, we found that networks with more distinct subpopulations (high network modularity) propagated signals more slowly than randomized networks, while strategic arrangement of subpopulations with respect to the inducer source (an electrode) can increase signal output and outperform otherwise homogeneous networks. Discussion: We developed this model to better understand our previous experimental results, but also to enable future designs wherein subpopulation composition, genetic circuits, and spatial configurations can be varied to tune performance. We suggest that this work may provide insight into the signaling which occurs in synthetically assembled systems as well as native microbial communities.

Chun, Kayla↗
DOE OSTI2021

Estimation of soil classes and their relationship to grapevine vigor in a Bordeaux vineyard: advancing the practical joint use of electromagnetic induction (EMI) and NDVI datasets for precision viticulture

Working within a vineyard in the Pessac Léognan Appellation of Bordeaux, France, this study documents the potential of using simple statistical methods with spatially-resolved and increasingly available electromagnetic induction (EMI) geophysical and normalized difference vegetation index (NDVI) datasets to accurately estimate Bordeaux vineyard soil classes and to quantitatively explore the relationship between vineyard soil types and grapevine vigor. First, co-located electrical tomographic tomography (ERT) and EMI datasets were compared to gain confidence about how the EMI method averaged soil properties over the grapevine rooting depth. Then, EMI data were used with core soil texture and soil-pit based interpretations of Bordeaux soil types (Brunisol, Redoxisol, Colluviosol and Calcosol) to estimate the spatial distribution of geophysically-identified Bordeaux soil classes. A strong relationship (r = 0.75, p < 0.01) was revealed between the geophysically-identified Bordeaux soil classes and NDVI (both 2 m resolution), showing that the highest grapevine vigor was associated with the Bordeaux soil classes having the largest clay fraction. The results suggest that within-block variability of grapevine vigor was largely controlled by variability in soil classes, and that carefully collected EMI and NDVI datasets can be exceedingly helpful for providing quantitative estimates of vineyard soil and vigor variability, as well as their covariation. The method is expected to be transferable to other viticultural regions, providing an approach to use easy-to-acquire, high resolution datasets to guide viticultural practices, including routine management and replanting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗
DOE OSTI2020

Revolutionizing Turbine Cooling with Micro-Architectures Enabled by Direct Metal Laser Sintering

The objective of this research effort was to explore innovative cooling architectures enabled by additive manufacturing techniques for improved turbine cooling performance. The ability to create complex internal geometries was leveraged to better distribute coolant as well as to integrate inherently unstable flow devices to enhance internal and external heat transfer. This was accomplished with a multi-faceted approach including analytical, experimental, and computational components. This final report documents progress during the total 4.25 year effort that was extended (at no additional cost) from the original 3-year cooperative agreement. Prior to starting the project, significant effort was invested in investigating innovative cooling designs from the open literature. When the project started, we quickly down-selected to 4 configurations: impingement jet fluidic oscillators (for internal leading edge cooling), reverse film cooling (for pressure surface), sweeping film cooling jets (for suction surface), and trailing edge slot cooling with microchannels. During the 2nd year, the reverse film cooling design was shown to be sub-optimal and work on that topic was halted. The three remaining technologies were integrated into a large scale nozzle guide vane and installed into a low-speed linear cascade facility for further interrogation. Initial results showed the sweeping film cooling jets (on the suction surface) to be the most promising technology compared to a baseline diffusion shaped film hole (777 design). The benefit was particularly evident at high blowing ratios (>1.5) when the 777 coolant film separated from the downstream surface. Benefits were also evident at elevated freestream turbulence levels. The leading edge cooling with unsteady jets was less effective in terms of peak or average cooling – however it was superior to round hole impingement cooling in terms of spatial uniformity. Finally, the trailing edge cooling design with micro-channels was scrapped in favor of a pinned arrangement with centerbody. During the 3rd year, the same three technologies were integrated into a transonic linear cascade for an assessment of compressibility effects. Again, the sweeping film cooling jets proved superior at high blowing ratios while the unsteady leading edge impingement and trailing edge pin-fins with centerbody designs yielded mixed results. Finally, during the 4th and final year, these 3 technologies were integrated into a direct metal laser sintered (DMLS) nozzle guide vane for testing in a high temperature, transonic annular vane cascade. The facility matches the flow temperature, Mach number, and coolant pressure ratios of an actual gas turbine. Before doing this, an additional series of tests were completed to validate a redesigned trailing edge cooling architecture. This final design included a centerbody with triangular pins between the centerbody and the vane external skin. The majority of these triangular pins are fabricated with a 30% gap to the centerbody – reducing their pressure drop considerably while still providing excellent heat transfer augmentation. This additional test campaign required an additional 3-month extension request to complete testing in the high temperature NGV test facility (Turbine Reacting Flow Rig – TuRFR). Testing of the DMLS vane in TuRFR pitted the 3 innovative cooling technologies against more traditional technologies (shaped 777 film hole, round impingement jet, and full pin arrangement in the trailing edge. The sweeping film cooling jets (on the suction surface) showed approximately 15% improvement in overall effectiveness compared to the 777 film hole while the trailing edge showed up to 20% improvement. The leading edge impingement was disappointing with a substantially lower effectiveness than traditional direct impingement. In summary, the study was successful in demonstrating that DMLS-enabled cooling technologies can yield significant gains in cooling performance.

03 NATURAL GAS↗
DOE OSTI2025

Alpha-relaxation by scattering and medium-range atomic correlation in simple liquids

The relaxation dynamics of liquid and glass can be studied by inelastic x-ray or neutron scattering through the intermediate scattering function F(Q, t), where Q is the momentum transfer of scattering. Because of the time-consuming nature of these measurements, F(Q, t) is usually measured only at the first peak of the structure function S(Q), and its principal decay time is referred to as the α-relaxation time τ α . τ α is generally considered to describe the relaxation time of the bulk, which is related to viscosity and is controlled by the atomic cage around an atom. Here, through simulations on metallic liquids, we show that the α-relaxation time determined by scattering experiments does not purely reflect viscosity but is influenced by changes in spatial cooperativity. We also demonstrate that atomic caging is not exerted by the nearest neighbors but involves more cooperative atomic dynamics of the atomic medium-range order.

Glass transitions↗
DOE OSTI2022

Ultrafast perturbation of magnetic domains by optical pumping in a ferromagnetic multilayer

Ultrafast optical pumping of systems with spatially nonuniform magnetic textures is known to cause far-from-equilibrium spin transport effects, such as the broadening of domain-walls. Here, we study the dynamics of labyrinth domain networks in ferromagnetic CoFe/Ni multilayers subject to a femtosecond optical pump and find an ultrafast domain dilation by 6% within 1.6 ps. This surprising result is based on the unambiguous determination of a harmonically-related shift of ultrafast magnetic X-ray diffraction for the first- and third-order rings. Domain dilation is plausible from conservation of momentum arguments, whereby inelastic scattering from a hot, quasi-ballistic, radial current transfers momentum to the magnetic domains. Our results suggest a potentially rich variety of unexpected physical phenomena associated with far-from-equilibrium inelastic electron-magnon scattering processes in the presence of spin textures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spatiotemporal dynamics of fast electron heating in solid-density matter via XFEL

Abstract High-intensity, short-pulse lasers are crucial for generating energetic electrons that produce high-energy-density (HED) states in matter, offering potential applications in igniting dense fusion fuels for fast ignition laser fusion. High-density targets heated by these electrons exhibit spatially non-uniform and highly transient conditions, which have been challenging to characterize due to limitations in diagnostics that provide simultaneous high spatial and temporal resolution. Here, we employ an X-ray Free Electron Laser (XFEL) to achieve spatiotemporally resolved measurements at sub-micron and femtosecond scales on a solid-density copper foil heated by laser-driven fast electrons. Our X-ray transmission imaging reveals the formation of a solid-density hot plasma localized to the laser spot size, surrounded by Fermi degenerate, warm dense matter within a picosecond, and the energy relaxation occurring within the hot plasma over tens of picoseconds. These results validate 2D particle-in-cell simulations incorporating atomic processes and provide insights into the energy transfer mechanisms beyond current simulation capabilities. This work significantly advances our understanding of rapid fast electron heating and energy relaxation in solid-density matter, serving as a key stepping stone towards efficient high-density plasma heating and furthering the fields of HED science and inertial fusion energy research using intense, short-pulse lasers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Engineering Electrical Conductivity in Stable Zirconium-Based PCN-222 MOFs with Permanent Mesoporosity

Electrically conductive metal–organic frameworks (MOFs) featuring zirconium-based nodes are of great interest for electrochemical and optoelectronic applications due to their exceptional thermal and chemical stability, although the number of such MOFs remains limited. Here in this paper, we demonstrate that electron deficient molecules such as transition metal(IV) bis(dicarbollide) (TM IV (C 2 B 9 H 11 ) 2 ; TM IV CB for short, where TM = Ni, Pd, and Pt) or C 60 , spatially infiltrated in the microporous channels of the zirconium-based porphyrinic MOFs PCN-222 and PCN-222-Zn, generate highly stable and electrically conductive frameworks in which the mesopores remain accessible to other guests. Solid-state density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations indicate incorporation of Ni IV CB and C 60 in the MOF micropores is energetically and structurally feasible. Interestingly, in contrast to the free species, strong donor–acceptor interactions between the MOF and Ni IV CB restrain it to a single conformation. Calculated electronic structures and charge-hopping conduction probabilities illustrate that efficient charge-transfer (CT) from photoexcited linkers to the guest molecules facilitates charge hopping in the framework, making the MOFs electrically conductive. The donor–acceptor conjugates also enhance exciton dissociation at their heterojunctions, fostering the formation of long-lived electron-trapped states with potential utility for photo- and electrochemical devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Spatially Resolved Investigation of Mixed Valence and Insulator-to-Metal Transition in an Organic Salt

Using scanning tunneling microscopy/spectroscopy (STM/STS), we investigate the evolution of electronic structures across the boundaries of 7,7,8,8-tetracyanoquinodimethane (TCNQ) and K-TCNQ assemblies on a weakly interacting substrate. Despite the semiconducting/insulating nature of TCNQ (TCNQ 0 ) and K-TCNQ (TCNQ -1 ), a continuum metallic-like density of states extending deep (~1.5 nm) into the TCNQ assembly is observed near domain boundary. We attribute the formation of these states to the abrupt change of molecular valence, which perturbs the electrostatics of the junction and creates local electric fields as evidenced by the band bending near domain boundary. To the best of our knowledge, this study provides the first microscopic understanding of the crucial physics occurring near domain boundaries of mixed valence in K-TCNQ, or broadly speaking charge transfer complexes, which highlights these boundaries as potential “weak” points to initiate the electric field induced insulator-to-metal transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Accurate Prediction of Voltage of Battery Electrode Materials Using Attention-Based Graph Neural Networks

Performing first-principles calculations to discover electrodes’ properties in the large chemical space is a challenging task. While machine learning (ML) has been applied to effectively accelerate those discoveries, most of the applied methods ignore the materials’ spatial information and only use predefined features: based only on chemical compositions. Here, we propose two attention-based graph convolutional neural network techniques to learn the average voltage of electrodes. Our proposed methods, which combine both atomic composition and atomic coordinates in 3D-space, improve the accuracy in voltage prediction significantly when compared to composition-based ML models. The first model directly learns the chemical reaction of electrodes and metal ions to predict their average voltage, whereas the second model combines electrodes’ ML predicted formation energy (E form ) to compute their average voltage. Our E form -based model demonstrates improved accuracy in transferability from our subset of learned Li ions to Na ions. Moreover, we predicted the theoretical voltage of 10 Na x MPO 4 F (M = Ti, Cr, Fe, Cu, Mn, Co, and Ni) fluorophosphate battery frameworks, which are unavailable in the Material Project database. It could be shown that we can expect average voltages higher than 3.1 V from those Na battery frameworks except from the NaTiPO 4 F and TiPO 4 F pair of electrodes, which offer an average voltage of 1.32 V.

25 ENERGY STORAGE↗
DOE OSTI2023

A streamlined tandem tip-based workflow for sensitive nanoscale phosphoproteomics

Effective phosphoproteome of nanoscale sample analysis remains a daunting task, primarily due to significant sample loss associated with non-specific surface adsorption during enrichment of low stoichiometric phosphopeptide. We develop a tandem tip phosphoproteomics sample preparation method that is capable of sample cleanup and enrichment without additional sample transfer, and its integration with our recently developed SOP (Surfactant-assisted One-Pot sample preparation) and iBASIL (improved Boosting to Amplify Signal with Isobaric Labeling) approaches provides a streamlined workflow enabling sensitive, high-throughput nanoscale phosphoproteome measurements. This approach significantly reduces both sample loss and processing time, allowing the identification of >3000 (>9500) phosphopeptides from 1 (10) µg of cell lysate using the label-free method without a spectral library. It also enables precise quantification of ~600 phosphopeptides from 100 sorted cells (single-cell level input for the enriched phosphopeptides) and ~700 phosphopeptides from human spleen tissue voxels with a spatial resolution of 200 µm (equivalent to ~100 cells) in a high-throughput manner. The new workflow opens avenues for phosphoproteome profiling of mass-limited samples at the low nanogram level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Probing exciton diffusion dynamics in photosynthetic supercomplexes via exciton–exciton annihilation

Photosynthesis converts solar energy into chemical energy through coordinated energy transfer between light-harvesting complexes and reaction centers (RCs). Understanding exciton motion, particularly the exciton diffusion length, is essential for optimizing energy efficiency in photosystems. In this work, we combine intensity-cycling transient absorption spectroscopy with kinetic Monte Carlo (kMC) simulation to investigate exciton motion in the C2S2 photosystem II supercomplex of spinach. Using exciton–exciton annihilation, revealed in the fifth-order response, we experimentally estimate an exciton diffusion length of 10.9 nm based on a 3D normal diffusion model, suggesting the ability of excitons to traverse the supercomplex. However, kMC simulations reveal that exciton motion is sub-diffusive because of spatial constraints and the strong RC traps. An anomalous diffusion model analysis of the experimental data yields a diffusion length of 9.7 nm, while the simulated diffusion length is 7.4 nm. The variable exciton residence time across subunits, partly influenced by their connectivity to the trap, indicates inhomogeneous annihilation probability and suggests how plants balance efficient light harvesting with photoprotection. We also explore the influence of specific assumptions in the annihilation simulation, which are challenging to access in more complex environments, such as the thylakoid membrane. Our study provides a framework for studying exciton dynamics using exciton–exciton annihilation, which can be extended to understand the light-harvesting efficiencies of larger, more complex photosynthetic assemblies.

Zhang, Kunyan (ORCID:000000026830409X)↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2021

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗