Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “link layer”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Physical inconsistencies in the representation of the ocean heat-carbon nexus in simple climate models

The Ocean Heat-Carbon Nexus, linking ocean heat and carbon uptake, is crucial for understanding climate responses to cumulative carbon dioxide (CO 2 ) emissions and to net zero CO 2 emissions. It results from a suite of processes involving the exchange of heat and carbon across the sea-air interface as well as their storage below the mixed layer and redistribution by the ocean large-scale circulation. The Ocean Heat and Carbon Nexus is assumed to be consistently represented across two modelling platforms used in the latest IPCC assessments: the Earth System Models (ESMs) and the Simple Climate Models (SCMs). However, our research shows significant deficiencies in state-of-the-art SCMs in replicating the ocean heat-carbon nexus of ESMs due to a crude treatment of the ocean thermal and carbon cycle coupling. With one SCM, we show that a more realistic heat-to-carbon uptake ratio exacerbates the projected warming by 0.1 °C in low overshoot scenarios and up to 0.2 °C in high overshoot scenarios. It is therefore critical to explore how SCMs’ physical inconsistencies, such as the representation of the ocean heat-carbon nexus, can affect future warming projections used in climate assessments, not just by SCMs in Working Group 3 but also by ESMs in Working Group 1 via SCM-driven emission-to-concentration translation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Topology-Dependent Performance of Free-Space Photonic Quantum Networks Under Noise

Photonic quantum communication enables secure and high-fidelity information transfer beyond classical limits, with direct relevance to emerging quantum networks operating in free-space environments. While physical-layer models of depolarizing noise, Gamma–Gamma turbulence statistics, entanglement swapping, and decoy-state QKD security bounds are individually well established, prior work typically treats these components in isolation or under fixed network assumptions. In this work, we develop a unified topology-aware analytical framework that simultaneously integrates free-space optical link budgets, turbulence-induced visibility degradation, depolarizing qubit noise, multi-hop entanglement cascade dynamics, teleportation fidelity thresholds, CHSH nonlocality certification, and asymptotic decoy-state secret key rate bounds across star, mesh, and ring graph structures. Rather than introducing new physical channel models, we demonstrate that identical physical links exhibit fundamentally different end-to-end performance once embedded within different network topologies. Mesh architectures minimize visibility cascade through hop-count reduction but incur quadratic hardware scaling. Star topologies minimize link count but concentrate noise and synchronization overhead at the hub. Ring configurations offer linear hardware scaling with multiplicative fidelity degradation. The results establish topology as a first-order design parameter in near-term free-space quantum networks operating without full quantum repeater infrastructures. While motivated by distributed multi-agent architectures, the framework applies broadly to terrestrial, airborne, and satellite-assisted photonic quantum communication systems.

QKD↗
DOE OSTI2026

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗
DOE OSTI2025

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spin-Canting-Induced Band Reconstruction in the Dirac Material Ca 1 – x Na x MnBi 2

The ternary A MnBi 2 ( A is alkaline as well as rare-earth atom) materials provide an arena for investigating the interplay between low-dimensional magnetism of the antiferromagnetic MnBi layers and the electronic states in the intercalated Bi layers, which harbor relativistic fermions. We report on a comprehensive study of the optical properties and magnetic torque response of Ca 1–x Na x MnBi 2 . Our findings give evidence for a spin canting occurring at T s ~50–100 K. Thus, with the support of first-principles calculations we establish a direct link between the spin canting and the reconstruction of the electronic band structure, having immediate implications for the spectral weight reshuffling in the optical response, signaling a partial gapping of the Fermi surface, and the dc transport properties below T s .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Promotion of superconductivity in magic-angle graphene multilayers

Graphene moiré superlattices show an abundance of correlated insulating, topological, and superconducting phases. Whereas the origins of strong correlations and nontrivial topology can be directly linked to flat bands, the nature of superconductivity remains enigmatic. We demonstrate that magic-angle devices made of twisted tri-, quadri-, and pentalayer graphene placed on monolayer tungsten diselenide exhibit flavor polarization and superconductivity. We also observe insulating states in the tril- and quadrilayer arising at finite electric displacement fields. As the number of layers increases, superconductivity emerges over an enhanced filling-factor range, and in the pentalayer it extends well beyond the filling of four electrons per moiré unit cell. Our results highlight the role of the interplay between flat and more dispersive bands in extending superconducting regions in graphene moiré superlattices.

Science & Technology - Other Topics↗
DOE OSTI2024

Processes Associated with Boundary Layer Cloud Ice Phase Precipitation in the High Southern Latitudes (Final Technical Report)

