Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Explanability”

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 145 records · Page 8

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗

Opposing Hysteresis Patterns in Flow and Outflow Macroscopic Fundamental Diagrams and Their Implications

Two key aggregated traffic models are the relationship between average network flow and density (known as the network or flow macroscopic fundamental diagram [flow-MFD]) and the relationship between trip completion and density (known as network exit function or the outflow-MFD [o-FMD]). The flow- and o-MFDs have been shown to be related by average network length and average trip distance under steady-state conditions. However, recent studies have demonstrated that these two relationships might have different patterns when traffic conditions are allowed to vary: the flow-MFD exhibits a clockwise hysteresis loop, while the o-MFD exhibits a counter-clockwise loop. One recent study attributes this behavior to the presence of bottlenecks within the network. The present paper demonstrates that this phenomenon may arise even without bottlenecks present and offers an alternative, but more general, explanation for these findings: a vehicle’s entire trip contributes to a network’s average flow, while only its end contributes to the trip completion rate. This lag can also be exaggerated by trips with different lengths, and it can lead to other patterns in the o-MFD such as figure-eight patterns. A simple arterial example is used to demonstrate this explanation and reveal the expected patterns, and they are also identified in real networks using empirical data. Then, simulations of a congestible ring network are used to unveil features that might increase or diminish the differences between the flow- and o-MFDs. Finally, more realistic simulations are used to confirm that these behaviors arise in real networks.

42 ENGINEERING↗

Chemical classification program synthesis using generative artificial intelligence

Accurately classifying chemical structures is essential for cheminformatics and bioinformatics, including tasks such as identifying bioactive compounds of interest, screening molecules for toxicity to humans, finding non-organic compounds with desirable material properties, or organizing large chemical libraries for drug discovery or environmental monitoring. However, manual classification is labor-intensive and difficult to scale to large chemical databases. Existing automated approaches either rely on manually constructed classification rules, or are deep learning methods that lack explainability. This work presents an approach that uses generative artificial intelligence to automatically write chemical classifier programs for classes in the Chemical Entities of Biological Interest (ChEBI) database. These programs can be used for efficient deterministic run-time classification of SMILES structures, with natural language explanations. The programs themselves constitute an explainable computable ontological model of chemical class nomenclature, which we call the ChEBI Chemical Class Program Ontology (C3PO). We validated our approach against the ChEBI database, and compared our results against deep learning models and a naive SMARTS pattern based classifier. C3PO outperforms the naive classifier, but does not reach the performance of state of the art deep learning methods. However, C3PO has a number of strengths that complement deep learning methods, including explainability and reduced data dependence. C3PO can be used alongside deep learning classifiers to provide an explanation of the classification, where both methods agree. The programs can be used as part of the ontology development process, and iteratively refined by expert human curators.

Artificial Intelligence↗

SoyFACE Fumigation Data Files

This data set is related to the SoyFACE experiments, which are open-air agricultural climate change experiments that have been conducted since 2001. The fumigation experiments take place at the SoyFACE farm and facility in Champaign County, Illinois during the growing season of each year, typically between June and October. - The "SoyFACE Plot Information 2001 to 2021" file contains information about each year of the SoyFACE experiments, including the fumigation treatment type (CO2, O3, or a combination treatment), the crop species, the plots (also referred to as 'rings' and labeled with numbers between 2 and 31) used in each experiment, important experiment dates, and the target concentration levels or 'setpoints' for CO2 and O3 in each experiment. - This data set includes files with minute readings of the fumigation levels ( "SoyFACE 1-Minute Fumigation Data Files" folder) from the SoyFACE experiments. The "Soyface 1-Minute Fumigation Data Files" folder contains sub-folders for each year of the experiments, each of which contains sub-folders for each ring used in that year's experiments. This data set also includes hourly data files for the fumigation experiments ( "SoyFACE Hourly Fumigation Data Files" folder) created from the 1-minute files, and hourly ambient/weather data files for each year of the experiments ( "Hourly Weather and Ambient Data Files" folder). The ambient CO2 and O3 data are collected at SoyFACE, and the weather data are collected from the SURFRAD and WARM weather stations located near the SoyFACE farm. - The "Fumigation Target Percentages" file shows how much of the time the CO2 and O3 fumigation levels are within a 10 or 20 percent margin of the target levels when the fumigation system is turned on. - The "Matlab Files" folder contains custom code (Aspray, E.K.) that was used to clean the "SoyFACE 1-Minute Fumigation Data" files and to generate the "SoyFACE Hourly Fumigation Data" and "Fumigation Target Percentages" files. Code information can be found in the "SoyFACE Hourly Fumigation Data Explanation" file. - Finally, the " * Explanation" files contain information about the column names, units of measurement, and other pertinent information for each data file. * NOTE: We have identified some files in the “SoyFACE 1-Minute Fumigation Data Files” folder in our SoyFACE data set submission that were not downloaded properly - the files were present in the folder, but the actual files were empty. V3 ensures that there are no longer any empty files in the data set.

