Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Oxygen”

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 37 records · Page 2

Insights into oxygen diffusion in rare earth disilicate environmental barrier coatings

Environmental barrier coatings (EBC) are crucial for the use of SiC-based ceramic matrix composites in high-temperature combustion environments, yet knowledge of oxygen diffusion in these coatings is limited. This study investigates oxygen diffusion dynamics in the β-RE 2 Si 2 O 7 system to minimize oxygen penetration in rare earth disilicates. We analyze defect formation energy under varying oxygen conditions, identifying key diffusion mechanisms. In oxygen-rich environments, the most favorable neutral interstitial oxygen diffuses along the [110] direction. In oxygen-poor conditions, neutral oxygen vacancies rotate around Y and Si atoms, exhibiting a diffusivity of 6.59×10 −22 m 2 /s at 1500 K for β-Yb 2 Si 2 O 7 . Under intermediate oxygen levels, charged interstitial oxygen diffuses via concerted interstitialcy along the [001] direction with a diffusivity of 6.21×10 −17 m 2 /s. Additionally, alloying rare earth Y with Er and Yb increases diffusion barriers, contributing to improved EBC performance in extreme environments. The insights gained provides valuable guidance for designing robust coatings tailored to withstand extreme operational environments.

36 MATERIALS SCIENCE

Substantial oxygen loss and chemical expansion in lithium-rich layered oxides at moderate delithiation

Delithiation of layered oxide electrodes triggers irreversible oxygen loss, one of the primary degradation modes in lithium-ion batteries. However, the delithiation-dependent mechanisms of oxygen loss remain poorly understood. Here we investigate the oxygen non-stoichiometry in Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes as a function of Li content by using cycling protocols with long open-circuit voltage steps at varying states of charge. Surprisingly, we observe substantial oxygen loss even at moderate delithiation, corresponding to 2.5, 4.0 and 7.6 ml O 2 per gram of Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ after resting at upper capacity cut-offs of 135, 200 and 265 mAh g −1 for 100 h. Our observations suggest an intrinsic oxygen instability consistent with predictions of high oxygen activity at intermediate potentials versus Li/Li + . In addition, we observe a large chemical expansion coefficient with respect to oxygen non-stoichiometry, which is about three times greater than those of classical oxygen-deficient materials such as fluorite and perovskite oxides. Furthermore, our work challenges the conventional wisdom that deep delithiation is a necessary condition for oxygen loss in layered oxide electrodes and highlights the importance of calendar ageing for investigating oxygen stability.

Atomistic models

Ultrafast Oxygen Conduction in Sillén Oxychlorides

Oxygen ion conductors are crucial for enhancing the efficiency of various clean energy technologies, including fuel cells, solid oxide air batteries, electrolyzers, membranes, sensors, and more. In this study, a structure-similarity analysis of ≈62k oxygen-containing compounds identified the MBi 2 O 4 X (M = rare-earth element, X = halogen element) family as promising candidates for fast oxygen transport. Among these, LaBi 2 O 4 Cl is found as an ultrafast oxygen conductor with an ultralow migration barrier of 0.1 eV based on ab initio studies. Its 2D layered structure, featuring a “triple fluorite” layer, supports diffusion of both oxygen vacancies and interstitials. In addition to vacancy diffusion with a 0.1 eV barrier, ab initio studies show interstitial diffusion exhibits a modest barrier of 0.6–0.8 eV. Frenkel pairs are found to be the dominant defects in intrinsic LaBi 2 O 4 Cl, facilitating significant vacancy-mediated oxygen diffusion at elevated temperatures. With 2.8% oxygen vacancies, LaBi 2 O 4 Cl is predicted to achieve a conductivity of 0.3 S/cm at 25 °C in a single crystal. Experimental synthesis and characterization of polycrystalline LaBi 2 O 4 Cl and Sr-doped LaBi 2 O 4 Cl revealed conductivity exceeding that of YSZ and LSGM below 400 °C, with lower activation energies, achieving a total conductivity of 0.1−0.2 mS/cm at 300 °C. Here, while these results confirm its potential of fast oxygen transport, we suggest further experimental optimization of LaBi 2 O 4 Cl, including aliovalent doping and microstructure refinement, could significantly enhance its performance, facilitating fast oxygen conduction approaching room temperature.