This is a final technical report for DOE Award No. DE-SC0018995. Using data collected during the MARCUS deployment of the 2nd ARM mobile facility, this project sought to address the following scientific questions: 1. What are the conditions that govern phase partitioning in supercooled boundary layer clouds over the high latitude Southern Oceans? 2. To what extent are supercooled mixed phase clouds responsible for the high cloud fractions over the Southern Ocean and do these mixed phase clouds contribute to the model biases in absorbed shortwave radiation? 3. Can the cloud properties and phase partitioning over the Southern Ocean be linked to the underlying aerosol environment? 4. How do the aerosol properties at the surface vary spatially and seasonally over the high latitude Southern Ocean and what processes control these features? 5. How well defined are the atmospheric compositional boundaries between Southern Ocean marine, Antarctic free tropospheric, and Antarctic surface airmasses, and how does this influence the chemistry and physics of the Southern Ocean atmosphere? In this project we published three peer review journal articles that address these questions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Photostability of Ambient-Processed, Conjugated Polymer Electrochromic Devices Encapsulated by Bioderived Barrier Films

Polymer-based electrochromic devices (ECDs) are a promising technology for enabling low-voltage, disposable displays, yet are currently limited by photo-oxidative bleaching of the active materials. Here, renewable barrier films composed of cellulose and chitin are presented as an alternative to poly(ethylene terephthalate) (PET) for encapsulating ECDs. To assess barrier film effectiveness, lateral ECDs composed of poly(3,4-propylenedioxythiophene-(CH 2 OEtHx) 2 ) (P(ProDOT)) active layers were constructed and encapsulated with a multilayer barrier consisting of chitin nanofibers and cellulose nanocrystals spray cast onto a cellulose acetate substrate (oxygen transmission rate (OTR) = 29 cm 3 m –2 day –1 ), a commercially available PET film (OTR = 8.5 cm 3 m –2 day –1 ), and a high-performance PET-Al 2 O 3 multilayer barrier film (OTR < 1 cm 3 m –2 day –1 ). The photodegradation of the P(ProDOT) active layer was determined by measuring the evolution of the colorimetric contrast (ΔE*) and switching speeds as a function of light exposure (100 mW cm –2 , AM 1.5 G light). Photodegradation was found to proceed at a similar rate for all encapsulated devices (roughly 10 times more slowly than unencapsulated devices), highlighting the opportunity for replacing petroleum packaging with bioderived barrier films. Analysis of the switching kinetics, the shifts in optical absorbance, and evidence of chemical degradation indicate that both photochemical breakdown of the electrolyte and cross-linking of the P(ProDOT) active material are key drivers for loss of device performance when oxygen flux to the active material is limited. Furthermore, pathways toward better understanding photodegradation are then proposed with sustainability in mind for future ECD design.

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

Origin of Phase Separation in Ni-Rich Layered Oxide Cathode Materials During Electrochemical Cycling

In intercalation materials, the kinetics and uniformity of mass transport across the nanocrystalline domains dictate the structural reversibility and transport capability at the macroscopic level. (De)intercalation-induced interlayer disintegrations exhibit anisotropic crystallite size change. Due to the anisotropic mass transport mechanism, separated phases are inherently crystallographically oriented. One such material is LiNi $1–x–y$ Mn $x$ Co $y$ O 2 , which plays a pivotal role in advanced Li-ion batteries but suffers from severe phase inhomogeneities under fast charge or electrochemical aging. Here, using operando synchrotron techniques, we probe the mechanistic origins of the compositional and orientational-dependent phase separations during the electrochemical cycling of LiNi 0.8 Mn 0.1 Co 0.1 O 2 by comprehensive analysis of both in-plane and out-of-plane reflections. In the H2/H3 phase regime, in-plane domain propagation occurs due to increased covalency despite the severe decay of interlayer crystallographic order, resulting in the change of crystalline domain shape from 3D spheres to two-dimensional nanosheets. The crystallographically selective XRD line splitting is linked to the geometry of the facets as mass transfers along the ab plane. In conclusion, this work provides mechanistic insights into crystallographic orientation-dependent phase inhomogeneity under fast charge and extended cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomic structure and phase assemblages in novel M-(N)-A-S-H materials

