Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field validation”

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

Higher-derivative relations between scalars and gluons

We extend the covariant color-kinematics duality introduced by Cheung and Mangan to effective field theories. We focus in particular on relations between the effective field theories of gluons only and of gluons coupled to bi-adjoint scalars. Maps are established between their respective equations of motion and between their tree-level scattering amplitudes. An additional rule for the replacement of flavor structures by kinematic factors realizes the map between higher-derivative amplitudes. As an example of new relations, the pure-gluon amplitudes of mass dimension up to eight, featuring insertions of the F 3 and F 4 operators which satisfy the traditional color-kinematics duality, can be generated at all multiplicities from just renormalizable amplitudes of gluons and bi-adjoint scalars. We also obtain closed-form expressions for the kinematic numerators of the dimension-six gluon effective field theory, which are valid in D space-time dimensions. Finally, we find strong evidence that this extended covariant color-kinematics duality relates the (DF) 2 +YM(+Φ 3 ) theories which, at low energies, generate infinite towers of operators satisfying the traditional color-kinematics duality, beyond aforementioned F 3 and F 4 ones.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska North Slope – Producer Responses and Operational Lessons Learned

The US Department of Energy and Hilcorp Alaska, LLC are jointly cosponsoring a field pilot project entitled “First Ever Field Pilot on Alaska’s North Slope to Validate the Use of Polymer Floods for Heavy Oil Enhanced Oil Recovery (EOR).” The objective of this pilot is to perform a field experiment to validate the use of an advanced polymer flooding technology to unlock the vast heavy oil resources on Alaska North Slope (ANS). Two pairs of horizontal injection and production wells have been deployed in an isolated fault block of the Schrader Bluff heavy oil reservoir at the Milne Point Field to conduct a polymer flood pilot which was designed to acquire field performance data to optimize polymer flood design in the Schrader Bluff heavy oil reservoirs on ANS. Hydrolyzed polyacrylamide (HPAM) polymer injection started on August 28, 2018 using a custom made polymer mixing and pumping unit. Polymer solution properties such as concentration, viscosity, and filter ratio, have been closely monitored and adjusted as necessary. Operational procedures have been improved to ensure adequate injectivity and polymer propagation through the formation. Oil and water production is tested frequently to monitor water cut performance, and produced water is sampled and analyzed for polymer content. Tracers were pumped pre and post polymer injection startup to compare polymer flood breakthrough timing and sweep efficiency. Polymer injectivity of the horizontal wells drilled in the Schrader Bluff NB sand is found to be sufficient to replace the reservoir production voidage although there were some declines as high viscosity polymer swept the near wellbore region. Production data show significant reduction in water cut and the oil production rate has stabilized. No polymer production has been confirmed from the two horizontal producers after 22 months of polymer injection into the two supporting horizontal injectors. Tracer data show very low sweep efficiencies with waterflood alone, whereas significant improvement has been achieved with polymer flood. Many operational lessons have been learned and improvements have been made regarding polymer mixing, pumping, and onsite quality control, particularly in a challenging Arctic environment. Adequate polymer injection rate can be achieved with diligent monitoring and quality assurance of the properties of the injected polymer solution.

Ning, Samson↗
DOE OSTI2023

Sensitivity of Optical Satellites to Estimate Windthrow Tree-Mortality in a Central Amazon Forest