Defects

Stationary Oxygen Vacancy Construction toward a Superior-Performance Ultrahigh Nickel Single-Crystal Cathode

Oxygen vacancies exert a complex and profound influence on the layered cathodes, especially those with ultrahigh nickel content. They can facilitate lithium-ion transport and enhance electronic conductivity, while aggressive oxygen vacancy formation causes structural degradation and electrolyte decomposition. Herein, taking ultrahigh nickel single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92) as a model material, we propose a pinning strategy to harness the benefits of oxygen vacancies while mitigating their detrimental effects. Through a carefully controlled thermal process, both oxygen vacancies and pinning atoms are successfully introduced into the surface region. The resulting anchored oxygen vacancies, capitalizing on their inherent advantages, improve conductivity and lithium-ion diffusion. Simultaneously, the neighboring pinning atoms effectively increase the migration barrier and suppress the adverse effects of these vacancies, including electrolyte decomposition and structural degradation during long-term electrochemical cycling. Consequently, oxygen vacancy-anchored single-crystal LiNi 0.92 Co 0.06 Mn 0.02 O 2 (SC-Ni92-OV) demonstrates significantly improved high-voltage electrochemical performance, with 86.16% capacity retention after 200 cycles at 4.6 V and 1 C in a half-cell and 90.71% after 300 cycles at 4.5 V and 1 C in a full cell. Furthermore, this study not only provides valuable insights into the chemistry of oxygen vacancy but also introduces a viable strategy for leveraging oxygen vacancies to achieve stable high-voltage performance in ultrahigh nickel single-crystal cathodes.

defects in solids

A Quasi-Ordered Mn-Rich Cathode with Highly Reversible Oxygen Anion Redox Chemistry

Anionic oxygen redox chemistry in Li-rich Mn-based layer oxide cathodes represents a transformative approach for boosting the energy density of next-generation lithium-ion batteries. However, conventional oxygen redox reactions often induce oxygen dimerization at high voltages, leading to irreversible lattice oxygen loss and a rapid voltage fade. Herein, we achieve highly reversible oxygen redox chemistry through a new quasi-ordered structural design that incorporates both intra- and interlayer cation disorder configurations. This unique structure significantly enhances lattice oxygen stability, effectively stabilizes oxidized oxygen, and inhibits the formation of peroxo- or superoxol-like species, thereby enabling anionic redox reactions to proceed reversibly even at deep delithiation states. The quasi-ordered design mitigates irreversible phase transitions and preserves the structural integrity throughout extended cycling. Consequently, the proposed cathode demonstrates exceptional cyclability with negligible capacity and voltage fade, retaining 99% capacity and 98% average voltage after long-term cycling. Finally this work provides fresh insights into addressing issues related to lattice oxygen instabilities and reforming strategies for developing long-life, high-energy-density anionic redox cathode materials for advanced batteries.

chemical structure

Influence of Oxygen Flow and Stoichiometry on Optical Properties and Damage Resistance of Hafnium Oxide Thin Films