This paper investigates the atomic structure and phase assemblages in new sodium-stabilized magnesium aluminosilicate hydrate (M-(N)-A-S-H) cementitious binders. Results indicate that in the absence of Ca{sup 2+}, Mg{sup 2+} promotes a binder atomic structure of Si-Al tetrahedral sheets and octahedral Mg sheets with hydrated Na{sup +} cations likely in the interlayer sites similar to trioctahedral micas (phyllosilicates). NMR studies verify the incorporation of Al in tetrahedral silicate sheets. XRD demonstrates the ability of these regions to nucleate and form zeolites (i.e., sodalite) as well as the formation of Mg-Al layered double hydroxide (LDH) phases (i.e., meixnerite), which is expected due to high concentrations of Mg and Al. TGA results indicate that M-(N)-A-S-H possesses chemically bound water and hydroxyl units similar to other Mg binders. These results evince the critical role of Mg to form unique atomic structures and durability-linked phases in low-calcium alkali-activated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimizing Non-Terrestrial Hybrid RF/FSO Links With Reinforcement Learning: Navigating Through Clouds

In the pursuit of ubiquitous broadband connectivity, there has been a significant shift towards the vertical expansion of communication networks into space, particularly through the exploitation of low Earth orbit (LEO) satellite constellations, which are favored for their relatively low latency. However, this approach faces many challenges that need to be addressed, including atmospheric turbulence, high path loss, and dynamic cloud formations. High-altitude pseudo-satellites (HAPS) have emerged as promising relaying layers between LEO satellites and ground stations, enhancing coverage, latency, and direct terrestrial user connectivity. While radio frequency (RF) bands suffer from congestion and limited bandwidth, free space optical (FSO) communications offer higher data rates, but are susceptible to misalignment and weather-induced signal degradation. To address these challenges, a hybrid RF/FSO approach has been proposed to take advantage of both technologies by dynamic switching between RF and FSO based on propagation channel conditions. This paper introduces a reinforcement learning-based algorithm designed to optimize the trajectory of HAPS, maneuver around cloudy areas, and seamlessly switch between the RF and FSO communication modes to maximize the achievable capacity. The proposed approach aims to maximize system performance by intelligently adapting to environmental conditions and offering a promising solution for next-generation space communication networks.

actor-critic algorithm↗
DOE OSTI2026

The Strength of Coupling to the Southern Ocean Modulates Tropical Eastern Pacific Variability and Forced Response

Despite rising global-mean temperatures, large parts of the Southern Ocean and tropical eastern Pacific Ocean have cooled during the satellite era. These regions may be linked by teleconnections, with Southern Ocean cooling contributing to tropical eastern Pacific cooling. We demonstrate that, on average, state-of-the-art Earth system models (ESMs) underestimate the magnitude of interaction between the Southern Ocean and tropical eastern Pacific Ocean. The strength of the teleconnection is shown to be mediated by the magnitude of the positive cloud–sea surface temperature (SST) feedback in the subtropical eastern Pacific Ocean and the strength of the wind–evaporation–SST (WES) feedback. We link excessive precipitation in the tropical Pacific south of the equator to the strength of the Southern Ocean–eastern Pacific teleconnection. This model bias, known as the double intertropical convergence zone (ITCZ), is shown to be related to erroneous convection south of the equator, weakened cross-equatorial trade winds, and unfavorable meteorological conditions for marine boundary layer subtropical clouds. We postulate there is a two-way interaction, in which a double-ITCZ occurs with weaker cloud–SST and WES feedbacks, which in turn impact local SSTs and amplify the double-ITCZ. Models with a stronger Southern Ocean to tropical Pacific teleconnection tend to exhibit more multidecadal variability in the Walker circulation, ITCZ, and west–east equatorial SST gradient, as well as greater delayed warming in the tropical eastern Pacific Ocean resulting from delayed Southern Ocean warming under greenhouse gas forcing. These results provide insight into why ESMs struggle to replicate observed tropical Pacific temperature trend patterns and point to ITCZ location as a key target for improvement in future model development.

Climate change↗
DOE OSTI2020

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Giant room-temperature nonlinearities in a monolayer Janus topological semiconductor

Abstract Nonlinear optical materials possess wide applications, ranging from terahertz and mid-infrared detection to energy harvesting. Recently, the correlations between nonlinear optical responses and certain topological properties, such as the Berry curvature and the quantum metric tensor, have attracted considerable interest. Here, we report giant room-temperature nonlinearities in non-centrosymmetric two-dimensional topological materials—the Janus transition metal dichalcogenides in the 1 T’ phase, synthesized by an advanced atomic-layer substitution method. High harmonic generation, terahertz emission spectroscopy, and second harmonic generation measurements consistently show orders-of-the-magnitude enhancement in terahertz-frequency nonlinearities in 1 T’ MoSSe (e.g., > 50 times higher than 2 H MoS 2 for 18 th order harmonic generation; > 20 times higher than 2 H MoS 2 for terahertz emission). We link this giant nonlinear optical response to topological band mixing and strong inversion symmetry breaking due to the Janus structure. Our work defines general protocols for designing materials with large nonlinearities and heralds the applications of topological materials in optoelectronics down to the monolayer limit.