Windthrow (i.e., trees broken and uprooted by wind) is a major natural disturbance in Amazon forests. Images from medium-resolution optical satellites combined with extensive field data have allowed researchers to assess patterns of windthrow tree-mortality and to monitor forest recovery over decades of succession in different regions. Although satellites with high spatial-resolution have become available in the last decade, they have not yet been employed for the quantification of windthrow tree-mortality. Here, we address how increasing the spatial resolution of satellites affects plot-to-landscape estimates of windthrow tree-mortality. We combined forest inventory data with Landsat 8 (30 m pixel), Sentinel 2 (10 m), and WorldView 2 (2 m) imagery over an old-growth forest in the Central Amazon that was disturbed by a single windthrow event in November 2015. Remote sensing estimates of windthrow tree-mortality were produced from Spectral Mixture Analysis and evaluated with forest inventory data (i.e., ground true) by using Generalized Linear Models. Field measured windthrow tree-mortality (3 transects and 30 subplots) crossing the entire disturbance gradient was 26.9 ± 11.1% (mean ± 95% CI). Although the three satellites produced reliable and statistically similar estimates (from 26.5% to 30.3%, p < 0.001), Landsat 8 had the most accurate results and efficiently captured field-observed variations in windthrow tree-mortality across the entire gradient of disturbance (Sentinel 2 and WorldView 2 produced the second and third best results, respectively). As expected, mean-associated uncertainties decreased systematically with increasing spatial resolution (i.e., from Landsat 8 to Sentinel 2 and WorldView 2). However, the overall quality of model fits showed the opposite pattern. We suggest that this reflects the influence of a relatively minor disturbance, such as defoliation and crown damage, and the fast growth of natural regeneration, which were not measured in the field nor can be captured by coarser resolution imagery. Our results validate the reliability of Landsat imagery for assessing plot-to-landscape patterns of windthrow tree-mortality in dense and heterogeneous tropical forests. Satellites with high spatial resolution can improve estimates of windthrow severity by allowing the quantification of crown damage and mortality of lower canopy and understory trees. However, this requires the validation of remote sensing metrics using field data at compatible scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microphysics of liquid water in sub-10 nm ultrafine aerosol particles

Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.

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

Superuniversality from disorder at two-dimensional topological phase transitions

We investigate the effects of quenched randomness on topological quantum phase transitions in strongly interacting two-dimensional systems. We focus first on transitions driven by the condensation of a subset of fractionalized quasiparticles (“anyons”) identified with “electric charge” excitations of a phase with intrinsic topological order. All other anyons have nontrivial mutual statistics with the condensed subset and hence become confined at the anyon condensation transition. Furthermore, using a combination of microscopically exact duality transformations and asymptotically exact real-space renormalization group techniques applied to these two-dimensional disordered gauge theories, we argue that the resulting critical scaling behavior is “superuniversal” across a wide range of such condensation transitions and is controlled by the same infinite-randomness fixed point as that of the 2D random transverse-field Ising model. We validate this claim using large-scale quantum Monte Carlo simulations that allow us to extract zero-temperature critical exponents and correlation functions in (2+1)D disordered interacting systems. We discuss generalizations of these results to a large class of ground-state and excited-state topological transitions in systems with intrinsic topological order as well as those where topological order is either protected or enriched by global symmetries. When the underlying topological order and the symmetry group are Abelian, our results provide prototypes for topological phase transitions between distinct many-body localized phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Emulation and detection of physical faults and cyber-attacks on building energy systems through real-time hardware-in-the-loop experiments