Hafnium oxide (HfO 2 ) is predominantly used as a high-index material in multi-layer dielectric coatings for high-peak- and high-average-power lasers, but laser damage often initiates within the HfO 2 layers despite their wide bandgap. Oxygen deficiency during deposition can introduce vacancy-related sub-bandgap states and absorptive defects, lowering damage resistance. This study investigates how oxygen flow during HfO 2 deposition with ion beam sputtering (IBS) affects its stoichiometry, defect formation, and nanosecond laser-induced damage threshold (LIDT) and whether single-layer trends predict multilayer performance. Single layers were deposited at varying oxygen flows, characterized for optical and structural properties, and tested for the LIDT at 1064 nm and 355 nm. Increasing oxygen flow drove the layer toward near-stoichiometric HfO 2 , reduced the refractive index, and altered the density of surface pinhole-like features. The single-layer LIDT at 355 nm increased with oxygen, whereas the 1064 nm LIDT was comparatively less sensitive to oxygen flow, consistent with the wavelength-dependent roles of absorptive precursors and microstructural defects. In contrast, a HfO 2 -based high-reflector (HR) showed a higher LIDT at lower oxygen flow, indicating that the family of damage precursors changes between single layers and multilayers; in stacks, structural properties such as stress, gas entrapment and thermal dissipation may outweigh the isolated absorptive defects found in single layers. These results demonstrate that the optimal oxygen flow condition depends on both LIDT wavelength and film architecture. We identified, for single layers, a 15–35 sccm window for maximizing the 1064 nm LIDT and a high-flow optimum (45 sccm) for the 355 nm LIDT and, for 355 nm HR stacks, a distinct lower-flow regime (~10 sccm).

Optics and optical instruments

Modeling Oxygen Partial Pressure in Solid Oxide Electrolysis Cells: The Microstructure Effect

Oxygen partial pressure is an important thermodynamic state variable that affects both the performance and degradation of solid oxide electrolysis cells. In this work, a 3D model developed from Virkar’s 1D model has been applied to reconstructed and synthetic microstructures of Ni-YSZ-GDC-LSCF cell. The effect of the microstructures, including the thickness of the YSZ and GDC layer, and the compositions of the hydrogen and oxygen electrode, on the distribution of oxygen partial pressure has been investigated. The results show that the maximum oxygen partial pressure may occur on the interface between oxygen electrode and GDC layer or between GDC layer and YSZ layer depending on the rate of oxygen ion exchange between GDC and YSZ. A thicker GDC layer lowers the maximum oxygen partial pressure in the cell, while a thicker YSZ layer lowers the maximum oxygen partial pressure in the hydrogen electrode. In addition, the Ni:YSZ ratio and porosity also affects the maximum partial pressure. These findings provide insights on mitigating degradation in solid oxide electrolysis cell by tuning the microstructures.

Lei, Yinkai

Electrode and Microstructure Dependence of Oxygen Diffusion in Ferroelectric Hafnium Zirconium Oxide Thin Films

Hafnia-based ferroelectrics hold promise to reduce energy demand for computing by enabling compute-in-memory and as non-volatile memories. The ferroelectric phase in this material system is, in part, stabilized by oxygen vacancies. While oxygen vacancies may be a necessity for phase stability, they limit device endurance through diffusion and accumulation into conducting channels. Herein, it is shown that oxygen diffusion is spatially variable within individual grains of ferroelectric hafnium zirconium oxide (HZO). Using 18 O tracers and finite difference modeling, it is shown that grain boundaries and regions near electrode interfaces allow for relatively rapid oxygen diffusion, with values as much as 10 4 larger than the grain cores. Further, the selection of electrode material affects the diffusion coefficients across all microstructural regions. HZO films in contact with TiN electrodes result in more oxygen-deficient HZO films and higher oxygen diffusion coefficients. Tungsten electrodes result in fewer vacancies and lower diffusion coefficients. Diffusion activation energy differences between the HZO with the two electrodes is reconciled by differing populations of charged and uncharged oxygen vacancies. This insight into the local vacancy populations and diffusion pathways provides a platform for designing hafnia-based films, deposition processes, and integration strategies to reduce vacancy gradients and improve performance.

36 MATERIALS SCIENCE

Elongated Fullerene Unlocks Oxygen-Resistant Perovskite Solar Minimodules