42 ENGINEERING↗
DOE OSTI2021

Three-Dimensional Free Gas Flow Focuses Basin-Wide Microbial Methane to Concentrated Methane Hydrate Reservoirs in Geological System

Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.

58 GEOSCIENCES↗
DOE OSTI2024

Unraveling the electronic structure and magnetic transition evolution across monolayer, bilayer, and multilayer ferromagnetic Fe 3 GeTe 2

Two-dimensional (2D) van der Waals (vdW) magnets have sparked widespread attention due to their potential in spintronic applications as well as in fundamental physics. Ferromagnetic vdW compound Fe 3 GeTe 2 (FGT) and its Ga variants have garnered significant interest due to their itinerant magnetism, correlated states, and high magnetic transition temperature. Experimental studies have demonstrated the tunability of FGT’s Curie temperature, T C , through adjustments in quintuple layer numbers (QL) and carrier concentrations, n. However, the underlying mechanism remains elusive. In this study, we employ molecular beam epitaxy (MBE) to synthesize 2D FGT films down to 1 QL with precise layer control, facilitating an exploration of the band structure and the evolution of itinerant carrier density. Angle-resolved photoemission spectroscopy (ARPES) reveals significant band structure changes at the ultra-thin limit, while first-principles calculations elucidate the band evolution from 1 QL to bulk, largely governed by interlayer coupling. Additionally, we find that n is intrinsically linked to the number of QL and temperature, with a critical value triggering the magnetic phase transition. Our findings underscore the pivotal role of band structure and itinerant electrons in governing magnetic phase transitions in such 2D vdW magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Holographic View of Mixed-State Symmetry-Protected Topological Phases in Open Quantum Systems

We establish a holographic duality between d -dimensional mixed-state symmetry-protected topological (mSPT) phases and ( d + 1 ) -dimensional subsystem symmetry-protected topological (SSPT) states. Specifically, we show that the reduced density matrix of the boundary layer of a ( d + 1 ) -dimensional SSPT state with subsystem symmetry S and global symmetry G corresponds to a d -dimensional mSPT phase with strong S and weak G symmetries. Conversely, we demonstrate that the wave function of an SSPT state can be constructed by replicating the density matrix of the corresponding lower-dimensional mSPT phase. This mapping links the density matrix in lower dimensions to the entanglement properties of higher-dimensional wave functions, providing an approach for analyzing nonlinear quantities and quantum information metrics in mixed-state systems. Our duality offers a new perspective for studying intrinsic mSPT phases that are unique to open quantum systems, without pure-state analogues. We show that strange correlators and twisted Rényi- N correlators can diagnose these nontrivial phases and we explore their connection to strange correlators in pure-state SSPT phases. Furthermore, we discuss several implications of this holographic duality, including a method for preparing intrinsic mSPT phases through the duality. Published by the American Physical Society 2025

Sun, Shijun (ORCID:0000000168880030)↗
DOE OSTI2023

Linking Stream Chemistry to Subsurface Redox Architecture

As drinking-water scarcity grows worldwide, we need to improve predictions of the quantity and quality of our water resources. An overarching problem for model improvement is that we do not know the geological structure of aquifers in sufficient detail. In this work, we demonstrate that mineral-water reactions imprint structure in the subsurface that impacts the flow and transport of some chemical species. Specifically, pyrite, a ubiquitous mineral, commonly oxidizes and depletes in the upper layers of the weathering profile in most humid watersheds, only remaining at depths of meters. We hypothesize that variations in concentrations (C) of pyrite-derived sulfate released into rivers as a function of discharge (q) reflect the rate-limiting step and depth of pyrite-oxidizing layers. We found that log C - log q behaviors thus differ in small and large watersheds in the Susquehanna River Basin as well as in selected watersheds in the Western United States. Although coal mining changes pyrite oxidation from closed to open system with respect to O 2 , patterns in stream chemistry as a function of discharge are consistent with deep and shallow pyrite oxidation zones in small and large watersheds respectively. Therefore, understanding the subsurface patterns of mineral reactions and how they affect the architecture of aquifers will elucidate patterns of changing river chemistry and our ability to manage water resources in the future under accelerated land use and climate change.

54 ENVIRONMENTAL SCIENCES↗