The increasing use of remote or mobile access, integrated wearable technologies, data exchange, and cloud-based data analytics in modern smart buildings is steering the building industry towards open communication technologies. The increased connectivity and accessibility could lead to more cyber-attacks in smart buildings. On the other hand, physical faults (e.g., HVAC -heating, ventilation, and air-conditioning faults) may have similar adverse impacts as those from the cyber-attacks on building energy systems, such as occupant discomfort, energy wastage, and equipment downtime. However, current physical behavior-based anomaly detection methods fail to differentiate between cyber-attacks and physical faults in building energy systems. Moreover, the challenge in collecting real-world threat data with ground truth has led researchers to rely on numerical models with user-defined assumptions, which may not accurately reflect real-world conditions due to the lack of in-situ experimental datasets. To address these challenges and gaps, this paper presents a flexible hardware-in-the-loop (HIL) testbed for generating cyber-attack and physical fault datasets and demonstrating threat detection algorithms in a real building automation system (BAS) environment. This testbed combines hardware (i.e., real BAS with local HVAC controllers and a physical network) with software (i.e., high-fidelity models to represent behaviors of building envelope and HVAC energy systems), enabling emulations of realistic threats. Five HIL experiments, including one baseline without any threats, two with physical faults, and two with cyber-attacks, were conducted to generate datasets containing detailed network traffic and system states. A joint classification framework, incorporating a network analyzer and a physical HVAC fault detector, was proposed to automatically detect cyber-physical abnormalities on BAS at both the network and the physical HVAC levels. The network analyzer comprises a conditional random fields (CRF) based command validator and a statistics-based detection strategy. The fault detector employs a weather and schedule-based pattern matching and feature-based principal component analysis (WPM-FPCA) method. Evaluation of the classification using four metrics from the multi-class confusion matrix revealed an average accuracy of 90.2%, recall of 89.7%, precision of 88.5% and F1-score of 89.2%. Finally, these results demonstrate that the proposed joint classification framework can effectively differentiate between specific types of cyber-attacks (e.g., device reinitialization attack, network Denial-of-Service attack) and physical faults (e.g., air handling unit operational fault, cooling coil valve stuck) in real time for improved building energy management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗
DOE OSTI2025

Liquid Sorption-Enhanced Haber–Bosch Process

The use of a liquid sorbent in a traditional Haber-Bosch process enables significant improvements in energy efficiency and potential cost savings for arguably the most important chemical process on the planet. The approach presented in this report employs an incompressible liquid sorbent that absorbs and releases ammonia (NH 3 ) under specific conditions. To achieve this, we investigate reactions of ammonia and pure phosphoric acid (H 3 PO 4 , PA), which rapidly neutralize to form an equilibrated solution of monoammonium phosphate (MAP) and diammonium phosphate (DAP) that functions as a reversible and regenerable sorbent. Through intimate contact of the gas-phase Haber-Bosch reaction mixture with this liquid absorbent, complete equilibrium uptake may be achieved in an appropriately sized separator, and facile separation occurs through the use of independent liquid and gas phases. Following depressurization and release of the ammonia product, only the incompressible fluid needs to be repressurized and returned to the reactor. This study documents proof-of-concept absorption and desorption experiments carried out in 75 mL batch reactors, predominantly charged with precise MAP and DAP mixtures that equilibrate at process-relevant temperatures and pressures. We then assemble the first thermodynamic relationships that underlie this advantaged separation strategy, validated by reactive force field (ReaxFF) interatomic potential simulations, and benchmarked with traditional separation routes via process modeling and technoeconomic analysis. The scale of energy consumption in the century-old Haber-Bosch process is massive, and the elegant liquid sorption approach reported here offers opportunities to enhance its energy efficiency for the next frontier of ammonia synthesis.

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

Backbone Stitching in Bottlebrush Copolymer Mesodomains and the Impact of Side Chain Crystallization

We synthesized bottlebrush statistical copolymers (BSCPs) having poly­(ethylene oxide) (PEO) and poly­(dimethylsiloxane) (PDMS) side chains attached to a polynorbornene backbone. Small-angle X-ray scattering analysis showed that for densely grafted BSCPs, the scattering length density gradually transitions between the PEO and PDMS domains. For loosely grafted BSCPs, the polymer backbone formed a distinct mesodomain, with a lower electron and mass density than both the PEO and PDMS domains. The bottlebrush backbone essentially “stitches” the PEO and PDMS side chains, looping back and forth from the PEO to PDMS domains with the backbone segments oriented normal to the domain interfaces. Self-consistent field theory (SCFT) calculations validated the stitching of the backbone driven by the microphase separation of PEO and PDMS, along with a strong segmental order of the side chains in the melt. The reduced birefringence upon PEO crystallization suggests the disruption of the strong segmental order by the crystallization. Both the static intrinsic and the form birefringences of the BSCPs decreased upon PEO crystallization. Solid-state NMR confirmed the rigidity of PEO crystallites and the bottlebrush backbone. Self-assembly of BSCPs containing polyhedral oligomeric silsesquioxane (POSS) pendent groups was also evaluated by X-ray scattering, showing the formation of lamellar microdomains that inhibited POSS crystallization.