State-of-the-art encapsulation cannot prevent the permeation of oxygen into perovskite solar cells (PSCs). Here, we report the finding that p-i-n structured PSCs degrade quickly under oxygen exposure. Oxidation of C60, rather than perovskites, dominates the degradation. Chemical absorption of oxygen by C60 is identified, which impairs its electron transport property. We modify the electron transport layers to address the instability under oxygen exposure. The elongated fullerene C70 is found to react with oxygen orders of magnitude slower than C60. In combination with a compact SnO2 buffer fabricated by atomic layer deposition, which can slow down the oxygen diffusion, the resulting unencapsulated PSCs with C70 retained 90% of their initial efficiency after 1-sun illumination in pure oxygen for 1,200 h at 70 degrees C, improving stability by hundreds of times. Testing of unencapsulated perovskite minimodules of different perovskite compositions with C70 gives extrapolated lifetimes of 17-41 years at 50 degrees C.

14 SOLAR ENERGY

Oxygen Vacancy Evolution at Li x V 2 O 5 /LiPON Solid State Electrochemical Interfaces Using Depth Resolved Cathodoluminescence Spectroscopy

The formation of oxygen vacancies at buried LiPON/ Li x V 2 O 5 interfaces has been observed on a near-nanometer scale and nondestructively using depth-resolved cathodoluminescence spectroscopy (DRCLS) and interfacial markers. Before electrochemical cycling, as-deposited LiPON/Li x V 2 O 5 exhibits a 1.6 eV defect optical emission, which density functional theory calculations identify as originating from oxygen vacancies. This defect appears first within a few nanometers of the buried LiPON/Li x V 2 O 5 interface without cycling, indicating that spontaneous O diffusion from the Li x V 2 O 5 lattice into LiPON may have caused these interface-localized oxygen vacancy defects. DRCLS measured the intensity and spatial distribution of this oxygen vacancy signal as a function of electrochemical cycling in a LiPON/Li x V 2 O 5 half-cell, showing oxygen vacancy signal increasing and moving deeper into the electrode with increased cycle number. Significant electrochemical irreversibility was also observed, with poor Coulombic efficiency and a 15% drop in capacity over 50 cycles. Theoretical simulations predict that the presence of oxygen vacancies increases the energy barrier for lithium diffusion significantly, indicating that this aggregation of oxygen vacancies could be another battery degradation mechanism accompanying lithiation induced phase changes.

36 MATERIALS SCIENCE

Modulating Surface Anionic Redox Chemistry toward Highly Stable Li-Rich Cathodes with Negligible Oxygen Loss

Low initial Coulombic efficiency and severe capacity/voltage fading during cycling caused by serious irreversible oxygen release, especially in the initial cycle, and resultantly induced unstable electrode/electrolyte interfacial chemistry, largely prohibit the commercial application of high-capacity Li-rich layered oxide cathodes (LLOs). Here, in this work, a dual reductive gas interface cotreatment strategy is applied to regulate the lattice oxygen redox activity and reversibility with a multiple defective structure design including Li/O/TM (TM = transition metal) vacancies and the intrinsic TM doping as well as a full-surface protective layer, which can suppress the irreversible TM migration and then undesirable phase transformation, resisting the corrosion of electrolyte during cycling effectively. Importantly, the introduced reversible SO 3 2- /SO 4 2- redox couple that provides extra capacity compensation could alleviate the distortion of oxygen-central octahedral structure and structural collapse caused by immoderate oxygen oxidation. Thus, the lattice oxygen redox chemistry is optimized, with negligible oxygen loss during the initial cycle. And the designed AS-LLO cathode with greatly enhanced structure stability shows high-capacity retentions of 99.2% at 0.3C after 100 cycles and 82.4% even after 1000 cycles at 5C. This work provides a guideline for manipulating the oxygen redox chemistry to achieve long-lifespan Li-rich layered oxide cathodes for high-energy-density lithium batteries.

25 ENERGY STORAGE