agricultural↗

Understanding oxidation of Fe-Cr-Al alloys through explainable artificial intelligence

Abstract The oxidation resistance of FeCrAl based on alloying composition and oxidizing conditions is predicted using a combinatorial experimental and artificial intelligence approach. A neural network (NN) classification model was trained on the experimental FeCrAl dataset produced at GE Research. Furthermore, using the SHapley Additive exPlanations (SHAP) explainable artificial intelligence (XAI) tool, we explore how the NN can showcase further material insights that are unavailable directly from a black-box model. We report that high Al and Cr content forms protective oxide layer, while Mo in FeCrAl creates thick unprotective oxide scale that is vulnerable to spallation due to thermal expansion. Graphical abstract

Materials Science↗

Parsimonious Inference Information-Theoretic Foundations for a Complete Theory of Machine Learning (CIS-LDRD Project 218313 Final Technical Report)

This work examines how we may cast machine learning within a complete Bayesian framework to quantify and suppress explanatory complexity from first principles. Our investigation into both the philosophy and mathematics of rational belief leads us to emphasize the critical role of Bayesian inference in learning well-justified predictions within a rigorous and complete extended logic. The Bayesian framework allows us to coherently account for evidence in the learned plausibility of potential explanations. As an extended logic, the Bayesian paradigm regards probability as a notion of degrees of truth. In order to satisfy critical properties of probability as a coherent measure, as well as maintain consistency with binary propositional logic, we arrive at Bayes' Theorem as the only justifiable mechanism to update our beliefs to account for empiracle evidence. Yet, in the machine learning paradigm, where explanations are unconstrained algorithmic abstractions, we arrive at a critical challenge: Bayesian inference requires prior belief. Conventional approaches fail to yield a consistent framework in which we could compare prior plausibility among the infinities of potential choices in learning architectures. The difficulty of articulating well-justified prior belief over abstract models is the provinence of memorization in traditional machine learning training practices. This becomes exceptionally problematic in the context of limited datasets, when we wish to learn justifiable predictions from only a small amount of data.

97 MATHEMATICS AND COMPUTING↗

Review of "ML for Surface Complexation Model Development"

The presentation given by Jadallah Zouabe on “ML for Surface Complexation Model Development”was extremely informative, but lacked some explanation of details that made the latter half hard to follow as someone outside of the field. He began his presentation with an anecdote to spilling milk or ink on the floor and seeing it spread, cleverly connecting that to the spread of contamination in the real world. He then explained the importance of being able to model contamination, as it allows specialists to see the potential effects on different areas and environments in the coming years. Heal so made it clear that the main issue trying to be addressed here is how expensive modelling is and that their goal is to make computationally cheaper models. The introduction was overall very well formulated and attention-grabbing, but it seemed to take up too much of the presentation time, as after this point the presentation lacked details and explanations.

54 ENVIRONMENTAL SCIENCES↗

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗

Understanding the Origin of the Hadron Mass within the Standard Model

Understanding the origin of the hadron mass, which constitutes 99% of our visible universe, is one of the central goals of nuclear physics. Although the Higgs mechanism provides mass for the fundamental building blocks of matter, it can only contribute less than 2% of the proton mass. The vast majority of the proton mass is believed to come from the strong force that tightly binds quarks and gluons (collectively called partons) together as described by Quantum Chromodynamics (QCD). The mass that emerges as a consequence of the strong interactions within QCD is commonly denoted as Emergent Hadronic Mass (EHM). Understanding how the nucleon mass emerges in QCD is a prerequisite to an explanation of how the Universe came into being, therefore it is of utmost importance and one of the key questions to be addressed by the future Electron-Ion Collider (EIC). When it comes to light mesons, particularly pions, the mass decomposition is drastically different. Since the pion is naturally massless in the chiral limit, the majority of its observed mass needs to come from other mechanisms within QCD. Any successful explanation for the EHM needs to be able to reconcile both the heavy proton mass and the very light pion mass (~15% of proton mass) simultaneously. The EHM theories have direct measurable implications on the description of the internal structure of the hadron, i.e., how the partons distribute inside the hadron. Precise measurement of the parton distribution functions (PDFs) will provide necessary experimental verifications and constraints of potential EHM mechanisms. We propose to carry out a comprehensive study of the poorly known pion PDFs at the AMBER experiment at CERN. The measurement will provide vital input to constrain the global analysis of the pion PDFs, which are still based on limited data obtained more than 30 years ago. The proposed pion measurement is the only direct measurement achievable within this decade, which could lead to a future major meson structure program parallel to the EIC’s proton structure measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗

X-Ray Generation for Diffraction Experiments on Z

Explanation of the concept and benefit of X-ray diffraction in the context of material sciences on Z. Explanation why we use ZBL and ZPW for these experiments and what the progress on X-ray generation for this purpose is.

Geissel, Matthias↗

Investigation of hydrometeorological influences on reservoir releases using explainable machine learning methods

Long short-term memory (LSTM) networks have demonstrated successful applications in accurately and efficiently predicting reservoir releases from hydrometeorological drivers including reservoir storage, inflow, precipitation, and temperature. However, due to its black-box nature and lack of process-based implementation, we are unsure whether LSTM makes good predictions for the right reasons. In this work, we use an explainable machine learning (ML) method, called SHapley Additive exPlanations (SHAP), to evaluate the variable importance and variable-wise temporal importance in the LSTM model prediction. In application to 30 reservoirs over the Upper Colorado River Basin, United States, we show that LSTM can accurately predict the reservoir releases with NSE ≥ 0.69 for all the considered reservoirs despite of their diverse storage sizes, functionality, elevations, etc. Additionally, SHAP indicates that storage and inflow are more influential than precipitation and temperature. Moreover, the storage and inflow show a relatively long-term influence on the release up to 7 days and this influence decreases as the lag time increases for most reservoirs. These findings from SHAP are consistent with our physical understanding. However, in a few reservoirs, SHAP gives some temporal importances that are difficult to interpret from a hydrological point of view, probably because of its ignorance of the variable interactions. SHAP is a useful tool for black-box ML model explanations, but the hydrological processes inferred from its results should be interpreted cautiously. More investigations of SHAP and its applications in hydrological modeling is needed and will be pursued in our future study.

54 ENVIRONMENTAL SCIENCES↗

Structure–Property Linkage in Alloys Using Graph Neural Network and Explainable Artificial Intelligence

Deep learning tools have recently shown significant potential for accelerating the prediction of microstructure–property linkage in materials. While deep neural networks like convolution neural networks (CNNs) can extract physics information from 3D microstructure images, they often require a large network architecture and substantial training time. In this research, we trained a graph neural network (GNN) using phase field generated microstructures of Ni-Al alloys to predict the evolution of mechanical properties. We found that a single GNN is capable of accurately predicting the strengthening of Ni-Al alloys with microstructures of varying sizes and dimensions, which cannot otherwise be done with a CNN. Additionally, GNN requires significantly less GPU utilization than CNN and offers more interpretable explanation of predictions using saliency analysis as features are manually defined in the graph. We also utilize explainable artificial intelligence tool Bayesian Inference to determine the coefficients in the power law equation that governs coarsening of precipitates. Overall, our work demonstrates the ability of the GNN to accurately and efficiently extract relevant information from material microstructures without having restrictions on microstructure size or dimension and offers an interpretable explanation.

Chemistry↗

Regulation of Solar Wind Electron Temperature Anisotropy by Collisions and Instabilities

Abstract Typical solar wind electrons are modeled as being composed of a dense but less energetic thermal “core” population plus a tenuous but energetic “halo” population with varying degrees of temperature anisotropies for both species. In this paper, we seek a fundamental explanation of how these solar wind core and halo electron temperature anisotropies are regulated by combined effects of collisions and instability excitations. The observed solar wind core/halo electron data in ( β ∥ , T ⊥ / T ∥ ) phase space show that their respective occurrence distributions are confined within an area enclosed by outer boundaries. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures and β ∥ is the ratio of parallel thermal energy to background magnetic field energy. While it is known that the boundary on the high- β ∥ side is constrained by the temperature anisotropy-driven plasma instability threshold conditions, the low- β ∥ boundary remains largely unexplained. The present paper provides a baseline explanation for the low- β ∥ boundary based upon the collisional relaxation process. By combining the instability and collisional dynamics it is shown that the observed distribution of the solar wind electrons in the ( β ∥ , T ⊥ / T ∥ ) phase space is adequately explained, both for the “core” and “halo” components.

Yoon, Peter H. (ORCID:0000000181343790)↗

Dark Photon Search at the Short-Baseline Near Detector

Neutrino physics has long been a key field in elementary particle physics, both enhancing our understanding of the Standard Model (SM) and raising new questions. Among these are the so-called "Short-Baseline Anomalies" observed by neutrino experiments, particularly the MiniBooNE experiment at Fermilab, which detected an excess of low-energy electron-like events. In recent years, beyond Standard Model (BSM) explanations have been proposed to address this anomaly, with a focus on neutrino beam-related processes. A novel interpretation involving a dark-sector explanation was recently suggested, introducing a vector portal that connects the SM and dark sectors through a new interaction mediated by a bosonic particle, the Dark Photon. This work investigates the production of dark-sector particles, specifically Dark Photons, in the Booster Neutrino Beam (BNB) at Fermilab and their potential detection at the Short-Baseline Near Detector (SBND). The BNB produces mesons which decay into Dark Photons, detectable via their decay into electron-positron pairs at SBND. By exploiting the temporal structure of the neutrino beam, we propose a method to isolate Dark Photon signals from neutrino backgrounds using time-delayed event detection. In this thesis, the sensitivity of SBND to Dark Photons is assessed using a three-year exposure, demonstrating that SBND has the potential to significantly improve current experimental constraints on Dark Photons. This analysis provides a promising avenue for future dark sector searches in neutrino experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. A systematic evaluation of the interaction between coastal catchment properties and carbon dioxide (CO2) emission by soil respiration is significant for assessing carbon dynamics and predicting the future trajectory of atmospheric CO2 concentrations in coastal TAIs. The soil CO2 efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation index (kNDVI), climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, there is no one-size-fits-all approach to modeling carbon fluxes in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

54 ENVIRONMENTAL SCIENCES↗

Deep Learning Estimation of Daily Ground–Level NO 2 Concentrations from Remote Sensing Data

The limited number of nitrogen dioxide (NO 2 ) surface measurements calls for the development of highly accurate approaches to estimating surface NO 2 concentrations. In this study, we leverage a new satellite instrument, the TROPOspheric Monitoring Instrument (TROPOMI), along with other predictor variables, to estimate daily surface NO 2 concentrations over Texas in 2019. We use the deep convolutional neural network (Deep-CNN), an advanced deep learning algorithm, to obtain estimates and achieve a correlation coefficient (R) of 0.91, an index of agreement (IOA) of 0.95, and a mean absolute bias (MAB) of 1.75 ppb in surface NO 2 estimation. Additionally, we leverage a novel approach, SHapley Additive exPlanations (SHAP), to describe how Deep-CNN understands each predictor variable. The SHAP results show that the Deep-CNN model has an advanced understanding of the dataset, revealing that TROPOMI closely captures levels of NO 2 . In addition, we show the superiority of our Deep-CNN model at estimating surface NO 2 over other well-known machine learning and regression models in the field, including the support vector machines (SVM), random forest (RF), and multiple linear regression (MLR). Although SVM and RF show strong capabilities at estimating surface NO 2 concentrations, their accuracy is inferior to that of the Deep-CNN model, ranking second and third in model accuracy in this study. The MLR, however, shows a poor ability at NO 2 estimation and ranks last among all models. Furthermore, testing the impact of sample size on model performance, we also show that, compared to other models, Deep-CNN needs more samples to trigger its strength at surface NO 2 estimation.

54 ENVIRONMENTAL SCIENCES↗

Contrasting Carbon–Water–Energy Dynamics in Perennial and Annual Bioenergy Agroecosystems Using Eddy Covariance and Interpretable Machine Learning

Understanding how agroecosystems respond to environmental variability is fundamental to predicting productivity and sustainability under a changing climate. We analyzed 55 site-years of high-frequency eddy covariance observations from five agroecosystems—two perennial grasses (miscanthus and switchgrass), two annual rotation systems (maize–soybean and sorghum–soybean), and a restored native prairie—to examine ecosystem-scale carbon, water, and energy fluxes. Using an interpretable machine-learning framework with regression tree ensembles, Shapley Additive Explanations, and Accumulated Local Effects, we quantified how environmental and temporal factors regulate gross primary productivity (GPP), evapotranspiration (ET), water-use efficiency, and the Bowen ratio. Perennials exhibited stronger physiological buffering and maintained fluxes across a broader range of temperature and moisture conditions, reflecting deeper rooting and persistent canopy cover. Annuals, in contrast, showed greater short-term variability and stronger coupling to atmospheric demand, with GPP and ET declining rapidly under low humidity or soil moisture. Differences in temperature sensitivity of Bowen ratio further revealed that perennials sustained proportionally greater sensible heat flux under cool conditions, whereas annuals exhibited constrained energy exchange when evaporative demand was low. Together, these results demonstrate that crop life cycle and canopy structure are fundamental determinants of ecosystem-scale carbon–water–energy coupling. By integrating long-term flux observations with interpretable machine learning, this study identifies the environmental drivers that shape agroecosystem function and highlights how conversion from annual to perennial feedstocks can enhance climatic resilience and alter land–atmosphere energy feedbacks. These findings provide a data-driven basis for improving crop and Earth-system models and for guiding bioenergy landscape design under future climate scenarios.

Accumulated Local Effects↗