Hu, Mingqiu↗
DOE OSTI2025

Chain Flexibility and Structure of a Polyimide Copolymer: Revisiting the Freely Rotating Chain Model

Poly(4–4′-oxydiphenylene-pyromellitimide)-based polyimides─trade name Kapton─have wide-ranging engineering applications owing to their thermal and mechanical stability, but little is known about underlying chain-level characteristics. While theoretical models have conceptualized Kapton as inflexible polycyclic rods separated by freely rotating diphenyl ether hinge groups, the model’s core predictions remain untested and subtleties of the relaxation behavior are missed, which atomistic modeling can resolve. To these ends, we examine all-atom Kapton structures in crystalline and glassy amorphous configurations using a DFT-validated class II force field. Constructing amorphous configurations is challenging, as the fused-ring-containing backbone has slow relaxation dynamics and scaling suggestive of entanglements even in oligomers. In conclusion, we find larger backbone rearrangements of the linear polycyclic segments about ether groups that are consistent with the rod-hinge picture on the monomer scale, whereas a ring rotation analysis suggests partially flexible rod-like segments and involves multiple facile rotational relaxation modes.

Liesen, Nicholas T. [Lawrence Livermore National L↗
DOE OSTI2024

On-demand nanoengineering of in-plane ferroelectric topologies

Hierarchical assemblies of ferroelectric nanodomains, so-called super-domains, can exhibit exotic morphologies that lead to distinct behaviours. Controlling these super-domains reliably is critical for realizing states with desired functional properties. Here we reveal the super-switching mechanism by using a biased atomic force microscopy tip, that is, the switching of the in-plane super-domains, of a model ferroelectric Pb 0.6 Sr 0.4 TiO 3 . We demonstrate that the writing process is dominated by a super-domain nucleation and stabilization process. A complex scanning-probe trajectory enables on-demand formation of intricate centre-divergent, centre-convergent and flux-closure polar structures. Correlative piezoresponse force microscopy and optical spectroscopy confirm the topological nature and tunability of the emergent structures. The precise and versatile nanolithography in a ferroic material and the stability of the generated structures, also validated by phase-field modelling, suggests potential for reliable multi-state nanodevice architectures and, thereby, an alternative route for the creation of tunable topological structures for applications in neuromorphic circuits.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Ferroelectric Al1 − xBxN–GaN heterostructures

This report demonstrates Al0.93B0.7N thin films grown epitaxially on n-type GaN (0002)/c-plane sapphire substrates by reactive magnetron sputtering at 300 °C. At 200 nm film thickness, the Al0.93B0.07N layers exhibit partially relaxed substrate-induced epitaxial strain, a 0.16° wide (0002) rocking curve, in-plane crystallographic registry, and sub-nanometer surface roughness. Electrically, the stack shows robust hysteresis over three frequency decades, a remanent polarization of ∼125 μC/cm2, a strongly frequency dependent coercive field, highly uniform dc leakage currents, and endurance >106 field cycles. This report validates possibilities for ferroelectric nitride integration into conventional III-nitride heterostructures with high crystalline fidelity, high electrical resistivity, and persistent hysteresis. Such materials are synthesizable at thermal budgets and temperatures compatible with back-end-of-the-line boundary conditions.

Physics↗
DOE OSTI2020

Adversarial super-resolution of climatological wind and solar data

Accurate and high-resolution data reflecting different climate scenarios are vital for policy makers when deciding on the development of future energy resources, electrical infrastructure, transportation networks, agriculture, and many other societally important systems. However, state-of-the-art long-term global climate simulations are unable to resolve the spatiotemporal characteristics necessary for resource assessment or operational planning. We introduce an adversarial deep learning approach to super resolve wind velocity and solar irradiance outputs from global climate models to scales sufficient for renewable energy resource assessment. Using adversarial training to improve the physical and perceptual performance of our networks, we demonstrate up to a 50 × resolution enhancement of wind and solar data. In validation studies, the inferred fields are robust to input noise, possess the correct small-scale properties of atmospheric turbulent flow and solar irradiance, and retain consistency at large scales with coarse data. An additional advantage of our fully convolutional architecture is that it allows for training on small domains and evaluation on arbitrarily-sized inputs, including global scale. We conclude with a super-resolution study of renewable energy resources based on climate scenario data from the Intergovernmental Panel on Climate Change’s Fifth Assessment Report.

14 SOLAR ENERGY↗
DOE OSTI2021

Scalar-tensor mixing from icosahedral inflation

Abstract We study the mixed scalar-tensor two-point function in icosahedral inflation. Within the regime of validity of the effective field theory, this has to be perturbatively small; in particular, much smaller than the scalar spectrum. However, it can be much bigger that the tensor spectrum itself. We discuss observational implications for the CMB temperature-polarization spectrum.

Astronomy & Astrophysics↗
DOE OSTI2025

Elastic cross section is entanglement entropy

We present universal relations between entanglement entropy, which quantifies the quantum correlation between subsystems, and the cross section, which is the primary observable for high-energy particle scattering, by employing a careful formulation of wave packets for the incoming particles. For 2-to-2 elastic scattering with no initial entanglement and subdividing the system along particle labels, we show that both the Rényi and Tsallis entropies in the final states are directly proportional to the elastic cross section in units of the transverse size for the initial wave packets, which is then interpreted as the elastic scattering probability. The relations do not depend on the underlying dynamics of the quantum field theory and are valid to all orders in coupling strengths. Furthermore, computing quantum correlations between momentum and nonkinematic data leads to entanglement entropies expressed as various semi-inclusive elastic cross sections. Our result gives rise to a novel “area law” for entanglement entropy in a two-body system. Published by the American Physical Society 2025

Low, Ian (ORCID:0000000275709597)↗
DOE OSTI2023

Scaling laws for two-dimensional dendritic crystal growth in a narrow channel

Here, we investigate analytically and computationally the dynamics of two-dimensional needle crystal growth from the melt in a narrow channel. Our analytical theory predicts that, in the low supersaturation limit, the growth velocity $\textit{V}$ decreases in time $\textit{t}$ as a power law $V ~t^{–2/3}$, which we validate by phase-field and dendritic-needle-network simulations. Simulations further reveal that, above a critical channel width $Λ ≈ 5l_D$, where $l_D$ is the diffusion length, needle crystals grow with a constant $V < V_s$, where $V_s$ is the free-growth needle crystal velocity, and approaches $V_s$ in the limit $Λ \gg l_D$.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing the Critical Nucleus Size in the Metal-Insulator Phase Transition of VO 2

In a first-order phase transition, critical nucleus size governs nucleation kinetics, but the direct experimental test of the theory and determination of the critical nucleation size have been achieved only recently in the case of ice formation in supercooled water. The widely known metal-insulator phase transition (MIT) in strongly correlated VO 2 is a first-order electronic phase transition coupled with a solid-solid structural transformation. It is unclear whether classical nucleation theory applies in such a complex case. Here we directly measure the critical nucleus size of the MIT by introducing size-controlled nanoscale nucleation seeds with focused ion irradiation at the surface of a deeply supercooled metal phase of VO 2 . The results compare favorably with classical nucleation theory and are further explained by phase-field modeling. This Letter validates the application of classical nucleation theory as a parametrizable model to describe phase transitions of strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

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↗