A Method for Rapid and Precise Triple Oxygen Isotope Measurements via High-Temperature Conversion to CO Followed by Nickel-Catalyzed CO to CO 2 Conversion and Laser Spectroscopy

Triple oxygen isotopic compositions ( 16 O, 17 O, 18 O) have conventionally been measured via isotope ratio mass spectrometry using O 2 as an analyte. Conversion of sample oxygen to O 2 typically utilizes fluorination chemistry or catalytic equilibration between CO 2 and O 2 . Recently, laser spectroscopy has become a viable alternative for triple oxygen isotope (Δ' 17 O) measurements due to its ease and rapid throughput. Laser spectrometers are currently available for Δ' 17 O analysis of either H 2 O or CO 2 as the analyte gas. So far, these instruments have been used to measure Δ' 17 O of water, carbonate (CO 2 liberated by acid digestion), and atmospheric CO 2 samples. Here, we present a new method for high-precision Δ' 17 O analysis of CO 2 via tunable infrared laser direct absorption spectroscopy that is compatible with a wider range of geochemically important materials. This approach involves converting sample oxygen to CO 2 in two steps. First, the sample oxygen is liberated and reduced to CO by high-temperature conversion at 1450 °C in the presence of excess elemental carbon. Then, CO is catalytically converted to CO 2 over hot nickel at 350 °C. The conversion process is rapid (10 to 30 min) and quantitative. Spectroscopic Δ' 17 O analysis of the resulting CO 2 takes approximately 45 min. By measuring several oxygen isotope standards, we demonstrate that the method is precise (1σ = 12 per meg for procedural replicates) and accurate (within 11 per meg of previously reported values). The method can be applied to most pyrolytic materials where quantitative oxygen conversion is attainable, such as sulfate, phosphate, nitrate, and oxide minerals, water, and organic molecules.

Ellis, Nicholas M. [University of California, Berk

Carbon-Mediated Oxygen Vacancy Creation at Hematite Interfaces

Nanoscale iron oxides (e.g., hematite (a-Fe 2 O 3 )) have unique properties, such as enhanced chemical reactivity and high surface area, when compared with their bulk counterparts. These nanoscale surfaces can be more reactive due to the presence of defects (e.g., oxygen vacancies). In this work, we probed the surface chemistry of bulk and nanoscale hematite via X-ray photoelectron spectroscopy, electron microscopy, and powder X-ray diffraction. Oxygen exposure and vacuum annealing experiments were conducted to add or remove oxygen vacancies and remove adventitious carbon. In the absence of the oxygen annealing step, vacuum annealing resulted in partial reduction of Fe(III) to Fe(II) on all hematite surfaces. This was a size-dependent effect with the extent of reduction increasing as the crystallite size decreased. In addition, the atomic concentrations of carbon increased on all iron oxide surfaces after vacuum annealing. Oxygen annealing almost completely removed carbon from sample surfaces, and no Fe(III) reduction was observed in the absence of carbon. Under these conditions, the results reveal that carbonaceous material enhances oxygen vacancy formation, which then facilitates the reduction of Fe(III) on hematite surfaces. We provide new insights into the mechanisms of Fe(III) reduction on both bulk and nanoscale hematite surfaces and establish the major role of carbon in oxygen vacancy formation.

Zengotita, Frances E. [University of Notre Dame, I

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato

Oxygen redox activities governing high-voltage charging reversibility of Ni-rich layered cathodes

A schematic illustration indicative of the improved structural reversibility through stabilizing oxidized oxygen by Al doping. The chemical reactions and phase transitions at high voltages are generally considered to determine the electrochemical properties of high-voltage layered cathodes such as Ni-rich rhombohedral oxides. Even if significantly higher SOCs (states-of-charge) are utilized above the capability of transition metal redox (primarily Ni and Co), the effect of oxygen redox on Ni-rich rhombohedral oxides still looks mysterious thereby necessitating research that can clarify the relationship between redox reactions and phase transitions. Here, we performed a comprehensive and comparative study of the cationic and anionic redox reactions, as well as the structural evolution of a series of commercial Ni-rich layered oxides with and without Al doping. We combined the results from X-ray spectroscopy, operando electrochemical mass spectrometry, and neutron diffraction with electrochemical properties and thereby revealed the different oxygen redox activities associated with structural and electrochemical degradations. We reveal that Al doping suppresses the irreversible oxygen release but enhances the lattice oxygen oxidization. With this modulated oxygen redox activity, the Ni-rich layered oxides’ notorious H2–H3 structural phase transition becomes highly reversible. Our findings disentangle the different oxygen redox activities during high-voltage cycling and clarify the role of dopants in the Ni-rich layered oxides in terms of structural and electrochemical stability, shedding light on the future direction of optimizing layered cathode materials for safer high energy-density secondary batteries.

Lee, Gi-Hyeok

Modifications in the charge trap landscape in Hf 0.5 Zr 0.5 O 2 as a function of oxygen vacancy concentration observed with photoemission electron microscopy

Oxygen vacancies in Hf x Zr (1-x) O 2 (HZO) both contribute to stabilization of the ferroelectric orthorhombic phase and promote leakage pathways that limit the endurance of devices based on the material. Here, for this reason, the defect states of oxygen vacancies were investigated using photoemission electron microscopy (PEEM) and photoluminescence spectroscopy (PL), as their concentration was varied via ex situ laser exposure. Following a controlled oxygen vacancy reduction via visible (2.54 eV) laser dosing of HZO, deep-ultraviolet (DUV, 5.82 eV) PEEM was used to spatially probe the resulting mid-gap defect states and work function. Work function was found to increase monotonically with the laser-induced reduction in oxygen vacancy concentration culminating in a total increase near 70 meV. The change implies a Fermi level shift toward the valence band as the total available electron-filled charge states are reduced with the removal of oxygen vacancies. A reduction in charge states is corroborated by the observed lessening of both photoemission and photoluminescence intensities after laser dosing. The deduced position of the Fermi level is within a band of near-conduction band defect states produced by oxygen vacancies that are linked to endurance limiting leakage currents. Together, these results directly identify the primary role of oxygen vacancies on the defect states in HZO while demonstrating that laser exposure can be used for their modification.

42 ENGINEERING

Irreversible oxygen poisoning: Modeling the initial water dissociation kinetics on δ-Pu(111) and δ-Pu(100) through density functional theory

In this work, the initial kinetics of water dissociation on two facets of δ-plutonium, δ-Pu(111) and δ-Pu(100), are explored through density functional theory in order to understand how water dissociation occurs on these facets. We explored the dissociation of water via the formation of hydroxyls, atomic hydrogen, and atomic oxygen species on each facet. We calculate low energetic barriers for water to split to adsorbed hydrogen and hydroxyl species at 0.19 eV for δ-Pu(111) and 0.07 eV for δ-Pu(100). The hydroxyl has a barrier of 0.64 and 0.37 eV to cleave the hydrogen–oxygen bond on δ-Pu(111) and δ-Pu(100), respectively. Due to the highly exergonic adsorption free energy of atomic oxygen of −2.10 eV, the metallic surfaces are found to be fully covered in oxygen, even with the inclusion of oxygen lateral interactions. When combined with the reaction thermodynamics, this free energy of adsorption forms a molecular oxygen desorption barrier greater than 9.51 eV (918 kJ/mol). These results, combined with simulated temperature programed desorption spectra, indicate that oxygen formed via water dissociation on δ-Pu (111) and (100) facets induces an irreversible poison that prevents further reaction of water directly on metallic plutonium surfaces, most likely due to the strong hybridization of the Pu and O valence states. Therefore, these results imply that another mechanism is responsible for the continuously experimentally measured water dissociation, which is most likely due to the Pu oxide.

36 MATERIALS SCIENCE

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling